|
|
||||||||
Article |
1 Seismological Laboratory and
Department of Geological Sciences
University of Nevada
Reno, Nevada
89557
pancha{at}seismo.unr.edu
jga{at}unr.edu
(A.P.,
J.G.A.)
2 Nevada Bureau of Mines and
Geology
Mail Stop 178
University of Nevada
Reno, Nevada
89557-0088
kreemer{at}unr.edu
(C.K.)
| Abstract |
|---|
|
|
|---|
5. The catalog was
compiled from 15 preexisting catalogs, supplemented by the review of 44 journal
articles. Throughout the catalog compilation, care was taken to obtain the
moment magnitude or a reasonable, and not inflated, equivalent. Seventy-six
percent of the moment release occurred during 10 earthquakes of magnitude
MW
6.79. The spatial distribution of earthquakes and
their moment release matches the geodetic pattern of deformation. All three are
concentrated in a
200-km-wide zone along the western boundary of the Great
Basin, with this zone widening to the north. Several techniques, ultimately
traceable to Kostrov and Brune, are used to translate the geodetic strain rates
into rates of seismic-moment release. The agreement between geodetic and
seismic- moment rate suggests that, within uncertainties, the rate of historic
earthquakes within the Basin and Range province, taken as a whole, provides a
reasonable estimate for the future rate of seismicity. These results support the
hypothesis that even a few years of detailed geodetic monitoring can provide a
good constraint on earthquake occurrence rate estimates for large-enough
regions. | Introduction |
|---|
|
|
|---|
Comparisons of seismicity and geology (e.g., Anderson, 1979; Doser and Smith, 1989) or comparisons of seismicity and geodesy (e.g., Ward, 1998a, 1998b; Shen Tu et al., 1998; Kreemer et al., 2000, 2002; Masson et al., 2005) have been conducted. This study improves on previous studies by including a decade of geodetic data, an improved seismicity catalog, and an attempt to characterize all of the active faults.
The adequacy of seismic catalogs to estimate average regional earthquake occurrence rates is governed by the catalog duration (Smith, 1976; Ward, 1998a) and the regional strain rate (Ward, 1998a). For seismicity rates estimated from the historical earthquake catalogs to be valid, the average recurrence interval is required to be shorter than the historical record. For an individual fault a complete earthquake cycle is required. Alternatively, for a region containing multiple faults, the historical seismicity record is required to be long enough to capture a statistical sample of all phases of the seismic cycle, including of course earthquakes, but different parts of the cycle can be represented by different faults. With this constraint, the catalog duration is almost always too short to give a reliable occurrence rate estimate for regions the size of an urban area, as desired for seismic-hazard analysis.
Fault slip rates may be used to estimate average regional earthquake occurrence rates (e.g., Brune, 1968; Wallace, 1970; Anderson, 1979; Molnar, 1979, Doser and Smith, 1982). For such geological estimates of seismicity to be reliable, all major faults must be recognized and characterized correctly. Where faulting is concentrated on a few major through-going structures, as in coastal California north or south of the Transverse Ranges, there is reason for confidence that this goal is close to being achieved. In the Basin and Range, fault characterization is incomplete.
Two conditions for geodesy to give reliable estimates of earthquake occurrence rates are necessary. First, the geodetic measurements should sample a large enough spatial scale so that they are not affected by nonlinear strain accumulation during the earthquake cycle on individual faults. Second, they should sample a long enough time interval that measurement uncertainties have a minimal effect on the estimated velocities. In the Great Basin, geodetic data meeting these conditions are obtained with as little as a few years of observations using the Global Positioning System (GPS), but the deformation is not uniquely assignable to specific faults. Agreement between the strike of active Quaternary faults and the azimuth of contemporary deformation lead Hammond and Thatcher (2004) to conclude that geodetic motion can be used to infer deformation over many earthquake cycles. Contemporary strain rates may be used to predict average earthquake moment rates using methods ultimately traceable to Kostrov (1974) or Brune (1968), assuming that all of the strain that is accumulated is ultimately released in earthquakes (e.g., Anderson, 1979; Ward, 1994; Working Group on California Earthquake Probabilities, 1995; Savage and Simpson, 1997; Shen-Tu et al., 1998; Ward, 1998a, 1998b).
The Basin and Range province extends from the rigid Sierra Nevada block in
the west to the Colorado Plateau in the east
(Fig. 1). The province is an
actively deforming region of Cenozoic extension and shear, dominated by normal
faulting throughout with strike-slip deformation superimposed primarily along
the western margin of the Great Basin
(Stewart, 1988). High heat flow,
negative Bouguer gravity anomalies, thin crust and lithosphere, and
high-attenuation low-velocity upper mantle characterize the region
(Catchings and Mooney, 1991;
Jones et al., 1992;
Chulick and Mooney, 2002;
Zandt et al., 1995).
From the geodynamic viewpoint, buoyancy forces within the crust and lithosphere
(Wernicke, 1992;
Sonder and Jones, 1999),
tractions applied to the plate edge
(Wernicke, 1992), and forces
applied to the base of the lithosphere
(Wernicke, 1992;
Sonder and Jones, 1999) are all
potentially contributing to drive the extension and shear observed in the
region. Topographic and geodetic data, along with plate-motion constraints,
indicate that extension is driven in part by gravitational potential energy, but
plate-boundary interaction stresses modify extension directions
(Flesch et al., 2000).
In addition to spreading, about 25% of the Pacific–North American
(PA–NA) relative plate motion (
12 mm/yr) is taken up by
displacement and deformation in the Basin and Range province
(Dokka and Travis, 1990;
Dixon et al., 2000;
Bennett et al., 2003).
|
Because deformation is distributed over a region nearly 1000 km wide, most major Basin and Range faults have recurrence times of several thousand years (Wallace, 1984; McCalpin and Nishenko, 1996; Lee et al., 2001; Dixon et al., 2003). We propose, however, that the deformation rate is fast enough that meaningful comparison of average seismicity rates estimated using seismic, geological, and geodetic data can still be accomplished for this region as a whole. The distributed seismicity of the region is assumed to be caused by a sufficient number of faults at different stages of the earthquake cycle to compensate for the long recurrence interval of the individual faults. This is reasonable because there are more than 400 major range-front faults distributed throughout the region (dePolo, 1998). Unfortunately many of the slip rates available for the region are based on reconnaissance techniques (dePolo, 1998) rather than on more detailed trenching. Geodetic strain measurements averaged over the region should be representative of the total geological strain across the entire region, considering that geodetic strain rates based on a decade of observations are broadly consistent with global plate models based on rates of seafloor spreading. Shen-Tu et al. (1998) and Humphreys and Weldon (1994) showed that fault slip rates in the western United States add up to be close to geological and geodetic estimates of the PA–NA motion. Although the geodesy data contain transients and localized areas of higher strain, these are small anomalies compared with the greater region.
To compare the three methods, we compile the best available GPS, geology, and seismicity data. GPS data from several studies in the Basin and Range area are combined for this study to model the present-day deformation field. Best estimates of slip rates on the most active faults characterized are obtained from input to the 1996 and 2002 U.S. Geological Survey (USGS) seismic hazard maps (Frankel et al., 1996, 2002; Haller et al., 2002).
To characterize the historical seismicity, this study compiles the most complete possible seismic-event catalog for the Basin and Range with strong emphasis on obtaining the most appropriate moment magnitude (MW) for each event. We utilize geodetic data from multiple stations in the Basin and Range to define the geodetic deformation on a fine spatial scale (5-km grid). This is in contrast to previous studies for the Basin and Range that utilized only a small number of data for this region (Ward, 1998a; Shen Tu et al., 1998). Besides the comparison of the entire region, we also compare the spatial distribution of earthquakes and their moment release with crustal deformation rates.
| Earthquake Data |
|---|
|
|
|---|
Seismic moment, M0, of an earthquake is defined in terms
of source parameters as the product of average shear modulus µ of
the crustal rock around the earthquake, area Ae of the fault
ruptured in the earthquake, and D, the average slip during the
earthquake (Brune, 1968), that
is,
|
| (1) |
We estimated seismic-moment rates from a new catalog of earthquakes from 1850
to the end of 1999, intended to be complete for magnitude M
5
(Fig. 1). Earthquakes within
the study region with M
4.8 in any of 15 pre- existing catalogs
were selected (Table 1). After
removing duplicate entries, the catalog contained a total of 800 earthquakes
since 1855, including 487 earthquakes with M
5.0
(www.seismo.unr.edu/htdocs/BandR.html).
By allowing a lower cutoff level of M 4.8, we sought to ensure that
magnitude differences arising from the use of different catalogs is accounted
for. In this way, no significant events should have been neglected from the
catalog.
|
This new catalog was supplemented by the review of 44 published journal articles to obtain MW values for many of the more prominent earthquakes (Slemmons, 1957; Tsai and Aki, 1966; Savage and Hastie, 1966, 1969; Bolt and Miller, 1975; Hanks et al., 1975; Hart et al., 1977; Langston and Butler, 1976; Hanks and Kanamori, 1978; Toppozada et al., 1981; Doser and Smith, 1982; Barrientos et al., 1985; Boatwright and Choy, 1985; Doser, 1985; Doser and Smith, 1985; Ekstrom and Dziewonski, 1985; Nabelek et al., 1985; Patton, 1985; Stein and Barrientos, 1985; Barker and Wallace, 1986; Doser, 1986; Sipkin, 1986; Ward and Barreintos, 1986; Doser, 1987; Doser and Kanamori, 1987; Barker and Doser, 1988; Doser, 1988; Pacheco and Nábelek, 1988; Patton and Doser, 1988; Doser, 1989a, 1989b; Doser and Smith, 1989; Westaway and Smith, 1989; Rogers et al., 1991; Smith and Arabasz, 1991; Kanamori et al., 1992; Beanland and Clark, 1993; Wells and Coppersmith, 1994; Caskey et al., 1996; Mason, 1996; dePolo and dePolo, 1999; dePolo et al., 2003; Ji et al., 2002; Ichinose et al., 2003). Only best estimates of the moment magnitude were considered from the literature. Some of the larger events consisted of a number of subevents. In these cases, the moments of each of the individual subevents were summed and the moment magnitude was calculated from that total.
The initial compiled catalog contained multiple entries from the same earthquake. Discrepancies were noted between catalog listings for single events. Timing differences of up to one day were observed, as well as differences in location, especially for early events. Some events were not listed in their primary catalogs. Hence, a small amount of subjectivity, based on the similarities of both the location and timing of events, was introduced for event association between catalogs. Errors were often noted in secondary sources. Where available, primary sources were therefore preferred.
Emphasis was placed on the accuracy of the magnitudes within each catalog listing, as discussed in the next section. The main objective was to include all the large earthquakes within the region and gain a reasonable estimate of MW or a reasonable equivalent. Preferred locations were those from primary catalogs, except where geological data have been used in relocations. In general, less emphasis was placed on the accuracy of the timing of events.
Large uncertainties surround some large historic events. Several catalogs include an earthquake in 1852 in western Nevada with M 7.3. The anecdotal evidence for this earthquake is not sufficient to assign a magnitude and location that is reliable enough to use in this study. The occurrence of the 1903 M 6.5 earthquake (Slemmons et al., 1959; Rogers et al., 1991) is based on geological mapping, field studies, and interviews with residents. Although omitted by many other catalogs, we include this event in our final catalog as its occurrence is sufficiently credible. Because of its smaller size, it does not contribute greatly to the final cumulative seismic-moment release of the region.
Earthquakes that are potentially related to volcanic processes, including
those associated with activity around Mammoth Lakes and Mt. St. Helens, were
included in the catalog. The Mt. St. Helens events were of magnitude M
5.3 or less and, hence, contribute little to the overall seismic- moment
release. The largest events near Mammoth Lakes (M
6) occur
outside the caldera, so that they may appropriately be considered to contribute
to the tectonic deformation of the region. Events in the locality of the Nevada
Test Site, occurring at times of known nuclear blasts, were removed from the
final database. Three additional events were removed, as they were located on
the Nevada Test Site and occurred either exactly on the hour or half hour. All
other events in the vicinity are assumed to be of tectonic origin or triggered
by the blasts.
Amplitude spectra of surface waves were used by
Patton (1985) to determine
seismic moments of western United States events with ML and
mb between 4.3 and 5.5. Whereas magnitude estimates from
other sources were M
4.8, the moment magnitudes estimated by
Patton (1985) were often less
than this. This discrepancy arises because of the high- frequency content of
these smaller events. Based on the rules for magnitude selection below, these
smaller moment magnitudes are used in the analyses.
Doser and Smith (1982) found that some events with ML < 4.8 can be modeled to give an estimate of MW > 4.8; these events are included in the new catalog. Thus, we may not have achieved completeness in the MW 4.8–5.5 range. However, these events would not have an impact on the main results of this study because of their small size.
Geological parameters given by Wells and Coppersmith (1994) and Mason (1996) were used to calculate moment- release values for significant earthquake events from equation (1) using a shear modulus of µ = 3 x 1011 dyne cm2. To determine Ae for normal faults, the dip was assumed to be 60°, whereas for strike-slip events, the dip was assumed to be 90°. Mason (1996) did not give estimates of the vertical depth of the faulting. Those given by Wells and Coppersmith (1994) were therefore used in calculation of the seismic moment from source parameters listed by Mason (1996). Best estimates of the fault length and displacement, as given by Mason, were used. In the case of the Wells and Coppersmith (1994) data, the surface fault length and average displacement were used. For four small earthquakes where average displacements were not available, maximum displacements were used. Equation (1) was also applied to other earthquake source parameter data listed within the literature.
These calculated moment values and those within the literature were converted
to moment magnitude estimates, so that they could be compared with other
magnitude estimates for each respective event using the relation first defined
by Kanamori (1977):
|
| (2) |
| Magnitude Assignment |
|---|
|
|
|---|
Careful selection criteria were established to retain the best MW estimate from the literature and catalog listings for the largest events. The Harvard long-period surface-wave estimates of the seismic moment have been consistent for the past 28 years; hence, these estimates were given primary preference. Second preference was given to other surface- wave moment tensor estimates followed by any other surface-wave estimates. Long-period surface waves are considered to be more representative of the average faulting process than body waves and thus were given precedence over moment estimates based on body waves. Within each of these MW categories, the catalog record with the minimum-magnitude listing was retained. After surface waves, body wave and/or Pnl wave, geological, and leveling data were then considered in the same category. In these cases, we needed to use judgment based on the quality of the data; because of concerns that each of these have the potential to underestimate the size of the earthquakes, we tended to favor larger values. Appendix 1 describes in detail how preferred MW values were selected for each of the 10 largest earthquake events in the catalog, and are listed in Table 2.
|
All magnitudes within the catalog were then treated as moment magnitudes. We
then estimated the seismic moment of each event using the relation
|
| (3) |
Figure 2 shows MW and coda magnitude estimates for 113 moderate-magnitude earthquakes from 1990 to 2000 in the western Great Basin. Moments are determined by Ichinose et al. (2003). Although these smaller events have relatively little influence on the study results, we use these results as an indication that our use of network magnitude estimates in the place of MW does not cause major bias.
|
| Strain Rate Field Model from GPS Velocities |
|---|
|
|
|---|
50
km) and western (
200 km) margins of the Basin and Range, coinciding with
regions of modern seismicity, and with little deformation in between
(Dixon et al., 1995;
Thatcher et al., 1999;
Dixon et al., 2000;
Svarc et al., 2002b;
Bennett et al., 2003;
Hammond and Thatcher, 2004).
Concentration of deformation in the westernmost 200 km of the Basin and
Range and along the eastern boundary may be related to rheological weakness of
the lithosphere
(Thatcher, 2003). Margin
parallel velocities at a latitude of 39–40° N increase from
1 to
2 mm/yr at 117.7° W, representative of the relatively more stable interior
of the Basin and Range, to
12 mm/yr at 120° W
(Bennett et al., 2003;
Hammond and Thatcher, 2004).
Strain rates also increase from north to south along the western boundary of the
region
(Bennett et al., 2003)
because of narrowing of the high- deformation zone from the northern Walker Lane
to the ECSZ in the south. For the purpose of this study it is appropriate to model the present-day deformation field by means of a continuous strain rate field based on GPS velocity observations. For this, we apply a spline interpolation technique (Haines and Holt, 1993; Holt et al., 2000) in which model velocities are fitted to observed GPS velocities in a least-squares sense, using the full data covariance matrix. After a continuous velocity gradient tensor field model has been obtained, we calculated model velocities for 0.05° (5 km) grid intervals to do the analysis in sections: Spatial Distribution of Moment and Geodetic Moment Rates.
We combine GPS velocities from multiple studies in the Great Basin area. Some of these studies are published (Freymueller et al., 1999; McClusky et al., 2001; Oldow et al., 2001; Svarc et al., 2002a, 2002b; Mazzotti et al., 2003; Hammond and Thatcher, 2004; Savage et al., 2004), others have been made available otherwise: Eastern Basin-Range and Yellowstone Hotspot GPS Network (EBRY) (R. Smith, personal comm., 2003); Southern California Earthquake Center (SCEC) Crustal Motion Map v.3; Basin and Range Geodetic Network (BARGEN) (R. Bennett, personal comm., 2003). In addition we have used several unpublished USGS GPS results from campaign-style surveys: these consist of the Hawthorne Profile, (part of) the Mammoth network, and networks presented by Hammond and Thatcher (2003) and Hammond et al. (2004).
Each of the geodetic studies used has realized a unique frame of reference (in general, by assuming North America [NA] to be stable). To be consistent, we have applied Helmert transformations (when possible) to translate velocities into the global GPSVEL velocity solution (Lavallée et al., 2001) which is in the ITRF2000 frame. We subsequently rotate all velocities into a North America fixed reference frame, using the ITRF2000 pole obtained from the GPSVEL solution: 3.0° S, 83.1° W, 0.198°/Ma. For some studies we have multiplied standard errors in the velocity components with a factor: EBRY data times 10, BARGEN data times 10, (Savage et al., 2004) times 2, (McClusky et al., 2001) times 2, and (Oldow et al., 2001) times 20.
| Analysis and Results |
|---|
|
|
|---|
4.8 are complete since 1954, M
5.5
earthquakes are complete since 1932, M
6.0 earthquakes are
complete since 1901, and M
7.0 earthquakes are complete for the
entire duration of the catalog back to 1855.
|
|
7.0, 0.54 earthquakes
per year with MW
6.0, and 6.6 earthquakes per year with
MW
5.0. These b-value curves are sensitive to
the magnitudes assigned to each earthquake. Of the total moment, 76% was
released during 10 earthquakes of magnitude
MW
6.79
(Table 2), and 90% was released
in the 38 events of MW
6.1. This confirms the
observation that small events do not significantly release accumulating strain
(Brune, 1968;
Anderson, 1979;
Anderson and Luco, 1983;
Shen-Tu et al., 1998).
Spatial Distribution of Moment
Figure 5 shows the
boundaries of four domains (A–D) that we use to compare deformation and
seismicity. The four domains are each 300 km wide.
Figure 6 shows the magnitude of
crustal velocity, as determined from inversion of the geodetic model, along
profiles. These velocities give a smoother picture of the deformation field than
the calculated strain-rate field. A northward broadening of the zone of high
deformations is evident, as observed by
Bennett et al. (2003).
|
|
|
|
Figure 7a–d compares these deformation profiles with the spatial distribution of earthquake numbers and of moment release within each domain. The seismic moment of each earthquake is more completely represented as a tensor. Here we use the magnitude of the maximum eigenvalue. Although tensor information is available for the 10 largest earthquakes, which release 76% of the total seismic moment, use of tensors increases the number of degrees of freedom. Therefore a much longer observation time is required to obtain a reliable comparison with regional components of the geodetic strain, considering that there is randomness in fault orientations.
The spatial patterns of seismic activity, seismic moment release, and geodetic deformation are similar along all of the profiles and show a northward widening of the deformation zone along the western edge of the province. One way to quantify the similarity of the profiles in Figure 7 is to tabulate the distances from the southwest boundary to the point along the profile within which 75% of the total of each activity measure occurs (Table 3). For profiles A, C, and D these widths agree within 20%. Across domain B (Fig. 7b), 75% of the earthquake numbers and the geodetic deformation occur within a zone about 90–113 km wide, but the moment release is concentrated by the 1872 Owens Valley event (Table 2), the largest event in the catalog. Removal of the 1872 MW 7.58 event results in improved correspondence between the three curves (Fig. 7b; Table 3). The 1915 Pleasant Valley (domain D) and 1872 Owens Valley (domain B) events occurred prior to when seismic instrumentation was capable of observing aftershocks. If those aftershocks could be included the distribution of earthquake numbers would change.
|
Plots similar to those in Figure 7 for the eastern domain (Fig. 5) are shown in Figure 8. Scales along the ordinate axes are normalized by the largest values in the study region for this latitude range. This allows comparison with profiles in Figure 7. The figure shows that activity along the eastern half of the Great Basin, across the Wasatch Mountains (Fig. 5), is significantly smaller than in the west. The greatest increase on all three rates in Figure 8 occurs near the Wasatch Front at 112° W. About 13% of the earthquakes and 4% of the seismic moment are concentrated east of 112.2° W. Less than about 15% of the geodetic deformation occurs there. Malservisi et al. (2003), in a study of similar GPS data, find that it is not possible to distinguish between two models: one being elastic strain accommodation on multiple faults with generally low strain rates (except the Wasatch) and the other being postseismic creep on the Wasatch fault.
Figures 7 and
8 suggest that spatial
distribution of moderate (M
5) earthquakes and moment release are
correlated in this region. If true, this implies that earthquake numbers and
moment release could be used to constrain the geodetic deformation field and,
conversely, be predicted from geodetic strain rates.
Kreemer et al. (2002)
have suggested that this is true on a global scale, whereas
Masson et al. (2005)
find this relation does not hold at a regional scale for Iran.
Seismic Scalar Moment Rate
The statistical procedure used to estimate the historical seismic-moment rate
is illustrated in Figure 9,
showing cumulative seismic moment as a function of time. As lower- bound
magnitudes of moderate earthquakes were considered, this yields lower-bound
estimates of the moment release rate. A least-squares fit to the points in
Figure 9 (1 point for each
year with an earthquake) has a slope of 6.85 x 1025 dyne cm/yr.
Lines on Figure 9 show a
nonunique, but plausible rationale for moment rates as low as
6.15 x 1025, or as high as 8.84 x
1025 dyne cm/yr.
|
To quantify uncertainties associated with these seismic- moment rates, we repeated the procedure shown in Figure 9; (1) using upper-bound estimates of the magnitudes of the smaller events, (2) using smallest and largest moment estimates for the 10 largest events, and (3) using Monte Carlo realizations of the moments of the 10 largest events randomly distributed between the smallest and the largest values. As the estimate of moment rates is least sensitive to the range of uncertainty of the 790 smaller events in the catalog (of the total of 800), we focused on the uncertainties associated with the 10 largest events. Results are given in Table 4.
|
The procedure in Figure 9 was automated and repeated for randomly chosen moments of the 10 largest earthquakes for the Monte Carlo approach. The moment release for each of these events was randomly selected, assuming a constant probability density between minimum and maximum Mw estimates. The maximum and minimum Mw values were selected based on the most reliable and appropriate estimates of Mw from the literature (Appendix 1, Table 2). Moment release for the 790 smaller earthquake events was held constant at the favored values. A total of 50,000 Monte Carlo realizations was generated. The distributions of rates from these realizations are shown in Figure 10 and summarized in Table 4. Considering one standard deviation about the mean values and the full set of fitting approaches, the most likely historical moment rate ranges from 5.8 x 1025 to 11.3 x 1025 dyne cm/yr determined from averaging since 1871. The extremes are 4.3 x 1025 to 14.9 x 1025 dyne cm/yr.
|
Geodetic Moment Rates
The range of moment rates, determined previously, can be compared with moment
rates that can be estimated from the geodetic deformation rates. To do this we
need models that relate the deformation rates to moment rates, which is
nonunique. Acknowledging the nonuniqueness and uncertainty involved with
converting surface strain to a scalar moment rate, this study utilizes four
methods (Anderson, 1979;
Ward, 1998a,
1998b;
Working Group on California Earthquake Probabilities, 1995;
Savage and Simpson, 1997) to
help quantify the moment rate from geodesy and its associated errors. These are
listed in Table 5.
Table 5 shows that for uniaxial
strain (e.g., in the x2 direction) the four methods are very
similar. Anderson (1979) is
different from the other three because he proposed an adjustment for inefficient
fault orientations (parameter k), whereas the other three methods
estimate the minimum rate.
|
The moment rate estimate from geodetic strain rates is proportional to the chosen seismogenic thickness. To choose the most appropriate thickness, depth distributions of microearthquakes from catalog listings from Nevada and Utah were utilized. Cumulative frequency plots of depth distribution of microearthquakes (Fig. 11) from general catalogs show that 98% of events occur at depths less than 15 km for the entire Utah (1962–1999) region and 17 km for the Nevada region (1990–1999). Time intervals for data presented in Figure 11 were selected based on station coverage. Previous studies have investigated the seismogenic depth for Utah in more detail. Based on well-constrained focal depths, 15 km is the preferred source depth used by Wong et al. (2001) for Utah. Arabasz et al. (1992) found a north–south dependence in focal depth distributions along the Wasatch Front with 99% of events occurring at 17 km in the south, to 11 km in the north. Considering the preceding studies and Figure 11, we assume all deformation occurs seismically above an average brittle-ductile transition depth of W = 15 km for the entire study region. The uncertainty introduced by this assumption is probably under 20%.
|
We predict the moment rate for the Basin and Range province from the regional geodetic strain rate tensor models introduced earlier (Table 5), using a shear modulus of µ = 3 x 1011 dyne cm2, W = 15 km, and k = 0.75 (Anderson, 1979). Resulting moments from geodesy are in the range from 3.91 x 1025 to 6.93 x 1025 dyne cm/yr.
Geological Moment Rates
To determine the geological moment release rate we utilize fault parameters
used as input to the 1996 and 2002 USGS seismic hazard maps
(Frankel et al., 1996,
2002;
Haller et al., 2002).
Data for California come from the 1996 model, whereas all other data for the
study region come from the 2002 model
(Haller et al., 2002).
Although improvements in geological fault characterizations are ongoing, the
USGS database represents an important landmark and plays a key role
in national building codes. The resulting geological moment rate for the region
is 2.55 x 1025 dyne cm/yr.
Comparison of Moment Rates
Moment rates from the seismic, geodetic, and geological methods are compared
in Table 6 and
Figure 12. As discussed
previously these moment rates are intended to display the full range of
uncertainties. Thus we believe that the results in
Figure 12 are a robust
comparison of these differing techniques which we later use to infer the
long-term seismicity rate. Within uncertainties, seismic and geodetic rates are
in agreement. Geological rates are much lower than the seismicity and geodetic
rates. This is not surprising considering the limiting paleoseismic data, which
are necessarily based only on faults that have been well characterized, a
minority of all the faults in the Basin and Range.
|
|
| Maximum Magnitude for the Great Basin |
|---|
|
|
|---|
Anderson and Luco (1983) related three functional forms of the Gutenburg-Richter curves to the moment rate. These models depend on (1) the rate of occurrence at a reference magnitude, (2) a b-value, and (3) the maximum magnitude Mmax. Major differences among the models occur in the way they are truncated as they approach Mmax. The integral of the area under these curves defines the moment release rate, based on the specified frequency-magnitude relationship. Note that these models are not physical laws to which earthquake statistics must ultimately conform. However, assuming that the frequency-magnitude distribution characterized by the seismic catalog is representative of the region, we apply these curves to the current catalog (1) to determine which of the Anderson and Luco (1983) models is most representative of the Basin and Range seismicity and (2) to determine what values of Mmax come out of fitting each of these curves to the observations.
We compare the historical cumulative magnitude distribution (Fig. 4b) with the Anderson and Luco (1983) models (Fig. 13). The curves are normalized using a moment release rate of 11.3 x 1025 dyne cm/yr, corresponding to the upper limit of our most likely historical moment rate range obtained from the historical seismicity. Increasing Mmax, while continuing to match a moment release rate of 11.3 x 1025 dyne cm/yr, results in lowering of the curves. The comparison of these models to the historical b-value curve indicates that model 1 best matches the data when Mmax 7.58, model 2 matches when Mmax 8.0, and model 3 matches when Mmax 8.2. Lowering the moment release rate to 5.8 x 1025 dyne cm/yr lowers the curve requiring Mmax 6.8 for model 1 to satisfy the historical seismicity, whereas models 2 and 3 fail to fit the high-magnitude end of the curve.
|
Models 2 and 3 appear to match the cumulative moment curve better than model 1 in Figure 13. However, at low- occurrence rates the shape of the curve determined by data is poorly constrained given the duration of observations. We feel that there is not sufficient data to determine which model best characterizes Basin and Range seismicity rates.
With consideration of moment rates from historical seismicity and geodesy,
the Anderson and Luco (1983)
models imply that there is no reason to expect an earthquake in the Great Basin
with magnitude greater than Mmax
8.2. This is
inconsistent with the suggestions of
Wernicke (1995) and
Kagan (1999), through separate
lines of research, that the region could experience earthquakes of magnitudes
much greater than eight.
Kagan (1999) proposed a
universal magnitude-frequency distribution in which the b-values and
maximum moment is the same for all continental regions. He proposed that for
shallow earthquakes the universal value of the effective maximum moment
magnitude is of the order 8.5–9.0. As shown previously, this magnitude is
inconsistent with the historical earthquake rates. An earthquake of this size
also runs into difficulty on physical grounds. Consider for instance that an
event with Mw 8.5 has, by
equation (3), a seismic moment of
7 x 1028 dyne cm. To maximize the fault area and shear modulus,
we consider values larger than those used previously: µ = 4
x 1011 dyne cm2 and a dip of 45° where the
seismogenic thickness is 17 km, yielding W = 24 km. From the
definition of seismic moment
(equation 1), the product of
fault length (L) and mean slip (D) for an event of this size
will be about LD = 7300 m km. A rupture on a fault of 300 km
length (greater than any fault in the Great Basin) would need an average slip of
24 m, which is larger than any observed historical rupture anyplace in the
world. We thus conclude that the Kagan model is not reasonable, both on
statistical and physical grounds.
Wernicke (1995) speculates that seismogenic low-angle normal faults that form the base of the entire seismogenic zone play an important role in accommodating Basin and Range extension. He further speculates that these faults have longer recurrence intervals than steeply dipping faults because they fail in infrequent, extremely large magnitude events. Wernicke suggests that the hypothesized ability of low-angle normal faults at the base of the crust to generate large events with greater magnitudes has not been tested globally because of the short historical records. Although low-angle normal faults have been imaged in the Basin and Range province (Allmendinger, 1983; Abbott et al., 2001), there is no evidence there or elsewhere that such structures experience brittle failure. A common expectation is that the extension of the region occurs by creep on these faults below the brittle-ductile transition. If the Wernicke speculation is correct, then the stress accumulation would roughly double the seismic-moment rate from the estimates given previously, and the statistical argument would allow Mmax to increase to about 8.2. However, evidence for accommodating extension by a creep mechanism (Buck et al., 2003) at the strain rates present in the Great Basin suggests that this model is unlikely.
| Discussion and Conclusions |
|---|
|
|
|---|
Uncertainties in estimates of moment rates do not arise only on the seismicity side. Geodetic and seismic measurements sample different aspects of the deformation field. Seismicity and geological estimates serve only as a record of brittle deformation, whereas geodesy encompasses both seismic and aseismic strain accumulation. Recently, significant aseismic deformation has been observed below Lake Tahoe, California (Smith et al., 2004), demonstrating the importance of aseismic strain in the Great Basin. Geodetic rates cannot uniquely determine slip at depth (Savage and Simpson, 1997) and may only give a measure of the instantaneous strain transients, which may not be preserved throughout the earthquake cycle (Savage and Lisowski, 1998; Shen-Tu et al., 1999). Recent geodetic observations have shown that historical moment release has occurred where contemporary strain is accumulating (Hammond and Thatcher, 2004). Whether the geodetic signature is representative of long- term deformation is argued by Thatcher (1995) to depend on fault characteristics and spacing, the extent of the cyclic zone, and whether the geodetic network completely spans the entire deforming zone. Geodetic measurements may therefore be sensitive to the duration and sampling of deformation with agreement also dependent on whether the sample period covered includes earthquakes that contributed significant post and coseismic displacements to the observed geodetic velocities. Apart from these uncertainties due to Earth processes, additional uncertainties are introduced in the processing decisions in converting individual GPS measurements to regional strain-rate models. These uncertainties become more important as the size of the region decreases. The final geodetic uncertainty arises in converting strain rates to seismic moment rates.
We suggest that it is possible to expand on and quantify the suggestions of
Smith (1976) and
Ward (1998a) that the adequacy
of seismic catalogs to estimate earthquake recurrence rates depends on the area
of the region, catalog duration, and regional strain rates. Where Smith looked
at the product of time and length of the seismic zone, and Ward looked at the
product of time and strain rate, we suggest that it is useful to define a
catalog adequacy parameter
|
| (4) |
= 7.25 x 105 km2,
and, for the region as a whole
V/L = 1.3 x 10–8/yr, where the
slip rate V = 13 mm/yr is the relative velocity of the Sierra
Nevada block relative to stable North America, and L = 1000 km
is the average width of the region. These parameters yield Z
1.38
km2. For the domains A–D
(Fig. 5), using the same
procedure, Z
0.1–0.25, and since within these areas the
historical and geodetic methods do not agree well
(Fig. 7), these values of
Z are apparently too small. Based on these Basin and Range results, it
seems reasonable to expect that in other regions with Z
1.5
km2, historical seismicity and geodesy will agree within
uncertainties, although of course more testing is needed to confirm this
hypothesis. The spatial consistency of the distribution of small earthquakes, deformation, and moment release shown in Figure 7 is interesting. Given the shortness of the catalog duration and variability in absolute rates, spatial similarity between geodetic deformation and seismicity is not necessarily expected. Even in domain B (Fig. 7b), where the moment release is highly concentrated in the Owens Valley because of the 1872 MW 7.74 Owens Valley earthquake, the spatial distribution of moderate earthquakes follows a broader curve similar to the deformation. Considering the uncertainties in estimates of the spatial distribution of the seismic hazard from historical seismicity, geodesy, and geology, consistency of more than one of these techniques is a key factor in providing confidence for how seismic hazards are localized. We suggest that it is worthwhile to investigate conditions for the similarity of seismicity and geodesy to hold and consider that it is possible that some criteria based on a catalog adequacy, or related, parameter might be possible.
| Appendix 1: Magnitudes of the Ten Largest Events |
|---|
|
|
|---|
Pleasant Valley 1915—10 3 6 53
Doser (1988) gives a
body-wave estimate of the earthquake size of MW 6.82.
However, Doser and Smith (1989)
model this event as two subevents and obtain a best estimate of
MW 6.89. The lower bound on this best estimate of
MW 6.87 gives us the minimum estimate of the size of
this event. A number of geological estimates were also available.
Doser (1988) and
Doser and Smith (1989) quote a
geological MW 7.15 from geological data using average
displacements, a focal depth of 16 km and a dip of 60°.
Wells and Coppersmith (1994)
also obtain a geological estimate of MW 7.15, whereas our
calculations using the Wells and Coppersmith data give an estimate of
MW 7.14.
Because of the low magnitude of the body-wave estimate we retain the Wells and Coppersmith (1994) value of MW 7.15 for the final catalog, which is also the maximum estimate.
Cedar Mountain 1932—12 21 6 10
Only body-wave and geological estimates were available for this event.
Doser (1986,
1988) and
Doser and Smith (1989) modeled
this earthquake as two subevents. We choose the moment calculated by summing the
individual moments of the two subevents as the best estimate of the body wave
estimate of MW 6.80. The lower bound of this estimate of
MW 6.73 gives us our minimum estimate.
Geological estimates were also available. Doser (1988) and Doser and Smith (1989) quote a geological estimate from Slemmons (1957) of MW 7.09. Using data from Wells and Coppersmith (1994), we obtain a value of MW 7.10, which was the maximum estimate for this event. We prefer the geological estimates over the lower body-wave estimate obtained by Doser and Smith because the geological data are more representative of the extent of faulting than the short- wavelength body-wave data.
Stillwater 1954—8 24 5 51
An average body-wave modeling estimate of MW 6.50 is
given by Doser (1986,
1987) and
Doser and Smith (1989),
combining three subevents.
Doser (1987) estimated
MW 6.55 from Pnl.
Barker and Doser (1988) invert
for the seismic moment using teleseismic body waves, and Pnl waves
separately, as well as both waveform types combined.
Doser (1986,
1987) also gives estimates using
body wave and Pnl waves, modeling the event as two subevents. A
combined teleseismic-Pnl determination from inversion of both waveform
types giving MW 6.27
(Barker and Doser, 1988) is
preferred compared with separate inversions. Geological estimates were also
given by Mason (1996),
Doser (1987), and
Doser and Smith (1989) and from
the
Wells and Coppersmith (1994)
data of MW 6.76, MW 6.41,
MW 6.41, and MW 6.52, respectively.
We prefer to use a value of MW 6.76, obtained using best estimates of geological parameters from Mason, whereas our maximum estimate of MW 6.77 comes from upper-bound estimates of these values. The minimum bound on the moment magnitude obtained from the literature was from Pnl modeling of MW 6.16 (Barker and Doser, 1988).
Fairview Peak 1954—12 16 11 7
There are no fewer than 40 estimates of the magnitude for this earthquake.
Two surface-wave estimates of MW 7.35 and
MW 7.12, using different models of the focal mechanism, were
given by
Doser and Kanamori (1987), with
two different focal mechanisms for both the Fairview Peak and Dixie Valley
events. Doser and Smith (1989)
quote the MW 7.12 as their best surface-wave estimate. We
also adopt this value as our preferred surface-wave magnitude estimate for this
earthquake.
The largest estimated moment magnitude from the literature was from geological estimates obtained from Mason (1996) of MW 7.36. Doser (1986) estimates MW 6.91 from body-wave modeling of two subevents, but the later study by Doser and Smith (1989) increases the minimum estimate to MW 6.95, which is what we use for the minimum estimate.
Dixie Valley 1954—12 16 11 11
A surface-wave estimate of MW 7.06 was given by
Doser and Smith (1989).
Doser and Kanamori (1987) gave
surface-wave estimates of MW 7.15 and MW
6.92 according to the two focal mechanisms that they employ for both the
Fairview Peak and Dixie Valley events. In consistency with the Fairview Peak
event, we choose to retain the best estimate of
Doser and Smith (1989)
MW 7.06.
Body-wave estimates were given by Doser (1986) and Doser and Smith (1989). They both model the event as two subevents. Geological estimates of the moment magnitude were available from Caskey et al. (1996), Mason (1996) and Wells and Coppersmith (1994). Doser and Kanamori (1987) quoted the geological estimate from Slemmons (1957) and Savage and Hastie (1969). Doser and Smith (1989) also quoted an estimate of the geological moment magnitude from the literature. From body-wave modeling, the minimum estimate of MW 6.53 was estimated by Doser and Smith (1989), while the maximum of MW 7.24 comes from geological estimates (Caskey et al., 1996).
Hebgen Lake 1959—8 18 6 37
Moment magnitude estimates were available from surface-wave and body-wave
modeling as well as from geodetic and geological observations.
Doser (1985) modeled this event as two subevents, with the second subevent 5 sec after the first. Doser and Smith (1989) provided a surface-wave estimate of MW 7.32 combining two subevents, whereas Doser and Kanamori (1987) calculated a surface-wave estimate of MW 7.39. We retained the lowest surface-wave estimate of MW 7.32 from Doser and Smith (1989) as our best estimate from the literature.
Geological data from Savage and Hastie (1969) are used to derive an MW 7.42. This was the largest value of the moment magnitude obtained from the literature. A minimum value of MW 7.01 was obtained from Doser and Smith (1989) from geological estimates, which includes deformation for the whole sequence. We therefore prefer to use a minimum magnitude of MW 7.25 derived from body-wave modeling, which included two subevents (Doser, 1985; Doser and Smith, 1989).
Borah Peak 1983—10 28 14 6
Harvard CMT estimate of MW 6.93 available for
this event was retained as the final estimate of the earthquake size according
to our selection criteria.
Moment tensor estimates of the moment magnitude were obtained from Ekstrom and Dziewonski (1985). Body- wave CMT estimates were also obtained from Nabelek et al. (1985). Estimates for geodetic leveling data were obtained from Barrientos et al. (1985), Stein and Barrientos (1985), Ward and Barrientos (1986), and Doser and Smith (1989). Body-wave modeling estimates were provided from Barker and Doser (1988) and Doser and Smith (1989). Geological estimates were obtained from Wells and Coppersmith (1994), Doser and Smith (1989), and Mason (1996).
Geological parameters documented by Mason (1996) give a maximum estimate of MW 7.25. This was the largest estimate obtained from the literature. Both Ekstrom and Dziewonski (1985) and Nabelek et al. (1985) give an estimate of MW 6.93 from body-wave centroid moment tensor solutions and body-wave inversions, respectively. A minimum estimate of MW 6.82 is documented by Doser and Smith (1989) from body-wave modeling.
Landers 1992—6 28 11 57
From the data provided by
Wells and Coppersmith (1994),
we calculate a minimum of MW 7.19. The maximum estimate is
from
Kanamori et al. (1992)
who give an estimate of MW 7.30 from surface-wave inversion.
However, a Harvard CMT moment magnitude estimate was available. As
this estimate of MW 7.29 is favored above all others, this
value of the earthquake size was retained in the final catalog.
Hector Mine 1999—10 16 9 46
Harvard CMT estimate of MW 7.12 was available
for this event and retained as the final estimate of the earthquake size
according to our selection criteria. This was the maximum estimate for this
event a minimum estimate of MW 7.10 is given by the Berkeley
(BK) and Council of the National Seismic System (CNSS)
catalogs. The maximum estimate of MW 7.14 is from
Ji et al. (2002).
| Appendix |
|---|
|
|
|---|
Manuscript received August 18, 2004
Abbott, R., J. N. Louie, S. J. Caskey, and S. Pullammanappallil (2001). Geophysical confirmation of low-angle normal slip on the historically active Dixie Valley Fault, Nevada, J. Geophys. Res. 106,4169 –4181.[CrossRef]
Allmendinger, R. W., J. W. Sharp, D. Von Tish, L. Serpa, L. Brown, S. Kaufman, J. E. Oliver, and R. B. Smith
(1983). Cenozoic and Mesozoic structure of the eastern Basin and
Range Province, Utah, from COCORP seismic-reflection data,
Geology 11,532
–536.
Anderson, J. G. (1979). Estimating the seismicity from
geological structure for seismic-risk studies, Bull. Seism. Soc.
Am. 69,135
–158.
Anderson, J. G. (2003). Strong-motion seismology, in International Handbook of Earthquake and Engineering Seismology, W. H. K. Lee, H. Kanamori, P. C. Jennings, and C. Kisslinger (Editors), Vol.81B , Academic Press, Amsterdam,937 –905.[CrossRef]
Anderson, J. G., and J. E. Luco (1983). Consequences of
slip rate constraints on earthquake occurrence relations, Bull.
Seism. Soc. Am. 73,471
–496.
Arabasz, W. J., J. C. Pechmann, and E. D. Brown (1992). Observational seismology and the evaluation of earthquake hazards and risk in the Wasatch Front area, Utah, Assessment of Regional Earthquake Hazards and Risk along the Wasatch Front, Utah, P. L. Gori and W. W. Hays (Editors), in U.S. Geol. Surv. Profess. Pap., 1500-A-J, D1–D36.
Barker, J. S., and T. C. Wallace (1986). A note on the
teleseismic body waves from the 23 November 1984 Round Valley, California,
earthquake, Bull. Seism. Soc. Am.76
,883
–888.
Barker, M. R., and D. I. Doser (1988). Joint inversion of regional and teleseismic earthquake waveforms, J. Geophys. Res. 93,2037 –2045.
Barrientos, S. E., S. N. Ward, J. Gonzalez-Ruiz, and R. S. Stein (1985). Inversion for moment as a function of depth from geodetic observations and long period body waves of the Borah Peak, Idaho, earthquake, Workshop XXVIII on the Borah Peak earthquake, U.S. Geol. Surv. Open-File Rept. 85-290,485 –518.
Beanland, S., and M. M. Clark (1993). Late Quaternary history of the Owens Valley fault zone, eastern California and surface rupture associated with the 1872 earthquakes (abstracts with programs), Geol. Soc. Am. 25,7 .
Bennett, R. A., B. P. Wernicke, N. A. Niemi, and M. Friedrich (2003). Contemporary strain rates in the northern Basin and Range province from GPS data, Tectonics22 , 1008, doi10.1029/2001TC1355 .
Boatwright, J., and G. L. Choy (1985). Teleseismic estimates of the energy radiated by shallow earthquakes, U.S. Geol. Surv. Open-File Rept. 85-290,409 –448.
Bolt, B. A., and R. D. Miller (1975). Catalogue of earthquakes in Northern California and adjoining areas, 1 January 1910–31 December 1972, Seismographic Stations, University of California, Berkeley, California.
Brune, J. N. (1968). Seismic moment, seismicity, and slip rate along major fault zones, J. Geophys. Res.73 ,777 –784.[CrossRef]
Buck, W. R., L. L. Lavier, and A. Babeyko (2003). A numerical model of lithospheric extension producing fault-bounded basins and ranges, Int. Geol. Rev.45 ,712 –723.
Caskey, S. J., S. G. Wesnousky, P. Zhang, and D. B. Slemmons
(1996). Surface faulting of the 1954 Fairview Peak (Ms 7.2) and
Dixie Valley (Ms 6.8) Earthquakes, Central Nevada, Bull. Seism. Soc.
Am. 86,761
– 787.
Catchings, R. D., and W. D. Mooney (1991). Basin and Range and upper mantle structure northwest to central Nevada, J. Geophys. Res. 96,6247 –6267.
Chulick, G. S., and W. D. Mooney (2002). Seismic
structure of the crust and uppermost mantle of North America and adjacent
oceanic basins: a synthesis, Bull. Seism. Soc. Am.92
,2478
–2492.
dePolo, C. M. (1998). A reconnaissance technique for estimating the slip rates of normal-slip faults in the Great Basin, and application to faults in Nevada, U.S.A., Ph.D. Dissertation, University of Nevada, Reno,381 p.
dePolo, D., and C. M. dePolo (1999). Earthquakes in Nevada, 1852–1998, Nevada Bureau of Mines and Geology, Map 119.
dePolo, C. M., A. R. Ramelli, R. H. Hess, and J. G. Anderson (2003). Reevaluation of pre-1900 earthquakes in Western Nevada, NBMG Open-File Report, 2003-3, 208 pp.
Dixon, T. H., M. Miller, F. Farina, H. Wang, and D. Johnson (2000). Present-day motion of the Sierra Nevada block and some tectonic implications for the Basin and Range Province, North American Cordillera, Tectonics19 ,1 –24.[CrossRef][Web of Science][GeoRef]
Dixon, T. H., E. Norabuena, and L. Hotaling (2003).
Paleoseismology and Global Positioning System; earthquake-cycle effects and
geodetic versus geologic fault slip rates in the Eastern California shear zone,
Geology 31,55
–58.
Dixon, T. H., S. Robaudo, J. Lee, and M. Reheis, (1995). Constraints on present-day Basin and Range deformation from space geodesy, Tectonics 14,755 –772.[CrossRef][Web of Science][GeoRef]
Dokka, R. K., and C. J. Travis (1990). Role of the Eastern California Shear Zone in accommodating Pacific-North American plate motion, Geophys. Res. Lett.17 ,1323 –1326.[Web of Science][GeoRef]
Doser, D. I. (1985). Source parameters and faulting processes of the 1959 Hebgen Lake, Montana, earthquake sequence, J. Geophys. Res. 90,4537 –4555.[CrossRef]
Doser, D. I. (1986). Earthquake processes in the Rainbow Mountain– Fairview Peak–Dixie Valley, Nevada, region 1954–1959, J. Geophys. Res.91 ,12,572 –12,586.
Doser, D. I. (1987). Modeling the Pnl waveforms of the Fairview Peak– Dixie Valley, Nevada, U.S.A. earthquake sequence (1954–1959), Phys. Earth Planet. Interiors48 ,64 –72.[CrossRef]
Doser, D. I. (1988). Source parameters of earthquakes in the Nevada seismic zone, 1915–1943, J. Geophys. Res. 93,15,001 –15,015.
Doser, D. I. (1989a). Extensional tectonics in northern Utah–southern Idaho, U.S.A., and the 1934 Hansel Valley sequence, Phys. Earth Planet. Interiors54 ,120 –134.[CrossRef]
Doser, D. I. (1989b). Source parameters of Montana
earthquakes (1925– 1964) and tectonic deformation in the northern
Intermountain seismic belt, Bull. Seism. Soc. Am.79
,31
–50.
Doser, D. I., and H. Kanamori (1987). Long-period
surface waves of four western United States earthquakes recorded by the Pasadena
strainmeter, Bull. Seism. Soc. Am.77
,236
–243.
Doser, D. I., and R. B. Smith (1982). Seismic moment
rates in the Utah region, Bull. Seism. Soc. Am.72
,525
–551.
Doser, D. I., and R. B. Smith (1985). Source parameters
of the 28 October 1983 Borah Peak, Idaho, Earthquake from Body Wave Analysis,
Bull. Seism. Soc. Am.75
,1041
–1051.
Doser, D. I., and R. B. Smith (1989). An assessment of
source parameters of earthquakes in the Cordillera of the western United States,
Bull. Seism. Soc. Am.79
,1383
–1409.
Ekstrom, G., and A. M. Dziewonski (1985).
Centroid-moment tensor solutions for 35 earthquakes in western North America
(1977–1983), Bull. Seism. Soc. Am.75
,23
–40.
Flesch, L. M., W. E. Holt, A. J. Haines, and B. Shen-Tu
(2000). Dynamics of the Pacific-North America plate boundary in the
Western United States, Science287
,834
–836.
Frankel, A., C. Mueller, T. Barnhard, D. Perkins, E. Leyendecker, N. Dickman, S. Hanson, and M. Hopper. (1996). National seismic-hazard maps: documentation June 1996, U.S. Geol. Surv. Open-File Rept. 96-532,110 p.
Frankel, A. D., M. D. Petersen, C. S. Mueller, K. M. Haller, R. L. Wheeler, E. V. Leyendecker, R. L. Wesson, S. C. Harmsen, C. H. Cramer, D. M. Perkins, and K. S. Rukstales (2002). Documentation for the 2002 update of the national seismic hazard maps, U.S. Geol. Surv. Open-File Rept. 02-420.
Freymueller, J. T., M. H. Murray, P. Segall, and D. Castillo (1999). Kinematics of the Pacific North America plate boundary zone, northern California, J. Geophys. Res.104 ,7419 –7441.[CrossRef]
Haines, A. J., and W. E. Holt (1993). A procedure for obtaining the complete horizontal motions within zones of distributed deformation from the inversion of strain rate data, J. Geophys. Res. 98,12,057 –12,082.
Haller, K. M., R. L. Wheerler, and K. S. Rukstales (2002). Document of changes in fault parameters for the 2002 National Seismic Hazard Maps–Conterminous United States except California, U.S. Geol. Surv. Open-File Rept. 02-467.
Hammond, W. C., and W. Thatcher (2003). Crustal deformation of the northern Basin and Range from measurements with the Global Positioning System (abstract), EOS Trans. AGU84 , no. 46, Fall Meet. Suppl.,G32B-07 .
Hammond, W. C., and W. Thatcher (2004), Contemporary tectonic deformation of the Basin and Range province, western United States: 10 years of observation with the Global Positioning System, J. Geophys. Res. 109,B08403 , doi10.1029/2003JB002746 .[CrossRef]
Hammond, W. C., W. Thatcher, and G. Blewitt (2004). Crustal deformation across the Sierra Nevada-northern Walker Lane, Basin and Range transition, western United States Measured with GPS, 2000–2004 (abstract), EOS Trans., AGU85 , no.47, Fall Meet. Suppl.,G31D-07 .
Hanks, T. C., J. A. Hileman, and W. Thatcher (1975).
Seismic moments of the larger earthquakes of the Southern California Region,
Geol. Soc. Am. Bull.86
,1131
–1139.
Hanks, T. C., and H. Kanamori (1978). A moment magnitude scale, J. Geophys. Res.84 ,2348 –2350.
Hart, R. S., R. Butler, and H. Kanamori (1977).
Surface-wave constraints on the August 1, 1975, Oroville earthquake,
Bull. Seism. Soc. Am.67
, no. 1,1
–7.
Hauksson, E., L. M. Jones, K. Hutton, and D. Eberhart-Phillips (1993). The 1992 Landers earthquake sequence: seismological observations, J. Geophys. Res.98 ,19,835 –19,858.
Holt, W. E., B. Shen-Tu, J. Haines, and J. Jackson (2000). On the determination of self-consistent strain rate fields within zones of distributed deformation, in The History and Dynamics of Global Plate Motions, M. A. Richards, R. G. Gordon, and R. D. van der Hilst, (Editors), American Geophysical Union, Washington, D.C.,113 –141.
Humphreys, E. D., and R. J. Weldon (1994). Deformation across the western United States: a local estimate of Pacific-North American transform deformation, J. Geophys. Res.99 ,19,975 –20,010.[CrossRef]
Ichinose, G. A., J. G. Anderson, K. D. Smith, and Y. Zeng
(2003). Source parameters of Eastern California and Western Nevada
Earthquakes from regional moment tensor inversion, Bull. Seism. Soc.
Am. 93,61
–84.
Ji, C., D. J. Wald, and D. V. Helmbereger (2002). Source
description of the 1999 Hector Mine, California, Earthquake, Part II: Complexity
of slip history, Bull. Seism. Soc. Am.92
,1208
–1226.
Jones, C. H., B. P. Wernicke, G. L. Farmer, J. D. Walker, D. S. Coleman, L. W. McKenna, and F. V. Perry (1992). Variations across and along a major continental rift: An interdisciplinary study of the Basin and Range Province, Western USA, Tectonophysics 213,57 –96.[CrossRef][Web of Science][GeoRef]
Kagan, Y. Y. (1999). Universality of the seismic moment-frequency relation, Pure Appl. Geophys.155 ,537 –573.[CrossRef]
Kanamori, H. (1977). The energy release in great earthquakes, J. Geophys. Res.82 ,2981 –2987.[CrossRef]
Kanamori, H. (1983). Magnitude scale and quantification of earthquakes, Tectonophysics93 ,185 –199.[CrossRef][Web of Science][GeoRef]
Kanamori, H., H.-K. Thio, D. Dreger, and E. Hauksson (1992). Initial investigation of the Landers, California, earthquake of 28 June 1992 using TERRAscope, Geophys. Res. Lett. 19,2267 –2270.[Web of Science][GeoRef]
Kostrov, V. V. (1974). Seismic moment and energy of earthquakes, and seismic flow of rock, Earth Phys.1 ,23 –40.
Kreemer, C., W. E. Holt, S. Goes, and R. Govers (2000). Active deformation in eastern Indonesia and the Philippines from GPS and seismicity data, J. Geophys. Res.105 ,663 –680.[CrossRef]
Kreemer, C., W. E. Holt, and A. J. Haines (2002). The global moment rate distribution within plate boundary zones, in Plate Boundary Zones Geodynamics Series 30, S. Stein and J. T. Freymueller (Editors), doi 10/1029/030GD10, American Geophysical Union, Washington, D.C.,173 –190.
Langston, C. A., and R. Butler (1976). Focal mechanism
of the August 1, 1975, Oroville earthquake, Bull. Seism. Soc.
Am. 66, no. 4,1111
– 1120.
Lavallée, D., G. Blewitt, P. J. Clarke, K. Nurutdinov, W. E. Holt, C. Kreemer, C. M. Meertens, W. S. Shiver, S. Stein, S. Zerbini, L. Bastos, and H.-G. Kahle (2001). GPSVEL project: towards a dense global GPS velocity field, AG Scientific Assembly, Budapest.
Lee, J., J. Spencer, and L. Owen (2001). Holocene slip
rates along the Owens Valley fault, California; implications for the recent
evolution of the Eastern California shear zone, Geology29
,819
–822.
Malservisi, R., T. H. Dixon, P. C. La Femina, and K. P. Furlong (2003). Holocene slip rate of the Wasatch fault zone, Utah, from geodetic data: earthquake cycle effects, Geophys. Res. Lett. 30, no. 13,1673 , doi10.1029/2003GL0174078 .[CrossRef]
Mason, D. B. (1996). Earthquake magnitude potential of
the Intermountain Seismic Belt, USA, from surface-parameter scaling of Late
Quaternary faults, Bull. Seism. Soc. Am.86
,1487
–1506.
Masson, F., J. Chery, D. Hatzfeld, J. Martinod, P. Vernant, F. Tavakoli, and M. Ghafory-Ashtiani (2005). Seismic versus aseismic deformation in Iran inferred from earthquakes and geodetic data, Geophys. J. Int.160 ,217 –226.[CrossRef]
Mazzotti, S., H. Dragert, J. Henton, M. Schmidt, R. Hyndman, T. S. James, Y. Lu, and M. Craymer (2003). Current tectonics of northern Cascadia from a decade of GPS measurements, J. Geophys. Res.108 , 2554, doi2510.1029/2003JB002653 .[CrossRef]
McCaplin, J. P., and S. P. Nishenko (1996). Holocene paleoseismicity, temporal clustering, and probabilities of future large (M > 7) earthquakes on the Wasatch fault zone, Utah, J. Geophys. Res. 101,6233 – 6253.[CrossRef]
McClusky, S. C., S. C. Bjornstad, B. H. Hager, R. W. King, B. J. Meade, M. M. Miller, F. C. Monastero, and B. J. Souter (2001). Present day kinematics of the Eastern California Shear Zone from a geodetically constrained block model, Geophys. Res. Lett. 28,3369 –3372.[CrossRef][Web of Science][GeoRef]
Molnar, P. (1979). Earthquake recurrence intervals and
plate tectonics, Bull. Seism. Soc. Am.69
,115
–133.
Nabelek, J. L., H. Eyidogan, and M. N. Tiksoz (1985). Source parameters of the Borah Peak, Idaho, earthquake of October 28, 1983 from body- wave inversion, EOS Trans. AGU66 ,308 .
Oldow, J. S., C. L. V. Aiken, J. L. Hare, J. F. Ferguson, and R. F. Hardyman
(2001). Active displacement transfer and differential block motion
within the central Walker Lane, western Great Basin,
Geology 29,19
–22.
Pacheco, J., and J. Nábelek
(1988). Source mechanisms of three moderate Californian earthquakes
of July 1986, Bull. Seism. Soc. Am.78
, 1907–
1929.
Patton, H. J. (1985). P-wave fault plane solutions and the generation of surface wave by earthquakes in the Western United States, Geophys. Res. Lett.112 ,518 –521.
Patton, H. J., and D. I. Doser (1988). Inversion of regional Pnl and surface- wave data for source parameters of the Borah Peak aftershock, Geophys. Res. Lett.5 ,459 –462.
Rogers, A. M., S. C. Harmsen, E. J. Corbett, K. F. Priestly, and D. dePolo (1991). The seismicity of Nevada and some adjacent parts of the Great Basin, in Neotectonics of North America, D. B. Slemmons, E. R. Engdahl, M. D. Zoback, and D. D. Blackwell (Editors), Geological Society of America, Boulder, Colorado, Decade Map Volume1 , 153– 184.
Savage, J. C., and L. M. Hastie (1966). Deformation associated with dip- slip faulting, J. Geophys. Res.71 ,4897 –4904.
Savage, J. C., and L. M. Hastie (1969). A dislocation
model for the Fairview Peak, Nevada, earthquake, Bull. Seism. Soc.
Am. 59,1937
–1948.
Savage, J. C., and M. Lisowski (1998). Viscoelastic coupling model of the San Andreas Fault along the big bend, southern California, J. Geophys. Res.103 ,7281 –7292.[CrossRef]
Savage, J. C., and R. W. Simpson (1997). Surface strain
accumulation and the seismic moment tensor, Bull. Seism. Soc.
Am. 87,1345
–1353.
Savage, J. C., W. J. Gan, W. H. Prescott, and J. L. Svarc (2004). Strain accumulation across the Coast Ranges at the latitude of San Francisco, J. Geophys. Res.109 , 3413, doi3410.1029/2003JB002612 .[CrossRef]
Shen-Tu, B., W. E. Holt, and A. J. Haines (1998). Contemporary kinematics of the Western United States determined from earthquake moment tensors, very long baseline interferometry, and GPS observations, J. Geophys. Res.103 ,18,087 –18,117.[CrossRef]
Shen-Tu, B., W. E. Holt, and A. J. Haines (1999). Deformation kinematics in the western United States determined from Quaternary fault slip rates and geodetic data, J. Geophys. Res.104 ,28,927 –28,955.[CrossRef]
Shimazaki, K., and T. Nakata (1980). Time-predictable recurrence model for large earthquakes, Geophys. Res. Lett. 7,279 –282.[Web of Science][GeoRef]
Sipkin, S. A. (1986). Interpretation of non-double-couple earthquake mechanisms derived from moment tensor inversion, J. Geophys. Res.91 ,531 –547.
Slemmons, D. B. (1957). Geological effects of the Dixie
Valley–Fairview Peak, Nevada, earthquakes of December 16, 1954,
Bull. Seism. Soc. Am.47
,353
–375.
Slemmons, D. B., K. V. Steinbugge, D. Tocher, G. B. Oakshott, and V. P. Gianella
(1959). Wonder, Nevada, Earthquake of 1903, Bull. Seism.
Soc. Am., 49,251
–265.
Smith, S. W. (1976). Determination of maximum earthquake magnitude, Geophys. Res. Lett.3 ,351 –354.[Web of Science][GeoRef]
Smith, R. B., and W. J. Arabasz (1991). Seismicity of the intermountain seismic belt, in Neotectonics of North America, D. B. Slemmons, E. R. Engdahl, M. D. Zoback, and D. D. Blackwell (Editors), Geological Society of America, Boulder, Colorado, Decade Map Volume1 , 185– 228.
Smith, K. D., D. von Seggrenm, G. Blewitt, L. Preston, J. G Anderson, B. P. Wernicke, and J. L. Davis
(2004). Evidence for deep magma injection beneath Lake Tahoe,
Nevada-California, Science305
, 1277–
1280.
Sonder, L. J., and C. H. Jones (1999). Western United States extension; how the west was widened, Ann. Rev. Earth Planet. Sci., 27,417 – 462.[CrossRef][Web of Science]
Stein, R. S., and S. E. Barrientos (1985). Planar high-angle faulting in the Basin and Range: geodetic analysis of the 1983 Borah Peak, Idaho earthquake, J. Geophys. Res.90 ,11,355 –11,366.
Stepp, J. C. (1972). Analysis of completeness of the earthquake sample in the Puget Sound area and its effect on statistical estimates of earthquake hazard, Int. Conf. Microzon. Proc. 2,897 –910.
Stewart, J. H. (1988). Tectonic of the Walker Lane Belt, Western Great Basin; Mesozoic and Cenozoic deformation in a zone of shear, in Metamorphism and Crustal Evolution of the Western United States, W. G. Ernst (Editor), Vol. 7,683 –713.
Svarc, J. L., J. C. Savage, W. H. Prescott, and M. H. Murray (2002a). Strain accumulation and rotation in western Oregon and southwestern Washington, J. Geophys. Res.107 , 2087, 7, doi2010.1029/2001JB000625 .[CrossRef]
Svarc, J. L., J. C. Savage, W. H. Prescott, and A. R. Ramelli (2002b). Strain accumulation and rotation in western Nevada, 1993–2000, J. Geophys. Res.107 , 2090, 13, doi2010.1029/2001JB000579 .[CrossRef]
Thatcher, W. (1995). Microplate versus continuum descriptions of active tectonic deformation, J. Geophys. Res. 100,3885 –3894.[CrossRef][Web of Science]
Thatcher, W. (2003). GPS constraints on the kinematics of continental deformation, Int. Geol. Rev.45 ,191 –212.[CrossRef]
Thatcher, W., G. R. Foulger, B. R. Julian, J. Svarc, E. Quilty, and G. W. Bawden
(1999). Present-day deformation across the Basin and Range
Province, Western United States, Science283
,1714
–1718.
Toppozada, T. R., C. R. Real, and D. L. Parke (1981). Preparation on isoseismal maps and summaries of reported effects for pre-1900 California earthquakes, annual technical report, U.S. Geol. Surv. Open- File Rept. 81-11.
Tsai, Y-B., and K. Aki (1966). Source mechanism of the Truckee, California, earthquake of September 12, 1966, Bull. Seism. Soc. Am. 60,1199 –1208.
Utsu, T. (2003). Relationships between magnitude scales, in International Handbook of Earthquake and Engineering Seismology W. H. K. Lee, H. Kanamori, P. C. Jennings, and C. Kisslinger (Editors), Vol.81A , Academic Press, Amsterdam,733 –746.
Wallace, R. E. (1970). Earthquake recurrence intervals
on the San Andreas fault, Bull. Geol. Soc. Am.81
,2875
–2890.
Wallace, R. E. (1984). Patterns and timing of Late Quaternary faulting in the Great Basin Province and relation to some regional tectonic features, J. Geophys. Res.89 ,5763 –5769.[Web of Science]
Wang, S. C., K. C. McNally, and R. J. Geller (1982). Seismic strain release along the Middle America Trench, Geophys. Res. Lett. 9,182 –185.[Web of Science][GeoRef]
Ward, S. N. (1994). A multidisciplinary approach to
seismic hazard in Southern California, Bull. Seism. Soc.
Am. 84,1293
–1309.
Ward, S. N. (1998a). On the consistency of earthquake moment rates, geological fault data, and space geodetic strain: the United States, Geophys. J. Int.134 ,172 –186.[CrossRef]
Ward, S. N. (1998b). On the consistency of earthquake moment rates, geological fault data, and space geodetic strain: Europe, Geophys. J. Int.135 ,1011 –1018.[CrossRef]
Ward, S. N., and S. E. Barrientos (1986). An inversion for slip distribution and fault shape from geodetic observations of the 1983, Borah Peak, Idaho earthquake, J. Geophys. Res.91 ,4909 –4919.[CrossRef]
Wells, D. L., and K. J. Coppersmith (1994). New
empirical relationships among magnitude, rupture length, rupture width, rupture
area, and surface displacement, Bull. Seism. Soc. Am.84
,974
–1002.
Wernicke, B. (1992). Cenozoic extensional tectonics of the U.S. Cordillera in The Cordilleran Orogen: Conterminous U.S., The Geology of North America, Vol. G-3, P. W. Lipman, B. C. Burchfiel, and M. L. Zoback (Editors), Geological Society of America, Boulder, Colorado, 553– 581.
Wernicke, B. (1995). Low-angle normal faults and seismicity: a review, J. Geophys. Res.100 ,20,159 –20,174.[CrossRef]
Westaway, R., and R. B. Smith (1989), Source parameters
of the Cache Valley (Logan), Utah earthquake of 30 August 1962, Bull.
Seism. Soc. Am. 79, no. 5,1410
–1425.
Wong, I., W. Silva, S. Olig, P. Thomas, D. Wright, F. Ashland, N. Gregor, J. Peachman, M. Dober, G. Christenson, and R. Gerth (2001). Earthquake scenario and probabilistic ground shaking maps for the Salt Lake City, Utah, metropolitan area, Miscellaneous Publication 02-5, Utah Geological Survey, p 57.
Working Group on California Earthquake Probabilities
(1995). Seismic hazards in southern California: probable
earthquakes, 1994–2024, Bull. Seism. Soc. Am.85
,379
–439.
Zandt, G., S. C. Myers, and T. C. Wallace (1995). Crust and mantle structure across the Basin and Range–Colorado Plateau boundary at 37°N latitude and implications for Cenozoic extension mechanism, J. Geophys. Res.100 ,10,529 –10,548.[CrossRef][Web of Science]
This article has been cited by other articles:
![]() |
D. H. von Seggern, K. D. Smith, and L. A. Preston Seismic Spatial-Temporal Character and Effects of a Deep (25-30 km) Magma Intrusion below North Lake Tahoe, California-Nevada Bulletin of the Seismological Society of America, June 1, 2008; 98(3): 1508 - 1526. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Stein Approaches to continental intraplate earthquake issues Geological Society of America Special Papers, January 1, 2007; 425(0): 1 - 16. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. G. Anderson and Y. Miyata Ranking States by Seismic Activity Seismological Research Letters, November 1, 2006; 77(6): 672 - 676. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |