FuncLab: MATLAB Receiver Function Analysis
Contents
1. 1 000000e 01 beginning time of seismogram E 1 000000e 02 end time of seismogram IFTYPE TIME SERIES FILE defines the type of file LEVEN TRUE defines the sampling of the seismogram as even this is a MUST DELTA 2 500000e 02 sample rate of the seismogram Updated February 29 2012 32 IDEP UNKNOWN type of dependent variable not necessary for our purposes although we assume it is a velocity seismogram DEPMIN 1 740033e 02 minimum amplitude in seismogram DEPMAX 1 473583e 02 maximum amplitude in seismogram DEPMEN 8 144172e 01 mean amplitude in seismogram OMARKER 721 03 number of points in the time series TIMARKER 0 P arrival time and marker of predicted P wave this needs to be set as your zero time for the seismograms KZDATE AUG 14 227 2008 origin date of earthquake KZTIME 00 18 41 200 origin time of earthquake IZTYPE BEGIN TIME defines the reference time for this file needs to be set to BEGIN TIME KSTNM GO8A station code CMPAZ 0 000000e 00 component of motion azimuth needs to be 0 for vertical BAZ 180 for radial and radial RFs and BAZ 90 for transverse and transverse RFs CMPINC 0 000000e 00 component of motion inclination needs to be 0 for vertical and 90 for radial or transverse including RFs STLA 4 529040e 01 station latitude STLO 1 189595e 02 station longitude STEL 1 318000e 03 station elevation in meters ST2. HkMetadataStrings and HkMetadataDoubles arrays to store metadata for each H xk stack that will be saved to the Project mat file Several visualizations to aid the user in the analysis are included in the package iene Updated February 29 2012 25 References Abers G A X X Hu and L R Sykes 1995 Source scaling of earthquakes in the Shumagin region Alaska time domain inversion of regional waveforms Geophys J Int 123 41 58 doi 10 1111 1365 1246X 1995 tb06660 x Bostock M G 1996 Ps conversions from the upper mantle transition zone beneath the Canadian landmass J Geophys Res 101 8393 8402 Clayton R W and R A Wiggins 1976 Source shape estimation and deconvolution of teleseismic body waves Geophys J R Astron Soc 47 151 177 Dziewonski A M and D L Anderson 1981 Preliminary reference Earth model Phys Earth Planet Int 25 297 356 Gurrola H G E Baker and J B Minster 1995 Simultaneous time domain deconvolution with application to the computation of receiver functions Geophys J Int 120 537 543 doi 10 1111 j 1365 1246X 1995 tb01837 x Eagar K C and M J Fouch 2012 FuncLab A MATLAB interactive toolbox for handling receiver function datasets Seis Res Lett in press Eagar K C M J Fouch and D E James 2010 Receiver function imaging of upper mantle complexity beneath the Pacific Northwest United States Earth Planet Sci Lett 297 141
3. e hh Aan hi HVN ihe 55 11 faba Aig 54 06 aR i o en is KVN j 53 Wall 51 95 3 50 89 W 49 84 hood a vind in BAN RARE vith 3 seule 48 78 aA iy i 4 if a i f a ti ly its hag OM Mi AN i WN fi WIA Py oe 0 20 40 60 80 0 20 40 60 80 0 20 40 60 80 0 20 40 60 80 Time sec Time sec Time sec Time sec Time sec w ii in i y Y Updated February 29 2012 16 090 lt Student Version gt TA 107A Stacked Records File Edit View Insert Tools Desktop Window Help TA I07A 1 he A afl 0 75 0 5 0 25 0 0 25 Vert 0 5 0 75 1 Rad Trans Rad RFs Trans RFs Time on The RF Moveout Record Section item plots receiver function wiggle traces of active records vs ray parameter and the RF Moveout Image plot the vertical radial and transverse seismograms and radial and transverse receiver functions The RF Moveout Record Section item plots receiver function wiggle traces of active records vs ray parameter and the RF Moveout Image stacks records into ray parameter bins and displays the amplitudes by a color scale eoo lt Student Version gt TA I07A Receiver Functions eoe lt Student Version gt TA 107A Receiver Functions File Edit View Insert Tools Desktop Window Help File Edit View Insert Tools Desktop Window Help a DSes kh BAVO8H4 a 08 oO gt DGH kh 828 93D84 8208 ao 2 S
4. 0 1 0 08 0 06 70 02 gg 0 02 0 04 0 06 0 08 Amplitude The Dataset Statistics item is not a true plot or visualization but displays a textbox summarizing various statistics such as the number of total active records about the dataset Another data coverage visualization is the Ray Parameter Backazimuth Diagram which plots these metadata on a polar plot Finally two map options a Station Map and a Global Event Map are available if the MATLAB Mapping Toolbox is installed eoo lt Student Version gt Data Coverage File Edit View Insert Tools a OGaSsS k Desktop Window Help qS ODEA a 16 ao Updated February 29 2012 18 5 Project Directory Structure and Description A FuncLab project defines an entire dataset to be analyzed It can contain multiple data directories as defined in section 3 A project directory contains all the raw data used by FuncLab as well as files and directories it creates when a user initiates a new project A new project is created within the FuncLab GUI by importing at least one data directory of the structure defined in section 3 The following shows the structure of the automatically generated project directory Updated February 29 2012 19 FUNCLAB_ PROJECT top directory containing project files no naming restriction Datastats Tablemenu RAWDATA RecordMetadataStrings RFDATA directory containing imported data as describ
5. 1s located at http geophysics asu edu funclab Download the latest release and untar it in the directory in which you would like to place the MATLAB M files For example tar xvf funclabl 5 3 tar The directory funclab1 5 3 contains the files needed to run all the functions within the FuncLab toolbox Appendix A contains a list of all the M files Note that data for mapping functions requires the Mapping Toolbox may be downloaded from outside data sources Links to these data and other relevant websites are provided in Appendix B of this manual and on the FuncLab website These data include DEM or topography such as ETOPOI ETOPO2v2c and GLOBE DEM as well as GSHHS coastline data I recommend placing these under funclab1 5 3 map_data In order to use the FuncLab functions the user will need to add the paths to the files in the MATLAB session To perform this task the user must change to the main funclab directory and then run the setup funclab script cd funclab1l 5 3 setup funclab The setup funclab script automatically adds all the directories to the search paths in MATLAB In addition this script checks for the Mapping and Signal Processing Toolboxes on the machine If these are installed certain functions will be available for use However if they are not installed some functions may not be available If the user has administrative rights on the machine running FuncLab they may choose to insert the full pathname of the set
6. 5 500000e 00 earthquake body wave magnitude One note of caution the CMPINC specifications in the description assume that the seismograms have been rotated in the direction of the earthquake along the plane of the surface Rotations into true radial and transverse along the back projected ray path will have a different inclination The iterative deconvolution code computes a receiver function by first performing a cross correlation on the vertical and radial seismogram Then based on the lag time convolves a spike train that represents that lag time with a Gaussian function The Gaussian parameter determines the width of the Gaussian wavelet and essentially acts as a low pass filter for the receiver function After the first pulse is generated on the predicted receiver function the receiver function is convolved with the vertical seismogram to get a predicted radial seismogram The difference between the predicted Updated February 29 2012 34 radial and the data is then computed to get a residual As long as the residual is greater than a given threshold the process continues The first pulse 1s stripped and the next pulse is generated again by the process of cross correlation convolution with a Gaussian wavelet to produce a predicted receiver function and convolution of the receiver function with the vertical seismogram More details can be found in Ligorria and Ammon 1999 Updated February 29 2012 35
7. TOMOMODELS directory in the MAT file NWUS_P09b_1 2 mat This file stores 6 variables used to define the model Depths Lons Lats ModelP ModelS and TomoType The first 5 are 3x3 matrices where the first second and third dimensions define x y and z longitudes latitudes and depth of the 3D volume The ModelP and ModelS variables define the P and S wave perturbations or velocities respectively at each point in the volume TomoType either defines the model as a perturbation model or a velocity model to handle differencing the volume from the 1D velocity model used in the ray tracing As with velocity models other tomography models can be built to be used in FuncLab These new models should be place under the funclab1 5 3 TOMOMODELS and have the extension mat Updated February 29 2012 D2 7 Add on Packages These are special code packages that are optional with FuncLab running as the base For instance common conversion point CCP stacking of receiver functions to look at mantle discontinuities is a common approach Similarly H stacking of receiver functions for examining Moho depth and crustal Vp Vs is another common analysis FuncLab includes add on packages for both of these analyses Although these add on packages are integrated with FuncLab it is also separate from the FuncLab main code This means that a core add on package does NOT have to be altered to change things associated with the CCP package Here w
8. The Logfile txt 1s a text file that is used to store a log of all the processes run within FuncLab on a specific project It contains a description of the process parameters used and time stamps This Updated February 29 2012 NO 0 is generated automatically and can be useful when retrieving a record of the different runs for each process The Project mat file is a MATLAB formatted data file that contains all the MATLAB variables used to reference the raw seismic data and metadata within a project It 1s essentially the key database of the project There are also variables contained in Project mat that are used to communicate with the different GUIs of FuncLab and keep them talking to one another Other directories will be added under RFPROJECT once certain add on packages are initiated with a project These new directories are CCPDATA GCCPDATA and HKDATA Each of these contains MAT files storing variables created in MATLAB specific for the types of analyses performed For example the CCPDATA directory contains CCPData mat which stores variables that define the CCP grid and ray tracing information for each receiver function used in the stacks It also contains a separate MAT file for each CCP bin e g Bin 1 000360 _040080 mat that stores variables defining the CCP stack for that bin This is similar for GCCPDATA The HKDATA directory contains a separate MAT file for each H x stack e g TA FIOA_ 000360 _040080 mat The
9. m deg2km m deg2rad m gg2gc m km2deg m latlon_ from m rad2deg m utilities colvector m ignoreNaN m makeColorMap m mmpolar m skm2srad m srad2skm m Updated February 29 2012 30 APPENDIX B WEBSITE REFERENCES Seismic Analysis Code SAC http www iris edu software sac manual html Standing Order for Data SOD http www seis sc edu sod Receiver Function Codes Charles Ammon Penn State http eqseis geosc psu edu cammon HTML RftnDocs thecodes01 html Generic Mapping Tools GMT http Awww soest hawai edu gmt Earthscope Automated Receiver Survey EARS http ears iris washington edu MATLAB Mapping Data e Coastlines O gshhs 2 1 1 http www ngdc noaa gov mg g shorelines data gshhs version2 1 1 gshhs _ 2 1 1 zip e Topography O l arc minute DEM ETOPO1 ice c 12 Ice Surface http www ngdc noaa gov mg global relief ETOPO1 data ice_surface ce IL_registered binary etopol_ ice _c_12 zip ETOPO1 bed c 12 Bedrock Surface http www ngdc noaa gov mg global relief ETOPO 1 data bedrock cell_r egistered binary etopol bed_c_12 zip 2 arc minute DEM ETOPO2v2c 12 MSB http www ngdc noaa gov mg global relief ETOPO2 ETOPO2v2 2006 ETOPO2v2c raw_binary ETOPO2v2c 12 MSB zip 30 arc second 1 km grid DEM GLOBE DEM http www ngdc noaa gov mgg topo DATATILES elev all10g zip Updated February 29 2012 31 APPENDIX C DESCRIPTION O
10. side 1 D ray tracing creates the RayMatrix XX and MidPoints XX matrices that store the ray coordinates at each depth increment and half depth increment respectively The XX is included in the matrix name to account for the possibility that more than one matrix may need to be generated since each matrix is allowed to store up to 25 000 records to save memory in MATLAB Both variables are stored in the CCPData mat MAT file under the created subdirectory in the project directory CCP The correction for 3 D heterogeneity is an optional step and requires the user to choose a tomography model which can be defined or imported by the user as explained in the user manual The MidPoint XX matrix will be populated with timing corrections as described in Eagar et al 2010 and are stored in the TimeCorrections XX matrix also saved to CCPData mat Time to depth conversion saves receiver function amplitudes in RayMatrix XX that are now associated with each ray coordinate and depth increment Each CCP bin saved in the CCPGrid matrix 1s assigned the amplitudes associated with it at each depth increment CCP stacking results in the creation of two new arrays to store the CCP metadata CCPMetadataStrings and CCPMetadataDoubles which are saved to the Project mat file The stacking function also adds a row to the Tablemenu variable in the Project mat file to control the drop down menu in the main FuncLab GUI and allow the user to view records of CCP bins in the same m
11. software sac manual html It can read and handle SAC files that are either little endian or big endian byte order automatically It also requires a specific directory structure for the seismic data import and the creation of a new FuncLab project next section This directory structure originated from the output of publicly available data download using the Standing Order for Data SOD program http www seis sc edu sod The choice of pre processing methods to obtain receiver functions 1s independent of FuncLab and is left up to the user Our choice of deconvolution schemes is an iterative deconvolution receiver function code written by Charles Ammon at Penn State University http eqseis geosc psu edu cammon HTML RftnDocs rftn01 html FuncLab is thus written assuming that the receiver function output will be similar to that for this iterative deconvolution package The link to download this code 1s also posted on the FuncLab website Also posted are some generic C shell scripts used to run the preprocessing of data in this directory structure and output filenames appropriately to help ensure that receiver functions are produced in a way in which FuncLab can read them We also note that the EarthScope Automated Receiver Survey EARS data product http ears iris washington edu provides receiver functions that can be used in FuncLab using the script EARS2FL csh also included in this release For a more complete description of SAC file
12. that S to P receiver function analysis could be incorporated via a community developed add on These tools were developed on the Mac OS X platform beginning in 2006 with MATLAB 7 1 0 21 R14 student version Most of the developed tools were created to be compatible with this version but were not tested to be backwards compatible with earlier versions of MATLAB Figure dimensions may also vary from platform to platform but reflect development on Red Hat Linux 5 3 running MATLAB 7 9 0 R2009b Testing has also been conducted on Mac OS X 10 6 7 running MATLAB 7 10 0 R2010a without problems There are some tools that are used within the toolkit that were not originally developed with FuncLab but are borrowed from free open source forums on the web They are l ignoreNaN m used to handle NaN values within matrices 2 makeColorMap m used to adjust color maps in figures 3 mmpolar m creates a polar plot of ray parameter vs backazimuth Updated February 29 2012 3 2 Download and Setup This section introduces the step by step instructions for downloading and setting up the M files in the FuncLab package for your MATLAB session For the purposes of this manual UNIX and MATLAB command line inputs are typed in courier font file names and MATLAB variables are in italics and directory names are bolded The FuncLab website which includes the source code a quick start guide this manual sample datasets and links to needed external files
13. 153 Kennett B L N and E R Engdahl 1991 Traveltimes for global earthquake location and phase identification Geophys J Int 105 429 465 Kennett B L N Engdahl E R and R Buland 1995 Constraints on seismic velocities in the Earth from travel times Geophys J Int 122 108 124 Langston C A 1979 Structure under Mount Rainier Washington inferred from teleseismic body waves J Geophys Res 84 4749 4762 Ligorria J P and C J Ammon 1999 Iterative deconvolution and receiver function estimation Bull Seism Soc Am 89 1395 1400 Park J and V Levin 2000 Receiver functions from multiple taper spectral correlation estimates Bull Seism Soc Am 90 1507 1520 Sheehan A F G A Abers C H Jones and A L Lerner Lam 1995 Crustal thickness variations across the Colorado Rocky Mountains from teleseismic receiver functions J Geophys Res 100 20 391 20 404 Vinnik L P 1977 Detection of waves converted from P to SV in the mantle Phys Earth Planet Int 15 39 45 Updated February 29 2012 26 West J D M J Fouch J B Roth and L T Elkins Tanton 2009 Vertical Mantle Flow Associated with a Lithospheric Drip Beneath the Great Basin Nature Geoscience 2 439 444 Updated February 29 2012 27 APPENDIX A LIST OF M FILES setup funclab m funclab m addons ccp sccp hk ccp_ 1D_raytracing m ccp assign m ccp binmap m ccp_ binstack m ccp_
14. CC STD Station Map CC SWFL Global Event Map CC VALT Ray Parameter Backazimuth Diagram j The Seismograms Record Section and Stacked Records items plot the vertical radial and transverse seismograms and radial and transverse receiver functions Note that the record section plot is more effective when switching to Event Tables and highlighting Updated February 29 2012 15 an event rather than a particular station A quick method of QC and trace editing your data maybe to look at these plots by event then by right clicking the event table you can select either Turn On or Turn Off to edit all records in that event eoo File Edit View Insert Tools lt Student Version gt Station TA I07A Event 2007 197 01 13 22 Seismograms Desktop Window Help Station TA I07A Event 2007 197 01 13 22 2000 4 Vertical o 2000 4 4000 1000 4 Radial 1000 4 1000 4 500 4 Transverse 0 500 4 1000 4 Radial RF Transverse RF Time s eoo lt Student Version gt 2007 076 22 43 09 Record Section File Edit View Insert Tools Desktop Window Help Record Section 2007 076 22 43 09 Vertical Radial Transverse Radial RFs Transverse RFs k kii hi 1D X ah iv Ih Ki yh n y ps al HUN has i M TA mu i mo f j f iy Pi i Hn i W ly i ih nil M ae iy ait an EASY Hy A a inch ahi ht ik h ie mea a 57 22 my Nis HIM Me ir vag Koy ini bi ei Me Ni ek Ae i 56 17 PA A Ai in
15. DP 0 000000e 00 station depth in meters EVLA 1 640000e 01 earthquake latitude EVLO 1 469000e 02 earthquake longitude EVEL 0 000000e 00 earthquake elevation in meters EVDP 5 300000e 04 earthquake depth in meters DIST 9 048033e 03 distance in kilometers from station to earthquake Updated February 29 2012 33 AZ 4 540545e 01 azimuth from earthquake to station BAZ 2 844766e 02 backazimuth from station to earthquake GCARC 8 138271e 01 distance in degrees from station to earthquake LOVROK TRUE set to TRUE to overwrite the file USERO 2 500000e 00 Gaussian parameter used to compute iterative deconvolution RF see description of process below USERI 3 028038e 02 predicted P wave ray parameter in sec rad USERS 0 000000e 00 predefines the trace edit status this is OPTIONAL but set to 0 if it is bad or I if it is good and you want to define it as an active RF USER9Y 6 000000e 01 fit or match of the predicted waveform from iterative deconvolution see description of process below NVHDR 6 header version number current version is 6 LPSPOL TRUE set to TRUE if component has positive polarity LCALDA TRUE set to TRUE if DIST AZ BAZ and GCARC are calculated from station earthquake coordinates KCMPNM BHZ component code should be set to BHZ for vertical BHR for radial and radial RFs and BHT for transverse and transverse RFs KNETWK TA network code MAG
16. F SAC FILES SAC files for the receiver functions and seismograms may be produced in a variety of ways but regardless of the method of preprocessing some important formatting standards must be followed This appendix describes the format and headers needed in each of the files FuncLab imports 5 separate SAC files for a given calculated receiver function These include the vertical seismogram radial seismogram transverse seismogram radial receiver function and transverse receiver function Itis preferred that all 5 files are time synched to the zero time of the receiver functions This makes plots in FuncLab make a lot more sense and prevents issues when interpreting the time axis Naming of the SAC files is an important first step All SAC files have a naming convention as follows STA 2 5 1 2 where STA is the station name 2 5 represents the Gaussian lowpass filter or any type of lowpass filter used 1 means that we used iterative deconvolution to calculate the receiver function and z is the seismogram component of motion The station name and the extension are the most important items of the name The extensions define the type of file as follows Z vertical seismogram r radial seismogram t transverse seismogram eqr radial receiver function eqt transverse receiver function The SAC header definitions and descriptions for the seismograms are below NPTS 4401 number of points in seismogram
17. FuncLab MATLAB Receiver Function Analysis Toolbox User Manual Version 1 5 3 Author Kevin C Eagar Email funclab dtm ciw edu Contents 1 ABOUT FUNCLAB 2 DOWNLOAD AND SETUP 3 RECEIVER FUNCTION OVERVIEW PREPARATION AND DATA DIRECTORY STRUCTURE 3 1 RECEIVER FUNCTION OVERVIEW 3 2 RECEIVER FUNCTION FORMAT AND PREPARATION 3 3 RECEIVER FUNCTION RAW DATA DIRECTORY STRUCTURE 4 GRAPHICAL USER INTERFACE 4 1 STARTING A NEW PROJECT 4 2 LOADING AN EXISTING PROJECT 4 3 PROJECT MAT FILE AND FUNCLAB VARIABLES 4 4 MAIN FUNCLAB GUI 4 5 TRACE EDITING 4 6 VIEWING DATA PLOTS AND VISUALIZATIONS 5 PROJECT DIRECTORY STRUCTURE AND DESCRIPTION 6 MODEL BUILDING AND SETUP 6 1 VELOCITY MODELS 6 2 TOMOGRAPHY MODELS 7 ADD ON PACKAGES APPENDIX A LIST OF M FILES APPENDIX B WEBSITE REFERENCES APPENDIX C DESCRIPTION OF SAC FILES Updated February 29 2012 2 1 About FuncLab FuncLab comprises a set of tools built within the MATLAB environment to analyze receiver functions Eagar and Fouch 2012 Receiver functions are seismic time series computed by deconvolving the radial from the vertical seismogram to isolate the Earth s radial S wave response While P to S conversions arrive later relative to the P wave S to P conversions atrive earlier relative to the S wave Both mode conversions can be utilized in receiver functions for different purposes In this first release of FuncLab tools are available for P to S receiver function analysis We note
18. FuncLab at this point will create the new project directory with some sub directories copy the data directories into the project directory and populate the new Project mat file that will handle most of the metadata for interaction with the GUI eoo lt Student Version gt FuncLab 1 5 2 File Editing View Export Add ons New FuncLab Project Setup Please type the full pathname for your NEW project top directory NOTE This directory must not already exist Users keagar research HLP_RFs Please select up to 10 datasets to add to this new project Browse Users keagar research HLP_RFs data_0 6 TA_ALL_07 Browse Browse Browse Browse Browse Browse Browse Browse Browse 4 2 Loading an Existing Project If you are working from an existing FuncLab project you can load the data into FuncLab by select File Load Project from the top menu Updated February 29 2012 eo Ji Editing View Export Add ons OO Select Your Project File About FuncLab Filter keagar research HLP_RFs New_PNW_Mantle maf Directories s New_PNW_Mantle G s New_PNW_Mantle s New_PNW_Mantle CCPDATA New Project Load Project FuncLab Preferences Exit FuncLab s New_PNW_Mantle GCCPDATA s New_PNW_Mantle RAWDATA Files of type MAT files mat Selection rresearch HLP_RFs New_PNW_Mantle Project mal Open Filter Cancel Select the Pro
19. XC_ALL_07_a to Users keagar research HLP_RFs Ne Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA XC_ALL_07 Moved Data Directory Users keagar research HLP_RFs data_0 6 XC_ALL_08 to Users keagar research HLP_RFs New_ Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA XC_ALL_0E 04 30 2011 18 55 50 Saved New Trace Editing Preferences 04 30 2011 19 30 33 Saved New FuncLab Preferences 04 30 2011 19 32 04 Saved New FuncLab Preferences 05 01 2011 06 26 22 Manual Trace Editing TA BO5A 05 01 2011 06 31 44 Manual Trace Editing XC ID001 05 01 2011 06 35 11 Manual Trace Editing XC ID002 05 01 2011 06 36 03 Saved New Trace Editing Preferences N5 N7 9N11 17 57 11 Manual Trace Faitinn XC OR1NG Pe 4 4 Trace Editing From the main FuncLab GUI the Trace Editor can be accessed through the Editing menu which contains options for manual and automated trace editing The first step that should be performed is setting the trace editing preferences These can be accessed through Editing Trace Editing Preferences These control aspects of the GUI display such as receiver function component radial or transverse time axis parameters gain number of receiver functions in the view and sample rate for resampling Other preferences such as the maximum amplitude and the minimum receiver function fit can be specified as control parameters in a
20. amp a fa a ES Radial RF Transverse RF Trace Count N A oo 1 N o 2 2 o 10 P35s ilo 0 06 im gt 7 10 F40 04 oe Au Git at Omer apt Ng A sh i PRC 40 02 S k a rs T H Amplitude 30 T D 3 a G z E 40 E r 7 0 02 5 50 C Pacs an m lss P4 0 04 l 0 06 70 i Lao 0 08 al 90 0 1 PEBCe i 0 08 Ray Parameter sec km Ray Parameter sec km 0 04 0 05 0 06 0 07 0 08 0 04 0 05 0 06 0 07 0 08 Ray Parameter sec km Ray Parameter sec km RF Backazimuth Record Section and RF Backazimuth Image are similar except that records are displayed vs backazimuth Updated February 29 2012 17 eoo lt Student Version gt TA 107A Receiver Functions Time sec a USH R 90 Radial RF 135 180 a 225 Event Backazimuth deg 270 315 File Edit View Insert Tools Desktop Window Help ALD eR 0a ag Transverse RF 360 ie 45 90 135 180 225 270 315 360 Event Backazimuth deg eoo lt Student Version gt TA 107A Receiver Functions File Edit View Insert Tools Desktop Window Help gt OGaS k AYDER a O88 Oo Time sec 0 07 Ray Parameter sec km Ray Parameter sec km Trace Count gt N A oo 6 1 WS o o zs o Em de F F30 p40 p50 p60 F70 F80 30 0 08
21. anner as receiver function records An lt Student Version gt CCP Stacking Export Plotting Analysis CCP Parameters CCP Processes Bootstrap Resamples o99c Bootstrap Percentage Mise 1D Ray Tracing Latitude Max M a Latitude Min Tac 3D Heterogeneity Corrections Longitude Max 41 Longitude Min 12 Depth Max km aoc Depth Inc km Time to Depth Migration Bin Width km 5c Bin Spacing km 5C CCP Assignment Backazimuth Min c Backazimuth Max EE Ray Parameter Min s km o 04 Ray Parameter Max s km 0 08 Vertical Overlap bins Min Hits Per Bin ad COP Sucking 1D Velocity Model AK135 vel Hd Tomography Model NWUS_P09b_1 2 mat The GCCP GUI is quite similar to the CCP GUI with only a few minor differences such as the choice of the Gaussian width parameter and weights which are Updated February 29 2012 24 stored in the RayMatrix XX variable In all other variables and functions CCP is replaced with GCCP M files for this add on are found in the subdirectory addons gccp Algorithms for H x stacking are included in m files found in the subdirectory addons hk To perform a single station stack the user selects the H x Stack Single Station button All active records from the station selected from the station tables in the main FuncLab GUI will be stacked based on the parameters selected on the left hand side of the H Stacking Analysis GUI Stacking one or all stations will create the
22. bootstrap m ccp Callback m ccp_corrections m ccp ex pp m ccp ex stacks m ccp fieldnames m ccp migration m ccp readamps m ccp record fields m ccp setup grid m ccp _stacking m ccp m init_ccp m gccp 1D raytracing m gccp assign m gccp binmap m gccp_ binstack m gccp_bootstrap m gccp Callback m gccp corrections m gccp ex pp m gccp ex volume m gccp_ fieldnames m gccp migration m gccp readamps m gccp record fields m gccp setup grid m gccp_ stacking m gccp m init_gccp m Updated February 29 2012 28 hk allstack m hk Callback m hk ex_grids m hk ex_solutions m hk record fields m hk _singlestack m hk view baz image m hk view _hkgrid m hk view _moveout_image m hk view_stacked_records m hk m init hk m custom flab errordlg kce m movegul kce m questdlg kce m uigetdir_kce m uigetfile kce m fl add data gui m fl add new data m fl_datainfo m fl_datastats gui m fl_default_prefs m fl_edit logfile m fl init _logfile m fl log gui m fl pref fieldnames m fl preferences gui m fl readseismograms m fl record fields m fl_traceedit gui m fl view baz image m fl view baz m fl view data coverage m fl view event map m fl view_moveout_image m fl view_moveout m fl view record section m fl view _record m fl view stacked _records m fl view station map m te default _prefs m te pref fieldnames m te preferences gui m te view record m Updated February 29 2012 29 map_utilites azimuth
23. e describe how to develop separate add on packages for additional receiver function processing There is a specific directory structure needed for an add on as well as a specific location relative to the FuncLab directories In funclab m the path of the funclab command is searched for That path should be the same path that the flab directory is housed This is the directory that contains the entire core FuncLab code Place the add on directory in the addons subdirectory For example the cep directory contains all the codes for the CCP add on All add on code must be in a single subdirectory FuncLab will not recognize any subdirectories beneath the main one such as addons gt ccp gt working As a point of organization FuncLab includes m files with the prefix of the add on directory e g ccp_assign m The next component is an m file under this add on directory called init_ add on m Replace add on with the name of the add on In the CCP add on case it would be called init_ccp m This file need only contain 1 line that looks like this uitmenuCMAddons Label Common Conversion Point Stacking Tag ccp_m Enable off Callback ccp_CalLlback The string after Label is the name to be listed under the Add ons toolbar in the main FuncLab GUI Tag should just be a unique name In FuncLab the add on name is a prefix and the letter m which stands for main Enable should ALWAYS be off E
24. ed 04 30 2011 09 15 31 Last Updated 06 16 2011 21 49 11 04 30 2011 09 15 31 Initiated FuncLab Project 1 Copied Data Directory Users keagar research HLP_RFs data_0 6 BK_ALL_04 to Users keagar research HLP_RFs New Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA BK_ALL_04 Copied Data Directory Users keagar research HLP_RFs data_0 6 BK_ALL_05 to Users keagar research HLP_RFs New Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA BK_ALL_0O amp Copied Data Directory Users keagar research HLP_RFs data_0 6 BK_ALL_06 to Users keagar research HLP_RFs New Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA BK_ALL_0O Copied Data Directory Users keagar research HLP_RFs data_0 6 BK_ALL_07 to Users keagar research HLP_RFs New Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA BK_ALL_07 Copied Data Directory Users keagar research HLP_RFs data_0 6 BK_ALL_07_a to Users keagar research HLP_RFs Ne 0 Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA BK_ALL_C Copied Data Directory Users keagar research HLP_RFs data_0 6 BK_ALL_08 to Users keagar research HLP_RFs Nev Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA BK_ALL_C Copied Data Directory User
25. ed in Figure 4 2 RecordMetadataDoubles FuncLabPreferences Project mat TraceEditPreferences project file containing metadata and FuncLab organization variables CCPMetadataDoubles CCPMetadataStrings GCCPMetadataDoubles GCCPMetadataStrings Logfile txt text file containing log of processes performed on project CCPDATA HkMetadataDoubles directory containing files created HkMetadataStrings using CCP Stacking GUI Add on CCPGrid CCPData mat n MAT file for CCP data created after 1 D Ray Tracing DepthAxis MidPoints XX Bin_1_ 000360 040080 mat RayMatrix XX CCP bin containing stacked data Bin_X_BMNBMX_RMNRMX mat Bin_ BootAxis subsequent CCP bins RRFBootMean RRFBootSigma GCCPDATA directory containing files created using GCCP Stacking GUI Add on HKDATA directory containing files created using H k Stacking GUI Add on TRFBootSigma 3 D 2 3 3 5 3 x TA_J08A_000360_ 040080 mat H k station analysis containing stacked data NET_STA_BMNBMX_RMNRMX mat HkEventList TA_KO9A kAxis can have multiple stations RFPROJECT this is a generic name given for the purposes of illustration is the top directory name which the user specifies at the onset of a new project The RAWDATA directory contains all the data directories that a user imports into the project Notice that these are of the same structure as described in the previous section
26. elected or the entire dataset Editing Auto Trace Edit ALL The second criteria is a minimum fit or measure of uncertainty in the receiver function expected to be set in the USER9 SAC header 8 0 0 8600 600 Fie eel View Export Add ons File E View Export Add ons File EA View Export Add ons Manual Trace Edit S Auto Trace Edit Selected TA Auto Trace Edit ALL TA Trace Edit Import TA Manual Trace Edit S Auto Trace Edit Selected TAST OEC N TA Trace Edit Import TA Manual Trace Edit S Auto Trace Edit Selected TA Auto Trace Edit ALL TA TA TA Trace Editing Preferences TAa be a a S TA Trace Editing Preferences TA Trace Editing Preferences 4 5 Viewing Data Plots and Visualizations A number of data plotting functions are located under the View menu in the main FuncLab GUI and are divided into three sections The top section displays plots related to the selected record in the center panel Options in the middle section display plots of active records in the selected table in the left hand panel The bottom section includes visualizations that involve active and inactive records within the dataset eoo File Editing Export Add ons Seismograms Station Tat Record Section Stacked Records BK HATC RF Moveout Record Section RF Moveout Image RF Backazimuth Record Section RF Backazimuth Image Dataset Statistics
27. ff are shaded gray eoc File EIEN View Export Add ons 68090 lt Student Version gt Manual Trace Edit S Auto Trace Edit Selected EK Auto Trace Edit ALL E BK Trace Edit Import 10 Receiver Function Trace Editing GUI Save Edits Cancel TA IO7A OOMOOODD BK Trace Editing Preferences DIF Lonanana Time sec 30 Sg ee ace anaes aera ey ean radar A za cares p 9 5 10 15 20 25 30 35 Trace Number Furthermore right clicking the waveform brings up a menu to start another GUI for visualizing the seismograms and radial receiver function change the record status and displays event and record information useful in determining the appropriate status of the record Figure 7 The record status is changed in the Status column of RecordMetadataDoubles and saved to Project mat when the Save Edits button 1s selected Updated February 29 2012 14 eoo lt Student Version gt Station TA I07A Event 2006 236 21 50 36 Seismograms File Edit View Insert Tools Desktop Window Help M Status Close Mag 5 9 mb Dist 53 59 deg Depth 43 00 km 2000 4 1000 4 Vertical 0 i ae 1000 5 2000 4 1000 4 500 4 Radial 0 500 4 1000 4 500 4 Transverse o 500 4 Radial RF 0 55 TT T T T T T T T T 1 10 0 10 20 30 40 50 60 70 80 90 Time s Automated trace editing can be performed on a selected table Editing Auto Trace Edit S
28. first and second set of numbers in the bin and H stack files refers to the backazimuthal range and ray parameter range in the stacking respectively Updated February 29 2012 21 6 Model Building and Setup 6 1 Velocity Models FuncLab is currently only designed to handle 1D Earth models in a very simple way This version of the package contains three commonly used 1D models PREM Dziewonski and Anderson 1981 ASP91 Kennett and Engdahl 1991 and AK135 Kennet et al 1995 These are contained under the funclab1 5 1 VELMODELS directory within the text files PREM vel IASP91 vel and AK 35 vel respectively Each contains three columns where column 1 is depth column 2 is Vp and column 3 is Vs Each row represents a new depth shell beginning at the depth defined in column 1 Other models can easily be built to use in FuncLab In order for FuncLab to recognize these however it is necessary to place the new velocity models in text files under the funclab1 5 3 VELMODELS directory and with the extension vel 6 2 Tomography Models Certain analysis packages such as CCP and GCCP stacking can utilize seismic tomography models for timing corrections with receiver function time to depth conversions or migrations This version of FuncLab 1s distributed with a P wave tomography model of the Pacific Northwest published in West et al 2009 and used for CCP stacking analysis in Eagar et al 2010 Itis placed under the funclab1 5 3
29. formats and header descriptions see Appendix C Updated February 29 2012 6 3 3 Receiver Function Raw Data Directory Structure The raw data directories should be structured in the follow manner RFDATA top directory containing dataset no naming restriction can have multiple datasets in a project Event_2006_319 11_14 13 event directory Event _YYYY_JJJ_HH_MM_SS TA_JO8A_0 6 i eqr radial receiver function file NET_STA_GAUSSIAN i eqr TA_JO8A_0 6 i eqt transverse receiver function file NET_STA_GAUSSIAN i eqt TA_JO8A_0 6 i r radial seismogram file NET_STA_GAUSSIAN i r TA_JO8A_0 6 i t transverse seismogram file NET_STA_GAUSSIAN i t TA_JO8A_0 6 i z vertical seismogram file NET_STA_GAUSSIAN i z TA_KO9QA can have multiple stations stations must have set of 5 files z Event_ can have multiple event directories Receiver function files begin with the station name code contains the Gaussian parameter used to compute the receiver function and ends with i eqr radial or i eqt transverse Files for both the radial and transverse receiver functions must be included for completeness Also it is helpful to include the radial transverse and vertical seismograms to exploit some of the seismogram viewing features within the FuncLab toolbox All five files are required for a single event station pair rrr rrr rr rrr rrr rrr rrr Updated Februa
30. ject mat file that is contained one level into the existing Project directory and select Open The main FuncLab GUI will open with your project and you are ready to begin 4 3 Project mat File and FuncLab Variables FuncLab saves the metadata for each station earthquake pair of waveform files in the Project mat MAT file Metadata for the records are organized into a cell array variable called RecordMetadataStrings and a double array variable called RecordMetadataDoubles where each column contains different information such as station name event origin time waveform sample rate etc 4 3 Main FuncLab GUI The main FuncLab GUI allows the user to interact with the entire dataset The GUI includes four top menus File Editing View and Add ons A drop down menu in the top left defines the manner in which receiver function metadata is organized for visualization in the GUI we refer to these as tables although this is not a formal relational database system Records common to a particular station are viewed in Station Tables The stations are then listed and selectable in the left hand panel Below this panel is listed the total number of tables displayed The records within the Station Table are listed by the event origin time in the center panel The total number of records and number of active records are listed below the panel Information or metadata about the selected record in the center panel are then listed in the right hand pa
31. meter axis parameters and mapping parameters These are accessed through the File FuncLab Preferences menu item Adjustments of map boundaries and axis settings from the defaults are likely required depending on the imaging target of the receiver function study 600 Revert to Defaults lt Student Version gt FuncLab Preferences Time axis increment Ray Parameter axis increment 10 0 01 Backazimuth axis increment 45 Begin Time 0 End Time 90 Time series resample rate 0 5 Phase arrivals P35s P35pPs P35sPs 1 D Velocity Model AK135 vel A Maps Locations All a Station Map Projection mercator Event Map Projection eqdazim Maps Background Coastlines low res Station Map North Latitude 46 Station Map South Latitude 40 Station Map West Longitude 123 Station Map East Longitude 113 Event Map Center Latitude 43 Event Map Center Longitude 119 Save Cancel Updated February 29 2012 12 A log file keeps information about each process performed on the dataset The log file can be viewed from the main FuncLab GUI using the menu File Log File This is a text file that 1s initiated when the project is first created Dates and times of processes are logged as well as the name of the process and certain preferences set to run the process This is not editable in the GUI and is intended for user guidance eoo lt Student Version gt FuncLab Log File Project Directory Users keagar research HLP_RFs New_PNW_Mantle Creat
32. nable is turned on eventually but must be initialized as off The Callback is followed by add on Callback and is important This refers to a function or m file by the same name This brings us to the third required component a callback function that transfers the FuncLab GUI handles to your add on GUI This is essential to connecting the core package with the add on This callback function must be named add on Callback m and must be in the add on directory As before replace add on with its proper name In the CCP example it is called ccp_Callback m It must contain 4 lines of code as follows function ccp_CalLback Ccbo eventdata handles h guidataCgcbf Project getCh project_t string Updated February 29 2012 23 ccpCh Replace ccp Callback with the name of the new callback function and replace ccp h with the m file that begins the add on GUI The input h is the FuncLab GUI handles that MUST pass to the Add on GUI for things to work properly The final step is to create this m file that calls the Add on GUI It can have any structure The only difference is that when you want to use handles that are associated with the main FuncLab GUI you must refer to them in the input structure array This is clearer using the examples provided with this manual The CCP GUI contains editable parameters on the left hand side and functional pushbuttons for data processing on the right hand
33. ncy domain or the time domain The most common method is done in the frequency domain but requires water level stabilization to eliminate spectral holes that result from division by zero in the frequency domain Clayton and Wiggins 1976 Although the algorithm is computationally fast it has the disadvantage of requiring the user to manually search for the optimal water level while trying to minimize the loss of spectral information Side lobes on Pds arrivals are also an adverse consequence of spectral leakage from this regularization method An alternative frequency domain method is the multi taper method which reduces the loss of spectral information but requires greater computation Park and Levin 2000 Time domain methods include Updated February 29 2012 5 matrix inversions e g Abers et al 1995 Gurrola et al 1995 Sheehan et al 1995 which also require regularization in the form of damping parameters and the iterative deconvolution method of Ligorria and Ammon 1999 which is more appropriately characterized as a predictive convolution that produces truly causal signals with no loss of spectral information 3 2 Receiver Function Format and Preparation In order for FuncLab to function in a more general sense with publicly available data certain choices were necessary regarding data formats and directory structures The current version is only capable of importing seismic data that is in the SAC format http www iris edu
34. nel These include station and event information time series information and the path of the SAC files for the record Updated February 29 2012 10 File Editing View Export Add ons Stations 310 Events 678 Active 108 The user can also view the records organized by event origin by selecting the Event Tables in the drop down menu Other options in the drop down menu become available when using add on analysis tools discussed below File Editing View Export Add ons File Editing View Export Add ons Y Station Tables ee rs 2004 001 04 39 12 004 00 w be Events 1597 Stations 25 Active 0 Message Event Tables loaded into workspace Updated February 29 2012 11 The user can also view the records organized by event origin by selecting the Event Tables in the drop down menu Other options in the drop down menu become available when using add on analysis tools discussed below The File menu contains options that control FuncLab and projects as a whole such as starting or loading projects adding data and changing preferences New data may be added to an existing project by selecting File Add Data from the top menu 00 About FuncLab New Project Load Project Add Data FuncLab Preferences Log File Exit FuncLab Editing View Export Add ons E FuncLab preferences control certain aspects of all the plotting functions such as time axis parameters backazimuth and ray para
35. nsitive to the local earth structure beneath a station We refer the reader to Langston 1979 and Vinnik 1977 for further details on method theory and development There are several choices in source normalization procedures to isolate the SV energy including the choices of coordinate system and deconvolution method Regarding the choice of the coordinate system Langston 1979 rotated the horizontal component seismograms along the free surface to point in the radial and tangential directions from the backazimuth of the earthquake ZRT coordinate system Vinnik 1977 performed a more complete rotation to the LQR coordinate system where the radial component includes the raypath of the incident wave in the radial direction assuming a slowness vector from a 1 D Earth model Bostock 1996 improved on the method of Vinnik 1977 by maximizing the SV wave energy using slowness as a function of time from the radial and vertical components covariance matrix over a given time window Although these successively improve upon each other it is still difficult to achieve full P SV wavefield separation and all methods represent an approximation For the purposes of most imaging however this approximation 1s acceptable as the signal of the Pds is visible with any of these rotation choices The second choice in source normalization involves the deconvolution method of the source from the seismograms Deconvolution can either be performed in the freque
36. ry 29 2012 4 Graphical User Interface To begin the user invokes the main FuncLab GUI using the funclab command in MATLAB From this interface the user may start a new project or load an existing one AAC lt Student Version gt FuncLab 1 5 File Editing View Export Add ons Welcome to FuncLab Start a new project File gt New Project or load an existing project File gt Load Project 4 1 Starting a new project To begin by starting a new project where the data has already been pre processed outside of FuncLab select File gt New Project from the top menu oo l Editing View Export Add ons About FuncLab New Project Load Project FuncLab Preferences Exit FuncLab The New FuncLab Project Setup interface will appear From here the user will first determine the directory for the new project This directory must not exist as FuncLab will create and populate it with the necessary files and sub directories for the project Up to 10 separate data directories formatted in the required structure as discussed above can be imported at the initiation of the project The user may either enter these directories Updated February 29 2012 8 manually in the specified boxes or browse for these directories on the computer After entering as many data directories as the user wishes to import for this project up to 10 but more can be added once the project is created the user hits the Start button
37. s keagar research HLP_RFs data_0 6 CC_ALL_04 to Users keagar research HLP_RFs Ne Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA CC_ALL_ Copied Data Directory Users keagar research HLP_RFs data_0 6 CC_ALL_05 to Users keagar research HLP_RFs Ne Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA CC_ALL Copied Data Directory Users keagar research HLP_RFs data_0 6 CC_ALL_06 to Users keagar research HLP_RFs Ne Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA CC_ALL_ Copied Data Directory Users keagar research HLP_RFs data_0 6 CC_ALL_07 to Users keagar research HLP_RFs Ne 20 Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA CC_ALL_ 04 30 2011 09 15 31 Created and Saved New Trace Editing Preferences 04 30 2011 17 33 43 Add Data to FuncLab Project Moved Data Directory Users keagar research HLP_RFs data_0 6 XC_ALL_06 to Users keagar research HLP_RFs New_ Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA XC_ALL_0 Moved Data Directory Users keagar research HLP_RFs data_0 6 XC_ALL_07 to Users keagar research HLP_RFs New Imported Data To Station and Event Tables Users keagar research HLP_RFs New_PNW_Mantle RAWDATA XC_ALL_07 Moved Data Directory Users keagar research HLP_RFs data_0 6
38. up funclab script to the startup m file located under the MATLAB application directory See MATLAB documentation for startup m for specifics see http www mathworks com help techdoc ref startup html This will then enable setup funclab to run automatically at the start of MATLAB FuncLab is now ready to be used on your machine Updated February 29 2012 4 3 Receiver Function Overview Preparation and Raw Data Directory Structure 3 1 Receiver Function Overview P to S converted waves Pds that arrive directly after the P wave provide information about velocity and density discontinuities in the Earth below seismic stations These phases may be investigated directly on the horizontal components if the recorded earthquakes have simple impulsive source time functions However many earthquake source wavelets are complicated and generally unknown Since the Pds arrivals lie within the coda of the P wave and have amplitudes an order of magnitude smaller stacking many records together to enhance the signal to noise ratio is necessary Both Langston 1979 and Vinnik 1977 recognized the need for an approach to compare converted phases from events with different source signals and developed similar methods for isolating the SV from the P wave signal and removing the influence of the source signal from the seismograms The receiver function therefore represents the Earth s impulse response containing time and amplitude information that is se
39. utomated trace editing e200 View Export Add ons Manual Trace Edit SS S Auto Trace Edit Selected bm TA Auto Trace Edit ALL L Trace Edit Import T TA Trace Editing Preferences Updated February 29 2012 13 Opening the Manual Trace Edit GUI will read in the receiver function SAC files of the records in the selected table in the left hand panel The Trace Editing GUI displays the receiver functions as wiggle traces with time increasing down the y axis and each record one unit wide across the x axis The scrollbar at the bottom allows the user to scan through records in the table to display them in the window The record status is changed in the Status column of RecordMetadataDoubles and saved to Project mat at the end of the trace editing session Each record has a checkbox above it which is used to select the status of the record check means on and no check means off The ability to change a record s status allows the user to specify good records that will be used in future analyses 1 e stacking and visualization processes and ignore bad records Records with status off remain part of the dataset and are not deleted from the project The color of the wiggle traces is another indicator of the record s status New data that have not yet been evaluated through the Trace Editor are colored red whereas data that has been looked at in the GUI before are colored blue All records with status o
Download Pdf Manuals
Related Search