MATLAB tools for modeling and plotting of shear-wave splitting data
- 01_stereoplots: Generate stereoplots of splitting results based on output of SplitLab and StackSplit.
- 02_data_handling: Read and visualize the shear-wave splitting data available from: https://doi.org/10.5445/IR/1000091427
- 03_modeling: Model shear-wave splitting results using synthetic splitting parameters computed for single-layer, two-layer and dipping layer scenarios.
If you make use of the content in this repository please acknowledge our paper published in GJI and/or my PhD thesis in whose framework several of the scripts and functions were developed:
-
Grund, M. & Ritter, J.R.R. (2020), Shear-wave splitting beneath Fennoscandia - Evidence for dipping structures and laterally varying multilayer anisotropy, Geophysical Journal International, 223, 1525–1547, https://doi.org/10.1093/gji/ggaa388
-
Grund, M. (2019), Exploring geodynamics at different depths with shear wave splitting, Dissertation, Karlsruhe Institute of Technology (KIT), https://doi.org/10.5445/IR/1000091425
Content of the directories:
Generate stereoplots of splitting results based on output of SplitLab and StackSplit.
- color-code bars of splits with respect to backazimuth (see below for supported colormaps)
- shade full background or specific backazimuthal wedges
- produce publication-ready plots in different formats (pdf, png, jpg...)
Supported colormaps
-
standard MATLAB colormaps: parula, winter, summer, copper....
-
MatPlotLib 2.0 Colormaps: Perceptually Uniform and Beautiful (see and download here: https://de.mathworks.com/matlabcentral/fileexchange/62729-matplotlib-2-0-colormaps-perceptually-uniform-and-beautiful)
-
Scientific colour-maps by F. Crameri (Zenodo. http://doi.org/10.5281/zenodo.1243862) (see and download here: http://www.fabiocrameri.ch/colourmaps.php)
Run function SWS_Analysis_BASICS_stereoplot
in the MATLAB command window to generate the stereoplots shown in the figure below using the provided test data set of station VAF.
Read and visualize the shear-wave splitting data available from: https://doi.org/10.5445/IR/1000091427
- generate individual stereoplots
SWS_read_evstruct('stereo')
- generate stereoplots of all 266 analyzed stations
SWS_read_evstruct('stereoall')
- generate histogram
SWS_read_evstruct('histo')
- generate output structs separated in split and null measurements
[SPLITS,NULLS]=SWS_read_evstruct
Model shear-wave splitting results using synthetic splitting parameters computed for single-layer, two-layer and dipping layer scenarios. Results plots are generated as well.
MSAT package of Walker & Wookey (2012) (https://github.com/andreww/MSAT) is required!
- reproduce the modeling results presented in our paper Grund & Ritter (2020, https://doi.org/10.1093/gji/ggaa388) using the shear-wave splitting data available from: https://doi.org/10.5445/IR/1000091427
- use your own measurements (SplitLab/Stacksplit output format, custom formatting is also supported)
Descriptions on how to set up different models and how to handle individual data sets can be found in the Supporting Information of Grund & Ritter (2020, https://doi.org/10.1093/gji/ggaa388).