ITERATIVE THREE‐DIMENSIONAL SOLUTION OF GRAVITY ANOMALY DATA USING A DIGITAL COMPUTER

Abstract

Given gridded gravity anomaly values and certain limiting restrictions on the mass distribution, a three‐dimensional structural model can be calculated automatically from gravity anomaly data by successive approximations. The causative body is optionally assumed to be either flat‐topped, flat‐bottomed, or symmetrical about a horizontal plane. Density contrast and the position of a horizontal plane delimiting either the top, base, or midsection of the causative body must be specified. A first approximation of structure is obtained by means of the Bouguer slab relationship. The gravity field of this first model is calculated and at each grid point the ratio of observed to calculated gravity is used to modify the first structural model, thus leading to a second approximation of structure. The process is iterated until a satisfactory agreement between observed and calculated gravity is obtained.