Prerpints.org logo
Preprint
Article

Singular points of the tremor of the earth's surface

Altmetrics

Downloads

67

Views

31

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

16 August 2023

Posted:

17 August 2023

You are already at the latest version

Alerts
Abstract
A method for studying properties of the earth's surface tremor, measured by means of GPS, is proposed. Four characteristics of tremor are considered: the entropy of the distribution of wavelet coefficients, the Donoho-Johnston wavelet index, and two estimates of the spectral slope. The anomalous areas of tremor are determined by estimating the probability densities of extreme values of the studied properties. The criteria for abnormal tremor behavior are based on the proximity to, or the difference between, tremor properties and white noise. The greatest deviation from the properties of white noise is characterized by entropy minima and spectral slope and DJ index maxima. This behavior of the tremor is called "active". The "passive" tremor behavior is characterized by the maximum proximity to the properties of white noise. The principal components approach provides weighted averaged density maps of these 2 variants of extreme distributions of parameters in a sliding time window of 3 years are considered. Singular points are the points of maximum average densities. The method is applied to the analysis of daily time series from a GPS network in the West USA in 2009-2022. It turned out that singular points of tremor form well-defined clusters.
Keywords: 
Subject: Environmental and Earth Sciences  -   Geophysics and Geology

1. Introduction

The thin structure of the time series of displacements of the earth’s surface, measured by means of GPS, is the subject of research by a large number of specialists. The maximum likelihood approach for estimating parameters of GPS time series models was used in [1,2,3,4]. The structure of the noise and estimates of errors were considered in [5,6,7]. Parametric models of GPS time series for the analysis of phase correlations were used in [8]. Seasonal effects in GPS tine series connected with hydrological loading and tectonic movements of the crust were studied in [9,10,11,12,13,14,15,16] by different methods including multidimensional statistic approach. Influence of ionospheric and tropospheric delays on the sensitivity of the GPS solutions were considered in [17]. Hypotheses about origin of “abnormal” harmonics in GPS time series were regarded in [18].
The coherence of the earth’s surface tremor was analyzed in [19,20]. In [21,22], fractal analysis of high-frequency GPS time series was used to assess seismic hazard in Japan and California.
This article presents a further development of the methods for analyzing the earth’s surface tremor, proposed in [23], and analyzes the high-frequency noise structure in the time series of the GPS station network in the Western United States during the time interval 2009-2022. The US West is one of the most seismically active regions. A dense network of GPS stations makes it possible to study in detail the features of the earth’s surface oscillations and try to compare them with global seismic activity. Note that earlier in [24], a study was made of the properties of seismic noise in Southern California in comparison with the seismic regime.

2. Materials

GPS ground displacement data taken from the Nevada Geodetic Laboratory website at: http://geodesy.unr.edu/gps_timeseries/tenv3/IGS14/.
These data represent displacements of the earth’s surface in three orthogonal directions, which we will further denote as E, N and U (Up) for displacements in the west-east, north-south and vertical directions, with a time step of 1 day [25].
In what follows, we will consider a regular grid of 100×150 nodes, covering an area with latitude from 32.5°N to 51°N and longitude from 115°W to 128.5°W. Figure 1(a) shows the positions of GPS points for measuring displacements of the earth’s surface. Figure 1(b) shows the distribution density of GPS stations obtained using the Gaussian kernel estimate [26] at the nodes of this regular grid with a smoothing radius of 0.25° (about 28 km).
We will consider estimates of the properties of GPS time series in a moving time window 365 days long with a shift of 7 days. For each time window, operable stations are determined from the condition that the number of gaps in the considered time window does not exceed 5% of the time window length, that is, that this number is not more than 18 values. Missing values are filled in with artificial data, which are found from the value to the left and right of the missing time interval of the same length as the length of the gap [19]. Figure 2 shows a graph of the number of working stations within a sliding time window of 365 days for 14 years of observations, 2009-2022.

3. Description of the time series statistics used

For each station in each time window, 4 surface tremor statistics are computed for the 3 ground displacement components.
Spectral exponent. Let x ( t ) ,   t = 1 , , , N be a sample of the GPS time series within a window of sample lengths N . Denote by S x x ( ω j ) the estimate of the signal x ( t ) power spectrum ω j , - a sequence of discrete values of frequencies with a step Δ ω , Δ ω = 2 π / ( N f ⋅ Δ s ) , j = 1 , ... , N f , Δ s , - a time step. Here N f is the minimum length, which is greater than or equal to the sample length N : N f = 2 n ≥ N , the index j in formula (6) varies up to N f / 2 due to the symmetry of the power spectrum of real signals with respect to the Nyquist frequency ω N f / 2 + 1 : S x x ( ω N f / 2 + 1 − k ) = S x x ( ω N f / 2 + 1 + k ) ,   k = 1 , ... , N f / 2 − 1 . We introduce the spectral exponent α by the formula:
log ( S x x ( T j ) ) = α ⋅ log ( T j ) ) + d + η j ,   j = 1 , ... , N f / 2
where period T j = 2 π / ω j , the value α is determined by the least squares method: ∑ j = 1 N f / 2 η j 2 → min d , α , η j is a sequence of random variables with zero mean. In each time window, the power spectrum was computed after transition to increments using an autoregressive model (maximum entropy estimate) [27] with an autoregressive order of N / 10 , i.e., 36 for a window of length 365 samples. The transition to increments before calculating the power spectrum aims to free from the dominant effect of low frequencies in the daily time series of GPS.
Wavelet-based entropy and wavelet-based spectral slope. The minimum normalized wavelet information entropy for the sample x ( t ) ,   t = 1 , , , N is determined by the formula:
E n = − ∑ k = 1 N p k ⋅ log ( p k ) / log ( N ) ,       p k = c k 2 / ∑ j = 1 N c j 2
Here c k denoted the coefficients of the orthogonal wavelet decomposition. Daubechies bases were used with the number of vanishing moments from 1 to 10 [28]. We chose such a basis from this set for which the value (2) is minimal for a given time window and for the considered component of the displacement of the earth’s surface. By formula (2) 0 ≤ E n ≤ 1 .
For the orthogonal wavelet which is found from minimum of entropy, the wavelet power spectrum can be calculated as the average values of the squares of the wavelet coefficients at each level of detail of the decomposition:
S k = ∑ j = 1 M k | c j ( k ) | 2 / M k
In formula (3) c j ( k ) are the wavelet coefficients at the detail level with number k ,     k =   1 , ... k max , index j =   1 , ... M k , where M k is the number of wavelet coefficients at the detail level with number k , which corresponds to the frequency band [ 1 2 ( k + 1 ) Δ   s ,   1 2 k Δ   s ] with the central period T k = 2 Δ   s 2 − k + 2 − ( k + 1 ) [28]. The maximum detail level number k max is chosen such that it contains at least a given number L min of wavelet coefficients. Recall that the number of wavelet coefficients at the level of detail decreases by a factor of 2 when the level number increases by one. The value L min = 4 was used in the calculations, which gives the value k max = 6 for the window length N = 365 . Values S k = S ( T k ) , are similar to conventional power spectra. The regression model is used to calculate the wavelet-based spectral slope:
log ( S ( T k ) ) = β ⋅ log ( T k ) ) + c + ε k ,   k = 1 , ... , k max
where ε k is a sequence of random variables with zero mean. The parameter β in formula (4) is the wavelet-based spectral exponent, the value of which can be found by the least squares method: ∑ k = 1 k max ε k 2 → min β , c .
Donoho-Johnstone index (DJ-index). The threshold T D J is defined by the formula [28,29]:
T D J = σ 2 ⋅ ln N
The threshold (5) separates rather large (informative) in their absolute values wavelet coefficients from other coefficients which are considered to be noisy. Thus, we can consider the dimensionless signal characteristics γ , 0 ≤ γ ≤ 1 , as the ratio of the number of the largest wavelet coefficients, for which the inequality | c k | >   T D J is fulfilled, to the number N of all wavelet coefficients.
The value σ in the formula (5) is the noise standard deviation estimate under the assumption that the noise is most concentrated in the 1st detail level of orthogonal wavelet decomposition. The estimate of standard deviation σ should be robust with respect to outliers in the values of the coefficients at the first level. To provide this, a median estimate of the standard deviation for a normal random variable is used:
σ = m e d { | c k ( 1 ) |   ,   k = 1 , ... , N / 2 } / 0.6745
For each station in each time window, 4 surface tremor statistics are computed for the 3 ground displacement components where сk(1) are wavelet coefficients at the first level of detail. The estimate of the standard deviation σ from formula (6) determines the value (5) as a threshold for extracting noise wavelet coefficients. The quantity (5) is known in wavelet analysis as the Donoho–Johnstone threshold, and the expression for this quantity is based on the formula for the asymptotic probability of the maximum deviations of Gaussian white noise [28,29]. The DJ index γ values can be interpreted as a measure of the non-stationarity of seismic noise. For stationary Gaussian white noise, the index γ is zero. In [30,31,32] DJ-index was used as one of the main statistics for investigating properties of global seismic noise.
Summarizing, we will briefly describe in the list the characteristics of the values of the tremor statistics for stationary white Gaussian noise, which is considered as a standard of a completely random oscillation:
1. The “usual” spectral exponent α is the coefficient of regression between the logarithm of the power spectrum values and the logarithm of the period (“spectral slope”). For white noise it is equal to zero (white noise has a “flat spectrum”).
2. Wavelet-based spectral exponent β - the regression coefficient between the logarithm of the mean values of the squared wavelet coefficients for each level of detail and the logarithm of the period of the central frequency corresponding to the level of detail of the wavelet decomposition (“spectral slope”). For white noise is zero.
3. Donoho-Johnstone threshold (DJ-index) γ - the ratio of the number of “large” in absolute magnitude orthogonal wavelet coefficients to the total number of coefficients. For white noise is zero.
4. Minimum entropy E n of the distribution of squared orthogonal wavelet coefficients in the enumeration of wavelet basis functions in the class of Daubechies wavelets with the number of vanishing moments from 1 to 10. For white noise, the entropy is maximum.
Since we use 4 statistics for 3 components of the ground displacement vector, we get 12 statistics for each type of tremor. Let’s put them in one matrix:
[ ( α E , max β E , max γ E , max E n E , min α N , max β N , max γ N , max E n N , min α U , max β U , max γ U , max E n U , min ) ( α E , min β E , min γ E , min E n E , max α N , min β N , min γ N , min E n N , max α U , min β U , min γ U , min E n U , max ) ]
On the left side of the matrix (7) there are the designations of the characteristics of active tremor, on the right side of the matrix (7) there are the designations of the characteristics of passive tremor. Recall that we are considering a regular grid of nodes with a size of 100 nodes in longitude and 150 nodes in latitude. For each time window with a length of 365 days with a shift of 7 days, the values of all 12 sequences of the distribution of tremor properties at the nodes of the regular grid are determined.

4. Extreme value probability density maps

Let U there be any value from the matrix (7), the estimates of which were obtained in a moving time window. For each grid node ( i , j ) and for each time window labeled t at the right end of the window, we find the 10 closest seismic working stations, which gives 10 values of U . Let’s take their median value U i j ( t ) . The values U i j ( t ) form an “elementary” map corresponding to a time window of 365 days. Consider the values of the parameter as a function of two-dimensional vectors z i j = ( x i , y j ) of longitudes and latitudes of nodes ( i , j ) in an explicit form: U i j ( t ) ≡ U ( t ) ( z i j ) . For each “elementary map” with a discrete time index t , we find the coordinates of the nodes z   m n ( t ) = ( x m ( t ) , y n ( t ) ) at which the extreme value U is reached with respect to all other nodes of the regular grid.
As noted above, the choice of the minimum for entropies and the maximum for spectral slopes is due to the considerations that the regions of active tremor should be characterized by the values of statistics that reflect the maximum deviation from the properties of white noise. For passive tremor, the opposite is true: maxima for entropies and minima for spectral slopes. A cloud of 2D vectors z   m n ( t ) considered within a certain time interval t ∈ [ t 0 , t 1 ] composes some random set. Let us estimate two-dimensional probability density function for each node of the regular grid. For this, the Parzen–Rosenblatt estimate with the Gaussian kernel function [26] will be applied:
p ( z i j | t   0 , t   1 ) = 1 2 π h 2 ( t   1 − t   0 + 1 ) ∑ t = t   0 t   1 ∑ m n exp ( − | z i j − z   m n ( t ) | 2 2 h 2 )
Here, h is the averaging radius, integer indices t   0 , t   1 enumerate the maps of each time window. The ( t   1 − t   0 + 1 ) equals the number of time windows in the considered time interval. The width of the smoothing band h = 0.25 ∘ was used, which is approximately equal to 28 km.
In our analysis, time windows of 2 lengths are used: a “short” window of 365 days in length, which is used to obtain estimates of tremor parameters, and a “long” window of 1095 days (3 years) in length, in which a sequence of estimates of the distribution densities of extreme values of statistics is obtained tremor. Estimates of characteristics from matrix (7) are calculated using “short” time windows 365 days long. Next, the probability densities of extreme values are calculated from the values that fall within the current “long” time window of 1095 days, which are also shifted by 7 days. Each “long” window includes 105 “short” time windows.
As a result of such assessments, 24 maps of the distribution density of extreme values of tremor parameters were constructed, presented in Figure 3, on which 12 maps for active tremor are in the upper half, and 12 maps for passive tremor are in the lower half.
Purely visually, a correlation of probability density maps of the distribution of extreme values of tremor statistics is noticeable. The principal component method [33] is using to extract the most common spatial characteristics. In each “long” 3-year time window, a correlation matrix of the size 12×12 between the probability densities of the distribution of extreme values of tremor statistics is calculated. The eigenvector corresponding to the maximum eigenvalue of the correlation matrix is found and the squares of its components were used as weights to calculate the average weighted probability density in each “long” window. The coordinates of the point that realizes the maxima of average density were found. These operations performed makes it possible to obtain the trajectory of geographic coordinates of points that realize the maximum values of the average probability density. These points are called singular tremor points.
The results of weighting the probability density maps of the distribution of extreme properties of tremor by the method of principal components are shown in Figure 4. On the map 4(a), red circles indicate the positions of singular points of active tremor, which are combined into 5 clusters. The clusters are numbered in descending order of the latitude of their centers of mass; the cluster numbers are also shown. On the map 4(b), blue circles indicate the positions of singular points of passive tremor, which form 3 clusters. On Figure 4(c) and (d) there are plots of histograms of the distribution of the number of singular points in each cluster for each type of tremor.
Figure 5 shows graphs of changes in the squares of the components of sequences of correlation matrices for probability density maps of the distribution of extreme values, respectively, for active (12 upper graphs) and passive (12 lower graphs). As noted above, these values are used as weight functions to obtain weighted average probability density maps. The result of such a weighing operation using the principal component method is shown in Figure 4. The graphs in Figure 5 show the average values of the weights calculated for all “long” time windows. For each tremor variant, 4 maximum average values of weights are highlighted in red - they select those components from matrix (7) that make the greatest contribution to obtaining the weighted average sum of probability densities.
The trajectories of changes in geographic latitudes and longitudes of singular points over time are shown in Figure 6(a) and (b). Along the trajectories of change of coordinates there are marks of belonging to clusters of stations. Figure 6(c) shows how the distances between the singular points of active and passive tremor change as the time window of the length 3 years moves from left to right. From this graph it can be seen that the concentration of singular points of both types of tremor in those regions where their positions are very close to each other, namely, cluster #2 of active singular points and cluster #1 of passive singular points mainly occurs in different time windows. The exception is the time intervals with labels of the right ends of the “long” 3-year windows 2012-2012.5 and 2015.5-2016, when the distances from singular points of different types are less than 100 km, but do not reach zero values.
Let λ ˜ 1 be the normalized maximum eigenvalue of the correlation matrix:
λ ˜ 1 = λ 1 / ∑ k = 1 q λ k
In formula (9) λ 1 , ... , λ q are the eigenvalues of the correlation matrix, sorted in descending order: λ k ≥ λ k + 1 ,   ∑ k = 1 q λ k = q . According to the principal component method [33], the value λ ˜ 1 is equal to the share of the total variance attributable to the first principal component and can be interpreted as a measure of synchronization of changes in the scalar components of the multivariate time series of the probability density values of the tremor statistics from the matrix in formula (7) in nodes regular grid covering the region under study.
Figure 7 shows graphs of changes in the eigenvalues of the correlation matrix for active and passive tremor. It can be seen from these plots that the active tremor (Figure 7(a)) is much more synchronized than the passive tremor (Figure 7(b)). Moreover, for active tremor, there is a very strong synchronization for the positions of the right end of the 3-year time window from 2013 to 2015. Given the length of the window, this synchronization takes a period of time from 2010 to 2015.
As was shown in [23,24], after a couple of the strongest earthquakes in the world, Maule in Chile on 2010.02.27 and Tohoku in Japan on 2011.11.03, there was an explosive increase in the radius of spatial correlations of global seismic noise properties. It is possible that the increase in synchronization of active tremor in the Western United States in 2010-2015 is also a response to this pair of mega earthquakes. Note that a sharp jump in the magnitude of global correlations of the earth’s surface tremor in 2010-2011 was shown in [19], and [34] showed that before the 2011 Tohoku earthquake in Japan, long-term correlations arose between seismic noise in California and Japan.

Conclusions

A method has been proposed for studying the properties of the earth’s surface tremor, measured by means of GPS, based on the use of oscillation statistics that describe the deviation or similarity to the properties of Gaussian white noise, considered as a standard of a stationary random signal. Based on the application of the principal component method to the estimation of the probability densities of the distribution of extreme properties of statistics in a moving time window, concepts of active and passive tremor and their singular points are introduced, and the trajectories of singular points are determined. The concept of tremor synchronization measure is introduced as the maximum eigenvalue of the correlation density matrix of distribution of extreme values of tremor statistics. The method is applied to the data of daily GPS time series on a network of stations located in the Western United States during the time interval 2009-2022. It was shown that singular points of both types of tremor form well-defined clusters. Estimates of measures of synchronization of active and passive are given. A hypothesis has been put forward that strong synchronization of active tremor in 2010-2015 caused by a jump in the global synchronization of the tremor of the earth’s surface and an explosive increase in the correlations of global seismic noise in 2010-2011.

Funding

This research received no external funding

Data Availability Statement

The open access data from the source: http://geodesy.unr.edu/gps_timeseries/tenv3/IGS14/ were used.

Acknowledgments

The work was carried out within the framework of the state assignment of the Institute of Physics of the Earth of the Russian Academy of Sciences (topic FMWU-2022-0018).

References

  1. Roncagliolo, P.A. , Garcнa J.G., Mercader P.I., Fuhrmann D.R., Muravchik C.H. Maximum-likelihood attitude estimation using GPS signals. Digital Signal Processing 2007, 17, 1089–1100. [Google Scholar] [CrossRef]
  2. Wang, F.; Li, H.; Lu, M. GNSS Spoofing Detection and Mitigation Based on Maximum Likelihood Estimation. Sensors 2017, 17, 1532. [Google Scholar] [CrossRef] [PubMed]
  3. Parkinson, B.W. Global Positioning System: Theory and Applications; AIAA: Reston, VA, USA, 1996. [Google Scholar]
  4. 4. Jing Xu, Jingsha He, Yuqiang Zhang, Fei Xu, and Fangbo Cai. A Distance-Based Maximum Likelihood Estimation Method for Sensor Localization in Wireless Sensor Networks. International Journal of Distributed Sensor Networks, 2080. [CrossRef]
  5. Agnew, D. The time domain behavior of power law noises, Geophys. Res. Lett. 1992, 19, 333–336. [Google Scholar] [CrossRef]
  6. Amiri-Simkooei, A. R., C. C. J. M. Tiberius, and P. J. G. Teunissen. Assessment of noise in GPS coordinate time series: Methodology and results, J. Geophys. Res. 2007, 112, B07413. [Google Scholar] [CrossRef]
  7. Bos, M. S., R. M. S. Fernandes, S. D. P. Williams, and L. Bastos, Fast error analysis of continuous GPS observations, J. Geod. 2008, 82, 157–166. [Google Scholar] [CrossRef]
  8. Kermarrec, G. , Schon, S. On modelling GPS phase correlations: A parametric model. Acta Geod Geophys. 2018, 53, 139–156. [Google Scholar] [CrossRef]
  9. Liu, B.; Xing, X.; Tan, J.; Xia, Q. Modeling Seasonal Variations in Vertical GPS Coordinate Time Series Using Independent Component Analysis and Varying Coefficient Regression. Sensors 2020, 20, 5627. [Google Scholar] [CrossRef]
  10. Tesmer, V.; Steigenberger, P.; van Dam, T.; Mayer-Gurr, T. Vertical deformations from homogeneously processed GRACE and global GPS long-term series. J. Geod. 2011, 85, 291–310. [Google Scholar] [CrossRef]
  11. Yan, J.; Dong, D.; Burgmann, R.; Materna, K.; Tan, W.; Peng, Y.; Chen, J. Separation of Sources of Seasonal Uplift in China Using Independent Component Analysis of GNSS Time Series. J. Geophys. Res. Solid Earth 2019, 124, 11951–11971. [Google Scholar] [CrossRef]
  12. Pan, Y.; Shen, W.-B.; Ding, H.; Hwang, C.; Li, J.; Zhang, T. The Quasi-Biennial Vertical Oscillations at Global GPS Stations: Identification by Ensemble Empirical Mode Decomposition. Sensors 2015, 15, 26096–26114. [Google Scholar] [CrossRef]
  13. Liu, B.; Dai, W.; Liu, N. Extracting seasonal deformations of the Nepal Himalaya region from vertical GPS position time series using Independent Component Analysis. Adv. Space Res. 2017, 60, 2910–2917. [Google Scholar] [CrossRef]
  14. Santamarнa-Gуmez, A.; Memin, A. Geodetic secular velocity errors due to interannual surface loading deformation. Geophys. J. Int. 2015, 202, 763–767. [Google Scholar] [CrossRef]
  15. Fu, Y.; Argus, D.F.; Freymueller, J.T.; Heflin, M.B. Horizontal motion in elastic response to seasonal loading of rain water in the Amazon Basin and monsoon water in Southeast Asia observed by GPS and inferred from GRACE. Geophys. Res. Lett. 2013, 40, 6048–6053. [Google Scholar] [CrossRef]
  16. Chanard, K.; Avouac, J.P.; Ramillien, G.; Genrich, J. Modeling deformation induced by seasonal variations of continental water in the Himalaya region: Sensitivity to Earth elastic structure. J. Geophys. Res. Solid Earth 2014, 119, 5097–5113. [Google Scholar] [CrossRef]
  17. He, X.; Montillet, J.P.; Fernandes, R.; Bos, M.; Yu, K.; Hua, X.; Jiang, W. Review of current GPS methodologies for producing accurate time series and their error sources. J. Geodyn. 2017, 106, 12–29. [Google Scholar] [CrossRef]
  18. Ray, J.; Altamimi, Z.; Collilieux, X.; van Dam, T. Anomalous harmonics in the spectra of GPS position estimates. GPS Solutions 2008, 12, 55–64. [Google Scholar] [CrossRef]
  19. Lyubushin, A. Global coherence of GPS-measured high-frequency surface tremor motions. GPS Solut. 2018, 22, 116. [Google Scholar] [CrossRef]
  20. Lyubushin, A. Field of coherence of GPS-measured earth tremors. GPS Solut. 2019, 23, 120. [Google Scholar] [CrossRef]
  21. Filatov, D.M.; Lyubushin, A.A. Fractal analysis of GPS time series for early detection of disastrous seismic events. Phys. A 2017, 469, 718–730. [Google Scholar] [CrossRef]
  22. Filatov, D.M.; Lyubushin, A.A. Precursory Analysis of GPS Time Series for Seismic Hazard Assessment. Pure Appl. Geophys. 2019, 177, 509–530. [Google Scholar] [CrossRef]
  23. Lyubushin, A. Identification of Areas of Anomalous Tremor of the Earth’s Surface on the Japanese Islands According to GPS Data. Appl. Sci. 2022, 12, 7297. [Google Scholar] [CrossRef]
  24. Lyubushin, A.A. Seismic Noise Wavelet-Based Entropy in Southern California. Journal of Seismology 2021, 25, 25–39. [Google Scholar] [CrossRef]
  25. Blewitt, G.; Hammond, W.C.; Kreemer, C. Harnessing the GPS data explosion for interdisciplinary science. Eos 2018, 99, 485. [Google Scholar] [CrossRef]
  26. Duda, R.O.; Hart, P.E.; Stork, D.G. Pattern Classification; Wiley: Hoboken, NJ, USA, 2000. [Google Scholar]
  27. Marple, S.L., Jr. Digital Spectral Analysis with Applications; Prentice-Hall, Inc.: Englewood Cliffs, NJ, USA, 1987. [Google Scholar]
  28. Mallat, S. A Wavelet Tour of Signal Processing, 2nd ed.Academic Press: Cambridge, MA, USA.
  29. Donoho, D.L. and Johnstone, I. M., Adapting to unknown smoothness via wavelet shrinkage, J. Am. Stat. Assoc., 1995, 90, 1200–1224. [Google Scholar]
  30. Lyubushin, A. Low-Frequency Seismic Noise Properties in the Japanese Islands, Entropy 2021, 23, 474. 23. [CrossRef]
  31. Lyubushin, A. Investigation of the Global Seismic Noise Properties in Connection to Strong Earthquakes. Front. Earth Sci. 2022, 10, 905663. [Google Scholar] [CrossRef]
  32. Lyubushin, A. Spatial Correlations of Global Seismic Noise Properties. Appl. Sci. 2023, 13, 6958. [Google Scholar] [CrossRef]
  33. Jolliffe, I.T. Principal Component Analysis; Springer: Berlin/Heidelberg, Germany, 1986. [Google Scholar]
  34. Lyubushin, A.A. Long-range coherence between seismic noise properties in Japan and California before and after Tohoku mega-earthquake. Acta Geodaetica et Geophysica 2017, 52, 467–478. [Google Scholar] [CrossRef]
Preprints 82560 g001
Figure 1. (a) - network of GPS stations in the US West, 2009-2022; (b) is the density of GPS stations, estimated using a Gaussian kernel with a smoothing radius of 0.25 degrees on a regular grid of nodes with a size of 100×150.
Figure 1. (a) - network of GPS stations in the US West, 2009-2022; (b) is the density of GPS stations, estimated using a Gaussian kernel with a smoothing radius of 0.25 degrees on a regular grid of nodes with a size of 100×150.
Preprints 82560 g001
Preprints 82560 g002
Figure 2. The number of operational GPS stations in a window of the length 365 days with mutual shift 7 days.
Figure 2. The number of operational GPS stations in a window of the length 365 days with mutual shift 7 days.
Preprints 82560 g002
Preprints 82560 g003
Figure 3. Average probability densities of distribution of extreme values of 4 statistics of “active” and “passive” tremor of the earth’s surface for 3 components of the displacement vector. Because each property value is estimated over a 365-day window offset by 7 days, the average probability density estimate is obtained by averaging 548 probability density maps, each derived from 105 tremor property value maps.
Figure 3. Average probability densities of distribution of extreme values of 4 statistics of “active” and “passive” tremor of the earth’s surface for 3 components of the displacement vector. Because each property value is estimated over a 365-day window offset by 7 days, the average probability density estimate is obtained by averaging 548 probability density maps, each derived from 105 tremor property value maps.
Preprints 82560 g003
Preprints 82560 g004
Figure 4. Figure (4a) and (4b) are the results of averaging the weighted average probability densities of extreme values of 12 surface tremor parameters over a moving time window of 3 years for active and passive tremor. The red circles in (a) represent the positions of active tremor singular points forming 6 numbered clusters, the blue circles in (b) represent the positions of passive tremor singular points forming 3 numbered clusters. Figure (4c) and (4d) are histograms of the numbers of singular points for “active” and “passive” tremors belonging to different clusters.
Figure 4. Figure (4a) and (4b) are the results of averaging the weighted average probability densities of extreme values of 12 surface tremor parameters over a moving time window of 3 years for active and passive tremor. The red circles in (a) represent the positions of active tremor singular points forming 6 numbered clusters, the blue circles in (b) represent the positions of passive tremor singular points forming 3 numbered clusters. Figure (4c) and (4d) are histograms of the numbers of singular points for “active” and “passive” tremors belonging to different clusters.
Preprints 82560 g004
Preprints 82560 g005
Figure 5. Plots of changes in the weights of each of the 12 probability densities in “long” time windows of 1095 days with a shift of 7 days for “active” and “passive” tremor. The numbers indicate the average values of the weights, the 4 largest averages out of 12 variants for each type of tremor are highlighted in red.
Figure 5. Plots of changes in the weights of each of the 12 probability densities in “long” time windows of 1095 days with a shift of 7 days for “active” and “passive” tremor. The numbers indicate the average values of the weights, the 4 largest averages out of 12 variants for each type of tremor are highlighted in red.
Preprints 82560 g005
Preprints 82560 g006
Figure 6. Figure (6a) and (6b) - time sequences of longitudes and latitudes of singular points of tremor with a shift in the time window of 3 years for “active” (red) and “passive” (blue) tremor with cluster membership marks. Figure (6c) is the distance between the singular points of active and passive tremor in each time window of 3 years.
Figure 6. Figure (6a) and (6b) - time sequences of longitudes and latitudes of singular points of tremor with a shift in the time window of 3 years for “active” (red) and “passive” (blue) tremor with cluster membership marks. Figure (6c) is the distance between the singular points of active and passive tremor in each time window of 3 years.
Preprints 82560 g006
Preprints 82560 g007
Figure 7. Figure (7a) and (7b) - time sequences of values of the normalized maximum eigenvalue of the correlation matrix with a shift in the time window of 3 years for “active” (a) and “passive” (b) tremor.
Figure 7. Figure (7a) and (7b) - time sequences of values of the normalized maximum eigenvalue of the correlation matrix with a shift in the time window of 3 years for “active” (a) and “passive” (b) tremor.
Preprints 82560 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated