[cig-commits] r11880 - in seismo/3D/automeasure/latex: . figures/japan
carltape at geodynamics.org
carltape at geodynamics.org
Thu May 1 07:21:15 PDT 2008
Author: carltape
Date: 2008-05-01 07:21:15 -0700 (Thu, 01 May 2008)
New Revision: 11880
Added:
seismo/3D/automeasure/latex/figures/japan/SHR_BO_200511211536A.pdf
seismo/3D/automeasure/latex/figures/japan/stats_T06.pdf
Removed:
seismo/3D/automeasure/latex/figures/japan/T06_rs-crop.pdf
Modified:
seismo/3D/automeasure/latex/AM-allcitations.bib
seismo/3D/automeasure/latex/appendix.tex
seismo/3D/automeasure/latex/figures/japan/051502B_T06_rs.pdf
seismo/3D/automeasure/latex/figures/japan/091502B_T06_rs.pdf
seismo/3D/automeasure/latex/figures/japan/200511211536A_T06_rs.pdf
seismo/3D/automeasure/latex/figures/japan/200511211536A_T24_rs.pdf
seismo/3D/automeasure/latex/figures/japan/ERM_II_051502B.pdf
seismo/3D/automeasure/latex/figures/japan/KIS_BO_091502B.pdf
seismo/3D/automeasure/latex/figures_paper.tex
seismo/3D/automeasure/latex/flexwin_paper.pdf
seismo/3D/automeasure/latex/results.tex
Log:
Checked in Min's modifications to flexwin paper.
Modified: seismo/3D/automeasure/latex/AM-allcitations.bib
===================================================================
--- seismo/3D/automeasure/latex/AM-allcitations.bib 2008-05-01 04:11:44 UTC (rev 11879)
+++ seismo/3D/automeasure/latex/AM-allcitations.bib 2008-05-01 14:21:15 UTC (rev 11880)
@@ -223,15 +223,25 @@
Volume = {96},
Year = {2006}}
+ at article{KennettEngdahl1991,
+ Author = {B. L. N. {Kennett} and E. R. {Engdahl}},
+ Date-Added = {2008-04-30 01:46:39 +0100},
+ Date-Modified = {2008-04-30 01:47:12 +0100},
+ Journal = {Geophys. J. Int.},
+ Pages = {429--465},
+ Title = {{Traveltimes for global earthquake location and phase identification}},
+ Volume = {105},
+ Year = {1991}}
+
@article{KomatitschEtal2004,
- Author = {Komatitsch, D. and Liu, Q. and Tromp, J. and S\"uss, P. and Stidham, C. and Shaw J.H.},
- Date-Added = {2008-03-22 19:38:39 +0100},
- Date-Modified = {2008-03-22 19:41:12 +0100},
- Journal = {B. Seismol. Soc. Am.},
- Pages = {187--206},
- Title = {{Simulations of ground motion in the Los Angeles basin based upon the spectral-element method}},
- Volume = {94},
- Year = {2004}}
+ Author = {Komatitsch, D. and Liu, Q. and Tromp, J. and S\"uss, P. and Stidham, C. and Shaw J.H.},
+ Date-Added = {2008-03-22 19:38:39 +0100},
+ Date-Modified = {2008-03-22 19:41:12 +0100},
+ Journal = {B. Seismol. Soc. Am.},
+ Pages = {187--206},
+ Title = {{Simulations of ground motion in the Los Angeles basin based upon the spectral-element method}},
+ Volume = {94},
+ Year = {2004}}
@article{ChenEtal2007,
Author = {Chen, M. and Tromp, J. and Helmberger, D. and Kanamori, H.},
Modified: seismo/3D/automeasure/latex/appendix.tex
===================================================================
--- seismo/3D/automeasure/latex/appendix.tex 2008-05-01 04:11:44 UTC (rev 11879)
+++ seismo/3D/automeasure/latex/appendix.tex 2008-05-01 14:21:15 UTC (rev 11880)
@@ -48,9 +48,40 @@
\pagebreak
\subsection{Japan scenario\label{ap:user_japan}}
+Below $t_P$ and $t_S$ denote the start of the time windows for P wave and S wave, predicted by 1-D IASPEI91 model \citep{KennettEngdahl1991}. The finish time for the surface-wave time window is indicated
+by $t_{R1}$.
-MIN: fill this in please.
+For the \trange{24}{120} data, we use
+%
+\begin{align}
+w_E(t) & =
+ \begin{cases}
+ 10 w_E & \text{$t < t_P$} \\
+ w_E & \text{$t_P \le t \leq t_{R1}$} \\
+ 10 w_E & \text{$t > t_{R1}$}
+ \end{cases}
+\end{align}
+For the \trange{6}{30} data, we avoid selecting surface-wave arrivals, and therefore we make the modifications
+%
+\begin{align}
+w_E(t) & =
+ \begin{cases}
+ 10 w_E & \text{$t < t_P$} \\
+ w_E & \text{$t_P \le t \leq t_S$} \\
+ 10 w_E & \text{$t > t_S$}
+ \end{cases}
+\end{align}
+
+For data in both period ranges, we use
+\begin{align}
+\Delta\tau_0(t) & =
+ \begin{cases}
+ 0.08 \text{$t_P$} & \text{$h \leq$ 70~km} \\
+ max(0.06 \text{$t_P$}, 1.4\tau_0) & \text{70~km $< h <$ 300~km} \\
+ max(0.06 \text{$t_P$}, 1.7\tau_0) & \text{$h \geq$ 300~km}
+ \end{cases}
+\end{align}
%--------------------------
\pagebreak
Modified: seismo/3D/automeasure/latex/figures/japan/051502B_T06_rs.pdf
===================================================================
(Binary files differ)
Modified: seismo/3D/automeasure/latex/figures/japan/091502B_T06_rs.pdf
===================================================================
(Binary files differ)
Modified: seismo/3D/automeasure/latex/figures/japan/200511211536A_T06_rs.pdf
===================================================================
(Binary files differ)
Modified: seismo/3D/automeasure/latex/figures/japan/200511211536A_T24_rs.pdf
===================================================================
(Binary files differ)
Modified: seismo/3D/automeasure/latex/figures/japan/ERM_II_051502B.pdf
===================================================================
(Binary files differ)
Modified: seismo/3D/automeasure/latex/figures/japan/KIS_BO_091502B.pdf
===================================================================
(Binary files differ)
Added: seismo/3D/automeasure/latex/figures/japan/SHR_BO_200511211536A.pdf
===================================================================
(Binary files differ)
Property changes on: seismo/3D/automeasure/latex/figures/japan/SHR_BO_200511211536A.pdf
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Deleted: seismo/3D/automeasure/latex/figures/japan/T06_rs-crop.pdf
===================================================================
(Binary files differ)
Added: seismo/3D/automeasure/latex/figures/japan/stats_T06.pdf
===================================================================
(Binary files differ)
Property changes on: seismo/3D/automeasure/latex/figures/japan/stats_T06.pdf
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Modified: seismo/3D/automeasure/latex/figures_paper.tex
===================================================================
--- seismo/3D/automeasure/latex/figures_paper.tex 2008-05-01 04:11:44 UTC (rev 11879)
+++ seismo/3D/automeasure/latex/figures_paper.tex 2008-05-01 14:21:15 UTC (rev 11880)
@@ -74,7 +74,7 @@
$T_{0,1}$ & 50, 150 & 24, 120 & 6, 30 & 6, 40 & 2, 40 \\
$r_{P,A}$ & 3.5, 3.0 & 3.5, 3.0 & 3.5, 3.0 & 3.5, 3.0 & 3.5, 2.5 \\
$r_0$ & 2.5 & 1.5 & 3.0 & 2.5 & 4.0 \\
-$w_E$ & 0.08 & 0.12 & 0.12 & 0.22 & 0.07 \\
+$w_E$ & 0.08 & 0.10 & 0.12 & 0.22 & 0.07 \\
$CC_0$ & 0.85 & 0.70 & 0.70 & 0.74 & 0.85 \\
$\Delta\tau_0$ & 15 & 12.0 & 3.0 & 3.0 & 2.0 \\
$\Delta\ln{A}_0$& 1.0 & 1.0 & 1.0 & 1.5 & 1.0 \\
@@ -262,9 +262,9 @@
(a)~Event and station map: event is 051502B indicated by the beach ball with the
CMT focal mechanism, station ERM is marked as red triangles and all the other stations
which recorded this event are marked by grey triangles.
-(b)~Results for station ERM for the period range \trange{6}{30}.
+(b)~Results for station ERM for the period range \trange{24}{120}.
Vertical (Z), radial (R), and transverse (T) records of data (black, left column) and synthetics (red, left column), as well as the STA/LTA records (right column) used to produce the window picks.
-(c)~Results for station ERM for the period range \trange{24}{120}.
+(c)~Results for station ERM for the period range \trange{6}{30}.
}
\end{figure}
\clearpage
@@ -278,25 +278,25 @@
(a)~Event and station map: event is 091502B indicated by the beach ball with the
CMT focal mechanism, station KIS is marked as red triangles and all the other stations
which recorded this event are marked by grey triangles.
-(b)~Results for station KIS for the period range \trange{6}{30}.
+(b)~Results for station KIS for the period range \trange{24}{120}.
Vertical (Z), radial (R), and transverse (T) records of data (black, left column) and synthetics (red, left column), as well as the STA/LTA records (right column) used to produce the window picks.
-(c)~Results for station KIS for the period range \trange{24}{120}.
+(c)~Results for station KIS for the period range \trange{6}{30}.
}
\end{figure}
\clearpage
\begin{figure}
%\center
-\includegraphics[width=5.7in]{figures/japan/SHR_BO_20051121536A}
-\caption{\label{fg:SHR_BO_20051121536A}
+\includegraphics[width=5.7in]{figures/japan/SHR_BO_200511211536A}
+\caption{\label{fg:SHR_BO_200511211536A}
Window selection results for event 20051121536A from Table~\ref{tb:events} recorded at station SHR ($44.06$\degN, $144.99$\degE, $\Delta=17.47$\deg).
(a)~Event and station map: event 20051121536A is indicated by the beach ball with the
CMT focal mechanism, station SHR is marked as red triangles and all the other stations
which recorded this event are marked by grey triangles.
-(b)~Results for station SHR for the period range \trange{6}{30}.
+(b)~Results for station SHR for the period range \trange{24}{120}.
Vertical (Z), radial (R), and transverse (T) records of data (black, left column) and synthetics (red, left column), as well as the STA/LTA records (right column) used to produce the window picks.
-(c)~Results for station SHR for the period range \trange{24}{120}.
-Note that corresponding high-frequency band-passed filtered version (b) has shorter record length (600~s).
+(c)~Results for station SHR for the period range \trange{6}{30}.
+Note that corresponding low-frequency band-passed filtered version (b) has longer record length (800~s).
}
\end{figure}
@@ -324,7 +324,7 @@
\clearpage
\begin{figure}
%\center
-\includegraphics[width=6in]{figures/japan/T06_rs-crop}
+\includegraphics[width=6in]{figures/japan/stats_T06}
\caption{\label{fg:T06_rs}
Summary statistics of windowing results for events 051502B, 200511211536A and 091502B in Table~\ref{tb:events}, for the period range \trange{6}{30}.
}
Modified: seismo/3D/automeasure/latex/flexwin_paper.pdf
===================================================================
(Binary files differ)
Modified: seismo/3D/automeasure/latex/results.tex
===================================================================
--- seismo/3D/automeasure/latex/results.tex 2008-05-01 04:11:44 UTC (rev 11879)
+++ seismo/3D/automeasure/latex/results.tex 2008-05-01 14:21:15 UTC (rev 11880)
@@ -174,7 +174,7 @@
\trange{6}{30}, using all the records, and \trange{24}{120}, using the broadband records only.
Figures~\ref{fg:ERM_II_051502B}, \ref{fg:KIS_BO_091502B} and
-\ref{fg:SHR_BO_20051121536A} show windowing results for three events at
+\ref{fg:SHR_BO_200511211536A} show windowing results for three events at
different depths (Table~\ref{tb:events}):
051502B, 22.4~km deep, Taiwan; 091502B, 589.4~km deep,
Northeastern China; 200511211536A, 155~km deep, Kyushu, Japan.
@@ -190,13 +190,19 @@
to predict the short-period surface waves.
For the shallow event beneath Taiwan
-recorded at station ERM (Figure~\ref{fg:ERM_II_051502B}), the synthetics are a poorer fit to the data in the short-period range (\trange{6}{30}), than in the long-period range (\trange{24}{120}).
-For \trange{6}{30}, the synthetics capture only the major $P$-wave
-arrival on the vertical component (Figure~\ref{fg:ERM_II_051502B}b), which is also the only arrival picked by the windowing algorithm.
-For \trange{24}{120}, the observed and synthetic seismograms are similar in shape, and the algorithm picks not only the long-period $P$- and $S$-wave arrivals, but also the
-Rayleigh-wave on the vertical and radial components and the Love-wave
-on the transverse component (Figure~\ref{fg:ERM_II_051502B}c). As the depth sensitivity of surface waves
-is frequency dependent, the good fits at long periods (\trange{24}{120})
+recorded at station ERM (Figure~\ref{fg:ERM_II_051502B}),
+the observed and synthetic seismograms in
+\trange{24}{120} are similar in shape, and the algorithm picks not only the long-period $P$- and $S$-wave arrivals, but also the
+Rayleigh-wave on the vertical and radial components and the
+Love-wave on the transverse component
+(Figure~\ref{fg:ERM_II_051502B}b). In the short-period range
+(\trange{6}{30}), the synthetics are a poorer
+fit to the data than in the long-period range (\trange{24}{120}).
+The synthetics capture only the major $P$-wave
+arrival on the vertical component (Figure~\ref{fg:ERM_II_051502B}c), which is also the only arrival picked by the windowing algorithm.
+As the depth sensitivity of surface waves
+is frequency dependent, the good fits at long periods
+(\trange{24}{120})
indicate that the starting model structure below the Moho
is adequate for tomography.
@@ -206,17 +212,17 @@
the seismograms from this event are
very simple, and contain only two major body-wave arrivals ($P$ and $S$). The windowing algorithm's similarity criterion clearly comes into play here, causing it
not to pick the short-period $S$ arrival
-on the vertical component (Figure~\ref{fg:KIS_BO_091502B}b) as
+on the vertical component (Figure~\ref{fg:KIS_BO_091502B}c) as
the distorted $S$-wave waveform of the data is quite
different from the Gaussian shaped synthetics.
The long-period $S$-wave arrival on the
same component is selected due to higher data-synthetic
-waveform similarity (Figure~\ref{fg:KIS_BO_091502B}c).
+waveform similarity (Figure~\ref{fg:KIS_BO_091502B}b).
-The records of the intermediate-depth event (20051121536A) recorded by station SHR (Figure~\ref{fg:SHR_BO_20051121536A}) contain more seismic phases
+The records of the intermediate-depth event (200511211536A) recorded by station SHR (Figure~\ref{fg:SHR_BO_200511211536A}) contain more seismic phases
than the previous two examples.
On the vertical component of the short-period
- seismogram (Figure~\ref{fg:SHR_BO_20051121536A}b), the $P$-wave
+ seismogram (Figure~\ref{fg:SHR_BO_200511211536A}c), the $P$-wave
arrives at $\sim$230~s, immediately followed by $pP$ and $sPn$,
and the $S$-wave arrives at $\sim$420~s, followed by $sS$ and
$PcP$. The windowing algorithm selects separate windows for the $P$, $S$ and $sS$ arrivals
@@ -226,8 +232,7 @@
intermediate-depth event are not as clearly defined as those of the shallow event (Figure~\ref{fg:ERM_II_051502B}).
Summary plots of window picks for event 200511211536A in the two
-period ranges \trange{6}{30} and \trange{24}{120} are shown in Figures
-Figure~\ref{fg:200511211536A_T06_rs} and~\ref{fg:200511211536A_T24_rs} respectively.
+period ranges \trange{6}{30} and \trange{24}{120} are shown in Figures~\ref{fg:200511211536A_T06_rs} and~\ref{fg:200511211536A_T24_rs} respectively.
On the short-period window record sections, the windows picked by the algorithm form two main branches
that correspond to $P$ and $S$ arrivals. Sparse $P$ arrivals are visible even on the transverse component (Figure~\ref{fg:200511211536A_T06_rs}g).
The number and width of windows for each trace varies with epicentral distance.
More information about the cig-commits
mailing list