Gold Exploration in Two and Three Dimensions: Improved and Correlative Insights from Microscopy and X-Ray Computed Tomography

: In this contribution, we highlight a correlative approach in which three-dimensional structural / positional data are combined with two dimensional chemical and mineralogical data to understand a complex orogenic gold mineralization system; we use the Kirk Range (southern Malawi) as a case study. Three dimensional structures and semi-quantitative mineral distributions were evaluated using X-ray Computed Tomography (XCT) and this was augmented with textural, mineralogical and chemical imaging using Scanning Electron Microscopy (SEM) and optical microscopy as well as ﬁre assay. Our results detail the utility of the correlative approach both for quantifying gold concentrations in core samples (which is often nuggety and may thus be misrepresented by quarter- or half-core assays), and for understanding the spatial distribution of gold and associated structures and microstructures in 3D space. This approach overlays complementary datasets from 2D and 3D analytical protocols, thereby allowing a better and more comprehensive understanding on the distribution and structures controlling gold mineralization. Combining 3D XCT analyses with conventional 2D microscopies derive the full value out of a given exploration drilling program and it provides an excellent tool for understanding gold mineralization. Understanding the spatial distribution of gold and associated structures and microstructures in 3D space holds vast potential for exploration practitioners, especially if the correlative approach can be automated and if the resultant spatially-constrained microstructural information can be fed directly into commercially available geological modelling software. The extra layers of information provided by using correlative 2D and 3D microscopies o ﬀ er an exciting new tool to enhance and optimize mineral exploration workﬂows, given that modern exploration e ﬀ orts are targeting increasingly complex and low-grade ore deposits.


Introduction
Gold exploration continues to attract a vast proportion of the global minerals exploration budget (e.g.,~50% of the $10.1 billion exploration budget for 2018 [1]). Much of this gold has been won from 'orogenic' classes of gold mineralization [2,3], which generally form in metamorphic terranes associated with accretionary or collisional orogenies [4]. These geological settings are typically characterized by a high level of structural complexity, in which gold mineralization may be multi-stage or may dimensional structures and semi-quantitative mineral distributions were evaluated using XCT and this was augmented with textural, mineralogical and chemical imaging using Scanning Electron Microscopy (SEM) and optical microscopy. Gold grades within the core samples were quantified using standard fire assay protocols [26].

High Resolution X-Ray Computed Tomography (XCT)
In order to understand the distribution of gold mineralization in mineralized drill core samples, four segments of the core sample (45 mm diameter; lengths respectively: 15, 3, 5 and 7 cm) were scanned using the Electric Phoenix VTomeX L240 micro X-ray Computed Tomography scanner system at the Central Analytical Facility (CAF), Stellenbosch University (South Africa), manufactured by General Electric Sensing and Inspection Technologies/Phoenix X-ray (Wunstorff, Germany). The XCT facility is described in more detail in [27]. To optimize scanning of gold-sulphide ores, the scanning parameters were set to 100 kV and 100 µA and the X-ray beam was filtered using a 0.5 mm Cu filter. Images were acquired using a dwell time of 500 ms. Each scan took about 4 h and the resulting model had a voxel size of 35 µm. Furthermore, a small subsample (dimensions 5 mm in diameter) cut from the original core sample JJ01 was also scanned at a high resolution using phoenix Nanotom S nanoCT scanner, manufactured by Phoenix X-ray (Wunstorff, Germany), the scanning parameters were set to 150 kV and 200 µA and the X-ray beam was filtered using a 0.5 mm Cu filter. Images were acquired using a dwell time of 500 ms. The resulting model had a voxel size of 10 µm.
These models were reconstructed using a filtered back projection algorithm built into the GE Datos 2.1 reconstruction package [27]. Subsequent data evaluation and visualization was conducted using Volume Graphics VGStudio Max 3.3, Heidelberg, Germany [28].
2.1.1. The XCT Thresholding Procedure for Gold in the Core Samples.
Thresholding is a process of material segmentation in X-ray computed tomography [16], in which materials or features of interest are mathematically separated from a background by distinguishing XCT model voxels according to threshold grey-scale values. To distinguish between voxels that represent gold, voxels that represent sulphides and voxels that represent matrix silicate minerals, we employed the global segmentation technique reported in [24]. Briefly, the global segmentation operation evaluates the grey-scale histogram of the entire core sample, identifies the main peaks in the histogram and defines the vertical segmentation line at the grey-scale value that marks the dip between the adjacent peaks. However, this method has a limitation of point spread function (PSF) [29]. Therefore, this threshold value is subsequently refined using a local threshold gradient search, in which the VGStudio Max 3.3 software identifies the sharpest change in grey value around the initial region of interest [24]. Figure 1 documents the grey value histograms for samples JJ01, JJ02, JJ04 and JJ05, wherein the automated global and local segmentation thresholding techniques have defined three main peaks. Based on the mineral densities and model mineral X-ray attenuation coefficients for the minerals that we have identified in our core samples (Table 1), we assume that the left peak denotes the less dense bulk silicate mineral matrix, the middle peak represents the sulphides, and the right peak denotes the much denser gold. Gold grains were defined using the minimum grey value threshold of 59,868. The peak of gold shifts a bit among the samples, but we still use the same threshold value for all the samples because the variation is not much. Sulphides were defined using the minimum threshold of CT number 49,774. These histogram-derived threshold values represent mathematical transformations of the original X-ray data, which were collected across 16 bit grey values. Because of the closeness between X-ray attenuation coefficients between the different sulphide minerals (Figure 1e), the XCT grey value alone could not distinguish between the different sulphide phases. The validity of these peak assignments was confirmed by comparing the grey value distributions in XCT image slices to correlated polished surfaces that were evaluated using optical microscopy and scanning electron microscopy (SEM). Notably, the SEM images revealed that many of the gold grains were much smaller than the 35 µm XCT scanning resolution. This partial volume effect is one reason for the broadness of the gold peaks (Figure 1), and although these partial volumes still allow identification of where gold is preferentially sited (by raising the voxel grey value), it does mean that the XCT-based approach should be acknowledged as being semi-quantitative. Using the thresholded data, we could visualize mineral grain size distribution and textures in 3D space. Furthermore, the thresholded data enabled us to generate volumes for each of the mineral phases (which were then used to calculate masses by considering known densities (Table 1)).
Minerals 2019, 9, x FOR PEER REVIEW 4 of 20 data, we could visualize mineral grain size distribution and textures in 3D space. Furthermore, the thresholded data enabled us to generate volumes for each of the mineral phases (which were then used to calculate masses by considering known densities (Table 1)). Histograms of the core samples from grey value mapping used in thresholding. The histograms distinguish three dominant peaks. The left peak denotes the less dense bulk silicate mineral matrix, the middle peak the sulphides and the right peak denotes the much denser gold mineral. Samples JJ01, JJ02, JJ04 and JJ05, x axis denotes grey values and y axis frequency. (e): X-ray attenuation coefficient at 100 keV [30]. The graph indicates that at 100 KeV energy, gold will have the highest attenuation by a large margin. The XCT technique could not resolve different sulphides from one another as their density range is narrow, but clearly distinguishes gold as high atomic number phases distinct from the sulphides.  The histograms distinguish three dominant peaks. The left peak denotes the less dense bulk silicate mineral matrix, the middle peak the sulphides and the right peak denotes the much denser gold mineral. Samples JJ01, JJ02, JJ04 and JJ05, x axis denotes grey values and y axis frequency. (e): X-ray attenuation coefficient at 100 keV [30]. The graph indicates that at 100 KeV energy, gold will have the highest attenuation by a large margin. The XCT technique could not resolve different sulphides from one another as their density range is narrow, but clearly distinguishes gold as high atomic number phases distinct from the sulphides. Table 1. Density range and X-ray attenuation coefficients of minerals present in the core samples. Density data are collected from [31] and X-ray attenuation coefficients are calculated using an Excel sheet developed by [32]. To empirically verify that the small peak is indeed gold and is not confused with other minerals, gold grains identified in XCT (brightest and most dense spots) were intersected and proved in reflected light microscopy. We cut the core at an appropriate location and made thick sections. Reflected light microscopy identifies ubiquitous gold grains within the core samples, which proves that the dense minerals identified in XCT are indeed gold ( Figure 2). Furthermore, in Figure 2c, the gold grain identified in XCT is proved in the reflected light microscopy in Figure 2d. This empirically proves that the dense minerals found in our samples are gold. In addition, we extracted grey values representing 100% gold, grey values representing 100% sulphides and grey value representing 100% silicates. Extracted XCT grey values for gold brightest spot grains range in 61,586-65,000, falling within the small peak on the histogram; those for sulphides range in 50,000-59,000, falling within the sulphide peak; and those for silicates range in 1600 and 49,771, falling within the silicate peak. This shows that our segmentation method is correct. However, grey values obtained from different scanners may differ, which would make comparison of the data difficult. To standardize this, grey values need to be correlated with Hounsfield units. This can be achieved by plotting the linear attenuation coefficient of the minerals in the core samples against grey values to provide a linear fit and a regression equation. The linear attenuation coefficient for each mineral at 100 keV energy can be converted to CT numbers' Hounsfield units (HU) using the equation HU = (µ material − µ water)/(µ water) × 1000, where µ material = linear attenuation coefficient of material and µ water = linear attenuation coefficient of water. This can make the data comparable and has been demonstrated by the authors in [33]. Minerals 2019, 9, x FOR PEER REVIEW 6 of 20 Figure 2. X-ray computed tomography and reflected light images indicating gold grains identified in XCT and gold grains identified in reflected light microscopy. (a) The drill core was cut on the rectangle shown to make thick sections. Grey colour in the drill core is the silicate matrix, the yellow colour indicates sulphides and other colours represent gold, distributed according to gold grain volume (warm colours = large grains, cold colours = smaller grains). (b) Reflected light microscopy shows ubiquitous presence of gold, which proves that the dense minerals we see in XCT are indeed gold. (c,d) Gold grain in XCT, same grain in reflected light microscopy. In Figure 2c, the yellow colour indicates sulphides, dark colour indicates silicate matrix and red colour indicates gold. Figure 2c is coloured according to density (warm colours = most dense gold grains, yellow colours = medium density sulphides, dark/grey colour = less dense silicate matrix).

Thin and Polished Section Photo Montages
Based on the gold distribution information provided by XCT data, we physically cut the cores samples at appropriate positions to intersect the identified gold and made thin and thick sections. The appropriate surfaces were polished and subjected to reflected and transmitted light microscopy. Thin and thick sections were described using a Nikon Eclipse E 200 petrological microscope at Stellenbosch University; digital photomicrographs were captured and stitched to make image mosaic of the thin and thick sections using PT GUI software. The 2D optical microscopy image was compared with the 3D XCT image for a specific region of the samples. The XCT data and microscopy data were aligned manually by visually aligning prominent features using commercially available Volume Graphics software [28]. Correlative microscopies in other geoscience studies (e.g., fluorescence microscopy + SEM [34]) have proven to be more precise and can expose the physical and chemical character of ore minerals at a better spatial resolution and analytical precision [35].  The drill core was cut on the rectangle shown to make thick sections. Grey colour in the drill core is the silicate matrix, the yellow colour indicates sulphides and other colours represent gold, distributed according to gold grain volume (warm colours = large grains, cold colours = smaller grains). (b) Reflected light microscopy shows ubiquitous presence of gold, which proves that the dense minerals we see in XCT are indeed gold. (c,d) Gold grain in XCT, same grain in reflected light microscopy. In Figure 2c, the yellow colour indicates sulphides, dark colour indicates silicate matrix and red colour indicates gold. Figure 2c is coloured according to density (warm colours = most dense gold grains, yellow colours = medium density sulphides, dark/grey colour = less dense silicate matrix).

Thin and Polished Section Photo Montages
Based on the gold distribution information provided by XCT data, we physically cut the cores samples at appropriate positions to intersect the identified gold and made thin and thick sections. The appropriate surfaces were polished and subjected to reflected and transmitted light microscopy. Thin and thick sections were described using a Nikon Eclipse E 200 petrological microscope at Stellenbosch University; digital photomicrographs were captured and stitched to make image mosaic of the thin and thick sections using PT GUI software. The 2D optical microscopy image was compared with the 3D XCT image for a specific region of the samples. The XCT data and microscopy data were aligned manually by visually aligning prominent features using commercially available Volume Graphics software [28]. Correlative microscopies in other geoscience studies (e.g., fluorescence microscopy + SEM [34]) have proven to be more precise and can expose the physical and chemical character of ore minerals at a better spatial resolution and analytical precision [35].

SEM Mapping
Scanning Electron Microscopy combined with an Electron Dispersive Spectrometer (SEM-EDS) was utilized to provide further insights into the mineralogy, textures and distribution of gold in the evaluated samples. Polished blocks, previously evaluated using reflected light optical microscopy, were coated with carbon to prevent localized charging and any resultant distortions or reflections of the electron beam. The SEM/EDS analyses were conducted at the Electron Microscopy Unit of the Central Analytical Facility (Stellenbosch University, South Africa) using a Zeiss EVO MA15VP Scanning Electron Microscope. SEM/EDS analyses were carried out in a nitrogen atmosphere at pressures from 65 to 96 Pa, voltages from 20 to 30 kV, working distance for EDS analyses of 8.5 mm and magnifications ranging between 1000× and 5000×.

Fire Assay
Fire assay analyses were conducted on subsamples of each of the cores that were scanned using XCT. Each subsample sent for fire assay was sectioned both physically and in the corresponding XCT reconstructed model. This ensured that the assayed gold grades could be accurately compared to calculated gold grades, which were derived from the XCT gold volumes generated after the thresholding procedure described in Section 2.1.1. The remaining core sample was reserved for future analytical work. Fire assay analyses were conducted at the ALS mineral services laboratory (Johannesburg, South Africa). Gold was assayed using fusion process at temperatures of 1200 • C to obtain the lead button that was then subjected to cupellation. Full details for this analytical technique are reported in reference [26].

Thresholding and Quantification using XCT Data
The histogram of the measured grey values ( Figure 1) provides a clear indication that the core sample contains mineral material from three markedly different density (or more correctly; X-ray attenuation coefficient) fractions viz. the silicate matrix, the sulphide paragenesis, and a distinctly high density material that we assume to largely represent gold. Other mineral phases that could be present with such high X-ray attenuation coefficients include galena, barite and native silver; however, none of these minerals were identified in the follow-up petrographic investigations.
The silicate mineralogy is quantitatively dominant throughout the core samples evaluated and have grey values ranging between 1600 and 49,771. Detailed interrogation of the XCT data revealed that we can see foliation within the silicates ( Figure 3a) and pick out some garnet porphyroblasts ( Figure 3b). Garnet porphyroblasts (density 3.93 gcm 3 ) have relatively high grey values in the range 42,000 to 47,000, which are close to the grey value range expected for sulphides (Table 1). Thin sections made from the core prove the presence of garnet ( Figure 3c). Due to the closeness of garnet grey values and those of sulphides, there could be a possible overlap of garnet with sulphides' peak. Immediately adjacent to the bright sulphides vein, there is a change in grey values, which may be an indicator of hydrothermal alteration (Figure 3b), followed by a dark silicate matrix with grey values less than 20,000. From the silicate matrix, we can see three phases, garnets, hydrothermally alteration adjacent to the vein and the dark silicate minerals that cannot be differentiated due to close density ranges. From microscopy, these dark silicate minerals are quartz, feldspars, amphiboles and biotite ( Table 1).
The sulphide mineral paragenesis is represented in the XCT data by grey values ranging between 49,774 and 59,868. From reflected light microscopy, this assemblage includes minor amounts of chalcopyrite, sphalerite and pyrrhotite, with pyrite as the dominant (>80% modal proportion among sulphide minerals) sulphide phase. The density for these sulphides range between 4.1 and 5.01 g/cm 3 (Table 1), and it is thus, not possible to discern between the different sulphide minerals using XCT alone (highlighting the need for correlative microscopies in order to derive additional detailed mineralogical information (Section 4.2)). In the core samples evaluated, the grey values representing the sulphides are generally distributed along planar features which reflect the main vein network, and which in turn is structurally controlled in response to the regional stresses experienced by the geology during its deformational history. By considering the number of voxels ascribed to sulphide minerals through our thresholding approach and by deriving a weighted average density for the sulphide paragenesis, the mass of sulphide could be estimated for each of the core samples. These estimated masses ranged between 1.15 and 5.3 g ( Table 2).
for the sulphide paragenesis, the mass of sulphide could be estimated for each of the core samples. These estimated masses ranged between 1.15 and 5.3 g ( Table 2).
Gold in the core samples was identified by using a minimum grey value threshold value of 59,868. This represents the brightest pixels in the XCT images (Figure 4), since gold has the highest X-ray attenuation coefficient of all the minerals identified in the cores (Figure 1e). The presence of gold in our samples was confirmed using reflected light microscopy, Scanning Electron Microscopy (Figures 2 and 3d), and traditional fire assay. Voxels representing gold were predominantly found in association with the sulphides in the vein network, although a small proportion of these voxels were also located in the silicate host rock, particularly in close proximity to the vein network (Figure 4a, b). The identified gold grains ranged in size between less than 0.14 to 0.7 mm in diameter, and a calculation based on the volume of the gold voxels and the known density of gold (19.32 g.cm −3 ) estimates that the mass of gold in each sample ranges between 0.095 and 0.25 g, corresponding to a calculated gold grade of between 0.03 and 0.14 ppm (Table 2).  Gold in the core samples was identified by using a minimum grey value threshold value of 59,868. This represents the brightest pixels in the XCT images (Figure 4), since gold has the highest X-ray attenuation coefficient of all the minerals identified in the cores (Figure 1e). The presence of gold in our samples was confirmed using reflected light microscopy, Scanning Electron Microscopy (Figures 2 and 3d), and traditional fire assay. Voxels representing gold were predominantly found in association with the sulphides in the vein network, although a small proportion of these voxels were also located in the silicate host rock, particularly in close proximity to the vein network (Figure 4a,b). The identified gold grains ranged in size between less than 0.14 to 0.7 mm in diameter, and a calculation based on the volume of the gold voxels and the known density of gold (19.32 g.cm −3 ) estimates that the mass of gold in each sample ranges between 0.095 and 0.25 g, corresponding to a calculated gold grade of between 0.03 and 0.14 ppm (Table 2). Table 2. Distribution of gold and sulphides in analysed drill core samples. Average density of sulphides was calculated by averaging the densities of individual sulphides identified in 2D microscopy.

Sample
Weight of Core (g)

Fire Assay to Calibrate XCT Threshold Values
Fire assay analyses were conducted to confirm the presence of gold in our samples and to calibrate the calculated gold grades derived from XCT data. The fire assay analyses were conducted on quarter cores derived from the original scanned cores, and these data compared to calculated gold grades for the corresponding quarter of the core, which was sectioned digitally from the 3D XCT model ( Figure 5). Fire assay data reveal that the gold grades for the respective cores in ppm are 0.03 (JJ1), 0.02 (JJ2), 0.01(JJ4) and 0 (JJ5), ( Table 3). The fire assay shows a strong and positive linear relationship with the gold grades calculated from our XCT thresholding protocol (Figure 6a: y = 0.5888x − 0.0011; r 2 = 0.5). The linear best-fit line has a slope that is less than the 1:1 line, indicative that the XCT protocol generally overestimates the mass of gold present in the core samples. Fire assay analyses were conducted to confirm the presence of gold in our samples and to calibrate the calculated gold grades derived from XCT data. The fire assay analyses were conducted on quarter cores derived from the original scanned cores, and these data compared to calculated gold grades for the corresponding quarter of the core, which was sectioned digitally from the 3D XCT model ( Figure 5). Fire assay data reveal that the gold grades for the respective cores in ppm are 0.03 (JJ1), 0.02 (JJ2), 0.01(JJ4) and 0 (JJ5), ( Table 3). The fire assay shows a strong and positive linear relationship with the gold grades calculated from our XCT thresholding protocol (Figure 6a: y = 0.5888x − 0.0011; r 2 = 0.5). The linear best-fit line has a slope that is less than the 1:1 line, indicative that the XCT protocol generally overestimates the mass of gold present in the core samples. Figure 5. (a,b) Quarter cores that were sent for fire assay analysis cut from the whole cores. The grey colour in the drill core is the silicate matrix and other colours represent gold, distributed according to gold grain volume (warm colours = large grains, cold colours = smaller grains). Samples JJ01 and JJ02.
The equation of the best-fit line does however provide a useful internal calibration for the XCT data. Figure 7 compares the distribution of gold from the XCT model data using the original threshold value derived from the global segmentation protocol 59,868, to a threshold value that has been selected to ensure that the XCT calculated gold grades fall on the 1:1 line with fire assay gold grades. This exercise reveals that calibrating XCT values based on fire assay to the threshold grey value in our case, of 63,000, other dense minerals disappear, and gold is only present in the vein and not necessarily in the wall rock ( Figure 7b) and is ~1:1 with fire assay (Figure 6c, Table 3), disregarding the JJ05 sample in which the assayed gold grades were below the detection limit. If we take grey values below 63,000, we have more gold grains making XCT to generally overestimate the mass of gold present in the core sample (Figure 7a).
A final test was conducted on the calibrated XCT core models in order to comment on the in situ nugget effect that may be induced by only assaying a quarter core. Quarter cores have relatively high gold values than the whole core due to the nugget effect (Table 3), except sample JJ04. The nugget effect is attributed to the heterogeneous distribution of gold grains in the core samples. Some regions of the core are rich in gold and others are not. This is so because vein systems typically contain both Small grain large grain Figure 5. (a,b) Quarter cores that were sent for fire assay analysis cut from the whole cores. The grey colour in the drill core is the silicate matrix and other colours represent gold, distributed according to gold grain volume (warm colours = large grains, cold colours = smaller grains). Samples JJ01 and JJ02.
The equation of the best-fit line does however provide a useful internal calibration for the XCT data. Figure 7 compares the distribution of gold from the XCT model data using the original threshold value derived from the global segmentation protocol 59,868, to a threshold value that has been selected to ensure that the XCT calculated gold grades fall on the 1:1 line with fire assay gold grades. This exercise reveals that calibrating XCT values based on fire assay to the threshold grey value in our case, of 63,000, other dense minerals disappear, and gold is only present in the vein and not necessarily in the wall rock ( Figure 7b) and is~1:1 with fire assay (Figure 6c, Table 3), disregarding the JJ05 sample in which the assayed gold grades were below the detection limit. If we take grey values below 63,000, we have more gold grains making XCT to generally overestimate the mass of gold present in the core sample (Figure 7a).
A final test was conducted on the calibrated XCT core models in order to comment on the in situ nugget effect that may be induced by only assaying a quarter core. Quarter cores have relatively high gold values than the whole core due to the nugget effect (Table 3), except sample JJ04. The nugget effect is attributed to the heterogeneous distribution of gold grains in the core samples. Some regions of the core are rich in gold and others are not. This is so because vein systems typically contain both barren and mineralized segments with gold grade varying laterally, vertically and across the core. Furthermore, finer grained gold mineralization usually has a high geological nugget effect [36]. Hence, some of the quarter cores have high gold values compared to the whole core (Table 3). barren and mineralized segments with gold grade varying laterally, vertically and across the core. Furthermore, finer grained gold mineralization usually has a high geological nugget effect [36]. Hence, some of the quarter cores have high gold values compared to the whole core (Table 3).   Figure 6. Correlation between XCT and Fire assay data. (a) Correlation between XCT and fire assay of quarter core before calibration. (b) 1:1 linear best fit line of XCT and fire assay. (c) Correlation between XCT and fire assay of quarter core after calibration. Error bar ± 0.001.

Textural and Mineralogical Correlations between 2D and 3D Study Approaches
The preceding section highlights the utility of XCT towards understanding and semiquantitatively analysing gold distributions in three dimensional reconstructions of exploration borehole core samples. When used independently however, this approach fails to discern between different sulphide minerals in the assemblage and poorly resolves microstructural features that may have important implications for the exploration efforts. It is thus imperative to combine 3D XCT analyses with conventional 2D microscopies in order to derive the full value out of a given exploration drilling program. Figure 8 compares a montage of 2D optical (reflected light) microscopy images with the corresponding image slice extracted from an XCT model. The correlation was conducted manually by visually aligning prominent features using the commercially available Volume Graphic software. An artefact of this imperfect overlap is depicted in Figure 8a, in which a large grain Small grain Figure 7. (a,b): Distribution of gold grains in the core after segmenting using a threshold value of 59,868 and 63,000 grey values respectively, at 59,868 XCT generally overestimate the mass of gold present in the core sample. The grey colour is the silicate matrix and other colours represent gold, distributed according to gold grain volume (warm colours = large grains, cold colours = smaller grains). Sample JJ01.

Textural and Mineralogical Correlations between 2D and 3D Study Approaches
The preceding section highlights the utility of XCT towards understanding and semi-quantitatively analysing gold distributions in three dimensional reconstructions of exploration borehole core samples. When used independently however, this approach fails to discern between different sulphide minerals in the assemblage and poorly resolves microstructural features that may have important implications for the exploration efforts. It is thus imperative to combine 3D XCT analyses with conventional 2D microscopies in order to derive the full value out of a given exploration drilling program. Figure 8 compares a montage of 2D optical (reflected light) microscopy images with the corresponding image slice extracted from an XCT model. The correlation was conducted manually by visually aligning prominent features using the commercially available Volume Graphic software. An artefact of this imperfect overlap is depicted in Figure 8a, in which a notable gold grain is intersected in the image slice but was either slightly below the polished surface or was above the polished surface and was subsequently cut and polished away. Indeed, when preparing polished sections of finely disseminated gold mineralization, intersection of discrete gold grains can only be fortuitous if cutting and polishing is not guided by prior XCT analysis. Figure 8b further contrasts the level of mineralogical detail that is achieved in the XCT image slice versus the petrographic montage. The XCT model provides a three-dimensional overview of the sulphide spatial distribution, reveals low-level textural information (e.g., mineral stretching along the shear plane is clearly evident), and enables a truer representation of the mineral size distributions (2D approaches provide apparent mineral sizes). However, it clearly lacks in its ability to discern between different sulphide mineral phases. In comparison, the 2D reflected light petrographic montage enables that the individual sulphide mineral phases can be identified, their paragenetic relationships can be more deeply investigated, and additional microstructural details (e.g., mica fish, rotated porphyroblasts) can be used to augment and enrich the understanding of the structural controls on the mineralization. Some disadvantages of the 2D petrographic study approaches include that they are destructive to the original core sample, the samples require preparation (i.e., cutting and polishing may be time-consuming), selection of sub-samples for thin/thick sectioning may introduce unconscious bias into the ultimate dataset, true mineral sizes cannot be discerned, and the 3D kinematics associated with the structural controls cannot be fully studied. These results (summarized in Table 4) clearly highlight that there is extra information that can be obtained from microscopy that cannot be obtained from XCT, and vice versa. Hence, from an optimized exploration perspective, it is important that these techniques should be used correlatively. Table 4. Summary on information that can be obtained from combined approach.

2D Petrography 3D Computed Tomography Combined and 2D and 3D Correlative Approach
Apparent size True size and orientation True size The correlation was conducted manually by visually aligning prominent features using Volume Graphic software. A notable gold grain is intersected in the XCT image slice but was either slightly below the polished surface or was above the polished surface and was subsequently cut and polished away. The highresolution X-ray computed tomography and microscopy analyses indicate that gold and accessory sulphides are not randomly distributed but rather follow specific orientation. They are distributed along shear fracture planes rather than disseminated throughout the cores. This suggests that gold was precipitated along the shear plane together with the sulphides during the ore-forming processes. The bright gold grain is not seen in reflected light microscopy, but it is easily seen in XCT image. (b) Zoomed in image of sulphides from XCT indicating that they cannot be differentiated due to a narrow density range. (c) The sulphide mineral assemblage comprises pyrite, pyrrhotite and chalcopyrite determined from detailed 2D examination of polished sections. Thick section cut from sample JJ04.

Discussion
The use of XCT in combination with microscopy provides an excellent tool for understanding gold mineralization, which is often nuggety and difficult to study because of its ppm-level grades. Although the application of this combined approach is not yet widespread, this approach overlays complementary datasets from 2D and 3D analytical protocols, thereby allowing a better and more comprehensive understanding on the distribution and structures controlling gold mineralization. In order to investigate the feasibility of the combined approach and its practical utility, we consider an example from the Manondo-Choma gold prospect (southern Malawi). The correlation was conducted manually by visually aligning prominent features using Volume Graphic software. A notable gold grain is intersected in the XCT image slice but was either slightly below the polished surface or was above the polished surface and was subsequently cut and polished away. The high-resolution X-ray computed tomography and microscopy analyses indicate that gold and accessory sulphides are not randomly distributed but rather follow specific orientation. They are distributed along shear fracture planes rather than disseminated throughout the cores. This suggests that gold was precipitated along the shear plane together with the sulphides during the ore-forming processes. The bright gold grain is not seen in reflected light microscopy, but it is easily seen in XCT image. (b) Zoomed in image of sulphides from XCT indicating that they cannot be differentiated due to a narrow density range. (c) The sulphide mineral assemblage comprises pyrite, pyrrhotite and chalcopyrite determined from detailed 2D examination of polished sections. Thick section cut from sample JJ04.

Discussion
The use of XCT in combination with microscopy provides an excellent tool for understanding gold mineralization, which is often nuggety and difficult to study because of its ppm-level grades. Although the application of this combined approach is not yet widespread, this approach overlays complementary datasets from 2D and 3D analytical protocols, thereby allowing a better and more comprehensive understanding on the distribution and structures controlling gold mineralization. In order to investigate the feasibility of the combined approach and its practical utility, we consider an example from the Manondo-Choma gold prospect (southern Malawi).

Geological Context
The Kirk Range is a medium-to high-metamorphic grade mobile belt which is related to the Southern Irumide orogeny and which is located in south-central Malawi [37]. It is situated approximately  [38]. These lithological units are characterized by polyphase deformation showing evidence of reworking possibly by the pan-African orogeny (800-550 Ma) [39].
Gold mineralization at the Manondo-Choma prospect has recently been described by the authors in [38]. Briefly, the gold is hosted predominantly in biotite schists and gneisses and is associated with a sulphide mineral paragenesis. Regional structures as well as mineralized quartz-sulphide vein sets exhibit a NE-SW trend, which is consistent with the deformation fabrics associated with the southern Irumide orogeny [40][41][42]. Textural and microstructural observations reveal that the vein sets have experienced subsequent deformation and overprinting, which has likely led to transposition of the primary gold mineralization on a regional scale [38].

Correlated 2D and 3D Insights into Gold Mineralization at the Manondo-Choma Prospect
Four segments of drill core were selected for detailed scanning and modelling using XCT and correlative reflected and transmitted light microscopies. The XCT reconstructions reveal that the gold is hosted primarily in sub-vertical vein structures that strike in a NE-SW direction (Figure 9b). The fine-grained gold flakes are sub-rounded in texture, range in size between 0.01 and 0.04 mm, and are typically elongated parallel to the direction of the strike of the host vein ( Figure 9c). Gold grades for the different core segments range between 0.01 and 0.03 ppm from fire assay and 0.01-0.031 ppm calculated from XCT (Table 3; after corrections to account for nugget effects inherent in gold mineralization (see Section 3.1.1.)) and show good agreement with the gold grades measured using fire assay analytical protocols. The gold is strongly associated with sulphide minerals which also orientate preferentially in the NE-SW strike direction. Minor gold is also associated with the wall rock, perhaps reflecting a second-order lithological or geochemical control on the siting of the mineralization [38]. Although the XCT image slices show subtle differences in the grey values among different sulphide minerals, these grey values could not be used to unequivocally determine sulphide speciation and did not allow detailed comment on the textural and temporal relationships within the ore paragenesis.
Thus, to augment the gold spatial, size and shape distribution data, correlative 2D petrographic analyses were conducted on selected sub-samples of the core segments. These sub-samples were best selected based on prior evaluation of the 3D XCT models, whereby gold-rich regions of the core segment were targeted for the further detailed petrographic work. Reflected light microscopy revealed that the sulphide mineral assemblage is dominated by pyrite (~80% modal proportion), with the remaining 20% of the sulphides represented by chalcopyrite, sphalerite and pyrrhotite. Early formed pyrrhotite showed the strongest mineral stretching, which was parallel to both the stretching identified in the gold flakes and to the orientation of the main vein structure. Transmitted light microscopy enabled that microtextures and shear sense indicators could be identified within the silicate mineralogy associated with the vein structure. For example, 2D microscopy indicates the presence of shear in the silicate mineralogy within the vein structure defined by grain size reduction (Figure 9e) and display flow structures that indicate ductile deformation of the core (Figure 9f,g). In standard 2D petrography, this information can give insights into whether the mineralized vein experienced a dextral or sinistral sense of shear (i.e., X-Y directionality). However, by considering this information in light of the known orientation of the vein from prior XCT work, the kinematics can further constrained with reference to the Z-direction. Furthermore, the correlative approach allows the microstructures observed in polished sections to be confidently positioned in X, Y, Z space (using the known depth of sub-sampling down an orientated drill string). To advance exploration efforts, these orientated and positioned microstructural data can be pulled directly into implicit modelling software as an additional layer of information which could be readily analysed and iterated to finalize geological or grade interpretations. allows the microstructures observed in polished sections to be confidently positioned in X,Y,Z space (using the known depth of sub-sampling down an orientated drill string). To advance exploration efforts, these orientated and positioned microstructural data can be pulled directly into implicit modelling software as an additional layer of information which could be readily analysed and iterated to finalize geological or grade interpretations. The outcomes from this correlative study conducted at the Manondo-Choma gold prospect thus highlight that the gold is intimately associated with sulphide mineralization in NE-SW trending SE-dipping quartz-sulphide vein structures that show a dextral sense of shear. These structures were identified by considering the complementary 2D and 3D datasets, and the orientations are found to match those of the regional NE-SW magnetic trends identified using geophysics [32] (Figure 9a), and which are tentatively related to the southern Irumide orogeny. In addition to the value that the correlative approach unlocks for exploration at the front end of the mine value chain, it also reveals vital information that can assist during later processing and beneficiation of the ore assemblage. For example, XCT scanning provides definitive information on the size, distribution and mineralogical association of the gold, all of which can control the liberation efficiencies during minerals processing. Similarly, the distribution and semi-quantitative evaluation of the sulphide minerals (Table 2), offers a unique methodology to predict acid mine drainage capacity associated with an ore prospect, especially given that sulphide concentrations are difficult to measure using standard XRF quantification approaches [43]. Related to this, the complementary petrographic information related to sulphide speciation further has important implications for the presence of metal 'sweeteners' (e.g., Cu, Zn) that may add additional economic value to an ore prospect, and for the presence of any deleterious elements (e.g., As) that may require special consideration to ensure that the environmental impacts of downstream mining are minimized.

Recommendations for Further Development
The use of correlative microscopies in minerals exploration has not yet reached its full maturity, yet the example from the Manondo-Choma gold prospect clearly indicates its high potential for incorporation in routine minerals exploration workflows. This is especially pertinent for gold exploration endeavours given the nuggety distribution of the ore within mineralized vein sets. As bench-scale XCT machines become increasingly commercially available, it is anticipated that the approach will be adopted more extensively among exploration practitioners. However, in order to optimize the approach, the following recommendations apply: -Need for automated alignment algorithms based on optical image analysis methodologies to match 2D montages to their positioning in XYZ space. This will help as it will require minimal initial user input. The manual approach can take 4 h. -Need to develop software interfaces that allow XYZ positioned and orientated microstructural and mineralogical data (mapped in 3D space in XCT software) directly into common implicit and explicit geological modelling software. -Need to optimize voxel sizes to match fine grained gold mineralization. -Need to optimize the scanning parameters, or resolution of the sensors to better discern between different sulphide moieties using XCT technologies (and between gold and other dense minerals with high X-ray attenuation coefficients such as barite, galena, etc.). -Opportunities to develop other computed tomography techniques to provide detailed mineralogical, chemical and valence distributions in 3D space (synchrotron XRD computed tomography [44], synchrotron XRF CT [45], and synchrotron XANES CT respectively [46]).

Conclusions
We have highlighted a correlative approach in which 3D XCT data are combined with 2D mineralogical and chemical data to understand a complex orogenic gold mineralization. The use of XCT in combination with microscopy provides an excellent tool for understanding gold mineralization. Contextualizing 2D structural features in a 3D volume enables a deeper understanding of the structures that control mineralization and characterize the distribution of gold. Given that modern exploration efforts are targeting increasingly complex and low-grade ore deposits, the extra layers of information provided by using correlative 2D and 3D microscopies offer an exciting new tool to enhance and optimize minerals exploration workflows. It is anticipated that the approach will be adopted more extensively among exploration practitioners. Although the approach is especially well suited to gold mineralization (nuggety textures, high density and X-ray absorption coefficients of Au), it can also be applied to a range of other mineralization styles (e.g., disseminated sulphide mineralization, tungstate mineralization in skarns, etc.).