CyberShake Source Filtering

As of February 2015, CyberShake runs show a damping effect at spectral acceleration periods close to the simulation period. If we perform a 0.5 Hz CyberShake run, we see a damping factor of 2 or more when we look at the 0.5 Hz/2 second hazard curves.

We are interested in examining the filtering we are using to see if we can improve this so that hazard curves at the simulation period do not show damping. This would be a big benefit, as a 1 Hz CyberShake run could produce undamped 1 second hazard curves, which are of broad interest.

Background

CyberShake only performs filtering when the source impulse is created for the SGTs. Neither the SGTs nor the seismograms are directly filtered. The source is produced as part of the Graves SGT code; since it is location-independent, we calculate it once for a given DT and NT and use it for all corresponding CyberShake runs. We decided to run some experiments by changing the filtering of the source. By default, the source is low-pass filtered with a 4th-order Butterworth filter at the frequency of the CyberShake run. For this experiment, we propose filtering the source at 2 times the frequency of the CyberShake run.

Spectral content of two sources for a 0.5 Hz simulation (NT=20000, DT=0.01), filtered with a 4th-order low-pass Butterworth filter at 0.5 Hz (blue) and 1.0 Hz (green).

Experiment 1

We calculated a hazard curve for USC, ERF 35 (1000m fault surfaces), Graves & Pitarka 2010, with the AWP-ODC-GPU code and CVM-S4.26 using the 1 Hz filtered source and compared it to a run with all the same parameters, except for using a 0.5 Hz filtered source.

2 second SA. 0.5 Hz filtered source results are in magenta, 1 Hz filtered source results in black

3 second SA. 0.5 Hz filtered source results are in magenta, 1 Hz filtered source results in black

5 second SA. 0.5 Hz filtered source results are in magenta, 1 Hz filtered source results in black

10 second SA. 0.5 Hz filtered source results are in magenta, 1 Hz filtered source results in black

Experiment 2

We decided to run another test to compare 4 conditions:

1. 0.5 Hz simulation, 0.5 Hz filtered source -
2. 0.5 Hz simulation, 1 Hz filtered source -
3. 1 Hz simulation, 1 Hz filtered source -
4. 1 Hz simulation, 2 Hz filtered source -

Results are compared below.

2 second SA. 0.5 Hz sim/0.5 Hz filtered in magenta, 0.5 Hz sim/1 Hz filtered source results in black, 1 Hz sim/1 Hz filtered in red

3 second SA. 0.5 Hz sim/0.5 Hz filtered in magenta, 0.5 Hz sim/1 Hz filtered source results in black, 1 Hz sim/1 Hz filtered in red

5 second SA. 0.5 Hz sim/0.5 Hz filtered in magenta, 0.5 Hz sim/1 Hz filtered source results in black, 1 Hz sim/1 Hz filtered in red

10 second SA. 0.5 Hz sim/0.5 Hz filtered in magenta, 0.5 Hz sim/1 Hz filtered source results in black, 1 Hz sim/1 Hz filtered in red

We also plotted seismograms for the seismogram with the largest 3 second SA value, seen below:

Here are timeseries for the three sources. The 0.5Hz sim/1Hz filtering and 1Hz sim/1Hz filtering sources are practically identical, which is why you can only see 1 on the plot.