Files of the global stacks and vespagrams for users wanting to analyze them or generate plots.
This page contains supplementary material including very large .zip files containing all of the files needed should a user want to analyze the stacks or vespagrams or generate plots themselves. Also, some additional figures and animations are available here that likely would interest a more limited audience relative to the images on the main globalstacks page.
Global stacks files
Each zip file has 360 BHZ global stack files, one for each half degree wide distance bin. Zip files are organized by the narrow band-pass filter applied.
Columns within each file are:
- time from origin. 0-5400 s for BHZ, 0-10800 s for LH?
- length of STA window (s). The STA/LTA functions are delayed by this amount.
- number of contributing traces
- number of traces rejected by QCing
- raw stack of the filtered & QCed traces
- raw stack of the filtered & QCed traces using square-root stacking
- stack of STA/LTA functions
- stack of STA/LTA functions using square-root stacking (preferred for plotting)
- stack of envelope functions of the filtered & QCed traces
Make your own plots: a MATLAB starter script for simple plotting of the stacks PlotGlobalStacksBasic.m. Some published global body wave stacking papers are given below in the Citation section for reference.
Make your own stacks with your own data & processing parameters. This quick and dirty stacking exercise uses Fetch scripts and MATLAB to let users collect & process data and make stacks at a single distance. quick_n_dirty_globalstack.zip
Long period vertical, radial and transverse
|0.1 – 0.3 Hz||LHZ_0.1to0.3Hz (60 MB)||LHR_0.1to0.3Hz||LHT_0.1to0.3Hz|
|0.05 – 0.1 Hz||LHZ_0.05to0.1Hz||LHR_0.05to0.1Hz||LHT_0.05to0.1Hz|
|0.03 – 0.05 Hz||LHZ_0.03to0.05Hz||LHR_0.03to0.05Hz||LHT_0.03to0.05Hz|
|0.005 – 0.1 Hz||LHZ_0.005to0.1Hz||LHR_0.005to0.1Hz||LHT_0.005to0.1Hz|
|0.005 – 0.033 Hz||LHZ_0.005to0.033Hz||LHR_0.005to0.033Hz||LHT_0.005to0.033Hz|
Vespagrams (slowness stacks)
Vespagrams are an array analysis tool where traces from stations assumed to be along a common backazimuth and within a narrow distance range are stacked at different slownesses.
Animation of vespagrams made using USArray data in different frequency bands sweeping through different distances
The vespagrams were generated using all available BHZ data from the BK, CI, TA, US, and UW networks. Data were analyzed in marching 5 degree wide distance bins every 0.5 degrees. Within each distance bin and in different frequency bands, data were culled after calculating basic QC metrics, then each trace was shifted in time according to the average delay from a multi-channel cross-correlation analysis windowed on the first arriving P-wave. Then each vespagram was normalized by the peak amplitude and then further ‘median-normalized’ by iteratively raising the entire trace to the power 0.9 until the median value exceeded a threshold. This final step during event processing ensures a common noise level before stacking vespagrams from all events. Finally, vespagrams from each event within each frequency band are stacked. The dominant phases in the resulting stacks still have a very high signal-to-noise ratio, so a nonlinear color scale is used in the animation in order to highlight weaker phases.
Each zip file has vespagrams for 0.5 degree wide distance bins. The files are 743400 lines long, 126 slownesses (-10:0.2:15 s/deg), 5900 times (-500:1:5400 s around the origin time), and have four columns: linear stacking, square-root stacking, linear stacking using a higher threshold (higher background noise levels/lower SNR for peak arrivals), square-root stacking using a higher threshold.
Vespfiles_1_sec.zip (1.7 GB)
Comparison of results with those from Astiz et al. 1996.
Ray paths animation
Animation showing paths of different seismic phases using an earlier global stacks image.
To cite the IRIS DMC Data Products effort:
- Trabant, C., A. R. Hutko, M. Bahavar, R. Karstens, T. Ahern, and R. Aster (2012), Data Products at the IRIS DMC: Stepping Stones for Research and Other Applications, Seismological Research Letters, 83(5), 846–854, https://doi.org/10.1785/0220120032.
To cite the IRIS DMC Global Stacks data product:
- IRIS DMC (2014), Data Services Products: globalstacks Global stacks of millions of seismograms, https://doi.org/10.17611/DP/GS.1.
Astiz, L., P. Earle and P. Shearer, Global stacking of broadband seismograms, Seis. Res. Lett., 67, 8-18, 1996.
Paul S. Earle, Polarization of the Earth’s teleseismic wavefield, Geophys. J. Int. (1999) 139 (1): 1-8 “https://doi.org/10.1046/j.1365-246X.1999.00908.x”:http://gji.oxfordjournals.org/content/139/1/1.full
Rost, S., Thorne, S., and Garnero, E. Imaging global seismic phase arrivals by stacking array processed short-period data, 77(6),697-707, Seismological Research Letters, 2006. 10.1785/gssrl.77.6.697
Shearer, P.M., Imaging global body-wave phases by stacking long-period seismograms, J. Geophys. Res., 96, 20,353-20,364, 1991. doi: 10.1029/91JB00421
Shearer, P. M., Rychert, C. A. and Liu, Q. (2011), On the visibility of the inner-core shear wave phase PKJKP at long periods. Geophysical Journal International, 185: 1379–1383. doi: 10.1111/j.1365-246X.2011.05011.x