Ambient Noise Tomography, Green’s Function and Earthquakes


Green’s function is well-known, among others, in the application of ambient noise tomography methodologies that may demonstrate the potential of hydrocarbon entrapment in the study area. Here it is also shown to be of key importance in identifying the fractal dimension in the unified scaling law for earthquakes as well as in studying an explicit relationship of a future strong earthquake epicenter to the average earthquake potential score. Such studies are now in progress.

Share and Cite:

Varotsos, P. and Skordas, E. (2022) Ambient Noise Tomography, Green’s Function and Earthquakes. International Journal of Geosciences, 13, 1082-1088. doi: 10.4236/ijg.2022.1312055.

1. Introduction

Seismic waves from earthquakes and/or artificial sources travel through the Earth carrying information about the Earth’s subsurface structure and properties. Aki [1] and Claerbout [2] contributed a lot to the development of the technique for seismic applications. A signal at location A can be cross correlated with a signal at location B to reproduce a virtual source-receiver pair. It has been demonstrated (e.g., [3]) that this cross correlation can reproduce the surface waves of the Earth’s impulse response, or the Green’s function, as if triggered by a point source. By studying the dispersion relation of these surface waves between multiple pairs of stations, surface wave tomography is possible [4].

It is the scope of this work to reveal the decisive importance of Green’s function in the study of earthquakes and in particular along the following two directions that are described in detail below, in Sections 2 and 3: First, the determination of the fractal dimension of the earthquake epicenters projected onto the Earth’s surface, e.g. Fig. 7d of the work by Lei and Kusunose [5]. Second, the interrelation between the epicenter of a future strong earthquake and the average earthquakes potential score (EPS) maps.

Before proceeding, we emphasize that scale invariance methods such as fractal dimension, correlation dimension and multifractal spectrum have potential applications in emerging fields of science, engineering, and seismology. Numerous studies have been undergoing to study the applications in seismology, in particular, see references [6] - [24].

2. Additional Usefulness of Green’s Function

2.1. Fractal Dimension of the Earthquake Epicenters Projected onto the Earth’s Surface

We have recently [25] shown that an estimate of the epicenter location of a future strong earthquake can be obtained by combining a new analysis of seismicity termed natural time analysis [26] [27] [28] [29] with earthquake networks based on similar activity patterns [30] and earthquake nowcasting [31] - [38]. This is based on the construction of average EPS maps.

A simple model was developed to understand the process of self-consistently averaging EPS e.g., see Fig. 5 of Reference [25]. To keep this model simple, we assume that all EQs occur practically in the center O of circular region of radius R. After a long enough period without a large EQ, the EPS takes the value of unity for all points inside and, since we assumed that no other EQ takes place, the EPS will be zero outside it. Hence, the value of 〈EPS〉 at a point P lying at a distance d away from O simplifies to:

EPS = A ( d , R ) π R 2 (1)

where A(d,R) is the overlapping area of the two intersecting circles, equaled:

A ( d , R ) = 2 R 2 arccos ( d 2 R ) d R 2 ( d 2 ) 2 (2)

Hence, the value of “EPS” in polar coordinates (ρ, θ) with center at O equaled to:

f ( ρ ) = 2 π arccos ( ρ 2 R ) ρ π R 1 ( ρ 2 R ) 2 (3)

for ρ < 2R and zero otherwise. If we now calculate the mean value of “EPS”:

m ( R , R ) = 0 R 0 2 π f ( ρ ) ρ d θ d ρ π R 2 = 2 0 R f ( ρ ) ρ d ρ R 2 (4)

in a circular region of radius R' > 2R, we obtain that:

m ( R , R ) = ( R R ) 2 (5)

by virtue of the integrals:

0 2 R arccos ( ρ 2 R ) ρ d ρ = π R 2 (6)

0 2 R ρ 2 d ρ 1 ( ρ 2 R ) 2 = π R 3 2 (7)

The validity of the latter Equation is verified by assuming that EQs occur according to the aforementioned simple model and, using the computer programs used for the calculation of the actual 〈EPS〉 maps shown in Fig. 4 of Reference [25], the numerically found mn(R,R') has the form:

m n ( R , R ) = ( R R ) d f (8)

When the quantity mn(R,R') was studied for the 〈EPS〉 maps of Fig. 4 of Reference [25] we find, df = 1.32. This value of df differed only slightly from the value df ≈ 1.2, which Bak et al. [39] found to describe the fractal dimension of the location of epicenters projected onto the surface of the Earth in a unified scaling law obeyed by the distribution of waiting times between EQs occurring in California and ranging from tens of seconds to tens of years.

2.2. On the Relation between Average EPS Maps and the Epicenter of a Future Strong Earthquake

We now proceed to another usefulness of Green’s function. A self-consistent method of producing average EPS maps, also written 〈EPS〉 maps, using a radius R has been suggested and applied to the Eastern Mediterranean area in [25]. To construct such a map, one first estimates EPS for disks of radius R at the points ( x i j , y i j ) of a lattice to obtain EPSij and then averages for each point ( x i 0 j 0 , y i 0 j 0 ) the estimated EPS values within the same radius R, i.e.,

EPS ( x i 0 j 0 , y i 0 j 0 ) = 1 N i , j d ( x i 0 j 0 , y i 0 j 0 ; x i j , y i j ) R EPS i j (9)

where the summation is restricted to the lattice points whose distance d ( x i 0 j 0 , y i 0 j 0 ; x i j , y i j ) from the observation point is smaller than or equal to R, and N stands for the number of lattice points included in the sum.

It has been shown that the study of 〈EPS〉 close to the epicenters of forthcoming strong EQs exhibits a logarithmic dependence on R, reminiscent of the Green’s function of the Poisson equation in two dimensions, while the mean value EPS ¯ of 〈EPS〉 over all the lattice points scales with R as a power law with an exponent df, i.e., EPS ¯ R d f , see Section 3 and Equation (12) of Reference [25].

A clear relation between such made 〈EPS〉 maps and the epicenter of an impending strong EQ has been observed in the respective regional studies of Reference [25].

3. Conclusions

Green’s function, beyond its usefulness in the application of ambient noise tomography methodologies that may demonstrate the potential of hydrocarbon entrapment in the study area, is shown to be of usefulness in the study of earthquakes along the following two directions:

1) In the identification of the fractal dimension of the earthquake epicenters projected onto the Earth’s surface in the unified scaling law for earthquakes obtained by Bak et al. [39], and

2) The two-dimensional Green’s function has recently been recovered when investigating the relationship of a future strong earthquake epicenter to the average earthquake potential score.

Such additional studies are now in progress in various areas due to their importance.

This importance is further strengthened by the fact that very recent aspects (see for example Reference [40]) are focused on the spatio-temporal variations of the correlation fractal dimension for earthquakes with magnitude M equal or greater than 2.5 in southern and Baja California to ascertain the incidence of seismic precursors before strong earthquakes [40].

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.


[1] Aki, K. (1957) Space and Time Spectra of Stationary Stochastic Waves, with Special Reference to Microtremors. Bulletin of the Earthquake Research Institute, 35, 415-456.
[2] Claerbout, J.F. (1968) Synthesis of a Layered Medium from ITS Acoustic Transmission Response. Geophysics, 33, 264-269.
[3] Snieder, R. and Wapenaar, K. (2010) Imaging with Ambient Noise. Physics Today, 63, 44-49.
[4] Shapiro, N.M., Campillo, M., Stehly, L. and Ritzwoller, M.H. (2005) High-Resolution Surface-Wave Tomography from Ambient Seismic Noise. Science, 307, 1615-1618.
[5] Lei, X. and Kusunose, K. (1999) Fractal Structure and Characteristic Scale in the Distributions of Earthquake Epicentres, Active Faults and Rivers in Japan. Geophysical Journal International, 139, 754-762.
[6] Ogata, Y., Jones, L.M. and Toda, S. (2003) When and Where the Aftershock Activity Was Depressed: Contrasting Decay Patterns of the Proximate Large Earthquakes in Southern California. Journal of Geophysical Research: Solid Earth, 108, Article No. 2318.
[7] Radulian, M. and Trifu, C.I. (1991) Would It Have Been Possible to Predict the 30 August 1986 Vrancea Earthquake? Bulletin of the Seismological Society of America, 81, 2498-2503.
[8] Matcharashvili, T., Chelidze, T. and Javakhishvili, Z. (2000) Nonlinear Analysis of Magnitude and Interevent Time Interval Sequences for Earthquakes of the Caucasian Region. Nonlinear Processes in Geophysics, 7, 9-20.
[9] Legrand, D. (2002) Fractal Dimensions of Small, Intermediate, and Large Earthquakes. Bulletin of the Seismological Society of America, 92, 3318-3320.
[10] Nakaya, S. (2005) Fractal Properties of Seismicity in Regions Affected by Large, Shallow Earthquakes in Western Japan: Implications for Fault Formation Processes Based on a Binary Fractal Fracture Network Model. Journal of Geophysical Research: Solid Earth, 110, B01310.
[11] Roy, P.N.S. and Ram, A. (2006) A Correlation Integral Approach to the Study of 26 January 2001 Bhuj Earthquake, Gujarat, India. Journal of Geodynamics, 41, 385-399.
[12] Roy, P.N.S. and Padhi, A. (2007) Multifractal Analysis of Earthquakes in the South-Eastern Iran-Bam Region. Pure and Applied Geophysics, 164, 2271-2290.
[13] Roy, P.N.S. and Mondal, S.K. (2012) Identification of Active Seismicity by Fractal Analysis for Understanding the Recent Geodynamics of Central Himalaya. Journal of the Geological Society of India, 79, 353-360.
[14] Roy, P.N.S. and Mondal, S.K. (2012) Multifractal Analysis of Earthquakes in Kumaun Himalaya and Its Surrounding Region. Journal of Earth System Science, 121, 1033-1047.
[15] Kumar, A., Rai, S.S., Joshi, A., Mittal, H., Sachdeva, R., Kumar, R. and Ghangas, V. (2013) The b-Value and Fractal Dimension of Local Seismicity around Koyna Dam (India). Earthquake Science, 26, 99-105.
[16] Roy, P.N.S., Chowdhury, S., Sarkar, P. and Mondal, S.K. (2015) Fractal Study of Seismicity in Order to Characterize the Various Tectonic Blocks of North-East Himalaya, India. Natural Hazards, 77, 5-18.
[17] Naimi-Ghassabian, N., Khatib, M.M., Nazari, H. and Heyhat, M.R. (2018) Regional Variations and Earthquake Frequency-Magnitude Distribution and Fractal Dimension in the North of Central-East Iran Blocks (NCEIB). Arabian Journal of Geosciences, 11, 243-264.
[18] Hamdache, M., Henares, J., Peláez, J.A. and Damerdji, Y. (2019) Fractal Analysis of Earthquake Sequences in the Ibero-Maghrebian Region. Pure and Applied Geophysics, 176, 1397-1416.
[19] Hui, C., Cheng, C., Ning, L. and Yang, J. (2020) Multifractal Characteristics of Seismogenic Systems and b Values in the Taiwan Seismic Region. International Journal of Geo-Information, 9, Article No. 384.
[20] Sukmono, S. (2001) New Evidence on the Fractal Pattern of Sumatra Fault Seismicity and Its Possible Application to Earthquake Prediction. Bulletin of the Seismological Society of America, 91, 870-874.
[21] Okubo, P.G. (1987) Fractal Geometry in the San Andreas Fault System. Journal of Geophysical Research: Solid Earth, 92, 345-355.
[22] Lee. H.-K. and Schwarch, H.P. (1995) Fractal Clustering of Fault Activity in California. Geology, 23, 377-380.<0377:FCOFAI>2.3.CO;2
[23] Robertson, M.C., Sammis, C.G., Sahimi, M. and Martin, A.J. (1995) Fractal Analysis of Three-Dimensional Spatial Distributions of Earthquakes with a Percolation Interpretation. Journal of Geophysical Research: Solid Earth, 100, 609-620.
[24] Ghosal, A., Ghosh, U. and Kayal, J.R. (2012) A Detailed b-Value and Fractal Dimension Study of the March 1999 Chamoli Earthquake (Ms 6.6) Aftershock Sequence in Western Himalaya. Geomatics, Natural Hazards and Risk, 3, 271-278.
[25] Varotsos, P.K., Perez-Oregon, J., Skordas, E.S. and Sarlis, N.V. (2021) Estimating the Epicenter of an Impending Strong Earthquake by Combining the Seismicity Order Parameter Variability Analysis with Earthquake Networks and Nowcasting: Application in Eastern Mediterranean. Applied Sciences, 11, Article No. 10093.
[26] Varotsos, P.A., Sarlis, N.V. and Skordas, E.S. (2002) Long-Range Correlations in the Electric Signals That Precede Rupture. Physical Review E, 66, Article ID: 011902.
[27] Varotsos, P.A., Sarlis, N.V. and Skordas, E.S. (2011) Natural Time Analysis: The New View of Time. Precursory Seismic Electric Signals, Earthquakes and Other Complex Time-Series. Springer, Berlin.
[28] Varotsos, P.A., Sarlis, N.V., Skordas, E.S., Uyeda, S. and Kamogawa, M. (2011) Natural Time Analysis of Critical Phenomena. Proceedings of the National Academy of Sciences of the United States of America, 108, 11361-11364.
[29] Sarlis, N.V., Skordas, E.S., Varotsos, P.A., Nagao, T., Kamogawa, M., Tanaka, H. and Uyeda, S. (2013) Minimum of the Order Parameter Fluctuations of Seismicity before Major Earthquakes in Japan. Proceedings of the National Academy of Sciences of the United States of America, 110, 13734-13738.
[30] Tenenbaum, J.N., Havlin, S. and Stanley, H.E. (2012) Earthquake Networks Based on Similar Activity Patterns. Physical Review E, 86, Article ID: 046107.
[31] Rundle, J.B., Turcotte, D.L., Donnellan, A., Grant Ludwig, L., Luginbuhl, M. and Gong, G. (2016) Nowcasting Earthquakes. Earth and Space Science, 3, 480-486.
[32] Rundle, J.B., Luginbuhl, M., Giguere, A. and Turcotte, D.L. (2018) Natural Time, Nowcasting and the Physics of Earthquakes: Estimation of Seismic Risk to Global Megacities. Pure and Applied Geophysics, 175, 647-660.
[33] Luginbuhl, M., Rundle, J.B. and Turcotte, D.L. (2018) Statistical Physics Models for Aftershocks and Induced Seismicity. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 377, Article ID: 20170397.
[34] Luginbuhl, M., Rundle, J.B. and Turcotte, D.L. (2018) Natural Time and Nowcasting Earthquakes: Are Large Global Earthquakes Temporally Clustered? Pure and Applied Geophysics, 175, 661-670.
[35] Rundle, J.B., Giguere, A., Turcotte, D.L., Crutchfield, J.P. and Donnellan, A. (2019) Global Seismic Nowcasting with Shannon Information Entropy. Earth and Space Science, 6, 191-197.
[36] Rundle, J.B. and Donnellan, A. (2020) Nowcasting Earthquakes in Southern California with Machine Learning: Bursts, Swarms, and Aftershocks May Be Related to Levels of Regional Tectonic Stress. Earth and Space Science, 7, e2020EA001097.
[37] Rundle, J., Stein, S., Donnellan, A., Turcotte, D.L., Klein, W. and Saylor, C. (2021) The Complex Dynamics of Earthquake Fault Systems: New Approaches to Forecasting and Nowcasting of Earthquakes. Reports on Progress in Physics, 84, Article ID: 076801.
[38] Rundle, J.B., Donnellan, A., Fox, G. and Crutchfield, J.P. (2022) Nowcasting Earthquakes by Visualizing the Earthquake Cycle with Machine Learning: A Comparison of Two Methods. Surveys in Geophysics, 43, 483-501.
[39] Bak, P., Christensen, K., Danon, L. and Scanlon, T. (2002) Unified Scaling Law for Earthquakes. Physical Review Letters, 88, Article ID: 178501.
[40] Chandriyan, H., Reddy, R. and Roy, R.N.S. (2022) Numerical Precursory Study on Strong Earthquakes in Southern and Baja California. Geosystems and Geoenvironment, 1, Article 100066.

Copyright © 2023 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.