A global optimizing approach for waveform inversion of receiver functions

A global optimizing approach is developed and implemented to retrieve one-dimensional crustal structure by waveform inversion of teleseismic receiver functions. The global optimization for the inversion is performed using a Differential Evolution (DE) algorithm. This modeling approach allows the user to perturb, within a preset range of reasonable bounds, multiple parameters such as V"p, V"p/V"s, thickness and anisotropy of each layer to fit the receiver function waveforms. Compared with linear modeling methods, the global optimal solution can be achieved with fewer model parameters (e.g., a small number of layers) and hence eliminate potential artifacts in the final model. Receiver function bins with small ray parameter intervals are used in the inversion, which can reduce distortion caused by modeling a single receiver function stacked from many recordings spread over a wide range of epicenter distance. The efficacy of this global optimizing approach is demonstrated with synthetic datasets and real receiver functions from the permanent seismic station BJT.

[1]  Clifford H. Thurber,et al.  Joint inversion of gravity and arrival time data from Parkfield: New constraints on structure and hypocenter locations near the SAFOD drill site , 2004 .

[2]  Brian J. Mitchell,et al.  Crustal thickness variations in the Aegean region and implications for the extension of continental crust , 2005 .

[3]  Bohuslav Růžek,et al.  Determination of the Earthquake Hypocenter: A Challenge for the Differential Evolution Algorithm , 2005 .

[4]  Mikhail K. Kaban,et al.  Receiver function tomography of the central Tien Shan , 2004 .

[5]  M. Sambridge.  Geophysical inversion with a neighbourhood algorithm—I. Searching a parameter space , 1999 .

[6]  Zhang Ling.  Inversion of crustal velocity model and earthquake relocation by using Differential evolution algorithm , 2006 .

[7]  Jonathan M. Lees,et al.  Crust and upper mantle of Kamchatka from teleseismic receiver functions , 2002 .

[8]  L. P. Vinnik,et al.  Detection of waves converted from P to SV in the mantle , 1977 .

[9]  S. Roecker,et al.  Receiver functions for the Tien Shan Analog Broadband Network: Contrasts in the evolution of structures across the Talasso‐Fergana Fault , 1993 .

[10]  Vadim Levin,et al.  Seismic anisotropy indicators in Western Tibet: Shear wave splitting and receiver function analysis , 2008 .

[11]  Vadim Levin,et al.  P-SH Conversions in Layered Media with Hexagonally Symmetric Anisotropy: A CookBook , 1998 .

[12]  Vadim Levin,et al.  Crustal anisotropy in the Ural Mountains Foredeep from teleseismic receiver functions , 1997 .

[13]  G. Zandt,et al.  Neighbourhood inversion of teleseismic Ps conversions for anisotropy and layer dip , 2003 .

[14]  Liang Zhao,et al.  A detailed receiver function image of the sedimentary structure in the Bohai Bay Basin , 2005 .

[15]  Kwang-Hee Kim,et al.  A preliminary study of crustal structure in Taiwan region using receiver function analysis , 2004 .

[16]  G. Randall,et al.  On the nonuniqueness of receiver function inversions , 1990 .

[17]  Charles A. Langston,et al.  Structure under Mount Rainier, Washington, inferred from teleseismic body waves , 1979 .

[18]  Liu Qi,et al.  The crustal thickness and Poisson's ratio beneath the Capital Circle Region , 2009 .

[19]  B. Růžek,et al.  Differential Evolution Algorithm in the Earthquake Hypocenter Location , 2001 .

[20]  R. Storn,et al.  Differential Evolution: A Practical Approach to Global Optimization (Natural Computing Series) , 2005 .

[21]  Charles J. Ammon,et al.  The isolation of receiver effects from teleseismic P waveforms , 1991, Bulletin of the Seismological Society of America.

[22]  Thomas J. Owens.  Crustal structure of the Adirondacks determined from broadband teleseismic waveform modeling , 1987 .

[23]  Rainer Storn,et al.  Differential Evolution – A Simple and Efficient Heuristic for global Optimization over Continuous Spaces , 1997, J. Glob. Optim..

[24]  Q Wu,et al.  THE CRUSTAL STRUCTURE OF QINGHAI-XIZANG PLATEAU INFERRED FROM BROADBAND TELESEISMIC WAVEFORM , 1998 .

[25]  Vadim Levin,et al.  P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation , 1997 .