Airborne Magnetic Data Research Articles

3‐D Magnetic Unstructured Inversion

AbstractMagnetic inversion methods based on structured grids have been used extensively in geomagnetic research. Nevertheless, structured grids limit the modeling capability and accuracy of models with arbitrary geometries. To address these challenges, 3‐D magnetic numerical forward modeling and inversion methods using an unstructured tetrahedral grid based on a partial differential equation (PDE) framework are proposed. The methods are derived from Maxwell's partial equations and constructed using the finite element method. Arbitrary undulating topographies, complex geometries, and demagnetization effects can be represented exactly, leading to high‐accuracy forward modeling and inversion solutions. The proposed methods are suitable for both local‐ and global‐scale magnetic data interpretation, for example, in mineral exploration and tectonic research. A synthetic example and real airborne magnetic survey data collected from Mount Iliamna, Alaska, USA, are tested. The results demonstrate that the novel methods significantly improve the capability and accuracy of magnetic data interpretation.

3D Inversion of Airborne Magnetic Data and Validation of it by Using Sourk Mine Drilling Data

The airborne magnetic method is one of the most significantly applied approaches in mineral exploration. This magnetic method has been widely used in recent years due to its extensive coverage, lower cost than other methods and coverage of impassable areas and forests. For many years in Iran, the magnetic method has been utilized by authorities to survey ore deposits and to prepare structural and geological maps. This paper presents part of the Dehshir area exploration data of airborne surveys by helicopter with flight line spacing of 250 meters and flight height of 40 meters conducted by Geological Survey and Mineral Exploration of Iran. Furthermore, among the potential areas in the Dehshir area is the studied hereby area of Block C2 of Sourk mine. This region in Yazd province with a 4.1 km2 area is located at 95 km northwest of Yazd city. The purpose of this paper is to present a three-dimensional model of magnetic data and its validation using drilling data. The three-dimensional modeling of magnetic data has been processed by Inversion method. Next, the locations of boreholes were proposed for drilling. Then, the boreholes data is used to evaluate the geophysical method model by comparing the three-dimensional modeling of magnetic data and three-dimensional modeling of the proposed drilled boreholes data.

Structural and lithological interpretation of aero-geophysical data in parts of the Lower Benue Trough and Obudu Plateau, southeast Nigeria

High resolution airborne magnetic and radiometric data were used to map probable sites for polymetallic-magmatic hydrothermal deposits in the study area. The first and second vertical derivatives, horizontal gradient, analytic signal and tilt angle derivative were used to enhance magnetic data in order to determine the location of short wavelength anomalies, magnetic sources and geologic boundaries. Also, to understand and detect radiometric anomalies controlling mineralization in host rocks, enhancement operations involving potassium (K), equivalent thorium (Th), equivalent uranium (U) concentrations and ternary imageries were carried out. All these operations enabled the delineation of responses associated with mineralization, lithology and geologic structures. Generally, lineaments trend in the NE-SW, NNE-SSW, NW-SE and E-W directions, and serve as potential pathway for hydrothermal fluid migration and mineralization. Radiometric anomalies from spectrometric imageries are caused by felsic minerals, hydrothermal rock alterations, widespread shales, residual clay, oxides and accessory minerals that constitute regolith of the area. The coincident high magnetic and potassium intensities sites in Anambra Bain (AB), Abakaliki Anticlinorium (AA), Ikom-Mamfe Rift (IMR) and Obudu Plateau (OP) are probable sites for polymetallic-magmatic hydrothermal deposits. Generally, the joint magnetic and radiometric results were able to establish link between lithology, geologic structures and hydrothermal alteration patterns.

Heat flow analyses over Bornu Basin and its environs, Northeast Nigeria, using airborne magnetic and radiometric data: implication for geothermal energy prospecting

Heat flow analyses over Bornu Basin and its environs, Northeast Nigeria, using airborne magnetic and radiometric data: implication for geothermal energy prospecting

Precambrian tectonic inheritance control of the NE Brazilian continental margin revealed by Curie point depth estimation

An investigation of Curie point depths (CPD) based on spectral analysis of airborne magnetic data was carried out in the NE Brazilian continental margin. The studied region represents a narrow hyper-extended margin with three sedimentary basins. Regional geothermal gradient and heat flow were also calculated. CPD results were integrated with interpretation of 2D deep seismic data and with estimated isostatic Moho depths. The results reveal that the narrow hyper-extended crust is 150 km wide in the southern sector and 80 km wide in the north, with a narrow ocean-continental transition (OCT) zone that varies from 50 km wide in the south sector to 30 to 20 km wide in the north. The CPD isotherm showed the strong influence of the three main continental blocks of Borborema ́s Shield in the tectonic evolution of the three marginal basins. The CPD analysis corroborated models provided by gravimetric data and successfully demonstrated the sharp control of basement compartments on the thermal properties of the marginal basins domains

Using self-organizing maps in airborne geophysical data for mapping mafic dyke swarms in NE Brazil

Mapping widespread dyke swarms shed light on the tectonic processes that culminate in continental fragmentation, especially in the early phases of the crustal extension and magma emplacement. Airborne magnetic data are an effective geophysical tool to clarify the extent of dyke swarms in continental-scale areas due to the expressive magnetic contrast between magmatic bodies and host rocks. However, many geological features display similar magnetic patterns, making the qualitative interpretation of magnetic anomaly maps quite subjective and ambiguous. To improve and optimize the predictive mapping of dyke swarms, this research presents a multivariate analysis of airborne geophysical surveys, applying a Self-Organizing Maps (SOM) approach using two magnetic and three gamma-spectrometric variables. This semiautomatic technique has allowed an integrated spatial analysis and the classification of the lithological units based on their magnetic and gamma-spectrometric signatures. We applied the SOM method to investigate a set of mafic dyke swarms that intruded in the Neoproterozoic Borborema Province (BP), Parnaíba Basin (PB) and the São Francisco Craton (SFC) in NE Brazil. These dykes are part of a large magmatic event associated with the continental breakup, which formed the Equatorial Atlantic Ocean in the Early Cretaceous, denominated EQUAMP. First, the SOM hyperparameters were defined by running the algorithm in a representative area in the central part of the BP. This training test area worked as a SOM template through which all data from the area were processed. The SOM analysis identified seven different populations, according to responses found in the five geophysical input variables. Two of these populations were associated with the mafic dykes, reducing the ambiguity of the magnetic anomaly interpretation. These populations represent high SOM quantization error, which means that these groups are the most anomalous data, evidenced in airborne magnetic data. These results were checked during fieldwork, revealing that dyke swarms occur more widely than were previously known, throughout the BP intruding the SFC, and showing some occurrences embedded in the Paleozoic sedimentary infill of the eastern border of the PB.

Stable Downward Continuation of Airborne Potential Field Geophysical Data: an Investigation of Stabilizer Family

Attenuation of the signal received from sources causing anomalies and reduction of data resolution are the negative features of airborne measurements. Using a stable downward continuation method is a practical way to address these shortcomings. In this study, we investigated the efficiency of various stabilizers in achieving stable downward continued data. The purpose of this study is to select the most appropriate stabilizer(s) for this operation. We examined the various stabilizing functions by introducing them into the Tikhonov regularization problem. The results of research on synthetic airborne gravity and magnetic data showed that βL1 (the other definition of L1 norm) and SM (the smoothest model) stabilizers have the potential to be used in the stable implementation of the downward continuation method. These stabilizers performed better than the other in the three comparisons, including visual, quantitative (RMS error), and graphical comparisons. Also, by examining the airborne magnetic data related to the Esfordi district in Yazd province (Iran), it was found that in general the βL1 stabilizer is more suitable than the other stabilizing functions studied in this research.

Model-Based Probabilistic Inversion Using Magnetic Data: A Case Study on the Kevitsa Deposit

Structural geological models are widely used to represent relevant geological interfaces and property distributions in the subsurface. Considering the inherent uncertainty of these models, the non-uniqueness of geophysical inverse problems, and the growing availability of data, there is a need for methods that integrate different types of data consistently and consider the uncertainties quantitatively. Probabilistic inference provides a suitable tool for this purpose. Using a Bayesian framework, geological modeling can be considered as an integral part of the inversion and thereby naturally constrain geophysical inversion procedures. This integration prevents geologically unrealistic results and provides the opportunity to include geological and geophysical information in the inversion. This information can be from different sources and is added to the framework through likelihood functions. We applied this methodology to the structurally complex Kevitsa deposit in Finland. We started with an interpretation-based 3D geological model and defined the uncertainties in our geological model through probability density functions. Airborne magnetic data and geological interpretations of borehole data were used to define geophysical and geological likelihoods, respectively. The geophysical data were linked to the uncertain structural parameters through the rock properties. The result of the inverse problem was an ensemble of realized models. These structural models and their uncertainties are visualized using information entropy, which allows for quantitative analysis. Our results show that with our methodology, we can use well-defined likelihood functions to add meaningful information to our initial model without requiring a computationally-heavy full grid inversion, discrepancies between model and data are spotted more easily, and the complementary strength of different types of data can be integrated into one framework.

Towards an understanding of mineral systems – Contributions from magnetotelluric data from the Fennoscandian Shield in northern Sweden

Magnetotelluric data from 104 new stations in northern Sweden together with previously available MaSca project data are analysed with respect to a hypothesised genetic relation between mineral occurrences and commonly reported deep crustal electrically conductive structures in Precambrian Shield areas. A total of 165 stations from the north-eastern part of the Norrbotten ore province are modelled from an area characterised as having an Archaean basement with later reworking during the Paleoproterozoic. The data cover an area of 200 × 250 km2 and provide a significant improvement of spatial resolution compared to previously available data. The derived 3D model reveals the presence of strong crustal conductors with conductance of more than 3000 S at depth of tens of kilometres within a generally resistive crust. A subdivision into a central more conductive domain bounded towards north and south by resistive domains is revealed. The boundaries between these three domains are oriented NW-SE. The southern boundary coincides roughly with the boundary between Archaean and Paleoproterozoic crust defined by εNd data. The northern boundary coincides with an interpreted 2.4–2.3 Ga coastline of a marine rift basin, which is spatially defined by the occurrence of the Kovo group rocks. The location of major conductors within the central domain correlates with the presence of mapped Paleoproterozoic metasupracrustal belts. A significant part of the middle crust conductors is elongated in directions that coincide with major deformation zones that have been mapped from airborne magnetic data and geological fieldwork. Some of these conductors have near-surface expression where they spatially correlate with the location of known mineralisation. In these cases, the conductive structures may mark the pathways for metal-bearing hydrothermal fluids. Interpretations involving enhanced quantities of magnetite or graphite in the middle crust are discussed but the magnetotelluric data alone are not conclusive with respect to cause of observed high conductivities.

Low-latitude reduction-to-the-pole and upward continuation between arbitrary surfaces based on the partial differential equation framework

SUMMARY Magnetic surveys conducted in complex conditions, such as low magnetic latitudes, uneven observation surfaces, or above high-susceptibility sources, pose significant challenges for obtaining stable solutions for reduction-to-the-pole (RTP) and upward-continuation processing on arbitrary surfaces. To tackle these challenges, in this study, we propose constructing an equivalent-susceptibility model based on the partial differential equation (PDE) framework in the space domain. A multilayer equivalent-susceptibility method was used for RTP and upward-continuation operations, thus allowing for application on undulating observation surfaces and strong self-demagnetization effect in a non-uniform mesh. A novel positivity constraint is introduced to improve the accuracy and efficiency of the inversion. We analysed the effect of the depth-weighting function in the inversion of equivalent susceptibility for RTP and upward-continuation reproduction. Iterative and direct solvers were utilized and compared in solving the large, sparse, non-symmetric and ill-conditioned system of linear equations produced by PDE-based equivalent-source construction. Two synthetic models were used to illustrate the efficiency and accuracy of the proposed method in processing both ground and airborne magnetic data. Aeromagnetic and ground data collected in Brazil at a low magnetic latitude region were used to validate the proposed method for processing RTP and upward-continuation operations on magnetic data sets with strong self-demagnetization.

Application of Airborne Electromagnetics and Magnetics for Mineral Exploration in the Baishiquan–Hongliujing Area, Northwest China

Airborne electromagnetics is an effective and efficient exploration tool in shallow mineral exploration for its high efficiency and low cost. In 2016, airborne electromagnetic and airborne magnetic surveys have been carried out at the border of Xinjiang Uygur Autonomous Region and Gansu Province, the Northwest China. With an integrated system, the airborne electromagnetics and airborne magnetic data were collected simultaneously by AreoTEM-IV system from Aeroquest International Limited in Vancouver, BC, Canada, and the CS3 Cesium Vapor magnetometer from Scintrex in Concord, ON, Canada. About 3149 line-km of both data with 250 m line space were acquired. After data processing, the comprehensive analysis and interpretation of resistivity and magnetic anomalies has been carried out to infer lithological structure and outline the potential ore deposits. Verified by the ground surveys, seven outlined anomalies are consistent with the known ore sites, and one new gold deposit and several mineralization clues were found. The prospective reserves of gold are expected to exceed 10 tons. Besides, some prospecting target areas were outlined as the possible locations of copper–nickel deposits. The successful case shows the airborne magnetic data accords with geological structures, and the airborne electromagnetic method is effective in finding metal mineral resources, which can help to quickly identify potential ore targets with no surface outcrop.

Inversion of magnetic data using deep neural networks

As a novel tool, deep learning is used to solve complex problems in the real world and has been successfully applied to invert for seismic and electromagnetic data. In this study, we propose to use deep neural networks (DNNs) to recover the distribution of the physical properties of buried magnetic orebodies from the surface and airborne magnetic anomaly data. This approach is based on data training instead of prior-knowledge assumptions used in traditional inversion methods. Once our generalized network is established, the computing time to predict new magnetic data using this method can be significantly reduced. By implementing the forward modeling of different types of synthetic physical property models, we obtained enough datasets to train a DNN model so that the network can establish a nonlinear mapping directly from magnetic anomaly data to physical properties. The pre-trained network can be used to estimate the distribution of magnetization intensity from new input magnetic data. Two DNN structures were employed to test the feasibility and generalization of the proposed method by implementing Experiments in serval two-dimensional (2D) synthetic examples. Compared with the conventional method, the predicted distribution of magnetization intensity obtained by using our method is more concentrated and has better resolution to determine the boundary of the magnetic body. In a field example of Galinge iron ore deposits in China, the magnetization distribution of concealed iron orebodies inverted by the proposed approach was in good agreement with that detected by borehole data. DNNs have good nonlinear inversion capability and exhibit some excellent merits of strong learning ability, wide coverage and strong adaptability, which is a considerable application prospect in geophysical exploration.

Assessment of depth to magnetic sources using high resolution aeromagnetic data of some parts of the Lower Benue Trough and adjoining areas, Southeast Nigeria

High resolution airborne magnetic data acquired between 2005 and 2010 were used to determine depth to shallow and deep magnetic sources in some parts of Southeastern Nigeria. Various depth estimation methods such as standard Euler deconvolution (SED), source parameter imaging (SPI), spectral depth analysis (SDA) and two dimensional (2-D) forward modeling were applied. Results obtained from SED, SPI and models of profiles 1 and 2 indicate that the Abakaliki Anticlinorium (AA) and Ikom-Mamfe Rift (IMR) regions are dominated by short wavelength magnetic anomalies caused by extensive tectonic events. The SED map showed depth to shallow and deep magnetic sources ranging from ~ 16.6 to ~ 338.3 m and ~ 394.3 to ~ 5748.1 m respectively. Likewise, depth estimates from the SPI map varies from ~ 147.1 to ~ 554.2 m (shallow magnetic sources) and ~ 644.2 to ~ 6141.6 m (deep magnetic sources). The result obtained from SDA revealed depths to deep magnetic basement in the range of ~ 769 to ~ 6666 m with an average of ~ 3449 m. Also, it showed that depth to shallow magnetic sources vary between ~ 119 and ~ 434 m with mean of ~ 269 m. The 2-D forward modelling showed maximum depth values of ~ 4700, ~4600 and ~ 6500 m in the models of profiles 1, 2 and 3 within the Anambra Basin (AB), Afikpo Syncline (AS) and Calabar Flank (CF) respectively. Generally, from all the various methods applied the results indicate that AB, AS and CF are dominated by long wavelength anomalies. The 2-D models indicated that the basement framework is undulant. Also, depth estimates involving the various methods used in this study correlate strongly with each other in the AB, AS and CF geological regions.

Integrating of Remote Sensing and Airborne Magnetic Data to Outline the Geologic Structural Lineaments That Controlled Mineralization Deposits for the Area around Gabal El-Niteishat, Central Eastern Desert, Egypt

The Gabal (G.) El-Niteishat area lies in the Central Eastern Desert of Egypt which is known for various mineral resources and geological structures. Umm Gheig, Umm Naggat, Umm Shaddad, Wadi (W.) Zeidun and Sigdit represent some important regions that contain mineral deposits in the study area. Various filters such as first vertical derivative (FVD), horizontal gradient magnitude (HGM), tilt derivative (TDR) and near-surface were applied to the airborne magnetic data for the study area to deduce the structural lineaments and magnetic source edges which were controlled by the presence of mineral deposits. Processed Landsat ETM+ images are used for delineating the rock unit boundaries that are exposed in the study area such as serpentinite, metagabbro, metavolcanics and metasediments. Also, band ratios, principal component analysis (PCA) and false-color composite image (Crosta alteration image) were applied to get specific results about the alteration zones. The structural lineaments analysis illustrated that the common trends that affected the study area were NW-SE, NE-SW, E-W and N-S. Integration of remote sensing and airborne magnetic data exhibited the relation between mineralization and structural lineaments.

Curie depth and heat flow analysis of Ikot Ekpene and environs, Eastern Niger Delta Basin, using airborne magnetic data

Aeromagnetic dataset over Ikot Ekpene and environs, Eastern Niger Delta Basin, was processed to compute the basement depth, Curie isotherm depth, geothermal gradient and heat flow within the area in order to investigate the depth to magnetic sources, geothermal prospect and the hydrocarbon potential of the place. The adopted computational method transformed the spatial data into frequency domain and provided a relationship between radially average power spectrum of the magnetic anomalies and the depths to respective sources. The results of the analysis showed that the depths to centroids and top boundaries range from 7.84 to 13.38 km and 0.233 to 0.459 km respectively. Curie depths within the basin undulate and vary between 15.42 and 26.49 km. The geothermal gradients range between 20.758 and 35.649 ⁰C/km while the corresponding heat flow is about 51.896 mWm⁻² within east of Ikono, north of Mbak and west of Abak Areas and 89.124 mWm⁻² within Amawum, Ndoro, Isiala, Ogbuebule and east of Uyo Areas. Based on the computed sedimentary thicknesses, high geothermal gradients and delineated major faults and fractures which could serve as migratory pathway for hydrocarbon or hydrothermal fluid, some parts of the study area have been demarcated for geothermal prospect and detail hydrocarbon exploration. Â

Identification of intrusive lithologies in volcanic terrains in British Columbia by machine learning using random forests: The value of using a soft classifier

Identifying the location of intrusions is a key component in exploration for porphyry Cu ± Mo ± Au deposits. In typical porphyry terrains, in the absence of outcrop, intrusions can be difficult to discriminate from the compositionally similar volcanic and volcanoclastic sedimentary rocks in which they are emplaced. The ability to produce lithological maps at an early exploration stage can significantly reduce costs by assisting in planning and prioritization of detailed mapping and sampling. Additionally, a data-driven strategy provides opportunity for the discovery of intrusions not identified during conventional mapping and interpretation. We used random forests (RF), a supervised machine-learning algorithm, to classify rock types throughout the Kliyul porphyry prospect in British Columbia, Canada. Rock types determined at geochemical sampling sites were used as training data. Airborne magnetic and radiometric data, geochemistry, and topographic data were used in classification. Results were validated using First Quantum Minerals’ geologic map, which includes additional detail from targeted location and transect mapping. The petrophysical and compositional similarity of rock types resulted in a noisy classification. Intrusions, particularly the more discrete, were inconsistently predicted, likely due to their limited extent relative to data sampling intervals. Closer examination of class membership probabilities (CMPs) identified locations where the probability of an intrusion being present was elevated significantly above the background. Indeed, a large proportion of mapped intrusions correspond to areas of elevated probability and, importantly, areas were highlighted as potential intrusions that were not identified in geologic mapping. The RF classification produced a reasonable lithological map, if lacking in resolution, but more significantly, great benefit comes from the insights drawn from the RF CMPs. Mapping the spatial distribution of elevated intrusion CMP, a soft classifier approach, produced a map product that can target intrusions and prioritize detailed mapping for mineral exploration.

Integrated Geophysical Investigation for Shallow-Scale Massive (Pb-Zn) Sulphide and Barite Exploration in the Abakaliki and Obubra Mining Districts (AOMD), Southeastern Nigeria

This contribution presents the first multi-scale geophysical investigation conducted in the Abakaliki and Obubra mining districts (AOMD), southeastern Nigeria, for the exploration of shallowly emplaced (< 50 m) massive sulfide and barite mineral deposits. Raw airborne magnetic data covering the AOMD was acquired and subjected to diverse enhancement routines to facilitate the identification of areas of subtle magnetic anomalies for subsequent ground truthing. The mineralization is mainly made up of galena, sphalerite, and barite. The contrasting physical properties between these minerals and their host rocks make resistivity, induced polarization (IP), spontaneous potential (SP), and magnetic methods ideal to study and delineate the mineralization and possibly locate the ore bodies. Wenner DC resistivity and IP data acquired along the same parallel traverses are inverted using Res2dinv to produce 2D tomography models of resistivity and chargeability. The resistivity and IP models were successful in delineating zones of sulfide and barite mineralization. The ground SP data collected along the same traverses as the resistivity and IP data were interpreted as line graphs. The 1D SP models show a number of central negative anomalies that could be related to sulfide mineralization. These anomalies show fair-to-good spatial continuity and are in close agreement with the result obtained in the resistivity and chargeability model. This study suggests the presence of exploitable mineral deposits in the AOMD and demonstrates the superior interpretative ability of combined geophysical methods over a single geophysical approach.

Interpolation artifacts as a result of spatial aliasing: A case study of the airborne magnetic data set of southeastern Minas Gerais, Brazil

Linear features at an acute angle with the flight direction are imaged as a series of aligned circular anomalies in the images of Area 15 aeromagnetic survey, which covered part of the Brazilian southeastern region. These features are interpolation artifacts, a recurring problem found in airborne magnetic images that cause problems for qualitative and quantitative geophysical-geologic interpretation. This imaging problem is attributed to spatial aliasing. By running simulations of magnetic data on a synthetic model, we have physically demonstrated that the interpolation artifacts from Area 15 are due to inappropriate survey design. Besides the most common expression of artifacts, we described a geologically noncoherent linear pattern as a new type of artifact. Supported by spectral analyses, we found that the Area 15 aliased spectrum is similar to geologic high-frequency magnetic features, which constitutes a motive for unearthing the correct geophysical signal. Thus, we made use of four techniques for removing the artifacts. The trend enforcement method partially improved the images, whereas the inverse interpolation method was ineffective, apparently because Area 15 data are severely aliased. The constrained coherence diffusion and multitrend gridding methods were able to significantly reduce the presence of artifacts. Despite the high-frequency attenuation, these tools adequately enhanced the magnetic trends and minimized the artifacts. Therefore, the improved images are better suited for reliable geologic interpretation.

Joint inversion of potential-fields data over the DO-27 kimberlite pipe using a Gaussian mixture model prior

We have carried out petrophysically and geologically guided inversions (PGIs) to jointly invert airborne and ground-based gravity data and airborne magnetic data to recover a quasi-geology model of the DO-27 kimberlite pipe in the Tli Kwi Cho (also referred to as TKC) cluster. DO-27 is composed of three main kimberlite rock types in contact with each other and embedded in a granitic host rock covered by a thin layer of glacial till. The pyroclastic kimberlite (PK), which is diamondiferous, and the volcanoclastic kimberlite (VK) have anomalously low density, due to their high porosity, and weak magnetic susceptibility. They are indistinguishable from each other based upon their potential-field responses. The hypabyssal kimberlite (HK), which is not diamondiferous, has been identified as highly magnetic and remanent. Quantitative petrophysical signatures for each rock unit are obtained from sample measurements, such as the increasing density of the PK/VK unit with depth and the remanent magnetization of the HK unit, and are represented as a Gaussian mixture model (GMM). This GMM guides the PGI toward generating a 3D quasi-geology model with physical properties that satisfies the geophysical data sets and the petrophysical signatures. Density and magnetization models recovered individually yield volumes that have physical property combinations that do not conform to any known petrophysical characteristics of the rocks in the area. A multiphysics PGI addresses this problem by using the GMM as a coupling term, but it puts a volume of the PK/VK unit at a location that is incompatible with geologic information from drillholes. To conform to that geologic knowledge, a fourth unit is introduced, PK-minor, which is petrophysically and geographically distinct from the main PK/VK unit. This inversion produces a quasi-geology model that presents good structural locations of the diamondiferous PK unit and can be used to provide a resource estimate or decide the locations of future drillholes.

Inversion of magnetic and frequency-domain electromagnetic data for investigating lithologies associated with gold mineralization in the Canadian Malartic area, Québec, Canada

Canadian Malartic is an Archean low-grade bulk tonnage native gold deposit. The deposit is mostly located in altered clastic metasedimentary rocks, mafic–ultramafic dykes, and monzodioritic porphyry intrusions. Airborne magnetic and frequency-domain electromagnetic (EM) data were inverted to reconstruct the geological units associated with the mineralization, especially the intrusive masses. The 3-D inversion of magnetic data, which used a tetrahedral mesh to a depth of 2.4 km, shows that mafic volcanic rocks and iron formation rocks extend to depth in the area, more so than diabase dykes. The magnetic inversion also shows that the diorite and monzodiorite rocks of the Lac Fournière A pluton are dipping toward the south on its northern edge at the contact with the metasedimentary rocks. The 1-D inversion of the frequency-domain EM data, for both electrical conductivity and magnetic susceptibility, is able to reconstruct geological structures to a depth of approximately 100 m, providing more details and information about these features. The intrusive masses such as diabase dykes, diorite and monzodiorite rocks, and mafic volcanic rocks are reconstructed as electrically conductive structures in the inversion results. The metasedimentary rocks are resistive, and the overburden is conductive in most of the area. The geophysical data and inversion results suggest the presence of some features (such as diabase dykes and monzodiorite rocks) that are not yet present on some parts of the geology map. A comparison of the EM-derived susceptibility and the magnetic-derived susceptibility over the iron formations can reveal the effect of remanent magnetization.