Quick
Search: 
 
advanced search
 GSW Home    GeoRef Home    My GSW Alerts    Contact GSW    About GSW    Journals List    Help 
Seismological  Research Letters Signup for GSW Email News
JOURNAL HOME HELP CONTACT PUBLISHER SUBSCRIBE ARCHIVE SEARCH TABLE OF CONTENTS

Seismological Research Letters; November 2007; v. 78; no. 6; p. 649-662; DOI: 10.1785/gssrl.78.6.649
© 2007 Seismological Society of America
This Article
Right arrow Extract
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by O'Connell, D. R. H.
Right arrow Articles by Liu, P.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation

Probabilistic Ground-motion Assessment of Balanced Rocks in the Mojave Desert

Daniel R. H. O'Connell and Roland LaForge
William Lettis and Associates, Inc.

Pengcheng Liu
Bureau of Reclamation


    INTRODUCTION
 TOP
 INTRODUCTION
 EFFECT OF CENSORING...
 SIMPLIFIED PSHA ANALYSIS
 MONTE CARLO SIMULATIONS OF...
 3D SITE EFFECTS
 VARIATION OF PEAK SPECTRAL...
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
Brune (1996, 1999) proposed using balanced rocks located at a variety of distances from faults, some with relatively high slip rates (e.g., the San Andreas fault), as proxy paleo-ground-motion indicators of peak acceleration. Shi et al. (1996) derived simplified relationships between peak horizontal accelerations (PHA) and toppling probabilities to convert pseudostatic field measurements of rock pedestal stability to limits on PHA consistent with rock stability and persistence. Anderson and Brune (1999) concluded that the persistence of balanced rocks 10-30 km from the San Andreas fault was inconsistent with the ergodic assumption implicit in probabilistic seismic hazard analyses (PSHA) derived using the approach of Cornell (1968). The existence of a modest population of balanced rocks on the Mojave side of the San Andreas fault affords an opportunity to test the hypothesis of Anderson and Brune (1999) that the persistence of these balanced rocks is inconsistent with site-specific PSHA and the ergodic assumption. Through Monte Carlo simulations of balanced-rock life cycles we show that the existing empirical ground-motion relation that explicitly accounts for site velocity (Boore et al. 1997) is consistent with PSHA assumptions (Cornell 1968) and the persistence of balanced rocks east of the San Andreas fault in the Mojave Desert (Brune 1996, 1999).

In this paper, discussion is not limited to balanced rocks, because a wide variety of balanced objects can provide useful information on ground motions. For example, various types of monuments have the advantage that the duration that the monument has remained standing is documented in the historical record, and the fragility of the monument can be well-established. The duration of balanced-object persistence is a necessary quantity to make statistical inferences about the rate of exceedance or nonexceedance of ground-motion amplitudes.

The relative paucity of well-constrained age estimates for balanced rocks is a significant limitation on the statistical inferences that can be derived from balanced rocks on ground-motion amplitude rates. However, Bell et al. (1998) provide sufficient age-dating constraints to make the Mojave Desert balanced-rock sites a test bed for statistical testing.

Two classes of information must be estimated from balanced rocks to allow direct comparisons with PSHA results, ground-motion rate information, and ground-motion amplitude information. Fragility relationships between toppling and ground-motion parameters provide a means to estimate ground-motion amplitude information. Housner (1963) indicated that without estimates of the prior population of balanced objects, it's not possible to infer ground-motion amplitude statistics solely from a surviving population of such objects after the occurrence of an earthquake. However, this doesn't mean that statistical information about ground-motion characteristics cannot be obtained from a surviving population of balanced objects. To do so, we need to employ statistical approaches that explicitly account for censoring. In statistics, censoring occurs when the value of an observation is only partially known.

Balanced rocks are represented by two types of censoring. The population of surviving balanced rocks represents a left censoring of ground-motion amplitudes. The surviving balanced rocks are left-censored because they experienced ground-motion accelerations, velocities, and durations, etc. less than the thresholds required to cause toppling, but we don't know the amplitudes of the differences between the critical values required for toppling and the actual values of the ground-motion parameters that the surviving balanced rocks were subjected to. Ground-motion rate information is censored because the total population of balanced rocks through time is unknown. While we are completely ignorant of the total population of balanced rocks prior to an earthquake, we know that the number of observed surviving balanced rocks lies within the interval of 0-100% of the pre-earthquake population over some period of time. This is a form of interval censoring, but clearly indicates we need supplementary balanced-rock rate information to constrain ground-motion rates.

The duration that a balanced object exists with a quantifiable fragility behavior provides constraints on ground-motion rates. Since balanced rocks are the product of differential erosion, their fragility behavior evolves with time toward states of increased fragility and eventual static collapse. In this paper we first consider a simplified static fragility curve, then we evaluate cases of continuous increases of fragility with age. Specific estimates of durations that balanced rocks have existed with particular fragility characteristics are not available from the Mojave balanced-rock sites, although it is generally assumed that they have persisted for at least 10,000-20,000 years (Bell et al. 1998; Anooshehpoor, personal communication). The combination of amplitude constraints (fragility curves), time constraints (time period that balanced objects have not toppled), and the rate that balanced rocks are created are required to extract information about the nonexceedance of ground-motion amplitudes as a function of frequency or annual exceedance probability (AEP). Without balanced-rock persistence-duration information and balanced-rock genesis-rate information, there is not sufficient information to compare against nonexceedance rate estimates from PSHAs for balanced-rock sites.

Knowledge of the population of toppled objects would provide a second type of censored information (binomial censoring) that represents positive exceedance information. As discussed in O'Connell (2005), the combination of exceedance and nonexceedance information can be sufficient in itself to provide statistically useful constraints on amplitude frequency (flood frequency in the specific examples in O'Connell 2005). Again a time-duration context is needed to convert exceedances into a useful binomial-censored dataset. For instance, if n toppled objects with fragility information could be associated with toppling occurring over a time interval T, this information can be used with appropriate uncertainties in n and T to estimate peak amplitude frequency (O'Connell 2005). Toro and Cornell (2006) have begun exploring application of the approach in O'Connell et al. (2002) to the statistical analysis of balanced rocks. However, in this paper we take a different approach to illustrate some important statistical concepts and use Monte Carlo simulations to investigate the sensitivity of statistical inferences on various balanced-rock parameters, particularly balanced-rock survival duration and simplified fragility behavior.

It has been clearly established that the toppling of balanced objects often involves nonlinear dynamics (Yim et al. 1980a,b; Zhang and Makris 2000; Makris and Zhang 2001). For example, Yim et al. (1980a) found that small changes in vertical motions could cause objects not to topple at higher horizontal accelerations, although the same objects toppled at lower horizontal accelerations. Further, balanced-object toppling can be a function of ground-motion parameters other than PHA, including duration, spectral accelerations, peak velocities, etc. However, for the purposes of the statistical formulation focus of this paper we only consider toppling as a function of PHA. Even for PHA, balanced-object fragility curves can be complex and have multiple modes (Yim et al. 1980a,b; Zhang and Makris 2000; Makris and Zhang 2001), but even if fragility curves are complex, their shapes can be explicitly accounted for in Bayesian parametric and nonparametric frequency estimation procedures (O'Connell et al. 2002; O'Connell 2005) or Monte Carlo simulations. However, this paper does not dwell on the details of balanced-rock fragility dynamics but rather focuses on statistical inference approaches that are general enough to incorporate any balanced-object fragility behavior that is appropriate.

We devise statistical tests that are relatively insensitive to or explicitly account for the multimodal, nonlinear nature of PHA-toppling relationships, to test if the persistence of the balanced rocks in Brune (1996, 1999) is consistent with site-specific PSHA and the ergodic assumption. We conduct a site-specific PSHA for the Mojave sites using site-specific shear-wave velocity estimates from Abbott et al. (2001) to construct 30-m-average shear-wave velocities (Vs30) for the ground-motion prediction relation of Boore et al. (1997) that explicitly accounts for Vs30, instead of using the lumped stiff-soil/rock and soft-soil classifications implemented in Abrahamson and Silva (1997) and Sadigh et al. (1997). While we show that mean PHA hazard curves are consistent with the existence of balanced rocks in Brune (1996, 1999), it's not clear if the balanced rocks represent mean, modal, or some other statistic correlated with PHA. Consequently, we delve into the possible statistical relationships between balanced rocks and PHA. We evaluate the impact of censoring (survival versus toppling) on ground-motion statistics. Then, assuming that the current existence of the balanced rocks implies persistence of the balanced rocks over the Quaternary, we show that persistence of balanced rocks in their current locations is consistent with site-specific ground motion associated with repeated characteristic earthquakes on the San Andreas fault, a result that directly contradicts the conclusions of Anderson and Brune (1999).


    EFFECT OF CENSORING (SURVIVAL/TOPPLING) ON PHA STATISTICS
 TOP
 INTRODUCTION
 EFFECT OF CENSORING...
 SIMPLIFIED PSHA ANALYSIS
 MONTE CARLO SIMULATIONS OF...
 3D SITE EFFECTS
 VARIATION OF PEAK SPECTRAL...
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
Each class of balanced rock (toppled or surviving) has its own population of ground-motion amplitudes with means and variances that are different from the means and variances associated with the total population of input ground-motion amplitudes (PHA is used to represent ground-motion amplitude for simplicity here). We construct several simple simulations to illustrate the resulting biases in the estimates of the mean PHA that result from using sample means and variances of either the surviving or toppled population of rocks. For the purposes of illustration we use a natural log PHA standard deviation, ln({sigma}), of 0.55, consistent with empirical ground-motion relations (Abrahamson and Silva 1997; Sadigh et al. 1997, and Boore et al. 1997) for M ~ 6, but larger than empirical ln({sigma}) for larger (M > 7) earthquakes. We use the ground-motion relation of Boore et al. (1997) that allows specification of Vs30 to illustrate the effects of site velocity on mean PHA for sites with Vs30 of 760 m/s, 1,250 m/s, and 2,250 m/s. The scenario strike-slip earthquake magnitude is M 7.8 and the site is located 20 km from the closest portion of the surface trace of the fault (figure 1).

We construct a balanced-rock fragility curve using a mean toppling PHA of 0.25 and a natural-log-PHA-toppling standard deviation of 0.55 (figure 2), and sample from it using a uniform random number generator with 10,000 trials to obtain simulations of outcomes for the three Vs30 scenarios in table 1. Specifically, in each trial we generate a random PHA and then generate another uniform random number between 0 and 1 to sample from figure 2 to determine if the rock topples or survives. As expected, the standard deviations of surviving and toppling rocks are smaller than the actual PHA standard deviation because the fragility curve splits (censors) the populations into two distinct groups. The surviving mean PHAs are systematically smaller than the actual ground motion PHAs because the fragility curve censors most of the upper tail of the ground motion-distribution for surviving rocks. Similarly, the mean toppling PHAs are larger than the mean ground-motion PHAs because the fragility curve censors most of the lower tail of the ground-motion distribution. A greater proportion of rocks survive as Vs30 increases because PHA decreases with increasing Vs30. Thus, the assumption that balanced-rock mean PHA corresponds to actual ground-motion PHA results in significant biases of 36% for a nominal National Earthquake Hazards Reduction Program (NEHRP) B site condition (Vs30 of 760 m/s), 29% for a Vs30 of 1,250 m/s, and 22% for a Vs30 of 2,250 m/s.


View this table:
[in this window]
[in a new window]

 
TABLE 1 PHA Censoring Results

 


Figure 1
View larger version (77K):
[in this window]
[in a new window]

 
Figure 1. Plan view of shaded topography, the San Andreas fault (white line), and balanced rock sites in the Mojave Desert from Brune (1999) (white circles).

 


Figure 2
View larger version (21K):
[in this window]
[in a new window]

 
Figure 2. Balanced object fragility distribution used for the censoring simulations.

 
Table 1 shows that for a known population of balanced objects observed at many distinct independent sites (sites located at sufficient distances from each other so that ground motions become uncorrelated at the periods of interest), the mean PHA associated with the surviving balanced objects underestimates the actual mean ground-motion PHA by > 20%, with the bias increasing with increasing mean PHA. Restricting balanced-object analyses to a single site does not provide multiple independent observations from a single earthquake, but multiple highly correlated observations. It's necessary to obtain observations from sites sufficiently separated that ground-motion coherence is negligible (statistical independence) to make inferences about ground-motion variability from a single earthquake. A collection of spatially independent sites that have been exposed to repeated earthquakes provides the spatial independence necessary to evaluate the ergodic assumption from repeated earthquakes. Thus, the Mojave balanced rock sites east of the San Andreas fault (Brune 1996, 1999) provide a critical mass (eight independent sites) of balanced rocks to conduct statistical investigations of probabilistic ground-motion predictions.


Figure 3
View larger version (14K):
[in this window]
[in a new window]

 
Figure 3. Balanced rock shear-velocity profiles modified from Abbott et al. (2001) for the Mojave Desert sites. Velocities are 30-m depth averages. Note the balanced rocks themselves are not located on 5-12 m thicknesses of grus; the profiles represent site velocities adjacent to the rock pedestals.

 

    SIMPLIFIED PSHA ANALYSIS
 TOP
 INTRODUCTION
 EFFECT OF CENSORING...
 SIMPLIFIED PSHA ANALYSIS
 MONTE CARLO SIMULATIONS OF...
 3D SITE EFFECTS
 VARIATION OF PEAK SPECTRAL...
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
We calculate a 1%-in-a-100-year PSHA analysis using the approach of Cornell (1968) that accounts for the Mojave site Vs30 estimates from Abbott et al. (2001) to compare to Brune's (1999) PSHA with his rock toppling PHAs estimated for the Mojave sites (figure 1). Abbott et al. (2001) showed that Mojave balanced rock sites (figure 3) all have Vs30 values that place the sites in NEHRP class B (Vs30 > 760 m/s). Since the rock pedestals are located outside the shallow low-velocity grus regions where the velocities were estimated in Abbott et al. (2001), as indicated in the photographs in Brune (1999), a Vs30 = 1,250 m/s, corresponding to the lower sub-grus velocity of the two site-velocity profiles in Abbott et al. (2001) was used in the PSHA analyses to represent typical site conditions beneath the rock pedestals. Rodriguez-Marek et al. (2001) found that competent rock sites had lower PHA dispersion (ln({sigma}) = 0.4) relative to stiff soil sites (ln({sigma}) = 0.6). However, for the purposes of illustration the standard ln({sigma}) of 0.55 from Boore et al. (1997) is used to calculate PHA hazard.

Following Cao et al. (2003) and Petersen et al. (1996), the Mojave segment of the San Andreas fault was assigned slip-rate scenarios of 23, 30, and 37 mm/yr and a characteristic magnitude of M 7.4. Earthquakes within the Mojave Desert region shown in figure 4 were used to establish background earthquake recurrence for a maximum background earthquake magnitude of 6.5 (figure 5). The existence of the balanced rocks is consistent with repeated M 7.4 earthquakes on the Mojave segment of the San Andreas fault and M 5.5+ background earthquakes in the region, particularly when the censoring bias is used to adjust the toppling accelerations (figure 6). As noted by Anooshehpoor et al. (2004), dynamic toppling accelerations are typically 30% higher than quasistatic toppling accelerations. Consequently, the combined effects of censoring and toppling dynamics further shift the balanced-rock toppling PHA for comparison to the PSHA results as indicated in figure 6. This implies that these balanced rocks persisted over about 48 characteristic M 7.4 1857-type earthquakes on the San Andreas fault for a ~ 10,000-year exposure period, consistent with the inferred age of the rocks (Bell et al. 1998).

The naïve conclusion is that the agreement between site-specific PSHA and balanced-rock toppling-acceleration estimates supports using the ergodic assumption in PSHA. Figure 6 shows there is no inconsistency, as indicated by Brune (1999) and Anderson and Brune (1999), between balanced rocks and PSHA methodology, once actual site conditions are used in a PSHA. However, a realistic comparison of balanced-rock constraints and PSHA predictions requires explicit consideration of the residence time and life-cycle characteristics of balanced rocks.


Figure 4
View larger version (160K):
[in this window]
[in a new window]

 
Figure 4. Shaded topography with the triangle enclosing the region of the Mojave Desert used to calculate background earthquake recurrence. Red circles are epicenters, with the smallest circle corresponding to M 3.5-4, and the largest circle to M 5.0-5.5. Seismicity from NCEDC (U.C. Berkeley)1930-April 2004. Constant light gray regions are areas below mean sea level.

 


Figure 5
View larger version (19K):
[in this window]
[in a new window]

 
Figure 5. Incremental recurrence estimates for the Mojave Desert background zone derived from the data in Figure 4.

 


Figure 6
View larger version (25K):
[in this window]
[in a new window]

 
Figure 6. Mean site-specific 1% in 100 years PHA results using Vs30=1250 m/s and only the San Andreas fault contribution are the diamonds connected by the solid line. The stars are the original quasi-static PHA estimates for the balanced rocks from Brune (1999), the crosses indicate the 29% increase to account for censoring bias, and the boxes indicate the estimates when censoring biases and dynamic toppling PHA biases from Anoonshehpoor et al. (2004) are included. The labeled dotted and dashed curves show the total PHA hazard from the San Andreas fault and background seismicity as smaller magnitudes are included in the background PHA hazard calculations.

 

    MONTE CARLO SIMULATIONS OF BALANCED-ROCK LIFE CYCLES
 TOP
 INTRODUCTION
 EFFECT OF CENSORING...
 SIMPLIFIED PSHA ANALYSIS
 MONTE CARLO SIMULATIONS OF...
 3D SITE EFFECTS
 VARIATION OF PEAK SPECTRAL...
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
We assume that the balanced rocks in the Mojave Desert have persisted over many cycles of large San Andreas fault earthquakes, based on the estimated ages of the balanced rocks (Bell et al. 1998) and the relatively short return period of 210 years for M 7.4 Mojave-segment San Andreas earthquakes (Petersen et al. 1996; Cao et al. 2003). The Monte Carlo simulation approach models the life cycles of a collection of balanced rocks, starting with an initial population of 20 balanced rocks, by subjecting the balanced-rock population to 4,000 cycles of San Andreas Mojave-segment M 7.4 earthquakes, and finds the combinations of ground motion and balanced-rock life-cycle parameters that allow a nonzero number of rocks to persist through all 4,000 earthquake cycles. To do this we must specify the lifetime of balanced rocks (rate they are destroyed due to static instability associated with erosive processes), their time-dependent fragility behavior, and the rate at which new balanced rocks are created.

The Mojave Desert balanced rocks are a product of fracture-controlled differential erosion of granitic boulders (Bell et al. 1998). Even in the absence of earthquake loadings, it's necessary to continually create new balanced rocks through differential erosion to replace rocks that have toppled as they reach their static stability limits. For instance, to maintain a population of about 20 balanced rocks that have a fragility of 20,000 years, it's necessary to create a new balanced rock about every 1,000 years as older rocks are retired through erosion to their static stability limit. So to use the 44 balanced rocks from the Mojave Desert in Brune (1999) in a probabilistic analysis it's necessary to postulate the average steady-state rate that new balanced rocks are created. The existence of the balanced rocks combined with earthquake loadings places a lower bound on the rate that new balanced rocks are created and is investigated numerically in the rock life-cycle simulations.


Figure 7
View larger version (20K):
[in this window]
[in a new window]

 
Figure 7. Three scenario time-dependent balanced rock fragility curves. Solid line is linear, long-dashed curve is exponential, and short-dashed curve is double exponential. The double exponential curve spends nearly 80% of its lifetime at < 0.5 g and 54% of its lifetime at < 0.25 g (thin vertical dotted lines). In contrast, the exponential curve spends about 62% of its lifetime at < 0.5 g and only 35% of its lifetime at < 0.25 g (thick vertical dotted lines).

 
We simulate the persistence or extinguishing of a population of balanced rocks over 800,000-1,200,000-year durations of 4,000 Mojave-segment M 7.4 San Andreas earthquakes at a distance of 20 km from the San Andreas fault by varying ground-motion and balanced-rock parameters. The balanced-rock and earthquake parameters considered are the erosion rate that creates new balanced rocks in terms of an effective birth interval in years and the time interval in years over which rocks degrade statically toward collapse, as listed in table 2, and the three time-dependent fragility-curve scenarios in figure 7.


View this table:
[in this window]
[in a new window]

 
TABLE 2 Balanced-Rock Life-Cycle Parameter Values

 
We calculate the most appropriate site-specific estimate of mean PHA for a M 7.4 strike-slip earthquake at 20 km using Boore et al. (1997) with a site Vs30=1250 m/s to establish the minimum mean PHA used in the simulations (table 3). The earthquake parameters include the mean San Andreas Mojave-segment M 7.4 return period, mean PHA, and mean PHA ln({sigma}), as listed in table 3. The simulations vary all the balanced-rock and earthquake parameters over the ranges in tables 2 and 3. Balanced rocks are statically created and destroyed based on the parameters in table 2. For each earthquake, PHAs are randomly generated consistent with the ground-motion parameters in table 3, and the balanced rocks either topple or survive depending on whether the PHAs exceed the fragility PHA limits in figure 7 as a function of time within the balanced rocks' life cycles. If during any of the 4,000 earthquake cycles no balanced rocks survive, the simulation is terminated and flagged as a nonpersistent balanced-rock outcome. Only simulation parameter combinations that maintained a nonzero number of surviving balanced rocks through all 4,000 earthquake cycles were flagged as persistent.


View this table:
[in this window]
[in a new window]

 
TABLE 3 Ground-Motion Parameter Ranges

 
The simulation results depend strongly on the shape of the time-dependent fragility curve and birth interval (rate) that new balanced rocks are created. A linear fragility curve (table 4) allows the largest range of ground-motion and balanced-rock parameters to be consistent with the persistence of balanced rocks in the Mojave Desert. Linear fragility allows the largest mean PHA, largest ln({sigma}), longest rock lifetimes (20,000 years), and slowest rock birth rates (two new balanced rocks created over 20,000 years) consistent with the existence of the balanced rocks in the Mojave Desert (Brune 1996, 1999). As fragility develops more rapidly with age (table 5), the largest PHA (0.4 g) becomes inconsistent with the existence of the balanced rocks in the Mojave Desert and shortest lifetimes, highest birth rates, and longest return periods are required to allow existence of balanced rocks for large mean PHA and large ln({sigma}) scenarios (table 5). The scenario where fragility develops most rapidly (table 6) excludes all mean PHAs larger than 0.3 g and only allows the longest lifetimes and largest ln({sigma}) for mean PHAs of 0.2 g or less. The one attribute that is independent of fragility-curve scenario is the persistence of balanced rocks in the Mojave Desert from Brune (1996, 1999) consistent with the site-specific estimates of mean PHA and ln({sigma}) from Boore et al. (1997). The double-exponential fragility-curve scenario does require the highest birth rate (table 6) to be consistent with a lifetime of 20,000 years and the ground-motion parameters of Boore et al. (1997). However, the birth rate in table 6 would result in a peak population of eight independent rock sites in the absence of earthquake loadings, and there are eight separate sites in Brune (1999), so it's clear the balanced-rock birth rates used in table 6 are reasonable and certainly not excessive.


View this table:
[in this window]
[in a new window]

 
TABLE 4 Linear Fragility Curve: Persistent Balanced-Rock Parameters

 


View this table:
[in this window]
[in a new window]

 
TABLE 5 Exponential Fragility Curve: Persistent Balanced-Rock Parameters

 


View this table:
[in this window]
[in a new window]

 
TABLE 6 Double Exponential Fragility Curve: Persistent Balanced-Rock Parameters

 


Figure 8
View larger version (24K):
[in this window]
[in a new window]

 
Figure 8. Plan view of receiver positions and surface velocities. The grey region is rock with a surface Vs=1250 m/s and the white region is grus with a surface Vs=300 m/s. Dots are receiver positions in 43 rings of 90 receivers at 4-degree spacing. A radius spacing of 5 m was used within 50 m of the rock-grus contact. Velocity-density relations are shown in Figure 9 and velocity model depth parameters in table 7 and figure 10.

 
If the rock site ln({sigma}) = 0.4 results of Rodriguez-Marek et al. (2001) from the 1994 M 6.7 Northridge and 1989 M 7.0 Loma Prieta earthquake prove appropriate in general, tables 4, 5, 6 indicate that even larger mean PHAs than predicted by Boore et al. (1997) for a Vs30 = 1,250 m/s are allowed by the persistence of the balanced rocks in the Mojave Desert in Brune (1999). A ln({sigma}) of 0.4 consistent with Rodriguez-Marek et al. (2001) allows mean PHAs 50% larger than Boore et al. (1997) to be consistent with the persistence of balanced rocks in the Mojave Desert. This becomes a progressively more important issue as annual exceedance probability (AEP) become smaller in PSHA investigations.


    3D SITE EFFECTS
 TOP
 INTRODUCTION
 EFFECT OF CENSORING...
 SIMPLIFIED PSHA ANALYSIS
 MONTE CARLO SIMULATIONS OF...
 3D SITE EFFECTS
 VARIATION OF PEAK SPECTRAL...
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
In addition to balanced rocks, important engineered structures, such as concrete dams, tend to be situated on the stiffest materials in a region. Using 3D finite-difference ground-motion simulations with correlated-random velocity fluctuations in the top 2 km of the crust, O'Connell (1999) found systematic reduction of mean PHA and PHA dispersion with increasing local site velocities. Here a deterministic 3D-velocity structure more relevant to the conditions at the Mojave Desert balanced-rock sites is used to compare ground-motion responses near a lateral discontinuity between rock and grus. A 3D cylindrically symmetric model in plan view (figure 8) has a 1-km-radius core of rock with Vs = 1,250 m/s at the surface surrounded by a progressively increasing depth of grus with Vs = 300 m/s at the surface (figures 9 and 10), with a Vs vertical velocity gradient of 0.5 m/s/m in both the rock and grus portions of the model located above 2.5 km depth. A 10-m node spacing is used in the shallow (< 0.2) km portion of the model that was increased to 30-m node spacing at depths > 0.2 km using the 3D finite-difference approach of Liu and Archuleta (2002). The model is 8 km wide and 3.5 km deep. The 3D variable-spacing viscoelastic finite-difference method of Liu and Archuleta (forthcoming) is used with the attenuation parameters indicated in table 7 to simulate ground motions to 5.5 Hz associated with a uniform-amplitude, plane SH-wave incident 10 degrees from vertical incidence in the homogenous depth region of the model (> 2.5 km).


View this table:
[in this window]
[in a new window]

 
TABLE 7 3D Model Attenuation Parameters

 
Acceleration seismograms were extracted from the SH-wave-polarized horizontal component from 90 sites located at 4-degree intervals in each of 43 receiver rings at a range of radii from the center of the model (figure 8). A 5-m-radius increment was used for the receiver rings within 50 m of the rock-grus contact. The thinnest grus, 10 m, was next to the contact. Grus thickness increased in 10-m increments at increasing distances from the contact (figure 10). Acceleration response spectra were calculated for a wide range of periods using 5% damping. There is a substantial contrast in short period (0.2 s) and longer period (1 s) responses near the rock-grus contact (figure 10).


Figure 9
View larger version (14K):
[in this window]
[in a new window]

 
Figure 9. Relationships between P-wave velocity and Vp/Vs and density in the 3D velocity model.

 


Figure 10
View larger version (26K):
[in this window]
[in a new window]

 
Figure 10. Mean spectral accelerations (5% damping) as a function of position relative to the rock/grus boundary for periods of 0.2 s (A) and 1.0 s (B) and corresponding ln({sigma}) from the 90 receivers in each ring for periods of 0.2 s (C) and 1.0 s (D). Dotted vertical lines indicate position of rock/grus boundary. The dashed curves in (A) and (B) show the thicknesses of grus as a function of distance from the rock/grus boundary.

 
The maximum 0.2-s PSA response occurs 300 m from the grus-rock contact (figure 10), although the ratio of the depth to the rock divided by the wavelength equals 0.83, not one-quarter of a wavelength at this position. The second largest 0.2-s PSA amplification occurs close to the grus/rock contact (figure 10), although the ratio of the depth to the rock divided by the wavelength equals 0.33 and not one-quarter at this position. These deviations from one-quarter wavelength amplification reflect the amplification complexity introduced by slightly nonvertical incidence, a first-order lateral velocity discontinuity at the free surface, and an irregular increasing thickness of grus with increasing distance from the contact that results from the interference of direct S waves with edge waves produced at the rock/grus contact and along the dipping interface between rock and grus. In contrast, the maximum amplification at 1 s occurs at the position in the grus corresponding to a quarter wavelength (~ 700 m from the rock-grus contact in figure 10). This quantitatively illustrates that even a quarter-wavelength approximation is inadequate at short periods to explain amplifications near significant lateral velocity contrasts, much less simpler amplification approximations that employ Vs30. The PSA ln({sigma}) results in figure 10 nearly mimic the empirical ground-motion findings of Rodriguez-Marek et al. (2001) that thin-soil short-period ln({sigma}) for soil is about 0.15-0.2 larger than ln({sigma}) for rock.


Figure 11
View larger version (20K):
[in this window]
[in a new window]

 
Figure 11. Mean response spectral acceleration ratios of 0.25 s response to 1.0 s response as a function of distance from the rock/grus boundary. Dotted vertical line indicates position of rock/grus boundary. Dotted lines indicate 25% and 50% amplification positions in the grus relative to rock sites 10 m from the rock/grus boundary.

 
To facilitate comparison with the site amplification investigations of Stirling et al. (2002), the ratio of 0.25 s to 1 s PSA responses are plotted as a function of position relative to the rock grus boundary (figure 11). Within 7.5-20 m of the rock, where the grus is thinnest, the 3D-synthetic plane-SH-wave PSA amplifications at 0.25 s are 25%-50% (figure 11), consistent with the same short-period amplifications observed by Stirling et al. (2002) adjacent to the Mojave Desert balanced rocks in 7-15-m thicknesses of grus (as indicated by the velocity profiles of Abbott et al. 2001). Thus, the balanced-rock sites that are nunataks of rock sticking out of a sea of grus are not likely to be amplified at short periods (< 0.25 s) as suggested by Stirling et al. (2002) precisely because they are not located in grus.


    VARIATION OF PEAK SPECTRAL RESPONSE PERIOD WITH MAGNITUDE
 TOP
 INTRODUCTION
 EFFECT OF CENSORING...
 SIMPLIFIED PSHA ANALYSIS
 MONTE CARLO SIMULATIONS OF...
 3D SITE EFFECTS
 VARIATION OF PEAK SPECTRAL...
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
Somerville (2003) noted that as earthquake magnitudes increase, rise times increase for constant mean stress drops, resulting in a shifting of peak spectral accelerations to longer periods with increasing magnitude. O'Connell and Ake (2007) use isochrone analyses to show that near-fault ground motions from strike-slip earthquakes will have peak spectral acceleration responses at longer periods than dip-slip earthquakes. The available near-fault strike-slip strong-motion data (table 8 and figure 12) are certainly more consistent with maximum spectral accelerations in the 0.35-0.7 s period range than the 0.2-0.25 s predictions from Abrahamson and Silva (1997), Boore et al. (1997), and Sadigh et al. (1997). Anooshehpoor et al. (2004) demonstrated that rock-toppling probabilities decrease as ground motions are less enriched in short-period energy. Certainly, if the available near-fault strike-slip strong-motion data are relatively depleted of short-period energy for strike-slip earthquakes compared to many existing ground-motion relations, it's reasonable to assume that ground motions in the 15-35 km distance range would be depleted of short-period energy also relative to the spectral-shape predictions of Abrahamson and Silva (1997), Boore et al. (1997), and Sadigh et al. (1997).


View this table:
[in this window]
[in a new window]

 
TABLE 8 Near-Fault Strike-Slip Ground Motion Used in Figure 12

 

    DISCUSSION
 TOP
 INTRODUCTION
 EFFECT OF CENSORING...
 SIMPLIFIED PSHA ANALYSIS
 MONTE CARLO SIMULATIONS OF...
 3D SITE EFFECTS
 VARIATION OF PEAK SPECTRAL...
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
The Monte Carlo balanced-rock life-cycle simulations are quite sensitive to the assumed shape of the toppling fragility curve, particularly the rate of degradation of stability over the lifetime of erosion to ultimate static instability. These simulations show that it's possible to make some statistical inferences with prior knowledge of the life-cycle characteristics of a population of balanced objects by statistically simulating balanced-object life cycles over long durations relative to inter-event recurrence times. However, some prior population information is still required to make inferences about PHA recurrence rates, although not exclusively in the form indicated by Housner (1963). As demonstrated here, quantitative rates for the static creation and destruction of balanced objects provide sufficient prior information to make quantitative estimates of PHA recurrence. This is information in addition to the durations that balanced objects exist with quantifiable fragility characteristics. Geomorphic investigations beyond the scope of Bell et al. (1998) may quantify the static components of balanced-rock creation and destruction.

Boore and Joyner (1997) advocate moving beyond Vs30 to using a quarter-wavelength approximation to quantify site effects. The 3D-synthetic plane-SH-wave simulation results indicate this approach is likely to work well at sites not located close to large lateral velocity contrasts. However, the 3D site-response results (figures 10 and 11) and Rodriguez-Marek et al. (2001) provide a first-order quantitative indication of the perils of extrapolating thin-soil responses to rock sites (e.g., not accounting for the actual site conditions beneath the balanced rocks or your site of interest). Both mean PHA and ln({sigma}) are larger at soil sites adjacent to rock-pedestal sites (figures 10 and 11). The balanced rocks in the Mojave Desert in Brune (1999) are consistent with the results of Boore et al. (1997) and Rodriguez-Marek et al. (2001) that indicate that ground shaking is reduced as a function of increasing site shear-wave velocity. However, well-constrained rock lifetimes, birth rates, and time-dependent fragility characteristics are required to further quantify ground-motion constraints based on balanced-rock observations.


Figure 12
View larger version (13K):
[in this window]
[in a new window]

 
Figure 12. Median fault-normal mean 5 %-damped acceleration response spectra from nine near-fault (0.6-3.5 km) M > 6 strike-slip earthquake recordings is the solid curve and the mean is the dotted curve. Shadded region is the period band where Abrahamson and Silva (1997), Boore et al. (1997), and Sadigh et al. (1997) predict the largest spectral acceleration will occur.

 
The difficulties of incorporating balanced-rock information into PSHA analyses in many ways parallels the evolution of incorporating paleohydrologic data in flood frequency analyses, with some significant differences. The statistical methodology to incorporate paleohydrologic information in flood frequency analyses (Russell 1982; Stedinger and Cohn 1986) was developed nearly concurrently with the beginning of intensive efforts to collect and analyze paleohydrologic field data, as summarized in Baker et al. (2002). In contrast, balanced-rock data have not been rigorously incorporated into PSHA analyses, but have instead been presented in simple comparisons to PSHA hazard curves (Brune 1999; Brune et al. 2006). Quantifying age-dating uncertainties is a fundamental requirement in using paleohydrologic data in flood frequency analyses and balanced rocks in PSHA calculations. Relatively simple binary erosion/stability interpretations are most often used with paleohydrologic bound data (Levish 2002), although more rigorous quantitative approaches have been used (Ostenaa and O'Connell 2005). In contrast, considerable effort has been devoted to quantifying balanced-object fragility (Yim et al. 1980a,b; Shi et al. 1996; Zhang and Makris 2000; Anooshehpoor et al. 2000; Makris and Zhang 2001; Anooshehpoor and Brune 2002; Anooshehpoor et al. 2004; Purvance 2004). However, there has not been a rigorous statistical implementation of balanced-object fragility complexities that exist in dynamic toppling behavior, although complex behavior like multimodal toppling states could easily be incorporated into the discrete fragility probability distributions used in flood frequency analyses by O'Connell et al. (2002) and O'Connell (2005). For instance, instead of using the simple mean toppling PHA in the binary toppling/survival fragility curves (figure 7) in the Monte Carlo simulations, we could implement probability distributions within discrete PHA bins that incorporate appropriate toppling probabilities in each bin. Both balanced-rock and paleohydrologic approaches started first with simple 1D interpretations of hydrologic models for step-back-water modeling (O'Connor and Webb 1988; Webb and Jarrett 2002) or seismic wave propagation and site response (Stirling et al. 2002). However, as demonstrated in Denlinger et al. (2002) and O'Connell et al. (2002), rigorous 2D hydraulic modeling is necessary to obtain unbiased peak-discharge/paleohydrologic bound probabilities, and as figures 10 and 11 illustrate, it's necessary to account for realistic site conditions immediately beneath balanced rocks. Paleohydrologic-bound information has often been erroneously dismissed as irrelevant because of perceived large discharge estimation uncertainties or as unusable in flood frequency analyses for quantifying the frequencies of low AEP extreme floods (Hosking and Wallis 1997; National Research Council 1999). This led to the Bayesian flood frequency implementations of O'Connell et al. (2002) and O'Connell (2005) to explicitly include complex fragility and age-dating uncertainties to obtain rigorous quantitative statistical estimates of flood frequency. Similar statistical advances are needed to extract quantitative PSHA information from balanced rocks. Promising approaches include the Monte Carlo approaches used here or Bayesian approaches (Toro and Cornell 2006).


    CONCLUSIONS
 TOP
 INTRODUCTION
 EFFECT OF CENSORING...
 SIMPLIFIED PSHA ANALYSIS
 MONTE CARLO SIMULATIONS OF...
 3D SITE EFFECTS
 VARIATION OF PEAK SPECTRAL...
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
In the absence of information about the prior static population of balanced objects in the manner envisioned by Housner (1963), we've shown it's sufficient to quantify four types of information to estimate ground-motion recurrence from balanced objects:

  1. The total duration that a balanced object persists (the clock starts when erosion is sufficient that the object can be toppled);
  2. balanced-object fragility in response to ground shaking as a function of time since the onset time that erosion or some other process created a detached pedestal capable of toppling;
  3. pseudostatic rate that balanced objects are created and destroyed independent of earthquake loadings; and
  4. a sufficient sample size of balanced rocks (more than just one or two).

Further, to test PSHA performance, it's necessary to obtain a spatially diverse sampling at sufficient site separations to render ground motions uncorrelated. Particularly near faults, peak ground motion amplitudes can become strongly correlated at site distance separations that are smaller than correlation distances associated with coherent patches of seismic radiation from faults (asperities). For example, mean asperity diameters for M > 6.5 earthquakes are > 6 km (Somerville et al. 1999). Increased peak ground motion amplitudes are likely to be strongly correlated among sites separated by less than 1-2 km that are located near asperities of M > 6.5 earthquakes. Thus, large populations of rocks distributed over large areas, as in Brune (1999) or Brune et al. (2006), have more statistical information (less correlation and potential for bias) than samples from only a few distinct sites.

Through Monte Carlo simulations of balanced-rock life cycles we show that the existing empirical ground-motion relation that explicitly accounts for site velocity (Boore et al. 1997) is consistent with PSHA assumptions (Cornell 1968) and the persistence of balanced rocks east of the San Andreas fault in the Mojave Desert (Brune 1996, 1999). Instead of demonstrating that there is a fundamental inconsistency between PSHA assumptions and the persistence of balanced rocks (Anderson and Brune 1999), the Monte Carlo simulations show quite the opposite. That is, once one accounts for what is physically known about site response in empirical data (Boore et al. 1994, 1997), through theoretical investigations (Boore and Joyner 1997; O'Connell 1999), and 3D wave propagation (figures 10 and 11), there is consistency between the persistence of balanced rocks (Brune 1996, 1999) and PSHA assumptions (Cornell 1968). Empirical strong-motion data (Rodriguez-Marek et al. 2001) indicate that sites that are truly rock sites don't shake as hard as thin soil sites, and that rock sites have smaller peak amplitude dispersion than stiff soil sites; stiff soil sites form the bulk of the data often used in "rock" ground motion relations (e.g., Abrahamson and Silva 1997; Sadigh et al. 1997; Campbell 1997; Spudich et al. 1999). There is complete consistency between observation (Stirling et al. 2002) and simulation (figures 10 and 11) that site acceleration responses on thin soft grus adjacent to rock pedestals (Abbott et al. 2001) are amplified 25%-50% at 4 Hz (Stirling et al. 2002) relative to the responses on the exclusively rock-pedestal sites. Nearly all normal faults cataloged in the western United States have hangingwall basins consisting of several hundred meters to several kilometers of low-velocity sedimentary fill (Zoback 1983) that are likely to strongly amplify hangingwall ground motions relative to footwall sites that typically are located on competent rock over a broad frequency band (O'Connell et al. 2003, 2007). Thus, it is clear that it probably won't be possible to quantify potential differences between hangingwall and footwall ground motions during normal faulting, which could be substantial (Oglesby et al. 1998, 2000; Brune and Anooshehpoor 1999; Bouchon et al. 2000; Brune 2003; Shi et al. 2003; Schürch and Becker 2005; O'Connell et al. 2007), using balanced objects if 3D crustal velocity structure and shallow site-response effects are ignored.

Finally, it would be a grave error to use rock-response results in most built environments because many urban areas are constructed on stiff soils; most ground motion relations (e.g., Abrahamson and Silva 1997; Sadigh et al. 1997; Campbell 1997; Spudich et al. 1999) appear to have appropriate means and log deviations for stiff-soil sites (Rodriguez-Marek et al. 2001). Possible revisions to estimating ground-motion relations that incorporate balanced-rock information will need to rigorously account for the site velocity effects to avoid introducing biases in estimated stiff-soil sites that could result in systematic underestimation of ground-motion amplitudes.


    ACKNOWLEDGMENTS
 TOP
 INTRODUCTION
 EFFECT OF CENSORING...
 SIMPLIFIED PSHA ANALYSIS
 MONTE CARLO SIMULATIONS OF...
 3D SITE EFFECTS
 VARIATION OF PEAK SPECTRAL...
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank Jon Ake and David Boore for helpful reviews. This research was supported by the Bureau of Reclamation Dam Safety Research Program as part of the SITER project.


    REFERENCES
 TOP
 INTRODUCTION
 EFFECT OF CENSORING...
 SIMPLIFIED PSHA ANALYSIS
 MONTE CARLO SIMULATIONS OF...
 3D SITE EFFECTS
 VARIATION OF PEAK SPECTRAL...
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 

Abbott, R. E., J. N. Louie, J. N. Brune, and R. Anooshehpoor (2001). Analysis of Shallow Site Response to LARSE-2 Blasts at Precarious Rock Sites near the San Andreas Fault. Final report to Southern California Earthquake Center.

Abrahamson, N. A., and W. J. Silva (1997). Empirical response spectral attenuation relations for shallow crustal earthquakes. Seismological Research Letters 68, 94-109.[GeoRef]

Anderson, J. G., and J. N. Brune (1999). Probabilistic seismic hazard analysis without the ergodic assumption. Seismological Research Letters 70, 19-28.[GeoRef]

Anooshehpoor, A., and J. N. Brune (2002). Verification of precarious methodology using shake table tests of rocks and rock models. Journal of Soil Dynamics and Earthquake Engineering 22,917 -922.[CrossRef]

Anooshehpoor, A., J. N. Brune, and Y. Teng (2004). Methodology for obtaining constraints on ground motion from precariously balanced rocks. Bulletin of the Seismological Society of America 94,285 -303, doi:10.1785/0120020242 .[Abstract/Free Full Text][CrossRef][ISI][GeoRef]

Anooshehpoor, A., T. H. Heaton, B. Shi, and J. N. Brune (2000). Reply to comment on "Estimates of ground accelerations at Point Reyes, California during the 1906 San Francisco earthquake" by J. Zhang and N. Makris. Bulletin of the Seismological Society of America 90,1,349 -1,351.[Free Full Text][CrossRef][ISI]

Baker, V. R., R. H. Webb, and P. K. House (2002). The scientific and societal value of paleoflood hydrology, In Paleoflood Hydrology, ed. P. K. House, D. R. Levish, R. H. Webb, and V. R. Baker, 1-19. Washington, DC: American Geophysical Union Monograph, 385 pps.

Bell, J. W., J. N. Brune, L. Tanzhuo, M. Zreda, and J. C. Yount (1998). Dating precariously balanced rocks in seismically active parts of California and Nevada. Geology 26,496 -498.

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

Boore, D. M., and W. B. Joyner (1997). Site amplifications for generic rock sites. Bulletin of the Seismological Society of America 87,327 -341.[Abstract/Free Full Text][ISI][GeoRef]

Boore, D. M., W. B. Joyner, and T. E. Fumal (1997). Equations for estimating horizontal response spectra and peak acceleration from western North American earthquakes: A summary of recent work. Seismological Research Letters 68,128 -153.[GeoRef]

Bouchon, M., S. Gaffet, C. Cornou, M. Dietrich, J. P. Glot, F. Courboulex, A. Caserta, G. Cultrera, F. Marra, and R. Guiguet (2000). Observations of vertical ground accelerations exceeding gravity during the 1997 Umbria-Marche (central Italy) earthquakes. Journal of Seismology 4, 517-523, DOI 10.1023/A:1026579908823. http://www.springerlink.com/content/v04k1155122445j2/

Brune, J. N. (1996). Precariously balanced rocks and ground motion maps for southern California. Bulletin of the Seismological Society of America 86, 43-54.[ISI][GeoRef]

____. (1999). Precarious rocks along the Mojave segment of the San Andreas fault, California: Constraints on ground motion for great earthquakes. Seismological Research Letters 70, 29-33.[GeoRef]

____. (2003). Precarious rock evidence for low near-source accelerations for trans-tensional strike-slip earthquakes. Physics of the Earth and Planetary Interiors 137,229 -239.[CrossRef][ISI][GeoRef]

Brune, J. N., and A. Anooshehpoor (1999). Dynamic geometrical effects on strong ground motion in a normal fault model. Journal of Geophysical Research 104,809 -815.[CrossRef][ISI]

Brune, J. N., A. Anooshehpoor, M. D. Purvance, and R. J. Brune (2006). Band of precariously balanced rocks between the Elsinore and San Jacinto, California, fault zones: Constraints on ground motion for large earthquakes. Geology 34,137 -140.[Abstract/Free Full Text][CrossRef][ISI][GeoRef]

Campbell, K. C. (1997). Empirical near-source attenuation relationships for horizontal and vertical components of peak ground acceleration, peak ground velocity, and pseudo-absolute acceleration response spectra. Seismological Research Letters 68,154 -179.[GeoRef]

Cao, T., W. A. Bryant, B. Rowshandel, D. Branum, and C. J. Wills (2003). The Revised 2002 California Probabilistic Seismic Hazard Assessment Maps. California Geological Survey, http://www.consrv.ca.gov/CGS/rghm/psha/.

Cornell, C. A. (1968). Engineering seismic risk analysis. Bulletin of the Seismological Society of America 58,1,583 -1,606.[Abstract/Free Full Text][ISI][GeoRef]

Denlinger, R. P., D. R. H. O'Connell, and P. K. House (2002). Robust determination of stage and discharge: An example from an extreme flood on the Verde River, Arizona. In Ancient Floods, Modern Hazards: Principles and Applications of Paleoflood Hydrology, ed. P. K. House, R. W. Webb, V. R. Baker, and D. R. Levish, 127-146. Water Science and Application Series, vol. 5. Washington, DC: American Geophysical Union.

Hosking, J. R. M., and J. R. Wallis (1997). Regional Frequency Analysis—An Approach based on L-Moments. Cambridge: Cambridge University Press, 224 pps.

Housner, G. W. (1963). The behavior of inverted pendulum structures during earthquakes. Bulletin of the Seismological Society of America 53,404 -417.

Levish, D. R. (2002). Paleohydrologic bounds—Non-exceedance information for flood hazard assessment. In Ancient Floods, Modern Hazards: Principles and Applications of Paleoflood Hydrology, ed. P. K. House, R. W. Webb, V. R. Baker, and D. R. Levish, 175-190. Water Science and Application Series, vol. 5. Washington, DC: American Geophysical Union.

Liu, P.-C., and R. Archuleta (2002). The effect of a low-velocity surface layer on the simulated ground motion. Seismological Research Letters 73, 267.[GeoRef]

____. (forthcoming). Efficient modeling of Q for 3D numerical simulation of wave propagation. Bulletin of the Seismological Society of America.

Makris, N., and J. Zhang (2001). Rocking response of anchored blocks under pulse-type motions. ASCE Journal of Engineering Mechanics 127 (5), 484-493.[CrossRef]

National Research Council (1999). Improving American River Flood Frequency Analyses. Washington, DC: National Academy Press, 120 pps.

O'Connell, D. R. H. (1999). Replication of apparent nonlinear seismic response with linear wave propagation models. Science 283,2,045 -2,050, DOI:10.1126/science.283.5410.2045 .[Abstract/Free Full Text][CrossRef][ISI][Medline][GeoRef]

____. (2005). Nonparametric Bayesian flood frequency estimation. Journal of Hydrology 313, 79-96, doi:10.1016/j.jhydrol.2005.02.005 .[CrossRef][ISI][GeoRef]

O'Connell, D. R. H., and J. P. Ake (2007). Earthquake ground motion estimation. In Earthquakes: (Hazards and Disasters), ed. C. Rodrigue and E. Rovai.2 vols. New York: Routledge.

O'Connell, D. R. H., S. Ma, and R. Archuleta (2007). Influence of fault dip and near-fault crustal heterogeneity on normal-faulting rupture dynamics and ground motions. Bulletin of the Seismological Society of America 97 (6).

O'Connell, D. R. H., D. A. Ostenaa, D. R. Levish, and R. E. Klinger (2002). Bayesian flood frequency analysis with paleohydrologic bound data. Water Resources Research 38 (5),1,058 , doi:10.1029/2000WR000028 .[CrossRef]

O'Connell, D. R. H., C. K. Wood, D. A. Ostenaa, L. V. Block, and R. C. LaForge (2003). Ground Motion Evaluation for Jackson Lake Dam, Minidoka Project, Wyoming. Seismotectonic report 2003-2, Bureau of Reclamation, Denver, Colorado, 493 pps.; http://www.usbr.gov/pn/programs/srao_misc/jackson/.

O'Connor, J. E., and R. H. Webb (1988). Hydraulic modeling for paleoflood analysis. In Flood Geomorphology, ed. V. R. Baker, R.C. Kochel and P. C. Patton,393 -402. New York: Wiley.

Oglesby, D. D., R. J. Archuleta, and S. B. Nielsen (1998). Earthquakes on dipping faults: The effects of broken symmetry. Science 280,1,055 -1,059.[Abstract/Free Full Text][CrossRef][ISI][Medline][GeoRef]

Oglesby, D. D., R. J. Archuleta, and S. B. Nielsen (2000). The three-dimensional dynamics of dipping faults. Bulletin of the Seismological Society of America 90,616 -628.[Abstract/Free Full Text][CrossRef][ISI][GeoRef]

Ostenaa, D. A., and D. R. H. O'Connell (2005). Big Lost River Flood Hazard Study. Idaho National Laboratory, Idaho, U.S. Bureau of Reclamation seismotectonic report 2005-2, Denver, Colorado, 212 pps., six appendices.

Petersen, M. D., W. A. Bryant, C. H. Cramer, T. Cao, M. S. Reichle, A. D. Frankel, J. J. Lienkaemper, P. A. McCrory, and D. P. Schwartz (1996). Probabilistic Seismic-hazard Assessment for the State of California. California Division of Mines and Geology Open-File Report 96-08, USGS Open-File Report 96-706.

Purvance, M. (2004). Parameters contributing to the toppling of precariously balanced rocks. In Proceedings and Abstracts of the SCEC 14th Annual Meeting,148 . September 19-23, 2004.