CyberShake Source Filtering

From SCECpedia
Jump to navigationJump to search

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.

Fig. 1: Spectral content of a source for 0.5 hz (NT=20000, DT=0.01), filtered with a 4th-order low-pass Butterworth filter at 0.5 Hz
Fig 2: Spectral content of a source for 0.5 hz (NT=20000, DT=0.01), filtered with a 4th-order low-pass Butterworth filter at 1.0 Hz

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 Figure 2 source and compared it to a run with all the same parameters, except for using the Figure 1 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 wave propagation, 0.5 Hz filtered source
  2. 0.5 Hz wave propagation, 1 Hz filtered source
  3. 1 Hz wave propagation, 1 Hz filtered source
  4. 1 Hz wave propagation, 2 Hz filtered source

Results 1-3 are compared below; we are still running test #4.