[cig-commits] r14613 - seismo/3D/ADJOINT_TOMO/flexwin_paper/latex

alessia at geodynamics.org alessia at geodynamics.org
Tue Apr 7 08:06:44 PDT 2009


Author: alessia
Date: 2009-04-07 08:06:44 -0700 (Tue, 07 Apr 2009)
New Revision: 14613

Added:
   seismo/3D/ADJOINT_TOMO/flexwin_paper/latex/cig_one_page.tex
Modified:
   seismo/3D/ADJOINT_TOMO/flexwin_paper/latex/AM-allcitations.bib
Log:
Added cig one-pager

Modified: seismo/3D/ADJOINT_TOMO/flexwin_paper/latex/AM-allcitations.bib
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin_paper/latex/AM-allcitations.bib	2009-04-07 05:03:16 UTC (rev 14612)
+++ seismo/3D/ADJOINT_TOMO/flexwin_paper/latex/AM-allcitations.bib	2009-04-07 15:06:44 UTC (rev 14613)
@@ -6208,6 +6208,14 @@
 	Volume = 168,
 	Year = 2007}
 
+ at article{TapeEtal2009,
+	Author = {Tape, C. and Liu, Q. and Maggi, A. and Tromp, J.},
+	Title = {{Adjoint tomography of the southern California crust}},
+	Year = 2009,
+	note = {(submitted)}
+	}
+
+
 @article{TarantolaValette1982,
 	Author = {Tarantola, A. and Valette, B.},
 	Journal = rog,

Added: seismo/3D/ADJOINT_TOMO/flexwin_paper/latex/cig_one_page.tex
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin_paper/latex/cig_one_page.tex	                        (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin_paper/latex/cig_one_page.tex	2009-04-07 15:06:44 UTC (rev 14613)
@@ -0,0 +1,109 @@
+%
+%=====================================================
+%
+%     Carl Tape
+%     cig_one_page.tex
+%     last printed xxx
+% 
+%=====================================================
+
+\documentclass[10pt,titlepage,fleqn]{article}
+
+\usepackage{amsmath}
+\usepackage{amssymb}
+\usepackage{latexsym}
+\usepackage[round]{natbib}
+\usepackage{xspace}
+\usepackage{graphicx}
+%\usepackage{epsfig}
+
+\usepackage{pifont}   % search for \ding
+
+%\usepackage{fancyhdr}
+%\pagestyle{fancy}
+
+%=====================================================
+%       SPACING COMMANDS (Latex Companion, p. 52)
+%=====================================================
+
+\usepackage{setspace}
+
+\renewcommand{\baselinestretch}{1.0}
+
+%\textwidth 460pt
+%\textheight 690pt
+\textwidth 570pt
+\textheight 790pt
+\oddsidemargin 0pt
+\evensidemargin 0pt
+
+% see Latex Companion, p. 85
+\voffset     -70pt
+\topmargin     0pt
+\headsep      20pt
+\headheight   15pt
+\headheight    0pt
+\footskip     30pt
+\hoffset     -50pt
+
+\include{def_base}
+
+\graphicspath{
+  {./figures/fig/}
+}
+
+%=====================================================
+\begin{document}
+%=====================================================
+
+{\bf \large
+\noindent An automated time-window selection algorithm for seismic tomography}
+
+\vspace{0.2cm}
+\noindent Alessia Maggi$^1$, Carl Tape$^2$, Min Chen$^3$, Daniel Chao$^2$, and Jeroen Tromp$^4$
+
+\vspace{0.2cm}
+\noindent $^1$Institut de Physique du Globe, Universit\'{e} de Strasbourg, Strasbourg, France \\ 
+\noindent $^2$Seismological Laboratory, California Institute of Technology, Pasadena, California, USA \\
+\noindent $^3$Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA\\ \noindent $^4$Department of Geosciences, Princeton University, Princeton, New Jersey, USA
+
+\bigskip\noindent
+
+We have developed an open source algorithm {\tt FLEXWIN} \citep{MaggiEtal2009} to automatically select time windows on pairs of observed and synthetic seismograms.  The algorithm was designed specifically to accommodate synthetic seismograms produced from 3D wavefield simulations such as {\tt SPECFEM3D} \citep{KomatitschEtal2002}, which capture complex phases that do not necessarily exist in 1D simulations.  Our motivation is to use the algorithm for iterative tomographic inversions, in which the synthetic seismograms change from one iteration to the next, and contain more and more 3D-specific phases \citep{TapeEtal2007,TapeEtal2009}.  
+
+Our strategy is to select time windows on the synthetic seismogram within
+which the waveform contains a distinct arrival, then require an adequate correspondence between observed and synthetic waveforms within these windows.  This selection paradigm is general, and can be applied to synthetic seismograms regardless of how they have been obtained. It is clear, however, that a synthetic seismogram obtained by 3D propagation (using e.g.~{\tt SPECFEM3D}) through a good 3D Earth model will provide a better fit to the observed seismogram over a greater proportion of its length than will be the case for a more approximate synthetic seismogram.
+
+In order to isolate changes in amplitude or frequency content potentially associated with distinct seismic arrivals, we have chosen a tool used in automated phase detection for earthquake location, the the short-term average / long-term average ratio, and we have adapted it to define time windows around distinct seismic phases (see Figure~\ref{fg:onepage}).
+
+\begin{figure}[h]
+\center \includegraphics[width=6in]{cig_onepager.pdf}
+\caption[Figure : ] 
+{
+Window selection results on the vertical component for a intermediate-depth event in Northern Peru recorded at station ABKT
+($\Delta=127$\deg). (a) Top: observed and synthetic seismograms
+(black and red traces); bottom: STA/LTA timeseries. Windows chosen by the algorithm are shown using light blue shading. The phases contained these windows are: (1) $PP$, (2) $PS + SP$, (3) $SS$, (4) $SSS$, (5) $S5$, (6) $S6$, (7) fundamental mode Rayleigh  wave. (b) Ray paths corresponding to the body-wave phases present in the data windows
+in (a).
+\label{fg:onepage}
+}
+\end{figure}
+
+We have successfully applied our windowing algorithm to diverse seismological scenarios: local and near regional tomography in Southern
+California, regional subduction-zone tomography in Japan, and global tomography \citep{MaggiEtal2009}.
+
+
+%=====================================================
+
+\vspace{0.2cm}
+
+\small
+\begin{spacing}{1.0}
+\bibliographystyle{plainnat}
+\bibliography{AM-allcitations}
+\end{spacing}
+
+
+
+%=====================================================
+\end{document}
+%=====================================================



More information about the CIG-COMMITS mailing list