Geomagnetic inverse modeling to estimate the volume of magnetization in areas with topographic variations

Abstract

A three-dimensional geomagnetic modeling is proposed to obtain information about a magnetization over topographic area from magnetic total force from land data survey. The model is divided into three-dimensional mesh grid of rectangular prisms that cover a finite depth below topographic elevation. The model parameter is assigned to each prism describing magnetization. A set of linear equations was formulated in terms of volume magnetization for each prism. The inverse calculation is obtained using the solution of over-determined problems. In this inverse calculation over three-dimensional mesh grid, the weighted matrix will be used to determine distribution of volume magnetization in the subsurface and keep the zero value of volume magnetization above the surface. A computer program has been tested with artificial topographic data and applied to digital elevation model of Bandung. The result of the inversion shows the comparison between synthetic model and inverse model to provide several analysis of algorithm performance. The high magnetization value can be characterize using inverse model up to 77% of volume area with average success recovery of model parameter value up to 56%.