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.)]
2007/02/24 02:36:22 -6.90 -80.32 23
The following compares this source inversion to the USGS Rapid Moment Tensor Solution and to the Harvard CMT solutions, if they are available.
SLU Moment Tensor Solution
2007/02/24 02:36:22
Best Fitting Double Couple
Mo = 2.43e+25 dyne-cm
Mw = 6.19
Z = 25 km
Plane Strike Dip Rake
NP1 162 73 100
NP2 310 20 60
Principal Axes:
Axis Value Plunge Azimuth
T 2.43e+25 61 87
N 0.00e+00 10 338
P -2.43e+25 27 243
Moment Tensor: (dyne-cm)
Component Value
Mxx -3.84e+24
Mxy -7.37e+24
Mxz 5.00e+24
Myy -9.67e+24
Myz 1.91e+25
Mzz 1.35e+25
--------------
############----------
##-----##############-------
--------#################-----
----------###################-----
-----------#####################----
------------######################----
--------------######################----
--------------#######################---
----------------########### #########---
----------------########### T #########---
-----------------########## #########---
-----------------#######################--
----- ---------#####################--
----- P ----------####################--
---- -----------###################-
------------------#################-
------------------################
-----------------#############
-----------------###########
---------------#######
-------------#
Harvard Convention
Moment Tensor:
R T F
1.35e+25 5.00e+24 -1.91e+25
5.00e+24 -3.84e+24 7.37e+24
-1.91e+25 7.37e+24 -9.67e+24
|
07/02/24 02:36:22.08
NEAR COAST OF NORTHERN PERU
Epicenter: -6.901 -80.322
MW 6.3
USGS MOMENT TENSOR SOLUTION
Depth 23 No. of sta: 69
Moment Tensor; Scale 10**18 Nm
Mrr= 1.95 Mtt= 0.13
Mpp=-2.08 Mrt= 0.46
Mrp=-2.82 Mtp=-0.15
Principal axes:
T Val= 3.47 Plg=62 Azm= 73
N 0.06 7 175
P -3.53 27 269
Best Double Couple:Mo=3.5*10**18
NP1:Strike= 15 Dip=19 Slip= 111
NP2: 173 73 83
-######
------#########--
-------############--
---------##############--
-----------###############---
------------################---
------------#################--
-------------####### #######---
--- -------####### T #######---
--- P -------####### #######---
--- -------#################---
--------------###############----
-------------###############---
--------------#############----
-------------############----
------------#########----
----------#######----
---------###-----
#####--
|
The following broadband stations passed the QC and were used for the source inversion. AAM ACSO AGMN AHID AMTX ANMO ASCN BBSR BFO BLA BMO BOZ BRAL BW06 CBKS CBN CNNC COLA COWI DGMT DUG DWPF ECSD EFI EGAK EGMT ELK ESK EYMN FFC GLMI GOGA HAWA HDIL HKT HLID HRV HWUT ISCO JCT JFWS KDAK KIP KSU1 KVTX LAO LBNH LONY LRAL MCWV MIAR MSO MVCO NATX NCB NEW NHSC NLWA PAB PFO PKME PMSA POHA QSPA RCBR RLMT RSSD SBA SCIA SSPA TPNV TRIS TUC TZTN WMOK WRAK WUAZ WVOR
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.
For the grid search, the observed traces and Green's functions are read in an cut using the following commands
Phase Gsac Command Comment P cut A -30 A CUTH = 95+0.3*DEPTH SH cut T1 -60 T1 CUTH = 95+0.6*DEPTH SV cut T0 -60 T0 CUTH = 95+0.6*DEPTH where the 95 is a maximum filter duration, DEPTH is in km, and the CUTH in secThe cut windows attempt to include the P, pP, sP, pS, S and sS arrivals. However, oen 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.0100 2 lp c 0.0400 2 int br c 0.12 0.2 n 4 p 2The traces were next integrated to ground displacment in meters. Finally the observed data are interpolated to ahve the same sampling at the Green's functions.
The source inversion is a multipass operation since a lower frequency filter band is used for larger earthquakes and since a search is made over depth. Up to three passed of the outer loop are made, after which the moment magnitude is determined and filter settings readjusted. The inner loop over depth samples all depths from 0 to 800 km with 5 km increments in depth to 50 km, followed by 10 km depth sampling for the remaining range.
The following filter ranges are used according to the moment magnitude Mw:
FILTER_BAND FH(s) 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.)] |
For this data set the favored solution is
WVFMTD96 20.0 166. 82. 90. 6.09 0.665 0.102E-05 0.639 0.816 0.551E-06 16.3
The following figures show the sensitivity of the goodness of fit parameter so source depth, the waveform comparison as a function of epicentral distance in degrees and the source to station azimuth
|
| Goodness of fit as a function of source depth. The measure is 1 - SUM (o -p)2 / SUM o2. A value of 1.0 is the best fit. The best double couple mechanism for the solution depth is plotted above goodness of fit value to indicate how the mefhanism may change with depth. |
| P-wave Z component |
|
| Comparison of the observed traces (red) and solution predicted traces (blue) ordered in terms of increasing epicentral distance. Each pair of traces is annotated with the station name, epicentral distance in degrees, source to station azimuth in degrees. Each pair of traces is plotted with the same scale and the peak amplitudes are indicated at the lect of each trace. Finally the time shift between the P-wave first arrival picked and the the theoretical P-wave first arrival in the predicted trace is indicated, with a positive sign indicating that the predicted trace has been shifted to the right by the given number of seconds. as a function of source to station azimuth in degrees (D). The purpose of this display is to highlight the azimuthal dependence on the first motion. The traces are annotated with the station name at the top. |
| SH-wave T component |
|
| Comparison of the observed traces (red) and solution predicted traces (blue) ordered in terms of increasing epicentral distance. Each pair of traces is annotated with the station name, epicentral distance in degrees, source to station azimuth in degrees. Each pair of traces is plotted with the same scale and the peak amplitudes are indicated at the lect of each trace. Finally the time shift between the P-wave first arrival picked and the the theoretical P-wave first arrival in the predicted trace is indicated, with a positive sign indicating that the predicted trace has been shifted to the right by the given number of seconds. as a function of source to station azimuth in degrees (D). The purpose of this display is to highlight the azimuthal dependence on the first motion. The traces are annotated with the station name at the top. |
| SV-wave R component |
|
| Comparison of the observed traces (red) and solution predicted traces (blue) ordered in terms of increasing epicentral distance. Each pair of traces is annotated with the station name, epicentral distance in degrees, source to station azimuth in degrees. Each pair of traces is plotted with the same scale and the peak amplitudes are indicated at the lect of each trace. Finally the time shift between the P-wave first arrival picked and the the theoretical P-wave first arrival in the predicted trace is indicated, with a positive sign indicating that the predicted trace has been shifted to the right by the given number of seconds. as a function of source to station azimuth in degrees (D). The purpose of this display is to highlight the azimuthal dependence on the first motion. The traces are annotated with the station name at the top. |
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.
For the grid search, the observed traces and Green's functions are read in an cut using the following commands
Phase Gsac Command Comment P cut A -30 A CUTH = 95+0.3*DEPTH SH cut T1 -60 T1 CUTH = 95+0.6*DEPTH SV cut T0 -60 T0 CUTH = 95+0.6*DEPTH where the 95 is a maximum filter duration, DEPTH is in km, and the CUTH in secThe cut windows attempt to include the P, pP, sP, pS, S and sS arrivals. However, oen 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.0100 2 lp c 0.0400 2 int br c 0.12 0.2 n 4 p 2The traces were next integrated to ground displacment in meters. Finally the observed data are interpolated to ahve the same sampling at the Green's functions.
The source inversion is a multipass operation since a lower frequency filter band is used for larger earthquakes and since a search is made over depth. Up to three passed of the outer loop are made, after which the moment magnitude is determined and filter settings readjusted. The inner loop over depth samples all depths from 0 to 800 km with 5 km increments in depth to 50 km, followed by 10 km depth sampling for the remaining range.
The following filter ranges are used according to the moment magnitude Mw:
FILTER_BAND FH(s) 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.)] |
For this data set the favored solution is
WVFGRD96 25.0 310 20 60 6.19 0.6771
The following figures show the sensitivity of the goodness of fit parameter so source depth, the waveform comparison as a function of epicentral distance in degrees and the source to station azimuth
|
| Goodness of fit as a function of source depth. The measure is 1 - SUM (o -p)2 / SUM o2. A value of 1.0 is the best fit. The best double couple mechanism for the solution depth is plotted above goodness of fit value to indicate how the mefhanism may change with depth. |
| P-wave Z component |
|
| Comparison of the observed traces (red) and solution predicted traces (blue) ordered in terms of increasing epicentral distance. Each pair of traces is annotated with the station name, epicentral distance in degrees, source to station azimuth in degrees. Each pair of traces is plotted with the same scale and the peak amplitudes are indicated at the lect of each trace. Finally the time shift between the P-wave first arrival picked and the the theoretical P-wave first arrival in the predicted trace is indicated, with a positive sign indicating that the predicted trace has been shifted to the right by the given number of seconds. as a function of source to station azimuth in degrees (D). The purpose of this display is to highlight the azimuthal dependence on the first motion. The traces are annotated with the station name at the top. |
| SH-wave T component |
|
| Comparison of the observed traces (red) and solution predicted traces (blue) ordered in terms of increasing epicentral distance. Each pair of traces is annotated with the station name, epicentral distance in degrees, source to station azimuth in degrees. Each pair of traces is plotted with the same scale and the peak amplitudes are indicated at the lect of each trace. Finally the time shift between the P-wave first arrival picked and the the theoretical P-wave first arrival in the predicted trace is indicated, with a positive sign indicating that the predicted trace has been shifted to the right by the given number of seconds. as a function of source to station azimuth in degrees (D). The purpose of this display is to highlight the azimuthal dependence on the first motion. The traces are annotated with the station name at the top. |
| SV-wave R component |
|
| Comparison of the observed traces (red) and solution predicted traces (blue) ordered in terms of increasing epicentral distance. Each pair of traces is annotated with the station name, epicentral distance in degrees, source to station azimuth in degrees. Each pair of traces is plotted with the same scale and the peak amplitudes are indicated at the lect of each trace. Finally the time shift between the P-wave first arrival picked and the the theoretical P-wave first arrival in the predicted trace is indicated, with a positive sign indicating that the predicted trace has been shifted to the right by the given number of seconds. as a function of source to station azimuth in degrees (D). The purpose of this display is to highlight the azimuthal dependence on the first motion. The traces are annotated with the station name at the top. |