An inversion scheme to trace three-dimensional density interfaces from thei
r gridded gravity anomalies is developed. The scheme requires the gravity a
nomaly values, profile and station spacings as the input, besides densities
of the underlying and overlying materials and the mean or undisturbed dept
h to the interface. The scheme calculates the initial values of the depths
to the interface below the internal grid points, and modifies them iterativ
ely until a best fit is achieved between the observed and calculated anomal
y values.
The gravity anomalies of the density interface are generated by equating th
e material below the interface to a series of juxtaposing rectangular block
s, one centered below each of the internal gridded anomaly points. It is as
sumed that the interface flattens out to its undisturbed depth well ahead o
f the boundaries of the sampled area. The values of initial thicknesses of
the prisms are calculated by equating the anomaly at a station to the algeb
raic sum of products of vertical gradients of gravity effects of the prisms
and their thicknesses. The differences between the observed and calculated
anomalies are then used to improve the thicknesses of the prisms by equati
ng these differences to the algebraic sum of the products of the vertical g
radients of gravity effects of prisms and increments to their thicknesses.
The computer program, named GRAV3DIN coded in FORTRAN 77 and used to invert
gravity anomalies of three-dimensional density interfaces, is presented. (
C) 1999 Elsevier Science Ltd. All rights reserved.