3D Gravity Cross-Correlation Imaging for Large Scale Data Analysis: Application to the Crustal Structure of Iran

Document Type : Research


1 Ph.D. Student, Department of Mining and Metallurgical Engineering, Yazd University, Yazd, Iran

2 Associate Professor. Department of Mining and Metallurgical Engineering, Yazd University, Yazd, Iran


We propose the 3D gravity cross-correlation method to large scale data analyses as a fast analysis method to image the underground mass distribution. This method presents the cross-correlation product of the observed gravity anomaly (or its vertical gradient) and the calculated field due to an elementary mass contrast source. The cross-correlation product of the domain is used to highlight the zones of the highest probability of mass concentrations. First, some synthetic examples demonstrate the reliability and resolution of the method. The synthetic models discover different parameters of investigation space as space dimensions and densities. Tests with synthetic bodies show that the resultant correlation coefficients of the approach can delineate causative bodies in the subsurface. Finally, terrestrial gravity anomaly data of Iran is used to study the crustal structure and the Moho depth of Iran. The result is in a good agreement compared with other research studies of the domain. This technique took about five minutes to calculate the 3D gravity cross-correlation of the whole terrestrial gravity data set of Iran (25,937 data) a computer. Hence, it can easily be used repeatedly to monitor changes of gravity field.


Main Subjects

Inácio, P. and Gunter, B.C., 2010, A Sensitivity Study into Strapdown Airborne Gravimetry. Aerosp. Eng., vol. M.Sc., P. 125.
Alberts, B.A., 2009, Regional gravity field modeling using airborne gravity data. http:// repository.tudelft.nl/assets/uuid: 9c6c7ce7- 7227-4bab-8699- db32697c92a8/Alberts_2009_ phdthesis. pdf.
Braile, L.W., Keller, G.R. and Peeples, W.J., 1974, Inversion of gravity data for two-dimensional density distributions. J. Geophys. Res., 79(14), 2017–2021.
Bear, G.W., Al-Shukri, H.J. and Rudman, A.J., 1995, Linear inversion of gravity data for 3-D density distributions. Geophysics J., 60(5), 1354–1364.
Li, Y. and Oldenburg, D.W., 1998, 3-D inversion of gravity data. Geophysics, 63(1), 109–119, doi: 10.1190/1.1444302.
Talwani, M. and Ewing, M., 1960, Rapid computation of gravitational attraction of three-dimensional bodies of arbitrary shape. Geophysics, 25(1), 203–225.
Cordell, L. and Henderson, R.G., 1968, Iterative three-dimensional solution of gravity anomaly data using a digital computer. Geophysics, 33(4), 596–601.
Oldenburg, D.W., 1974, The inversion and interpretation of gravity anomalies. Geophysics, 39(4), 526–536.
Gómez-Ortiz, D. and Agarwal, B.N.P., 2005, 3DINVER. M: a MATLAB program to invert the gravity anomaly over a 3D horizontal density interface by Parker--Oldenburg’s algorithm. Comput. Geosci. J., 31(4), 513–520.
Chakravarthi, V. and Sundararajan, N., 2007, 3D gravity inversion of basement relief—A depth-dependent density approach. Geophysics, 72(2), I23-I32.
Patella, D., 1997, Introduction to ground surface self-potential tomography. Geophys. Prospect., 45(4), 653–681.
Mauriello, P. and Patella, D., 1999a, Principles of probability tomography for natural-source electromagnetic induction fields. Geophysics, 64(5), 1403–1417.
Mauriello, P. and Patella, D., 1999b, Resistivity anomaly imaging by probability tomography. Geophys. Prospect. J., 47(3), 411–429.
Mauriello, P. and Patella, D., 2001, Gravity probability tomography: a new tool for buried mass distribution imaging. Geophys. Prospect. J., 49(1), 1–12.
Alaia, R., Patella, D. and Mauriello, P., 2009, Imaging multipole self-potential sources by 3D probability tomography. Prog. Electromagn. Res. J., 14, 311–339.
Guo, L., Meng, X. and Shi, L., 2010, 3D correlation imaging of the vertical gradient of gravity data.  Geophys. Eng. J., 8(1), 6–12.
Plouff, D., 1976, Gravity and magnetic fields of polygonal prisms and application to magnetic terrain corrections. Geophysics J., 41(4), 727–741.
Berberian, M. and King, G.C.P., 1981, Towards a paleogeography and tectonic evolution of Iran: Reply. Can. J. Earth Sci., 18(11), 1764–1766.
Berberian, M., 1983, The southern Caspian: a compressional depression floored by a trapped, modified oceanic crust. Can. J. Earth Sci., 20(2), 163–183.
Stoneley, R., 1981, The geology of the Kuh-e Dalneshin area of southern Iran, and its bearing on the evolution of southern Tethys. J. Geol. Soc. London., 138(5), 509–526.
Richards, J.P., Wilkinson, D. and Ullrich, T., 2006, Geology of the Sari Gunay epithermal gold deposit, northwest Iran. Econ. Geol. J., 101(8), 1455–1496.
Bird, P., 1978, Finite element modeling of lithosphere deformation: the Zagros collision orogeny. Tectonophysics, 50(2–3), 307–336.
Page, W.D., Alt, J.N., Cluff, L.S. and Plafker, G., 1979, Evidence for the recurrence of large-magnitude earthquakes along the Makran coast of Iran and Pakistan. Tectonophysics, 52 (1–4), 533–547.
Dehegani, G.A. and Makris, J., 1983, The gravity field and crustal structure of Iran. Geodynamic project (Geotraverse) in Iran. Neues Jahrbuch Geologie und Palaeontologie, Abhandlungen.
Mokhtari, M., Farahbod, A. M., Lindholm, C., Alahyarkhani, M. and Bungum,  H., 2004, An approach to a comprehensive Moho depth map and crust and upper mantle velocity model for Iran. Iran. Int. J. Sci., 5(2), 223–244.
Taghizadeh-Farahmand, F., Afsari, N. and Sodoudi, F., 2015, Crustal thickness of Iran inferred from converted waves. Pure Appl. Geophys., 172(2), 309–331.
Mousavi, N. and Ebbing, J., 2018, Basement characterization and crustal structure beneath the Arabia--Eurasia collision (Iran): a combined gravity and magnetic study. Tectonophysics, 731, 155–171.
Shad Manaman, N., Shomali, H. and Koyi, H., New constraints on upper-mantle S-velocity structure and crustal thickness of the Iranian plateau using partitioned waveform inversion. Geophys. J. Int., 184(1), 247–267.
Radjaee, A., Rham, D., Mokhtari, M., Tatar, M., Priestley, K. and Hatzfeld, D., 2010, Variation of Moho depth in the central part of the Alborz Mountains, northern Iran. Geophys. J. Int., 181(1), 173–184.
Sodoudi, F., Yuan, X., Kind, R., Heit, B. and Sadidkhouy, A., 2009,  Evidence for a missing crustal root and a thin lithosphere beneath the Central Alborz by receiver function studies. Geophys. J. Int., 177(2), 733–742.
Afsari, N., Sodoudi, F., Farahmand, F.T. and Ghassemi, M.R., 2011, Crustal structure of northwest Zagros (Kermanshah) and Central Iran (Yazd and Isfahan) using teleseismic PS converted phases. J. Seismol., 15(2), 341–353.
Abdollahi, S., Ardestani,V.E.,  Zeyen, H. and Shomali, Z., 2018, Crustal and upper mantle structures of Makran subduction zone, SE Iran by combined surface wave velocity analysis and gravity modeling. Tectonophysics, 747, 191-210.
Abdollahi, S., Zeyen, H., Ardestani, V.E. and Shomali, Z., 2019, 3D joint inversion of gravity data and Rayleigh wave group velocities to resolve shear-wave velocity and density structure in the Makran subduction zone, southeast Iran. Journal of Asian Earth Sciences, 173, 275-290.
Abedi, M. and Oskooi, B., 2015, A combined magnetometry and gravity study across Zagros orogeny in Iran. Tectonophysics, 664, 164–175.