USGS Science⌕ Search

SEARCH · USGS Science

Results for “Bulletin of the Geological Society of America”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39Linked to original sources

Structure and seismic hazard of the Ventura Avenue anticline and Ventura fault, California: Prospect for large, multisegment ruptures in the Western Transverse Ranges

The Ventura Avenue anticline is one of the fastest uplifting structures in southern California, rising at ∼5 mm/yr. We use well data and seismic reflection profiles to show that the anticline is underlain by the Ventura fault, which extends to seismogenic depth. Fault offset increases with depth, implying that the Ventura Avenue anticline is a fault‐propagation fold. A decrease in the uplift rate since ∼30±10 ka is consistent with the Ventura fault breaking through to the surface at that time and implies that the fault has a recent dip‐slip rate of ∼4.4–6.9 mm/yr. To the west, the Ventura fault and fold trend continues offshore as the Pitas Point fault and its associated hanging wall anticline. The Ventura–Pitas Point fault appears to flatten at about 7.5 km depth to a detachment, called the Sisar decollement, then step down on a blind thrust fault to the north. Other regional faults, including the San Cayetano and Red Mountain faults, link with this system at depth. We suggest that below 7.5 km, these faults may form a nearly continuous surface, posing the threat of large, multisegment earthquakes. Holocene marine terraces on the Ventura Avenue anticline suggest that it grows in discrete events with 5–10 m of uplift, with the latest event having occurred ∼800 years ago ( Rockwell, 2011 ). Uplift this large would require large earthquakes ( M w 7.7–8.1) involving the entire Ventura/Pitas Point system and possibly more structures along strike, such as the San Cayetano fault. Because of the local geography and geology, such events would be associated with significant ground shaking amplification and regional tsunamis.

California↗

2001 Bhuj, India, earthquake engineering seismoscope recordings and Eastern North America ground-motion attenuation relations

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

Bulletin of the Seismological Society of America↗

Near real-time monitoring of seismic events and status of portable digital recorders using satellite telemetry

Near real-time monitoring of seismic events and status of portable 16-bit digital recorders has been established for arrays near Parkfield, Mammoth Lakes, and San Francisco, California. This monitoring system provides near real-time seismic event identification (rough location and magnitude) and a cost-effective means to maintain arrays at near 100% operational level. Principal objectives in the design of this system have been portability and low-cost telemetry. The system has been developed to use portable digital seismic recorders ( GEOS —General Earthquake Observation System) and portable data collection platforms ( DCP 's) for the Geostationary Operational Environmental Satellite (GEOS) telemetry system. Data are transferred asynchronously from the GEOS seismic system through a microprocessor-controlled interface every 10 min. The interface stores, determines priority, converts, and synchronously transfers these data to a Sutron Corp. model 8004 DCP for transmission through the GEOS satellite telemetry system. Event parameters include trigger time, peak amplitude, time of peak amplitude, and event duration. Instrument configuration parameters, transmitted at system start-up time and every 24 hr, include recording parameters, trigger parameters, GEOS software version, clock reference, and location parameter. Instrument status includes battery voltage, number of events, and percentage of tape usage. These data are transmitted as appropriate to the U.S. Geological Survey satellite downlink and computers located in Menlo Park, California, where they are processed and displayed.

Bulletin of the Seismological Society of America↗

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington↗

Holocene fault reactivation in the eastern Cascades, Washington

Significant uncertainty remains concerning how and where crustal shortening occurs throughout the eastern Cascade Range in Washington State. Using light detection and ranging (lidar) imagery, we identified an ∼ 5 ‐ km ‐ long "> ∼5‐km‐long lineament in Swakane canyon near Wenatchee, roughly coincident with a strand of the Entiat fault. Topographic profiles across the lineament reveal a southwest‐side‐up break in slope, with an average of 2–3 m of vertical separation of the hillslope surface. We consider a range of possible origins for this feature, including differential erosion across a fault‐line scarp, slope failure (sackung or landslide), and surface deformation across an active fault strand. Based on trenching, radiocarbon and luminescence dating, and ground‐penetrating radar (GPR) across the lineament, we conclude that warped saprolite observed in the shallow subsurface is most consistent with southwest‐side‐up folding caused by blind reverse faulting at depth. Following this reasoning, dating of overlying colluvial deposits suggests that at least one Holocene earthquake occurred on this strand of the southern Entiat fault, with an approximate vertical separation of ≥ 1    m "> m ≥1 m . GPR reveals up to 4 m of cumulative vertical separation of the saprolite, suggesting a history of multiple earthquakes on the structure. Taken in context with other potential fault‐related lineaments along the Entiat fault, our interpretation of Holocene earthquakes in Swakane canyon could suggest reactivation of longer sections of the Entiat fault, as well as of other bedrock faults in the eastern Cascades. Although active erosion and slow strain rates lead to a subdued geomorphic expression of recent deformation, we conclude that the reactivated Entiat fault represents a seismogenic structure that should be considered in regional seismic hazard analyses. The difficulty of recognizing low‐slip‐rate structures in forested and mountainous terrain underscores the importance of additional lidar surveys and geological and geophysical studies for fully understanding seismic hazard in regions with infrequent but potentially large earthquakes.

Washington↗

Relaxing segmentation on the Wasatch Fault Zone: Impact on seismic hazard

The multisegment Wasatch fault zone is a well-studied normal fault in the western United States that has paleoseismic evidence of recurrent Holocene surface-faulting earthquakes. Along the 270-km-long central part of the fault, four primary structural complexities provide possible along-strike limits to these ruptures and form the basis for models of fault segmentation. Here, we assess the impact that the Wasatch fault segmentation model has on seismic hazard by evaluating the time-independent long-term rate of ruptures on the fault that satisfy fault slip rates and paleoseismic event rates, adapting standard inverse theory used in the Uniform California Earthquake Rupture Forecast 3 (UCERF3), and implementing a segmentation constraint where ruptures across primary structural complexities are penalized. We define three models with varying degrees of rupture penalization: (1) segmented (ruptures confined to individual segments), (2) penalized (multi-segment ruptures allowed, but penalized), and (3) unsegmented (all ruptures allowed). Seismic-hazard results show that on average, hazard is highest for the segmented model, where seismic moment is accommodated by frequent moderate (moment magnitude, M w 6.2–6.8) earthquakes. The unsegmented model yields the lowest average seismic hazard because part of the seismic moment is accommodated by large (M w 6.9–7.9), but infrequent ruptures. We compare these results to model differences derived from other inputs such as slip rate and magnitude scaling relationships and conclude that segmentation exerts a primary control on seismic hazard. This study demonstrates the need for additional geologic constraints on rupture extent and methods by which these observations can be included in hazard-modeling efforts.

Utah↗

Late Holocene slip rate of the Mojave section of the San Andreas Fault near Palmdale, California

The geologic slip rate on the Mojave section of the San Andreas fault is poorly constrained, despite its importance for understanding earthquake hazard, apparent discrepancies between geologic and geodetic slip rates along this fault section, and long‐term fault interactions in southern California. Here, we use surficial geologic mapping, excavations, and radiocarbon and luminescence dating to quantify the displacements and ages of late Holocene landforms offset by the fault at three sites. At the Ranch Center site, the slip rate is determined using the base of a fan marking incision and deflection of an ephemeral channel. At the adjacent Key Slide site, the margin of a landslide deposited on indigenous fire hearths provides a minimum rate. At the X‐12 site, the slip rate is determined from a channel that incised into a broad fan surface, and is deflected and beheaded by the fault. We use maximum–minimum bounds on both the displacement and age of each offset feature to calculate slip rate for each site independently. Overlap of the three independent rate ranges yields a rate of 33–39 mm/yr over the last 3 ka, under the assumption that the sites share a common history, given their proximity. Considered in sequence, site‐level epistemic uncertainties in the data permit but do not require a rate increase since ∼1200 cal B.P. Modest rate changes can be explained by aleatory variability in earthquake timing and magnitude; larger changes could suggest a shared regional variation with the Garlock and other faults. The new late Holocene slip rates are consistent with geodetic model estimates that include a viscoelastic crust and earthquake cycle effects. The geologic slip rates also provide average slip over dozens of earthquake cycles—a key constraint for long‐term earthquake rupture forecasts.

California↗

Magnetic field variations in Alaska: Recording space weather events on seismic stations in Alaska

Seismometers are highly sensitive instruments to not only ground motion but also many other nonseismic noise sources (e.g., temperature, pressure, and magnetic field variations). We show that the Alaska component of the Transportable Array is particularly susceptible to recording magnetic storms and other space weather events because the sensors used in this network are unshielded and magnetic flux variations are stronger at higher latitudes. We also show that vertical‐component seismic records across Alaska are directly recording magnetic field variations between 40 and 800 s period as opposed to actual ground motion during geomagnetic events with sensitivities ranging from 0.004 to 0.48    ( m / s 2 ) / T "> 0.48 ( m / s 2 ) / T 0.48 (m/s2)/T ⁠ . These sensitivities were found on a day where the root mean square variation in the magnetic field was 225 nT. Using a method developed by Forbriger (2007 , his section 3.1), we show that improving vertical seismic resolution of an unshielded sensor by as much as 10 dB in the 100–400 s period band using magnetic data from a collocated three‐component magnetometer is possible. However, due to large spatial variations in Earth’s magnetic field, this methodology becomes increasingly ineffective as the distance between the seismometer and magnetometer increases (no more than 200 km separation). A potential solution to this issue may be to incorporate relatively low‐cost magnetometers as an additional environmental data stream at high‐latitude seismic stations. We demonstrate that the Bartington Mag‐690 sensors currently deployed at Global Seismographic Network sites are not only acceptable for performing corrections to seismic data, but are also capable of recording many magnetic field signals with similar signal‐to‐noise ratios, in the 20–1000 s period band, as the observatory grade magnetometers operated by the U.S. Geological Survey Geomagnetism Program. This approach would densify magnetic field observations and could also contribute to space weather monitoring by supplementing highly calibrated magnetometers with additional sensors.

Alaska↗

Finite-fault source inversion using teleseismic P waves: Simple parameterization and rapid analysis

We examine the ability of teleseismic P waves to provide a timely image of the rupture history for large earthquakes using a simple, 2D finite‐fault source parameterization. We analyze the broadband displacement waveforms recorded for the 2010 M w ∼7 Darfield (New Zealand) and El Mayor‐Cucapah (Baja California) earthquakes using a single planar fault with a fixed rake. Both of these earthquakes were observed to have complicated fault geometries following detailed source studies conducted by other investigators using various data types. Our kinematic, finite‐fault analysis of the events yields rupture models that similarly identify the principal areas of large coseismic slip along the fault. The results also indicate that the amount of stabilization required to spatially smooth the slip across the fault and minimize the seismic moment is related to the amplitudes of the observed P waveforms and can be estimated from the absolute values of the elements of the coefficient matrix. This empirical relationship persists for earthquakes of different magnitudes and is consistent with the stabilization constraint obtained from the L‐curve in Tikhonov regularization. We use the relation to estimate the smoothing parameters for the 2011 M w 7.1 East Turkey, 2012 M w 8.6 Northern Sumatra, and 2011 M w 9.0 Tohoku, Japan, earthquakes and invert the teleseismic P waves in a single step to recover timely, preliminary slip models that identify the principal source features observed in finite‐fault solutions obtained by the U.S. Geological Survey National Earthquake Information Center (USGS/NEIC) from the analysis of body‐ and surface‐wave data. These results indicate that smoothing constraints can be estimated a priori to derive a preliminary, first‐order image of the coseismic slip using teleseismic records.

Darfield, El Mayor-Cucapah↗

Length-scale-dependent relationships between VS30 and topographic slopes in southern California

In engineering seismology, the time‐averaged shear‐wave velocity ( ⁠ V S ⁠ ) of the upper 30 m of the crust ( ⁠V S30⁠ ) is the primary parameter used in ground‐motion models to predict seismic site effects. ⁠V S30⁠ is typically derived from in situ recordings of V S ⁠ , although proxy‐based approaches (using geologic and/or geomorphometric classifications) are provisionally adopted when measurement‐based ⁠V S30⁠ are sparse or not readily available. Despite the acceptance of proxy approaches, there are limited studies that examine the empirical relationships between ⁠V S30⁠ and topographic attributes measured from various length scales and different resolutions of the digital elevation model. In this study, we examine the relationships between compiled ⁠V S30⁠ measurements from 218 sites in southern California and topographic metrics of slope and relief measured over various length scales. We find that the correlations between topographic metrics and ⁠V S30⁠ are weak but statistically significant. The correlations are improved when topographic slopes and relief are measured over length scales longer than typical hillslopes and ⁠V S30⁠ sites are separated by different geologic groups. This is likely because ⁠V S30⁠ ⁠ , especially on the rock sites, is better reflected in topographic metrics that capture large‐scale topographic relief, as well as landscape positions such as hilltops and valley bottoms.

California↗

The seismic response of the Los Angeles basin, California

Using strong-motion data recorded in the Los Angeles region from the 1992 ( M w 7.3) Landers earthquake, we have tested the accuracy of existing three-dimensional (3D) velocity models on the simulation of long-period (≧2 sec) ground motions in the Los Angeles basin and surrounding San Fernando and San Gabriel Valleys. First, the overall pattern and degree of long-period excitation of the basins were identified in the observations. Within the Los Angeles basin, the recorded amplitudes are about three to four times larger than at sites outside the basins; amplitudes within the San Fernando and San Gabriel Valleys are nearly a factor of 3 greater than surrounding bedrock sites. Then, using a 3D finite-difference numerical modeling approach, we analyzed how variations in 3D earth structure affect simulated waveforms, amplitudes, and the fit to the observed patterns of amplification. Significant differences exist in the 3D velocity models of southern California that we tested (Magistrale et al. , 1996; Graves, 1996a; Hauksson and Haase, 1997). Major differences in the models include the velocity of the assumed background models; the depth of the Los Angeles basin; and the depth, location, and geometry of smaller basins. The largest disparities in the response of the models are seen for the San Fernando Valley and the deepest portion of the Los Angeles basin. These arise in large part from variations in the structure of the basins, particularly the effective depth extent, which is mainly due to alternative assumptions about the nature of the basin sediment fill. The general ground-motion characteristics are matched by the 3D model simulations, validating the use of 3D modeling with geologically based velocity-structure models. However, significant shortcomings exist in the overall patterns of amplification and the duration of the long-period response. The successes and limitations of the models for reproducing the recorded ground motions as discussed provide the basis and direction for necessary improvements to earth structure models, whether geologically or tomographically derived. The differences in the response of the earth models tested also translate to variable success in the ability to successfully model the data and add uncertainty to estimates of the basin response given input “scenario” earthquake source models.

California↗

Contemporary seismicity in and around the Yakima-Fold-and-Thrust Belt in eastern Washington

We examined characteristics of routinely cataloged seismicity from 1970 to the present in and around the Yakima fold‐and‐thrust belt (YFTB) in eastern Washington to determine if the characteristics of contemporary seismicity provide clues about regional‐scale active tectonics or about more localized, near‐surface processes. We employed new structural and hydrologic models of the Columbia River basalts (CRB) and found that one‐third to one‐half of the cataloged earthquakes occur within the CRB and that these CRB earthquakes exhibit significantly more clustered, and swarmlike, behavior than those outside. These results and inferences from published studies led us to hypothesize that clustered seismicity is likely associated with hydrologic changes in the CRB, which hosts the regional aquifer system. While some general features of the regional groundwater system support this hypothesis, seismicity patterns and mapped long‐term changes in groundwater levels and present‐day irrigation neither support nor refute it. Regional tectonic processes and crustal‐scale structures likely influence the distribution of earthquakes both outside and within the CRB as well. We based this inference on qualitatively assessed alignments between the dominant northwest trends in the geologic structure and the seismicity generally and between specific faults and characteristics of the 2009 Wooded Island swarm and aseismic slip, which is the only cluster studied in detail and the most vigorous since regional monitoring began.

Washington↗

Including foreshocks and aftershocks in time-independent probabilistic seismic hazard analyses

Time‐independent probabilistic seismic‐hazard analysis treats each source as being temporally and spatially independent; hence foreshocks and aftershocks, which are both spatially and temporally dependent on the mainshock, are removed from earthquake catalogs. Yet, intuitively, these earthquakes should be considered part of the seismic hazard, capable of producing damaging ground motions. In this study, I consider the mainshock and its dependents as a time‐independent cluster, each cluster being temporally and spatially independent from any other. The cluster has a recurrence time of the mainshock; and, by considering the earthquakes in the cluster as a union of events, dependent events have an opportunity to contribute to seismic ground motions and hazard. Based on the methods of the U.S. Geological Survey for a high‐hazard site, the inclusion of dependent events causes ground motions that are exceeded at probability levels of engineering interest to increase by about 10% but could be as high as 20% if variations in aftershock productivity can be accounted for reliably.

Bulletin of the Seismological Society of America↗

Maximum magnitude ( M max ) in the central and eastern United States for the 2014 U.S. Geological Survey Hazard Model

Probabilistic seismic‐hazard assessment (PSHA) requires an estimate of M max , the moment magnitude M of the largest earthquake that could occur within a specified area. Sparse seismicity hinders M max estimation in the central and eastern United States (CEUS) and tectonically similar regions worldwide (stable continental regions [SCRs]). A new global catalog of moderate‐to‐large SCR earthquakes is analyzed with minimal assumptions about enigmatic geologic controls on SCR M max . An earlier observation that SCR earthquakes of M 7.0 and larger occur in young (250–23 Ma) passive continental margins and associated rifts but not in cratons is not strongly supported by the new catalog. SCR earthquakes of M 7.5 and larger are slightly more numerous and reach slightly higher M in young passive margins and rifts than in cratons. However, overall histograms of M from young margins and rifts and from cratons are statistically indistinguishable. This conclusion is robust under uncertainties in M , the locations of SCR boundaries, and which of two available global SCR catalogs is used. The conclusion stems largely from recent findings that (1) large southeast Asian earthquakes once thought to be SCR were in actively deforming crust and (2) long escarpments in cratonic Australia were formed by prehistoric faulting. The 2014 seismic‐hazard model of the U.S. Geological Survey represents CEUS M max as four‐point probability distributions. The distributions have weighted averages of M 7.0 in cratons and M 7.4 in passive margins and rifts. These weighted averages are consistent with M max estimates of other SCR PSHAs of the CEUS, southeastern Canada, Australia, and India.

Bulletin of the Seismological Society of America↗

Coda duration magnitudes in central California: An empirical approach

A new empirical coda magnitude M D is presented for the central California seismic network (CALNET), in which station corrections are introduced to reduce the influence of hypocenter and station distributions. The new magnitude scale is obtained by relating the lapse-time τ^ to local magnitudes M L using data from 55 earthquakes with 1.1 ≦ M L ≤ 5.6 from June 1977 through 1981. Magnitude estimates are significantly improved by introducing corrections for the seismograph instrument attenuation setting α ij and a site correction δ j . The magnitude of the i th event 〈 M D 〉, is the median value of individual M D ij = -1.03 + 2.10log 10 τˆ ij + 0.00268 τˆ ij + α ij + δ j , where τ^ ij is the lapse-time measured at the j th station. Local site effects influence the duration of the coda and cause stations to either overestimate or underestimate 〈 M D 〉, typically up to 0.5 magnitude units. The distribution of the site corrections is spatially correlated, changing coherently across major tectonic and geologic features such as the San Andreas fault. It appears that the site corrections are influenced by the physical attenuation properties near the station, seemingly contradicting the concept that the coda is related to back-scattering and attenuation in a volume with dimensions on the order of the hypocenter to receiver distance.

California↗

MMI attenuation and historical earthquakes in the basin and range province of western North America

Earthquakes in central Nevada (1932-1959) were used to develop a modified Mercalli intensity (MMI) attenuation model for estimating moment magnitude M for earthquakes in the Basin and Range province of interior western North America. M is 7.4-7.5 for the 26 March 1872 Owens Valley, California, earthquake, in agreement with Beanland and Clark's (1994) M 7.6 that was estimated from geologic field observations. M is 7.5 for the 3 May 1887 Sonora, Mexico, earthquake, in agreement with Natali and Sbar's (1982) M 7.4 and Suter's (2006) M 7.5, both estimated from geologic field observations. MMI at sites in California for earthquakes in the Nevada Basin and Range apparently are not much affected by the Sierra Nevada except at sites near the Sierra Nevada where MMI is reduced. This reduction in MMI is consistent with a shadow zone produced by the root of the Sierra Nevada. In contrast, MMI assignments for earthquakes located in the eastern Sierra Nevada near the west margin of the Basin and Range are greater than predicted at sites in California. These higher MMI values may result from critical reflections due to layering near the base of the Sierra Nevada.

Bulletin of the Seismological Society of America↗

How to determine an earthquake rate budget? A comparison of geodetic, geologic, and seismologic moment rates

How much tectonic geodetic strain is accommodated as earthquakes? Can estimates of tectonic geodetic strain be used in probabilistic seismic hazard analyses? Understanding the balance of moment accumulation (i.e., geodetic moment) and accommodation (i.e., geologic and seismologic moment) is a necessary step toward leveraging corrected geodetic strain rates to quantify earthquake budgets. We investigate these questions using updated geodetic, geologic, and seismologic datasets compiled for the western United States (WUS) portion of the 2023 U.S. National Seismic Hazard Model (NSHM23‐WUS) to provide new views of moment rates across the WUS. We compare moment rates using a floating footprint approach across the extent of the NSHM23‐WUS fault sections database. We make these comparisons using “moment rate deviation (MRD),” which is defined as a surplus or deficit of geodetic moment rate relative to the sum of geologic and seismologic moment rates. We find that, in general, higher strain rate areas of the WUS (i.e., plate boundary) have smaller values of MRD, meaning that geodetic moment rates approximate the earthquake record preserved by geologic and seismologic moment rates well. In contrast, the lower strain rate, eastern portions of our analysis domain (i.e., intraplate) have very high values of MRD, with geodetic moment rates exceeding measured geologic and seismologic moment rates. We find that geodetic moment rates appear to be representative of the earthquake budget in relatively well‐characterized, interplate regions, and geodetic estimates of off‐fault deformation could be considered in seismic hazard analyses where MRD is minimized.

western United States↗

Normal-faulting slip maxima and stress-drop variability: A geological perspective

We present an empirical estimate of maximum slip in continental normal-faulting earthquakes and present evidence that stress drop in intraplate extensional environments is dependent on fault maturity. A survey of reported slip in historical earthquakes globally and in latest Quaternary paleoearthquakes in the Western Cordillera of the United States indicates maximum vertical displacements as large as 6–6.5 m. A difference in the ratio of maximum-to-mean displacements between data sets of prehistoric and historical earthquakes, together with constraints on bias in estimates of mean paleodisplacement, suggest that applying a correction factor of 1.4±0.3 to the largest observed displacement along a paleorupture may provide a reasonable estimate of the maximum displacement. Adjusting the largest paleodisplacements in our regional data set (∼6 m) by a factor of 1.4 yields a possible upper-bound vertical displacement for the Western Cordillera of about 8.4 m, although a smaller correction factor may be more appropriate for the longest ruptures. Because maximum slip is highly localized along strike, if such large displacements occur, they are extremely rare. Static stress drop in surface-rupturing earthquakes in the Western Cordillera, as represented by maximum reported displacement as a fraction of modeled rupture length, appears to be larger on normal faults with low cumulative geologic displacement (<2 km) and larger in regions such as the Rocky Mountains, where immature, low-throw faults are concentrated. This conclusion is consistent with a growing recognition that structural development influences stress drop and indicates that this influence is significant enough to be evident among faults within a single intraplate environment.

Bulletin of the Seismological Society of America↗