UCERF3

The 2015 Uniform California Earthquake Rupture Forecast, Version 3, or UCERF3, is the latest official earthquake rupture forecast (ERF) for the state of California, superseding UCERF2. It provides authoritative estimates of the likelihood and severity of potentially damaging earthquake ruptures in the long- and near-term. Combining this with ground motion models produces estimates of the severity of ground shaking that can be expected during a given period (seismic hazard), and of the threat to the built environment (seismic risk). This information is used to inform engineering design and building codes, planning for disaster, and evaluating whether earthquake insurance premiums are sufficient for the prospective losses. A variety of hazard metrics can be calculated with UCERF3; a typical metric is the likelihood of a magnitude M 6.7 earthquake (the size of the 1994 Northridge earthquake) in the 30 years (typical life of a mortgage) since 2014.

UCERF3 was prepared by the Working Group on California Earthquake Probabilities (WGCEP), a collaboration between the United States Geological Survey (USGS), the California Geological Survey (CGS), and the Southern California Earthquake Center (SCEC), with significant funding from the California Earthquake Authority (CEA).



Highlights
A major achievement of UCERF3 is use of a new methodology that can model multifault ruptures such as have been observed in recent earthquakes. This allows seismicity to be distributed in a more realistic manner, which has corrected a problem with prior studies that overpredicted earthquakes of moderate size (between magnitude 6.5 and 7.0). The rate of earthquakes of magnitude (M ) 6.7 and greater (over the entire state) is now believed to be about one in 6.3 years, instead of one in 4.8 years. On the other hand, earthquakes of magnitude 8 and larger are now expected about every 494 years (down from 617). Otherwise the overall expectations of seismicity are generally in line with earlier results. (See Table A for a summary of the overall rates.)

The fault model database has been revised and expanded to cover over 350 fault sections, up from about 200 for UCERF2, and new attributes added to better characterize the faults. Various technical improvements have also been made.


 * {| class="wikitable"

! M≥ ! 6.0 ! 6.7 ! 7.0 ! 7.5 ! 7.7 ! 8.0
 * + Table A: Estimated probabilities (minimum,, and maximum) of an earthquake of the given magnitude in the next thirty years for different regions of California1
 * -!Region
 * All CA
 * 100%   100%
 * 97%   100%
 * 77%    100%
 * 17%     85%
 * 3%     71%
 * 0%     32%
 * N. CA
 * 100%   100%
 * 84%    100%
 * 55%     96%
 * 8%     60%
 * 1%     45%
 * 0%     25%
 * S. CA
 * 100%   100%
 * 77%    100%
 * 44%     97%
 * 9%     79%
 * 2%     68%
 * 0%     32%
 * SF
 * 89%    100%
 * 52%     94%
 * 27%     84%
 * 5%     43%
 * 0%     32%
 * 0%     21%
 * LA
 * 84%    100%
 * 28%     92%
 * 17%     87%
 * 5%     77%
 * 1%     68%
 * 0%     32%
 * }
 * 1. From Table 7 in . "M" is moment magnitude (p. 512).
 * 5%     77%
 * 1%     68%
 * 0%     32%
 * }
 * 1. From Table 7 in . "M" is moment magnitude (p. 512).

Of the six main faults evaluated in previous studies the Southern San Andreas fault remains the most likely to experience an M ≥ 6.7 earthquake in the next 30 years. The largest increase in such likelihood is on the Calaveras fault (see main faults map for location), where the mean (most likely) value is now set at 25%. The old value, of 8%, is less than the minimum now expected (10%). The previous under-estimate is believed to be due mostly to not modeling multifault ruptures, which limited the size of many ruptures.

The largest probability decrease is on the San Jacinto fault, which went from 32% to 9%. Again this is due to multifault rupturing, but here the effect is fewer earthquakes, but they are more likely to be bigger (M ≥ 7.7)

Table B

 * {| class=wikitable

!Fault2 !Section maps3 !QFFDB fault#4 !Length5 !Notable Earthquakes !Min.6 ||Mean ||Max.
 * + Table B: Aggregate probabilities of an M ≥ 6.7 earthquake within 30 years (and change from UCERF2)1

Parkfield
 * San Andreas Fault south
 * San Andreas Fault south

Cholame

Carrizo

Big Bend

Mojave N

Mojave S

San Bernardino N

San Bernardino S

San Gorgonio Pass

N. Branch Mill Cr

Coachella

1f

1g

1h

1i

1j

546 km 339 miles 1857 Fort Tejon earthquake Offshore
 * 17% (−6%)|| (−7%)|| 93% (−1%)
 * San Andreas Fault north
 * San Andreas Fault north

North Coast

Peninsula

Santa Cruz Mts

Creeping Section

1a

1b

1c

1d

1e

472 km 293 miles 1906 San Francisco earthquake Rodgers Creek
 * 1% (−5%) || (+12%) || 73% (+33%)
 * Hayward/ Rodgers Creek Fault
 * Hayward/ Rodgers Creek Fault

Hayward North

Hayward South

55a

55b

55c

32 150 km 93 miles 1868 Hayward earthquake
 * 14% (−2%) || (0%) || 54% (−14%)

North
 * Calaveras Fault
 * Calaveras Fault

Central

South

54a

54b

54c

54d

123 km 76 miles 1911 Calaveras earthquake

1979 Coyote Lake earthquake

1984 Morgan Hill earthquake

2007 Alum Rock earthquake San Bernardino
 * 10% (+8%) || (+17%) || 54% (+32%)
 * San Jacinto Fault Zone
 * San Jacinto Fault Zone

San Jacinto Valley

Stepovers

Anza

Clark

Coyote Creek

Borrego

Superstition Mtn

125a

125b

125c

125d

125e

125f

125g

309 km 192 miles 1918 San Jacinto earthquake East
 * 0% (−14%) || (23%) || 35% (−20%)
 * Garlock Fault
 * Garlock Fault

Central

West

69a

69b

69c

254 km 158 miles Whittier
 * 0% (−3%) || (+2%) || 37% (+24%)
 * Elsinore Fault Zone
 * Elsinore Fault Zone
 * Elsinore Fault Zone

Glen Ivy

Stepovers

Temecula

Julian

Coyote Mountains

126a

126b

126c

126d

126e

126f

126g

249 km 217 miles 1910 Elsinore earthquake
 * 1% (−4%)|| (−6%) || 17% (−8%)


 * }


 * Notes.
 * 1. Adapted from Table 6 in . Values are aggregated from the fault sections comprising each fault. Some sections have higher individual probabilities; see Table 4 in . "M" is moment magnitude (p. 512).
 * 2. These are the six faults for which UCERF2 had enough data to do stress-renewal modeling. The Hayward fault zone and Rodgers Creek fault are treated as a single fault; the San Andreas fault is treated as two sections.
 * 3. UCEF3 fault sections, with links to "participation" maps for each section (outlined in black), showing the rate (in color) that section participates in ruptures with other sections. Participation maps for all fault sections available at http://pubs.usgs.gov/of/2013/1165/data/UCERF3_SupplementalFiles/UCERF3.3/Model/FaultParticipation/ Some faults have had sections added or split since UCERF2.
 * 4. USGS Quaternary Fault and Fold Database fault numbers, with links to summary reports. QFFDB maps are no longer available.
 * 5. Lengths from UCERF-2, Table 4; may vary from QFFDB values.
 * 6. Min. and Max. probabilities correspond to the least and most likely alternatives in the logic tree; the Mean is a weighted average.
 * 7. Slip-rates not included due to variation across sections and deformation models. See figure C21 (below) for an illustration.

Methodology
California earthquakes result from the Pacific Plate, heading approximately northwest, sliding past the North American continent. This requires accommodation of 34 to 48 millimeters (about one and a half inches) of slippage per year, with some of that taken up in portions of the Basin and Range Province to the east of California. This slippage is accommodated by ruptures (earthquakes) and aseismic creep on the various faults, with the frequency of ruptures dependent (in part) on how the slippage is distributed across the various faults.

Modeling
Like its predecessor, UCERF3 determines this based on four layers of modeling:
 * 1) The fault models (FM 3.1 and 3.2) describe the physical geometry of the larger and more active faults.
 * 2) Deformation models determine the slip rates and related factors for each fault section, how much strain accumulates before a fault ruptures, and how much energy is then released. Four deformation models are used, reflecting different approaches to handling earthquake dynamics.
 * 3) The earthquake rate model (ERM) fits together all this data to estimate the long-term rate of rupturing.
 * 4) The probability model estimates how close (ready) each fault segment is to rupturing given how much stress has accumulated since its last rupture.

The first three layers of modeling are used to determine the long-term, or Time Independent, estimates of the magnitude, location, and frequency of potentially damaging earthquakes in California. The Time Dependent model is based on the theory of elastic rebound, that after an earthquake releases tectonic stress there will be some time before sufficient stress accumulates to cause another earthquake. In theory, this should produce some regularity in the earthquakes on a given fault, and knowing the date of the last rupture is a clue to how soon the next one can be expected. In practice this is not so clear, in part because slip rates vary, and also because fault segments influence each other, so a rupture on one segment triggers rupturing on adjacent segments. One of the achievements of UCERF3 is to better handle such multifault ruptures.

The various alternatives (see diagram), taken in different combinations, form a logic tree of 1440 branches for the Time Independent model, and, when the four probability models are factored in, 5760 branches for the Time Dependent model. Each branch was evaluated and weighted according to its relative probability and importance. The UCERF3 results are an average of all these weighted alternatives.

"The Grand Inversion"
In UCERF2 each fault was modeled separately, as if ruptures do not extend to other faults. This assumption of fault segmentation was suspected as the cause of UCERF2 predicting nearly twice as many earthquakes in the M 6.5 to 7.0 range then actually observed, and is contrary to the multifault rupturing seen in many earthquakes.

UCERF3 subdivides each fault section (as modeled by the Fault Models) into subsections (2606 segments for FM 3.1, and 2665 for FM 3.2), then considers ruptures of multiple segments regardless of which parent fault they belong to. After removing those ruptures considered implausible there are 253,706 possibilities to consider for FM 3.1, and 305,709 for FM 3.2. This compares to less than 8,000 ruptures considered in UCERF2, and reflects the high connectivity of California's fault system.

A significant achievement of UCERF is development of system-level approach called the "grand inversion". This uses a supercomputer to solve a system of linear equations that simultaneously satisfies multiple constraints such as known slip rates, etc. The result is a model (set of values) that best fits the available data. In balancing these various factors it also provides an estimate of how much seismicity is not accounted for in the fault model, possibly in faults not yet discovered. The amount of slip occurring on unidentified faults has been estimated at between 5 and about 20 mm/yr depending on the location (generally higher in the LA area) and deformation model, with one model reaching 30 mm/yr just north of LA.

Assessment
While UCERF3 represents a considerable improvement over UCERF2, and the best available science to-date for estimating California's earthquake hazard, the authors caution that it remains an approximation of the natural system. There are a number of assumptions in the Time Independent model, while the final (Time Dependent) model explicitly "assumes elastic rebound dominates other known and suspected processes that are not included in the model." Among the known processes not included is spatiotemporal clustering.

There are a number of sources of uncertainty, such as insufficient knowledge of fault geometry (especially at depth) and slip rates, and there is considerable challenge in how to balance the various elements of the model to achieve the best fit with the available observations. For example, there is difficulty fitting paleoseismic data and slip rates on the southern San Andreas Fault, resulting in estimates of seismicity that run about 25% less than seen in the paleoseismic data. The data does fit if a certain constraint (the regional Magnitude-Frequency Distribution) is relaxed, but this brings back the problem over-predicting moderate events.

An important result is that the generally accepted Gutenberg-Richter (GR) relationship (that the distribution of earthquakes shows a certain relationship between magnitude and frequency) is inconsistent with certain parts of the current UCERF3 model. The model implies that achieving GR consistency would require certain changes in seismological understanding that "fall outside the current bounds of consensus-level acceptability". Whether the Gutenberg-Richter relation is inapplicable at the scale of individual faults, or some basis of the model is incorrect, "will be equally profound scientifically, and quite consequential with respect to hazard."