Recovering Magnetization of Rock Formations by Jointly Inverting Airborne Gravity Gradiometry and Total Magnetic Intensity Data

Conventional 3D magnetic inversion methods are based on the assumption that there is no remanent magnetization, and the inversion is run for magnetic susceptibility only. This approach is well-suited to targeting mineralization; however, it ignores the situation where the direction of magnetization of the rock formations is different from the direction of the induced magnetic field. We present a novel method of recovering a spatial distribution of magnetization vector within the rock formation based on joint inversion of airborne gravity gradiometry (AGG) and total magnetic intensity (TMI) data for a shared earth model. Increasing the number of inversion parameters (the scalar components of magnetization vector) results in a higher degree of non-uniqueness of the inverse problem. This increase of non-uniqueness rate can be remedied by joint inversion based on (1) Gramian constraints or (2) joint focusing stabilizers. The Gramian constraints enforce shared earth structure through a correlation of the model gradients. The joint focusing stabilizers also enforce the structural similarity and are implemented using minimum support or minimum gradient support approaches. Both novel approaches are applied to the interpretation of the airborne data collected over the Thunderbird V-Ti-Fe deposit in Ontario, Canada. By combining the complementary AGG and TMI data, we generate jointly inverted shared earth models that provide a congruent image of the rock formations hosting the mineral deposit.