Long-term fault slip rates, distributed deformation rates, and forecast of seismicity in the western United States from joint fitting of community geologic, geodetic, and stress direction data sets

The long-term average velocity field of the western United States is computed with a kinematic finite element code. Community data sets include fault traces, geologic offset rates, geodetic velocities, principal stress directions, and Euler poles. There is an irreducible minimum amount of distributed permanent deformation, which accommodates one third of Pacific-North America relative motion in California. Much of this may be due to slip on faults not included in the model. All data sets are fit at a common RMS level of 1.8 datum standard deviations. Experiments with alternate weights, fault sets, and Euler poles define a suite of acceptable community models. In pseudoprospective tests, fault offset rates are compared to 126 additional published rates not used in the computation: 44% are consistent; another 48% have discrepancies under 1 mm a(-1), and 8% have larger discrepancies. Updated models are then computed. Novel predictions include dextral slip at 2-3 mm a(-1) in the Brothers fault zone, two alternative solutions for the Mendocino triple junction, slower slip on some trains of the San Andreas fault than in recent hazard models, and clockwise rotation of some domains in the eastern California shear zone. Long-term seismicity is computed by assigning each fault and finite element the seismicity parameters (coupled thickness, corner magnitude, and spectral slope) of the most comparable type of plate boundary. This long-term seismicity forecast is retrospectively compared to instrumental seismicity. The western United States has been 37% below its long-term average seismicity during 1977-2008, primarily because of (temporary) reduced activity in the Cascadia subduction zone and San Andreas fault system.