A Quality Control Method for Teleseismic P‐Wave Receiver Functions

We introduce a new semiautomated quality control (QC) method for teleseismic P ‐wave receiver functions consisting of 13 editing procedures. These procedures are categorized into three groups based on (1) deconvolution attributes, (2) characteristics of each receiver function, and (3) the statistical attributes of a station gather. We evaluated the efficacy of individual procedure and combination of multiple procedures by applying them to receiver functions estimated from the EarthScope Ozark–Illinois–Indiana–Kentucky experiment for the Illinois basin area, and 10 years of USArray Transportable Array data for the contiguous United States. The method developed in this study significantly improves the efficiency and accuracy of receiver function QC. Our approach has three useful properties: (1) unlike a manual editing method, our approach produces reproducible results defined by a finite set of tunable parameters; (2) it provides an effective and efficient tool for processing large datasets containing millions of receiver functions; (3) the QC procedures are less dependent on phase similarities comparing to the automated scheme used by the EarthScope Automated Receiver Survey while producing comparable common conversion point stacking images.

[1]  Gary L. Pavlis,et al.  Array processing of teleseismic body waves with the USArray , 2010, Comput. Geosci..

[2]  G. Zandt,et al.  Pervasive lower-crustal seismic anisotropy in Southern California: Evidence for underplated schists and active tectonics , 2011 .

[3]  C. Langston.  An integrated study of crustal structure and regional wave propagation for southeastern Missouri , 1994 .

[4]  P. Molnar,et al.  Receiver functions in the western United States, with implications for upper mantle structure and dynamics , 2003 .

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

[6]  Michael H. Ritzwoller,et al.  Joint inversion of surface wave dispersion and receiver functions: a Bayesian Monte-Carlo approach , 2013 .

[7]  Hiroo Kanamori,et al.  Moho depth variation in southern California from teleseismic receiver functions , 2000 .

[8]  C. Langston,et al.  The Vertical Component P-Wave Receiver Function , 2001 .

[9]  Gary L. Pavlis,et al.  Three-dimensional, wavefield imaging of broadband seismic array data , 2011, Comput. Geosci..

[10]  Matthew J. Fouch,et al.  FuncLab: A MATLAB Interactive Toolbox for Handling Receiver Function Datasets , 2012 .

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

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

[13]  D. James,et al.  Receiver function imaging of upper mantle complexity beneath the Pacific Northwest, United States , 2010 .

[14]  Generalized iterative deconvolution for receiver function estimation , 2016 .

[15]  Thomas J. Owens,et al.  Automated Receiver Function Processing , 2005 .

[16]  M. Ritzwoller,et al.  Crustal layering in northeastern Tibet: a case study based on joint inversion of receiver functions and surface wave dispersion , 2015 .

[17]  Utkarsh Ayachit,et al.  The ParaView Guide: A Parallel Visualization Application , 2015 .

[18]  Stephen S. Gao,et al.  Determining crustal structure beneath seismic stations overlying a low‐velocity sedimentary layer using receiver functions , 2015 .

[19]  B. Kennett,et al.  Traveltimes for global earthquake location and phase identification , 1991 .

[20]  S. Sapkota,et al.  Imaging the Indian subcontinent beneath the Himalaya , 2005, Nature.

[21]  Anne F. Sheehan,et al.  Structure of the Bighorn Mountain region, Wyoming, from teleseismic receiver function analysis: Implications for the kinematics of Laramide shortening , 2014 .

[22]  H. P. Crotwell.  High data volume seismology: Surviving the avalanche , 2007 .

[23]  C. Haberland,et al.  Seismic imaging of subducting continental lower crust beneath the Pamir , 2013 .

[24]  A. Sheehan,et al.  Mantle discontinuity structure from midpoint stacks of converted P to S waves across the Yellowstone hotspot track , 1997 .

[25]  V. Schlindwein.  On the use of teleseismic receiver functions for studying the crustal structure of Iceland , 2006 .

[26]  H. Gilbert,et al.  Crustal signatures of the tectonic development of the North American midcontinent , 2016 .

[27]  G. Pavlis.  Three-dimensional wavefield imaging of data from the USArray: New constraints on the geometry of the Farallon slab , 2011 .

[28]  S. Kay,et al.  Structure of the crust and the lithosphere beneath the southern Puna plateau from teleseismic receiver functions , 2014 .

[29]  A. Velasco,et al.  Rocky Mountain evolution: Tying Continental Dynamics of the Rocky Mountains and Deep Probe seismic experiments with receiver functions , 2009 .

[30]  J. Lawrence,et al.  Characterizing the Main Himalayan Thrust in the Garhwal Himalaya, India with receiver function CCP stacking , 2013 .

[31]  Shaping wavelet effects in scattered wave imaging of P to S conversion data , 2015 .

[32]  Hersh Gilbert,et al.  Lithospheric and upper mantle structure of central Chile and Argentina , 2006 .

[33]  Wenjin Zhao,et al.  Receiver function imaging of crustal suture, steep subduction, and mantle wedge in the eastern India–Tibet continental collision zone , 2015 .

[34]  W. Yeck,et al.  Short Note Sequential H-κ Stacking to Obtain Accurate Crustal Thicknesses beneath Sedimentary Basins , 2013 .

[35]  K. Mahan,et al.  A method for mapping crustal deformation and anisotropy with receiver functions and first results from USArray , 2014 .

[36]  Stephen S. Gao,et al.  Mantle transition zone discontinuities beneath the contiguous United States , 2014 .