...
| Version | Num sources | SRF generation | VM domain | VM generation | Simulation specs | Recording stations | Estimated core-hours | Actual core hours | Notes |
|---|---|---|---|---|---|---|---|---|---|
| v17.8 | 15 (dominant in Canterbury) | 3 hypo and 2 slip dist per source | automated based on PGV>5cm/s; 15kmRup, 5km land cutoff Default depth and duration scaling | 0.4km regular grid, Vs_min=500 | Transition freq = 0.25Hz | 19,604 (virtual + Geonet stations) | ~3k Fitzroy | First implementation; Focus on running workflow and comparison with empirically-derived hazard curves | |
| v17.9 + v18p4 | ~ (South Island) | hypo every 20 km along strike, 3 slip dist per source | as v17.8 | as for v17.8 | as for v17.8 | 19,604 (virtual + Geonet stations) | -- hours on Kupe | Focus on extending number of sources and srf uncertianties | |
v18.5 | 251 North Island faults | Same as v17.9 | as for v17.9
| as for v17.9 | as for v17.9 | 19,604 (virtual + Geonet stations) | -- hours on Kupe |
| Running for the North Island sources to merge with 18p4 (and 17p9) to have a nation-wide hazard results |
v18p6
| 483 sources across NZ | Magnitude-dependant number of realizations one slip distribution per hypocentre | Pgv threshold of 2 cm/s | as for v17.9 | as for v17.9 | 27,481 (virtual + Geonet stations) | 150 k on Kupe |
| |
| v ?? |
Ideas for future implementations (no particular order):
...
- Bootstrap sampling to understand how many ruptures are needed for a given source
- Source uncertainties (currently slip and hypo; but need to add uncertainty in G&P parametrization).
- Velocity model uncertainties (random pertubations).
- Explicit modelling of subduction zone sources in Cybershake
- Neural Net for GMM trained with CS and validation results in order to use for distributed seismicity
Verification process:
Fault rupture (SRF) files:
- Source plot on the map
- Length-dependent hypocenter location (done)
- Different slip distribution realizations (done)
Velocity model (VM) files:
- Velocity model plot on the map (done)
- Rotation and reduction in size to minimize the ocean coverage (done)
- Generate a code to check VM files (done)
- Check if given VM parameter (folder path) exists
- .p, .s, .d and params_vel.py files must exist
- warns if model_params/coords/bounds etc don't exist using params_vel sufx in filename
- params_vel matching hh/xlen vs nx for x, y, z
- file size = nx * ny * nz * 4 bytes, checked for .p, .s, .d
- if numpy available: checks first xz slice for >0 and not NaN in .p, .s, .d
...