Skip to main content

Spectral-infinite-element simulations of gravity anomalies

Author(s): Gharti, Hom Nath; Tromp, Jeroen; Zampini, Stefano

Download
To refer to this page use: http://arks.princeton.edu/ark:/88435/pr1k35md15
Full metadata record
DC FieldValueLanguage
dc.contributor.authorGharti, Hom Nath-
dc.contributor.authorTromp, Jeroen-
dc.contributor.authorZampini, Stefano-
dc.date.accessioned2022-01-25T14:59:31Z-
dc.date.available2022-01-25T14:59:31Z-
dc.date.issued2018-08-07en_US
dc.identifier.citationGharti, Hom Nath, Jeroen Tromp, and Stefano Zampini. "Spectral-infinite-element simulations of gravity anomalies." Geophysical Journal International 215, no. 2 (2018): 1098-1117. doi:10.1093/GJI/GGY324.en_US
dc.identifier.issn0956-540X-
dc.identifier.urihttp://arks.princeton.edu/ark:/88435/pr1k35md15-
dc.description.abstractGravity anomalies induced by density heterogeneities are governed by Poisson’s equation. Most existing methods for modelling such anomalies rely on its integral solution. In this approach, for each observation point, an integral over the entire density distribution needs to be carried out, and the computational cost is proportional to the number of observation points. Frequently, such methods are sensitive to high density contrasts due to inaccurate resolution of the volume integral. We introduce a new approach which directly solves a discretized form of the Poisson/Laplace equation. The main challenge in our approach involves the unbounded nature of the problem, because the potential exists in all of space. To circumvent this challenge, we combine a mapped infinite-element approach with a spectral-element method. Spectral elements represent the domain of interest, and a single layer of infinite elements captures outer space. To solve the weak form of the Poisson/Laplace equation, we use Gauss–Legendre–Lobatto (GLL) quadrature in spectral elements inside the domain of interest. Outside the domain, we use Gauss–Radau quadrature in the infinite direction, and GLL quadrature in the other directions. We illustrate the efficiency and accuracy of our method by comparing calculated gravity anomalies for various density heterogeneities with corresponding analytical solutions. Finally, we consider a complex 3-D model of an ore mine, which consists of both positive and negative density anomalies.en_US
dc.format.extent1098 - 1117en_US
dc.language.isoen_USen_US
dc.relation.ispartofGeophysical Journal Internationalen_US
dc.rightsFinal published version. Article is made available in OAR by the publisher's permission or policy.en_US
dc.titleSpectral-infinite-element simulations of gravity anomaliesen_US
dc.typeJournal Articleen_US
dc.identifier.doidoi:10.1093/GJI/GGY324-
dc.identifier.eissn1365-246X-
pu.type.symplectichttp://www.symplectic.co.uk/publications/atom-terms/1.0/journal-articleen_US

Files in This Item:
File Description SizeFormat 
Spectral-infinite-element_simulations_gravity_anomalies.pdf7.48 MBAdobe PDFView/Download


Items in OAR@Princeton are protected by copyright, with all rights reserved, unless otherwise indicated.