USGS Science⌕ Search

Geology topics

W. H. Bakun

Publications and source records attributed to W. H. Bakun.

At least 37 records · Page 2Linked to original sources

Association of earthquakes and faults in the San Francisco Bay area using Bayesian inference

Bayesian inference provides a method to use seismic intensity data or instrumental locations, together with geologic and seismologic data, to make quantitative estimates of the probabilities that specific past earthquakes are associated with specific faults. Probability density functions are constructed for the location of each earthquake, and these are combined with prior probabilities through Bayes' theorem to estimate the probability that an earthquake is associated with a specific fault. Results using this method are presented here for large, preinstrumental, historical earthquakes and for recent earthquakes with instrumental locations in the San Francisco Bay region. The probabilities for individual earthquakes can be summed to construct a probabilistic frequency–magnitude relationship for a fault segment. Other applications of the technique include the estimation of the probability of background earthquakes, that is, earthquakes not associated with known or considered faults, and the estimation of the fraction of the total seismic moment associated with earthquakes less than the characteristic magnitude. Results for the San Francisco Bay region suggest that potentially damaging earthquakes with magnitudes less than the characteristic magnitudes should be expected. Comparisons of earthquake locations and the surface traces of active faults as determined from geologic data show significant disparities, indicating that a complete understanding of the relationship between earthquakes and faults remains elusive.

California↗

Estimating locations and magnitudes of earthquakes in eastern North America from Modified Mercalli intensities

We use 28 calibration events (3.7 < or = M < or = 7.3) from Texas to the Grand Banks, Newfoundland, to develop a Modified Mercalli intensity (MMI) model and associated site corrections for estimating source parameters of historical earthquakes in eastern North America. The model, MMI = 1.41 + 1.68 XM - 0.00345X Delta - 2.08log (Delta), where Delta is the distance in kilometers from the epicenter and M is moment magnitude, provides unbiased estimates of M and its uncertainty, and, if site corrections are used, of source location. The model can be used for the analysis of historical earthquakes with only a few MMI assignments. We use this model, MMI site corrections, and Bakun and Wentworth's (1997 technique to estimate M and the epicenter for three important historical earthquakes. The intensity magnitude M1 is 6.1 for the 18 November 1755 earthquake near Cape Ann, Massachusetts; 6.0 for the 5 January 1843 earthquake near Marked Tree, Arkansas; and 6.0 for the 31 October 1895 earthquake. The 1895 event probably occurred in southern Illinois, about 100 km north of the site of significant ground failure effects near Charleston, Missouri.

Bulletin of the Seismological Society of America↗

Empirical modified Mercalli intensity site corrections for towns in eastern North America

Modified Mercalli intensity (MMI) assignments for earthquakes in eastern North America (ENA) were used by Bakun et al. (2003) and Bakun and Hopper (in press) to develop models for estimating the location and moment magnitude M of earthquakes in ENA from MMI observations. The MMI empirical site corrections developed and used by Bakun et al. (2003) and Bakun and Hopper (in press) are listed in this Open-file Report.

Open-File Report↗

Differences in attenuation among the stable continental regions

There are systematic differences in the attenuation of damaging earthquake ground motions between different stable continental regions (SCRs). Seismic intensity and weak-motion data show that the attenuation in seismic waves for eastern North America (ENA) is less than for India, Africa, Australia, and northwest Europe. If ENA ground-motion attenuation relations are used in seismic hazard models for other SCRs, as is commonly done, then the estimated ground motions and resulting hazard may be too large. If an attenuation model that averages observations from ENA and the other SCRs is used to estimate the magnitudes of large historical earthquakes in ENA, as is the case for recent estimates of M for the 1811-1812 New Madrid, Missouri and the 1886 Charleston, South Carolina events, then the magnitude estimates for these events will be too large, as will be the resulting hazard.

Geophysical Research Letters↗

A bilinear source-scaling model for M-log a observations of continental earthquakes

The Wells and Coppersmith (1994) M-log A data set for continental earthquakes (where M is moment magnitude and A is fault area) and the regression lines derived from it are widely used in seismic hazard analysis for estimating M, given A. Their relations are well determined, whether for the full data set of all mechanism types or for the subset of strike-slip earthquakes. Because the coefficient of the log A term is essentially 1 in both their relations, they are equivalent to constant stress-drop scaling, at least for M ??? 7, where most of the data lie. For M > 7, however, both relations increasingly underestimate the observations with increasing M. This feature, at least for strike-slip earthquakes, is strongly suggestive of L-model scaling at large M. Using constant stress-drop scaling (???? = 26.7 bars) for M ??? 6.63 and L-model scaling (average fault slip u?? = ??L, where L is fault length and ?? = 2.19 &times 10-5) at larger M, we obtain the relations M = log A + 3.98 ?? 0.03, A ??? 537 km2 and M = 4/3 log A + 3.07 ?? 0.04, A > 537 km2. These prediction equations of our bilinear model fit the Wells and Coppersmith (1994) data set well in their respective ranges of validity, the transition magnitude corresponding to A = 537 km2 being M = 6.71.

Bulletin of the Seismological Society of America↗

The December 1872 Washington state earthquake

The largest historical earthquake in eastern Washington occurred on 15 December 1872. We used Modified Mercalli intensity (MMI) assignments for 12 twentieth-century earthquakes to determine attenuation relations for different regions in the Pacific Northwest. MMI attenuation for propagation paths east and west of the Cascade Mountains differs significantly only for epicentral distances greater than about 225 km. We used these attenuation relations and the MMI assignments for the 15 December 1872 earthquake to conclude that its epicentral region was east of the Cascade Mountains near Lake Chelan, Washington, and most probably near the south end of Lake Chelan. The intensity magnitude, M I , is 6.8 and moment magnitude, M , is 6.5-7.0 at the 95% confidence level .

Washington↗

Modified Mercalli Intensities (MMI) for some earthquakes in eastern North America (ENA) and empirical MMI site corrections for towns in ENA

Modified Mercalli Intensity (MMI) assignments for earthquakes in eastern North America (ENA) were used by Bakun et al. (submitted) to develop a model for eastern North America for estimating the location and moment magnitude M of earthquakes from MMI observations. MMI assignments for most of the earthquakes considered by Bakun et al. (submitted) are published. MMI assignments for 6 other earthquakes used by Bakun et al. (submitted) are listed in this report: November 18, 1755 near Cape Ann, Massachusetts; January 5, 1843 near Marked Tree, Arkansas; October 31, 1895 in southern Illinois; November 18, 1929 on the Grand Banks, Newfoundland; September 26, 1990 in southeast Missouri; and May 4, 1991 near Risco, Missouri. MMI empirical site corrections developed and used by Bakun et al. (submitted) are also listed in this report.

Arkansas;Illinois;Massachusetts;Missouri;Newfoundl↗

Modified Mercalli intensities (MMI) for large earthquakes near New Madrid, Missouri, in 1811-1812 and near Charleston, South Carolina, in 1886

Large historical earthquakes occurred in the eastern United States on December 16, 1811 near New Madrid, MO, on January 23, 1812 near New Madrid, MO, on February 7, 1812 near New Madrid, MO, and on September 1, 1886 near Charleston, SC. Modified Mercalli Intensity (MMI) assignments for these earthquakes were used by Bakun et al. (submitted) to estimate the location and moment magnitude M of these earthquakes from MMI observations. The MMI assignments used by Bakun et al. (submitted) are listed in this report.

Missouri, South Carolina↗

Comments on potential geologic and seismic hazards affecting Mare Island, Solano County, California

This report was prepared in response to a written request from the City of Vallejo, California, to the U.S. Geological Survey (USGS). By letter of October 4, 2002, the City requested that the USGS "provide advice to the City’s LNG Health and Safety Committee on its review of a potential liquid natural gas project" on the southern portion of Mare Island. The City specifically requested that the USGS advise the committee on potential hazards including fault rupture, earthquake ground motion, soil failure during earthquakes, tsunami and seiche, and landslides. The City requested that the USGS: (1) comment on these hazards, (2) describe its degree of confidence in its opinions, and (3) describe the scope of additional studies that will be needed if the City enters into an agreement with project sponsors. Advice was also requested on the selection of the safe shutdown and operating basis earthquakes as specified in the NFPA 59A standard (NFPA, 2001). This review of published reports and other publicly available information indicates that all of the hazards on which the USGS was asked to comment should be considered for the proposed project on the southern portion of Mare Island. Available information differs greatly for each of these potential hazards, and adequate understanding for design will require detailed site-specific investigations.

California↗

Seismicity of California's north coast

At least three moment magnitude (M) 7 earthquakes occurred along California's north coast in the second half of the nineteenth century. The M 7.3 earthquake on 23 November 1873 occurred near the California-Oregon coast and likely was located on the Cascadia subduction zone or within the North American plate. The M 7.0 + earthquake on 9 May 1878 was located about 75 km offshore near the Mendocino fault. The surface-wave magnitude (M(s)) 7.0 earthquake on 16 April 1899 was located about 150 km offshore within the Gorda plate. There were at least three M 7 north-coast earthquakes in the 35 years before 1906, two M 7 earthquakes in the 20 years after 1906, no M 7 earthquakes from 1923 until 1980, and four M 7 earthquakes since 1980. The relative seismic quiescence after 1906 for M 7 earthquakes along California's north coast mimics the post-1906 seismic quiescence in the San Francisco Bay area for M 6 earthquakes. The post-1906 relative quiescence did not extend to lower magnitudes in either area. The 18 April 1906 earthquake apparently influenced the rate of occurrence of M 7 north-coast earthquakes as it apparently influenced the rate of M 6 earthquakes in the San Francisco Bay area. The relative seismic quiescence along the California north-coast region after 1906 should be taken into account when evaluating seismic hazards in northwest California.

California↗

Seismic activity of the San Francisco Bay region

Moment magnitude M with objective confidence-level uncertainties are estimated for felt San Francisco Bay region earthquakes using Bakun and Went-worth's (1997) analysis strategy for seismic intensity observations. The frequency-magnitude distribution is well described for M ≧ 5.5 events since 1850 by a Gutenberg-Richter relation with a b -value of 0.90. The seismic moment rate Σ M 0 /yr since 1836 is 2.68 × 10 18 N-m/yr (95% confidence range = 1.29 × 10 18 N-m/yr to 4.07 × 10 18 N-m/yr); the seismic moment rate since 1850 is nearly the same. Σ M 0 /yr in the 56 years before 1906 is about 10 times that in the 70 years after 1906. In contrast, Σ M 0 /yr since 1977 is about equal that in the 56 years before 1906. 80% (1σ = 14%) of the plate-motion moment accumulation rate is available for release in earthquakes. The historical Σ M 0 /yr and the portion of the plate-motion moment accumulation rate available for release in earthquakes are used in a seismic cycle model to estimate the rate of seismic activity in the twenty-first century. High and low rates of future seismic activity are both permissible given the range of possible seismic-cycle recurrence times T and the uncertainties in the historical Σ M 0 and in the percentage of plate motion available for release in earthquakes. If the historical seismic moment rate is not greater than the estimated 2.68 × 10 18 N-m/yr and the percentage of the plate-motion moment accumulation available for release in earthquakes is not less than the estimated 80%, then for all T , the rate of seismic moment release from now until the next 1906-sized shock will be comparable to the rate from 1836 to 1905 when M 6 1/2 shocks occurred every 15 to 20 years.

California↗

Modified Mercalli intensities for some California north-coast earthquakes

Modified Mercalli Intensity (MMI) data are used by Bakun (in preparation) to estimate the location and moment magnitude M of California north-coast earthquakes from MMI observations only. The MMI data and site corrections used by Bakun (in preparation) are listed in this Open-file Report. This report is also accessible at http://quake.wr.usgs.gov/~bakun/.

Open-File Report↗

Corrections to "Estimating earthquake location and magnitude from seismic intensity data"

The confidence parameters in Table 5 and Figures 7 -12 in Bakun and Wentworth (1997) were calculated using the site corrections for MI(4) rather than Bakun and Wentworth's (1997) preferred relation MI(3). The conclusions of Bakun and Wentworth (1997) are not changed by these corrections. The corrected Table 5 and Figures are in this report and at http://quake.wr.usgs.gov/~bakun/.

Open-File Report↗

Review of Seismic Hazard Issues Associated with Auburn Dam Project, Sierra Nevada Foothills, California

Summary -- The U.S. Geological Survey was requested by the U.S. Department of the Interior to review the design values and the issue of reservoir-induced seismicity for a concrete gravity dam near the site of the previously-proposed Auburn Dam in the western foothills of the Sierra Nevada, central California. The dam is being planned as a flood-control-only dam with the possibility of conversion to a permanent water-storage facility. As a basis for planning studies the U.S. Army Corps of Engineers is using the same design values approved by the Secretary of the Interior in 1979 for the original Auburn Dam. These values were a maximum displacement of 9 inches on a fault intersecting the dam foundation, a maximum earthquake at the site of magnitude 6.5, a peak horizontal acceleration of 0.64 g, and a peak vertical acceleration of 0.39 g. In light of geological and seismological investigations conducted in the western Sierran foothills since 1979 and advances in the understanding of how earthquakes are caused and how faults behave, we have developed the following conclusions and recommendations: Maximum Displacement. Neither the pre-1979 nor the recent observations of faults in the Sierran foothills precisely define the maximum displacement per event on a fault intersecting the dam foundation. Available field data and our current understanding of surface faulting indicate a range of values for the maximum displacement. This may require the consideration of a design value larger than 9 inches. We recommend reevaluation of the design displacement using current seismic hazard methods that incorporate uncertainty into the estimate of this design value. Maximum Earthquake Magnitude. There are no data to indicate that a significant change is necessary in the use of an M 6.5 maximum earthquake to estimate design ground motions at the dam site. However, there is a basis for estimating a range of maximum magnitudes using recent field information and new statistical fault relations. We recommend reevaluating the maximum earthquake magnitude using current seismic hazard methodology. Design Ground Motions. A large number of strong-motion records have been acquired and significant advances in understanding of ground motion have been achieved since the original evaluations. The design value for peak horizontal acceleration (0.64 g) is larger than the median of one recent study and smaller than the median value of another. The value for peak vertical acceleration (0.39 g) is somewhat smaller than median values of two recent studies. We recommend a reevaluation of the design ground motions that takes into account new ground motion data with particular attention to rock sites at small source distances. Reservoir-Induced Seismicity. The potential for reservoir-induced seismicity must be considered for the Auburn Darn project. A reservoir-induced earthquake is not expected to be larger than the maximum naturally occurring earthquake. However, the probability of an earthquake may be enhanced by reservoir impoundment. A flood-control-only project may involve a lower probability of significant induced seismicity than a multipurpose water-storage dam. There is a need to better understand and quantify the likelihood of this hazard. A methodology should be developed to quantify the potential for reservoir induced seismicity using seismicity data from the Sierran foothills, new worldwide observations of induced and triggered seismicity, and current understanding of the earthquake process. Reevaluation of Design Parameters. The reevaluation of the maximum displacement, maximum magnitude earthquake, and design ground motions can be made using available field observations from the Sierran foothills, updated statistical relations for faulting and ground motions, and current computational seismic hazard methodologies that incorporate uncertainty into the analysis. The reevaluation does not require significant new geological field studies.

Open-File Report↗

Early warning system for aftershocks

A prototype early warning system to provide San Francisco and Oakland, California a few tens-of-seconds warning of incoming strong ground shaking from already-occurred M ≧ 3.7 aftershocks of the magnitude 7.1 17 October 1989 Loma Prieta earthquake was operational on 28 October 1989. The prototype system consisted of four components: ground motion sensors in the epicentral area, a central receiver, a radio repeater, and radio receivers. One of the radio receivers was deployed at the California Department of Transportation ( CALTRANS ) headquarters at the damaged Cypress Street section of the I-880 freeway in Oakland, California on 28 October 1989 and provided about 20 sec of warning before shaking from the M 4.5 Loma Prieta aftershock that occurred on 2 November 1989 at 0550 UTC. In its first 6 months of operation, the system generated triggers for all 12 M > 3.7 aftershocks for which trigger documentation is preserved, did not trigger on any M ≦ 3.6 aftershocks, and produced one false trigger as a result of a now-corrected single point of failure design flaw. Because the prototype system demonstrated that potentially useful warnings of strong shaking from aftershocks are feasible, the USGS has completed a portable early warning system for aftershocks that can be deployed anywhere.

California↗