A reweighted damped singular spectrum analysis method for robust seismic noise suppression

(Multichannel)Singular spectrum analysis is considered as one of the most effective methods for seismic incoherent noise suppression.It utilizes the low-rank feature of seismic signal and regards the noise suppression as a low-rank reconstruction problem.However,in some cases the seismic geophones receive some erratic disturbances and the amplitudes are dramatically larger than other receivers.The presence of this kind of noise,called erratic noise,makes singular spectrum analysis(SSA)reconstruction unstable and has undesirable effects on the final results.We robustify the low-rank reconstruction of seismic data by a reweighted damped SSA(RD-SSA)method.It incorporates the damped SSA,an improved version of SSA,into a reweighted framework.The damping operator is used to weaken the artificial disturbance introduced by the low-rank projection of both erratic and random noise.The central idea of the RD-SSA method is to iteratively approximate the observed data with the quadratic norm for the first iteration and the Tukeys bisquare norm for the rest iterations.The RD-SSA method can suppress seismic incoherent noise and keep the reconstruction process robust to the erratic disturbance.The feasibility of RD-SSA is validated via both synthetic and field data examples.

Frequency dependent polarization analysis of ambient seismic noise recorded at a broadband seismometer in the central United States

We present a new approach to polarization analysis of seismic noise recorded by three-component seismometers. It is based on statistical analysis of frequency-dependent particle motion properties determined from a large number of time windows via eigenanalysis of the 3-by-3, Hermitian, spectral covariance matrix. We applied the algorithm to continuous data recorded in 2009 by the seismic station SLM, located in central North America. A rich variety of noise sources was observed. At low frequencies （〈0.05 Hz） we observed a tilt-related signal that showed some elliptical motion in the horizontal plane. In the microseism band of 0.05-0.25 Hz, we observed Rayleigh energy arriving from the northeast, but with three distinct peaks instead of the classic single and double frequency peaks. At intermediate frequencies of 0.5-2.0 Hz, the noise was dominated by non-fundamental-mode Rayleigh energy, most likely P and Lg waves. At the highest frequencies （〉3 Hz）, Rayleigh-type energy was again dominant in the form of Rg waves created by nearby cultural activities. Analysis of the time dependence of noise power shows that a frequency range of at least 0.02-1.0 Hz （much larger than the microseism band） is sensitive to annual, meteorologically induced sources of noise.

Comparison of ground truth location of earthquake from InSAR and from ambient seismic noise: A case study of the 1998 Zhangbei earthquake

Because ambient seismic noise provides estimated Green’s function （EGF） between two sites with high accuracy, Rayleigh wave propagation along the path connecting the two sites is well resolved. Therefore, earthquakes which are close to one seismic station can be well located with calibration extracting from EGF. We test two algorithms in locating the 1998 Zhangbei earthquake, one algorithm is waveform-based, and the other is traveltime-based. We first compute EGF between station ZHB （a station about 40 km away from the epicenter） and five IC/IRIS stations. With the waveform-based approach, we calculate 1D synthetic single-force Green’s functions between ZHB and other four stations, and obtain traveltime corrections by correlating synthetic Green’s functions with EGFs in period band of 10–30 s. Then we locate the earthquake by minimizing the differential travel times between observed earthquake waveform and the 1D synthetic earthquake waveforms computed with focal mechanism provided by Global CMT after traveltime correction from EGFs. This waveform-based approach yields a location which error is about 13 km away from the location observed with InSAR. With the traveltime-based approach, we begin with measuring group velocity from EGFs as well as group arrival time on observed earthquake waveforms, and then locate the earthquake by minimizing the difference between observed group arrival time and arrival time measured on EGFs. This traveltime-based approach yields accuracy of 3 km, Therefore it is feasible to achieve GT5 （ground truth location with accuracy 5 km） with ambient seismic noises. The less accuracy of the waveform-based approach was mainly caused by uncertainty of focal mechanism.

Seismic Noise Suppression for Ground-Based Investigation of an Inertial Sensor by Suspending the Electrode Cage

Performance test of a high precise accelerometer or an inertial sensor on the ground is inevitably limited by the seismic noise. A torsion pendulum has been used to investigate the performances of an electrostatic accelerometer, where the test mass is suspended by a fiber to compensate for its weight, and this scheme demonstrates an advantage, compared with the high-voltage levitation scheme, in which the effect of the seismic noise can be suppressed for a few orders of magnitude in low frequencies. In this work, the capacitive electrode cage is proposed to be suspended by another pendulum, and theoretical analysis shows that the effects of the seismic noise can be further suppressed for more than one order by suspending the electrode cage.

Local-scale cross-correlation of seismic noise from the Calico fault experiment

Most studies of seismic noise cross-correlation （NCC） have focused on regional/continental scale imaging using empirical surface-wave Green＇s functions extracted from primary （0.05-0.08 Hz） and secondary （0.1-0.16 Hz） microseisms. In this work, we present the NCC results at higher frequencies （〉0.5 Hz） from 6 months seismic noise recorded by a local array （ - 4 km aperture） deployed along the Calico fault in the Mojave Desert, California. Both fast and slow propagating waves are observed from the NCC record-sections. We compare the NCCs from sensor pairs that share a common sensor with the records of a borehole shot located very close to this common sensor. The result shows a good match of the slow surface-wave arrivals, indicating that the NCC method is able to recover unbiased surface-wave Green＇s functions at local scales. The strong body-wave NCC component is caused by the P waves generated offshore California. Along a SW-NE profile across the fault, we observe apparent P-wave arrivals and their reflections, which can be explained by a low-velocity- zone （LVZ） along the Calico fault. We calculat6 the LVZ width to be - 2.3 kin, and the P-wave velocity reduction within the LVZ to be -35 %. These estimates are consistent with other evidence for a relatively wide LVZ along the Calico fault.

Research on Real-time Monitoring of Abnormal Seismic Noise

The noise data in vertical component records of 85 seismic stations in Fujian Province during 2012 is used as the research object in this paper. The noise data is divided into fiveminute segments to calculate the power spectra. The high reference line and low reference line of station are then identified by drawing a probability density function graph( PDF)using the power spectral probability density function. Moreover, according to the anomalies of PDF graphs in 85 seismic stations,the abnormal noise is divided into four categories: dropped packet, low noise, high noise, and median noise anomalies.Afterwards,four selection methods are found by the high or low noise reference line of the stations,and the system of real-time monitoring of seismic noise is formed by combining the four selection methods. Noise records of 85 seismic stations in Fujian Province in July2013 are selected for verification,and the results show that the anomalous noise-recognition system could reach a 90% success rate at most stations and the effect of selection are very good. Therefore,it could be applied to the seismic noise real-time monitoring in stations.

Rayleigh Wave Group Velocity Maps of East Guangdong and Its Surrounding Regions from Ambient Seismic Noise Tomography

We perform Rayleigh wave tomography in east Guangdong and its surrounding regions by applying the ambient noise method to broadband data recorded at 26 stations from Guangdong,Fujian and Jiangxi Digital Seismic Networks.Cross-correlations of verticalcomponent ambient noise data are computed in one-day segments and stacked over seven months from March to September,2011.Then Rayleigh wave group dispersion curves are measured using the frequency-time analysis method.Group velocity maps at periods from5s to 15s are inverted.The resulting group velocity maps generally show good correlation with tectonic features,reflecting the velocity variations in the shallow crust.The basin areas are clearly resolved with lower group velocities at the short periods due to thick sedimentary layers,and the mountain areas with higher group velocities due to thin sedimentary layers.The variations of group velocity on the map can draw out the distribution of basins and mountains in study areas.The geothermal field can change the group velocity obviously,and lower group velocities are always found in high geothermal areas.The velocity maps indicate that a low-velocity layer may be found in the study areas.

Rayleigh waves from correlation of seismic noise in Great Island of Tierra del Fuego,Argentina:Constraints on upper crustal structure

In this study, the ambient seismic noise cross-correlation technique is applied to estimate the upper structure of the crust beneath Great Island of Tierra del Fuego（TdF）, Argentina, by the analysis of shortperiod Rayleigh wave group velocities. The island, situated in the southernmost South America, is a key area of investigation among the interaction between the South American and Scotia plates and is considered as a very seismically active one. Through cross-correlating the vertical components of ambient seismic noise registered at four broadband stations in TdF, we were able to extract Rayleigh waves which were used to estimate group velocities in the period band of 2.5-16 s using a timefrequency analysis. Although ambient noise sources are distributed in homogeneously, robust empirical Green＇s functions could be recovered from the cross-correlation of 12 months of ambient noise, The observed group velocities were inverted considering a non-linear iterative damped least-squares inversion procedure and several 1-D shear wave velocity models of the upper crust were obtained.According to the inversion results, the S-wave velocity ranges between 1.75 and 3,7 km/s in the first10 km of crust, depending on the pair of stations considered. These results are in agreement to the major known surface and sub-surface geological and tectonic features known in the area. This study represents the first ambient seismic noise analysis in TdF in order to constraint the upper crust beneath this region.It can also be considered as a successful feasibility study for future analyses with a denser station deployment for a more detailed imaging of structure.

Analysis of the Influence of Wind Turbine Noise on Seismic Recordings at Two Wind Parks in Germany

We recorded continuous seismic noise close and around two wind parks with different geological site conditions in order to investigate the effects of wind turbines on seismic signals in the frequency interval between 1 and 10 Hz. In wind park Fraureuth-Beiersdorf (5 turbines, 11 km south from Zwickau, Saxony), we used small seismic networks, with 3 to 5 stations for 1 - 2 weeks. In Heinde (2 turbines, close to Hildesheim, Lower Saxony), we recorded 1 week with one station around 1 km away from the wind turbines and some additional stations for several hours only. With the spectrogram analysis of the data, we clearly identify the diurnal variation on the spectral amplitude separately from the noise generated by the turbines. The turbine noise appears at certain frequency bands around 2.2, 2.7, 3.3, 4.5, 5.2 and 6.6 Hz. A linear relation between the spectral amplitudes of these frequency bands and the wind velocity or rotation velocity of the turbines is clearly identified. The seismic signals produced by the operation of the wind turbines are not peaks at single frequencies, but look more like frequency bands with increased noise amplitudes. They could be identified up to at least 10 km in the case of Fraureuth-Beiersdorf. These bands depend on numerous parameters, i.e. wind turbine height, weight and construction, number of turbines, geology, etc. In both wind parks we also recorded along profiles with increasing distances from the wind turbines. With the analysis of these data, we propose an amplitude attenuation model for the wind park Fraureuth-Beiersdorf to a distance of 9 km, and for Heinde to a distance of 4 km. The attenuation models for both wind parks are quite different, depending probably on the local geology and topography.

Crustal S-wave velocity structure of the Yellowstone region using a seismic ambient noise method

The Yellowstone volcano is one of the largest active volcanoes in the world, and its potential hazards demand detailed seismological and geodetic studies. Previous studies with travel time tomography and receiver functions have revealed a low-velocity layer in the crust beneath the Yellowstone volcano, suggesting the presence of a magma chamber at depth. We use ambient seismic noise from regional seismic stations to retrieve short-period surface waves and then study the shallow shear velocity structure of the Yellowstone region by surface wave dispersion analysis. We first obtained a crustal model of the area outside of the Yellowstone volcano and then constructed an absolute shear wave velocity structure in combination with receiver function results for the crust beneath the Yellowstone volcano. The velocity model shows a low-velocity layer with shear velocity at around 1.3 km/s, suggesting that a large-scale magma chamber exists at shallow levels within the crust of the Yellowstone volcanic region.

Processing Seismic Ambient Noise Data to Obtain Reliable Surface Wave Dispersion Measurements

Ambient noise tomography is a rapidly emerging field of seismological research. This paper presents the current status of ambient noise data processing and its development history over the past several years, with the intention to explain and justify this development through salient examples. The ambient noise data processing procedure can be divided into four principal phases： ① single station data preparation; ② cross- correlation and temporal stacking; ③ measurements of dispersion curves （ performed with frequency-time analysis for both group and phase speeds） ; ④ quality control, including SNR analysis and selection of the acceptable measurements. In addition, we provide a specific solution for a better use of the seismic station data to ambient noise study.

Using Ambient Noise to Study the Seismic Velocity Changes Caused by the Rise and Fall of the Water Level in the Zipingpu Reservoir Region

We study the feature of media changes beneath the Zipingpu reservoir and discuss the process of permeation with the water level rise and fall of the reservoir from January 2005 to January 2008 from ambient noise cross correlation by using continuous seismic data recorded by the stations of Zipingpu seismic network and YZP station. A moving-window cross-spectrum technique has been used to calculate the relative seismic velocity changes between station pairs. Results revealed an obvious relationship between relative seismic velocity, and the water level changes with a time delay that may be caused by permeation during three main impoundments and two large scale disemboguements. Impoundment generates a fast and large impact on the superficial layer, and the changes of seismic velocity is the result of increased pressure and permeation during the impoundment. At the first impoundment, the main effect factor is pressure. During the next two process of impoundment, permeation becomes the main effect factor, affecting the fault at a depth of about 8kin.

Crustal velocity structures beneath North China revealed by ambient noise tomography

We collected continuous noise waveform data from January 2007 to February 2008 recorded by 190 broadband and 10 very broadband stations of the North China Seismic Array. The study region is divided into grid with interval 0.25°×0.25°, and group velocity distribution maps between 4 s and 30 s are obtained using ambient noise tomography method. The lateral resolution is estimated to be 20-50 km for most of the study area. We construct a 3-D S wave velocity model by inverting the pure path dispersion curve at each grid using a genetic algorithm with smoothing constraint. The crustal structure observed in the model includes sedimentary basins such as North China basin, Yanqing-Huailai basin and Datong basin. A well-defined low velocity zone is observed in the Beijing-Tianjin-Tangshan region in 22-30 km depth range, which may be related to the upwelling of hot mantle material. The high velocity zone near Datong, Shuozhou and Qingshuihe within the depth range of 1-23 km reveals stable characteristics of Ordos block. The Taihangshan front fault extends to 12 km depth at least.

Stationary phase approximation in the ambient noise method revisited

The method of extracting Green＇s function between stations from cross correlation has proven to be effective theoretically and experimentally. It has been widely applied to surface wave tomography of the crust and upmost mantle. However, there are still controversies about why this method works. Snieder employed stationary phase approximation in evaluating contribution to cross correlation function from scatterers in the whole space, and concluded that it is the constructive interference of waves emitted by the scatterers near the receiver line that leads to the emergence of Green＇s function. His derivation demonstrates that cross correlation function is just the convolution of noise power spectrum and the Green＇s function. However, his derivation ignores influence from the two stationary points at infinities, therefore it may fail when attenuation is absent. In order to obtain accurate noise-correlation function due to scatters over the whole space, we compute the total contribution with numerical integration in polar coordinates. Our numerical computation of cross correlation function indicates that the incomplete stationary phase approximation introduces remarkable errors to the cross correlation function, in both amplitude and phase, when the frequency is low with reasonable quality factor Q. Our results argue that the dis- tance between stations has to be beyond several wavelengths in order to reduce the influence of this inaccuracy on the applications of ambient noise method, and only the station pairs whose distances are above several （〉5） wavelengths can be used.

Effect of uneven noise source and/or station distribution on estimating the azimuth anisotropy of surface waves

With the development of the dense array,the surface wave velocity and azimuthal anisotropy under the array can be directly obtained by beamforming the noise cross-correlation functions(NCFs). However, the retrieval of the Green’s function by cross-correlating the seismic noise requires that the noise source has a uniform distribution. For the case with uneven noise source, the azimuthal dependence on the sources in the expression for the spatial coherence function, which corresponds to the NCF in the time domain,has the same form as the azimuthal dependence of the surface wave velocity in weakly anisotropic media. Therefore, the uneven noise source will affect the surface wave anisotropy extraction. In this study, three passive seismic methods, i.e.,beamforming, SPAC(spatial autocorrelation), and NCF, are compared to demonstrate that an uneven source distribution and uneven station distribution have equivalent effects on the outcome from each method. A beamforming method is proposed to directly extract the velocity and azimuthal anisotropy of surface waves. The effect of uneven noise source and/or station distribution on estimating the azimuth anisotropy of surface waves was investigated using data from the ChinArray Phase Ⅱ. A method for correcting the apparent anisotropy in beamforming results caused by an uneven station distribution is suggested.

Ambient noise surface wave tomography of the Makran subduction zone,south-east Iran:Implications for crustal and uppermost mantle structures

Seismic ambient noise of surface wave tomography was applied to estimate Rayleigh wave empirical Green＇s functions （EGFs） and then to study crust and uppermost mantle structure beneath the Makran region in south-east 1mn. 12 months of continuous data from January 2009 through January 2010, recorded at broadband seismic stations, were analyzed. Group velocities of the fundamental mode Rayleigh wave dispersion curves were obtained from the empirical Green＇s functions. Multiple- filter analysis was used to plot group velocity variations at periods from 10 to 50 s. Using group velocity dispersion curves, 1-D Vs velocity models were calculated between several station pairs. The final results demonstrate signifi- cant agreement to known geological and tectonic features. Our tomography maps display low-velocity anomaly with SW-NE trend, comparable with volcanic arc settings of the Makran region which may be attributable to the geometry of Arabian Plate subducting beneath the overriding the Lut block. The northward subducting Arabian Plate is deter- mined by high-velocity anomaly along the Straits of Hor- muz. At short periods （〈20 s）, there is a sharp transition boundary between low- and high-velocity transition zone with the NW trending at the western edge of Makran which is attributable to the Minab fault system.

Comparison of four techniques for estimating temporal change of seismic velocity with passive image interferometry

Passive image interferometry （PII） is becoming a powerful tool for detecting the temporal variations in the Earth＇s structure, which applies coda wave interferometry to the waveforrns from the cross-correlation of seismic ambient noise. There are four techniques for estimating temporal change of seismic velocity with PII： moving-window cross-correlation technique （MWCCT）, moving-window cross-spectrum technique （MWCST）, stretching technique （ST） and moving-window stretching technique （MWST）. In this paper, we use the continuous seismic records from a typical station pair near the Wenchuan Ms8.0 earthquake fault zone and generate three sets of waveforms by stacking cross-correlation function of ambient noise with different numbers of days, and then apply four techniques to processing the three sets of waveforms and compare their results. Our results indicate that the techniques based on moving-window （MWCCT, MWCST and MWST） are superior in detecting the change of seismic velocity, and the MWCST can give a better estimate of velocity change than the other moving-window techniques due to measurement error. We also investigate the clock errors and their influences on measuring velocity change. We find that when the clock errors are not very large, they have limited impact on the estimate of the velocity change with the moving-window techniques.

Detecting remotely triggered temporal changes around the Parkfield section of the San Andreas fault

Detecting temporal changes in fault zone properties at seismogenic depth have been a long-sought goal in the seismological community for many decades. Recent studies based on waveform analysis of repeating earthquakes have found clear temporal changes in the shallow crust and around active fault zones associated with the occurrences of large nearby and teleseismic earthquakes. However, repeating earthquakes only occur in certain locations and their occurrence times cannot be controlled, which may result in inadequate sampling of the interested regions or time periods. Recent developments in passive imaging via auto- and cross-correlation of ambient seismic wavefields （e.g., seismic noise, earthquake coda waves） provide an ideal source for continuous monitoring of temporal changes around active fault zones. Here we conduct a systematic search of temporal changes along the Parkfield section of the San Andreas fault by cross-correlating relatively high-frequency （0.4-1.3 Hz） ambient noise signals recorded by 10 borehole stations in the High Resolution Seismic Network. After using stretch/compressed method to measure the delay time and the decorrelation-index between the daily noise cross-correlation functions （NCCFs）, we find clear temporal changes in the median seismic velocity and decorrelation-index associated with the 2004 M6.0 Parkfield earthquake. We also apply the same procedure to the seismic data around five regional/teleseismic events that have triggered non-volcanic tremor in the same region, but failed to find any clear temporal changes in the daily NCCFs. The fact that our current technique can detect temporal changes from the nearby but not regional and teleseismic events, suggests that temporal changes associated with distance sources are very subtle or localized so that they could not be detected within the resolution of the current technique （-0.2%）.

Site survey and assessment for the planned seismogeodynamic monitoring network in the Republic of Armenia

We present results of a detailed analysis of data obtained from seismic geodynamic field studies conducted at proposed sites for the development of advanced seismic monitoring stations in the Republic of Armenia.These studies aim to determine the background seismic and geodynamic noise level around such sites.Finally,based on the received data and international standards,nine survey points were classified into respective classes according to their noise level.We also calculated minimum significant earthquake magnitude detectable by the proposed seismic network in different regions of Armenia and mapped it based on recorded regional earthquakes.The resulting map indicates that the proposed seismic monitoring network will provide homogenous initial data for the various seismically active regions of the territory of the Republic of Armenia.

Monitoring Media Velocity Variations with Coda Wave Interferometry

Multiply scattered waves are sensitive to media changes owing to the effect of repeated sampling,superposition and amplification. Based on this characteristic,small-medium changes could be detected by using coda wave interferometry. In recent years,coda wave interferometry has been widely used in estimating velocity variation with high precision in areas such as seismology and non-destructive testing. This paper systematically presents the principle and research status of coda wave interferometry,and especial focus is placed on the research of media velocity variations by using repeating earthquakes,artificial sources,and ambient noise. Applications of coda wave interferometry can contribute to the more subtle understanding of dynamic evolution process in the medium.