Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress

Size: px
Start display at page:

Download "Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress"

Transcription

1 JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi: /2010jb007586, 2010 Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress Jeanne Hardebeck 1 Received 23 March 2010; revised 27 July 2010; accepted 24 August 2010; published 3 December [1] It has been proposed that the crustal stress field contains small length scale heterogeneity of much larger amplitude than the uniform background stress. This model predicts that earthquake focal mechanisms should reflect the loading stress rather than the uniform background stress. So, if the heterogeneous stress hypothesis is correct, focal mechanisms before and after a large earthquake should align with the tectonic loading and the earthquake induced static stress perturbation, respectively. However, I show that the off fault triggered aftershocks of the 1992 M7.3 Landers, California, earthquake align with the same stress field as the pre Landers mechanisms. The aftershocks occurred on faults that were well oriented for failure in the pre Landers stress field and then loaded by the Landers induced static stress change. Aftershocks in regions experiencing a 0.05 to 5 MPa coseismic differential stress change align with the modeled Landers induced static stress change, implying that they were triggered by the stress perturbation. Contrary to the heterogeneous stress hypothesis, these triggered aftershocks are also well aligned with the pre Landers stress field obtained from inverting the pre Landers focal mechanisms. Therefore, the inverted pre Landers stress must represent the persistent background stress field. Earthquake focal mechanisms provide an unbiased sample of the spatially coherent background stress field, which is large relative to any small scale stress heterogeneity. The counterexample provided by the Landers earthquake is strong evidence that the heterogeneous stress model is not widely applicable. Citation: Hardebeck, J. (2010), Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress, J. Geophys. Res., 115,, doi: /2010jb Introduction [2] Some have argued for small scale spatial heterogeneity in the crustal stress field with large amplitude compared to the uniform background stress [e.g., Smith, 2006; Smith and Dieterich, 2010; see also edu/caltechetd:etd ]. There are reasonable arguments both for and against highly heterogeneous stress in the crust. Earthquake slip models exhibit large changes in displacement over 1 km distances, implying large localized stress changes that could produce heterogeneity [e.g., Wald and Heaton, 1994; Aagaard and Heaton, 2008]. However, finite slip models are nonunique because slip at depth is usually underdetermined [e.g., Liu et al., 2006; Page et al., 2009], so it is unclear how well resolved apparently sharp changes in displacement actually are. Large changes in surface displacement over 1 km along strike distances [McGill and Rubin, 1999] have also been cited as evidence for localized high stress [e.g., Smith, 2006; Smith and Dieterich, 2010]. However, large changes in displacement only imply large stresses if the deformation is elastic. Given the relatively low strength of near surface 1 U.S. Geological Survey, Menlo Park, California, USA. This paper is not subject to U.S. copyright. Published in 2010 by the American Geophysical Union. materials, it is much more likely that the deformation reflects permanent deformation of near surface materials, as the original authors concluded [McGill and Rubin, 1999]. [3] Borehole breakouts indicate that stress orientations in the upper few kilometers can sometimes change significantly over length scales as small as 10 m [e.g., Wilde and Stock, 1997; Hickman and Zoback, 2004]. The breakout heterogeneity, if present, does not necessarily extrapolate to seismogenic depths, however, where greater confining stresses may dampen the relative amplitude of fluctuations. At depth, heterogeneous focal mechanisms have been observed for aftershocks along parts of the rupture surfaces of large earthquakes, including the 1989 M6.9 Loma Prieta [Michael et al., 1990] and 1992 M7.3 Landers [Hauksson, 1994] main shocks in California, implying local stress heterogeneity caused by the main shock. A similar argument, based on the apparent variability of the full southern California focal mechanism catalog, has also been made for heterogeneous stress in the entire upper crust [Rivera and Kanamori, 2002]. However, when differences in tectonic regime and the sizable focal mechanism uncertainties are accounted for, the southern California mechanisms appear much more uniform [Hardebeck, 2006]. [4] One of the predictions of the heterogeneous stress model is that earthquake focal mechanisms during the interseismic period provide a biased sample of the stress field 1of10

2 [e.g., Smith and Dieterich, 2010]. In this model, earthquakes should preferentially occur in volumes of crust where the stress field is aligned with the tectonic loading, as the tectonic loading would tend to push the well oriented faults in those volumes closer to failure, much more than faults in other volumes. Focal mechanisms should reflect the orientation of the tectonic loading, and stress inversion of focal mechanisms [e.g., Michael, 1984] during interseismic periods would return the loading tensor. Similarly, after a large earthquake, the focal mechanisms of triggered aftershocks would reflect the stress perturbations from the earthquake, because faults in volumes where stress was aligned with the earthquake induced stress changes would be preferentially pushed toward failure. Smith and Dieterich [2010] hypothesized that apparent stress rotations due to large earthquakes [e.g., Michael, 1987; Hauksson, 1994] may be artifacts from sampling different spatial volumes in a heterogeneous stress field. [5] The prediction that focal mechanisms should reflect the loading stress, rather than the background stress, provides an opportunity to test the heterogeneous stress hypothesis. Large earthquakes are known to change the stress field [e.g., King et al., 1994], producing a temporary loading stress that differs from the usual tectonic loading. Therefore, if the heterogeneous stress hypothesis is correct, the pre earthquake events and the triggered aftershocks should reflect different stress fields. This can be tested by inverting the pre earthquake focal mechanisms for the apparent prestress and checking whether the aftershocks are aligned with this stress field. [6] The heterogeneous stress model predicts that the apparent prestress found from inversion of pre earthquake focal mechanisms is actually the tectonic loading stress. The aftershock focal mechanisms should not be preferentially aligned with the prestress, as they should be independent of the tectonic loading, assuming that the cumulative tectonic loading over the time period of the aftershock catalog is negligible. Alternatively, if the homogeneous part of the stress field is larger than the local variations, inversion of the pre earthquake focal mechanisms for prestress would result in the uniform background stress. In this case, the aftershock focal mechanisms should be consistent with a combination of the background and earthquake induced stress fields. The aftershocks should therefore be consistent with the prestress, assuming they are far enough away from the main shock that the stress change is sufficiently small compared to the background stress that a true rotation of the stress field could not occur. Therefore, the key test is whether or not the focal mechanisms of triggered aftershocks align with the inverted pre earthquake stress field. 2. Data [7] The 1992 M7.3 Landers, California, earthquake triggered aftershocks over a large region, and most aftershocks are consistent with the Landers induced static stress change field [King et al., 1994; Hardebeck et al., 1998]. The Landers aftershocks were well recorded by the Southern California Seismic Network [Hauksson et al., 1993] and provide an ideal data set to test whether the triggered aftershocks were also aligned with the apparent pre Landers stress field. Southern California exhibits complex fault structure, associated with the large scale bend in the San Andreas Fault, and so it is a good candidate location for a heterogeneous stress field. The complex fault geometry also ensures that there are a variety of available fault orientations, so changes in preferred failure orientation should be detectable. [8] I use a focal mechanism catalog for southern California generated by the method of Hardebeck and Shearer [2002], selecting only quality A C mechanisms without multiples. The pre Landers catalog is complicated by the occurrence of the M6.2 Joshua Tree earthquake 2 months prior to the Landers main shock, and the post Landers catalog is complicated by the occurrence of the M6.5 Big Bear aftershock a few hours after Landers. To simplify this situation, the Joshua Tree, Landers, and Big Bear earthquakes are considered together as the main shock. The pre Landers catalog (Figure 1a) uses event from 1984 until the time of the Joshua Tree earthquake, and the post Landers catalog (Figure 1b) starts just after the Big Bear earthquake. The focal mechanisms are shown by the P axes and T axes in the insets to Figures 1a and 1b. [9] Two stress fields are required for the test: the inverted pre Landers stress field and the Landers induced static stress change field. The pre Landers focal mechanisms were inverted for a spatially varying stress field using the 2 D damped stress inversion method of Hardebeck and Michael [2006] (Figure 1c). This inversion assumes that the stress field does not change significantly over the 8 years of the pre Landers catalog and that the stress orientations do not change with depth. These assumptions are justified by Hardebeck and Hauksson s [2001] observation that stress orientations in southern California do not change substantially with depth and do not change significantly with time except in areas with recent large earthquakes. To test the sensitivity to the grid configuration, I tried four binning schemes, two with 8 km grid spacing and two with 10 km spacing, with each pair shifted relative to each other by half the grid spacing. I allow 8 10 km scale variations in the pre Landers stress field because my tests target stress heterogeneity on much smaller scales of 10 m to 1 km. [10] The earthquake induced static stress changes (Figure 1d) were computed using DLC software codes (R. Simpson, personal communication) based on the Okada [1992] subroutines for dislocations in an elastic half space. Because the Joshua Tree, Landers, and Big Bear earthquakes are being considered together as the main shock, the static stress changes induced by all three were modeled. Wald and Heaton s [1994] source model was used for the Landers earthquake, and uniform slip models were constructed for Joshua Tree and Big Bear earthquakes (Table 1) on the basis of their seismic moments and apparent fault orientations defined by seismicity [Hauksson et al., 1993]. [11] A subset of the focal mechanism catalog is selected such that the Landers induced stress change at each event hypocenter is large relative to the tectonic loading over the time scale of the aftershock data set, but small enough that it could not cause a true rotation of the uniform background stress. The magnitude of the stress change is measured as the differential stress (s1 s3) of the stress change tensor, where s1 and s3 are the maximum and minimum compressive stresses, respectively. The upper cutoff of the differential stress change is chosen to be 5 MPa, which is small compared to most estimates of the magnitude of the crustal 2of10

3 deviatoric stress, and so should not cause a true stress rotation. Additionally, all aftershocks within 5 km of the modeled fault planes are removed, because the stress changes this close to the fault are poorly constrained due to the coarse slip distribution model, and it is undesirable for the data set to be dominated by events with high but poorly modeled stress changes. The lower cutoff differential stress is chosen to be 5 kpa, based on the observation of Jackson et al. [1997] that most of the region under consideration exhibits <0.15 microstrain/yr, or <4.5 kpa/yr, of tectonic Figure 1 3of10

4 Table 1. Parameters for the Uniform Slip Source Models for the Joshua Tree and Big Bear Main Shocks End Point 1 End Point 2 Longitude (deg) Latitude (deg) Longitude (deg) Latitude (deg) Dip (deg) Width (km) Rake (deg) Slip (m) Joshua Tree Big Bear loading. The aftershock catalog ends 1 year after the Big Bear earthquake, so the cumulative tectonic loading is less than the earthquake induced loading over most of the region. Hardebeck et al. [1998] found significant aftershock triggering at Landers for stress changes of 10 kpa to 1 MPa, so this data set should encompass most of the off fault triggered aftershocks. 3. Method [12] I found the misfit of each focal mechanism in the catalog with respect to two stress tensors: the inverted stress field prior to the main shock (Figure 1c) and the static stress change due to the main shock (Figure 1d). The pre Landers stress is taken at the nearest stress inversion grid point, whereas the Landers induced stress change is computed at each hypocenter. Misfit is calculated as the angle between the rake vector and the shear stress vector projected on the fault plane, as is often done in stress inversions [e.g., Michael, 1984]. [13] The misfits of pre Landers events to the prestress are expected to be low, because the stress field was inverted from these events. These events are independent of the Landers stress change, so their misfits to the static stress change field should somewhat uniformly fill the range of possible misfit values, If the heterogenous stress hypothesis is correct, aftershocks are controlled only by the main shock stress change and are independent of the inverted prestress. So aftershocks should have low misfit relative to the modeled static stress change and somewhat uniformly fill the range of possible misfit values relative to the inverted prestress. However, if stress is more homogeneous, the aftershocks are controlled by both the stress change and the prestress and should exhibit low misfit relative to both. Any background events misidentified as aftershocks will have misfits in the same range as the pre Landers events. Therefore, simple scatterplots or histograms of the misfit of each aftershock relative to the two stress fields provide a qualitative evaluation of the heterogeneous stress hypothesis. [14] I also performed a quantitative test of whether or not the pre Landers and post Landers focal mechanisms differ significantly in their misfits to the two stress tensors. In this quantitative test, I considered the magnitude of the Landersinduced stress change and divided the catalog into nine subcatalogs based on the differential (s1 s3) stress change at the event hypocenter. For each subcatalog, I subtracted the average misfit of events before Landers from the average misfit of events after Landers to find the change in misfit. If the post Landers events are dominated by aftershocks, the misfit change should be negative with respect to the Landers induced static stress change. For heterogeneous stress, the misfit change relative to the prestress should be positive for catalogs dominated by aftershocks, which should be less dependent on the prestress than the pre Landers events. For homogeneous stress, the misfits of the aftershocks and pre Landers events to the prestress should be approximately the same. [15] The significance of the misfit change is assessed by comparison to the null hypothesis that the pre Landers and post Landers events sample the same stress field. Another way to state the null hypothesis is that it is unimportant which events come before or after Landers. This null hypothesis is implemented by reshuffling, in each of the nine subcatalogs, which events occur before and after Landers, keeping the original number of each. For the misfit with respect to the Landers induced static stress changes, 3000 reshufflings are performed. The misfit with respect to the background stress requires a new inversion for background stress for each reshuffled catalog, and I performed 1000 reshufflings for each of the four inversion grid configurations. 4. Results [16] The post Landers focal mechanisms are concentrated toward low misfits with respect to the Landers induced Figure 1. (a) Pre Landers earthquakes, from 1 January 1984 until the time of the 23 April 1992 M6.2 Joshua Tree earthquake, 2 months prior to the 28 June 1992 M7.3 Landers main shock. Color reflects the magnitude of the modeled Landers static stress change at the event hypocenter, expressed as the differential stress (s1 s3). Only the 2223 events with differential stress changes between 5 kpa and 5 MPa are plotted and used in the analysis. Inset shows the P axes (red) and T axes (blue) of the focal mechanisms for these events. (b) Post Landers earthquakes, from just after the 28 June 1992 M6.5 Big Bear aftershock until 28 June Only events with differential stress changes between 5 kpa and 5 MPa are used, and all aftershocks <5 km from the modeled main shock source are excluded, leaving 725 events in the catalog. Inset shows the P axes (red) and T axes (blue) of the focal mechanisms for these events. Blue lines represent the Landers surface rupture and model fault planes for Joshua Tree and Big Bear. SAF, San Andreas Fault. (c) Pre Landers stress field, inverted from the events in Figure 1a. Arrows show the orientation and magnitude of the horizontal principal axes of the normalized deviatoric stress tensor. The 8 km 8 km inversion grid has been decimated for plotting purposes. The inset shows stress orientation around Cajon Pass, the only location where the change in stress orientation is not fully captured by the decimated stress field. (d) Main shock static stress changes, from the combined Landers, Joshua Tree, and Big Bear earthquakes. Arrows show the orientation and magnitude of the horizontal principal axes of the normalized deviatoric stress tensor. 4of10

5 Figure 2. (a) Misfits of the pre Landers earthquakes, with respect to the inverted pre Landers stress field (Figure 1c) and modeled static stress change field (Figure 1d). Misfit is calculated as the angle between the rake vector and the shear stress vector projected on the fault plane. (b) Misfits of the post Landers earthquakes, with respect to the inverted prestress field and modeled static stress change field. Note that most post Landers earthquakes are consistent with both the prestress and the main shock static stress change, and very few appear to be independent of the prestress. stress change, indicating that the post Landers catalog is dominated by aftershocks that are consistent with triggering by the coseismic static stress changes (Figure 2b). The mechanisms of the post Landers events also exhibit low misfits with respect to the prestress and appear more consistent with the prestress than with the Landers induced stress changes. So, while the aftershock mechanisms are clearly responding to the Landers stress change, they remain aligned with the pre Landers stress field. For comparison, the misfits of the pre Landers events are low with respect to the prestress, as expected, and span the full misfit range with respect to the Landers stress changes (Figure 2a). The misfit of the post Landers events to the prestress appears slightly larger than the misfit of the pre Landers events, which is expected because the stress field was fit to the pre Landers focal mechanisms. [17] Post Landers events in regions experiencing a MPa stress change have an average misfit to the modeled Landers static stress change that is less than the pre Landers events (Figure 3a). This misfit decrease is outside of what would occur by chance, at >99% confidence, implying that the post Landers events with MPa stress change are dominated by aftershocks aligned with the static stress change. This stress threshold is an order of magnitude greater than the estimated tectonic loading for the duration of the catalog, so the aftershocks should not be affected by the tectonic loading stress. However, the aftershocks show only a modest increase of 5 38 in average misfit with respect to the prestress, compared to the pre Landers mechanisms used to fit the prestress (Figure 3b). For most subcatalogs, the average misfit change is indistinguishable from the reshuffled data at 95% confidence. For subcatalogs with stress changes of MPa and 1 5of10

6 Figure 3. (a) Difference in mechanism misfit before versus after Landers, with respect to the modeled Landers static stress change (Figure 1d), as a function of differential stress change. The misfit change is defined as the average misfit of events after Landers minus the average misfit of events before Landers. Black lines represent the results for the real data, with the horizontal error bars representing the range of differential stress values and the vertical error bars representing nodal plane uncertainty. Shading shows the range of results expected from random chance, found by reshuffling the data. The middle 90% of 3000 reshufflings are shown. The numbers at the top of each bin indicate the number of pre Landers and post Landers events in that bin. (b) Difference in mechanism misfit before versus after Landers, with respect to the inverted pre Landers stress field (Figure 1c), as a function of differential stress change. Black lines represent the real data, with the horizontal error bars representing the range of differential stress values and the vertical error bars representing nodal plane uncertainty and different stress inversion grid configurations. Shading shows the range of results expected from random chance, found by reshuffling the data. The middle 90% of 4000 reshuffling results are shown. 2 MPa, the average misfit change is different from that of the reshuffled data at 95% confidence, although, when the nodal plane ambiguity and inversion grid effects are accounted for, the 90% confidence ranges overlap. These two bins together contain only 7% of the aftershock catalog, so this marginally significant movement away from the pre Landers stress tensor is limited to a small number of aftershocks. The majority of triggered aftershocks are well aligned with the prestress, inconsistent with the predictions of the heterogeneous stress hypothesis. [18] I repeated this test, using only events within 10 km of the San Andreas Fault, to test whether stress heterogeneity along this major fault could have been missed in the analysis of the larger data set. The results are quite similar. The post Landers events experiencing a 0.1 MPa stress change have an average misfit to the modeled stress change that is less than the pre Landers events (Figure 4a). This misfit decrease is outside of what would occur by chance at >99% confidence, implying that the post Landers events along the San Andreas are dominated by aftershocks. These aftershocks show only a modest increase of 3 17 in average 6of10

7 Figure 4. As in Figure 3, but for earthquakes within 10 km of the San Andreas Fault. misfit with respect to the prestress, compared to the pre Landers mechanisms (Figure 4b). This average misfit change is indistinguishable from the reshuffled data at >95% confidence. Therefore, I conclude that the heterogeneous stress hypothesis can be rejected near the San Andreas Fault as well as for the rest of southern California. [19] Interestingly, the triggered aftershocks are simultaneously consistent with two stress fields: the pre Landers stress field and the Landers induced stress change field. One possible explanation is that aftershocks may preferentially occur in locations where the main shock stress change aligns with the background stress. This idea is the basis for modeling aftershock triggering using Coulomb stress changes on optimally oriented planes [e.g., King et al., 1994]. However, as McKenzie [1969] demonstrated, a single focal mechanism may be consistent with a range of different stress tensors if the fault is not assumed to be optimally oriented. Therefore, another possible explanation is that the stress fields are different and the triggered aftershocks represent the subset of focal mechanisms that are consistent with both stress fields. [20] The data suggest a combination of these two explanations. The background stress and the Landers induced stress change are in general not very similar, as measured by the normalized tensor dot product (Figure 5). The post Landers events preferentially occur where the two stress tensors are somewhat more similar. The distribution of the normalized tensor dot product is shifted toward higher values when sampled at the locations of the aftershocks compared to the locations of the pre Landers earthquakes or a uniform sample. However, many aftershocks occur in locations where the two stress tensors are dissimilar: about two thirds occur where the normalized tensor dot product is less than 0.5, and one third occurs where the tensor dot product is negative. Figure 6 shows a set of aftershocks that occurred where the two stress fields are quite dissimilar, with a normalized tensor dot product of 0.44, yet the orientation of the shear stress projected onto the fault planes is similar. 5. Discussion [21] The heterogeneous stress hypothesis predicts that triggered aftershocks should occur in volumes of the crust 7of10

8 Figure 5. The similarity between the pre Landers stress field and the Landers induced stress change field, measured by the normalized tensor dot product between the deviatoric parts of the stress tensors. The solid black line shows the cumulative distribution of the similarity sampled at each grid node used in the stress inversion. The dashed and shaded lines show the similarity sampled at the locations of the pre Landers and post Landers earthquakes, respectively. where stress is aligned with the Landers stress change, whereas the pre Landers events should occur in volumes where stress is aligned with the tectonic loading stress. The results, however, show that the triggered aftershocks are aligned with the same stress field as the pre Landers mechanisms. One might argue that, in a heterogeneous stress field, the aftershocks might preferentially align with the tectonic loading if years of loading have brought many faults in well oriented stress volumes near failure, while also perhaps moving faults in other stress volumes away from failure. However, the accumulated stress owing to past tectonic loading, which is a persistent spatially coherent stress field, should actually be considered part of the uniform background stress. If the background stress is large enough to have a significant impact on the orientations of faults near failure, the heterogeneous part of the stress field must be small in comparison, implying a relatively uniform stress field. [22] The simplest explanation of the focal mechanism catalog is that the prestress is the persistent spatially coherent background stress, which is large enough compared to any heterogeneity that it can be successfully retrieved using standard stress inversion techniques [e.g., Michael, 1984; Hardebeck and Michael, 2006]. The aftershocks are consistent with both the background stress and the Landers induced stress change, in part because they preferentially occur in regions where the two stress fields are more similar and in part because a single focal mechanism can be consistent with a range of different stress tensors. This requires failure of nonoptimally oriented planes, such as preexisting faults, which are not included in Smith and Dieterich s [2010] modeling. These faults failed because they were well aligned within the dominant background stress field and then received a small additional stress load from the Landers main shock. [23] The heterogeneous stress model, if it applies anywhere, should apply to the complex active fault system in southern California. Therefore, the counterexample provided by the Landers earthquake is strong evidence that the heterogeneous stress model is not widely applicable. The Landers counterexample also disproves the hypothesis that focal mechanisms are biased toward the stress loading direction and implies that focal mechanism stress inversion techniques [e.g., Michael, 1984], when used properly, correctly return the background stress orientations. Reported stress rotations due to large earthquakes, including in the near field of the Landers main shock [Hauksson, 1994], are therefore not sampling artifacts and imply true rotations of the background stress. [24] The only convincing focal mechanism evidence for heterogeneous stress comes from a few data sets, typically aftershock sequences, where the focal mechanism variability is too large to be explained by a uniform stress tensor. The best studied example is the aftershock sequence of the 1989 M6.9 Loma Prieta earthquake [e.g., Michael et al., 1990], which exhibited right lateral, left lateral, thrust, and normal faulting mechanisms. The heterogeneity is restricted to the events along and directly above the main shock rupture surface, whereas the rest of the aftershocks are more homogeneous [Michael et al., 1990; Hardebeck, 2006]. Other earthquakes, such as the Landers main shock [Hauksson, 1994], exhibit aftershock heterogeneity along only some parts of the rupture surface. This suggests that, although large earthquakes can produce heterogeneous stress, the effects are very localized. The heterogeneity also cannot be long lived; otherwise stress heterogeneity would be apparent along all major faults. This suggests either that the aftershocks significantly reduce the stress heterogeneity or that some other stress relaxation mechanism is in effect. [25] The general homogeneity of the stress field also has important implications for the dynamics of earthquake rupture. Assuming that earthquake slip distributions are as heterogeneous as they appear in kinematic slip models, the initial stress distribution, the fault constitutive properties, or both must be heterogeneous [e.g., Aagaard and Heaton, 2008]. Low heterogeneity in the stress field implies that significant spatial variation is required in the fault constitutive properties. 6. Conclusions [26] I used earthquakes before and after the 1992 M7.3 Landers, California, earthquake to test the hypothesis that the crustal stress field contains large amplitude, small lengthscale heterogeneity. The heterogeneous stress hypothesis predicts that triggered aftershocks should occur in volumes of the crust where stress is aligned with the Landers stress change, whereas the pre Landers events should occur in volumes where stress is aligned with the tectonic loading stress. However, both the pre Landers and aftershock focal mechanisms are well aligned with the inverted prestress, indicating that they are in fact sampling the same stress field. The pre Landers stress field and the Landers stress change field are quite different over much of the region, ruling out coincidental alignment. The alignment of the aftershock 8of10

9 Figure 6. Examples of Landers aftershocks with focal mechanisms that are consistent with both the pre Landers stress field and the Landers induced static stress change. These mechanisms occurred in a region where the two stress tensors are dissimilar, with a normalized tensor dot product of The open red and blue symbols show the pre Landers stress and Landers stress change tensors, respectively. The black and shaded lines show the fault plane in the lower and upper hemispheres, respectively, and the rake is shown as a solid black or shaded circle. The solid red and blue circles show the projected shear stress direction of the two stress tensors onto the fault plane, with lighter colors indicating the upper hemisphere. focal mechanisms with the prestress is inconsistent with the prediction of the heterogeneous stress model that the pre Landers events and triggered aftershocks should sample crustal volumes with different stress orientations. Therefore, the inverted pre Landers stress must be the persistent background stress field. The aftershocks are consistent with both the background stress and the Landers induced stress change, in part because they preferentially occur in regions where the two stress fields are more similar and in part because a single focal mechanism can be consistent with a range of different stress tensors. If heterogeneous stress were common in the crust, it should be expected in southern California, given the complexity of the fault system. The Landers counterexample is therefore strong evidence that the heterogeneous stress model is not widely applicable. [27] Acknowledgments. IamgratefultoAndyMichaelandBrad Aagaard for constructive reviews of an earlier version of this manuscript, to the associate editor and two anonymous reviewers for their comments, and to Debbie Smith and Tom Heaton for many interesting discussions about stress heterogeneity. I thank the Southern California Seismic Network (SCSN) and the Southern California Earthquake Data Center (SCEDC) for collecting and providing access to the earthquake first motion data used to generate the focal mechanism catalog. 9of10

10 References Aagaard, B. T., and T. H. Heaton (2008), Constraining fault constitutive behavior with slip and stress heterogeneity, J. Geophys. Res., 113, B04301, doi: /2006jb Hardebeck, J. L. (2006), Homogeneity of small scale earthquake faulting, stress and fault strength, Bull. Seismol. Soc. Am., 96, Hardebeck, J. L., and E. Hauksson (2001), Crustal stress field in southern California and its implications for fault mechanics, J. Geophys. Res., 106(B10), 21,859 21,882. Hardebeck, J. L., and A. J. Michael (2006), Damped regional scale stress inversions: Methodology and examples for southern California and the Coalinga aftershock sequence, J. Geophys. Res., 111, B11310, doi: /2005jb Hardebeck, J. L., and P. M. Shearer (2002), A new method for determining first motion focal mechanisms, Bull. Seismol. Soc. Am., 92, Hardebeck, J. L., J. J. Nazareth, and E. Hauksson (1998), The static stress change triggering model: Constraints from two southern California aftershock sequences, J. Geophys. Res., 103(B10), 24,427 24,438. Hauksson, E. (1994), State of stress from focal mechanisms before and after the 1992 Landers earthquake sequence, Bull.Seismol.Soc.Am., 84, Hauksson, E., L. M. Jones, K. Hutton, and D. Eberhart Phillips (1993), The 1992 Landers earthquake sequence: Seismological observations, J. Geophys. Res., 98(B11), 19,835 19,858. Hickman, S., and M. Zoback (2004), Stress orientation and magnitudes in the SAFOD pilot hole, Geophys. Res. Lett., 31, L15S12, doi: / 2004GL Jackson, D. D., Z. K. Shen, D. Potter, X. Ge, and L. Sung (1997), Southern California deformation, Science, 277, King, G. C. P., R. S. Stein, and J. Lin (1994), Static stress changes and the triggering of earthquakes, Bull. Seismol. Soc. Am., 84, Liu, P., S. Custódio, and R. J. Archuleta (2006), Kinematic inversion of the 2004 M 6.0 Parkfield earthquake including an approximation to site effects, Bull. Seismol. Soc. Am., 96, S143 S158. McGill, S. F., and C. M. Rubin (1999), Surficial slip distribution on the central Emerson fault during the June 28, 1992, Landers earthquake, California, J. Geophys. Res., 104(B3), McKenzie, D. P. (1969), The relation between fault plane solutions for earthquakes and the directions of the principal stresses, Bull. Seismol. Soc. Am., 59, Michael, A. J. (1984), Determination of stress from slip data: Faults and folds, J. Geophys. Res., 89(B13), 11,517 11,526. Michael, A. J. (1987), Stress rotation during the Coalinga aftershock sequence (USA), J. Geophys. Res., 92(B8), Michael, A. J., W. L. Ellsworth, and D. H. Oppenheimer (1990), Coseismic stress changes induced by the 1989 Loma Prieta, California earthquake, Geophys. Res. Lett., 17(9), Okada, Y. (1992), Internal deformation due to shear and tensile faults in a half space, Bull. Seismol. Soc. Am., 82, Page, M. T., S. Custódio, R. J. Archuleta, and J. M. Carlson (2009), Constraining earthquake source inversions with GPS data: 1. Resolution based removal of artifacts, J. Geophys. Res., 114, B01314, doi: / 2007JB Rivera, L., and H. Kanamori (2002), Spatial heterogeneity of tectonic stress and friction in the crust, Geophys. Res. Lett., 29(6), 1088, doi: / 2001GL Smith, D. E. (2006), A new paradigm for interpreting stress inversions from focal mechanisms: How 3D stress heterogeneity biases the inversions toward the stress rate, Ph. D. dissertation, Calif. Inst. of Technol, Pasadena, Cal. Smith, D. E., and J. H. Dieterich (2010), Aftershock sequences modeled in 3 D stress heterogeneity and rate state seismicity equations: Implications for crustal stress estimation, Pure Appl. Geophys., 167, , doi: /s Wald, D. J., and T. H. Heaton (1994), Spatial and temporal distribution of slip for the 1992 Landers, California, earthquake, Bull. Seismol. Soc. Am., 84, Wilde, M., and J. Stock (1997), Compression directions in southern California (from Santa Barbara to Los Angeles basin) obtained from borehole breakouts, J. Geophys. Res., 102(B3), J. Hardebeck, U.S. Geological Survey, 345 Middlefield Rd., MS 977, Menlo Park, CA 94025, USA. (jhardebeck@usgs.gov) 10 of 10

by Deborah Elaine Smith * and Thomas H. Heaton Introduction

by Deborah Elaine Smith * and Thomas H. Heaton Introduction Bulletin of the Seismological Society of America, Vol. 11, No. 3, pp. 1396 1421, June 211, doi: 1.1785/12158 Models of Stochastic, Spatially Varying Stress in the Crust Compatible with Focal-Mechanism

More information

Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications

Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications Bulletin of the Seismological Society of America, Vol. 85, No. 5, pp. 1513-1517, October 1995 Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications by

More information

Stress Orientations Obtained from Earthquake Focal Mechanisms: What Are Appropriate Uncertainty Estimates?

Stress Orientations Obtained from Earthquake Focal Mechanisms: What Are Appropriate Uncertainty Estimates? Bulletin of the Seismological Society of America, 91, 2, pp. 250 262, April 2001 Stress Orientations Obtained from Earthquake Focal Mechanisms: What Are Appropriate Uncertainty Estimates? by Jeanne L.

More information

Stress orientations at intermediate angles to the San Andreas Fault, California

Stress orientations at intermediate angles to the San Andreas Fault, California JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2004jb003239, 2004 Stress orientations at intermediate angles to the San Andreas Fault, California Jeanne L. Hardebeck and Andrew J. Michael U.S.

More information

Earth Tides Can Trigger Shallow Thrust Fault Earthquakes

Earth Tides Can Trigger Shallow Thrust Fault Earthquakes Earth Tides Can Trigger Shallow Thrust Fault Earthquakes Elizabeth S. Cochran, 1 * John E. Vidale, 1 Sachiko Tanaka 2 1 Department of Earth and Space Sciences and Institute of Geophysics and Planetary

More information

Aftershock Sequences Modeled with 3-D Stress Heterogeneity and Rate-State Seismicity Equations: Implications for Crustal Stress Estimation

Aftershock Sequences Modeled with 3-D Stress Heterogeneity and Rate-State Seismicity Equations: Implications for Crustal Stress Estimation Pure Appl. Geophys. 167 (2010), 1067 1085 Ó 2010 The Author(s) This article is published with open access at Springerlink.com DOI 10.1007/s00024-010-0093-1 Pure and Applied Geophysics Aftershock Sequences

More information

In geophysics there is increasing interest in modeling spatially nonuniform stress,

In geophysics there is increasing interest in modeling spatially nonuniform stress, Chapter 1. Introduction I-1 In geophysics there is increasing interest in modeling spatially nonuniform stress, i.e., spatially heterogeneous stress, on 2D planar faults as a means of explaining a variety

More information

Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation.

Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. In the process zone, stress amplitudes are poorly determined and much

More information

Onto what planes should Coulomb stress perturbations be resolved?

Onto what planes should Coulomb stress perturbations be resolved? JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2004jb003356, 2005 Onto what planes should Coulomb stress perturbations be resolved? Sandy Steacy, Süleyman S. Nalbant, and John McCloskey School

More information

Regional tectonic stress near the San Andreas fault in central and southern California

Regional tectonic stress near the San Andreas fault in central and southern California GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L15S11, doi:10.1029/2003gl018918, 2004 Regional tectonic stress near the San Andreas fault in central and southern California J. Townend School of Earth Sciences,

More information

to: Interseismic strain accumulation and the earthquake potential on the southern San

to: Interseismic strain accumulation and the earthquake potential on the southern San Supplementary material to: Interseismic strain accumulation and the earthquake potential on the southern San Andreas fault system by Yuri Fialko Methods The San Bernardino-Coachella Valley segment of the

More information

COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES

COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES Madlazim Physics Department, Faculty Mathematics and Sciences of Surabaya State University (UNESA) Jl. Ketintang, Surabaya 60231, Indonesia. e-mail:

More information

Measurements in the Creeping Section of the Central San Andreas Fault

Measurements in the Creeping Section of the Central San Andreas Fault Measurements in the Creeping Section of the Central San Andreas Fault Introduction Duncan Agnew, Andy Michael We propose the PBO instrument, with GPS and borehole strainmeters, the creeping section of

More information

Predicted reversal and recovery of surface creep on the Hayward fault following the 1906 San Francisco earthquake

Predicted reversal and recovery of surface creep on the Hayward fault following the 1906 San Francisco earthquake GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L19305, doi:10.1029/2008gl035270, 2008 Predicted reversal and recovery of surface creep on the Hayward fault following the 1906 San Francisco earthquake D. A. Schmidt

More information

Does Aftershock Duration Scale With Mainshock Size?

Does Aftershock Duration Scale With Mainshock Size? GEOPHYSICAL RESEARCH LETTERS, VOL.???, NO., PAGES 1 16, Does Aftershock Duration Scale With Mainshock Size? A. Ziv A. Ziv, Ben-Gurion University of the Negev, Beer-Sheva, 84105, Israel. (e-mail: zival@bgu.ac.il)

More information

SUPPLEMENTAL INFORMATION

SUPPLEMENTAL INFORMATION GSA DATA REPOSITORY 2013310 A.M. Thomas et al. MOMENT TENSOR SOLUTIONS SUPPLEMENTAL INFORMATION Earthquake records were acquired from the Northern California Earthquake Data Center. Waveforms are corrected

More information

JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. 0, XXXX, doi: /2001jb000678, 2002

JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. 0, XXXX, doi: /2001jb000678, 2002 JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. 0, XXXX, doi:10.1029/2001jb000678, 2002 A strong correlation between induced peak dynamic Coulomb stress change from the 1992 M7.3 Landers, California, earthquake

More information

Homogeneous vs. realistic heterogeneous material-properties in subduction zone models: Coseismic and postseismic deformation

Homogeneous vs. realistic heterogeneous material-properties in subduction zone models: Coseismic and postseismic deformation Homogeneous vs. realistic heterogeneous material-properties in subduction zone models: Coseismic and postseismic deformation T. Masterlark 1, C. DeMets 2, H.F. Wang 2, O. S nchez 3, and J. Stock 4 1 US

More information

Stress triggering and earthquake probability estimates

Stress triggering and earthquake probability estimates JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2003jb002437, 2004 Stress triggering and earthquake probability estimates Jeanne L. Hardebeck 1 Institute for Geophysics and Planetary Physics, Scripps

More information

Earthquake Stress Drops in Southern California

Earthquake Stress Drops in Southern California Earthquake Stress Drops in Southern California Peter Shearer IGPP/SIO/U.C. San Diego September 11, 2009 Earthquake Research Institute Lots of data for big earthquakes (rupture dimensions, slip history,

More information

Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment

Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment Abstract D. Weatherley University of Queensland Coulomb stress change analysis has been applied in

More information

A mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements

A mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L15S13, doi:10.1029/2004gl019521, 2004 A mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements Jean Chéry Laboratoire Dynamique de la

More information

Seismic and aseismic processes in elastodynamic simulations of spontaneous fault slip

Seismic and aseismic processes in elastodynamic simulations of spontaneous fault slip Seismic and aseismic processes in elastodynamic simulations of spontaneous fault slip Most earthquake simulations study either one large seismic event with full inertial effects or long-term slip history

More information

Kinematics of the Southern California Fault System Constrained by GPS Measurements

Kinematics of the Southern California Fault System Constrained by GPS Measurements Title Page Kinematics of the Southern California Fault System Constrained by GPS Measurements Brendan Meade and Bradford Hager Three basic questions Large historical earthquakes One basic question How

More information

Automatic Moment Tensor Analyses, In-Situ Stress Estimation and Temporal Stress Changes at The Geysers EGS Demonstration Project

Automatic Moment Tensor Analyses, In-Situ Stress Estimation and Temporal Stress Changes at The Geysers EGS Demonstration Project PROCEEDINGS, 42nd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 13-15, 2017 SGP-TR-212 Automatic Moment Tensor Analyses, In-Situ Stress Estimation and

More information

arxiv: v1 [physics.geo-ph] 30 Apr 2013

arxiv: v1 [physics.geo-ph] 30 Apr 2013 arxiv:.v [physics.geo-ph] Apr Sensitivity of the Static Earthquake Triggering Mechanism to Elastic Heterogeneity and Main Event Slip Musa Maharramov musa@sep.stanford.edu September, Abstract This paper

More information

Synthetic Seismicity Models of Multiple Interacting Faults

Synthetic Seismicity Models of Multiple Interacting Faults Synthetic Seismicity Models of Multiple Interacting Faults Russell Robinson and Rafael Benites Institute of Geological & Nuclear Sciences, Box 30368, Lower Hutt, New Zealand (email: r.robinson@gns.cri.nz).

More information

Jack Loveless Department of Geosciences Smith College

Jack Loveless Department of Geosciences Smith College Geodetic constraints on fault interactions and stressing rates in southern California Jack Loveless Department of Geosciences Smith College jloveless@smith.edu Brendan Meade Department of Earth & Planetary

More information

Short Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood, California, Earthquake

Short Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood, California, Earthquake Bulletin of the Seismological Society of America, Vol. 100, No. 6, pp. 3269 3277, December 2010, doi: 10.1785/0120100087 Short Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood,

More information

High Resolution Imaging of Fault Zone Properties

High Resolution Imaging of Fault Zone Properties Annual Report on 1998-99 Studies, Southern California Earthquake Center High Resolution Imaging of Fault Zone Properties Yehuda Ben-Zion Department of Earth Sciences, University of Southern California

More information

FOCAL MECHANISM DETERMINATION OF LOCAL EARTHQUAKES IN MALAY PENINSULA

FOCAL MECHANISM DETERMINATION OF LOCAL EARTHQUAKES IN MALAY PENINSULA FOCAL MECHANISM DETERMINATION OF LOCAL EARTHQUAKES IN MALAY PENINSULA Siti Norbaizura MAT SAID Supervisor: Tatsuhiko HARA MEE10505 ABSTRACT Since November 30, 2007, small local earthquakes have been observed

More information

A TESTABLE FIVE-YEAR FORECAST OF MODERATE AND LARGE EARTHQUAKES. Yan Y. Kagan 1,David D. Jackson 1, and Yufang Rong 2

A TESTABLE FIVE-YEAR FORECAST OF MODERATE AND LARGE EARTHQUAKES. Yan Y. Kagan 1,David D. Jackson 1, and Yufang Rong 2 Printed: September 1, 2005 A TESTABLE FIVE-YEAR FORECAST OF MODERATE AND LARGE EARTHQUAKES IN SOUTHERN CALIFORNIA BASED ON SMOOTHED SEISMICITY Yan Y. Kagan 1,David D. Jackson 1, and Yufang Rong 2 1 Department

More information

Supporting Information for Break of slope in earthquake-size distribution reveals creep rate along the San Andreas fault system

Supporting Information for Break of slope in earthquake-size distribution reveals creep rate along the San Andreas fault system GEOPHYSICAL RESEARCH LETTERS Supporting Information for Break of slope in earthquake-size distribution reveals creep rate along the San Andreas fault system Inessa Vorobieva, 1,2 Peter Shebalin, 1,2 Clément

More information

Earthquake stress drop estimates: What are they telling us?

Earthquake stress drop estimates: What are they telling us? Earthquake stress drop estimates: What are they telling us? Peter Shearer IGPP/SIO/U.C. San Diego October 27, 2014 SCEC Community Stress Model Workshop Lots of data for big earthquakes (rupture dimensions,

More information

UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions

UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions UCERF3 Task R- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions Bruce E. Shaw Lamont Doherty Earth Observatory, Columbia University Statement of the Problem In UCERF Magnitude-Area

More information

Ground displacement in a fault zone in the presence of asperities

Ground displacement in a fault zone in the presence of asperities BOLLETTINO DI GEOFISICA TEORICA ED APPLICATA VOL. 40, N. 2, pp. 95-110; JUNE 2000 Ground displacement in a fault zone in the presence of asperities S. SANTINI (1),A.PIOMBO (2) and M. DRAGONI (2) (1) Istituto

More information

UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions

UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions Bruce E. Shaw Lamont Doherty Earth Observatory, Columbia University Statement of the Problem In UCERF2

More information

EARTHQUAKE LOCATIONS INDICATE PLATE BOUNDARIES EARTHQUAKE MECHANISMS SHOW MOTION

EARTHQUAKE LOCATIONS INDICATE PLATE BOUNDARIES EARTHQUAKE MECHANISMS SHOW MOTION 6-1 6: EARTHQUAKE FOCAL MECHANISMS AND PLATE MOTIONS Hebgen Lake, Montana 1959 Ms 7.5 1 Stein & Wysession, 2003 Owens Valley, California 1872 Mw ~7.5 EARTHQUAKE LOCATIONS INDICATE PLATE BOUNDARIES EARTHQUAKE

More information

Lecture 20: Slow Slip Events and Stress Transfer. GEOS 655 Tectonic Geodesy Jeff Freymueller

Lecture 20: Slow Slip Events and Stress Transfer. GEOS 655 Tectonic Geodesy Jeff Freymueller Lecture 20: Slow Slip Events and Stress Transfer GEOS 655 Tectonic Geodesy Jeff Freymueller Slow Slip Events From Kristine Larson What is a Slow Slip Event? Slip on a fault, like in an earthquake, BUT

More information

On the nucleation of creep and the interaction between creep and seismic slip on rate- and state-dependent faults

On the nucleation of creep and the interaction between creep and seismic slip on rate- and state-dependent faults Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L15303, doi:10.1029/2007gl030337, 2007 On the nucleation of creep and the interaction between creep and seismic slip on rate- and state-dependent

More information

Time-Dependent Viscoelastic Stress Transfer and Earthquake Triggering

Time-Dependent Viscoelastic Stress Transfer and Earthquake Triggering Time-Dependent Viscoelastic Stress Transfer and Earthquake Triggering Jian Lin Department of Geology and Geophysics, Woods Hole Oceanographic Institution, Woods Hole, MA 02543 USA, jlin@whoi.edu tel 508

More information

California foreshock sequences suggest aseismic triggering process

California foreshock sequences suggest aseismic triggering process GEOPHYSICAL RESEARCH LETTERS, VOL. 0, 1, doi:10.100/grl.50, 013 California foreshock sequences suggest aseismic triggering process Xiaowei Chen 1 and Peter M. Shearer 1 Received January 013; revised 9

More information

of other regional earthquakes (e.g. Zoback and Zoback, 1980). I also want to find out

of other regional earthquakes (e.g. Zoback and Zoback, 1980). I also want to find out 4. Focal Mechanism Solutions A way to investigate source properties of the 2001 sequence is to attempt finding well-constrained focal mechanism solutions to determine if they are consistent with those

More information

Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source

Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source Nevada Bureau of Mines and Geology Special Publication 36 Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source by Mendoza, C. 1 and Hartzell S. 2 1 Centro de Geociencias, Universidad

More information

Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013

Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013 Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013 Abstract Earthquakes do not fit into the class of models we discussed in Physics 219B. Earthquakes

More information

Fault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction

Fault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction Fault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction H. Sekiguchi Disaster Prevention Research Institute, Kyoto University, Japan Blank Line 9 pt Y. Kase Active Fault and Earthquake

More information

State of Stress from Focal Mechanisms Before and After the 1992 Landers Earthquake Sequence

State of Stress from Focal Mechanisms Before and After the 1992 Landers Earthquake Sequence Bulletin of the Seismological Society of America, Vol. 84, No. 3, pp. 917-934, June 1994 State of Stress from Focal Mechanisms Before and After the 1992 Landers Earthquake Sequence by Egill Hauksson Abstract

More information

Off-fault aftershocks of the 2005 West Off Fukuoka Prefecture Earthquake: Reactivation of a structural boundary?

Off-fault aftershocks of the 2005 West Off Fukuoka Prefecture Earthquake: Reactivation of a structural boundary? LETTER Earth Planets Space, 58, 81 86, 26 Off-fault aftershocks of the 25 West Off Fukuoka Prefecture Earthquake: Reactivation of a structural boundary? Kazutoshi Imanishi 1, Yasuto Kuwahara 1, and Yoshikatsu

More information

Today: Basic regional framework. Western U.S. setting Eastern California Shear Zone (ECSZ) 1992 Landers EQ 1999 Hector Mine EQ Fault structure

Today: Basic regional framework. Western U.S. setting Eastern California Shear Zone (ECSZ) 1992 Landers EQ 1999 Hector Mine EQ Fault structure Today: Basic regional framework Western U.S. setting Eastern California Shear Zone (ECSZ) 1992 Landers EQ 1999 Hector Mine EQ Fault structure 1 2 Mojave and Southern Basin and Range - distribution of strike-slip

More information

Testing the stress shadow hypothesis

Testing the stress shadow hypothesis JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2004jb003277, 2005 Testing the stress shadow hypothesis Karen R. Felzer and Emily E. Brodsky Department of Earth and Space Sciences, University of

More information

M 7.0 earthquake recurrence on the San Andreas fault from a stress renewal model

M 7.0 earthquake recurrence on the San Andreas fault from a stress renewal model Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2006jb004415, 2006 M 7.0 earthquake recurrence on the San Andreas fault from a stress renewal model Tom Parsons 1 Received

More information

On the origin of diverse aftershock mechanisms 1989 Loma Prieta earthquake. following the

On the origin of diverse aftershock mechanisms 1989 Loma Prieta earthquake. following the Geophys. J. Int. (1997) 128,557-570 On the origin of diverse aftershock mechanisms 1989 Loma Prieta earthquake following the D. Kilb,' M. Ellis,' J. Gomberg2 and S. Davis3 Center for Earthquake Research

More information

3D Finite Element Modeling of fault-slip triggering caused by porepressure

3D Finite Element Modeling of fault-slip triggering caused by porepressure 3D Finite Element Modeling of fault-slip triggering caused by porepressure changes Arsalan Sattari and David W. Eaton Department of Geoscience, University of Calgary Suary We present a 3D model using a

More information

The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets

The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets P. Banerjee, 1 F. F. Pollitz, 2 R. Bürgmann 3 * 1 Wadia Institute of Himalayan Geology, Dehra Dun, 248001, India. 2

More information

Earthquake Doublet Sequences: Evidence of Static Triggering in the Strong Convergent Zones of Taiwan

Earthquake Doublet Sequences: Evidence of Static Triggering in the Strong Convergent Zones of Taiwan Terr. Atmos. Ocean. Sci., Vol. 19, No. 6, 589-594, December 2008 doi: 10.3319/TAO.2008.19.6.589(PT) Earthquake Doublet Sequences: Evidence of Static Triggering in the Strong Convergent Zones of Taiwan

More information

Aftershocks and Pore Fluid Diffusion Following the 1992 Landers Earthquake

Aftershocks and Pore Fluid Diffusion Following the 1992 Landers Earthquake The University of San Francisco USF Scholarship: a digital repository @ Gleeson Library Geschke Center Nursing and Health Professions Faculty Research and Publications School of Nursing and Health Professions

More information

Focused Observation of the San Andreas/Calaveras Fault intersection in the region of San Juan Bautista, California

Focused Observation of the San Andreas/Calaveras Fault intersection in the region of San Juan Bautista, California Focused Observation of the San Andreas/Calaveras Fault intersection in the region of San Juan Bautista, California Clifford Thurber, Dept. of Geology and Geophysics, Univ. Wisconsin-Madison Mick Gladwin,

More information

Slab pull, slab weakening, and their relation to deep intra-slab seismicity

Slab pull, slab weakening, and their relation to deep intra-slab seismicity GEOPHYSICAL RESEARCH LETTERS, VOL. 32, L14305, doi:10.1029/2005gl022922, 2005 Slab pull, slab weakening, and their relation to deep intra-slab seismicity Susan L. Bilek Earth and Environmental Science

More information

Rate and State-Dependent Friction in Earthquake Simulation

Rate and State-Dependent Friction in Earthquake Simulation Rate and State-Dependent Friction in Earthquake Simulation Zac Meadows UC Davis - Department of Physics Summer 2012 REU September 3, 2012 Abstract To better understand the spatial and temporal complexity

More information

Velocity contrast along the Calaveras fault from analysis of fault zone head waves generated by repeating earthquakes

Velocity contrast along the Calaveras fault from analysis of fault zone head waves generated by repeating earthquakes Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L01303, doi:10.1029/2007gl031810, 2008 Velocity contrast along the Calaveras fault from analysis of fault zone head waves generated by

More information

ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL

ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL 1 Best Practices in Physics-based Fault Rupture Models for Seismic Hazard Assessment of Nuclear ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL Hideo AOCHI

More information

Foreshocks, Aftershocks, and Earthquake Probabilities: Accounting for the Landers Earthquake

Foreshocks, Aftershocks, and Earthquake Probabilities: Accounting for the Landers Earthquake Bulletin of the Seismological Society of America, Vol. 84, No. 3, pp. 892-899, June 1994 Foreshocks, Aftershocks, and Earthquake Probabilities: Accounting for the Landers Earthquake by Lucile M. Jones

More information

Bulletin of the Seismological Society of America, Vol. 81, No. 6, pp , December 1991

Bulletin of the Seismological Society of America, Vol. 81, No. 6, pp , December 1991 Bulletin of the Seismological Society of America, Vol. 81, No. 6, pp. 2310-2319, December 1991 AFTERSHOCK SEQUENCE OF THE 3 DECEMBER 1988 PASADENA EARTHQUAKE BY Kuo-FoNG MA AND HIROO KANAMORI ABSTRACT

More information

State of Stress in Seismic Gaps Along the SanJacinto Fault

State of Stress in Seismic Gaps Along the SanJacinto Fault ELEVEN State of Stress in Seismic Gaps Along the SanJacinto Fault Hirao Kanamori and Harold Magistrale NTRODUCTON Data from the Southern California Seismic Network have been extensively used to map spatial

More information

Potency-magnitude scaling relations for southern California earthquakes with 1.0 < M L < 7.0

Potency-magnitude scaling relations for southern California earthquakes with 1.0 < M L < 7.0 Geophys. J. Int. (2002) 148, F1 F5 FAST TRACK PAPER Potency-magnitude scaling relations for southern California earthquakes with 1.0 < M L < 7.0 Yehuda Ben-Zion 1, and Lupei Zhu 2 1 Department of Earth

More information

Spatial and temporal stress drop variations in small earthquakes near Parkfield, California

Spatial and temporal stress drop variations in small earthquakes near Parkfield, California Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jb004395, 2007 Spatial and temporal stress drop variations in small earthquakes near Parkfield, California Bettina

More information

Limitations of Earthquake Triggering Models*

Limitations of Earthquake Triggering Models* Limitations of Earthquake Triggering Models* Peter Shearer IGPP/SIO/U.C. San Diego September 16, 2009 Earthquake Research Institute * in Southern California Why do earthquakes cluster in time and space?

More information

External Grant Award Number 04HQGR0058 IMPROVED THREE-DIMENSIONAL VELOCITY MODELS AND EARTHQUAKE LOCATIONS FOR CALIFORNIA

External Grant Award Number 04HQGR0058 IMPROVED THREE-DIMENSIONAL VELOCITY MODELS AND EARTHQUAKE LOCATIONS FOR CALIFORNIA External Grant Award Number 04HQGR0058 IMPROVED THREE-DIMENSIONAL VELOCITY MODELS AND EARTHQUAKE LOCATIONS FOR CALIFORNIA Clifford H. Thurber University of Wisconsin-Madison 1215 W. Dayton St. Madison,

More information

GPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise

GPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise GPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise Strain Analysis Introduction Name: The earthquake cycle can be viewed as a process of slow strain accumulation

More information

! EN! EU! NE! EE.! ij! NN! NU! UE! UN! UU

! EN! EU! NE! EE.! ij! NN! NU! UE! UN! UU A-1 Appendix A. Equations for Translating Between Stress Matrices, Fault Parameters, and P-T Axes Coordinate Systems and Rotations We use the same right-handed coordinate system as Andy Michael s program,

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (21) 183, 433 45 doi: 1.1111/j.1365-246X.21.4745.x Quantifying focal mechanism heterogeneity for fault zones in central and southern California Iain W.

More information

Testing the stress shadow hypothesis

Testing the stress shadow hypothesis Testing the stress shadow hypothesis Karen Felzer Emily Brodsky March 14, 2005 1 Abstract A fundamental question in earthquake physics is whether aftershocks are predominantly triggered by static stress

More information

On rate-state and Coulomb failure models

On rate-state and Coulomb failure models JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 105, NO. B4, PAGES 7857 7871, APRIL 10, 2000 On rate-state and Coulomb failure models J. Gomberg U.S. Geological Survey, Center for Earthquake Research and Information,

More information

Basics of the modelling of the ground deformations produced by an earthquake. EO Summer School 2014 Frascati August 13 Pierre Briole

Basics of the modelling of the ground deformations produced by an earthquake. EO Summer School 2014 Frascati August 13 Pierre Briole Basics of the modelling of the ground deformations produced by an earthquake EO Summer School 2014 Frascati August 13 Pierre Briole Content Earthquakes and faults Examples of SAR interferograms of earthquakes

More information

Triggering of earthquakes during the 2000 Papua New Guinea earthquake sequence

Triggering of earthquakes during the 2000 Papua New Guinea earthquake sequence JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jb004480, 2007 Triggering of earthquakes during the 2000 Papua New Guinea earthquake sequence Sun-Cheon Park 1 and Jim Mori 1 Received 3 May

More information

The role of Coulomb stress changes for injection induced seismicity: The Basel enhanced geothermal system

The role of Coulomb stress changes for injection induced seismicity: The Basel enhanced geothermal system GEOPHYSICAL RESEARCH LETTERS, VOL. 4, 72 77, doi:1.129/212gl54147, 213 The role of Coulomb stress changes for injection induced seismicity: The Basel enhanced geothermal system Flaminia Catalli, 1 Men

More information

The Length to Which an Earthquake will go to Rupture. University of Nevada, Reno 89557

The Length to Which an Earthquake will go to Rupture. University of Nevada, Reno 89557 The Length to Which an Earthquake will go to Rupture Steven G. Wesnousky 1 and Glenn P. Biasi 2 1 Center of Neotectonic Studies and 2 Nevada Seismological Laboratory University of Nevada, Reno 89557 Abstract

More information

Di#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes

Di#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes Bull. Earthq. Res. Inst. Univ. Tokyo Vol. 2+,**0 pp.,/3,00 Di#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes Paul Somerville* and Arben Pitarka

More information

CHANGES OF COULOMB FAILURE STRESS DUE TO DISLOCATIONS DURING STIMULATION OF GPK2 AT SOULTZ-SOUS-FORÊTS

CHANGES OF COULOMB FAILURE STRESS DUE TO DISLOCATIONS DURING STIMULATION OF GPK2 AT SOULTZ-SOUS-FORÊTS PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 CHANGES OF COULOMB FAILURE STRESS DUE TO DISLOCATIONS

More information

Oceanic Detachment Faults Generate Compression in Extension

Oceanic Detachment Faults Generate Compression in Extension GSA Data Repository 2017309 Oceanic Detachment Faults Generate Compression in Extension R. Parnell-Turner 1, R. A. Sohn 1, C. Peirce 2, T. J. Reston 3, C. J. MacLeod 4, R. C. Searle 2 and N. M. Simão 2.

More information

Mechanical origin of aftershocks: Supplementary Information

Mechanical origin of aftershocks: Supplementary Information Mechanical origin of aftershocks: Supplementary Information E. Lippiello Department of Mathematics and Physics, Second University of Naples, Via Vivaldi 43, 81100 Caserta, Italy & Kavli Institute for Theoretical

More information

INGV. Giuseppe Pezzo. Istituto Nazionale di Geofisica e Vulcanologia, CNT, Roma. Sessione 1.1: Terremoti e le loro faglie

INGV. Giuseppe Pezzo. Istituto Nazionale di Geofisica e Vulcanologia, CNT, Roma. Sessione 1.1: Terremoti e le loro faglie Giuseppe Pezzo Istituto Nazionale di Geofisica e Vulcanologia, CNT, Roma giuseppe.pezzo@ingv.it The study of surface deformation is one of the most important topics to improve the knowledge of the deep

More information

Depth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution

Depth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L24301, doi:10.1029/2008gl036249, 2008 Depth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution

More information

Microseismic monitoring is a valuable

Microseismic monitoring is a valuable SPECIAL SECTION: M i c r o s e i s m i c moment tensors: A path to understanding frac growth ADAM BAIG and TED URBANCIC, Engineering Seismology Group Canada monitoring is a valuable tool in understanding

More information

Constraining fault constitutive behavior with slip and stress heterogeneity

Constraining fault constitutive behavior with slip and stress heterogeneity Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2006jb004793, 2008 Constraining fault constitutive behavior with slip and stress heterogeneity B. T. Aagaard 1 and T.

More information

Triggering of Aftershocks of the Japan 2011 Earthquake by Earth Tides

Triggering of Aftershocks of the Japan 2011 Earthquake by Earth Tides Triggering of Aftershocks of the Japan 2011 Earthquake by Earth Tides A. Datta and Kamal* Department of Earth Sciences, Indian Institute of Technology Roorkee, Roorkee-247667, Uttarakhand, India email:

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (2010) 182, 329 342 doi: 10.1111/j.1365-246X.2010.04609.x Spatial heterogeneity of tectonic stress and friction in the crust: new evidence from earthquake

More information

SOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION

SOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION SOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION Kimiyuki Asano 1 and Tomotaka Iwata 2 1 Assistant Professor, Disaster Prevention

More information

DATA REPOSITORY SUPPLEMENTARY MATERIAL. We analyzed focal mechanism solutions for aftershocks of the Loma Prieta

DATA REPOSITORY SUPPLEMENTARY MATERIAL. We analyzed focal mechanism solutions for aftershocks of the Loma Prieta GSA Data Repository item 2007141 DATA REPOSITORY SUPPLEMENTARY MATERIAL DR.1 Selection of aftershock sets and subsets We analyzed focal mechanism solutions for aftershocks of the Loma Prieta earthquake

More information

Short-Term Properties of Earthquake Catalogs and Models of Earthquake Source

Short-Term Properties of Earthquake Catalogs and Models of Earthquake Source Bulletin of the Seismological Society of America, Vol. 94, No. 4, pp. 1207 1228, August 2004 Short-Term Properties of Earthquake Catalogs and Models of Earthquake Source by Yan Y. Kagan Abstract I review

More information

Dynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement

Dynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement Emily Morton Geop 523: Theoretical Seismology 29 April 2011 Introduction Dynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement Earthquakes can be triggered by other earthquakes.

More information

Heterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault model

Heterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault model Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L21306, doi:10.1029/2008gl035614, 2008 Heterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault

More information

Resolving Stress Components and Earthquake Triggering

Resolving Stress Components and Earthquake Triggering Resolving Stress Components and Earthquake Triggering Earthquake Triggering Do certain events make an earthquake more likely to occur? Earthquakes Slow Slip Wastewater Fluids Dams The focus of this presentation

More information

RELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE

RELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE RELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE Paulino C. FEITIO* Supervisors: Nobuo HURUKAWA** MEE07165 Toshiaki YOKOI** ABSTRACT

More information

Probabilities for Jumping Fault Segment Stepovers

Probabilities for Jumping Fault Segment Stepovers GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:1.129/, Probabilities for Jumping Fault Segment Stepovers Bruce E. Shaw Lamont-Doherty Earth Observatory, Columbia University, New York James H. Dieterich

More information

Simulation of earthquake rupture process and strong ground motion

Simulation of earthquake rupture process and strong ground motion Simulation of earthquake rupture process and strong ground motion Takashi Miyatake (1) and Tomohiro Inoue (2) (1) Earthquake Research Institute, University of Tokyo, Yayoi, Bunkyo, Tokyo, 113-0032, Japan

More information

Seismic Activity near the Sunda and Andaman Trenches in the Sumatra Subduction Zone

Seismic Activity near the Sunda and Andaman Trenches in the Sumatra Subduction Zone IJMS 2017 vol. 4 (2): 49-54 International Journal of Multidisciplinary Studies (IJMS) Volume 4, Issue 2, 2017 DOI: http://doi.org/10.4038/ijms.v4i2.22 Seismic Activity near the Sunda and Andaman Trenches

More information

Mechanics of Earthquakes and Faulting

Mechanics of Earthquakes and Faulting Mechanics of Earthquakes and Faulting Lecture 18, 16 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Earthquake Magnitude and Moment Brune Stress Drop Seismic Spectra & Earthquake Scaling laws Scaling and

More information

Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics

Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Naoyuki Kato (1) and Tomowo Hirasawa (2) (1) Geological

More information

High-resolution image of Calaveras Fault seismicity

High-resolution image of Calaveras Fault seismicity JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. B9, 2186, doi:10.1029/2001jb000633, 2002 High-resolution image of Calaveras Fault seismicity David P. Schaff, 1 Götz H. R. Bokelmann, and Gregory C. Beroza

More information