Interfaces characterized by a density contrast are common models to interpret gravity data. Two key parameters of the interface inversion are the density contrast and constant offset, which are free parameters. We propose to iteratively estimate the two free parameters by a new method: the iterative-ratio (ITRAT) method. The method is based on the linear approximate invariance of depth and free parameters and on the availability of several depth constraints. The free parameters are estimated iteratively to update the interface depth. By processing the synthetic Moho model, we find that the ITRAT method generally converges after two iterations with good accuracy for the free-parameter estimates and Moho depth. Random tests for the constant density model show that the accuracy of ITRAT method is slightly affected by the noise of gravity data and by the number of constraints, but it is sensitive to the uncertainty of constraints. For the variable density model, the ITRAT method will estimate a group of constant density contrast and constant offset to obtain a better Moho model, which has the lower depth error than the model with a mean density contrast. The new method is finally demonstrated by mapping the Moho of the Santos basin (Brazil) according to the depth constraints and regional gravity anomalies.
Read full abstract