The 2014 Working Group on California Earthquake Probabilities (WGCEP14) present the time-independent component of the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3), which provides authoritative estimates of the magnitude, location, and time-averaged frequency of potentially damaging earthquakes in California. The primary achievements have been to relax fault segmentation and include multifault ruptures, both limitations of UCERF2. The rates of all earthquakes are solved for simultaneously and from a broader range of data, using a system-level inversion that is both conceptually simple and extensible. The inverse problem is large and underdetermined, so a range of models is sampled using an efficient simulated annealing algorithm. The approach is more derivative than prescriptive (e.g., magnitude-frequency distributions are no longer assumed), so new analysis tools were developed for exploring solutions. Epistemic uncertainties were also accounted for using 1440 alternative logic-tree branches, necessitating access to supercomputers. The most influential uncertainties include alternative deformation models (fault slip rates), a new smoothed seismicity algorithm, alternative values for the total rate of M w ≥ 5 events, and different scaling relationships, virtually all of which are new. As a notable first, three deformation models are based on kinematically consistent inversions of geodetic and geologic data, also providing slip-rate constraints on faults previously excluded due to lack of geologic data. The grand inversion constitutes a system-level framework for testing hypotheses and balancing the influence of different experts. For example, we demonstrate serious challenges with the Gutenberg-Richter hypothesis for individual faults. UCERF3 is still an approximation of the system, however, and the range of models is limited (e.g., constrained to stay close to UCERF2). Nevertheless, UCERF3 removes the apparent UCERF2 overprediction of M 6.5-7 earthquake rates and also includes types of multifault ruptures seen in nature. Although UCERF3 fits the data better than UCERF2 overall, there may be areas that warrant further site-specific investigation. Supporting products may be of general interest, and we list key assumptions and avenues for future model improvements. Manuscript OrganizationBecause of manuscript length and model complexity, we begin with an outline of this report to help readers navigate the various sections:

After the 1987 Whittier Narrows and 1994 Northridge earthquakes revealed that blind thrust faults represent a significant threat to metropolitan Los Angeles, a network of 250 continuously recording global positioning system (GPS) stations was deployed to monitor displacements associated with deep slip on both blind and surface faults. Here we augment this GPS data with interferometric synthetic aperture radar imagery to take into account the deformation associated with groundwater pumping and strike-slip faulting. After removing these non-tectonic signals, we are left with 4.4 mm yr-1 of uniaxial contraction across the Los Angeles basin, oriented N 36 degrees E (perpendicular to the major strike-slip faults in the area). This indicates that the contraction is primarily accommodated on thrust faults rather than on the northeast-trending strike-slip faults. We have found that widespread groundwater and oil pumping obscures and in some cases mimics the tectonic signals expected from the blind thrust faults. In the 40-km-long Santa Ana basin, groundwater withdrawal and re-injection produces 12 mm yr-1 of long-term subsidence, accompanied by an unprecedented seasonal oscillation of 55 mm in the vertical direction and 7 mm horizontally.

Two recent large earthquakes in the Mojave Desert, California-the magnitude 7.3 1992 Landers and magnitude 7.1 1999 Hector Mine earthquakes-have each been followed by elevated crustal strain rates over periods of months and years. Geodetic data collected after the Hector Mine earthquake exhibit a temporally decaying horizontal velocity field and a quadrant uplift pattern opposite to that expected for localized shear beneath the earthquake rupture. We interpret the origin of this accelerated crustal deformation to be vigorous flow in the upper mantle in response to the stress changes generated by the earthquake. Our results suggest that transient flow in the upper mantle is a fundamental component of the earthquake cycle and that the lower crust is a coherent stress guide coupling the upper crust with the upper mantle.

Two contrasting models of the earthquake deformation cycle on strike slip faults predict significant temporal declines in shear strain rate near the fault, accompanied by a progressive broadening of the zone of deformation adjacent to it. In the thin lithosphere model, transient deformation results from flow in the asthenosphere due to stress relaxation following faulting through most or all of the lithosphere. For an earth model with a thick elastic lithosphere (plate thickness » depth of seismic slip), transient motions are due to postearthquake aseisrnic slip below the coseismic fault plane. Data from the San Andreas fault indicate a long‐term temporal decrease in strain rate that persists for at least 30 years and may extend through the entire earthquake cycle. Observations support a cycle‐long rate decrease and a temporal spreading of the deformation profile only if movement cycles on the northern and southern locked sections of the fault are basically similar. If so, the usually lower strain rates and broader deformation zone currently observed on the southern San Andreas represent a later evolutionary stage of the northern locked section, where a great earthquake is a more recent occurrence. Although the data allow some extreme models to be discarded, no sufficiently strong constraints exist to decide between the thin and thick lithosphere models. Regardless of the appropriate model the geodetic observations themselves indicate that strain buildup is sufficiently nonlinear to cause significant departures from recurrence estimates based on linear strain accumulation and the time‐predictable model.

The distribution of deformation within the Basin and Range province was determined from 1992, 1996, and 1998 surveys of a dense, 800-kilometer-aperture, Global Positioning System network. Internal deformation generally follows the pattern of Holocene fault distribution and is concentrated near the western extremity of the province, with lesser amounts focused near the eastern boundary. Little net deformation occurs across the central 500 kilometers of the network in western Utah and eastern Nevada. Concentration of deformation adjacent to the rigid Sierra Nevada block indicates that external plate-driving forces play an important role in driving deformation, modulating the extensional stress field generated by internal buoyancy forces that are due to lateral density gradients and topography near the province boundaries.