USGS Science⌕ Search

SEARCH · USGS Science

Results for “Bulletin of the Seismological 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 397 records · Page 22Linked to original sources

Comparing stochastic point-source and finite-source ground-motion simulations: SMSIM and EXSIM

Comparisons of ground motions from two widely used point-source and finite-source ground-motion simulation programs (SMSIM and EXSIM) show that the following simple modifications in EXSIM will produce agreement in the motions from a small earthquake at a large distance for the two programs: (1) base the scaling of high frequencies on the integral of the squared Fourier acceleration spectrum; (2) do not truncate the time series from each subfault; (3) use the inverse of the subfault corner frequency for the duration of motions from each subfault; and (4) use a filter function to boost spectral amplitudes at frequencies near and less than the subfault corner frequencies. In addition, for SMSIM an effective distance is defined that accounts for geometrical spreading and anelastic attenuation from various parts of a finite fault. With these modifications, the Fourier and response spectra from SMSIM and EXSIM are similar to one another, even close to a large earthquake (M 7), when the motions are averaged over a random distribution of hypocenters. The modifications to EXSIM remove most of the differences in the Fourier spectra from simulations using pulsing and static subfaults; they also essentially eliminate any dependence of the EXSIM simulations on the number of subfaults. Simulations with the revised programs suggest that the results of Atkinson and Boore (2006), computed using an average stress parameter of 140 bars and the original version of EXSIM, are consistent with the revised EXSIM with a stress parameter near 250 bars.

Bulletin of the Seismological Society of America↗

Effects of 3D random correlated velocity perturbations on predicted ground motions

Three-dimensional, finite-difference simulations of a realistic finite-fault rupture on the southern Hayward fault are used to evaluate the effects of random, correlated velocity perturbations on predicted ground motions. Velocity perturbations are added to a three-dimensional (3D) regional seismic velocity model of the San Francisco Bay Area using a 3D von Karman random medium. Velocity correlation lengths of 5 and 10 km and standard deviations in the velocity of 5% and 10% are considered. The results show that significant deviations in predicted ground velocities are seen in the calculated frequency range (≤1 Hz) for standard deviations in velocity of 5% to 10%. These results have implications for the practical limits on the accuracy of scenario ground-motion calculations and on retrieval of source parameters using higher-frequency, strong-motion data.

Bulletin of the Seismological Society of America↗

The influence of maximum magnitude on seismic-hazard estimates in the Central and Eastern United States

I analyze the sensitivity of seismic-hazard estimates in the central and eastern United States (CEUS) to maximum magnitude (m max ) by exercising the U.S. Geological Survey (USGS) probabilistic hazard model with several m max alternatives. Seismicity-based sources control the hazard in most of the CEUS, but data seldom provide an objective basis for estimating m max . The USGS uses preferred m max values of moment magnitude 7.0 and 7.5 for the CEUS craton and extended margin, respectively, derived from data in stable continental regions worldwide. Other approaches, for example analysis of local seismicity or judgment about a source's seismogenic potential, often lead to much smaller m max . Alternative models span the m max ranges from the 1980s Electric Power Research Institute/Seismicity Owners Group (EPRI/SOG) analysis. Results are presented as haz-ard ratios relative to the USGS national seismic hazard maps. One alternative model specifies m max equal to moment magnitude 5.0 and 5.5 for the craton and margin, respectively, similar to EPRI/SOG for some sources. For 2% probability of exceedance in 50 years (about 0.0004 annual probability), the strong m max truncation produces hazard ratios equal to 0.35-0.60 for 0.2-sec spectral acceleration, and 0.15-0.35 for 1.0-sec spectral acceleration. Hazard-controlling earthquakes interact with m max in complex ways. There is a relatively weak dependence on probability level: hazardratios increase 0-15% for 0.002 annual exceedance probability and decrease 5-25% for 0.00001 annual exceedance probability. Although differences at some sites are tempered when faults are added, m max clearly accounts for some of the discrepancies that are seen in comparisons between USGS-based and EPRI/SOG-based hazard results.

Bulletin of the Seismological Society of America↗

On the composition of earth's short-period seismic noise field

In the classic microseismic band of 5-20 sec, seismic noise consists mainly of fundamental mode Rayleigh and Love waves; however, at shorter periods seismic noise also contains a significant amount of body-wave energy and higher mode surface waves. In this study we perform a global survey of Earth's short-period seismic noise field with the goal of quantifying the relative contributions of these propagation modes. We examined a year's worth of vertical component data from 18 seismic arrays of the International Monitoring System that were sited in a variety of geologic environments. The apertures of the arrays varied from 2 to 28 km, constraining the periods we analyzed to 0.25-2.5 sec. Using frequency-wavenumber analysis we identified the apparent velocity for each sample of noise and classified its mode of propagation. The dominant component was found to be L g , occurring in about 50% of the noise windows. Because L g does not propagate across ocean-continent boundaries, this energy is most likely created in shallow water areas near coastlines. The next most common component was P-wave energy, which accounted for about 28% of the noise windows. These were split between regional P waves (P n =P g at 6%), mantle bottoming P waves (14%), and core-sensitive waves (PKP at 8%). This energy is mostly generated in deep water away from coastlines, with a region of the North Pacific centered at 165?? W and 40?? N being especially prolific. The remainder of the energy arriving in the noise consisted of R g waves (28%), a large fraction of which may have a cultural origin. Hence, in contrast to the classic micro-seismic band of 5-20 sec, at shorter periods fundamental mode Rayleigh waves are the least significant component.

Bulletin of the Seismological Society of America↗

Evidence for a twelfth large earthquake on the southern hayward fault in the past 1900 years

We present age and stratigraphic evidence for an additional paleoearthquake at the Tyson Lagoon site. The acquisition of 19 additional radiocarbon dates and the inclusion of this additional event has resolved a large age discrepancy in our earlier earthquake chronology. The age of event E10 was previously poorly constrained, thus increasing the uncertainty in the mean recurrence interval (RI), a critical factor in seismic hazard evaluation. Reinspection of many trench logs revealed substantial evidence suggesting that an additional earthquake occurred between E10 and E9 within unit u45. Strata in older u45 are faulted in the main fault zone and overlain by scarp colluviums in two locations.We conclude that an additional surfacerupturing event (E9.5) occurred between E9 and E10. Since 91 A.D. (??40 yr, 1??), 11 paleoearthquakes preceded the M 6:8 earthquake in 1868, yielding a mean RI of 161 ?? 65 yr (1??, standard deviation of recurrence intervals). However, the standard error of the mean (SEM) is well determined at ??10 yr. Since ~1300 A.D., the mean rate has increased slightly, but is indistinguishable from the overall rate within the uncertainties. Recurrence for the 12-event sequence seems fairly regular: the coefficient of variation is 0.40, and it yields a 30-yr earthquake probability of 29%. The apparent regularity in timing implied by this earthquake chronology lends support for the use of time-dependent renewal models rather than assuming a random process to forecast earthquakes, at least for the southern Hayward fault.

Bulletin of the Seismological Society of America↗

Rupture directivity of moderate earthquakes in northern California

We invert peak ground velocity and acceleration (PGV and PGA) to estimate rupture direction and rupture velocity for 47 moderate earthquakes (3.5≥M≥5.4) in northern California. We correct sets of PGAs and PGVs recorded at stations less than 55–125 km, depending on source depth, for site amplification and source–receiver distance, then fit the residual peak motions to the unilateral directivity function of Ben-Menahem (1961). We independently invert PGA and PGV. The rupture direction can be determined using as few as seven peak motions if the station distribution is sufficient. The rupture velocity is unstable, however, if there are no takeoff angles within 30° of the rupture direction. Rupture velocities are generally subsonic (0.5β–0.9β); for stability, we limit the rupture velocity at v=0.92β, the Rayleigh wave speed. For 73 of 94 inversions, the rupture direction clearly identifies one of the nodal planes as the fault plane. The 35 strike-slip earthquakes have rupture directions that range from nearly horizontal (6 events) to directly updip (5 events); the other 24 rupture partly along strike and partly updip. Two strike-slip earthquakes rupture updip in one inversion and downdip in the other. All but 1 of the 11 thrust earthquakes rupture predominantly updip. We compare the rupture directions for 10 M≥4.0 earthquakes to the relative location of the mainshock and the first two weeks of aftershocks. Spatial distributions of 8 of 10 aftershock sequences agree well with the rupture directivity calculated for the mainshock.

California↗

Localized surface disruptions observed by InSAR during strong earthquakes in Java and Hawai'i

Interferometric Synthetic Aperture Radar data spanning strong earthquakes on the islands of Java and Hawai‘i in 2006 reveal patches of subsidence and incoherence indicative of localized ground failure. Interferograms spanning the 26 May 2006 Java earthquake suggest an area of about 7.5 km 2 of subsidence (~2 cm) and incoherence south of the city of Yogyakarta that correlates with significant damage to housing, high modeled peak ground accelerations, and poorly consolidated geologic deposits. The subsidence and incoherence is inferred to be a result of intense shaking and/or damage. At least five subsidence patches on the west side of the Island of Hawai‘i, ranging 0.3–2.2 km 2 in area and 3–8 cm in magnitude, occurred as a result of a pair of strong earthquakes on 15 October 2006. Although no felt reports or seismic data are available from the areas in Hawai‘i, the Java example suggests that the subsidence patches indicate areas of amplified earthquake shaking. Surprisingly, all subsidence areas in Hawai‘i were limited to recent, and supposedly stable, lava flows and may reflect geological conditions not detectable at the surface. In addition, two ‘a‘ā lava flows in Hawai‘i were partially incoherent in interferograms spanning the earthquakes, indicating surface disruption as a result of the earthquake shaking. Coearthquake incoherence of rubbly deposits, like ‘a‘ā flows, should be explored as a potential indicator of earthquake intensity and past strong seismic activity.

Bulletin of the Seismological Society of America↗

Fundamental questions of earthquake statistics, source behavior, and the estimation of earthquake probabilities from possible foreshocks

Estimates of the probability that an M L 4.8 earthquake, which occurred near the southern end of the San Andreas fault on 24 March 2009, would be followed by an M 7 mainshock over the following three days vary from 0.0009 using a Gutenberg–Richter model of aftershock statistics (Reasenberg and Jones, 1989) to 0.04 using a statistical model of foreshock behavior and long‐term estimates of large earthquake probabilities, including characteristic earthquakes (Agnew and Jones, 1991). I demonstrate that the disparity between the existing approaches depends on whether or not they conform to Gutenberg–Richter behavior. While Gutenberg–Richter behavior is well established over large regions, it could be violated on individual faults if they have characteristic earthquakes or over small areas if the spatial distribution of large‐event nucleations is disproportional to the rate of smaller events. I develop a new form of the aftershock model that includes characteristic behavior and combines the features of both models. This new model and the older foreshock model yield the same results when given the same inputs, but the new model has the advantage of producing probabilities for events of all magnitudes, rather than just for events larger than the initial one. Compared with the aftershock model, the new model has the advantage of taking into account long‐term earthquake probability models. Using consistent parameters, the probability of an M 7 mainshock on the southernmost San Andreas fault is 0.0001 for three days from long‐term models and the clustering probabilities following the M L 4.8 event are 0.00035 for a Gutenberg–Richter distribution and 0.013 for a characteristic‐earthquake magnitude–frequency distribution. Our decisions about the existence of characteristic earthquakes and how large earthquakes nucleate have a first‐order effect on the probabilities obtained from short‐term clustering models for these large events.

Bulletin of the Seismological Society of America↗

Tsunami simulations of the 1867 Virgin Island earthquake: Constraints on epicenter location and fault parameters

The 18 November 1867 Virgin Island earthquake and the tsunami that closely followed caused considerable loss of life and damage in several places in the northeast Caribbean region. The earthquake was likely a manifestation of the complex tectonic deformation of the Anegada Passage, which cuts across the Antilles island arc between the Virgin Islands and the Lesser Antilles. In this article, we attempt to characterize the 1867 earthquake with respect to fault orientation, rake, dip, fault dimensions, and first tsunami wave propagating phase, using tsunami simulations that employ high-resolution multibeam bathymetry. In addition, we present new geophysical and geological observations from the region of the suggested earthquake source. Results of our tsunami simulations based on relative amplitude comparison limit the earthquake source to be along the northern wall of the Virgin Islands basin, as suggested by Reid and Taber (1920) , or on the carbonate platform north of the basin, and not in the Virgin Islands basin, as commonly assumed. The numerical simulations suggest the 1867 fault was striking 120°–135° and had a mixed normal and left-lateral motion. First propagating wave phase analysis suggests a fault striking 300°–315° is also possible. The best-fitting rupture length was found to be relatively small (50 km), probably indicating the earthquake had a moment magnitude of ∼7.2. Detailed multibeam echo sounder surveys of the Anegada Passage bathymetry between St. Croix and St. Thomas reveal a scarp, which cuts the northern wall of the Virgin Islands basin. High-resolution seismic profiles further indicate it to be a reasonable fault candidate. However, the fault orientation and the orientation of other subparallel faults in the area are more compatible with right-lateral motion. For the other possible source region, no clear disruption in the bathymetry or seismic profiles was found on the carbonate platform north of the basin.

Bulletin of the Seismological Society of America↗

LiDAR-Assisted identification of an active fault near Truckee, California

We use high-resolution (1.5-2.4 points/m2) bare-earth airborne Light Detection and Ranging (LiDAR) imagery to identify, map, constrain, and visualize fault-related geomorphology in densely vegetated terrain surrounding Martis Creek Dam near Truckee, California. Bare-earth LiDAR imagery reveals a previously unrecognized and apparently youthful right-lateral strike-slip fault that exhibits laterally continuous tectonic geomorphic features over a 35-km-long zone. If these interpretations are correct, the fault, herein named the Polaris fault, may represent a significant seismic hazard to the greater Truckee-Lake Tahoe and Reno-Carson City regions. Three-dimensional modeling of an offset late Quaternary terrace riser indicates a minimum tectonic slip rate of 0.4 ?? 0.1 mm/yr.Mapped fault patterns are fairly typical of regional patterns elsewhere in the northern Walker Lane and are in strong coherence with moderate magnitude historical seismicity of the immediate area, as well as the current regional stress regime. Based on a range of surface-rupture lengths and depths to the base of the seismogenic zone, we estimate a maximum earthquake magnitude (M) for the Polaris fault to be between 6.4 and 6.9.

Bulletin of the Seismological Society of America↗

Short baseline variations in site response and wave-propagation effects and their structural causes: Four examples in and around the santa clara valley, California

Ground motion records of local and regional events from a portable array are used to investigate the structural causes of variations in ground motion over distances of a few hundred meters to a few kilometers in the sedimentary basin environment of the Santa Clara Valley, California, and its margins. Arrays of portable seismic stations are used to target four study areas with different ground motion patterns: (1) an edge of the alluvial basin extending up onto a marginal ridge (Blossom Hill), (2) a Cenozoic basin with a nearly flat bottom (Cupertino Basin), (3) a long, narrow Cenozoic basin with a steep V profile (Evergreen Basin), and (4) a line perpendicular to the trace of the Hayward fault. Average peak velocities on Blossom Hill from local earthquakes are a factor of 2.5 times higher than nearby valley sites. Three-dimensional (3D) modeling is used to conclude that the majority of the amplification is due to lower shear-wave velocities along a local fault zone (Shannon–Berrocal). Site amplification over the Cupertino Basin in the frequency band 0.5–4 Hz is generally low (less than 2.0 relative to a Mesozoic rock site) and spatially uniform. This response is attributed to the shallow, flat-bottomed shape of the basin and the uniform, flat-laying sedimentary fill. In contrast, site amplification in the Evergreen Basin generally exceeds 3.0 and is attributed to the deep, V-shaped geometry of the basin and younger sedimentary fill. 3D waveform modeling shows the elongated shape of the Evergreen Basin causes more efficient trapping of long-period waves for sources along the long axis of the basin. A low-velocity zone is postulated along the Hayward fault with a width between 100 and 200 m, based on elevated site response along the fault trace and 4.5-Hz fault zone guided waves on the horizontal components of stations near the fault.

California↗

Seismotectonics and fault structure of the California Central Coast

I present and interpret new earthquake relocations and focal mechanisms for the California Central Coast. The relocations improve upon catalog locations by using 3D seismic velocity models to account for lateral variations in structure and by using relative arrival times from waveform cross-correlation and double-difference methods to image seismicity features more sharply. Focal mechanisms are computed using ray tracing in the 3D velocity models. Seismicity alignments on the Hosgri fault confirm that it is vertical down to at least 12 km depth, and the focal mechanisms are consistent with right-lateral strike-slip motion on a vertical fault. A prominent, newly observed feature is an ~25 km long linear trend of seismicity running just offshore and parallel to the coastline in the region of Point Buchon, informally named the Shoreline fault. This seismicity trend is accompanied by a linear magnetic anomaly, and both the seismicity and the magnetic anomaly end where they obliquely meet the Hosgri fault. Focal mechanisms indicate that the Shoreline fault is a vertical strike-slip fault. Several seismicity lineations with vertical strike-slip mechanisms are observed in Estero Bay. Events greater than about 10 km depth in Estero Bay, however, exhibit reverse-faulting mechanisms, perhaps reflecting slip at the top of the remnant subducted slab. Strike-slip mechanisms are observed offshore along the Hosgri–San Simeon fault system and onshore along the West Huasna and Rinconada faults, while reverse mechanisms are generally confined to the region between these two systems. This suggests a model in which the reverse faulting is primarily due to restraining left-transfer of right-lateral slip.

California↗

Repeatable source, site, and path effects on the standard deviation for empirical ground-motion prediction models

In this study, we quantify the reduction in the standard deviation for empirical ground-motion prediction models by removing ergodic assumption.We partition the modeling error (residual) into five components, three of which represent the repeatable source-location-specific, site-specific, and path-specific deviations from the population mean. A variance estimation procedure of these error components is developed for use with a set of recordings from earthquakes not heavily clustered in space.With most source locations and propagation paths sampled only once, we opt to exploit the spatial correlation of residuals to estimate the variances associated with the path-specific and the source-location-specific deviations. The estimation procedure is applied to ground-motion amplitudes from 64 shallow earthquakes in Taiwan recorded at 285 sites with at least 10 recordings per site. The estimated variance components are used to quantify the reduction in aleatory variability that can be used in hazard analysis for a single site and for a single path. For peak ground acceleration and spectral accelerations at periods of 0.1, 0.3, 0.5, 1.0, and 3.0 s, we find that the singlesite standard deviations are 9%-14% smaller than the total standard deviation, whereas the single-path standard deviations are 39%-47% smaller.

Bulletin of the Seismological Society of America↗

Ground-motion modeling of Hayward fault scenario earthquakes, part I: Construction of the suite of scenarios

We construct kinematic earthquake rupture models for a suite of 39 M w 6.6-7.2 scenario earthquakes involving the Hayward, Calaveras, and Rodgers Creek faults. We use these rupture models in 3D ground-motion simulations as discussed in Part II (Aagaard et al. , 2010) to provide detailed estimates of the shaking for each scenario. We employ both geophysical constraints and empirical relations to provide realistic variation in the rupture dimensions, slip heterogeneity, hypocenters, rupture speeds, and rise times. The five rupture lengths include portions of the Hayward fault as well as combined rupture of the Hayward and Rodgers Creek faults and the Hayward and Calaveras faults. We vary rupture directivity using multiple hypocenters, typically three per rupture length, yielding north-to-south rupture, bilateral rupture, and south-to-north rupture. For each rupture length and hypocenter, we consider multiple random distributions of slip. We use two approaches to account for how aseismic creep might reduce coseismic slip. For one subset of scenarios, we follow the slip-predictable approach and reduce the nominal slip in creeping regions according to the creep rate and time since the most recent earthquake, whereas for another subset of scenarios we apply a vertical gradient to the nominal slip in creeping regions. The rupture models include local variations in rupture speed and use a ray-tracing algorithm to propagate the rupture front. Although we are not attempting to simulate the 1868 Hayward fault earthquake in detail, a few of the scenarios are designed to have source parameters that might be similar to this historical event.

Bulletin of the Seismological Society of America↗

Reply to “Comment on ‘Is There a Basis for Preferring Characteristic Earthquakes over a Gutenberg–Richter Distribution in Probabilistic Earthquake Forecasting?’ by Tom Parsons and Eric L. Geist” by Jens-Uwe Klügel

The focus of Parsons and Geist (2009) was to test whether the key observational data used in earthquake forecasting necessitate a characteristic earthquake rupture model. The point of our article was not to suggest that a specific form of the Gutenberg–Richter earthquake distribution is a perfect representation of reality. The uncertainties surrounding event slip estimates, paleoseismic event rates, and observed a and b values in catalog magnitude–frequency distributions are broad. So broad, in fact, that giving full weight to just one model of earthquake rupture behavior in formal forecasting is unjustified. Further, the characteristic earthquake model requires definition of rupture segments, which introduces a series of unquantifiable uncertainties that are seldom addressed in forecasts (e.g., Field et al., 2009).

Bulletin of the Seismological Society of America↗

Reply to “Comment on ‘Near-surface location, geometry, and velocities of the Santa Monica fault zone, Los Angeles, California’ by R. D. Catchings, G. Gandhok, M. R. Goldman, D. Okaya, M. J. Rymer, and G. W. Bawden” by T. L. Pratt and J. F. Dolan

In a comment on our 2008 paper (Catchings, Gandhok, et al. , 2008) on the Santa Monica fault in Los Angeles, California, Pratt and Dolan (2010) (herein referred to as P&D) cite numerous objections to our work, inferring that our study is flawed. However, as shown in our reply, their objections contradict their own published works, published works of others, and proven seismic methodologies. Rather than responding to each repeated invalid objection, we address their objections by topic in the subsequent sections. In Catchings, Gandhok, et al. (2008), we presented high-resolution seismic-reflection images that showed two near-surface faults in the upper 50 m beneath the grounds of the Wadsworth Veterans Administration Hospital (WVAH). Although P&D suggest we effectively duplicated their seismic acquisition, our survey was not a duplication of their efforts. Rather, we conducted a seismic-imaging survey over a similar profile as Pratt et al. (1998) but used a different data acquisition system and different data processing methods to evaluate methods of seismically imaging blind faults in the wake of the 17 January 1994 M 6.7 Northridge earthquake. We used an acquisition method that provides both tomographic seismic velocities and reflection images. Our combined-data approach allowed for shallower imaging (∼2.5 m minimum) than the ∼20-m minimum of Pratt et al. (1998), clearer images of the fault zone, and more accurate depth determinations (rather than time images). In processing the reflection images, we used prestack depth migration, which is generally accepted as the only proper imaging method for imaging subsurface structures with strong lateral velocity variations (Versteeg, 1993), a condition shown to exist at the WVAH site. We correlated our reflection images with refraction tomography images, borehole lithology, and velocity data, Interferometric Synthetic Aperture Radar images, and changes in groundwater depths. Except for some minor differences, our seismic-reflection images coincide with previously published seismic-reflection images by Dolan and Pratt (1997) and Pratt et al. (1998), and a paleoseismic study by Dolan et al. (2000). Principal differences among our interpretations and those of Pratt et al. (1998) relate to the upper 20 m and the south side of the fault, which Pratt et al. (1998) did not clearly image. In contrast, our seismic images included structures on both sides of the fault zone from about 2.5 m depth to about 100 m depth at WVAH, allowing us to interpret more details.

California↗

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↗

Surface-wave potential for triggering tectonic (nonvolcanic) tremor

Source processes commonly posed to explain instances of remote dynamic triggering of tectonic (nonvolcanic) tremor by surface waves include frictional failure and various modes of fluid activation. The relative potential for Love- and Rayleigh-wave dynamic stresses to trigger tectonic tremor through failure on critically stressed thrust and vertical strike-slip faults under the Coulomb-Griffith failure criteria as a function of incidence angle is anticorrelated over the 15- to 30-km-depth range that hosts tectonic tremor. Love-wave potential is high for strike-parallel incidence on low-angle reverse faults and null for strike-normal incidence; the opposite holds for Rayleigh waves. Love-wave potential is high for both strike-parallel and strike-normal incidence on vertical, strike-slip faults and minimal for ~45?? incidence angles. The opposite holds for Rayleigh waves. This pattern is consistent with documented instances of tremor triggered by Love waves incident on the Cascadia mega-thrust and the San Andreas fault (SAF) in central California resulting from shear failure on weak faults (apparent friction, ????? 0.2). However, documented instances of tremor triggered by surface waves with strike-parallel incidence along the Nankai megathrust beneath Shikoku, Japan, is associated primarily with Rayleigh waves. This is consistent with the tremor bursts resulting from mixed-mode failure (crack opening and shear failure) facilitated by near-lithostatic ambient pore pressure, low differential stress, with a moderate friction coefficient (?? ~ 0.6) on the Nankai subduction interface. Rayleigh-wave dilatational stress is relatively weak at tectonic tremor source depths and seems unlikely to contribute significantly to the triggering process, except perhaps for an indirect role on the SAF in sustaining tremor into the Rayleigh-wave coda that was initially triggered by Love waves.

Bulletin of the Seismological Society of America↗