Correlation based imaging

Size: px
Start display at page:

Download "Correlation based imaging"

Transcription

1 Correlation based imaging George Papanicolaou Stanford University International Conference on Applied Mathematics Heraklion, Crete September 17, 2013 G. Papanicolaou, ACMAC-Crete Correlation based imaging 1/33

2 Imaging with ambient noise Key words in imaging: Passive sensors, coherent imaging, noise source illumination, deterministic or random medium, auxiliary arrays. 1. Recovering information (travel times, the Green s function) from cross correlations of recorded noise signals. 2. Using stationary phase to understand the role of directivity in noise sources. Background velocity estimation. 3. Virtual source imaging with auxiliary arrays. Minimizing the effect of medium inhomogeneities with correlation-based imaging using auxiliary arrays. G. Papanicolaou, ACMAC-Crete Correlation based imaging 2/33

3 Passive sensor imaging with noise cross correlations 1. When there is enough ambient noise, cross correlations can determine wave energy inter-sensor travel times (singular support of time domain Green s function) 2. If sensor locations are known, a background speed of wave propagation can be estimated tomographically, as has been done very successfully with regional seismic data (San Diego group and Grenoble 2005) 3. If background speed is known then sensor locations can be estimated (passive topological self-awareness of sensor networks) 4. How can we detect and possibly image changes in the environment, such as geophysical changes or moving reflectors, with noise cross correlations? G. Papanicolaou, ACMAC-Crete Correlation based imaging 3/33

4 Southern California Seismic Network (a) Map of the 150 stations of the SCSN; (b) Shot record generated by cross correlations over one month of noise measurements; (c) Surface-wave group velocity map Figure from: A. Curtis, P. Gerstoft, H. Sato, R. Snieder, and K. Wapenaar, Seismic Interferometry -Turning Noise to Signal, The Leading Edge, September G. Papanicolaou, ACMAC-Crete Correlation based imaging 4/33

5 European Seismic Network (Alps) Left the seismic sensor network with 3000 selected paths. Right the Rayleigh surface wave velocity map obtained tomographically from the correlation-estimated travel times. Figures from: P. Gouiédard, L. Stehly, F. Brenguier, M. Campillo, Y. Colin de Verdière, E. Larose, L. Margerin, P. Roux, F. J. Sánchez-Sesma, N. M. Shapiro and R. L. Weaver, Cross-correlation of random fields: mathematical approach and applications, Geophysical Prospecting, 2008, 56, G. Papanicolaou, ACMAC-Crete Correlation based imaging 5/33

6 Hellenic Unified Seismic Network From: Ambient Noise Levels in Greece as Recorded at the Hellenic Unified Seismic Network by C. P. Evangelidis and N. S. Melis, Bulletin of the Seismological Society of America, Vol. 102, No. 6, pp , December G. Papanicolaou, ACMAC-Crete Correlation based imaging 6/33

7 The empirical cross correlation and geophysics The empirical cross correlation of the signals recorded at x 1 and x 2 for an integration time T is C T (τ, x 1, x 2 ) = 1 T T 0 u(t, x 1 )u(t + τ, x 2 )dt. (1) Here u(t, x j ) is the signal recordded at x j and T is chosen appropriately (a non-trivial statistical signal processing issue in applications). What is generating the wave energy recorded in the noise signals? Frequency band of microseisms from wave breaking: Hz. Dispersiveness of the underlying wave propagation. Field, regional and continental seismic data. What is the relevant geophysics? G. Papanicolaou, ACMAC-Crete Correlation based imaging 7/33

8 Simulated cross correlation of noise signals piece of signal recorded at x 1 t x 1 x 2 piece of signal recorded at x 2 t cross correlation x 1!x 2 t!500!400!300!200! When the spatial support of the noise sources (circles) extends over all space then the cross correlation function is symmetric. The positive and negative parts correspond to the Green s function between x 1 and x 2 and its anti-causal counterpart, respectively. For the scalar wave equation used here all the information about the Green s function is contained in the travel time. G. Papanicolaou, ACMAC-Crete Correlation based imaging 8/33

9 A model correlation-gf calculation The solution u(t, x) of the wave equation in an inhomogeneous medium: 1 2 u c 2 (x) t 2 xu = n ɛ (t, x). (2) The term n ɛ (t, x) models a random distribution of noise sources. It is a zero-mean stationary (in time) Gaussian process with autocorrelation function n ɛ (t 1, y 1 )n ɛ (t 2, y 2 ) = F ɛ (t 2 t 1 )θ(y 1 )δ(y 1 y 2 ). (3) We assume that the decoherence time of the noise sources is much smaller than typical travel times between sensors. If we denote with ɛ the (small) ratio of these two time scales, we can then write the time correlation function F ɛ in the form ( F ɛ t2 t ) 1 (t 2 t 1 ) = F, (4) ɛ where t 1 and t 2 are scaled relative to typical inter-sensor travel times. G. Papanicolaou, ACMAC-Crete Correlation based imaging 9/33

10 The model GF calculation continued The stationary solution of the wave equation has the integral representation u(t, x) = = t n ɛ (s, y)g(t s, x, y)dsdy n ɛ (t s, y)g(s, x, y)dsdy, (5) where G(t, x, y) is the time-dependent Green s function. It is the fundamental solution of the wave equation 1 2 G c 2 (x) t 2 xg = δ(t)δ(x y), (6) starting from G(0, x, y) = t G(0, x, y) = 0 (and continued on the negative time axis by G(t, x, y) = 0 t 0). G. Papanicolaou, ACMAC-Crete Correlation based imaging 10/33

11 Wave equation in a homogeneous medium and cross correlations When the medium is homogeneous with background velocity c 0 and the source distribution extends over all space, i.e. θ 1, the signal amplitude diverges because the contributions from noise sources far away from the sensors are not damped. So we consider the solution u of the damped wave equation with T a an absorption time: 1 ( 1 + ) 2u x u = n ɛ (t, x). (7) T a t c 2 0 In a three-dimensional open medium with dissipation and with noise sources over all space, θ 1, then x 1 x 2 [ ] c 0 Ta F ɛ G(τ, x 1, x 2 ) F ɛ G( τ, x 1, x 2 ), τ C(1) (τ, x 1, x 2 ) = c2 0 T a 4 e (8) where stands for the convolution in τ and G is the Green s function of the wave equation in a homogeneous medium without dissipation: G(t, x 1, x 2 ) = 1 ( 4π x 1 x 2 δ t x 1 x 2 ). c 0 G. Papanicolaou, ACMAC-Crete Correlation based imaging 11/33

12 The basic GF calculation continued If the decoherence time of the sources is much shorter than the travel time (i.e., ɛ 1), then F ɛ behaves like a Dirac delta function and we have τ C(1) (τ, x 1, x 2 ) e x 1 x 2 [ ] c 0 Ta G(τ, x 1, x 2 ) G( τ, x 1, x 2 ), up to a multiplicative constant. We see from this simple model calculation how it is theoretically possible to estimate the travel time τ(x 1, x 2 ) = x 1 x 2 /c 0 between x 1 and x 2 from the cross correlation, with an accuracy of the order of the decoherence time of the noise sources. With real geophysical data the connection with (elastic) wave propagation models is exceedingly complex and has not been carried out. G. Papanicolaou, ACMAC-Crete Correlation based imaging 12/33

13 What do we want from models? What can we expect from data? From models: How does the directivity of the wave energy affect cross correlation calculations? From models: How can we interpret the observed dispersiveness of the waves? From models: Do we really want to fit surface wave models to the data? From DATA: Can we identify other interesting travel times or peaks in the cross correlations? From DATA: What about absorption and multiple scattering? G. Papanicolaou, ACMAC-Crete Correlation based imaging 13/33

14 Using stationary phase to understand energy directivity y x 1 x 1 x 1 x 2 x 2 x 2 If the ray going through x 1 and x 2 (solid line) enters into the source region (left figure), then the travel time can be estimated from the cross correlation. If this is not the case, then the cross correlation does not have a peak at the travel time (middle figure). Right figure: The main contribution to the singular components of the cross correlation is from pairs of ray segments issuing from a source y going to x 1 and to x 2, respectively (solid and dashed lines, respectively), and starting in the same direction. G. Papanicolaou, ACMAC-Crete Correlation based imaging 14/33

15 Noise signal cross correlations piece of signal recorded at x 1 t x 1 x 2 piece of signal recorded at x 2 t cross correlation x 1!x 2 t!500!400!300!200! piece of signal recorded at x 1 t x 1 x 2 piece of signal recorded at x 2 t cross correlation x 1!x 2 t!500!400!300!200! When the distribution of noise sources is spatially localized then the cross correlation function is not symmetric G. Papanicolaou, ACMAC-Crete Correlation based imaging 15/33

16 Noise signal cross correlations, continued piece of signal recorded at x 1 x 1 t piece of signal recorded at x 2 x 2 t cross correlation x 1!x 2 t!500!400!300!200! When the distribution of noise sources is spatially localized then the coherent part of the cross correlation function can be difficult or even impossible to distinguish if the axis formed by the two sensors is perpendicular to the main direction of energy flux from the noise sources. G. Papanicolaou, ACMAC-Crete Correlation based imaging 16/33

17 Cross correlations of Long Beach array seismic data A 4 8 kilometer array with an average 100m spacing (From Chang and Biondi, SEP, Stanford 2013). Note the highway (405) crossing the array. G. Papanicolaou, ACMAC-Crete Correlation based imaging 17/33

18 Cross correlations of Long Beach array seismic data II Snapshots of a virtual source in the center of the array at 5s lag on left and +5s lag on right for frequencies Hz. Speed is 250m/s. (From Chang and Biondi 2013). Note the directivity of the wave energy. G. Papanicolaou, ACMAC-Crete Correlation based imaging 18/33

19 Cross correlations of Long Beach array seismic data III Snapshots of a virtual source in the center of the array at 5s lag on left and +5s lag on right for frequencies Hz. Speed is 500m/s. (From Chang and Biondi 2013). Note the relative isotropy of the wave energy. G. Papanicolaou, ACMAC-Crete Correlation based imaging 19/33

20 Imaging with passive auxiliary arrays schematic x s x r y Active sensor array for imaging a reflector at y through a scattering medium. Here x s, x r, and y denote source, receiver and reflector locations, respectively. x s x q x q y An auxiliary passive array for imaging through a scattering medium. x s is the source location (primary array), x q and x q are receiver locations (auxiliary array) below the scattering medium, and y is the reflector. G. Papanicolaou, ACMAC-Crete Correlation based imaging 20/33

21 How is such an imaging setup realized? Original setup in seismic exploration imaging 1 : Use passive arrays placed in drill holes to image below them. Typical wave propagation parameters: Speed of propagation: 3km/s, wavelength: 100m (30Hz peak frequency), bandwidth: broadband (more than 50%), synthetic array size: up to 10 10km in marine surveys, range (depth) 5 6km. Placement of auxiliary arrays at depths of 5 6km is being contemplated for enhanced imaging to well below 6km. Proposed radar application: Ground based radar arrays to image and track (small) satellites using passive auxiliary (synthetic) arrays on high-flying platforms. Potential elimination of atmospheric turbulence effects below the flying platform (at 15 20km hight). 1 A. Bakulin and R. Calvert, The virtual source method: Theory and case study, Geophysics, 71 (2006), pp. SI139-SI150 G. Papanicolaou, ACMAC-Crete Correlation based imaging 21/33

22 Travel-time migration imaging with the primary, active array Travel-time migration imaging using the array data p(t, x s, x r ), with x s the source location and x r the receiver location, recorded at the active array on top, is done with N r N s I( y S ) = p ( ) T( x s, y S ) + T( y S, x r ), x r ; x s, r=1 s=1 where the travel time, assuming a homogeneous medium, is given by T( x, y) = x y c 0, G. Papanicolaou, ACMAC-Crete Correlation based imaging 22/33

23 The correlation based imaging algorithm with the passive auxiliary array Step 1: Compute cross-correlations of the recorded signals p(t, x q ; x s ) on the auxiliary array receiver located at x q when the illumination is from the primary array at x s, through the inhomogeneous medium: C T (τ, x q, x q ) = 1 T T 0 N s p(t, x q ; x s )p(t+τ, x q ; x s )dt, q, q = 1,..., N q s=1 Step 2: Migrate the cross-correlations as if they were primary data (that is, actively generated) on the auxiliary array: I( y S ) = N q q,q =1 C T ( T( xq, y S ) + T( y S, x q ), x q, x q ). where, assuming a homogeneous medium below the auxiliary array, T( x, y) = x y c 0, G. Papanicolaou, ACMAC-Crete Correlation based imaging 23/33

24 Simulations (with geophysical parameters) 2 Imaging through random media with an auxiliary array: Left: The active, primary imaging array is at the top, the random medium below it, the auxiliary array is just below the random medium, and the object to be imaged is the black dot at the bottom. Middle: The image obtained with standard migration with active array at the top. Right: The image obtained when the auxiliary array is used. It is remarkable that the effects of the random medium have been taken out almost entirely. 2 Signal to Noise Ratio estimation in passive correlation-based imaging, J. Garnier, G. Papanicolaou, A. Semin and C. Tsogka, SIAM J. Imaging Sci G. Papanicolaou, ACMAC-Crete Correlation based imaging 24/33

25 Remarks on the simulations Use of a 2D finite element code (high resolution) that can handle about numerical domains wavelengths. This explains why the scatterer is about 180 wavelengths (1800m) deep. One dimensional layer codes can go much deeper but they are too special and mostly of theoretical interest today. 2D (3D) random effects are very important because they increase the lateral diversity of the illumination from the primary array. If this is properly quantified it can be exploited to reduce data acquisition costs by reducing the size of the primary array. Secondary array need not be horizontal and can be synthetic. G. Papanicolaou, ACMAC-Crete Correlation based imaging 25/33

26 Why does the auxiliary array imaging with correlations work so well? There are two important effects that we need to understand and exploit When imaging with ambient noise sources (spatially uncorrelated), the cross-correlations on the auxiliary passive array behave like the Green s function between the two points on the array. This is true only in a homogeneous or slowly varying medium. It is not true in a strongly scattering medium. Clearly the random medium helps make surface illumination more diverse when it gets to the auxiliary array. This is good for correlation based imaging. But it does not help if the medium is layered! So why does correlation based imaging work? It is time-reversal that makes the random medium effects go away. G. Papanicolaou, ACMAC-Crete Correlation based imaging 26/33

27 How does time reversal work? Basic fact: The cross correlation of two signals is the convolution of one with the time-reversed version of the other. Recall that C T (τ, x q, x q ) = 1 T T 0 N s p(t, x q ; x s )p(t+τ, x q ; x s )dt, q, q = 1,..., N q s=1 Interpretation: The cross correlation of signals on the auxiliary array is the same as having an impulsive point source source at one location x q, a time reversal mirror at the primary array on top, and an observation point x q of the time-reversed field at the other location on the auxiliary array. Main TR facts 3 : The TR field is peaked very near the source and is statistically stable (self-averaging) relative to the random medium 3 Theory and applications of time reversal and interferometric imaging, L. Borcea, G. Papanicolaou and C. Tsogka. Inverse Problems, vol 19, (2003), pp G. Papanicolaou, ACMAC-Crete Correlation based imaging 27/33

28 Why does it work, continued Conclusion: The cross-correlations computed on the auxiliary array are very different from the ones in a homogeneous medium. They do not peak at the travel time between the two points. They peak near zero lag, and only when the points are near each other. Effect: This is perfect for imaging below the auxiliary array, where the medium is (assumed) homogeneous, because there, which means at later lag times, the cross correlation will have peaks as we expect from the theory in homogeneous media. The effect of the random medium has thus been removed! This has been proven in detail for (i) randomly layered media and, (ii) for media for which the paraxial approximation can be used 4. The paraxial approximation is particularly useful for microwave propagation in a turbulent atmosphere. 4 Correlation based virtual source imaging in strongly scattering media, Josselin Garnier and George Papanicolaou. Inverse Problems, vol 28, (2012). G. Papanicolaou, ACMAC-Crete Correlation based imaging 28/33

29 Concluding remarks Passive sensor imaging with noise sources and in scattering media has both theoretical and applied potential. It is mathematically challenging and draws on many related research areas (wave propagation, random media, statistical signal processing, optimization,...). Passive auxiliary arrays can have a huge impact in imaging through inhomogeneities. It is the first time in imaging that the real benefits of time-reversal are being used. Theoretical issues, resolution theory for example, can often be reduced to already established results in active sensor imaging, both for distributed sensors and for arrays. G. Papanicolaou, ACMAC-Crete Correlation based imaging 29/33

30 Predicting volcanic eruptions The schematic of the seismic sensor deployment around the Piton de la Fournese volcano (by M. Campillo in Grenoble and collaborators, 2012). G. Papanicolaou, ACMAC-Crete Correlation based imaging 30/33

31 Echolocation in nature: Dolphins Ultrasonic echolocation by the dolphin is about ten times more accurate than what one expects from the basic physics of array imaging. (From: W. Au and J. Simmons, Echolocation in dolphins and bats, Physics Today, 2007) G. Papanicolaou, ACMAC-Crete Correlation based imaging 31/33

32 Echolocation in nature: Bats Ultrasonic echolocation by bats is much more accurate than what one expects from the basic physics of array imaging. (From: W. Au and J. Simmons, Echolocation in dolphins and bats, Physics Today, 2007) G. Papanicolaou, ACMAC-Crete Correlation based imaging 32/33

33 Basic resolution theory in array imaging Imaging formation key words: Coherent, array-data structure, active-passive, broadband-narrowband illumination-signals, homogeneous-inhomogeneous (random) medium. In small aperture regime, λ a L: Cross-range resolution: λl λl Range resolution: 2 a narrowband, and c 2 B broadband. For full (surrounding) arrays resolution is λ 2, narrowband and full recovery in broadband, in 3D. 2 λ (L/a) a. a y a λ L/a L L G. Papanicolaou, ACMAC-Crete Correlation based imaging 33/33

Imaging with Ambient Noise II

Imaging with Ambient Noise II Imaging with Ambient Noise II George Papanicolaou Stanford University Michigan State University, Department of Mathematics Richard E. Phillips Lecture Series April 21, 2009 With J. Garnier, University

More information

PASSIVE SENSOR IMAGING USING CROSS CORRELATIONS OF NOISY SIGNALS IN A SCATTERING MEDIUM

PASSIVE SENSOR IMAGING USING CROSS CORRELATIONS OF NOISY SIGNALS IN A SCATTERING MEDIUM PASSIVE SENSOR IMAGING USING CROSS CORRELATIONS OF NOISY SIGNALS IN A SCATTERING MEDIUM JOSSELIN GARNIER AND GEORGE PAPANICOLAOU Abstract. It is well known that the travel time or even the full Green s

More information

Imaging with Ambient Noise III

Imaging with Ambient Noise III Imaging with Ambient Noise III George Papanicolaou Stanford University http://math.stanford.edu/ papanico Department of Mathematics Richard E. Phillips Lectures April 23, 2009 With Liliana Borcea (Rice

More information

Ghost Imaging. Josselin Garnier (Université Paris Diderot)

Ghost Imaging. Josselin Garnier (Université Paris Diderot) Grenoble December, 014 Ghost Imaging Josselin Garnier Université Paris Diderot http:/www.josselin-garnier.org General topic: correlation-based imaging with noise sources. Particular application: Ghost

More information

(1) (2) (3) Main Menu. Summary. reciprocity of the correlational type (e.g., Wapenaar and Fokkema, 2006; Shuster, 2009):

(1) (2) (3) Main Menu. Summary. reciprocity of the correlational type (e.g., Wapenaar and Fokkema, 2006; Shuster, 2009): The far-field approximation in seismic interferometry Yingcai Zheng Yaofeng He, Modeling Imaging Laboratory, University of California, Santa Cruz, California, 95064. Summary Green s function retrieval

More information

RESOLUTION ENHANCEMENT FROM SCATTERING IN PASSIVE SENSOR IMAGING WITH CROSS CORRELATIONS. Josselin Garnier. George Papanicolaou

RESOLUTION ENHANCEMENT FROM SCATTERING IN PASSIVE SENSOR IMAGING WITH CROSS CORRELATIONS. Josselin Garnier. George Papanicolaou Inverse Problems and Imaging Volume 8, No. 3, 214, 645 683 doi:1.3934/ipi.214.8.645 RESOLUTION ENHANCEMENT FROM SCATTERING IN PASSIVE SENSOR IMAGING WITH CROSS CORRELATIONS Josselin Garnier Laboratoire

More information

Séminaire Laurent Schwartz

Séminaire Laurent Schwartz Séminaire Laurent Schwartz EDP et applications Année 211-212 Josselin Garnier Identification of Green s Functions Singularities by Cross Correlation of Ambient Noise Signals Séminaire Laurent Schwartz

More information

Distributed Sensor Imaging

Distributed Sensor Imaging Distributed Sensor Imaging George Papanicolaou Department of Mathematics - Stanford University URL: http://georgep.stanford.edu In collaboration with: Gregoire Derveaux, INRIA, France Chrysoula Tsogka,

More information

Séminaire Laurent Schwartz

Séminaire Laurent Schwartz Séminaire Laurent Schwartz EDP et applications Année 2013-2014 Josselin Garnier Multiscale analysis of wave propagation in random media. Application to correlation-based imaging Séminaire Laurent Schwartz

More information

CONSTRUCTION OF EMPIRICAL GREEN S FUNCTIONS FROM DIRECT WAVES, CODA WAVES, AND AMBIENT NOISE IN SE TIBET

CONSTRUCTION OF EMPIRICAL GREEN S FUNCTIONS FROM DIRECT WAVES, CODA WAVES, AND AMBIENT NOISE IN SE TIBET Proceedings of the Project Review, Geo-Mathematical Imaging Group (Purdue University, West Lafayette IN), Vol. 1 (2009) pp. 165-180. CONSTRUCTION OF EMPIRICAL GREEN S FUNCTIONS FROM DIRECT WAVES, CODA

More information

Effects of Surface Geology on Seismic Motion

Effects of Surface Geology on Seismic Motion 4 th IASPEI / IAEE International Symposium: Effects of Surface Geology on Seismic Motion August 23 26, 2011 University of California Santa Barbara TOMOGRAPHIC ESTIMATION OF SURFACE-WAVE GROUP VELOCITY

More information

Interferometric array imaging in clutter

Interferometric array imaging in clutter INSTITUTE OF PHYSICS PUBLISHING Inverse Problems 1 005 1419 1460 INVERSE PROBLEMS doi:10.1088/066-5611/1/4/015 Interferometric array imaging in clutter Liliana Borcea 1, George Papanicolaou and Chrysoula

More information

ACMAC s PrePrint Repository

ACMAC s PrePrint Repository ACMAC s rerint Repository Filtering deterministic layer effects in imaging Liliana Borcea and Fernando Gonzalez del Cueto and George apanicolaou and Chrysoula Tsogka Original Citation: Borcea, Liliana

More information

Coherent imaging without phases

Coherent imaging without phases Coherent imaging without phases Miguel Moscoso Joint work with Alexei Novikov Chrysoula Tsogka and George Papanicolaou Waves and Imaging in Random Media, September 2017 Outline 1 The phase retrieval problem

More information

Kinematics in Iterated Correlations of a Passive Acoustic Experiment

Kinematics in Iterated Correlations of a Passive Acoustic Experiment Kinematics in Iterated Correlations of a Passive Acoustic Experiment Sjoerd de Ridder and George Papanicolaou ABSTRACT Correlating ambient seismic noise can yield the inter-station Green s function, but

More information

Seismic Coda Waves. L. Margerin. CNRS, Toulouse, France

Seismic Coda Waves. L. Margerin. CNRS, Toulouse, France Mesoscopic Physics in Complex Media, 11 (21) DOI:1.151/iesc/21mpcm11 Owned by the authors, published by EDP Sciences, 21 Seismic Coda Waves L. Margerin CNRS, Toulouse, France In collaboration with B. Van

More information

Virtual Seismometers in the Subsurface of the Earth from Seismic Interferometry

Virtual Seismometers in the Subsurface of the Earth from Seismic Interferometry 1 Virtual Seismometers in the Subsurface of the Earth from Seismic Interferometry Andrew Curtis 1,2, Heather Nicolson 1,2,3, David Halliday 1,2, Jeannot Trampert 4, Brian Baptie 2,3 1 School of GeoSciences,

More information

Stationary phase approximation in the ambient noise method revisited

Stationary phase approximation in the ambient noise method revisited Earthq Sci (010)3: 45 431 45 Doi: 10.1007/s11589-010-0741-7 Stationary phase approximation in the ambient noise method revisited Zhongwen Zhan 1 and Sidao Ni, 1 Mengcheng National Geophysical Observatory,

More information

Approximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets. University of California, LLNL

Approximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets. University of California, LLNL 147th Annual Meeting: Acoustical Society of America New York City, May 25, 2004 Approximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets Speaker: James G. Berryman University

More information

Emergence rate of the time-domain Green s function from the ambient noise cross-correlation function

Emergence rate of the time-domain Green s function from the ambient noise cross-correlation function Emergence rate of the time-domain Green s function from the ambient noise cross-correlation function Karim G. Sabra, a Philippe Roux, and W. A. Kuperman Marine Physical Laboratory, Scripps Institution

More information

Signal to Noise Ratio estimation in passive correlation-based imaging

Signal to Noise Ratio estimation in passive correlation-based imaging Vol. xx, pp. x c xxxx Society for Industrial and Applied Mathematics x x Signal to Noise Ratio estimation in passive correlation-based imaging Josselin Garnier, George Papanicolaou, Adrien Semin, and Chrysoula

More information

Using SVD for improved interferometric Green s function retrieval

Using SVD for improved interferometric Green s function retrieval Using SVD for improved interferometric Green s function retrieval Gabriela Melo, Alison Malcolm, Dylan Mikesell 2,3, and Kasper van Wijk 3 Earth Resources Laboratory - Earth, Atmospheric, and Planetary

More information

Adaptive interferometric imaging in clutter and optimal illumination

Adaptive interferometric imaging in clutter and optimal illumination Adaptive interferometric imaging in clutter and optimal illumination Liliana Borcea and George Papanicolaou and Chrysoula Tsogka Computational and Applied Mathematics, MS 134, Rice University, 61 Main

More information

Filtering Deterministic Layer Effects in Imaging

Filtering Deterministic Layer Effects in Imaging SIAM REVIEW Vol. 54, No. 4, pp. 757 798 c 01 Society for Industrial and Applied Mathematics Filtering Deterministic Layer Effects in Imaging L. Borcea F. González del Cueto G. Papanicolaou C. Tsogka Abstract.

More information

Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a).

Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a). 7.1. Low-Coherence Interferometry (LCI) Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a). The light is split by the beam splitter (BS) and

More information

Earthquake ground motion prediction using the ambient seismic field

Earthquake ground motion prediction using the ambient seismic field Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L14304, doi:10.1029/2008gl034428, 2008 Earthquake ground motion prediction using the ambient seismic field Germán A. Prieto 1 and Gregory

More information

Using SVD for Improved Interferometric Green's Function Retrieval

Using SVD for Improved Interferometric Green's Function Retrieval Boise State University ScholarWorks CGISS Publications and Presentations Center for Geophysical Investigation of the Shallow Subsurface (CGISS) 9-1-2013 Using SVD for Improved Interferometric Green's Function

More information

A LOWER BOUND FOR PASSIVE SENSOR-NETWORK AUTO-LOCALIZATION

A LOWER BOUND FOR PASSIVE SENSOR-NETWORK AUTO-LOCALIZATION A LOWER BOUND FOR PASSIVE SENSOR-NETWORK AUTO-LOCALIZATION Rémy VINCENT, Mikael CARMONA CEA-Leti Department of Systems and Integration MINATEC Campus, Grenoble, FRANCE Olivier MICHEL, Jean-Louis LACOUME

More information

Coherent Interferometric Imaging in Clutter

Coherent Interferometric Imaging in Clutter Coherent Interferometric Imaging in Clutter Liliana Borcea, George Papanicolaou, Chrysoula Tsogka (December 9, 2005) ABSTRACT Coherent interferometry is an array imaging method in which we back propagate,

More information

Controlled source electromagnetic interferometry by multidimensional deconvolution: Spatial sampling aspects

Controlled source electromagnetic interferometry by multidimensional deconvolution: Spatial sampling aspects Controlled source electromagnetic interferometry by multidimensional deconvolution: Spatial sampling aspects Jürg Hunziker, Evert Slob & Kees Wapenaar TU Delft, The Netherlands ABSTRACT We review electromagnetic

More information

Summary. Introduction

Summary. Introduction : a philosophical view on seismic interferometry Kees Wapenaar*, Delft University of Technology, and Roel Snieder, Colorado School of Mines Summary We discuss the phenomenon of turning noise into signal

More information

Ambient Noise Tomography in the Western US using Data from the EarthScope/USArray Transportable Array

Ambient Noise Tomography in the Western US using Data from the EarthScope/USArray Transportable Array Ambient Noise Tomography in the Western US using Data from the EarthScope/USArray Transportable Array Michael H. Ritzwoller Center for Imaging the Earth s Interior Department of Physics University of Colorado

More information

Multifrequency Interferometric Imaging with Intensity-Only Measurements

Multifrequency Interferometric Imaging with Intensity-Only Measurements SIAM J. IMAGING SCIENCES Vol., No. 3, pp. 5 32 c 7 Society for Industrial and Applied Mathematics Multifrequency Interferometric Imaging with Intensity-Only Measurements Miguel Moscoso, Alexei Novikov,

More information

Th E Elastic Autofocusing SUMMARY

Th E Elastic Autofocusing SUMMARY Th E12 12 Elastic Autofocusing C.A. da Costa* (University of Edinburgh), M. avasi (University of Edinburgh), G. Meles (University of Edinburgh) & A. Curtis (University of Edinburgh) SUMMAY The method of

More information

Imagerie de la Terre profonde avec le bruit sismique. Michel Campillo (ISTERRE, Grenoble)

Imagerie de la Terre profonde avec le bruit sismique. Michel Campillo (ISTERRE, Grenoble) Imagerie de la Terre profonde avec le bruit sismique Michel Campillo (ISTERRE, Grenoble) Body waves in the ambient noise: microseisms (Gutenberg, Vinnik..) The origin of the noise in the period band 5-10s:

More information

Resolution Analysis of Passive Synthetic Aperture Imaging of Fast Moving Objects

Resolution Analysis of Passive Synthetic Aperture Imaging of Fast Moving Objects SIAM J. IMAGING SCIENCES Vol. 10, No. 2, pp. 665 710 c 2017 Society for Industrial and Applied Mathematics Resolution Analysis of Passive Synthetic Aperture Imaging of Fast Moving Objects L. Borcea, J.

More information

Local-scale cross-correlation of seismic noise from the Calico fault experiment

Local-scale cross-correlation of seismic noise from the Calico fault experiment DOI 10.1007/s11589-014-0074-z RESEARCH PAPER Local-scale cross-correlation of seismic noise from the Calico fault experiment Jian Zhang Peter Gerstoft Received: 27 November 2013 / Accepted: 28 February

More information

Robust seismic velocity change estimation using

Robust seismic velocity change estimation using submitted to Geophys. J. Int. Robust seismic velocity change estimation using ambient noise recordings E. Daskalakis,C.Evangelidis 2, J. Garnier 3, N. Melis 4,G.Papanicolaou 5 and C. Tsogka 6 Mathematics

More information

Seismic Noise Correlations. - RL Weaver, U Illinois, Physics

Seismic Noise Correlations. - RL Weaver, U Illinois, Physics Seismic Noise Correlations - RL Weaver, U Illinois, Physics Karinworkshop May 2011 Over the last several years, Seismology has focused growing attention on Ambient Seismic Noise and its Correlations. Citation

More information

27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies

27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies SHORT PERIOD SURFACE WAVE DISPERSION MEASUREMENTS FROM AMBIENT SEISMIC NOISE IN NORTH AFRICA, THE MIDDLE EAST, AND CENTRAL ASIA Michael H. Ritzwoller 1, Nikolai M. Shapiro 1, Michael E. Pasyanos 2, Gregory

More information

Edinburgh Research Explorer

Edinburgh Research Explorer Edinburgh Research Explorer Nonlinear Scattering-based Imaging in Elastic Media Citation for published version: Ravasi, M & Curtis, A 213, 'Nonlinear Scattering-based Imaging in Elastic Media' Paper presented

More information

Final Report for DOEI Project: Bottom Interaction in Long Range Acoustic Propagation

Final Report for DOEI Project: Bottom Interaction in Long Range Acoustic Propagation Final Report for DOEI Project: Bottom Interaction in Long Range Acoustic Propagation Ralph A. Stephen Woods Hole Oceanographic Institution 360 Woods Hole Road (MS#24) Woods Hole, MA 02543 phone: (508)

More information

RESOLUTION ANALYSIS OF PASSIVE SYNTHETIC APERTURE IMAGING OF FAST MOVING OBJECTS

RESOLUTION ANALYSIS OF PASSIVE SYNTHETIC APERTURE IMAGING OF FAST MOVING OBJECTS RESOLUTION ANALYSIS OF PASSIVE SYNTHETIC APERTURE IMAGING OF FAST MOVING OBJECTS L. BORCEA, J. GARNIER, G. PAPANICOLAOU, K. SOLNA, AND C. TSOGKA Abstract. In this paper we introduce and analyze a passive

More information

2D Imaging in Random Media with Active Sensor Array

2D Imaging in Random Media with Active Sensor Array 2D Imaging in Random Media with Active Sensor Array Paper Presentation: Imaging and Time Reversal in Random Media Xianyi Zeng rhodusz@stanford.edu icme, Stanford University Math 221 Presentation May 18

More information

Mandatory Assignment 2013 INF-GEO4310

Mandatory Assignment 2013 INF-GEO4310 Mandatory Assignment 2013 INF-GEO4310 Deadline for submission: 12-Nov-2013 e-mail the answers in one pdf file to vikashp@ifi.uio.no Part I: Multiple choice questions Multiple choice geometrical optics

More information

Correlating seismic wave velocity measurements with mining activities at Williams Mine

Correlating seismic wave velocity measurements with mining activities at Williams Mine Underground Mining Technology 2017 M Hudyma & Y Potvin (eds) 2017 Australian Centre for Geomechanics, Perth, ISBN 978-0-9924810-7-0 https://papers.acg.uwa.edu.au/p/1710_19_rebuli/ Correlating seismic wave

More information

c 2009 Society for Industrial and Applied Mathematics

c 2009 Society for Industrial and Applied Mathematics MULTISCALE MODEL. SIMUL. Vol. 7, No. 3, pp. 167 1301 c 009 Society for Industrial and Applied Mathematics FILTERING DETERMINISTIC LAYER EFFECTS IN IMAGING L. BORCEA, F. GONZÁLEZ DEL CUETO, G. AANICOLAOU,

More information

Investigation of long period amplifications in the Greater Bangkok basin by microtremor observations

Investigation of long period amplifications in the Greater Bangkok basin by microtremor observations Proceedings of the Tenth Pacific Conference on Earthquake Engineering Building an Earthquake-Resilient Pacific 6-8 November 2015, Sydney, Australia Investigation of long period amplifications in the Greater

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (2010) 183, 1648 1662 doi: 10.1111/j.1365-246X.2010.04825.x General representation theorem for perturbed media and application to Green s function retrieval

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (010) 18, 1509 1514 doi: 10.1111/j.1365-46X.010.04693. The relationship between noise correlation and the Green s function in the presence of degeneracy

More information

Imaging and time reversal in random media

Imaging and time reversal in random media 0 1 2 Imaging and time reversal in random media Liliana Borcea George Papanicolaou Chrysoula Tsogka James Berryman December 15, 2002 Abstract We present a general method for estimating the location of

More information

Title: Exact wavefield extrapolation for elastic reverse-time migration

Title: Exact wavefield extrapolation for elastic reverse-time migration Title: Exact wavefield extrapolation for elastic reverse-time migration ummary: A fundamental step of any wave equation migration algorithm is represented by the numerical projection of the recorded data

More information

An interferometric theory of source-receiver scattering and imaging

An interferometric theory of source-receiver scattering and imaging GEOPHYSICS, VOL. 75, NO. 6 NOVEMBER-DECEMBER 21; P. SA95 SA13, 7 FIGS. 1.119/1.3486453 An interferometric theory of source-receiver scattering and imaging David Halliday 1 and Andrew Curtis 2 ABSTRACT

More information

Deconvolution imaging condition for reverse-time migration

Deconvolution imaging condition for reverse-time migration Stanford Exploration Project, Report 112, November 11, 2002, pages 83 96 Deconvolution imaging condition for reverse-time migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT The reverse-time

More information

Imaging and monitoring with industrial seismic noise.

Imaging and monitoring with industrial seismic noise. Imaging and monitoring with industrial seismic noise. M. Campillo also : Boston, May 2016 Passive imaging: Long range correla@ons ()*&#"'!"#"$%"&'!!!!!!"! #! Source in A the signal recorded in B characterizes

More information

Cross-correlation of random fields: mathematical approach and applications

Cross-correlation of random fields: mathematical approach and applications Cross-correlation of random fields: mathematical approach and applications P. Gouédard 1,, L. Stehly 1,2, F. Brenguier 1,3, M. Campillo 1, Y. Colin de Verdière 4, E. Larose 1, L. Margerin 5, P. Roux 1,

More information

arxiv: v1 [physics.app-ph] 10 Aug 2018

arxiv: v1 [physics.app-ph] 10 Aug 2018 Retrieving the scattering strength using noise correlation in a reverberant medium arxiv:1808.03450v1 [physics.app-ph] 10 Aug 2018 Aida Hejazi Nooghabi, 1, 2 Lapo Boschi, 1 Philippe Roux, 3 and Julien

More information

SUMMARY ANGLE DECOMPOSITION INTRODUCTION. A conventional cross-correlation imaging condition for wave-equation migration is (Claerbout, 1985)

SUMMARY ANGLE DECOMPOSITION INTRODUCTION. A conventional cross-correlation imaging condition for wave-equation migration is (Claerbout, 1985) Comparison of angle decomposition methods for wave-equation migration Natalya Patrikeeva and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Angle domain common image gathers offer

More information

s=1 and the receivers at ( x r) Nr

s=1 and the receivers at ( x r) Nr ROLE OF SCATTERING IN VIRTUAL SOURCE ARRAY IMAGING JOSSELIN GARNIER AND GEORGE PAPANICOLAOU Abstract. We consider imaging in a scattering medium where the illumination goes through this medium but there

More information

Multicomponent imaging with distributed acoustic sensing Ivan Lim Chen Ning & Paul Sava, Center for Wave Phenomena, Colorado School of Mines

Multicomponent imaging with distributed acoustic sensing Ivan Lim Chen Ning & Paul Sava, Center for Wave Phenomena, Colorado School of Mines Multicomponent imaging with distributed acoustic sensing Ivan Lim Chen Ning & Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY The usage of Distributed Acoustic Sensing (DAS in geophysics

More information

MODELING AND INVERSION OF THE MICROTREMOR H/V SPECTRAL RATIO:

MODELING AND INVERSION OF THE MICROTREMOR H/V SPECTRAL RATIO: MODELING AND INVERSION OF THE MICROTREMOR H/V SPECTRAL RATIO: THE PHYSICAL BASIS BEHIND THE DIFFUSE FIELD APPROACH Francisco J. Sánchez-Sesma Instituto de Ingeniería, Universidad Nacional Autónoma de México

More information

Retrieving virtual reflection responses at drill-bit positions using seismic interferometry with drill-bit noise

Retrieving virtual reflection responses at drill-bit positions using seismic interferometry with drill-bit noise Geophysical Prospecting, 2015 doi: 10.1111/1365-2478.12292 Retrieving virtual reflection responses at drill-bit positions using seismic interferometry with drill-bit noise Yi Liu 1, Deyan Draganov 2, Kees

More information

Controlled source interferometry with noisy data Jürg Hunziker, Joost van der Neut, Evert Slob and Kees Wapenaar, Delft University of Technology

Controlled source interferometry with noisy data Jürg Hunziker, Joost van der Neut, Evert Slob and Kees Wapenaar, Delft University of Technology Jürg Hunziker, Joost van der Neut, Evert Slob and Kees Wapenaar, Delft University of Technology SUMMARY We investigate the effects of noise on controlled-source interferometry by multi-dimensional deconvolution

More information

Seismic interferometry versus spatial auto-correlation method on the regional coda of the NPE

Seismic interferometry versus spatial auto-correlation method on the regional coda of the NPE Seismic interferometry versus spatial auto-correlation method on the regional coda of the NPE Sjoerd de Ridder ABSTRACT A seismic recording of the non-proliferation experiment (NPE) contains the first

More information

Imaging and time reversal in random media

Imaging and time reversal in random media Imaging and time reversal in random media Liliana Borcea George Papanicolaou Chrysoula Tsogka James Berryman June 4, 2002 Abstract We present a general method for estimating the location of small, well-separated

More information

Seismic interferometry with antipodal station pairs

Seismic interferometry with antipodal station pairs GEOPHYSICAL RESEARCH LETTERS, VOL. 4, 1 5, doi:1.12/grl.597, 213 Seismic interferometry with antipodal station pairs Fan-Chi Lin 1 and Victor C. Tsai 1 Received 25 June 213; revised 19 August 213; accepted

More information

Seismic interferometry and estimation of the Green s function using Gaussian beams

Seismic interferometry and estimation of the Green s function using Gaussian beams Earthq Sci (1)3: 417 44 417 Doi: 1.17/s11589-1-74-8 Seismic interferometry and estimation of the Green s function using Gaussian beams Robert L. Nowack Department of Earth and Atmospheric Sciences, Purdue

More information

Modeling, migration, and inversion in the generalized source and receiver domain

Modeling, migration, and inversion in the generalized source and receiver domain Modeling, migration, and inversion in the generalized source and receiver domain Yaxun Tang ABSTRACT I extend the theory of Born modeling, migration and inversion to the generalized source and receiver

More information

Global P, PP, and PKP wave microseisms observed from distant storms

Global P, PP, and PKP wave microseisms observed from distant storms GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L23306, doi:10.1029/2008gl036111, 2008 Global P, PP, and PKP wave microseisms observed from distant storms Peter Gerstoft, 1 Peter M. Shearer, 1 Nick Harmon, 1 and

More information

Global propagation of body waves revealed by cross-correlation analysis of seismic hum

Global 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 information

Effects of Surface Geology on Seismic Motion

Effects of Surface Geology on Seismic Motion 4 th IASPEI / IAEE International Symposium: Effects of Surface Geology on Seismic Motion August 23 26, 2011 University of California Santa Barbara VELOCITY STRUCTURE INVERSIONS FROM HORIZONTAL TO VERTICAL

More information

Emergence of the Green s Functions from Noise and Passive Acoustic Remote Sensing of Ocean Dynamics

Emergence of the Green s Functions from Noise and Passive Acoustic Remote Sensing of Ocean Dynamics Emergence of the Green s Functions from Noise and Passive Acoustic Remote Sensing of Ocean Dynamics Oleg A. Godin CIRES/Univ. of Colorado and NOAA/OAR/Earth System Research Lab., R/PSD99, 325 Broadway,

More information

arxiv: v1 [math-ph] 15 Jul 2016

arxiv: v1 [math-ph] 15 Jul 2016 IMAGING IN RANDOM MEDIA WITH CONVEX OPTIMIZATION LILIANA BORCEA AND ILKER KOCYIGIT arxiv:67.465v math-ph 5 Jul 6 Abstract. We study an inverse problem for the wave equation where localized wave sources

More information

Th Guided Waves - Inversion and Attenuation

Th Guided Waves - Inversion and Attenuation Th-01-08 Guided Waves - Inversion and Attenuation D. Boiero* (WesternGeco), C. Strobbia (WesternGeco), L. Velasco (WesternGeco) & P. Vermeer (WesternGeco) SUMMARY Guided waves contain significant information

More information

arxiv:submit/ [physics.geo-ph] 18 Feb 2016

arxiv:submit/ [physics.geo-ph] 18 Feb 2016 ROBUST SEISMIC VELOCITY CHANGE ESTIMATION USING AMBIENT NOISE RECORDINGS E. DASKALAKIS, C. P. EVANGELIDIS, J. GARNIER, N. S. MELIS, G. PAPANICOLAOU AND C. TSOGKA arxiv:submit/48479 [physics.geo-ph] 8 Feb

More information

10. OPTICAL COHERENCE TOMOGRAPHY

10. OPTICAL COHERENCE TOMOGRAPHY 1. OPTICAL COHERENCE TOMOGRAPHY Optical coherence tomography (OCT) is a label-free (intrinsic contrast) technique that enables 3D imaging of tissues. The principle of its operation relies on low-coherence

More information

Finite difference elastic modeling of the topography and the weathering layer

Finite difference elastic modeling of the topography and the weathering layer Finite difference elastic modeling of the topography and the weathering layer Saul E. Guevara and Gary F. Margrave ABSTRACT Finite difference 2D elastic modeling is used to study characteristics of the

More information

LETTER Earth Planets Space, 64, 43 48, 2012

LETTER Earth Planets Space, 64, 43 48, 2012 LETTER Earth Planets Space, 64, 43 48, 1 Retrieval of long-wave tsunami Green s function from the cross-correlation of continuous ocean waves excited by far-field random noise sources on the basis of a

More information

Marchenko imaging. Kees Wapenaar 1, Jan Thorbecke 1, Joost van der Neut 1, Filippo Broggini 2,3, Evert Slob 1 and Roel Snieder 3

Marchenko imaging. Kees Wapenaar 1, Jan Thorbecke 1, Joost van der Neut 1, Filippo Broggini 2,3, Evert Slob 1 and Roel Snieder 3 CWP-796 Marchenko imaging Kees Wapenaar 1, Jan Thorbecke 1, Joost van der Neut 1, Filippo Broggini 2,3, Evert Slob 1 Roel Snieder 3 (1) Department of Geoscience Engineering, Delft University of Technology,

More information

San Jacinto Fault Zone and Sage Brush Flat High Frequency Experiments

San Jacinto Fault Zone and Sage Brush Flat High Frequency Experiments San Jacinto Fault Zone and Sage Brush Flat High Frequency Experiments Frank Vernon! Scripps Institution of Oceanography! University of California, San Diego!! 2015 Udine AUG! 13 March 2015! Southern California

More information

arxiv: v1 [physics.geo-ph] 22 Apr 2009

arxiv: v1 [physics.geo-ph] 22 Apr 2009 Stability of Monitoring Weak Changes in Multiply Scattering Media with Ambient Noise Correlation: Laboratory Experiments. Céline Hadziioannou, Eric Larose, Olivier Coutant, Philippe Roux and Michel Campillo

More information

Originally published as:

Originally published as: Originally published as: Ryberg, T. (2011): Body wave observations from cross correlations of ambient seismic noise: A case study from the Karoo, RSA. Geophysical Research Letters, 38, DOI: 10.1029/2011GL047665

More information

Analysis of ambient noise energy distribution and phase velocity bias in ambient noise tomography, with application to SE Tibet

Analysis of ambient noise energy distribution and phase velocity bias in ambient noise tomography, with application to SE Tibet Geophys. J. Int. (2009) doi: 10.1111/j.1365-246X.2009.04329.x Analysis of ambient noise energy distribution and phase velocity bias in ambient noise tomography, with application to SE Tibet Huajian Yao

More information

Scripps Institution of Oceanography Seminar University of California San Diego, October 7, 2005

Scripps Institution of Oceanography Seminar University of California San Diego, October 7, 2005 Scripps Institution of Oceanography Seminar University of California San Diego, October 7, 2005 TIME REVERSAL FOR RADAR IMAGING OF SMALL CONDUCTING AND/OR DIELECTRIC OBJECTS HIDDEN BY CLUTTER James G.

More information

Chapter 7. Seismic imaging. 7.1 Assumptions and vocabulary

Chapter 7. Seismic imaging. 7.1 Assumptions and vocabulary Chapter 7 Seismic imaging Much of the imaging procedure was already described in the previous chapters. An image, or a gradient update, is formed from the imaging condition by means of the incident and

More information

Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source

Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source Nevada Bureau of Mines and Geology Special Publication 36 Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source by Mendoza, C. 1 and Hartzell S. 2 1 Centro de Geociencias, Universidad

More information

Marchenko redatuming: advantages and limitations in complex media Summary Incomplete time reversal and one-sided focusing Introduction

Marchenko redatuming: advantages and limitations in complex media Summary Incomplete time reversal and one-sided focusing Introduction Marchenko redatuming: advantages and limitations in complex media Ivan Vasconcelos*, Schlumberger Gould Research, Dirk-Jan van Manen, ETH Zurich, Matteo Ravasi, University of Edinburgh, Kees Wapenaar and

More information

CBSE 12th Physics 2016 Unsolved Paper Delhi Board ARYAN INSTITUTE

CBSE 12th Physics 2016 Unsolved Paper Delhi Board ARYAN INSTITUTE CBSE 12th Physics 2016 Unsolved Paper Delhi Board CBSE 12th Physics 2016 Unsolved Paper Delhi Board TIME - 3HR. QUESTIONS - 26 THE MARKS ARE MENTIONED ON EACH QUESTION SECTION-A Q.1. A point charge +Q

More information

APPLICATION OF OPTIMAL BASIS FUNCTIONS IN FULL WAVEFORM INVERSION

APPLICATION OF OPTIMAL BASIS FUNCTIONS IN FULL WAVEFORM INVERSION METHODS AND APPLICATIONS OF ANALYSIS. c 2004 International Press Vol. 11, No. 3, pp. 345 352, September 2004 004 APPLICATION OF OPTIMAL BASIS FUNCTIONS IN FULL WAVEFORM INVERSION PING SHENG, GANG SUN,

More information

Correction of Ocean-Bottom Seismometer Instrumental Clock Errors Using Ambient Seismic Noise

Correction of Ocean-Bottom Seismometer Instrumental Clock Errors Using Ambient Seismic Noise Bulletin of the Seismological Society of America, Vol. 14, No. 3, pp. 1276 1288, June 214, doi: 1.1785/1213157 Correction of Ocean-Bottom Seismometer Instrumental Clock Errors Using Ambient Seismic Noise

More information

Phase velocities from seismic noise using beamforming and cross correlation in Costa Rica and Nicaragua

Phase velocities from seismic noise using beamforming and cross correlation in Costa Rica and Nicaragua Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L19303, doi:10.1029/2008gl035387, 2008 Phase velocities from seismic noise using beamforming and cross correlation in Costa Rica and Nicaragua

More information

Comparison of four techniques for estimating temporal change of seismic velocity with passive image interferometry

Comparison of four techniques for estimating temporal change of seismic velocity with passive image interferometry Earthq Sci (2010)23: 511 518 511 Doi: 10.1007/s11589-010-0749-z Comparison of four techniques for estimating temporal change of seismic velocity with passive image interferometry Zhikun Liu Jinli Huang

More information

Behavior and Sensitivity of Phase Arrival Times (PHASE)

Behavior and Sensitivity of Phase Arrival Times (PHASE) DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Behavior and Sensitivity of Phase Arrival Times (PHASE) Emmanuel Skarsoulis Foundation for Research and Technology Hellas

More information

Seismic Waves Propagation in Complex Media

Seismic Waves Propagation in Complex Media H4.SMR/1586-1 "7th Workshop on Three-Dimensional Modelling of Seismic Waves Generation and their Propagation" 5 October - 5 November 004 Seismic Waves Propagation in Comple Media Fabio ROMANELLI Dept.

More information

Mesoscopic Physics with Seismic waves «not nano but kilo»

Mesoscopic Physics with Seismic waves «not nano but kilo» Mesoscopic Physics with Seismic waves «not nano but kilo» Bart van Tiggelen Laboratoire de Physique et Modélisation des Milieux Condensés CNRS/University of Grenoble, France Collaborators: Michel Campillo

More information

ARTICLE IN PRESS. International Journal of Greenhouse Gas Control xxx (2012) xxx xxx. Contents lists available at SciVerse ScienceDirect

ARTICLE IN PRESS. International Journal of Greenhouse Gas Control xxx (2012) xxx xxx. Contents lists available at SciVerse ScienceDirect International Journal of Greenhouse Gas Control xxx (212) xxx xxx Contents lists available at SciVerse ScienceDirect International Journal of Greenhouse Gas Control j ourna l ho mepage: www.elsevier.com/locate/ijggc

More information

Source function estimation in extended full waveform inversion

Source function estimation in extended full waveform inversion Source function estimation in extended full waveform inversion Lei Fu The Rice Inversion Project (TRIP) May 6, 2014 Lei Fu (TRIP) Source estimation in EFWI May 6, 2014 1 Overview Objective Recover Earth

More information

Surface-wave tomography from ambient seismic noise of accelerograph networks in southern Korea

Surface-wave tomography from ambient seismic noise of accelerograph networks in southern Korea GEOPHYSICAL RESEARCH LETTERS, VOL. 33,, doi:10.1029/2006gl027044, 2006 Surface-wave tomography from ambient seismic noise of accelerograph networks in southern Korea Tae-Seob Kang 1 and Jin Soo Shin 1

More information

FloatSeis Technologies for Ultra-Deep Imaging Seismic Surveys

FloatSeis Technologies for Ultra-Deep Imaging Seismic Surveys FloatSeis Technologies for Ultra-Deep Imaging Seismic Surveys 25 th January, 2018 Aleksandr Nikitin a.nikitin@gwl-geo.com Geology Without Limits Overview 2011-2016 GWL Acquired over 43000 km 2D seismic

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (2014) 196, 1034 1042 Advance Access publication 2013 November 7 doi: 10.1093/gji/ggt434 Measuring of clock drift rates and static time offsets of ocean

More information