Skip to content

Latest commit

 

History

History
59 lines (36 loc) · 3.84 KB

README.md

File metadata and controls

59 lines (36 loc) · 3.84 KB
Raw

sws_tools

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:

01_stereoplots

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

  1. standard MATLAB colormaps: parula, winter, summer, copper....

  2. 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)

  3. 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.

PLOT_res_ALL

02_data_handling

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

PLOTS_merged

03_modeling

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!

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).

PLOTS_output