How Risk Selection in California Can Be Affected by Earthquake Source Modeling Assumptions: A Spatial Loss Correlation Perspective
Delphine FitzenzSeptember 27, 2018
There used to be several ways to ensure risk diversification in a California earthquake insurance portfolio. You could select risks on the Peninsula and risks in the East Bay; or select risks in Ventura and Orange counties; or risks in Santa Barbara and Los Angeles counties. Better yet, it was considered that selecting risks in the San Francisco Bay Area and in the Los Angeles region was a perfectly good way of achieving risk diversification. This practice was largely based on an understanding of the spatial correlation of expected loss between counties in California and selecting risks for counties which decreased loss correlations in the insured portfolio.
While California and the large-scale plate motions that it is subjected to have not changed in recent years, the way earthquake sources are modeled has. The two main areas scientists are trying to explore are: first, whether there are preferential locations in a fault network where ruptures are likely to start or stop. The second area examines what the relationship is, if any, between the timing of the latest events on a fault network and the timing of the next event that will overlap with those events. A third avenue of research that is relevant for California is the behavior of aseismic faults — faults that deform without making felt earthquakes, and what happens to them when large ruptures propagate in their direction.
RMS led a study to quantify the impact of these three major modeling assumptions on spatial loss correlations. The study used sixteen county portfolios made using the RMS Industry Exposure Database (2017), and two vintages of source model: the Uniform California Earthquake Rupture Forecast 2 and 3 (UCERF2 and UCERF3). One major conclusion was that new and different risk selection strategies would be required by the spatial loss correlation study to ensure portfolio diversification with the most recent United States Geological Survey (USGS) model (UCERF3) as compared to the previous versions of the model (i.e. UCERF1 or 2).
RMS shared this study with the broader earthquake engineering community, including those researchers who work on source models. The hope is that the study findings can inform them on the societal and commercial implications of their model assumptions, but also that it can help discriminate between source models that yield similar hazard maps at given return periods.
A Word on Loss Correlation
Common events for each portfolio drive up the value of the correlation coefficient. So, perfect correlation (+1.0) means that the list of events, event rates and secondary uncertainties between two portfolios will be identical.
Change in spatial correlation can come from 1) change in number of events that affect both portfolios (longer/shorter ruptures than before, aka segmentation assumption), but also from 2) change in the rate of events in common (especially for different time-dependent models).
Models up to 2008 (including UCERF2) had most ruptures start and stop at fault segment boundaries (i.e., discontinuities visible in the landscape), with a few exceptions along the San Andreas fault, for example. The most recent models (2014, UCERF3) have explored the opposite assumption: ruptures can start and stop anywhere in the fault network. In this latter case, how long and which faults the rupture is going to span depends only on a few criteria.
Instead of relying on past observed ruptures and generic assumptions about ruptures being confined to fault segments, scientists have decided to include all ruptures that could not be ruled out. That is a lot of ruptures, most of them very similar to each other, with a lot of overlap. Most of them are very long, and many involve more than one fault. Some others are shorter than individual fault segments. Since loss correlation between portfolios is driven by those events that impact both sets of exposure, one can see how this change in segmentation assumption can correlate previously uncorrelated exposures, or on the contrary reduce the correlation.
Effects of removing segment boundaries as constraints for rupture size and location:
If two counties are each close to a different fault segment and now a rupture can connect those segments, correlation is going to increase. e.g., increase in correlation between Santa Barbara and Los Angeles, for example.
On the contrary, if each end of a fault segment can rupture alone, two sets of exposure at each end of the same fault segment do not have to be correlated anymore. e.g., decrease in correlation between San Francisco county and Santa Clara county.
Pre-UCERF3, most ruptures were confined to fault segments and the time-dependent (TD) model was applied at the scale of fault segments and only for more active and better-known faults. The largest component of the time-dependent model was a renewal model, defined by a mean recurrence time and a deviation from periodic behavior (aka aperiodicity).
The aperiodicity is not well constrained because we do not often have the dates for many consecutive large events on a fault. There are two main ways around that: 1) to use data from many faults around the world, and 2) to use numerical models simulating earthquake catalogs in a complex fault network. A third way involves using all available earthquake geology data to infer probability distributions for aperiodicity given a chosen recurrence model.
Until 2008, the Working Group on California Earthquake Probabilities (WGCEP) used the first option. The global average pointed to an aperiodicity of 0.5 being the most probable, with 0.3 (slightly more periodic) and 0.7 (slightly less periodic) also possible choices.
Effect of high aperiodicity on few single-segment ruptures on spatial loss correlation — UCERF2 case:
Minimal effect in Southern California, slight decrease between the Peninsula counties, and slight increase in the East Bay counties (reflecting rate changes between time-independent and time-dependent perspective).
For UCERF3, the WGCEP chose the second option, using results emerging from simulations of the fault network geometry representing California. Those exhibited a trend as a function of magnitude: the larger the magnitude, the more periodic the behavior. For most magnitudes, the aperiodicity is lower in UCERF3 than in UCERF2 (ruptures are modeled as more periodic in UCERF3 than in UCERF2).
Once coupled to the unsegmented approach, one can see how the largest, multi-segment and multi-fault ruptures become the most periodic ones as well. This does not have much impact on ruptures whose fault sections do not have a known time since the last event. However, for ruptures where we have a good understanding of time since the last event, such as the Hayward Fault, this impact is pronounced.
The probability gain relative to the time-independent model is larger for the larger magnitude, multi-segment events involving the Hayward Fault than for those events confined to the Hayward Fault alone (by virtue of the smaller aperiodicity for larger events). The same thing (but with a negative sign) happens for events overlapping with the 1906 event on the northern San Andreas Fault, In other words, the effect of time-dependence is also larger for the larger events, but the effect is to decrease their rate (the time since last event is still small compared to the mean recurrence time, an effect also called stress shadow).
New aperiodicity scheme + TD for all fault events + unsegmented approach: UCERF3 time-dependent model relative to UCERF3 time-independent.
Loss correlations increase between the North Bay and the South Bay (through the East Bay), but they decrease between Marin County and the peninsula counties.
Participation of the Central San Andreas Fault to Large Ruptures Linking Up Northern California Faults and Southern California Faults
Some faults are mostly aseismic, meaning they are known to deform without making felt or damaging earthquake. The San Andreas Fault, south of the San Francisco Bay Area but north of Parkfield exhibits such behavior. While scientists agree that the measured aseismic deformation across this fault does not account for all plate motion, it was previously thought that only the asperities (i.e., the parts that are stuck) could produce/participate in earthquakes. Tohoku in 2011 has shown that aseismic portions of the interface can participate in seismic ruptures if the ruptures are fast and large enough when they get to the (usually) aseismic area. In UCERF3, the WGCEP has assigned 20 percent of the accumulated moment rate for seismic events going through that section of the San Andreas fault.
Large ruptures through the central “creeping” section of the San Andreas
Events from high M5 to M>8.0 are allowed to go through that section and therefore correlate San Francisco Bay Area exposure with Los Angeles region exposure
Between San Bernardino or Riverside in the south, and Marin, San Francisco, San Mateo or Santa Clara in the north, loss correlation is up to 15 percent in the time-independent perspective and up to 10 percent in the time-dependent perspective in UCERF3*, whereas it was zero percent in previous models.
* Residential industry exposure 2017.
To highlight and quantify those largely unstudied effects of source model assumptions on spatial loss correlations in California, RMS used RiskLink 16 (UCERF2 implementation) and RiskLink 17 (UCERF3 implementation), as well as hybrid versions between the two which kept all components uniform except the source model, and presented and published our findings at the National Conference on Earthquake Engineering in Los Angeles in June 2018. The article from the conference proceedings can be made available upon request. While the main trends are presented here, more cases and more details are discussed in the proceedings.
Notations: UCERF Uniform California Earthquake Rupture Forecast, RiskLink version 16 (based on UCERF2 in terms of California on-fault source model); RiskLink version 17 (based on UCERF3 in terms of California on-fault source model).
Delphine Fitzenz works on earthquake source modeling for risk products, with a particular emphasis on spatio-temporal patterns of large earthquakes. Since joining RMS in 2012 after 10+ years in academia, she has strived to bring the risk and the earthquake science communities closer together through her articles and by organizing special sessions at conferences.
Delphine holds an Eng. MSc and a research MSc in geophysics from University Louis Pasteur, Strasbourg, France, and a Ph.D. from ETH Zurich, Switzerland. Her research has been on two major themes: the role of fluids and fault structure on the behavior of fault systems, and constraints on earthquake recurrence probabilistic models from individual segments to system-wide behavior. She is currently on the Technical Advisory Group for the National Seismic Hazard update project in New Zealand and on the Editorial Board of the Bulletin of the Seismological Society of America (BSSA).