Citation: | Zhao S,Ma Z T,Sun Y C,Xue M,Wang X M. 2025. Detection of station clock error and polarity reversal based on surface wave phase measurement and P-wave arrival time. Acta Seismologica Sinica,47(3):356−373. DOI: 10.11939/jass.20230164 |
With the widespread deployment of global seismic networks, the amount of available seismic data has increased substantially, enhancing opportunities to investigate Earth’s dynamics. However, issues such as station clock errors and polarity reversals, which arise from timing system failures, deployment mistakes, or digitization errors, undermine the reliability of seismological research. Rapid and accurate detection and correction of these errors in terabyte-scale datasets are thus essential for effective data preprocessing.
Clock errors are primarily caused by asynchrony between internal and external station clocks. Existing detection methods include comparing P-wave arrival times at nearby stations, and using Green’s functions from ambient noise. The former requires high seismicity and clear phase records, and the latter is limited to stations within small-to-medium apertures due to signal attenuation. Neither method is suitable for global, long-term datasets spanning decades. Consequently, new approaches are needed for efficient and accurate detection of clock errors on a global scale. Polarity reversals, caused by incorrect instrument orientation, cable misconnection, or digitization symbol errors, also require reliable identification and correction to ensure data quality. Although the origins of clock errors and polarity reversals are different, both manifest as phase anomalies in waveform data, affecting measured travel times. Therefore, developing a unified method to detect both errors simultaneously is feasible and necessary.
The method in this paper utilizes the combination of seismic surface wave relative travel-time measurements and P-wave arrivals to efficiently detect both clock error and polarity reversal at stations. Firstly, globally observable long-period surface waves excited by events with MW≥5.5 are selected. The relative travel-times are obtained by measuring recordings of the same seismic event at different stations. The anomalies of the relative travel-time measurements at a station over a certain time period can reflect the clock errors or polarity reversals of the station. In order to further reduce the possibility of misjudgment in detection using surface wave relative travel-time measurements and to distinguish the clock error and polarity reversal, the proposed method combines the prediction of P-wave arrivals as a secondary check for the anomalous data.
The proposed method was applied to detect long-period surface waveform data from 2 010 stations of 30 networks during 2008−2012. The results show that twelve stations have clock error greater than 10 s in fourteen periods, and seven stations have polarity reversal in eight periods. Clock error detection technique of seismic wave data based on surface wave or P-wave techniques has their advantages and disadvantages. Surface waves have stronger energy and usually propagate over long distances, which can satisfy the detection of stations on a global scale. Due to the existence of interference factors during the measurement of surface waves, it is difficult to ensure absolute accuracy of the results. Although the detection accuracy of P-wave method is higher than that of the surface wave method, the lower signal-to-noise ratio of P-wave propagation over long distances could not satisfy the requirements of the method for a clear P-wave phase. The proposed method combines the two techniques to achieve better detection results.
The proposed method utilizes the cluster analysis for the efficient measurement of the relative travel-times of surface waves and it contains two checks. The first check by surface waves realizes the accurate and fast detection of a large number of data, and the secondary check by the P-wave arrival time ensures the accuracy of the final detection results. The application of the actual waveform data proves that the method can quickly and accurately detect the clock error and polarity reversal of the stations, which improves the quality of global surface wave velocity tomography and also provides a reference for the research on the inversion of source mechanism using relevant data. However, the proposed method based on seismic surface wave technology has its own limitations. Due to measurement errors and source parameter errors, the method can only detect large clock error. Additionally, the method is applicable to the clock error over long time and cannot effectively detect the anomalies of a single event or clock error of less than one day, which need to be further verified by combining with other data. Considering the limitations of the current method, further improvements may be worth considering in the future.
王俊,郑定昌,詹小艳,江昊琳,缪发军,高景春. 2013. 基于背景噪声互相关格林函数的单台时间误差估计[J]. 地震学报,35(6):888–901. doi: 10.3969/j.issn.0253-3782.2013.06.012
|
Wang J,Zheng D C,Zhan X Y,Jiang H L,Miao F J,Gao J C. 2013. Estimation of time error for a single station based on Green’s function by ambient noise cross-correlation[J]. Acta Seismologica Sinica,35(6):888–901 (in Chinese).
|
Abbas A,Zhu G H,Zi J P,Chen H,Yang H F. 2023. Evaluating and correcting short-term clock drift in data from temporary seismic deployments[J]. Earthquake Research Advance,3(2):100199. doi: 10.1016/j.eqrea.2022.100199
|
Anchieta M C,Wolfe C J,Pavlis G L,Vernon F L,Eakins J A,Solomon S C,Laske G,Collins J A. 2011. Seismicity around the Hawaiian islands recorded by the PLUME seismometer networks:Insight into faulting near Maui,Molokai,and Oahu[J]. Bull Seismol Soc Am,101(4):1742–1758. doi: 10.1785/0120100271
|
Bondár I,Storchak D. 2011. Improved location procedures at the international seismological centre[J]. Geophys J Int,186(3):1220–1244. doi: 10.1111/j.1365-246X.2011.05107.x
|
Butler R,Lay T,Creager K,Earl P,Fischer K,Gaherty J,Laske G,Leith B,Park J,Ritzwolle M,Tromp J,Wen L X. 2004. The global seismographic network surpasses its design goal[J]. Eos Trans Am Geophys Union,85(23):225-229.
|
Dziewonski A,Bloch S,Landisman M. 1969. A technique for the analysis of transient seismic signals[J]. Bull Seismol Soc Am,59(1):427–444. doi: 10.1785/BSSA0590010427
|
Ekström G,Nettles M,Abers G A. 2003. Glacial earthquakes[J]. Science,302(5645):622–624. doi: 10.1126/science.1088057
|
Ekström G,Dalton C A,Nettles M. 2006. Observations of time-dependent errors in long-period instrument gain at global seismic stations[J]. Seismol Res Lett,77(1):12–22. doi: 10.1785/gssrl.77.1.12
|
Ekström G,Nettles M. 2010. Performance of the GSN Station CASY-IU,1996−2009[R]. New York:Waveform Quality Center:1−15.
|
Ekström G. 2011. A global model of Love and Rayleigh surface wave dispersion and anisotropy,25−250 s[J]. Geophys J Int,187(3):1668–1686. doi: 10.1111/j.1365-246X.2011.05225.x
|
Gouédard P,Seher T,McGuire J J,Collins J A,van der Hilst R D. 2014. Correction of ocean-bottom seismometer instrumental clock errors using ambient seismic noise[J]. Bull Seismol Soc Am,104(3):1276–1288. doi: 10.1785/0120130157
|
Hartigan J A. 1975. Clustering Algorithms[M]. New York:John Wiley and Sons:351.
|
Houser C,Masters G,Shearer P,Laske G. 2008. Shear and compressional velocity models of the mantle from cluster analysis of long-period waveforms[J]. Geophys J Int,174(1):195–212. doi: 10.1111/j.1365-246X.2008.03763.x
|
Hurst A W,Bibby H M,Robinson R R. 2002. Earthquake focal mechanisms in the central Taupo volcanic zone and their relation to faulting and deformation[J]. N Z J Geol Geophys,45(4):527–536. doi: 10.1080/00288306.2002.9514989
|
Jin G,Gaherty J B. 2015. Surface wave phase-velocity tomography based on multichannel cross-correlation[J]. Geophys J Int,201(3):1383–1398. doi: 10.1093/gji/ggv079
|
Lentas K. 2021. Indications of seismic station phase reversals detected from parametric data in the ISC bulletin[J]. J Seismol,25(1):1–23. doi: 10.1007/s10950-020-09960-1
|
Ma Z,Masters G. 2014. A new global Rayleigh- and Love-wave group velocity dataset for constraining lithosphere properties[J]. Bull Seismol Soc Am,104(4):2007–2026. doi: 10.1785/0120130320
|
Ma Z T,Masters G. 2015. Effect of earthquake locations on Rayleigh wave azimuthal anisotropy models[J]. Geophys J Int,203(2):1319–1333. doi: 10.1093/gji/ggv369
|
Ma Z T,Masters G,Laske G,Pasyanos M. 2014. A comprehensive dispersion model of surface wave phase and group velocity for the globe[J]. Geophys J Int,199(1):113–135. doi: 10.1093/gji/ggu246
|
Niu F L,Li J. 2011. Component azimuths of the CEArray stations estimated from P-wave particle motion[J]. Earthquake Science,24(1):3–13. doi: 10.1007/s11589-011-0764-8
|
Petersen G M,Cesca S,Kriegerowski M,the AlpArray Working Group. 2019. Automated quality control for large seismic networks:Implementation and application to the AlpArray seismic network[J]. Seismol Res Lett,90(3):1177–1190. doi: 10.1785/0220180342
|
Pollitz F F,Stein R S,Sevilgen V,Bürgmann R. 2012. The 11 April 2012 east Indian Ocean earthquake triggered large aftershocks worldwide[J]. Nature,490(7419):250–253. doi: 10.1038/nature11504
|
Priestley K,McKenzie D,Ho T. 2018. A lithosphere-asthenosphere boundary:A global model derived from multimode surface‐wave tomography and petrology[G]//Lithospheric Discontinuities. Hoboken:John Wiley & Sons:111−123.
|
Roman D C,Neuberg J,Luckett R R. 2006. Assessing the likelihood of volcanic eruption through analysis of volcanotectonic earthquake fault-plane solutions[J]. Earth Planet Sci Lett,248(1/2):244–252.
|
Schneider J F,Aster R C,Powell L A,Meyer R P. 1988. Timing of portable seismographs from Omega navigation signals[J]. Bull Seismol Soc Am,78(2):1035–1035. doi: 10.1785/BSSA0780021035
|
Sens-Schönfelder C. 2008. Synchronizing seismic networks with ambient noise[J]. Geophys J Int,174(3):966–970. doi: 10.1111/j.1365-246X.2008.03842.x
|
Stehly L,Campillo M,Shapiro N M. 2007. Traveltime measurements from noise correlation:Stability and detection of instrumental time-shifts[J]. Geophys J Int,171(1):223–230. doi: 10.1111/j.1365-246X.2007.03492.x
|
Tape C,Liu Q Y,Maggi A,Tromp J. 2009. Adjoint tomography of the Southern California crust[J]. Science,325(5943):988–992. doi: 10.1126/science.1175298
|
Tian Y,Ning J Y,Yu C Q,Cai C,Tao K. 2013. Focal mechanism solutions of the 2008 Wenchuan earthquake sequence from P-wave polarities and SH/P amplitude ratios:New results and implications[J]. Earthquake Science,26(6):357–372. doi: 10.1007/s11589-014-0067-y
|
VanDecar J C,Crosson R S. 1990. Determination of teleseismic relative phase arrival times using multi-channel cross-correlation and least squares[J]. Bull Seismol Soc Am,80(1):150–169.
|
Xia Y J,Ni S D,Zeng X F,Xie J,Wang B S,Yuan S Y. 2015. Synchronizing intercontinental seismic networks using the 26 s persistent localized microseismic source[J]. Bull Seismol Soc Am,105(4):2101–2108. doi: 10.1785/0120140252
|
Yang W Z,Hauksson E,Shearer P M. 2012. Computing a large refined catalog of focal mechanisms for Southern California (1981−2010):Temporal stability of the style of faulting[J]. Bull Seismol Soc Am,102(3):1179–1194. doi: 10.1785/0120110311
|
Young B A,Chen K C,Huang B S,Chiu J M. 2017. Identification and elimination of data peculiarities in the strong‐motion downhole array in Taipei Basin[J]. Seismol Res Lett,88(1):82–95. doi: 10.1785/0220160164
|
Yukutake Y,Iio Y. 2017. Why do aftershocks occur? Relationship between mainshock rupture and aftershock sequence based on highly resolved hypocenter and focal mechanism distributions[J]. Earth Planets Space,69(1):68. doi: 10.1186/s40623-017-0650-2
|