Discriminating non-seismic long-period pulses and noise to improve earthquake source inversion

Broadband seismometers produce artifacts resembling long-period pulses (non-seismic pulses) that degrade centroid moment tensor (CMT) estimations based on waveform inversion of broadband seismic records in long-period bands (50–200 s). We propose a method to discriminate non-seismic pulses and long-period noise from seismic signals, which can be applied to automatic CMT inversion analysis. In this method, we calculate source amplitudes as peak-to-peak displacement amplitudes in individual long-period seismic records after each event has been corrected for medium attenuation and geometric spreading and then estimate the ratios of individual source amplitudes to the minimum source amplitude. Because source amplitude ratios for non-seismic pulses tend to be greater than those of the seismic signals, we use seismic records in CMT estimations only if their source amplitude ratios are lower than a threshold value (R). We tested this method using broadband seismic data from the Philippines and found that reprocessed inversion solutions using this method showed a clear improvement when using R = 11, although focal mechanism estimations were not entirely stable. To investigate the general applicability of this method, we analyzed broadband seismic data from F-net in Japan. Our analysis indicated that source amplitude ratios in F-net data ranged up to about 20, indicating that the threshold value may be dependent on station density. Given that F-net is one of the highest density networks in the world, we may assume that a threshold value between 10 and 20 is appropriate for application of our method for most regional broadband networks. Our synthetic tests indicated that source amplitude ratios can be as high as 103, although observed ratios are only within the range 10–20. This suggests that we happened to observe only events having focal mechanisms with source amplitude ratios of 10–20. Alternatively, these high source amplitude ratios can be explained by distortion of radiation patterns in the long-period band, which reduces maximum source amplitude ratios and affects CMT estimates.

[1]  Hiroyuki Fujiwara,et al.  Recent Progress of Seismic Observation Networks in Japan , 2004 .

[2]  Sadaki Hori,et al.  Broadband Seismic Observation , 2009 .

[3]  E. R. Engdahl,et al.  Constraints on seismic velocities in the Earth from traveltimes , 1995 .

[4]  John Ristau,et al.  Implementation of Routine Regional Moment Tensor Analysis in New Zealand , 2008 .

[5]  S. Sipkin.  Rapid determination of global moment-tensor solutions , 1994 .

[6]  Tatsuhiko Saito,et al.  Characterization of scattered seismic wavefields simulated in heterogeneous media with topography , 2010 .

[7]  Nobuaki Kubo,et al.  Performance analysis of GPS augmentation using Japanese Quasi-Zenith Satellite System , 2004 .

[8]  T. Okada,et al.  Performance of regional distance centroid moment tensor inversion applied to the 2004 mid‐Niigata prefecture earthquake, Japan , 2006 .

[9]  E. Görgün,et al.  Seismicity of the 24 May 2014 Mw 7.0 Aegean Sea earthquake sequence along the North Aegean Trough , 2015 .

[10]  Sunarjo,et al.  Centroid moment tensor catalogue for Indonesia , 2010 .

[11]  J. Bonita,et al.  Regional Moment Tensor Analysis in the Philippines: CMT Solutions in 2012–2013 , 2015 .

[12]  Ishmael C. Narag,et al.  Performance of Broadband Seismic Network of the Philippines , 2015 .

[13]  Y. Kagan.  Simplified algorithms for calculating double-couple rotation , 2007 .

[14]  Sunarjo,et al.  Source model of an earthquake doublet that occurred in a pull-apart basin along the Sumatran fault, Indonesia , 2010 .

[15]  John H. Woodhouse,et al.  Determination of earthquake source parameters from waveform data for studies of global and regional seismicity , 1981 .

[16]  Cheng‐Horng Lin,et al.  Very low frequency earthquakes along the Ryukyu subduction zone , 2012 .

[17]  T. Furumura,et al.  Distortion of the apparent S-wave radiation pattern in the high-frequency wavefield: Tottori-Ken Seibu, Japan, earthquake of 2000 , 2009 .

[18]  M. Bouchon.  Discrete wave number representation of elastic wave fields in three-space dimensions , 1979 .

[19]  Ishmael C. Narag,et al.  Development and Operation of a Regional Moment Tensor Analysis System in the Philippines: Contributions to the Understanding of Recent Damaging Earthquakes , 2015 .

[20]  Mizuho Ishida,et al.  Automated Seismic Moment Tensor Determination by Using On-line Broadband Seismic Waveforms , 1998 .

[21]  J. Zahradník,et al.  Long-Period Pulses in Broadband Records of Near Earthquakes , 2005 .

[22]  Barbara Romanowicz,et al.  Regional and far-regional earthquake locations and source parameters using sparse broadband networks: A test on the Ridgecrest sequence , 1998, Bulletin of the Seismological Society of America.

[23]  K. Aki,et al.  Location of seismic events and eruptive fissures on the Piton de la Fournaise volcano using seismic amplitudes , 2003 .

[24]  C. Frohlich.  Triangle diagrams: ternary graphs to display similarity and diversity of earthquake focal mechanisms , 1992 .

[25]  Göran Ekström,et al.  The global CMT project 2004–2010: Centroid-moment tensors for 13,017 earthquakes , 2012 .

[26]  Donald V. Helmberger,et al.  The 23:19 aftershock of the 15 October 1979 Imperial Valley earthquake: More evidence for an asperity , 1985 .

[27]  Laura Scognamiglio,et al.  Real-Time Determination of Seismic Moment Tensor for the Italian Region , 2009 .

[28]  J. Ristau.  Update of Regional Moment Tensor Analysis for Earthquakes in New Zealand and Adjacent Offshore Regions , 2013 .

[29]  Jiří Zahradník,et al.  Toward Understanding Subtle Instrumentation Effects Associated with Weak Seismic Events in the Near Field , 2010 .

[30]  Takashi Furumura,et al.  Scattering of high-frequency seismic waves caused by irregular surface topography and small-scale velocity inhomogeneity , 2015 .

[31]  E. Görgün,et al.  The 2011 October 23 Mw 7.2 Van-Erciş, Turkey, earthquake and its aftershocks , 2013 .

[32]  M. Ruiz,et al.  Broadband seismic monitoring of active volcanoes using deterministic and stochastic approaches , 2009 .

[33]  Paul Bodin,et al.  Broadband Sensor Nonlinearity during Moderate Shaking , 2008 .

[34]  S. Takemura,et al.  Frequency and distance changes in the apparent P-wave radiation pattern: effects of seismic wave scattering in the crust inferred from dense seismic observations and numerical simulations , 2015 .

[35]  Taku Urabe,et al.  GRiD MT (grid-based real-time determination of moment tensors) monitoring the long-period seismic wavefield , 2009 .

[36]  Hitoshi Kawakatsu,et al.  Automated near‐realtime CMT inversion , 1995 .

[37]  S. Aoi,et al.  Strong localized asperity of the 2011 Tohoku-Oki earthquake, Japan , 2012, Earth, Planets and Space.

[38]  Y. Y.,et al.  3-D rotation of double-couple earthquake sources , 2006 .

[39]  M. Nakano,et al.  Waveform inversion in the frequency domain for the simultaneous determination of earthquake source mechanism and moment function , 2008 .

[40]  D. Kalafat,et al.  Focal mechanism determinations of earthquakes along the North Anatolian fault, beneath the Sea of Marmara and the Aegean Sea , 2015, Earth, Planets and Space.