Cascadia Simulations
From SCECpedia
Jump to navigationJump to searchDaniel Roten, Kim Olsen have used AWP-ODC-GPU software to run simulations of large magnitude earthquakes in the Cascadia region.
Cascadia M9.2
- Daniel Roten recently ran a simulation of a M9.2 Cascadia earthquake on Summit, to compare ground motions against the previous Titan simulations, using 375 V100 GPUs for 3:30 hrs.
- These data are on the Alpine file system in
/gpfs/alpine/proj-shared/geo112/rotend/Cascadia.alpine/scen_13_a_S_66m_dm_summit_flt
A previous version of the data is on Rhea/Titan in the directory:
/lustre/atlas/proj-shared/geo112/rotend/Cascadia/scen_8_a_N_66m_dm/output_sfc
Metadatda
Daniel Roten put together a list of metadata for NVIDIA to work with the Cascadia simulation output.
Points in X: 9,720 Points in Y: 16,065 Resolution along X and Y: 66.666 m The projection is UTM for zone 10. There is no rotation, with the X axis pointing from west to east and the Y axis from south to north. The UTM coordinates for the X and Y axis are given in the two attached ASCII files. The sampling interval in the output is 40 ms (250 samples per second). The total simulation duration is 450 seconds. The time step used in the finite difference simulation was 4 ms, and 102,000 time steps were computed. The output files are binary, and contain ground velocities (m/s) as 32-bit floating point numbers. There is no file header. Each component is saved in separate files, and the data are aggregated in files containing 100 time steps. The naming convention is S%c_0_%07d, where the character stands for the component (X, Y or Z), the 7-digit number for the last time step in the file (e.g., 1000 for the first file, and 102,000 for the last one). Inside each binary file, X increases fastest, followed by Y, and time t increasing slowest.
Animations:
Daniel Roten created a sample visualization done from a similar simulation (same resolution, but different scenario), which you might want to share along with the data:
