Inversion of gravity anomalies of three-dimensional density interfaces

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.