Probabilistic inversion of Rayleigh-wave dispersion data: an application to Mt. Etna, Italy
Tóm tắt
We present a methodology for determining the elastic properties of the shallow crust from inversion of surface wave dispersion characteristics through a fully nonlinear procedure. Using volcanic tremor data recorded by a small-aperture seismic array on Mount Etna, we measured the surface waves dispersion curves with the multiple signal classification technique. The large number of measurements allows the determination of an a priori probability density function without the need of making any assumption about the uncertainties on the observations. Using this information, we successively conducted the inversion of phase velocities using a probabilistic approach. Using a wave-number integration method, we calculated the predicted dispersion function for thousands of 1-D models through a systematic grid search investigation of shear-wave velocities in individual layers. We joined this set of theoretical dispersion curves to the experimental probability density function (PDF), thus obtaining the desired structural model in terms of an a posteriori PDF of model parameters. This process allowed the representation of the objective function, showing the non-uniqueness of the solutions and providing a quantitative view of the uncertainties associated with the estimation of each parameter. We then compared the solution with the surface wave group velocities derived from diffuse noise Green’s functions calculated at pairs of widely spaced (~5–10 km) stations. In their gross features, results from the two different approaches are comparable, and are in turn consistent with the models presented in several earlier studies.
Tài liệu tham khảo
Aki J (1957) Space and time spectra of stationary stochastic waves with special reference to microtremors. Bull Earthq Res Inst Tokyo 25:415–457
Aki K, Richards PG (2002) Quantitative seismology, 2nd edn. University Science Book
Asten MW (1978) Geological control on the three-component spectra of Rayleigh wave microseism. Bull Seismol Soc Am 68:1623–1636
Bean C, Lokmer I, O’Brien G (2008) Influence of near-surface volcanic structure on long-period seismic signals and on moment tensor inversions: simulated examples from Mt. Etna. J Geophys Res 113:B08308
Bensen GD, Ritzwoller MH, Barmin MP, Levshin AL, Lin F, Moschetti MP, Shapiro NM, Yang Y (2007) Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements. Geoph J Int 169:1239–1260
Biondi BL, Kostov C (1989) High resolution velocity spectra using eigenstructure methods. Geophysics 54(7):832–842
Capon J (1969) High-resolution frequency-wavenumber spectrum analysis. Proc IEEE 57(8):1408–1418
Chiou S-J, Bolt BA (1993) Seismic wave slowness-vector estimation from broad-band array data. Geoph J Int 114(2):234–248
Corsaro RA, Miraglia L (2005) Dynamics of 2004-2005 Mt. Etna effusive eruption as inferred from petrologic monitoring. Geophys Res Lett 32:L13302
Cristiano L, Petrosino S, Saccorotti G, Ohrnberger M, Scarpa R (2010) Shear-wave velocity structure at Mt. Etna from inversion of Rayleigh-wave dispersion patterns (2s < T < 20s). Ann Geophys 53(2):69–78
Derode A, Larose E, Tanter M, De Rosny J, Tourim A, Campillo M, Fink M (2003) Recovering the Green’s function from field–field correlations in an open scattering medium. J Acoust Soc Am 113:2973–2976
Di Giulio G, Cornou C, Ohrnberger M, Wathelet M, Rovelli A (2006) Deriving wavefield characteristics and shear-velocity profiles from two-dimensional small-aperture arrays analysis of ambient vibrations in a small-size alluvial basin, Colfiorito, Italy. Bull Seismol Soc Am 96(5):1915–1933
Dziewonski AM, Bloch S, Landisman M (1969) A technique for the analysis of transient seismic signals. Bull Seismol Soc Am 59:427–444
Gardner GHF, Gardner LW, Gregory AR (1974) Formation velocity and density–the diagnostic basics for stratigraphics straps. Geophysics 39(6):770–780
Goldstein P, Archuleta RJ (1987) Array analysis of seismic signals. Geophys Res Lett 14(1):13–16
Haskell NA (1951) The dispersion of surface waves on multilayered media. Bull Seismol Soc Am 43(1):17–34
Herrmann RB (1987) Computer programs in seismology. Saint Louis University, Saint Louis, MO
Herrmann RB (2002) Computer programs seismology, an overview of synthetic seismogram computation. Saint Louis University, Missouri
Horike M (1985) Inversion of phase velocity of long-period microtremors to the S-wave velocity structure down to the basement in urbanized areas. J Phys Earth 33:59–96
Knopoff L (1964) A matrix method for elastic wave problems. Bull Seismol Soc Am 54(1):431–438
Larose E, Derode A, Corennec D, Margerin L, Campillo M (2005) Passive retrieval of Rayleigh waves in disoredered elastic media. Phys Rev E 72:046607
Lay e Wallace (1995) Modern global seismology. Academic Press
Molnar S, Dosso E, Cassidy JF (2010) Bayesian inversion of microtremor array dispersion data in southwestern British Columbia. Geophys J Int 183:923–940
Mosegaard K, Tarantola A (1995) Monte Carlo sampling of solutions to inverse problems. J Geoph Res 100(B7):12431–12447
Moser TJ, Van Eck T, Nolet G (1992) Hypocenter determination in strongly heterogeneous earth models using the shortest path method. J Geophys Res 97(B5):6563–6572
Murphy JR, Shah HK (1988) An analysis of the effects of site geology on the characteristics of near-field Rayleigh waves. Bull Seismol Soc Am 78:64–82
Nolet G (1981) Linearized inversion of (teleseismic) data. In: Cassinis R (ed) The solution of the inverse problem in geophysical interpretation. Plenum, pp 9–37
Park CB, Miller RD, Xia J (1999) Multichannel analysis of surface waves. Geophys 64(3):800–808
Petrosino S, La Rocca M, Del Pezzo E (1999) Shallow velocity model of the nothern flank of Stromboli Volcano, deduced by high velocity frequency surface wave dispersion. J Seismol 3:83–94
Ritzwoller MH, Shapiro NM, Levshin AL, Leahy GM (2001) The structure of the crust and upper mantle beneath Antartica and the surroundings oceans. J Geophys Res 106:30645–30670
Rost S, Thomas C (2002) Array seismology: methods and applications. Rev Geophys 40:3
Sabra KG, Gerstoft P, Roux P, Kuperman WA (2005) Extracting time-domain Green’s function estimated from ambient seismic noise. Geophys Res Lett 32:LO3310
Saccorotti G, Zuccarello L, Del Pezzo E, Ibanez J, Gresta S (2004) Quantitative analysis of the tremor wavefield at Etna Volcano, Italy. J Volcanol Geotherm Res 136:223–245
Sambridge M, Mosegaard K (2000) Monte Carlo methods in geophysical inverse problems. Rev Geophys 40:3
Satoh T, Kawase H, Matsushima SI (2001) Differences between site characteristics obtained from microtremors, S-waves, P-waves, and codas. Bull Seismol Soc Am 91:313–334
Scherbaum F, Hinzen KG, Ohrnberger M (2003) Determination of shallow shear wave velocity profiles in the Cologne/Germany area using ambient vibrations. Geophys J Int 152:597–612
Schmidt RO (1986) Multiple emitter location and signal parameter estimation. IEEE Trans Anten Propag 34(3):276–280
Shapiro NM, Campillo M (2004) Emergence of broadband Rayleigh waves from correlations of the ambient seismic noise. Geophys Res Lett 31:L07614
Shapiro NM, Ritzwoller MH (2002) Monte-Carlo inversion for a global shear-velocity model of the crust and upper mantle. Geophys J Int 151:88–105
Snieder R (2004) Extracting the Green’s function from the correlation of coda waves: a derivation based on stationary phase. Phys Rev E 69:046610
Tarantola A (1987) Inverse problem theory. Elsevier
Tarantola A (2005) Inverse problem theory and methods for model parameter estimation. Society for Industrial and Applied Mathematics, Philadelphia
Tokimatsu K (1995) Geotechnical site characteriszation using surface waves. In: Earthquake geotechnical engineering. Ishihara, Balkema, Rotterdam, pp 1333–1368
Wang H, Kaveh M (1985) Coherent signal–subspace processing for the detection and estimation of angles of arrival of multiple wide-band sources. IEEE Trans ASSP 33:823–831
Wapenaar K (2004) Retrieving the elastodynamic Green’s function of an arbitrary homogeneous medium by cross-correlation. Phys Rev Lett 93:254301
Wathelet M, Jongmans D, Ohrnberger M (2004) Surface waves inversion using a direct search algorithm and its application to ambient vibration measurements. Near Surf Geophys 2:211–221
Weaver RL, Lobkis OI (2001) On the emergence of the Green’s function in the correlations of a diffuse wavefield. J Acoust Soc Am 110:3011–3017
Wegler U, Seidl D (1997) Kinematic parameters of the tremor wavefield at Mt. Etna (Sicily). Geophys Res Lett 24(7):753–762
Xia J, Miller RD, Park CB (1999) Estimation of near-surface shear wave velocity by inversion of Rayleigh waves. Geophys 64(3):691–700