Location

2007/11/29 19:00:19 14.94 -61.24 146

Moment Tensor Comparison

The following compares this source inversion to the USGS Rapid Moment Tensor Solution and to the Harvard CMT solutions, if they are available.

SPECIAL NOTE

These waveforms were taken from the Orfeus Wilbur II for the purpose of testing the dataless SEED instrument responses. The Mednet (MN) stations CEL, TIP and TUE could not be used since the RESP files otained from the dataless SEED at the Orfeus data center could not parsed using evalresp.

The earthquake selected, the 2007/11/29 Martinique earthquake, has a source time function more complicated that the modeling assumption of a single, smoothed pulse. Do not focus on the RED dots, which indicate a devision of more that a factor of two, may not be indicative of the quality of the waveforms. Focus on large differences, e.g., CEYBHRSL or others. EIL is another station with an incorrect response.

Obviously this is just a demonstration of how teleseismic modeling can be used to test a network performance.

Data set

The following broadband stations passed the QC and were used for the source inversion. AIGLE ANTO APE AQU ARBF ARE0 ARG ARSA ARU BALST BERNI BFO BNI BOURR BRANT BRG BSD BSEG BUC1 BUD BUG CADS CALF CEL CEY CFF CHIF CLL COP CRAR CRES CRVS CSOR CSS CUC DAVA DAVOX DIX DPC DRGR DSB EBR ECAL ECH EDI EIL EJON EKB EMUR EMV EOSO ESCA ESK ESPZ FUORN FUR FUSIO GIMEL GKP GROS HASLI HEF HGN HLG HTL IAS IBBN IDI ISP ITM JAVC JAVS JMB JOF KBA KBS KEK KHC KIV KNDS KONO KPL KRUC KSP KWP LAST LIENZ LJU LLS LRW LTK MAHO MCH1 MMK MOA MOX MTE MUD MUGIO MUO NAO01 NE05 NKC OBKA OBN OGAG OGDI OJC OKC OPLO ORCD OUL PERS PESTR PKSM PLONS PSZ PVAQ RAF RDO RENF RGN ROBS RTC RUE RUSF SAOF SBD SENIN SFJD SFS SLE SMAI SOP SSB STU SUF SULZ SWN1 TAM TIP TIR TIRR TORNY TRI TUE VAF VDL VISS VLC VRAC VRI VSL VSU VYHS WAR WDD WET WILA WIMIS WIT WLF WOL WTSB WTTA ZKR ZST ZUR

Forward Synthetic

All observed and Greens function waveforms are corrected to instrument response to ground velocity in meters/sec for the passband of 0.004 - 5 Hz. The traces were then lowpass filtered at 0.25 Hz and interpolated to a sample rate of 1 second.

The moment tensor solution used has the parameters:

HS=146
STK=109
DIP=59
RAKE=-31
MW=7.4 

The Green's function closest to the desired depth was
DEP=1500 , where (DEP/10) is the computed depth.
The cut windows attempt to include the P, pP, sP, pS, S and sS arrivals. However, one must be very careful about the fact that PP may be included in some distance ranges.

The waveforms are then bandpass filtered by the application of the following high- and low-pass stages (an optional microseism filter):

hp c 0.0070 2
lp c 0.0125 2
int
br c 0.12 0.25 n 4 p 2
The traces were next integrated to ground displacment in meters. Finally the observed data are interpolated to have the same sampling at the Green's functions.
NOTE: this was done for speed. The proper sequence is to read traces, filter and then cut - gsac will be modified to introduce a command CUTWR to define the cut upon a write.

The following filter ranges are used according to the moment magnitude Mw:

   FILTER_BAND   1/FH(s)  1/FL(s)
       1           60     12      Mw < 6.4
       2          100     20      6.4 < Mw <= 6.9
       3          120     40      Mw > 6.9

The map displays the distribution of stations used for this source inversion.


Location of the earthquake (yellow star) and great circle path from the epicenter to each station (red) [created using GMT (Wessel, P., and W. H. F. Smith, New version of Generic Mapping Tools released, EOS Trans. AGU, 76 329, 1995.)]

Detailed Waveform Comparison

P-wave Z component
Observed (red) and predicted seismograms (blue) ordered by increasing epicentral distance. Each pair of traces is annotated with the wave type and a station identifier (station, network, and channel id's), epicentral distance in degrees, source-to-station azimuth in degrees. Each seismogram pair is plotted with the same scale and the peak amplitudes in meters are shown above to the left of each trace. Red circles flag seismograms with amplitude misfits of a factor of 2 or more.
SH-wave T component
Observed (red) and predicted seismograms (blue) ordered by increasing epicentral distance. Each pair of traces is annotated with the wave type and a station identifier (station, network, and channel id's), epicentral distance in degrees, source-to-station azimuth in degrees. Each seismogram pair is plotted with the same scale and the peak amplitudes in meters are shown above to the left of each trace. Red circles flag seismograms with amplitude misfits of a factor of 2 or more.
SV-wave R component
Observed (red) and predicted seismograms (blue) ordered by increasing epicentral distance. Each pair of traces is annotated with the wave type and a station identifier (station, network, and channel id's), epicentral distance in degrees, source-to-station azimuth in degrees. Each seismogram pair is plotted with the same scale and the peak amplitudes in meters are shown above to the left of each trace. Red circles flag seismograms with amplitude misfits of a factor of 2 or more.

Waveforms used

The following waveforms are used for the plot above. Note that the naming convention used by the generating program does not distinguish between LH, BH and HH channels. One can change the DOPZSTA, DOSHTSTA and DOSVRSTA scripts to preserve the component names.

AIGLECHBHR AIGLECHBHT AIGLECHBHZ ANTOIU00BHR ANTOIU00BHT ANTOIU00BHZ APEGEBHZ AQUMNBHR AQUMNBHT AQUMNBHZ ARBFFRBHR ARBFFRBHT ARBFFRBHZ ARE0NOBHR ARE0NOBHT ARE0NOBHZ ARGHLBHR ARGHLBHT ARGHLBHZ ARSAOEBHR ARSAOEBHT ARSAOEBHZ ARUII00BHR ARUII00BHT ARUII00BHZ BALSTCHBHR BALSTCHBHT BALSTCHBHZ BERNICHBHR BERNICHBHT BERNICHBHZ BFOII00BHR BFOII00BHT BFOII00BHZ BNIMNBHR BNIMNBHT BNIMNBHZ BOURRCHBHR BOURRCHBHT BOURRCHBHZ BRANTCHBHR BRANTCHBHT BRANTCHBHZ BRGGRBHR BRGGRBHT BRGGRBHZ BSDDKBHR BSDDKBHT BSDDKBHZ BSEGGRBHR BSEGGRBHT BSEGGRBHZ BUC1ROBHR BUC1ROBHT BUC1ROBHZ BUDHUBHR BUDHUBHT BUDHUBHZ BUGGRBHR BUGGRBHT BUGGRBHZ CADSSLBHR CADSSLBHT CADSSLBHZ CALFFRBHR CALFFRBHT CALFFRBHZ CELMNBHR CELMNBHT CELMNBHZ CEYSLBHR CEYSLBHT CEYSLBHZ CFFFRBHR CFFFRBHT CFFFRBHZ CHIFFRBHR CHIFFRBHT CHIFFRBHZ CLLGRBHR CLLGRBHT CLLGRBHZ COPDKBHR COPDKBHT COPDKBHZ CRARROBHR CRARROBHT CRARROBHZ CRESSLBHR CRESSLBHT CRESSLBHZ CRVSSKBHR CRVSSKBHT CRVSSKBHZ CSORCABHR CSORCABHT CSORCABHZ CSSGEBHR CSSGEBHT CSSGEBHZ CUCMNBHR CUCMNBHT CUCMNBHZ DAVAOEBHR DAVAOEBHT DAVAOEBHZ DAVOXCHBHR DAVOXCHBHT DAVOXCHBHZ DIXCHBHR DIXCHBHT DIXCHBHZ DPCCZBHR DPCCZBHT DPCCZBHZ DRGRROBHR DRGRROBHT DRGRROBHZ DSBGEBHR DSBGEBHT DSBGEBHZ EBRSSBHR EBRSSBHT EBRSSBHZ ECALESBHR ECALESBHT ECALESBHZ ECHG10BHR ECHG10BHT ECHG10BHZ EDIGBBHR EDIGBBHT EDIGBBHZ EILGEBHR EILGEBHT EILGEBHZ EJONESBHR EJONESBHT EJONESBHZ EKBBNBHR EKBBNBHT EKBBNBHZ EMURESBHR EMURESBHT EMURESBHZ EMVCHBHR EMVCHBHT EMVCHBHZ EOSOESBHR EOSOESBHT EOSOESBHZ ESCAFRBHR ESCAFRBHT ESCAFRBHZ ESKII00BHR ESKII00BHT ESKII00BHZ ESPZAIBHR ESPZAIBHT ESPZAIBHZ FUORNCHBHR FUORNCHBHT FUORNCHBHZ FURGRBHR FURGRBHT FURGRBHZ FUSIOCHBHR FUSIOCHBHT FUSIOCHBHZ GIMELCHBHR GIMELCHBHT GIMELCHBHZ GKPPLBHR GKPPLBHT GKPPLBHZ GROSSLBHR GROSSLBHT GROSSLBHZ HASLICHBHR HASLICHBHT HASLICHBHZ HEFHEBHR HEFHEBHT HEFHEBHZ HGNNL01BHR HGNNL01BHT HGNNL01BHZ HLGGEBHR HLGGEBHT HLGGEBHZ HTLGBBHR HTLGBBHT HTLGBBHZ IASROBHR IASROBHT IASROBHZ IBBNGEBHR IBBNGEBHT IBBNGEBHZ IDIMNBHR IDIMNBHT IDIMNBHZ ISPGEBHR ISPGEBHT ISPGEBHZ ITMHLBHR ITMHLBHT ITMHLBHZ JAVCCZBHR JAVCCZBHT JAVCCZBHZ JAVSSLBHR JAVSSLBHT JAVSSLBHZ JMBBSBHR JMBBSBHT JMBBSBHZ JOFHEBHR JOFHEBHT JOFHEBHZ KBAOEBHR KBAOEBHT KBAOEBHZ KBSGE00BHR KBSGE00BHT KBSGE00BHZ KEKHLBHR KEKHLBHT KEKHLBHZ KHCCZBHR KHCCZBHT KHCCZBHZ KIVII00BHR KIVII00BHT KIVII00BHZ KNDSSLBHR KNDSSLBHT KNDSSLBHZ KONOIU00BHR KONOIU00BHT KONOIU00BHZ KONOIU10BHR KONOIU10BHT KONOIU10BHZ KPLGBBHR KPLGBBHT KPLGBBHZ KRUCCZBHR KRUCCZBHT KRUCCZBHZ KSPPLBHR KSPPLBHT KSPPLBHZ KWPGEBHR KWPGEBHT KWPGEBHZ LASTGEBHR LASTGEBHT LASTGEBHZ LIENZCHBHR LIENZCHBHT LIENZCHBHZ LJUSLBHR LJUSLBHT LJUSLBHZ LLSCHBHR LLSCHBHT LLSCHBHZ LRWGBBHR LRWGBBHT LRWGBBHZ LTKHPBHR LTKHPBHT LTKHPBHZ MAHOGEBHR MAHOGEBHT MAHOGEBHZ MCH1GBBHR MCH1GBBHT MCH1GBBHZ MMKCHBHR MMKCHBHT MMKCHBHZ MOAOEBHR MOAOEBHT MOAOEBHZ MOXGRBHR MOXGRBHT MOXGRBHZ MTEGEBHR MTEGEBHT MTEGEBHZ MUDDKBHR MUDDKBHT MUDDKBHZ MUGIOCHBHR MUGIOCHBHT MUGIOCHBHZ MUOCHBHR MUOCHBHT MUOCHBHZ NAO01NOBHR NAO01NOBHT NAO01NOBHZ NE05NRBHR NE05NRBHT NE05NRBHZ NKCCZBHR NKCCZBHT NKCCZBHZ OBKAOEBHR OBKAOEBHT OBKAOEBHZ OBNII00BHR OBNII00BHT OBNII00BHZ OGAGFRBHR OGAGFRBHT OGAGFRBHZ OGDIFRBHR OGDIFRBHT OGDIFRBHZ OJCPLBHR OJCPLBHT OJCPLBHZ OKCCZBHR OKCCZBHT OKCCZBHZ OPLONLBHR OPLONLBHT OPLONLBHZ ORCDAIBHR ORCDAIBHT ORCDAIBHZ OULFNBHR OULFNBHT OULFNBHZ PERSSLBHR PERSSLBHT PERSSLBHZ PESTRPMBHR PESTRPMBHT PESTRPMBHZ PKSMHUBHR PKSMHUBHT PKSMHUBHZ PLONSCHBHR PLONSCHBHT PLONSCHBHZ PSZGEBHR PSZGEBHT PSZGEBHZ PVAQPMBHR PVAQPMBHT PVAQPMBHZ RAFHEBHR RAFHEBHT RAFHEBHZ RDOHLBHR RDOHLBHT RDOHLBHZ RENFFRBHR RENFFRBHT RENFFRBHZ RGNGEBHR RGNGEBHT RGNGEBHZ ROBSSLBHR ROBSSLBHT ROBSSLBHZ RTCMNBHR RTCMNBHT RTCMNBHZ RUEGEBHR RUEGEBHT RUEGEBHZ RUSFFRBHR RUSFFRBHT RUSFFRBHZ SAOFFRBHR SAOFFRBHT SAOFFRBHZ SBDBNBHR SBDBNBHT SBDBNBHZ SENINCHBHR SENINCHBHT SENINCHBHZ SFJDIU00BHR SFJDIU00BHT SFJDIU00BHZ SFJDIU10BHR SFJDIU10BHT SFJDIU10BHZ SFSGEBHR SFSGEBHT SFSGEBHZ SLECHBHR SLECHBHT SLECHBHZ SMAIAIBHR SMAIAIBHT SMAIAIBHZ SOPHUBHR SOPHUBHT SOPHUBHZ SSBG00BHR SSBG00BHT SSBG00BHZ STUGEBHR STUGEBHT STUGEBHZ SUFHEBHR SUFHEBHT SUFHEBHZ SULZCHBHR SULZCHBHT SULZCHBHZ SWN1GBBHR SWN1GBBHT SWN1GBBHZ TAMG00BHR TAMG00BHT TAMG00BHZ TIPMNBHR TIPMNBHT TIPMNBHZ TIRMNBHR TIRMNBHT TIRMNBHZ TIRRGEBHR TIRRGEBHT TIRRGEBHZ TORNYCHBHR TORNYCHBHT TORNYCHBHZ TRIMNBHR TRIMNBHT TRIMNBHZ TUEMNBHR TUEMNBHT TUEMNBHZ VAFHEBHR VAFHEBHT VAFHEBHZ VDLCHBHR VDLCHBHT VDLCHBHZ VISSSLBHR VISSSLBHT VISSSLBHZ VLCMNBHR VLCMNBHT VLCMNBHZ VRACCZBHR VRACCZBHT VRACCZBHZ VRIROBHR VRIROBHT VRIROBHZ VSLMNBHR VSLMNBHT VSLMNBHZ VSUGEBHR VSUGEBHT VSUGEBHZ VYHSSKBHR VYHSSKBHT VYHSSKBHZ WARPLBHR WARPLBHT WARPLBHZ WDDMNBHR WDDMNBHT WDDMNBHZ WETGRBHR WETGRBHT WETGRBHZ WILACHBHR WILACHBHT WILACHBHZ WIMISCHBHR WIMISCHBHT WIMISCHBHZ WITNLBHR WITNLBHT WITNLBHZ WLFGEBHR WLFGEBHT WLFGEBHZ WOLBNBHR WOLBNBHT WOLBNBHZ WTSBNL01BHR WTSBNL01BHT WTSBNL01BHZ WTTAOEBHR WTTAOEBHT WTTAOEBHZ ZKRGEBHR ZKRGEBHT ZKRGEBHZ ZSTSKBHR ZSTSKBHT ZSTSKBHZ ZURCHBHR ZURCHBHT ZURCHBHZ

Processing times

Starting Processing         : Fri Dec 7 21:26:03 UTC 2007
Unpacking SEED volume       : Fri Dec 7 21:26:04 UTC 2007
Starting deconvolution      : Fri Dec 7 21:26:04 UTC 2007
Starting trace rotation     : Fri Dec 7 22:04:34 UTC 2007
Starting distance selection : Fri Dec 7 22:09:07 UTC 2007
Starting trace QC           : Fri Dec 7 22:11:08 UTC 2007
Starting synthetic          : Fri Dec 7 22:11:09 UTC 2007
Starting documentation      : Fri Dec 7 22:28:46 UTC 2007

Last Changed Fri Dec 7 22:29:17 UTC 2007