Finite-frequency tomography
|
|
- Alvin Griffith
- 5 years ago
- Views:
Transcription
1 Finite-frequency tomography Lapo Boschi October 14, 2009 Application of Born theory (scattering theory) to the ray-theory solution (the forward problem) In earlier lectures we have seen how the tomographic inverse problem can be formulated on the basis of ray theory, i.e. the approximation that ω 1. This approximation leads to the introduction of the eikonal equation, on whose basis seismic ray paths are defined and travel-time or phase anomalies are written as integrals of Earth s structure along the ray-path only: in the ray-theory approximation, sensitivity kernels associated to such measurements are zero everywhere except along the ray path. The requirement that ω 1 limits the quality of our modeled seismic waves (whose frequency, however high, is finite), and therefore the quality of tomographic maps based on ray theory. The limit it imposes will be more severe as frequency decreases and we move away from a regime of very high frequency. How can the resolution limit caused by this flaw in the theory be quantified? as frequency decreases, and the concept of ray path becomes meaningless, what effect will Earth structure away from the ray path have on the seismograms, on the time and phase anomaly that we observe? how should we change our formulation of the inverse problem to account for non-ray-theoretical phenomena? One way to answer these questions is to solve numerically the equations of motion of the Earth, with no approximations. But this is very expensive and practically unaffordable. Another proposed approach involves perturbative theory, or Born theory, that physicists are already familiar with. To see how Born theory works in global seismology, let us once again go back to the Earth s equation of motion, written as a differential equation in displacement, which we dubbed equation (*) in the first lecture of this course. Let L be an operator such that eq. (*) can be simply written Lu = 0. (1) In earlier lectures we assumed that the effect of the earthquake excitation, necessary to have a nontrivial solution to (1), could be prescribed in the form of an initial/boundary condition. It is now more convenient to write it as a body force equivalent f, replacing (1) with Lu = f. (2) If f is impulsive, eq. (2) is called Green s problem ; its solution is denoted G and called Green s function 1. It is useful to write G as a function of the source and receiver positions, 1 The concept of Green s problem is nicely illustrated, for example, by Dahlen and Tromp, Theoretical 1
2 which we shall denote r S and r R, respectively; then G = G(r R, r S ). The complete Green s function here is a 3 3 tensor, and eq. (2) has to be solved for three different impulsive forcing terms f, each oriented in the direction of one of the axes. It can be proved that, once the Green s problem is solved, the displacement field u associated with any seismic source can be derived, by a simple convolution of a function describing the source in time and space, and the Green s function G. We can use the ray-theory approach to solve the Green s problem and find a Green s tensor G 0, valid at high frequencies. Let us then introduce a small perturbation δµ(r), δλ(r), δρ(r), in the parameters describing Earth s structure. µ, λ and ρ, are replaced by µ + δµ, λ+δλ and ρ+δρ in the analytical expression for L; after linearization, a perturbation operator δl can be defined such that the perturbed Green s problem can be written (L + δl)(g 0 + δg) = f, (3) where everything is known apart from the perturbation δg to the Green s function. From eq. (3), neglecting second order terms, LG 0 + δlg 0 + LδG = f, (4) and since LG 0 = f, we are left with LδG = δlg 0, (5) where the right hand side is known. Having already solved (in the ray-theory approximation) the Green s problem associated with the operator L, we are able to find δg by a simple convolution of G 0 with the new forcing term δlg 0. In the frequency domain, δg = G 0 (r R, x) δlg 0 (x, r S )d 3 x, (6) V with x denoting the integration variable. The response u of the perturbed Earth model (µ+δµ, etc.) at any location r, to any earthquake with arbitrary hypocenter r S, can now be found convolving G(r, r S ) + δg(r, r S ) with the appropriate source function. Equation (6) is often interpreted as follows: a wave travels, according to ray theory, from r S to x. Once it is hit, the point x becomes a secundary source and another wave travels from x to r R. The cumulative effect of all possible secundary sources is integrated. The importance of a secundary source (in principle, any point in the Earth can function as secundary source) depends on the properties of the operator δl. If δλ(x) = δµ(x) = δρ(x) = 0, then also δlg 0 (x, r S ) = 0. If δlg 0 (x, r S ) 0, we say that the point x acts as a scatterer for the incident wave. Born theory is often referred to as scattering theory. Born theory holds so long as perturbations δµ, δλ, δρ are small. Then, the perturbed Green s tensor G 0 + δg provides a higher order of accuracy than the ray-theory one, G 0. Scattering and tomography (the inverse problem) After replacing δlg 0 with its explicit expression, eq. (6) can be rewritten in terms of a scattering tensor S, δg = V G 0 (r R, x) S(x) G 0 (x, r S )d 3 x. (7) Global Seismology, Princeton Univ. Press 1998, section 4.1.7; a more lengthy treatment is given by Aki and Richards, Quantitative Seismology, chapters 2 and 4.
3 After some algebra, tensors S P, S S and S ρ are defined so that the scattering tensor can in turn be written S = S P ( δvp v P ) ( ) ( ) δvs δρ + S S + S ρ. (8) v S ρ Replacing (8) into (7), an expression for δg in terms of relative perturbations to compressional and shear velocities and density throughout the Earth is found. This is the main ingredient we need to set up a linear inverse problem. The next step, necessary to derive from (7) and (8) an equation relating seismic observable to Earth heterogeneity, is to write a quantity that we can observe directly from a seismogram in terms of δg. Let Γ(τ) denote the cross-correlation of a reference seismogram u 0 and a perturbed seismogram u 1 (the treatment that follows is going to hold independently for any component of u as defined above), Γ(τ) = t2 t 1 u 0 (t τ)u 1 (t)dt, (9) with the time window (t 1, t 2 ) chosen to isolate the phase of interest. It is reasonable to establish that the delay time δt between reference and perturbed phase equals the value of τ for which Γ(τ) is maximum. After expanding Γ(τ) in a Taylor series around τ = 0, and equating to zero the derivative of the Taylor series (up to second order) with respect to τ, an expression for the delay time in terms of the perturbed seismogram, and therefore δg, is found 2. Eqs. (7) and (8), with explicit expressions for S P, S S, S ρ, are plugged into the latter expression, and after some algebra kernels K P, K S and K ρ can be defined such that 3 δt = V [ ( ) ( ) ( )] δvp δvs δρ K P + K S K ρ d 3 x, (10) v P v S ρ and heteroengeities in Earth structure are now directly related to the delay time, i.e. something that we can pick from a seismogram. The functions K P (x), K S (x) and K ρ (x) have the peculiar banana-doughnut shape (with zero value on the infinite-frequency ray path) long discussed by Dahlen, Nolet and their co-workers at Princeton, and illustrated here in figures 1 and 2. Eq. (10) should be compared with eq. (2) of my lecture notes on body-wave tomography: they relate the same quantities, but the earlier equation rests on pure ray-theory approximation, and the integral there is to be performed only along the ray path; eq. (10) here should represent an improvement to that approximation. Some authors 4 have suggested that the application of finite-frequency tomography, as described here, to global body wave databases helps to enhance model resolution, imaging narrow features like plumes that have been invisible to traditional tomography (figure 3). However, the significance of such improvements is under debate (figure 4). Surface wave scattering The above treatment is naturally applied to the case of body waves, where it is possible to define travel time, and pick it on a seismogram. Surface waves require a separate formulation. The principles, however, are the same, and the results are similar. 2 Dahlen, Hung and Nolet, GJI 2000, vol. 141 page 157, section The following equation is equivalent to eq. (77) of Dahlen and co-authors. 4 e.g., Montelli et al., Science, 303 pages 338ff, 2004.
4 Figure 1: Finite-frequency kernels for a uniform background medium. From Hung, Dahlen and Nolet, Fréchet kernels for finite-frequency traveltimes II. Examples, Geophysical Journal International, 141, pages 175ff, 2000.
5 Figure 2: Finite-frequency kernels for a spherically symmetric reference Earth model. From Hung et al., 2000.
6 Figure 3: P-velocity in the mantle from scattering theory: ascending plumes. From Montelli et al., 2004 (supplementary online material).
7 Figure 4: P-velocity at several depths in the mantle, derived tomographically from ray theory (left) and the finite-frequency approach (right: model of figure 3). From Montelli et al., 2004.
8 Figure 5: Finite-frequency kernels relating the phase (measured from the seismogram) and the phase velocity of 150 s Love waves, for two different source-station geometries. From Boschi, Geophys. J. Int., 167, pages , 2006.
9 The problem can be collapsed to two dimensions only, like in the previous lecture (surface wave tomography). A (θ, φ)-dependent (not r-dependent) kernel relating δt(ω) and the corresponding phase velocity anomaly δc(θ, φ) is found at each frequency ω; let us call it K L (θ, φ; ω) for Love waves, or K R for Rayleigh waves; then eq. (9) from the previous lecture (surface wave tomography) is replaced by δt L (ω) = ω ray path K L(θ, φ; ω)δc L (θ, φ; ω)ds (11) (and an analogous expression for Rayleigh waves), and the inverse problem is set up exactly as in earlier lectures. Looking at K L (θ, φ; ω), K R (θ, φ; ω) we find them to be nonzero over very large portions of the Earth s surface, if compared to the body wave kernels, K P, etc. As surface waves are waves of lower frequency than body waves, it is to be expected that, in their case, limits of the high-frequency (ray theory) approximation be more evident, and scattering effects more relevant. Alternatively, one can leave explicit the dependence of surface wave phase on Earth structure at varying depth; an expression analogous to equation (17) of lecture 3 follows, with the kernels K SH, etc., now accounting for scattering effects as well. The former approach has been followed, for example, by R. Snieder and other authors from the Utrecht group 5, and likewise in the works of the ETH group illustrated in figures 5 and 6. In a more recent article 6, the Princeton group have worked out the latter, more general, approach. They show that known 2-D equations for phase velocity maps can be derived as a particular case of the 3-D formulation. Note that their unperturbed solution (that I have always called ray-theoretical) is now defined as a linear combination of normal modes. Hence, it does not involve ray paths. However, the requirement that heterogeneities be smooth is still needed. The rest of the Princeton treatment is substantially analogous to what I have outlined in this lecture, though much more detailed. 5 Snieder and Nolet, J. geophys. Res. 61 pages 55ff, Spetzler, Trampert and Snieder, Geophys. J. Int., 149 pages 755ff, Zhou, Dahlen and Nolet, Geophys. J. Int., 158, pages 142ff, 2004.
10 Figure 6: The 2D inverse problem: phase-velocity maps obtained from inversions based on numerical finite-frequency (top) and ray theory (bottom), from Rayleigh-wave measurements at 40 s, 75 s, 150 s and 250 s. From Peter et al., Geophys. Res. Lett., 35, L16315, Significant differences between ray-theory and finite-frequency results emerge at long periods.
Global surface-wave tomography
Global surface-wave tomography Lapo Boschi (lapo@erdw.ethz.ch) October 7, 2009 Love and Rayleigh waves, radial anisotropy Whenever an elastic medium is bounded by a free surface, coherent waves arise that
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2) 85, 385 42 doi:./j.365-246x.2.4945.x Tomographic errors from wave front healing: more than just afastbias Alison E. Malcolm and Jeannot Trampert 2
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2013) 194, 1042 1049 Advance Access publication 2013 May 3 doi: 10.1093/gji/ggt144 Finite-frequency sensitivity kernels for two-station surface wave
More informationTomographic Errors from Wavefront Healing: more than
Geophys. J. Int. (2) 42, Tomographic Errors from Wavefront Healing: more than just a fast bias Alison E. Malcolm, Jeannot Trampert 2 Department of Earth, Atmospheric and Planetary Sciences, Massachusetts
More informationFinite-frequency sensitivity of body waves to anisotropy based upon adjoint methods
Geophys. J. Int. (2007) 171, 368 389 doi: 10.1111/j.1365-246X.2007.03528.x Finite-frequency sensitivity of body waves to anisotropy based upon adjoint methods Anne Sieminski, 1 Qinya Liu, 2 Jeannot Trampert
More informationThe Basic Properties of Surface Waves
The Basic Properties of Surface Waves Lapo Boschi lapo@erdw.ethz.ch April 24, 202 Love and Rayleigh Waves Whenever an elastic medium is bounded by a free surface, coherent waves arise that travel along
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2010) doi: 10.1111/j.1365-246X.2010.04643.x Empirically determined finite frequency sensitivity kernels for surface waves Fan-Chi Lin and Michael H.
More informationTomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima
1 2 3 Auxiliary materials for Tomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima nuclear power plant area 4 5 6 7 8 9 Ping Tong 1,2, Dapeng Zhao 1 and Dinghui Yang 2 [1] {Department of Geophysics,
More informationOn the relevance of Born theory in global seismic tomography
GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:10.1029/, On the relevance of Born theory in global seismic tomography Lapo Boschi ETH Zürich, Switzerland. Thorsten W. Becker University of Southern California,
More informationDatabases of surface wave dispersion
submitted to Geophys. J. Int. Databases of surface wave dispersion Simona Carannante (1,2), Lapo Boschi (3) (1) Dipartimento di Scienze Fisiche, Università Federico II, Naples, Italy (2) Istituto Nazionale
More informationTraveltime sensitivity kernels: Banana-doughnuts or just plain bananas? a
Traveltime sensitivity kernels: Banana-doughnuts or just plain bananas? a a Published in SEP report, 103, 61-68 (2000) James Rickett 1 INTRODUCTION Estimating an accurate velocity function is one of the
More informationStructural sensitivities of finite-frequency seismic waves: a full-wave approach
Geophys. J. Int. (2006) 165, 981 990 doi: 10.1111/j.1365-246X.2006.02993.x Structural sensitivities of finite-frequency seismic waves: a full-wave approach Li Zhao and Thomas H. Jordan Department of Earth
More informationEarthscope Imaging Science & CIG Seismology Workshop
Earthscope Imaging Science & CIG Seismology Introduction to Direct Imaging Methods Alan Levander Department of Earth Science Rice University 1 Two classes of scattered wave imaging systems 1. Incoherent
More informationTraveltimes and amplitudes of seismic waves: a re-assessment
Traveltimes and amplitudes of seismic waves: a re-assessment 1 Guust Nolet, F.A. Dahlen and Raffaella Montelli Department of Geosciences, Princeton University, Princeton NJ. Abstract In this paper we give
More informationSeismic tomography, adjoint methods, time reversal, and banana-doughnut kernels
Seismic tomography, adjoint methods, time reversal, and banana-doughnut kernels Geophysical Journal International, 25, v. 16, p. 195 216 Jeroen Tromp, Carl Tape, and Qinya Liu Seismological Laboratory,
More informationSurface Wave Tomography
Surface Wave Tomography Alexandra Mauerberger gassner@gfz-potsdam.de GEO-DEEP9300 course at CEED, University of Oslo Nov. 6 - Nov. 10 2017 Outline Recap General Methods Techniques to Derive Tomography
More informationDownloaded 11/08/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Velocity tomography based on turning waves and finite-frequency sensitivity kernels Xiao-Bi Xie*, University of California at Santa Cruz, Jan Pajchel, Statoil, Norway Summary We present a velocity tomography
More informationSeismic tomography: Art or science? Frederik J Simons Princeton University
Seismic tomography: Art or science? Frederik J Simons Princeton University What s inside the Earth? www.usgs.gov www.glencoe.com Dalton, Nature 2003 Only seismic waves have actually been there, done that
More informationSeismogram Interpretation. Seismogram Interpretation
Travel times in the Earth Ray paths, phases and their name Wavefields in the Earth: SH waves, P-SV waves Seismic Tomography Receiver Functions Seismogram Example Long-period transverse displacement for
More informationSeismic tomography: Art or science?
Seismic tomography: Art or science? Frederik J Simons Princeton University What s inside the Earth? www.usgs.gov www.glencoe.com Dalton, Nature 2003 Only seismic waves have actually been there, done that
More informationANEWJOINTP AND S VELOCITY MODEL OF THE MANTLE PARAMETERIZED IN CUBIC B-SPLINES
ANEWJOINTP AND S VELOCITY MODEL OF THE MANTLE PARAMETERIZED IN CUBIC B-SPLINES Michael Antolik, Göran Ekström, Adam M. Dziewonski, Yu J. Gu, Jian-feng Pan, and Lapo Boschi Department of Earth and Planetary
More informationGlobal surface wave diffraction tomography
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. B12, 2335, doi:10.1029/2002jb001777, 2002 Global surface wave diffraction tomography Michael H. Ritzwoller, Nikolai M. Shapiro, Mikhail P. Barmin, and Anatoli
More informationSurface wave tomography for azimuthal anisotropy in a strongly reduced parameter space
Geophys. J. Int. (2008) doi: 10.1111/j.1365-246X.2008.03833.x Surface wave tomography for azimuthal anisotropy in a strongly reduced parameter space Mark P. Panning 1 and Guust Nolet 1,2 1 Department of
More informationNon-linear crustal corrections in high-resolution regional waveform seismic tomography
Geophys. J. Int. (27) 17, 46 467 doi: 1.1111/j.1365-246X.27.3399.x Non-linear crustal corrections in high-resolution regional waveform seismic tomography Federica Marone and Barbara Romanowicz Berkeley
More informationPEAT SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II
PEAT8002 - SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II Nick Rawlinson Research School of Earth Sciences Australian National University Waveform modelling P-wave first-motions
More information27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
CRUSTAL AND MANTLE STRUCTURE BENEATH EASTERN EURASIA FROM FINITE FREQUENCY SEISMIC TOMOGRAPHY (FFST) Ting Yang 2, Yang Shen 2, and Xiaoping Yang 1 Science Applications International Corporation 1 and University
More informationShear wave splitting in three-dimensional anisotropic media
Geophys. J. Int. (24) 59, 7 72 doi:./j.365-246x.24.2432.x Shear wave splitting in three-dimensional anisotropic media Sébastien Chevrot, Noémie Favier and Dimitri Komatitsch 2 Laboratoire de Dynamique
More informationChapter 6. Secondary source theories
Chapter 6 Secondary source theories So far we have concentrated on synthetic seismogram calculations for the case of one-dimensional Earth models in which the velocity varies only as a function of depth.
More informationMinor-arc and major-arc global surface wave diffraction tomography
Physics of the Earth and Planetary Interiors 149 (2005) 205 223 Minor-arc and major-arc global surface wave diffraction tomography Anatoli L. Levshin a,, Michael P. Barmin a, Michael H. Ritzwoller a, Jeannot
More informationGlobal anisotropic phase velocity maps for higher mode Love and Rayleigh waves
Geophys. J. Int. (8) 7, 6 3 doi:./j.365-46x.7.3685.x Global anisotropic phase velocity maps for higher mode Love and Rayleigh waves K. Visser, J. Trampert and B. L. N. Kennett Department of Earth Sciences,
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
FINITE-FREQUENCY SEISMIC TOMOGRAPHY OF BODY WAVES AND SURFACE WAVES FROM AMBIENT SEISMIC NOISE: CRUSTAL AND MANTLE STRUCTURE BENEATH EASTERN EURASIA Yong Ren 2, Wei Zhang 2, Ting Yang 3, Yang Shen 2,and
More informationSurface wave focusing effects: Numerical modeling and statistical observations
Physics of the Earth and Planetary Interiors 155 (2006) 191 200 Surface wave focusing effects: Numerical modeling and statistical observations B.G. Bukchin a,b, T.B. Yanovskaya c, J.-P. Montagner b,, A.Z.
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. 2014) 196, 1827 1843 Advance Access publication 2014 January 2 doi: 10.1093/gji/ggt496 Non-linear partial derivative and its De Wolf approximation for
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2) 82, 25 42 doi:./j.365-246x.2.4667.x Frequency-dependent effects on global S-wave traveltimes: wavefront-healing, scattering and attenuation Christophe
More informationWide-band coupling of Earth s normal modes due to anisotropic inner core structure
Geophys. J. Int. (28) 174, 919 929 doi: 1.1111/j.1365-246X.28.3824.x Wide-band coupling of Earth s normal modes due to anisotropic inner core structure J. C. E. Irving, A. Deuss and J. Andrews University
More informationSeismology and Seismic Imaging
Seismology and Seismic Imaging 4. Ray theory N. Rawlinson Research School of Earth Sciences, ANU Seismology lecture course p.1/23 The ray approximation Here, we consider the problem of how body waves (P
More informationExistence of finite rigidity layer at the base of the Earth s liquid outer core inferred from anomalous splitting of normal modes
LETTER Earth Planets Space, 54, 67 7, 22 Existence of finite rigidity layer at the base of the Earth s liquid outer core inferred from anomalous splitting of normal modes Seiji Tsuboi,2 and Masanori Saito
More informationSome aspects of seismic tomography
Some aspects of seismic tomography Peter Shearer IGPP/SIO/U.C. San Diego September 7, 2009 Earthquake Research Institute Part 1: Global Tomography P velocity perturbations 200 km 1000 km 2700 km MIT 2006
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (215) 23, 232 2 GJI Seismology doi: 1.193/gji/ggv13 Finite-frequency sensitivity kernels of seismic waves to fault zone structures A.A. Allam, 1 C. Tape
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2010) 180, 1187 1199 doi: 10.1111/j.1365-246X.2009.04468.x Assessment of tomographic mantle models using spectral element seismograms Ebru Bozdağ and
More informationTomography of core-mantle boundary and lowermost mantle coupled by geodynamics: joint models of shear and compressional velocity
1 2 3 4 5 6 7 8 9 10 11 12 Tomography of core-mantle boundary and lowermost mantle coupled by geodynamics: joint models of shear and compressional velocity Gaia Soldati 1, Lapo Boschi 2,3, Steve Della
More informationSupporting Information for An automatically updated S-wave model of the upper mantle and the depth extent of azimuthal anisotropy
GEOPHYSICAL RESEARCH LETTERS Supporting Information for An automatically updated S-wave model of the upper mantle and the depth extent of azimuthal anisotropy Eric Debayle 1, Fabien Dubuffet 1 and Stéphanie
More informationWhat Can Seismology Say About Hot Spots?
What Can Seismology Say About Hot Spots? Bruce R. Julian U. S. Geological Survey, Menlo Park, CA 94025 USA julian@usgs.gov G. R. Foulger Dept. of Geological Sciences, Univ. of Durham, Durham DH1 3LE, U.K.
More informationSeismic ray path variations in a 3D global velocity model
Physics of the Earth and Planetary Interiors 141 (2004) 153 166 Seismic ray path variations in a 3D global velocity model Dapeng Zhao, Jianshe Lei Geodynamics Research Center, Ehime University, Matsuyama
More informationBorehole Geophysics. Acoustic logging measurements
Acoustic logging measurements - Review of basic physics background - Concept of P- and S-wave measurements and logging tools - Tube waves - Seismic imaging - Synthetic seismograms - Field application examples
More informationTheoretical Seismology (GEOP 523) Seismic Sources
Theoretical Seismology GEOP 523 Seismic Sources Rick Aster April 25, 2011 Up to this point, we have been mostly concerned with homogeneous solutions to the wave equation, gx, t 2 g t 2 c2 2 g = 0. 1 In
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2011) doi: 10.1111/j.1365-246X.2011.05100.x Apparent anisotropy in inhomogeneous isotropic media Fan-Chi Lin and Michael H. Ritzwoller Center for Imaging
More informationValidation of first-order diffraction theory for the traveltimes and amplitudes of propagating waves
GEOPHYSICS, VO. 71, NO. 6 NOVEMBER-DECEMBER 26; P. T167 T177, 13 FIGS. 1.119/1.2358412 Validation of first-order diffraction theory for the traveltimes and amplitudes of propagating waves Jeroen Jocker
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) 188, 293 300 doi: 10.1111/j.1365-246X.2011.05256.x An analysis of SS precursors using spectral-element method seismograms L. Bai, Y. Zhang and
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (218) 213, 851 867 Advance Access publication 218 January 23 GJI Seismology doi: 1.193/gji/ggy2 Effects of induced stress on seismic forward modelling
More informationObservation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story
Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear
More informationFinite-frequency resolution limits of wave path traveltime tomography for smoothly varying velocity models
Geophys. J. Int. (3) 5, 669 676 Finite-frequency resolution limits of wave path traveltime tomography for smoothly varying velocity models Jianming Sheng and Gerard T. Schuster Department of Geology and
More informationWe briefly discuss two examples for solving wave propagation type problems with finite differences, the acoustic and the seismic problem.
Excerpt from GEOL557 Numerical Modeling of Earth Systems by Becker and Kaus 2016 1 Wave propagation Figure 1: Finite difference discretization of the 2D acoustic problem. We briefly discuss two examples
More informationAVO inversion in V (x, z) media
Stanford Exploration Project, Report 97, July 8, 998, pages 75 94 AVO inversion in V (x, z) media Yalei Sun and Wenjie Dong keywords:.5-d Kirchhoff integral, AVO inversion, fluid-line section ABSTRACT
More informationSURFACE WAVE GROUP VELOCITY MEASUREMENTS ACROSS EURASIA
SURFACE WAVE GROUP VELOCITY MEASUREMENTS ACROSS EURASIA A. L. Levshin, M. H. Ritzwoller, and L. I. Ratnikova Department of Physics, University of Colorado at Boulder -Contract Number F49620-95-1-0139 Sponsored
More information2 RITZWOLLER ET AL.: GLOBAL SURFACE WAVE DIFFRACTION TOMOGRAPHY Abstract. Our purpose is to determine the eect of replacing geometrical ray-theory in
JOURNAL OF GEOPHYSICAL RESEARCH, VOL., NO., PAGES 1??, Global surface wave diraction tomography Michael H. Ritzwoller, Nikolai M. Shapiro, Mikhail P. Barmin, Anatoli L. Levshin Center for Imaging the Earth's
More informationThe Local Discrete Wave-number Method to Simulate Wave Propagation in Irregular Layer
October 1-17, 8, Beijing, hina The Local Discrete Wave-number Method to Simulate Wave Propagation in Irregular Layer Hong Zhou 1 Xiaofei hen 1 Associate Professor, Dept. of Structural Engineering, Institute
More information620 C 2003 RAS. Geophys. J. Int. (2003) 152,
Geophys. J. Int. (2003) 152, 620 632 Surface wave tomography of North America and the Caribbean using global and regional broad-band networks: phase velocity maps and limitations of ray theory Stéphanie
More informationComposite memory variables for viscoelastic synthetic seismograms
Geophys. J. Int. (1995) 121,634-639 RESEARCH NOTE Composite memory variables for viscoelastic synthetic seismograms Tong Xu George A. McMechan Center for Lithospheric Studies, The University of Texas at
More informationStrain Green s Tensors, Reciprocity, and Their Applications to Seismic Source and Structure Studies
Bulletin of the Seismological Society of America, Vol. 96, No. 5, pp. 1753 1763, October 2006, doi: 10.1785/0120050253 Strain Green s Tensors, Reciprocity, and Their Applications to Seismic Source and
More informationDevelopment and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks
Development and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks Erika Tudisco Research Group: Stephen A. Hall Philippe
More informationObservations of long period Rayleigh wave ellipticity
Geophys. J. Int. (), doi:./j.-x...x Observations of long period Rayleigh wave ellipticity Ana M. G. Ferreira and John H. Woodhouse Department of Earth Sciences, University of Oxford, Parks Road, Oxford
More information3D IMAGING OF THE EARTH S MANTLE: FROM SLABS TO PLUMES
3D IMAGING OF THE EARTH S MANTLE: FROM SLABS TO PLUMES Barbara Romanowicz Department of Earth and Planetary Science, U. C. Berkeley Dr. Barbara Romanowicz, UC Berkeley (KITP Colloquium 9/11/02) 1 Cartoon
More informationRELATIVE SEISMIC ATTENUATION ESTIMATION
RELATIVE SEISMIC ATTENUATION ESTIMATION ADAM T. RINGLER 1. INTRODUCTION The end goal of any solid Earth study is to better constrain physical parameters of the Earth (e.g. viscosity, temperature, and mineral
More informationGlobal upper-mantle tomography with the automated multimode inversion of surface and S-wave forms
Geophys. J. Int. (28) 173, 55 518 doi: 1.1111/j.1365-246X.28.3721.x Global upper-mantle tomography with the automated multimode inversion of surface and S-wave forms Sergei Lebedev 1 and Rob D. van der
More informationThe Math, Science and Computation of Hydraulic Fracturing
The Math, Science and Computation of Hydraulic Fracturing 03/21/2013 LLNL-PRES-XXXXXX This work was performed under the auspices of the U.S. Department of Energy by under contract DE-AC52-07NA27344. Lawrence
More informationContinent-sized anomalous zones with low seismic velocity at the base of Earth s mantle
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO2733 Continent-sized anomalous zones with low seismic velocity at the base of Earth s mantle Edward J. Garnero 1, Allen K. McNamara 1, and Sang-Heon D. Shim 1
More informationDatabases of surface wave dispersion
ANNALS OF GEOPHYSICS, VOL. 48, N. 6, December 2005 Databases of surface wave dispersion Simona Carannante ( 1 ) and Lapo Boschi ( 2 ) ( 1 ) Istituto Nazionale di Geofisica e Vulcanologia, Sede di Bologna,
More informationMantle plume tomography
Chemical Geology 241 (2007) 248 263 www.elsevier.com/locate/chemgeo Mantle plume tomography Guust Nolet a, Richard Allen b, Dapeng Zhao c, a Department of Geosciences, Princeton University, Princeton,
More information2.1 Introduction. 16 Chapter2
Chapter 2 Surface wave tomography of North America and the Caribbean using global and regional broadband networks: phase velocity maps and limitations of ray theory Abstract We present phase velocity maps
More informationSupplementary Online Material for. Seismic evidence for a chemically distinct thermochemical reservoir in Earth s deep mantle beneath Hawaii
Supplementary Online Material for Seismic evidence for a chemically distinct thermochemical reservoir in Earth s deep mantle beneath Hawaii Authors: Chunpeng Zhao 1, Edward J. Garnero 1,*, Allen K. McNamara
More informationEstimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 24, 2248, doi:10.1029/2003gl018413, 2003 Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
More informationAPPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY
APPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY Timur TEZEL Supervisor: Takuo SHIBUTANI MEE07169 ABSTRACT In this study I tried to determine the shear wave velocity structure in the crust
More informationboundaries with additional record sections, as emphasized in Fig. S2. The observations at the
Data used to Sample African Anomaly. The great circle paths of the source-receiver combinations used in this study are shown in Fig. S1. The event information is given in Table S1. Abrupt Changes across
More informationNormal modes of the Earth
Journal of Physics: Conference Series Normal modes of the Earth To cite this article: Jean-Paul Montagner and Geneviève Roult 08 J. Phys.: Conf. Ser. 118 0104 View the article online for updates and enhancements.
More informationPhysics of the Earth and Planetary Interiors xxx (2004) xxx xxx. Minor-arc and major-arc global surface wave diffraction tomography
Physics of the Earth and Planetary Interiors xxx (2004) xxx xxx 3 4 Minor-arc and major-arc global surface wave diffraction tomography 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 Abstract
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
TRANSITION ZONE WAVE PROPAGATION: CHARACTERIZING TRAVEL-TIME AND AMPLITUDE INFORMATION Peter M. Shearer and Jesse F. Lawrence University of California San Diego, Institute of Geophysics and Planetary Physics
More informationOn the Earthquake-Source Numerical Implementation in the Seismic Wave Equation
Journal of Earthquake Engineering ISSN: 1363-2469 (Print) 1559-808X (Online) Journal homepage: http://www.tandfonline.com/loi/ueqe20 On the Earthquake-Source Numerical Implementation in the Seismic Wave
More informationVertical coherence in mantle heterogeneity from global seismic data
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl049281, 2011 Vertical coherence in mantle heterogeneity from global seismic data L. Boschi 1,2 and T. W. Becker 3 Received 11 August 2011; revised
More informationThe signature of attenuation and anisotropy on AVO and inversion sensitivities
The signature of attenuation and anisotropy on AVO and inversion sensitivities Shahpoor Moradi and Kristopher Innanen University of Calgary, Department of Geoscience, Calgary, Canada Summary We study the
More informationSOEE3250/5675/5115 Inverse Theory Lecture 10; notes by G. Houseman
SOEE3250/5675/5115 Inverse Theory Lecture 10; notes by G. Houseman Travel-time tomography Examples of regional lithospheric tomography CBP / SCP projects data acquisition: array / sources arrival time
More informationPHASE TIME INVERSION: A SIMPLE METHOD FOR REGIONAL WAVEFORM INVERSION. Charles A. Langston. University of Memphis
PHASE TIME INVERSION: A SIMPLE METHOD FOR REGIONAL WAVEFORM INVERSION Charles A. Langston University of Memphis Sponsored by National Nuclear Security Administration Office of Nonproliferation Research
More informationHigh Resolution Imaging of Fault Zone Properties
Annual Report on 1998-99 Studies, Southern California Earthquake Center High Resolution Imaging of Fault Zone Properties Yehuda Ben-Zion Department of Earth Sciences, University of Southern California
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. 1 189, 138 146 doi: 1./j.1365-46X.1.546.x Calculation of normal mode spectra in laterally heterogeneous earth models using an iterative direct solution
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/326/5949/112/dc1 Supporting Online Material for Global Surface Wave Tomography Using Seismic Hum Kiwamu Nishida,* Jean-Paul Montagner, Hitoshi Kawakatsu *To whom correspondence
More informationGlobal propagation of body waves revealed by cross-correlation analysis of seismic hum
GEOPHYSICAL RESEARCH LETTERS, VOL. 4, 1691 1696, doi:1.12/grl.5269, 213 Global propagation of body waves revealed by cross-correlation analysis of seismic hum K. Nishida 1 Received 1 January 213; revised
More informationPEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity
PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical
More informationData Repository: Seismic and Geodetic Evidence For Extensive, Long-Lived Fault Damage Zones
DR2009082 Data Repository: Seismic and Geodetic Evidence For Extensive, Long-Lived Fault Damage Zones Fault Zone Trapped Wave Data and Methods Fault zone trapped waves observed for 2 shots and 5 local
More informationComputational Seismology: An Introduction
Computational Seismology: An Introduction Aim of lecture: Understand why we need numerical methods to understand our world Learn about various numerical methods (finite differences, pseudospectal methods,
More informationCompensating for attenuation by inverse Q filtering. Carlos A. Montaña Dr. Gary F. Margrave
Compensating for attenuation by inverse Q filtering Carlos A. Montaña Dr. Gary F. Margrave Motivation Assess and compare the different methods of applying inverse Q filter Use Q filter as a reference to
More informationTime-lapse traveltime change of singly scattered acoustic waves
Geophys. J. Int. 2006 65, 485 500 doi: 0./j.365-246X.2006.02856.x Time-lapse traveltime change of singly scattered acoustic waves C. Pacheco and R. Snieder Center for Wave Phenomena, Department of Geophysics,
More informationAcoustic wave-equation-based earthquake location
submitted to Geophys. J. Int. Acoustic wave-equation-based earthquake location Ping Tong 1, Dinghui Yang 2, Qinya Liu 3, Xu Yang 4 and Jerry Harris 1 1 Department of Geophysics, Stanford University, California,
More information1. University of Ottawa 2. Dublin Institute for Advanced Studies 3. University of Texas at Austin
Supplementary Figures: Azimuthal seismic anisotropy in the Earth s upper mantle and the thickness of tectonic plates A.J. Schaeffer 1, S. Lebedev 2 and T.W. Becker 3 Geophysical Journal International July
More informationIntroduction to Seismic Imaging
Introduction to Seismic Imaging Alison Malcolm Department of Earth, Atmospheric and Planetary Sciences MIT August 20, 2010 Outline Introduction Why we image the Earth How data are collected Imaging vs
More informationSeismic Waves and Earthquakes A Mathematical Overview
Seismic Waves and Earthquakes A Mathematical Overview The destruction caused by earthquakes is caused by the seismic waves propagating inside and in particular, near the surface of earth. We describe the
More informationOn sensitivity kernels for wave-equation transmission tomography
Geophys. J. Int. (24) doi: 1.1111/j.1365-246X.24.259.x On sensitivity kernels for wave-equation transmission tomography Maarten V. de Hoop and obert D. van der Hilst Department of Earth, Atmospheric and
More informationA magnitude 7.4 earthquake struck 255 km (158 miles) southwest of Tonga, according to the US Geological Survey, but there were no reports of damage.
A magnitude 7.4 earthquake struck 255 km (158 miles) southwest of Tonga, according to the US Geological Survey, but there were no reports of damage. The earthquake hit at 5:19 am local time and was centered
More informationSensitivity of frequency-dependent traveltimes to laterally heterogeneous, anisotropic Earth structure
Geophys. J. Int. (1998) 133, 683 704 Sensitivity of frequency-dependent traveltimes to laterally heterogeneous, anisotropic Earth structure Li Zhao and Thomas H. Jordan Department of Earth, Atmospheric
More informationCrust and upper mantle P- and S-wave delay times at Eurasian seismic stations
Physics of the Earth and Planetary Interiors 123 (2001) 205 219 Crust and upper mantle P- and S-wave delay times at Eurasian seismic stations E.R. Engdahl, M.H. Ritzwoller Center for Imaging the Earth
More information1. A few words about EarthScope and USArray. 3. Tomography using noise and Aki s method
1. A few words about EarthScope and USArray 2. Surface-wave studies of the crust and mantle 3. Tomography using noise and Aki s method 4. Remarkable images of US crust (and basins)! Unlocking the Secrets
More information