Recent developments in large-scale geophysical inversions made it possible to invert the results of entire airborne geophysical surveys over large areas into 3D models of the subsurface. This paper presents the methods for and results of the interpretation of the data acquired by a multiphysics airborne geophysical survey in Saudi Arabia. The project involved the acquisition, processing, and interpretation of airborne electromagnetic, gravity, and magnetic geophysical data over an 8000 square kilometer area. All the collected data were carefully analyzed and inverted in 3D models of the corresponding physical properties of the subsurface, including 3D density, magnetization vector, and conductivity models. This paper summarizes the interpretation of all geophysical data sets collected during the field airborne survey. The goal of the paper is to demonstrate how the advanced 3D modeling and inversion methods can be effectively used for interpretation of multiphysics airborne survey data and to study and analyze the potential of the survey area for natural resource exploration in Saudi Arabia.

The Saudi Arabian Glass Earth Pilot Project is a geophysical exploration program to explore the upper crust of the Kingdom for minerals, groundwater, and geothermal resources as well as strictly academic investigations. The project began with over 8000 km² of green-field area. Airborne geophysics including electromagnetic (EM), magnetics, and gravity were used to develop several high priority targets for ground follow-up. Based on the results of airborne survey, a spectral induced polarization (SIP) survey was completed over one of the prospective targets. The field data were collected with a distributed array system, which has the potential for strong inductive coupling. This was examined in a synthetic study, and it was determined that with the geometries and conductivities in the field survey, the inductive coupling effect may be visible in the data. In this study, we also confirmed that time domain is vastly superior to frequency domain for avoiding inductive coupling, that measuring decays from 50 ms to 2 s allow discrimination of time constants from 1 ms to 5 s, and the relaxation parameter C is strongly coupled to intrinsic chargeability. We developed a method to fully include all 3D EM effects in the inversion of induced polarization (IP) data. The field SIP data were inverted using the generalized effective-medium theory of induced polarization (GEMTIP) in conjunction with an integral equation-based modeling and inversion methods. These methods can replicate all inductive coupling and EM effects, which removes one significant barrier to inversion of large bandwidth spectral IP data. The results of this inversion were interpreted and compared with results of drill hole set up in the survey area. The drill hole intersected significant mineralization which is currently being further investigated. The project can be considered a technical success, validating the methods and effective-medium inversion technique used for the project.