CyberShake Study 21.12

CyberShake 21.12 is a computational study to use a new ERF with CyberShake, generated from an RSQSim catalog. We plan to calculate results for 336 sites in Southern California using the RSQSim ERF, a minimum Vs of 500 m/s, and a frequency of 1 Hz. We will use the CVM-S4.26.M01 model, and the GPU implementation of AWP-ODC-SGT enhanced from the BBP verification testing. We will begin by generating all sets of SGTs, on Summit, then post-process them on a combination of Summit and Frontera.

Status

This study is in the pre-production phase. Production is scheduled to begin in mid-December, 2021.

Science Goals

The science goals for this study are:

• Calculate a regional CyberShake model using an alternative, RSQSim-derived ERF.
• Compare results from an RSQSim ERF to results using a UCERF2 ERF (Study 15.4).

Technical Goals

The technical goals for this study are:

• Perform a study using OLCF Summit as a key compute resource.
• Evaluate the performance of the new workflow submission host, shock-carc.
• Use Globus Online for staging of output data products.

ERF

The ERF was generated from an RSQSim catalog, with the following parameters:

• 715kyr catalog (the first 65k years of events were dropped, so that every fault's first event is excluded)
• 220,927 earthquakes with M6.5+
• All events have equal probability, 1/715k

Additional details are available on the catalog's metadata page.

Sites

We will run a list of 335 sites, taken from the site list that was used in other Southern California studies.

• CSV site list
• KML site list with names
• KML site list without names

Velocity Model

We will use CVM-S4.26.M01.

To better represent the near-surface layer, we will populate the velocity parameters for the surface point by querying the velocity model at a depth of (grid spacing)/4. For this study, the grid spacing is 100m, so we will query UCVM at a depth of 25m and use that value to populate the surface grid point. The rationale is that the media parameters at the surface grid point are supposed to represent the material properties for [0, 50m], and this is better represented by using the value at 25m than the value at 0m.

Output Data Products

We plan to produce the following data products:

• PSA: X and Y spectral acceleration at 44 periods (10, 9.5, 9, 8.5, 8, 7.5, 7, 6.5, 6, 5.5, 5, 4.8, 4.6, 4.4, 4.2, 4, 3.8, 3.6, 3.4, 3.2, 2, 2.8, 2.6, 2.4, 2.2, 2, 1.66667, 1.42857, 1.25, 1.11111, 1, .66667, .5, .4, .33333, .285714, .25, .22222, .2, .16667, .142857, .125, .11111, .1 sec)
• RotD: RotD50, the RotD50 azimuth, and RotD100 at 22 periods (1.0, 1.2, 1.4, 1.5, 1.6, 1.8, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.5, 4.0, 4.4, 5.0, 5.5, 6.0, 6.5, 7.5, 8.5, 10.0 sec)
• Durations: for X and Y components, energy integral, Arias intensity, cumulative absolute velocity (CAV), and for both velocity and acceleration, 5-75%, 5-95%, and 20-80%.

We plan to store the following data products in the database:

• PSA: none
• RotD: RotD50 and RotD100 at 10, 7.5, 5, 4, 3, and 2 sec.
• Durations: acceleration 5-75% and 5-95% for X and Y components