Difference between revisions of "CyberShake Study 22.12"
Line 33: | Line 33: | ||
| [[File:cvmsi_default_LAPD.png|thumb|300px]] | | [[File:cvmsi_default_LAPD.png|thumb|300px]] | ||
| [[File:cvmsi_elygtl_ely_z_700m_LAPD.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_LAPD.png|thumb|300px]] | ||
+ | | [[File:LAPD_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! s650 | ! s650 | ||
Line 41: | Line 42: | ||
| [[File:cvmsi_default_s764.png|thumb|300px]] | | [[File:cvmsi_default_s764.png|thumb|300px]] | ||
| [[File:cvmsi_elygtl_ely_z_700m_s764.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_s764.png|thumb|300px]] | ||
+ | | [[File:s764_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! s568 | ! s568 | ||
| [[File:cvmsi_default_s568.png|thumb|300px]] | | [[File:cvmsi_default_s568.png|thumb|300px]] | ||
| [[File:cvmsi_elygtl_ely_z_700m_s568.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_s568.png|thumb|300px]] | ||
+ | | [[File:s568_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! s035 | ! s035 | ||
| [[File:cvmsi_default_s035.png|thumb|300px]] | | [[File:cvmsi_default_s035.png|thumb|300px]] | ||
| [[File:cvmsi_elygtl_ely_z_700m_s035.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_s035.png|thumb|300px]] | ||
+ | | [[File:s035_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! RIDG | ! RIDG | ||
Line 57: | Line 61: | ||
| [[File:cvmsi_default_PERR.png|thumb|300px]] | | [[File:cvmsi_default_PERR.png|thumb|300px]] | ||
| [[File:cvmsi_elygtl_ely_z_700m_PERR.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_PERR.png|thumb|300px]] | ||
+ | | [[File:PERR_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! MRVY | ! MRVY | ||
| [[File:cvmsi_default_MRVY.png|thumb|300px]] | | [[File:cvmsi_default_MRVY.png|thumb|300px]] | ||
| [[File:cvmsi_elygtl_ely_z_700m_MRVY.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_MRVY.png|thumb|300px]] | ||
+ | | [[File:MRVY_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! s211 | ! s211 | ||
| [[File:cvmsi_default_s211.png|thumb|300px]] | | [[File:cvmsi_default_s211.png|thumb|300px]] | ||
| [[File:cvmsi_elygtl_ely_z_700m_s211.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_s211.png|thumb|300px]] | ||
+ | | [[File:s211_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! PIBU | ! PIBU | ||
Line 79: | Line 86: | ||
! s117 | ! s117 | ||
| [[File:cvmsi_default_s117.png|thumb|300px]] | | [[File:cvmsi_default_s117.png|thumb|300px]] | ||
− | |||
| [[File:cvmsi_elygtl_ely_z_700m_s117.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_s117.png|thumb|300px]] | ||
+ | | [[File:s117_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! USC | ! USC | ||
| [[File:cvmsi_default_USC.png|thumb|300px]] | | [[File:cvmsi_default_USC.png|thumb|300px]] | ||
− | |||
| [[File:cvmsi_elygtl_ely_z_700m_USC.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_USC.png|thumb|300px]] | ||
+ | | [[File:USC_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! WNGC | ! WNGC | ||
| [[File:cvmsi_default_WNGC.png|thumb|300px]] | | [[File:cvmsi_default_WNGC.png|thumb|300px]] | ||
− | |||
| [[File:cvmsi_elygtl_ely_z_700m_WNGC.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_WNGC.png|thumb|300px]] | ||
+ | | [[File:WNGC_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! PEDL | ! PEDL | ||
| [[File:cvmsi_default_PEDL.png|thumb|300px]] | | [[File:cvmsi_default_PEDL.png|thumb|300px]] | ||
− | |||
| [[File:cvmsi_elygtl_ely_z_700m_PEDL.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_PEDL.png|thumb|300px]] | ||
|- | |- | ||
! SBSM | ! SBSM | ||
| [[File:cvmsi_default_SBSM.png|thumb|300px]] | | [[File:cvmsi_default_SBSM.png|thumb|300px]] | ||
− | |||
| [[File:cvmsi_elygtl_ely_z_700m_SBSM.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_SBSM.png|thumb|300px]] | ||
+ | | [[File:SBSM_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
! SVD | ! SVD | ||
| [[File:cvmsi_default_SVD.png|thumb|300px]] | | [[File:cvmsi_default_SVD.png|thumb|300px]] | ||
− | |||
| [[File:cvmsi_elygtl_ely_z_700m_SVD.png|thumb|300px]] | | [[File:cvmsi_elygtl_ely_z_700m_SVD.png|thumb|300px]] | ||
+ | | [[File:SVD_ifless_profile.png|thumb|300px]] | ||
|- | |- | ||
|} | |} |
Revision as of 03:48, 4 October 2022
CyberShake Study 22.10 is a proposed study in Southern California which will include deterministic low-frequency (0-1 Hz) and stochastic high-frequency (1-50 Hz) simulations. We will use the Graves & Pitarka (2019) rupture generator and the high frequency modules from the SCEC Broadband Platform v22.4.
Contents
- 1 Status
- 2 Data Products
- 3 Science Goals
- 4 Technical Goals
- 5 Sites
- 6 Velocity Model
- 7 High-frequency codes
- 8 Spectral Content around 1 Hz
- 9 Verification
- 10 Updates and Enhancements
- 11 Output Data Products
- 12 Computational and Data Estimates
- 13 Lessons Learned
- 14 Stress Test
- 15 Performance Metrics
- 16 Production Checklist
- 17 Presentations, Posters, and Papers
Status
This study is in the planning phase. We estimate calculations will begin in October 2022.
Data Products
Science Goals
Technical Goals
Sites
Velocity Model
We are planning to use CVM-S4.26 with a GTL applied, and the CVM-S4 1D background model outside of the region boundaries.
We are investigating applying the Ely-Jordan GTL down to 700 m instead of the default of 350m. We extracted profiles for a series of southern California CyberShake sites, with no GTL, a GTL applied down to 350m, and a GTL applied down to 700m.
Sites (the CVM-S4.26 basins are outlined in red):
These are for sites outside of the CVM-S4.26 basins:
Site | No GTL | 700m GTL | Smaller value, extracted from mesh |
---|---|---|---|
LAPD | |||
s650 | |||
s764 | |||
s568 | |||
s035 | |||
RIDG | |||
PERR | |||
MRVY | |||
s211 | |||
PIBU |
These are for sites inside of the CVM-S4.26 basins:
Site | No GTL | 700m GTL | Smaller value |
---|---|---|---|
s117 | |||
USC | |||
WNGC | |||
PEDL | |||
SBSM | |||
SVD |
Proposed Algorithm
Our proposed algorithm for generating the velocity model is as follows:
- Set the surface mesh point to a depth of 25m.
- Query the CVM-S4.26.M01 model ('cvmsi' string in UCVM) for each grid point.
- Calculate the Ely taper at that point using 700m as the transition depth.
- Compare the values before and after the taper modification; at each grid point down to the transition depth, use the value from the method with the lower Vs value.
- Check values for Vp/Vs ratio, minimum Vs, Inf/NaNs, etc.
Cross-sections
Below are cross-sections of the TEST site velocity model at depths down to 600m for CVM-S4.26.M01, the taper to 700m, the lesser of the two values, and the % difference.
Depth | CVM-S4.26.M01 | Ely taper to 700m | Selecting smaller value | % difference, smaller value vs S4.26.M01 |
---|---|---|---|---|
0m (queried at 25m depth) | ||||
100m | ||||
200m | ||||
300m | ||||
400m | ||||
500m | ||||
600m |
Implementation details
To support the Ely taper approach where the smaller value is selected, the following algorithmic changes were made to the CyberShake mesh generation code. Since the UCVM C API doesn't support removing models, we couldn't just add the 'elygtl' model - it would then be included in every query, and we need to run queries with and without it to make the comparison. We also must include the Ely interpolator, since all that querying the elygtl model does is populates the gtl part of the properties with the Vs30 value, and it's the interpolator which takes this and the crustal model info and generates the taper.
- The decomposition is done by either X-parallel or Y-parallel stripes, which have a constant depth. Check the depth to see if it is shallower than the transition depth.
- If so, initialize the elygtl model using ucvm_elygtl_model_init(), if uninitialized. Set the id to UCVM_MAX_MODELS-1. There's no way to get model id info from UCVM, so we use UCVM_MAX_MODELS-1 as it is 29, we are very unlikely to load 29 models, and therefore unlikely to have an id conflict.
- Change the depth value of the points to query to the transition depth. This is because we need to know the crustal model Vs value at the transition depth, so that the Ely interpolator can match it.
- Query the crustal model again, with the modified depths.
- Modify the 'domain' parameter in the properties data structure to UCVM_DOMAIN_INTERP, to indicate that we will be using an interpolator.
- Change the depth value of the points to query back to the correct depth, so that the interpolator can work correctly on them.
- Query the elygtl using ucvm_elygtl_model_query() and the correct depth.
- Call the interpolator for each point independently, using the correct depth.
- For each point, choose to use either the original crustal data or the Ely taper data, depending on which has the lower Vs.
In the future we plan for this approach to be implemented in UCVM, and we can simplify the CyberShake query code.
High-frequency codes
For this study, we will use the Graves & Pitarka high frequency module (hb_high) from the Broadband Platform v22.4, hb_high_v6.0.5. We will use the following parameters. Parameters in bold have been changed for this study.
Parameter | Value |
---|---|
stress_average | 50 |
rayset | 2,1,2 |
siteamp | 1 |
nbu | 4 (not used) |
ifft | 0 (not used) |
flol | 0.02 (not used) |
fhil | 19.9 (not used) |
irand | Seed used for generating SRF |
tlen | Seismogram length, in sec |
dt | 0.01 |
fmax | 10 (not used) |
kappa | 0.04 |
qfexp | 0.6 |
mean_rvfac | 0.775 |
range_rvfac | 0.1 |
rvfac | Calculated using BBP hfsims_cfg.py code |
shal_rvfac | 0.6 |
deep_rvfac | 0.6 |
czero | 2 |
c_alpha | -99 |
sm | -1 |
vr | -1 |
vsmoho | 999.9 |
nlskip | -99 |
vpsig | 0 |
vshsig | 0 |
rhosig | 0 |
qssig | 0 |
icflag | 1 |
velname | -1 |
fa_sig1 | 0 |
fa_sig2 | 0 |
rvsig1 | 0.1 |
ipdur_model | 11 |
ispar_adjust | 1 |
targ_mag | -1 |
fault_area | -1 |
default_c0 | 57 |
default_c1 | 34 |
Spectral Content around 1 Hz
We investigated the spectral content of the Broadband CyberShake results in the 0.5-3 second range, to look for any discontinuities.
The plot below is from WNGC, Study 15.12 (run ID 4293).
Below is a plot of the hypocenters from the 706 rupture variations which meet the target.
These events have a different distribution than the rupture variations as a whole.
Fault | Percent of target RVs | Percent of all RVs |
---|---|---|
San Andreas | 60 | 44 |
Elsinore | 21 | 9 |
San Jacinto | 13 | 8 |
Other | 6 | 39 |
Additionally, 88% of the selected events have a magnitude greater than the average for their source. 4% are average, and 8% are lower.
Below is a spectral plot but which includes all rupture variations and the overall mean (in orange).
Additional Sites
We created spectral plots for 7 additional sites (STNI, SBSM, PAS, LGU, LBP, ALP, PLS), located here:
STNI | SBSM | PAS | LGU | LBP | ALP | PLS |
---|---|---|---|---|---|---|
Verification
Updates and Enhancements
Output Data Products
Computational and Data Estimates
Lessons Learned
Stress Test
Performance Metrics
Production Checklist
Science to-dos
- Perform validation for Northridge, Chino Hills, Whittier, Landers, Hector Mine, and North Palm Springs
- For each validation event, calculate BBP results, CS results, and CS results with BBP Vs30 values
- Check for spectral discontinuities around 1 Hz
Decide if we should stick with rvfrac=0.8 or allow it to vary
- Determine appropriate SGT, low-frequency, and high-frequency seismogram durations
- Update high-frequency Vs30 to use released Thompson values
Technical to-dos
Integrate refactoring of BBP codes into latest BBP releaseSwitch to using github repo version of CyberShake on SummitUpdate to latest UCVM (v22.7)- Switch to optimized version of rupture generator
- Test DirectSynth code with fixed memory leak from Frontera
- Switch to using Pegasus-supported interface to Globus transfers
- Test bundled glide-in jobs for SGT and DirectSynth jobs