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 that occurred between October of 1989 and December of 1990 (Figs. 1, 2, 5). Approximately 90% of these aftershocks have a magnitude less than 3, and almost all of those have magnitudes greater than 1.5. We visually subdivided the data into seventeen different sets on the basis of spatially distinct clusters of seismicity (distinguished by color in Figs. 1, 2, and 5) using the statistical analysis program JMP (SAS Institute Inc.), which includes features for rotating plots of three-dimensional data and for selecting out individual groups. We tested whether each data set records a homogeneous or a heterogeneous deformation, and if it was heterogeneous, we subdivided the set further into subsets each of which records a homogeneous deformation. We provide details of the protocol we used for assigning data to subsets in the following discussion. This process of subdivision resulted in a total of 33 separate subsets. Subsets within the sets (plotted with different symbols of the same color in Figs. 1, 2, and 5) were separated according to (1) the association of data with a prominent point maximum of P and/or T axes; and (2) the association of data with a prominent cluster of solutions, where the solutions are those for data sets chosen from all the data in a given subset by random selection with replacement. The protocol for assigning a datum to a particular solution is determined iteratively, as described below. (1) We make a Kamb-coutoured stereonet plot of the P and T axes for the focal mechanisms in each set, to reveal the preferred orientations of the axes (we used

GSA Data Repository item 2007141 Allmendinger's StereonetPPC v.6.0). The P and T axis for the data almost always cluster into one or more groups for each axis. If the P and T axes plot shows multiple maxima for P and/or for T, the data set is probably heterogeneous and reflects a mixture of different deformations. We define initial data subsets by selecting the data that are associated either with individual P and T axis point maxima, or at worst with P or T axis girdles. In some cases, this separation may be done conveniently using histograms of the trends and/or plunges of the P and/or T axes (we used the statistical program JMP (SAS Institute Inc.) for this purpose). The inversion solutions for the subsets are generally in the vicinity of, but not identical to, the P and T axis maxima. (2) Alternatively, or in addition, we define 100 to 200 data sets of the same size as the original set, selecting members of each new set from the original data set by random selection with replacement, and we invert each new set for a best-fit solution. To avoid false minima and find the best-fit solution, we use 5 random-restarts for each inversion, each restart starting from a set of randomly selected values for the initial model parameters. We then plot the resulting best-fit solutions for the principal instantaneous deformations d 1 and d 3 on a stereonet. If the solutions for these data sets concentrate around a single model, then that is taken as evidence that the data reflect a homogeneous deformation. If the solutions for d 1 and d 3 show multiple maxima, or the solutions are widely dispersed, then, either the data are inhomogeneous, or the solution is poorly resolved by the data. The point maxima for the P and T axes, and the point maxima for the d 1 and d 3 axis solutions to the randomly selected data sets provide initial trial solutions for the different deformations in the heterogeneous data set.

GSA Data Repository item 2007141 For this second procedure to identify the different solutions, an inhomogeneous data set should contain somewhat comparable numbers of events in each set, because the procedure relies on the random selection with replacement of the data to produce some data sets that will be best-fit by each of the possible solutions. The orientations of these solution-clusters on stereograms indicate the general location of the solutions for the subsets. These clusters usually are closely related to maxima in the P and T axis plots. We divide the data set into subsets, or refine the subsets defined in (1), such that each subset is most consistent with one of the solution clusters defined in (2), using the procedure described below. In a large data set that is overwhelmingly dominated by events associated with one solution but has a relatively small number of events that define a different solution, the probability is small of producing a data set by random selection with replacement that is dominated by the subordinate subset. Therefore, the solution to the subordinate subset is not likely to appear as a significant cluster in the inversion solutions for the randomly selected data sets. Similarly, if a data set contains more than two or three subsets with different solutions, the random selection of data might not reliably pick out the different solutions, because the probability of randomly selecting numerous data sets that are dominated by one of the homogeneous subsets may become too small. We have also found, especially for shear-plane/slickenline data gathered in the field, that orientation biases in the data set can lead to erroneous indications of solutions with this random selection method. For such data sets, criterion (1) provides a more reliable way of identifying subsets.

GSA Data Repository item 2007141 In refining the definition of the subsets, we determine the consistency of each individual datum with a particular subset solution according to the following procedure: i. Identify K initial trial solutions for K different subsets by locating the approximate average d 1 and d 3 pairs of orientations that represent each cluster of solutions either from the P and T axis maxima described in (1) or from the inversions of randomly selected data sets as described in (2); ii. Determine the misfit e (k) i of every datum (i = 1:N) in the whole data set to each of the (k = 1:K) trial solutions, and assign each datum to a subset α based on the solution for which the misfit for the datum is a minimum. Thus for a given subset α, if we renumber all the data in that subset j = 1:n α, then each of those data must satisfy the relation (α e ) (k ) j < e j j =1:n α k =1K]α[KK where the notation implies that k takes on any value from 1 to K except for α. The same relation holds for any specific subset numbered as α = 1:K. iii. Calculate the inversion solutions for each subset using only those data assigned to each subset. (k iv. Calculate the new misfits e ) i for each datum in the whole data set (i = 1:N), for each of the new subset solutions (k = 1:K) (α ) v. Calculate the means e j (α ) and the standard deviations SD e j of the misfits for just the data (j = 1:n α ) that are assigned to each subset (α = 1:K). vi. Assign the datum to a subset according to the following criteria: If the P and T axes for a datum k (belonging to the set 1:N) appear to be a part of a particular maximum α in the plot of all P and T axes for the set (which we define qualitatively as described in

GSA Data Repository item 2007141 (1) above), and if the misfit of that datum for the α th solution is within one standard deviation of the mean misfit for the α subset as determined in (iii), (α e ) (α ) k e i (α ) < SD e i i =1:n α where the angled brackets denote the mean and SD... denotes the standard deviation from the mean, then the datum is assigned to the α th subset, even if the misfit for a different solution β α is less than the misfit for the α th solution. If the P and T axes for a datum j are associated with a particular P or T axis maximum α, but its misfit for the α th solution exceeds the mean misfit for that subset by more than one standard deviation, (α e ) (α ) j e i (α > SD e ) i, i =1:n α and if the misfit of that datum for a different solution β α is within one standard deviation of the mean misfit for the β subset, i.e., (β e ) (β ) j e i (β < SD e ) i, i =1:n β then the datum is reassigned to the β subset; In all other cases, for example if a datum has no clear association with a particular P or T axis maximum, or if the misfits for the datum exceed the means for all solutions by more than one standard deviation, then it is assigned to the subset for which its misfit for the associated solution is a minimum. Steps (iii) through (vi) are iterated to a stable definition of the subsets. vii. Finally, the spatial clustering of events within individual hypocenter sets is in some cases a criterion used to test different possible subsets.

GSA Data Repository item 2007141 The identification of subsets defined by this protocol resulted in generally stable subsets each of which records a homogeneous deformation. Some data were consistent with more than one solution, however, and other data did not fit any of the solutions well, so we cannot rule out the possibility that this protocol results in an imperfect separation. DR.2 The FLTSLP Program We used our program FLTSLP (Twiss and Guenther, 2002; see "Acknowledgements" for availability) to search for the best-fit deformation model that accounts for the shear-plane/slip-line data given by the aftershock focal mechanism solutions. For a given shear-plane-normal, the orientation of the slip-line is uniquely determined from a given set of model parameters. For a given trial deformation model, the program uses a conjugate gradient technique (Press et al. 1989) to find the orientation of the model shear-plane-normal that, with the associated calculated slip-line, optimizes the misfit to each datum shear-plane/slip-line pair. The measure of misfit to each datum is the mean of the cosines of two angles, the first measured between the datum shearplane-normal and the model shear-plane-normal, and the second measured between the datum slip-line and model slip-line. For small angles, this is approximately the same as the cosine of the unique rotation angle that takes the model shear-plane/slip-line pair into the observed datum pair (Unruh et al. 1996). For each datum, therefore, neither the shear-plane nor the slip-line is held fixed in finding the best-fit model pair. The modelmisfit assigned to a given model is then the average of the misfits for all the data shearplane/slip-line pairs. The best-fit model is the set of five model parameters that minimizes the model-misfit and we determine this set of parameters by using a downhill simplex routine ( amoeba, Press et al. 1989). The program avoids false minima in the

GSA Data Repository item 2007141 five-dimensional model-misfit surface by restarting the search from at least 20 different randomly selected starting models. In the process of inverting seismic focal mechanisms, the inversion program always tests both possible nodal planes for the fit to the model, and it selects as the preferred nodal plane the one that yields a minimum misfit. Michael (1987) has done tests with his inversion algorithm, which does not include the relative instantaneous rotation and which uses different misfit criteria from ours, and has shown that the minimum misfit plane is very often not the correct choice. In such a case, the misfit will be smaller than the misfit to the correct shear plane. We have not run such tests with our program, and thus accept the resulting orientation maxima of the preferred nodal planes (see Section 3.2) as an indication, rather than a definitive determination, of the real dominant orientations of these planes. To determine the confidence limits on the model parameters, we invert 2000 bootstrap data sets, selecting each data set from the original data set by random selection with replacement. The inversion of each bootstrap data set uses the best-fit model as the starting model. To determine the 95% confidence limits we calculate the Mahalanobis distance of each bootstrap model from the mean bootstrap model. That distance is defined as the square root of the sum of the squares of the distance of each model parameter from same parameter of the mean model, where all parameters are normalized to values between 0 and 1. We ignore the 5% (100 models) with the largest distances, which leaves us with a cloud of models that defines the 95% confidence limits for each of the model parameters. We use the downhill simplex routine amoeba (Press et al. 1989) to find the mean bootstrap model by minimizing the L1 norm of the normalized

GSA Data Repository item 2007141 parameter distances (sum of absolute values of the distances) between the mean model and those bootstrap models within the 95% confidence limits. The calculation for the mean model and for the confidence limits is iterated to a stable set of bootstrap models that define the 95% confidence limits. DR.3 Fitting Planes to Hypocenters Alignments A plane in space is characterized by the equation n x L = 0 (DR3.1) where n is the unit normal to the plane, x is the vector to any point in the plane, the dot indicates the scalar product, and L is the plane location constant, which is the orthogonal distance to the plane from the coordinate origin. We determined the best-fit plane to a set of hypocenters by minimizing the L1 norm, which is a measure of the misfit defined by the sum of the absolute values of the perpendicular distances of all hypocenters from the plane. Compared with the more usual L2 norm, or least squares misfit, the L1 norm gives less weight to the outliers on the solution. The misfit for each hypocenter is thus the left side of Equation (DR3.1) where x is the vector to the hypocenter and L is the orthogonal distance from the origin to the best-fit plane. Two coefficients of determination, c 13 and c 23, define the fit of the plane to the data (1) m i (3) m i c 13 i i i m i (1) (2) m i (3) m i c 23 i i i m i (2) (DR3.2) where m i (α) is the perpendicular distance of the i th datum to the α th plane, and where α = 3 for the best-fit plane, α = 1 for the worst-fit plane normal to the best-fit plane through the mean point of the data, and α = 2 for the plane normal to the other two planes through the

GSA Data Repository item 2007141 mean point of the data. Thus the ranges of values for the two coefficients are 0 c 13 1; 0 c 23 c 13. The significance of different values of the coefficients is indicated in Table DR.3-1. The closer c 13 is to a value of 1, the closer the points are to lying in a plane, because then, in the first Equation (DR3.2), the length of the minimum misfit axis approaches zero, Σ m (3) i 0. The closer the ratio c 23 /c 13 is to 1, the closer the projection of hypocenters onto the best-fit plane approximates a circular distribution, because under those circumstances, the length of the intermediate misfit axis approaches the length of the maximum misfit axis Σ m i (2) Σ m i (1). The closer the ratio c 23 /c 13 is to 0, the closer the points are to defining a linear trend, because then length of the intermediate misfit axis approaches the length of the minimum misfit axis Σ m i (2) Σ m i (3). The local hypocenter alignment planes for the different sets were defined visually by rotating the data in a three-dimensional plot and noting orientations in which distinct planar alignments of the hypocenters occurred. The aligned hypocenters were then selected and the best-fit plane was calculated by minimizing the L1 norm of perpendicular distances from hypocenter to the plane. We used the three dimensional 'spinning plot' in the statistical program JMP (SAS Institute, Inc.) for this purpose. In a number of cases, the different alignment planes are defined predominantly by hypocenters from one particular subset, which are defined according to completely independent criteria. In these cases, we associate the different inversion solution for each subset with the planar structures defined predominantly by that subset. In other cases, there seems to be no preference of individual subsets to be associated with specific planar alignments, and in these cases, we infer that a partitioning into different deformations, as

GSA Data Repository item 2007141 defined by the different subsets, has occurred within the volume containing the set of hypocenters. We discuss this partitioning in 3.8 and 4.2. We report the strike and dip of these best-fit planes in Table DR2, where we assume the convention of the right-hand rule: Facing in the direction of the strike azimuth, the plane dips down to the right. We also report the coefficients of determination c 13 and c 23 and the ratio c 23 /c 13. DR.4. Kinematics of Strain Partitioning The triaxial strain is calculated by summing two partial plane-strains that have different orientations, each of which is the sum of strains due to slip on a conjugate set of shear planes. We express the shear on each shear plane in terms of local orthogonal φ coordinate axes ξ ) α for the two partial strains (f = 1, 2), and for the two conjugate shear φ planes associated with each partial strain (φ = 1, 2) (Fig. DR1). ξ ) 1 is normal to the φ shear plane, ξ ) φ 2 is parallel to the intersection line of the conjugate planes, and ξ ) 3 is parallel to the slip direction. Each partial plane-strain has its separate principal axes x ) k, which are the symmetry axes for the conjugate shear planes (Fig. DR1). The shearing on φ each shear plane defined in the ξ ) α coordinates is transformed into the x ) k coordinates so that the strain from two conjugate shear planes can be summed to give a partial strain. Each of the two partial strains is then transformed into the X K coordinate system so the partial strains can be summed to give the total strain. X K are the principal axes of the total triaxial strain (Fig. DR1). The details of this derivation follow. Refer to Figure DR1 for the relations among the different axes and angles used in the derivation. The displacement gradient associated with shear on a shear plane is most easily φ described in terms of coordinates defined such that ξ ) 1 is normal to the shear plane,

GSA Data Repository item 2007141 ξ 3 φ ) is parallel to the slip direction, and ξ 2 φ ) = ξ 3 φ ) ξ 1 φ ), where the right side of this equation is the vector product, or cross product. In such a coordinate system, the slip parallel to ξ 3 φ ) on a set of parallel shear planes produces a displacement gradient in the material given by 0 0 0 φ ) ν α φ ) ξ = 0 0 0 β S φ ) 0 0 (DR4.1) φ where ν ) α is the displacement field smoothed over the fault spacing. The strain from conjugate sets of shear planes (φ = 1, 2) with slip perpendicular to the intersection lines of the conjugate shears can be expressed in terms of the coordinates x ) k for each partial strain (f = 1, 2), defined so that x ) 1 bisects the acute angle between the normals to the two sets of shear planes (the obtuse angle between the planes themselves), x ) 2 is parallel to the intersection line of the planes, and x ) 3 bisects the obtuse angle between the normals to the conjugate sets of shear planes (i.e. the acute angle between the two planes themselves). For each of the partial strains (f = 1, 2), the components for the φ displacement gradient in the local coordinates ξ ) α for each of the conjugate shear planes ) (φ = 1, 2) given by Equation (DR4.1) are transformed into components in the x k coordinates and added together. The transformation is given by u k φ ) x l ) = ν φ ) α ξ Q φ ) φ ) φ ) αk Q βl β (DR4.2) where u k φ ) are the displacement components in the x k ) coordinates of the displacement on the φ th φ set of shear planes, and Q ) αk is the orthogonal transformation that relates

GSA Data Repository item 2007141 vector components in the ξ β φ ) coordinates to components of the same vector in the x k ) coordinates. If λ β φ ) and i k ) are the unit base vectors respectively in the ξ β φ ) and x k ) coordinates, then, φ Q ) φ αk λ ) ) α i k (DR4.3) where the right hand side is the scalar product. We define θ φ ) to be the angle between the normal to the (fφ) th φ shear plane (i.e. ξ ) 1 ) and x ) 1, where for φ = 1, the angle is a positive (right-handed) rotation about the positive x ) 2 axis, and for φ = 2 the angle is a negative rotation. The orthogonal transformations in Equation (DR4.3) then become, Q 1) αk = cosθ 1) 0 sinθ 1) 0 1 0 sinθ 1) 0 cosθ 1) Q 2) αk = cosθ 2) 0 sinθ 2) 0 1 0 (DR4.4) sinθ 2) 0 cosθ 2) The total displacement gradient for the partial strain is the sum of the expressions like Equation (DR4.2) for both sets of shear planes, u k ) N φ ) u k x = ) (DR4.5) ) l x l φ=1 We assume for simplicity that θ 1) = θ 2) = θ for f =1, 2 (DR4.6) which defines conjugate shear planes that are symmetrical about the x ) 1 and x ) 3 axes, with a comparable geometry for both partial strains. In this case, performing the transformations indicated in Equation (DR4.2) with Equations (DR4.4), and adding as prescribed by Equation (DR4.5) gives

GSA Data Repository item 2007141 Σ ) sinθ cosθ 0 Δ ) sin 2 θ ) u k ) x = 0 0 0 l Δ ) cos 2 θ 0 Σ ) sinθ cosθ (DR4.7) where Σ ) S 1) + S 2) Δ ) S 1) S 2) (DR4.8) This is a plane-strain with no deformation parallel to the x ) 2 axis, as indicated by the zeros in the second row and column. If the magnitude of the shear on both shear planes for both partial strains is the same, S 1) = S 2) = S whereby Σ ) = 2S, Δ ) = 0 for f =1, 2 (DR4.9) Thus the off-diagonal terms in Equation (DR4.7) become zero, the deformation is irrotational, and the x ) k axes are the principal axes of the strain. The simplifying assumptions in Equations (DR4.6) and (DR4.9) maximize the symmetry of the system and ensure that the coordinate axes x ) k are the principal axes of the associated strains, and that the strains are irrotational. Our data indicate that a triaxial strain is accommodated by a pair of plane strains (f = 1, 2), and that each plane strain is commonly accommodated on a conjugate set of shear planes (Table 2; Fig. 14). Thus each plane strain (i.e. for f = 1, 2) is defined by Equation (DR4.7), and the principal axes x ) k of each of the plane strains are related by a rotation about the common axis. In our data, that common axis is usually the d 3 axis (Fig. 14A, B; the X 3 axis in Fig. DR1A), which means that d (1) 3 is parallel to d (2) 3 ( x (1) 3 is parallel to x (2) 3 in Fig. DR1A), but it can also be the d 1 axis (Fig. 14C, D; X 1, x (1) 1, and x (2) 1 are all

GSA Data Repository item 2007141 parallel in Fig. DR1B). We can calculate the triaxial deformation by adopting a coordinate system X K defined such that the two sets of conjugate shear planes (four shear planes in all) are symmetrically arranged with respect to these axes. For the case in which d 3 is the common axis for the two sets of conjugate shear planes (Fig. 14A, B), we take X 3 to be parallel to d 3, and X 1 to bisect the angle between the maximum extension axes for the two plane-strains, i.e. d (1) 1 and d (2) 1 ( x (1) 1 and x (2) 1 in Fig. DR1A). Thus because each set of conjugate shear planes has its own local coordinate axes x ) k, we must transform the components of deformation given by Equation (DR4.7) in the each of the x ) k coordinates into the X K coordinate system, in order to add the two deformations. The components of the total displacement gradient tensor U K X L in the X K coordinates can therefore be calculated from P ) U K u = m ) X L x T ) mk T ) nl (DR4.10) n f =1 where we take P = 2, where u m ) x n ) is given by Equation (DR4.7) for each f = 1, 2, and where T ) mk is the orthogonal transformation that transforms components from the x k ) coordinates to the X K coordinates. If i m ) and I K are the unit base vectors in the x k ) and the X K coordinate systems, respectively, then the components of the two transformation tensors T ) mk are defined by the scalar product, T ) mk i ) m I K (DR4.11) We define α (f) to be the angle between X 1 and x 1 ), such that is a positive (right-hand) rotation by an angle α (1) about the X 3 axis carries X 1 to x 1 (1) (Fig. DR1A), and a negative

GSA Data Repository item 2007141 rotation by an angle α (2) about X 3 carries X 1 to x 1 (2). The transformation matrices for f = 1 and 2, then, are T (1) mk = cosα (1) sinα (1) 0 sinα (1) cosα (1) 0 0 0 1 T (2) mk = cosα (2) sinα (2) 0 sinα (2) cosα (2) 0 (DR4.12) 0 0 1 We make the simplifying assumption that the angle between X 1 and x 1 ) for both values of f is the same, which means the conjugate planes are symmetrically arrayed about X 1 : α (1) = α (2) = α (DR4.13) Using the assumptions in Equations (DR4.6), (DR4.9), and (DR4.13) as well as Equations (DR4.7) and (DR4.12) in Equation (DR4.10) then gives for the triaxial displacement gradient tensor, U K X L cos 2 α 0 0 = 2Σsinθ cosθ 0 sin 2 α 0 0 0 1 (DR4.14) Our simplifying assumptions in Equations (DR4.6), (DR4.9), and (DR4.13) have imposed symmetries on our model system that results in the triaxial instantaneous strain being irrotational and the X K axes being the principal axes of that strain. For the case in which d 1 is the common axis for the two partial plane-strains (Fig. 14C, D; Table 2), the derivation proceeds in a similar manner. Equations (DR4.10) and (DR4.11) remain the same, but we define the angle β (f) to be the angle between X 3 and x ) 3, where a positive (right-hand) rotation of angle β (1) (1) about the X 1 axis carries X 3 to x 3 (Fig. DR1B), and a negative rotation of angle β (2) about the X 1 axis carries X 3 to x (2) 3. Then for f = 1 and 2, the transformation matrices are

GSA Data Repository item 2007141 T (1) mk = 1 0 0 0 cosβ (1) sinβ (1) 0 sinβ (1) cosβ (1) T (2) mk = 1 0 0 0 cosβ (2) sinβ (2) (DR4.15) 0 sinβ (2) cosβ (2) We make the same simplifying assumptions as in Equations (DR4.6) and (DR4.9), and we assume that the angle between X 3 and x 3 ) for both values of f is the same (Fig. DR1B), β (1) = β (2) = β (DR4.16) which means the two sets of conjugate planes are symmetrically arrayed about X 3. With these assumptions, we find the triaxial deformation gradient from Equations (DR4.7) and (DR4.10): U K X L 1 0 0 = 2Σsinθ cosθ 0 sin 2 β 0 0 0 cos 2 β (DR4.17) Equations (DR4.14) and (DR4.17) with Equation (DR4.9) give the Equations (6) and (7) presented in 4.2 of the main text.

GSA Data Repository item 2007141 DATA REPOSITORY MATERIAL: REFERENCES CITED Michael, A.J., 1987, Use of focal mechanisms to determine stress: A control study: Journal of Geophysical Research, v. 92, p. 357-368. Press, W.H., Flannery, B.P., Teukolsky, S.A., and Vettering, W.T., 1989, Numerical Recipes, the Art of Scientific Computing (FORTRAN Version): New York, Cambridge University Press, 702 p. Unruh, J.R., Twiss, R.J. and Hauksson, E., 1996, Seismogenic deformation field in the Mojave block and implications for the tectonics of the eastern California shear zone: Journal Geophysical Research, v. 101(B4), p. 8335-8362.

TABLE DR1. INVERSION RESULTS FOR LOMA PRIETA SUBSETS Data Set No. of data d 1 trend d 1 plung e d 1 C.I. angle from mean d 2 trend d 2 plung e d 2 C.I. angle from mean d 3 trend d 3 plung e d 3 C.I. angle from mean D D C.I. from mean 1 S.D. Max for 95% C.I. W # W C.I. from mean Best-fit Mean 1 S.D. Max for 95% C.I. Misfit Bestfit Bestfit Mean Mean 1 S.D. Max. for 95% C.I. Bestfit Bestfit Mean Mean 1 S.D. Max for 95% C.I. Bestfit Bestfit Mean Mean 1 S.D. Max for 95% C.I. Bestfit Mean Bestfit Mean V Vertical deformatio n ratio Best-fit 1-32 29 073.86 1.04 344.83-43.11 162.75-46.87 0.54 0.62@ 4.6110-0.54 073.39 1.38 5.7 344.71-43.58 6.3 161.94-46.39 4.6 0.54 ±0.04 0.58 ±0.16 3.93 17.1 18.5 12.5 +0.12 / -0.13 +0.42 / -0.60 1-42 16 204.73 74.04 254.75-10.41 342.52 11.96 0.40 0.38&@ 2.8280 0.88 205.40 74.34 3.3 255.12-10.27 3.9 342.97 11.69 3.3 0.41 ±0.04 0.40 ±0.08 2.40 12.0 17.2 14.7 +0.15 / -0.12 +0.30 / -0.30 2-32 160 087.22 20.55 54.53-65.99 172.70-11.88 0.66* - 0.57&@ 087.12 20.54 4.2 54.37-65.99 4.2 172.60-11.89 3.3 0.66 ±0.04-0.55 ±0.08 13.2 14.8 10.1 +0.12 / -0.12 +0.34 / -0.22 4.5686 0.27 4.44 2-42 29 105.69-38.93 111.79 50.91 018.11 2.98 0.36 0.08 3.6093 0.29 105.73-38.93 3.7 111.81 50.92 5.7 018.14 2.98 5.1 0.37 ±0.04 0.09 ±0.09 3.27 16.5 24.8 23.9 +0.18 / -0.13 +0.29 / -0.32 3-1 20 310.52 7.83 026.83-59.84 044.88 28.92 0.70* 0.29&@ 3.2289-0.06 310.43 7.56 3.7 027.02-60.24 4.5 044.58 28.59 4.3 0.69 ±0.03 0.30 ±0.11 2.79 14.0 21.1 22.9 +0.10 / -0.12 +0.46 / -0.22 3-2 18 336.45 72.88 124.85 14.70 217.12 8.58 0.41-0.84& 1.9832 0.89 336.78 73.63 2.5 124.70 13.98 3.4 216.79 8.35 3.1 0.41 ±0.03-0.83 ±0.04 1.72 8.0 10.9 10.2 +0.13 / -0.09 +0.14 / -0.12 1

4-52 114 091.50 23.76 307.81 61.35 008.30-15.06 0.49 0.21 5.4446 0.08 91.33 24.28 4.0 307.69 60.75 4.8 008.42-15.31 4.0 0.49 ±0.03 0.24 ±0.08 5.23 12.2 17.8 14.7 +0.09 / -0.10 +0.32 / -0.31 4-62 31 085.38-88.58 248.01-1.35 158.02 0.42 0.80* 0.07@ 5.1099 1.00 80.54-86.62 6.4 248.87-3.31 6.9 158.90 0.68 3.7 0.77 ±0.05 0.07 ±0.10 4.66 26.0 27.2 14.5 +0.14 / -0.23 +0.27 / -0.41 5-52 18 278.74-21.61 130.21-65.08 013.49-11.80 0.54 0.07@ 5.1053 0.13 278.62-20.99 5.2 130.50-65.69 7.3 013.19-11.71 6.0 0.55 ±0.08 0.08 ±0.09 4.35 23.6 26.1 22.3 +0.27 / -0.34 +0.34 / -0.28 5-62 14 086.69-38.42 297.08-47.41 189.45-15.56 0.17 0.53&@ 6.8328 0.14 086.50-37.74 5.7 297.84-47.82 11.2 189.35-16.03 11.4 0.21 ±0.09 0.52 ±0.12 5.44 26.7 44.1 42.9 +0.35 / -0.15 +0.42 / -0.38 6-0 52 065.32-86.43 106.96 2.67 016.85 2.37 0.37 0.52 5.9876 0.99 051.30-87.81 11.4 108.01 1.20 13.9 017.97 1.83 12.3 0.38 ±0.05 0.51 ±0.15 5.38 36.1 52.6 46.2 +0.19 / -0.16 +0.37 / -0.55 7-02 50 003.72-40.59 110.94-19.06 039.93 43.28 0.41 0.65&@ 4.6365-0.01 003.79-40.93 3.5 111.46-19.28 4.5 040.36 42.80 4.2 0.41 ±0.03 0.65 ±0.08 4.30 12.4 16.4 14.9 +0.12 / -0.13 +0.24 / -0.28 8-51 65 179.03 74.49 282.53 3.71 013.53 15.04 0.69* -0.56& 3.9120 0.84 179.03 74.49 2.1 282.53 3.71 2.3 013.53 15.04 2. 0.69 ±0.03-0.56 ±0.04 3.78 6.7 8.3 8.3 +0.09 / -0.1 +0.15 / -0.15 8-61 40 289.53 12.45 143.57 75.08 201.33-8.09 0.44 0.27& 3.8773-0.04 289.53 12.45 2.3 143.57 75.08 3.3 201.33-8.09 2.8 0.44 ±0.03 0.27 ±0.06 3.70 8.5 11.1 10.3 +0.11 / -0.08 +0.23 / -0.21 9-01 54 297.26-45.13 111.15-44.71 024.18 3.06 0.53-0.71@ 7.7180 0.52 296.45-46.00 13.7 110.84-43.86 17.3 023.54 2.80 12.8 0.53 ±0.06-0.67 ±0.19 6.76 53 62.0 43.5 +0.19 / -0.21 +0.82 / -0.33 2

10-52 29 116.12 37.82 260.61 46.36 011.1 18.46 0.52 0.1 4.4088 0.29 115.84 37.80 4.5 261.97 46.95 7.0 011.69 17.50 6.1 0.53 ±0.05 0.11 ±0.07 3.95 17.1 21.3 20.1 +0.16 / - 0.18 +0.22 / -0.23 10-62 17 267.89 78.14 108.15 11.14 197.36-4.01 0.37-0.35 3.8335 0.95 267.86 78.14 3.2 108.15 11.14 7.4 197.36-4.01 7.1 0.38 ±0.05-0.35 ±0.09 3.26 12.2 35.0 34.8 +0.22 / -0.15 +0.42 / - 0.33 11-33 30 198.50-33.95 132.93 31.57 074.03-40.05 0.55-0.69& 3.3017-0.11 197.31-33.29 3.8 132.28 32.74 4.7 074.41-39.59 4.0 0.54 ±0.03-0.68 ±0.07 3.05 14.2 16.7 14.5 +0.07 / -0.08 +0.19 / -0.21 11-43 27 262.60 8.81 265.05-81.18 352.66 0.37 0.46-0.73&@ 263.68 7.53 4.8 270.50-82.41 6.8 353.80 0.89 6.3 0.48 ±0.05-0.73 ±0.07 15.2 19.1 18.8 +0.20 / -0.20 +0.25 / -0.24 5.9291-0.03 5.30 12-5 17 012.90-34.37 242.42-43.50 123.38-27.09 0.48 0.73&@ 3.0603 0.1 013.00-34.53 3.9 242.47-43.37 6.0 123.59-27.07 5.5 0.48 ±0.04 0.73 ±0.06 2.61 16.3 26.9 26.2 +0.15 / -0.12 +0.13 / -0.28 12-62 18 251.36-67.57 326.81 5.92 234.46 21.55 0.72-0.26 4.3781 0.68 252.39-68.03 5.2 327.78 5.81 6.1 235.53 21.11 6.0 0.69 ±0.06-0.26 ±0.12 3.63 22.5 24.1 22.4 +0.19 / -0.25 +0.39 / -0.31 12-72 17 050.06 40.76 250.77 47.34 329.22-10.46 0.58-0.32 4.3204 0.45 050.15 40.70 5.1 250.13 47.54 5.5 328.94-10.07 4.3 0.59 ±0.04-0.32 ±0.10 3.73 31.9 36.9 21.8 +0.18 / -0.19 +0.47 / -0.32 13-33 34 248.41 57.20 324.52-8.80 049.12 31.31 0.52-0.41&@ 247.96 56.81 3.8 324.27-8.80 3.9 048.78 31.71 3.6 0.53 ±0.03-0.44 ±0.09 10.7 11.3 9.9 +0.12 / -0.10 +0.25 / -0.26 4.3586 0.43 3.97 13-43 31 281.32-14.93 186.29-18.20 048.35-66.12 0.54 0.46&@ 5.2003-0.81 282.04-14.54 3.7 187.23-17.91 3.8 048.88-66.61 4.0 0.55 ±0.05 0.46 ±0.11 4.64 11.4 13.2 13.7 +0.19 / -0.16 +0.36 / -0.30 3

13-53 40 126.20-30.28 195.42 31.28 070.03 43.63 0.53 0.48& 3.0900-0.23 126.20-30.28 2.6 195.42 31.28 3.3 070.03 43.63 2.6 0.53 ±0.03 0.48 ±0.04 2.91 9.3 10.9 9.6 +0.08 / -0.11 +0.11 / -0.15 14 13 275.58-68.79 283.19 21.04 192.20 2.56 0.52 0.77& 5.5106 0.87 275.70-69.23 9.0 282.94 20.61 13.4 192.04 2.40 11.1 0.52 ±0.11 0.76 ±0.09 4.08 36.9 57.6 45.7 +0.22 / -0.38 +0.23 / -0.33 15-732 26 322.98-62.00 115.68-25.29 31.06 11.23 0.50-0.26@ 3.1399 0.74 322.83-61.96 3.1 115.77-25.37 3.6 31.13 11.14 2.8 0.51 ±0.03-0.24 ±0.08 2.81 9.8 13.8 10.3 +0.09 / - 0.10 +0.27 / -0.24 15-832 22 109.27 13.81 282.61 76.10 018.89 1.55 0.36* 0.23& 3.2971-0.1 109.27 13.81 2.5 280.79 76.04 3.8 018.79 1.98 3.3 0.38 ±0.03 0.23 ±0.05 3.00 8.6 11.5 11.1 +0.11 / - 0.09 +0.18 / -0.18 15-842 19 080.90 21.84 346.43 11.00 051.44-65.28 0.50-0.82& 2.3914-0.69 080.92 21.91 4.2 346.44 10.98 5.1 051.57-65.23 3.6 0.50 ±0.03-0.82 ±0.04 2.04 12.4 14.1 10.7 +0.09 / -0.08 +0.16 / -0.12 15-932 20 107.87-60.37 320.27-25.65 223.49-13.82 0.61-0.37&@ 108.14-60.10 3.5 320.61-25.88 3.7 223.70-13.93 3.0 0.62 ±0.06-0.35 ±0.07 13.3 12.9 10.7 +0.18 / -0.20 +0.24 / -0.22 16-11 27 159.40 53.98 121.71-29.91 222.59-18.16 0.38-0.66&@ 158.31 53.70 4.3 119.81-29.90 5.6 220.99-18.63 4.7 0.40 ±0.06-0.66 ±0.06 18.5 23.1 19.2 +0.25 / -0.14 +0.23 / -0.21 2.8533 0.72 2.46 4.0318 0.53 3.60 16-21 14 292.57 29.94 090.32 58.11 196.75 9.98 0.53 0.09 2.0074 0.25 293.37 29.94 3.2 090.04 57.90 4.0 197.24 10.50 3.8 0.55 ±0.04 0.10 ±0.09 1.60 12.7 13.9 11.7 +0.14 / -0.15 +0.35 / -0.21 17 36 125.42 18.41 354.66 62.99 042.02-19.05 0.51 0.54& 4.6033 0.01 124.69 18.79 5.3 352.36 63.20 5.8 041.18-18.39 5.1 0.51 ±0.04 0.54 ±0.14 4.16 27.5 36.1 28.2 +0.16 / -0.14 +0.46 / -0.49 4

Notes: Solutions are listed as both the best-fit value to the data and the mean of the bootstrap solutions that lie within the 95% confidence limits. The two solutions usually differ by less than 1 in the orientations of the principal instantaneous strain axes, and less than 0.05 in D and W. Confidence Interval: listed first in terms of 1 S.D. (standard deviation) from the mean of the bootstrap solutions, and second in terms of the maximum angle from the mean of those bootstrap solutions defining the 95% confidence limits. An asterisk in this column indicates a value of D that is significantly different from 0.5 at the 95% confidence level. # An ampersand '&' in this column indicates a value of W that is significantly different from 0 at the 95% confidence level. An 'at' sign '@' in this column indicates that the negative of this value and orientation of W is plotted in Figures 11B and 8C The misfit is the angle given by the inverse cosine of the mean of the misfit cosines for each datum. The misfit cosine for each datum is the mean of the cosines of the angles between the model and the datum shear-plane-normals, and between the model and the datum slip-directions. 5

TABLE DR2: PLANAR BEST-FITS (L1 NORM) TO HYPOCENTER ALIGNMENTS* Plane Association Best-Fit Plane Unit Normal n to Best-Fit Plane All data but Sets 1 and 3 Sets 2, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17 Plane Location Constant Strike Dip n 1 n 2 n 3 L [km] Coefficients of Determination Ratio Standard deviation [km] c 13 c 23 c 23 /c 13 strike of fault-zone segment Angle Between Strike of Best-Fit Plane and azimuth of plate motion 131.1 64.8 0.5952 0.682-0.425 4.24 0.9208 0.7436 0.81 1.776-17.3 Comment Southern Central Northern 2, 5, 7 142.5 76.2 0.7707 0.5906-0.2389 5.65 0.9007 0.616 0.68 1.316-5.9 4, 9, 11 125 62 0.5062 0.7237-0.469 4.8 0.9459 0.7143 0.76 0.291-23.4 12.62, 12.72, 13, 14, 15, 16, 17 136 65.1 0.652 0.6304-0.4213 2.79 0.8691 0.8261 0.95 0.356-12.4 1-ne Shallow Zone 295.5 88.1 0.4298 0.9023 0.0323 0.94 0.9759 0.8904 0.91 0.012-27.0-32.9 Mostly Subset 1.32 data; 3 data from Subset. 1.42; vertical. plane on NE of set 1-sw1 Shallow Zone 115.3 75.5-0.4135-0.8754 0.2504-0.71 0.9273 0.817 0.88 0.105-27.2-33.1 Subsets 1.32 & 1.42 data not in 1_ne; plane with steep dip SW 1-sw2 Shallow Zone 117.8 58.2 0.3958 0.752-0.5271 1.18 0.8732 0.9025 1.03 0.062-24.7-30.6 Subsets 1.32 & 1.42 data not in 1_ne; plane with moderate dip SW 2-ne Southern 141.9 88.5 0.7868 0.6166-0.0269 5.58 0.9319 0.7058 0.76 0.437-0.6-6.5 Subsets 2.32 & 2.42 events in the NE planar alignment 2-sw Southern 142.7 79.3 0.7818 0.5953-0.1856 5.07 0.9124 0.8656 0.95 0.165 0.2-5.7 All events in SW alignment plane 3-nw Off-plane set 297.9 44.6 0.3291 0.6206 0.7117-9.17 0.8038 0.5951 0.74 2.799-30.5 All data from the NW scattered events not in the SE cluster 3-se Off-plane set 326.8 72.5-0.7983-0.5217-0.3009 7.51 0.8257 0.637 0.77 0.446-1.6 Events in the SE cluster 4-62-ne Central 289.4 90 0.3314 0.9435 0.0004-4.77 0.9846 0.9616 0.98 0.005-15.6-39.0 Subset 4.62 data along NE plane 4-62-sw Central 291.7 89.9 0.3696 0.9292 0.002-6.11 0.9393 0.6875 0.73 0.015-13.3-36.7 Subset 4.62 data for short SW planar alignment 5-pl1 (deep) Southern 157.5 27.6 0.4279 0.1769-0.8863 10.99 0.9377 0.6259 0.67 0.046 15.0 9.1 Events below 9 km except for Datum No. 433 5-pl2 (shallow) Southern 142.2 76.8 0.7696 0.5965-0.2278 5.74 0.8506 0.7978 0.94 0.23-0.3-6.2 Events above 9 km except for Data Nos. 72, 433, 874 6-pl1 Shallow Zone 101.8 1.6 0.0058 0.0277-0.9996 4.86 0.9784 0.8881 0.91 0.004-40.7-46.6 Horizontal plane 6-pl2 Shallow Zone 290.3 28.6-0.1662-0.4489-0.878 3.38 0.9423 0.8552 0.91 0.021-32.2-38.1 Shallowly-NE-dipping plane 6-pl3 Shallow Zone 111 82.6 0.356 0.9255-0.1296 1.44 0.8735 0.9209 1.05 0.058-31.5-37.4 Steeply SW-dipping plane; N side of data set 6-pl4 Shallow Zone 111.8 80.8 0.366 0.9166-0.1606 0.13 0.8943 0.8051 0.90 0.02-30.7-36.6 Steeply SW-dipping plane; S side of data set 1

7-02 Southern 146.4 71 0.787 0.5237-0.3261 6.39 0.8072 0.6185 0.77 0.099 3.9-2.0 All subset data, refined data set 7-02-1 Southern 336.5 25.7 0.3972 0.1726 0.9014-4.67 0.958 0.7819 0.82 0.011 14.0 8.1 top subhorizontal alignment in [315, -01] view 7-02-2 Southern 322.3 26.1 0.3476 0.2686 0.8983-5.48 0.9644 0.8762 0.91 0.007-0.2-6.1 2nd from top subhorizontal alignment in [315, -01] view 7-02-3 Southern 308 17.3 0.1825 0.2338 0.955-7.01 0.9562 0.8655 0.91 0.005-14.5-20.4 3rd from top subhorizontal alignment in [315, -01] view lp7-02-4 Southern 329.1 11.8 0.175 0.1049 0.979-7.92 0.9518 0.772 0.81 0.011 6.6 0.7 4th from top subhorizontal alignment in [315, -01] view lp7-02-5 Southern 135.2 3.4-0.0426-0.0423 0.9982-9.52 0.9691 0.8717 0.90 0.002-7.3-13.2 5th from top subhorizontal alignment in [315, -01] view 8-51-pl1 Shallow Zone 287.9 55.4 0.2529 0.7834 0.5678 2.92 0.9681 0.8057 0.83 0.02-17.0-40.5 Shallow data defining NE-dipping plane 8-51-pl2 Shallow Zone 106.7 86.4 0.2875 0.9558-0.062 5.48 0.9436 0.8592 0.91 0.012-18.2-41.7 Data defining vertical SW boundary to Subset 8.51 8-61-pl3 Shallow Zone 109.9 57.7 0.2874 0.7953-0.5337 5.92 0.9299 0.853 0.92 0.032-15.0. -38.5 Data defining mod. SW dipping plane in Subset 8.61 9-01 Central 124.9 64 0.5142 0.7369-0.4388 4.1 0.8043 0.7135 0.89 0.07-0.1-23.5 10-52-pl1 Shallow Zone 315.4 83.6 0.7075 0.6978 0.112 0.93 0.8894 0.8355 0.94 0.064 0.3-13.0 Data selected to define vertical alignment plane 10-pl2 Shallow Zone 134.2 44 0.4837 0.4981-0.7196 3.39 0.9087 0.806 0.89 0.101-0.9-14.2 Data selected to define SWdipping plane at top of hypocenter distribution 10-pl3 Off-plane set 151.6 82.3 0.872 0.4708-0.1341-4.19 0.7569 0.6088 0.80 0.288 3.2 Data from NW cluster; steep offplane group 11-33 Central 117.3 79.7 0.4506 0.8747-0.1785 1.49 0.7332 0.3747 0.51 0.179-7.7-31.1 11-43 Central 128.7 60.3 0.5429 0.6784-0.4951 5.08 0.8649 0.6303 0.73 0.243 3.7-19.7 12-50 Off-plane set 153.1 81.8 0.8829 0.4472-0.1433-0.55 0.8481 0.6222 0.73 0.14 4.7 12-pl2 12-pl1 Northern Northern 13 Northern 13-14-pl1 Northern 13-14-pl2 Northern 13-pl3 Northern 13-pl4 Northern 13-pl5 Northern 14 Northern 142.3 88.9 0.7912 0.6113-0.0191-3.16 0.8408 0.4458 0.53 0.099 6.3-6.1 Data below 12 km along planar alignment from 12.62 and 12.72 139.4 63.8 0.6813 0.5843-0.4409 2.17 0.882 0.8282 0.94 0.061 3.4-9.0 Data above 11 km along planar alignment from 12.62 and 12.72 126.3 49.8 0.4523 0.6153-0.6457 8.33 0.722 0.5465 0.76 0.263-9.7-22.1 169.6 81.3 0.9723 0.1779-0.1517-7.77 0.8791 0.5299 0.60 0.073 33.6 21.2 Sets 13 (& 14) NS alignment on SE of the cluster 357.7 88.4 0.9988 0.0409 0.0282-13.32 0.8892 0.5989 0.67 0.066 41.7 29.3 Sets 13 (& 14) NS alignment on the SW of the cluster 168.7 85.4 0.9775 0.1952-0.0798-9.44 0.7273 0.5421 0.75 0.049 32.7 20.3 Set 13 NS alignment on the NE of the cluster 340.8 88.8 0.944 0.3292 0.021-10.79 0.8313 0.7951 0.96 0.023 24.8 12.4 Set 13 NS alignment on the NW of the cluster 77.6 83.6-0.2133 0.9705-0.1121 13.14 0.7011 0.5063 0.72 0.073-58.4-70.8 EW alignment on the N side of the cluster 129.6 54.1 0.5171 0.6241-0.5857 7.49 0.6936 0.2102 0.30 0.22-6.4-18.8 Set 14.0 data minus Data Nos. 103, 168, 388 2

15-pl1 Northern 183.1 8.9-0.1549 0.0083 0.9879-7.61 0.9127 0.8187 0.90 0.011 47.1 34.7 15-pl2 Northern 317.4 17.4 0.2198 0.2024 0.9543-9.16 0.962 0.911 0.95 0.006 1.4-11.0 15-pl3 Northern 314.1 34 0.3894 0.4015 0.829-5.78 0.9105 0.723 0.79 0.012-1.9-14.3 15-pl4 Northern 151.8 68.2-0.8181-0.4396 0.3709 2.5 0.891 0.7825 0.88 0.02 15.8 3.4 15-pl5 Northern 327.4 89.9 0.8423 0.5389 0.001-3.46 0.892 0.7044 0.79 0.035 11.4-1.0 16-pl1b Northern 131.2 14.2 0.1619 0.1853-0.9693 11.3 0.9613 0.8504 0.88 0.02-4.8-17.2 16-pl2b Northern 129.3 67.9 0.5862 0.7173-0.3766 4.03 0.9629 0.8516 0.88 0.018-6.7-19.1 17-pl1 Northern 226.1 14.3 0.1711-0.1779-0.9691 6.3 0.9722 0.9393 0.97 0.003-89.9 77.7 Sub-horizontal plane at top of data 17-pl2 Northern 179.1 66.8 0.9191 0.0137-0.3938-11.51 0.9422 0.9268 0.98 0.008 43.1 30.7 Moderate W dip 17-pl3 Northern 106.7 36.6 0.1716 0.5716-0.8024 16.68 0.8243 0.7219 0.88 0.037-29.3-41.7 Moderate S dip * Refer to Figure 8; see discussion in Section DR.4 Azimuth of plate motion vector is 148.4 3

1 TABLE DR3: SLIP LINE ORIENTATIONS ON MAIN FAULT SEGMENTS FOR SUBSET SOLUTIONS Shear_Plane Slip_Line Slip_Line Dip Trend Plunge Subset Shear_Plane Strike Slip_Line Rake (from strike azimuth) Normalized Slip Magnitude Southern 1-32 142.5 76.2 131.1 38.8 40.2 0.68 1-42 142.5 76.2 132.7 34.7 35.9 0.16 2-32 142.5 76.2 143.5 3.9 4.0 0.16 2-42 142.5 76.2 138.5 15.8 16.3 0.27 5-52 142.5 76.2 148.8 24.1 24.9 0.53 5-62 142.5 76.2 127.7 46.1 47.9 0.64 7-02 142.5 76.2 133.3 33.2 34.3 0.10 mean rake for planes 13.9 Plate slip vector (a) 142.5 76.2 148.4 22.3 23 Plate slip vector (b) 142.5 76.2 142.5 0.0 0.0 Central 6 125 62 193. 60.2 79.3 0.63 8-51 125 62 160.4 47.4 56.5 0.55 8-61 125 62 182.5 57.8 73.4 0.15 4-52 125 62 141.2 27.7 31.8 0.36 4-62 125 62 139. 24.5 28.0 0.35 9-01 125 62 253.4 55.8 110.4 0.07 11-33 125 62 50.1 61.2 82.8 0.14 11-43 125 62 130.3 9.9 11.2 0.11 mean rake for planes 45.4 Plate slip vector (a) 125 62 148.4 37 43. Plate slip vector (b) 125 62 142.5 29.6 34 Northern 10-52 136 65.1 142.8 14.2 15.7 0.39 10-62 136 65.1 156.3 36.8 41.3 0.16 12-5 136 65.1 65.6 63.8 81.5 0.24 12-62 136 65.1 233.3 64.9 93.1 0.62 12-72 136 65.1 153. 32.3 36.1 0.30 13-33 136 65.1 279.8 51.8 119.9 0.11 13-43 136 65.1 150.2 27.9 31.0 0.25 13-53 136 65.1 318.7 5.8 173.6 0.36 14 136 65.1 153.6 33.1 37.0 0.11 15-732 136 65.1 159. 40.1 45.2 0.18 15-832 136 65.1 143.2 15.1 16.7 0.32 15-842 136 65.1 171.1 51.1 59.1 0.41 15-932 136 65.1 255.2 62 103.2 0.34 16-11 136 65.1 347.5 48.4 124.5 0.20 16-21 136 65.1 157.9 38.8 43.7 0.26 17-02 136 65.1 245.9 63.7 98.7 0.19 mean rake for planes 62.2 Plate slip vector (a) 136 65.1 148.4 25.2 28. Plate slip vector (b) 136 65.1 142.5 13.6 15. 1

2 TABLE DR.3-1: INTERPRETATION OF COEFFICIENTS OF DETERMINATION FOR FITS OF PLANES TO HYPOCENTER DATA c 13 c 23 c 23 /c 13 Description 1 1 1 All data lie exactly on the best-fit plane distributed in a circular cloud <1 = c 13 1 Data do not lie exactly on a plane, but they project onto the best-fit plane as a circular cloud 1 0 0 Data lie exactly along a line perpendicular to the worst-fit plane 0 0 0 Data are distributed in a spherical cloud 2

TABLE DR4: SLIP LINE ORIENTATIONS ON LOCAL HYPOCENTER ALIGNMENT PLANES FOR SUBSET SOLUTIONS Alignment Plane Subset Solution Slip Line Trend Slip Line Rake Comment Shear Plane Strike Shear Plane Dip Slip Line Plunge Normalized slip Magnitude Southern 2-ne 2-32 141.9 88.5 142.1 6.2 6.2 0.16 2.32 & 2.42 events in the NE planar alignment 2-ne 2-42 141.9 88.5 141.1 28.3 28.3 0.36 2.32 & 2.42 events in the NE planar alignment 2-sw 2-32 142.7 79.3 143.7 5.1 5.2 0.16 All events in SW alignment plane 2-sw 2-42 142.7 79.3 138.7 20.2 20.6 0.29 All events in SW alignment plane 5-pl2 5-52 142.2 76.8 148.3 24.3 25.0 0.53 Events above 9 km except for Data Nos. 72, 433, 874 5 pl2 5-62 142.2 76.8 128.1 46.1 47.7 0.65 Events above 9 km except for Data Nos. 72, 433, 874 5-pl1 5-52 157.5 27.6 129.0 14.0 31.5 0.25 Events below 9 km except for Data No. 433 5-pl1 5-62 157.5 27.6 236.9 27.2 80.6 0.05 Events below 9 km except for Data No. 433 7-02 7-02 146.4 71 146.0 1.1 1.2 0.15 7-02-1 7-02 336.5 25.7 189.3 14.6 144.4 0.60 7-02-2 7-02 322.3 26.1 189.0 19.6 130.2 0.60 7-02-3 7-02 308.0 17.3 195.7 16.1 111.4 0.70 7-02-4 7-02 329.1 11.8 196.8 8.8 131.7 0.73 7-02-5 7-02 135.2 3.4 199.7 3.1 64.5 0.78 Central 4-62-ne 4-52 289.4 90 289.4 47.7 0.0 0.18 ne plane; 4.62 data: slip line for 4.52 solution 4-62-ne 4-62 289.4 90 289.4 3.4 0.0 0.40 ne pl; 4.62 data: slip line for 4.62 solution 4-62-sw 4-52 291.7 89.9 291.6 42.3 42.3 0.20 4.62 data for short sw planar alignment 4-62-sw 4-62 291.7 89.9 291.7 3.5 0.0 0.40 4.62 data for short sw planar alignment 9-1 9-01 124.9 64 271.0 48.9 123.1 0.06 11-33 11-33 117.3 79.7 140.8 65.5 67.6 0.05 11-43 11-43 128.7 60.3 126.7 3.5 4.0 0.11 Northern 12-pl1 12-62 139.4 63.8 230.6 63.8 90.5 0.62 Shallow data along planar alignment from 12.62 and 12.72 12-pl1 12-72 139.4 63.8 158.5 33.6 38.1 0.25 Shallow data along planar alignment from 12.62 and 12.72 12-pl2 12-62 142.3 88.9 212.9 88.8 89.6 0.48 Deep data along planar alignment from 12.62 and 12.72 12-pl2 12-72 142.3 88.9 143.9 55.3 55.3 0.31 Deep data along planar alignment from 12.62 and 12.72 15-pl1 15-732 183.1 8.9 162.0 3.2 21.3 0.18 Defined largely by Subset.732; partly by.932 1 1

15-pl1 15-932 183.1 8.9 246.5 8.0 63.7 0.48 15-pl2 15-732 317.4 17.4 81.6 14.5 123.0 0.30 Defined largely by Subset.732; partly by.932 15-pl2 15-932 317.4 17.4 277.8 11.3 40.9 0.18 15-pl3 15-842 314.1 34 254.4 30.2 64.1 0.52 Defined mostly by subset.842 15-pl4 15-832 151.8 68.2 153.7 4.7 5.1 0.38 Defined by subsets.832; partly by.932 15-pl4 15-932 151.8 68.2 189.4 56.7 64.2 0.26 15-pl5 15-932 327.4 89.9 327.5 38.3 38.3 0.13 Defined by Subset.932, with.732 and.832 15-pl5 15-732 327.4 89.9 327.5 49.9 49.9 0.22 15-pl5 15-832 327.4 89.9 327.4 10.2 0.0 0.36 13 13-33 126.3 49.8 255.9 42.3 118.1 0.40 13 13-43 126.3 49.8 112.3 16.0 21.1 0.15 13 13-53 126.3 49.8 301.9 5.2 173.2 0.55 13-14-pl1 13-33 169.6 81.3 160.1 47.2 47.9 0.21 Sets 13 (& 14) N S alignment on SE of cluster 13-14-pl1 13-43 169.6 81.3 182.7 56.0 57.0 0.44 13-14-pl1 13-53 169.6 81.3 144.5 70.2 72.1 0.21 13-14-pl2 13-33 357.7 88.4 1.3 65.9 66.0 0.25 Sets 13 (& 14) N S alignment on SW of cluster 13-14-pl2 13-43 357.7 88.4 354.0 66.6 66.7 0.55 13-14-pl2 13-53 357.7 88.4 000.3 58.8 58.8 0.25 13-pl3 13-33 168.7 85.4 161.5 57.3 57.6 0.23 Set 13 N S alignment on NE of cluster 13-pl3 13-43 168.7 85.4 176.4 59.0 59.3 0.49 13-pl3 13-53 168.7 85.4 155.8 70.2 70.7 0.23 13-pl4 13-33 340.8 88.8 344.1 69.9 69.9 0.26 Set 13 N S alignment on NW of cluster 13-pl4 13-43 340.8 88.8 338.8 59.3 59.3 0.57 13-pl4 13-53 340.8 88.8 346.0 77.0 77.1 0.22 13-pl5 13-33 077.6 83.6 258.4 7.4 172.6 0.24 E W alignment on the N side of cluster 13-pl5 13-43 077.6 83.6 261.6 31.7 148.1 0.14 13-pl5 13-53 077.6 83.6 257.9 2.3 177.7 0.40 14 14 129.6 54.1 171.7 42.8 57.0 0.11 Set 14 data minus Data Nos 103, 168, 388 16-pl1b 16-11 131.2 14.2 165.3 8.1 34.9 0.17 16-pl1b 16-21 131.2 14.2 272.3 9.0 140.2 0.19 16-pl2b 16-11 129.3 67.9 336.4 48.3 126.3 0.30 16-pl2b 16-21 129.3 67.9 158.8 50.5 56.4 0.26 17-pl1 17-02 226.1 14.3 084.8 9.0 140.4 0.46 Data Nos. 458, 648, 805, 982 omitted 17-pl2 17-02 179.1 66.8 179.7 1.3 1.4 0.17 17-pl3 17-02 106.7 36.6 244.6 26.5 131.6 0.11 2 2

3 Shallow Zone Off plane sets 1-ne 1-32 295.5 88.1 296.6 30.6 30.6 0.64 1-ne 1-42 295.5 88.1 298.7 59.3 59.4 0.44 1-sw1 1-32 115.3 75.5 106.4 31.0 32.1 0.75 1-sw1 1-42 115.3 75.5 103.8 37.6 39.1 0.26 1-sw2 1-32 117.8 58.2 100.1 26.1 31.2 0.79 Mostly 1.32 data; some 1.42 data 1-sw2 1-42 117.8 58.2 131.5 20.9 24.8 0.17 1.42 data that are not aligned with 1_ne plane 6-pl1 6 101.8 1.6 201.8 1.6 100.0 0.33 6-pl2 6 290.3 28.6 012.8 28.4 83.4 0.14 6-pl3 6 111.0 82.6 182.7 82.2 87.6 0.34 Data defining steep plane on NE of data set 6-pl4 6 111.8 80.8 183.3 80.3 86.9 0.37 Shallow data defining shallow NE dipping plane 8-51-pl1 8-51 287.9 55.4 012.0 55.3 86.6 0.21 Shallow data defining NE dipping plane 8-51-pl2 8-51 106.7 86.4 112.8 59.4 59.6 0.10 Data defining vertical SW boundary to 8_51 subset 8-61-pl3 8-61 109.9 57.7 260.0 38.2 132.9 0.28 Data defining moderately SW dipping plane in 8_61 subset 10-pl2 10-52 134.2 44 154.2 18.3 26.9 0.43 Data define SW dipping plane at top of hypocenter distribution 10-52-pl1 10-52 315.4 83.6 318.5 25.6 25.8 0.27 Data selected to define vertical alignment plane 10-52-pl1 10-62 315.4 83.6 326.0 58.7 59.3 0.42 10_52_vert data 10-pl2 10-62 134.2 44 216.8 43.8 84.7 0.29 10_sw_pl data 10-pl3 10-52 151.6 82.3 152.7 7.9 8.0 0.37 Data from NW cluster; steep off plane group 10-pl3 10-62 151.6 82.3 146.0 35.7 36.1 0.25 Data from NW cluster; steep off plane group 12_50 12-5 153.1 81.8 338.8 34.6 145.0 0.09 3_se 3-1 326.8 72.5 349.4 50.7 54.2 0.47 Events in the SE cluster 3_se 3-2 326.8 72.5 354.9 56.2 60.6 0.64 Events in the SE cluster 3_nw 3-1 297.9 44.6 039.7 44.0 98.5 0.16 All data from the NW scattered events not in the SE cluster 3_nw 3-2 297.9 44.6 022.0 44.4 85.8 0.86 All data from the NW scattered events not in the SE cluster 3