Chris H. Cramer

A Time-Dependent Probabilistic Seismic-Hazard Model for California

For the purpose of sensitivity testing and illuminating nonconsensus components of time-dependent models, the California Department of Conservation, Division of Mines and Geology (CDMG) has assembled a time-dependent version of its statewide probabilistic seismic hazard (PSH) model for California. The model incorporates available consensus information from within the earth-science community, except for a few faults or fault segments where consensus information is not available. For these latter faults, published information has been incorporated into the model. As in the 1996 CDMG/U.S. Geological Survey (USGS) model, the time-dependent models incorporate three multisegment ruptures: a 1906, an 1857, and a southern San Andreas earthquake. Sensitivity tests are presented to show the effect on hazard and expected damage estimates of (1) intrinsic (aleatory) sigma, (2) multisegment (cascade) vs. independent segment (no cascade) ruptures, and (3) time-dependence vs. time-independence. Results indicate that (1) differences in hazard and expected damage estimates between time-dependent and independent models increase with decreasing intrinsic sigma, (2) differences in hazard and expected damage estimates between full cascading and not cascading are insensitive to intrinsic sigma, (3) differences in hazard increase with increasing return period (decreasing probability of occurrence), and (4) differences in moment-rate budgets increase with decreasing intrinsic sigma and with the degree of cascading, but are within the expected uncertainty in PSH time-dependent modeling and do not always significantly affect hazard and expected damage estimates.

2001 Bhuj, India, Earthquake Engineering Seismoscope Recordings and Eastern North America Ground-Motion Attenuation Relations

Engineering seismoscope data collected at distances less than 300 km for the M 7.7 Bhuj, India, mainshock are compatible with ground-motion attenuation in eastern North America (ENA). The mainshock ground-motion data have been corrected to a common geological site condition using the factors of Joyner and Boore (2000) and a classification scheme of Quaternary or Tertiary sediments or rock. We then compare these data to ENA ground-motion attenuation relations. Despite uncertainties in recording method, geological site corrections, common tectonic setting, and the amount of regional seismic attenuation, the corrected Bhuj dataset agrees with the collective predictions by ENA ground-motion attenuation relations within a factor of 2. This level of agreement is within the dataset uncertainties and the normal variance for recorded earthquake ground motions.

Quantifying the Uncertainty in Site Amplification Modeling and Its Effects on Site-Specific Seismic-Hazard Estimation in the Upper Mississippi Embayment and Adjacent Areas

The Mississippi embayment, located in the central United States, and its thick deposits of sediments (over 1 km in places) have a large effect on earthquake ground motions. Several previous studies have addressed how these thick sediments might modify probabilistic seismic-hazard maps. The high seismic hazard associated with the New Madrid seismic zone makes it particularly important to quantify the uncertainty in modeling site amplification to better represent earthquake hazard in seismic-hazard maps. The methodology of the Memphis urban seismic-hazard- mapping project (Cramer et al. , 2004) is combined with the reference profile approach of Toro and Silva (2001) to better estimate seismic hazard in the Mississippi embayment. Improvements over previous approaches include using the 2002 national seismic-hazard model, fully probabilistic hazard calculations, calibration of site amplification with improved nonlinear soil-response estimates, and estimates of uncertainty. Comparisons are made with the results of several previous studies, and estimates of uncertainty inherent in site-amplification modeling for the upper Mississippi embayment are developed. I present new seismic-hazard maps for the upper Mississippi embayment with the effects of site geology incorporating these uncertainties.

Average QLg, QSn, and observation of Lg blockage in the Continental Margin of Nova Scotia

The term “Lg blockage” refers to the sudden disappearance of the Lg phase along a particular propagation path which is commonly seen at continental-oceanic transition zones. In this paper we present observational evidence of Lg blockage across the continental margin of Nova Scotia in eastern Canada. Regional Lg and Sn spectra from 91 events with epicentral distances between 100 and 1200 km and magnitudes between 2.5 and 4.7 are inverted simultaneously for the source spectrum, site amplification, and average attenuation. The vertical displacement spectra were estimated between 0.9 and 10.75 Hz. The assumptions include a fixed frequency-independent geometric spreading rate for Lg and a frequency-dependent spreading model for the Sn. Estimates for the apparent regional attenuations are QLg (f) = 615(±25) f0.35(±0.04) and QSn (f) = 404(±23) f0.45(±0.03). Results from this study provide an accurate parameterization of observed amplitude spectra and are valuable for representing wave propagation in the region. Based on the observation of a strong trade-off between Sn and Lg amplitudes which have different attenuation characteristics, we conclude any attenuation study based on measuring amplitude of a package of several different phases, without taking into consideration the propagation characteristics of individual waveforms at the region of study, may bias the estimation of average regional Q.

Discrepancy between earthquake rates implied by historic earthquakes and a consensus geologic source model for California

We examine the difference between expected earthquake rates inferred from the historical earthquake catalog and the geologic data that was used to develop the consensus seismic source characterization for the state of California [California Department of Conservation, Division of Mines and Geology (CDMG) and U.S. Geological Survey (USGS) Petersen et al., 1996; Frankel et al., 1996]. On average the historic earthquake catalog and the seismic source model both indicate about one M 6 or greater earthquake per year in the state of California. However, the overall earthquake rates of earthquakes with magnitudes ( M ) between 6 and 7 in this seismic source model are higher, by at least a factor of 2, than the mean historic earthquake rates for both southern and northern California. The earthquake rate discrepancy results from a seismic source model that includes earthquakes with characteristic (maximum) magnitudes that are primarily between M 6.4 and 7.1. Many of these faults are interpreted to accommodate high strain rates from geologic and geodetic data but have not ruptured in large earthquakes during historic time. Our sensitivity study indicates that the rate differences between magnitudes 6 and 7 can be reduced by adjusting the magnitude-frequency distribution of the source model to reflect more characteristic behavior, by decreasing the moment rate available for seismogenic slip along faults, by increasing the maximum magnitude of the earthquake on a fault, or by decreasing the maximum magnitude of the background seismicity. However, no single parameter can be adjusted, consistent with scientific consensus, to eliminate the earthquake rate discrepancy. Applying a combination of these parametric adjustments yields an alternative earthquake source model that is more compatible with the historic data. The 475-year return period hazard for peak ground and 1-sec spectral acceleration resulting from this alternative source model differs from the hazard resulting from the standard CDMG–USGS model by less than 10% across most of California but is higher (generally about 10% to 30%) within 20 km from some faults.

Estimating Earthquake Magnitudes from Reported Intensities in the Central and Eastern United States

A new macroseismic intensity prediction equation is derived for the central and eastern United States and is used to estimate the magnitudes of the 1811–1812 New Madrid, Missouri, and 1886 Charleston, South Carolina, earthquakes. This work improves upon previous derivations of intensity prediction equations by including additional intensity data, correcting magnitudes in the intensity datasets to moment magnitude, and accounting for the spatial and temporal population distributions. The new relation leads to moment magnitude estimates for the New Madrid earthquakes that are toward the lower range of previous studies. Depending on the intensity dataset to which the new macroseismic intensity prediction equation is applied, mean estimates for the 16 December 1811, 23 January 1812, and 7 February 1812 mainshocks, and 16 December 1811 dawn aftershock range from 6.9 to 7.1, 6.8 to 7.1, 7.3 to 7.6, and 6.3 to 6.5, respectively. One‐sigma uncertainties on any given estimate could be as high as 0.3–0.4 magnitude units. We also estimate a magnitude of 6.9±0.3 for the 1886 Charleston, South Carolina, earthquake. We find a greater range of magnitude estimates when also accounting for multiple macroseismic intensity prediction equations. The inability to accurately and precisely ascertain magnitude from intensities increases the uncertainty of the central United States earthquake hazard by nearly a factor of two. Relative to the 2008 national seismic hazard maps, our range of possible 1811–1812 New Madrid earthquake magnitudes increases the coefficient of variation of seismic hazard estimates for Memphis, Tennessee, by 35%–42% for ground motions expected to be exceeded with a 2% probability in 50 years and by 27%–35% for ground motions expected to be exceeded with a 10% probability in 50 years.

Active fault near-source zones within and bordering the State of California for the 1997 Uniform Building Code

The fault sources in the Project 97 probabilistic seismic hazard maps for the state of California were used to construct maps for defining near-source seismic coefficients, N a and N v , incorporated in the 1997 Uniform Building Code (ICBO 1997). The near-source factors are based on the distance from a known active fault that is classified as either Type A or Type B. To determine the near-source factor, four pieces of geologic information are required: (1) recognizing a fault and determining whether or not the fault has been active during the Holocene, (2) identifying the location of the fault at or beneath the ground surface, (3) estimating the slip rate of the fault, and (4) estimating the maximum earthquake magnitude for each fault segment. This paper describes the information used to produce the fault classifications and distances.

An Assessment of the Impact of the 2003 epri Ground-Motion Prediction Models on the usgs National Seismic-Hazard Maps

Ground-motion attenuation relations have an important impact on seismic hazard analyses. Ground-motion modeling is particularly sensitive to assumptions about wave-propagation attenuation (crustal Q and geometrical spreading), as well as source and site conditions. Studies of path attenuation from earthquakes in eastern North America (ena) provide insights into the appropriateness of specific attenuation relations. An Electric Power Research Institute (EPRI) (2003, 2004) study combines published ena ground-motion attenuation relations into four model forms: single-corner, double-corner, hybrid-empirical, and finite-fault. When substituted in the U.S. Geological Survey 2002 national seismic hazard maps for the five ena relations originally used in those hazard calculations, the EPRI (2003) relations predict similar ground motions and hazard at short periods ( 0.5 sec), relative to the 2002 national maps. A major reason for this difference is due to the crustal seismic-wave attenuation model assumed in a few of the ena relations combined into the EPRI (2003, 2004) models. Although appropriate differences in geometrical spreading models among ena relations can also be significant, a few ena relations have 1-Hz Q -values ( Q ) that are below the EPRI (1993) consensus range for Q when coupled with a geometrical spreading of R −0.5. The EPRI (2003, 2004) single-corner relation is strongly influenced by the inclusion of ena relations with assumed Q below the EPRI (1993) range, which explains much of the discrepancy in predictions at longer periods.

Seismic hazard estimation of northern Iran using smoothed seismicity

This article presents a seismic hazard assessment for northern Iran, where a smoothed seismicity approach has been used in combination with an updated seismic catalog and a ground motion prediction equation recently found to yield good fit with data. We evaluate the hazard over a geographical area including the seismic zones of Azerbaijan, the Alborz Mountain Range, and Kopeh-Dagh, as well as parts of other neighboring seismic zones that fall within our region of interest. In the chosen approach, seismic events are not assigned to specific faults but assumed to be potential seismogenic sources distributed within regular grid cells. After performing the corresponding magnitude conversions, we decluster both historical and instrumental seismicity catalogs to obtain earthquake rates based on the number of events within each cell, and smooth the results to account for the uncertainty in the spatial distribution of future earthquakes. Seismicity parameters are computed for each seismic zone separately, and for the entire region of interest as a single uniform seismotectonic region. In the analysis, we consider uncertainties in the ground motion prediction equation, the seismicity parameters, and combine the resulting models using a logic tree. The results are presented in terms of expected peak ground acceleration (PGA) maps and hazard curves at selected locations, considering exceedance probabilities of 2 and 10% in 50 years for rock site conditions. According to our results, the highest levels of hazard are observed west of the North Tabriz and east of the North Alborz faults, where expected PGA values are between about 0.5 and 1 g for 10 and 2% probability of exceedance in 50 years, respectively. We analyze our results in light of similar estimates available in the literature and offer our perspective on the differences observed. We find our results to be helpful in understanding seismic hazard for northern Iran, but recognize that additional efforts are necessary to obtain more robust estimates at specific areas of interest and different site conditions.

Viscosity of the Atlantic Ocean Bottom

Two profiles across the Mid-Atlantic Ridge were analyzed to determine viscosity values. Viscous creeping of bottom features due to gravitational stress was assumed, as was sea floor spreading at a rate of 2 centimeters per year. Values obtained agreed well with previous results obtained on the Fennoscandian Uplift, despite great differences in the horizontal dimensions of the bottom relief.