Carsten Laukamp

Airborne hyperspectral imaging of hydrothermal alteration zones in granitoids of the Eastern Fold Belt, Mount Isa Inlier, Australia

ABSTRACT Hyperspectral remote sensing data from the Eastern Fold Belt, Mount Isa Inlier, Australia were compared with petrographic and geochemical studies to map the spatial extension and compositional variations of Proterozoic granitoids and endoskarns as well as hydrothermal alteration patterns in adjoining metasedimentary successions. Detailed spatial analysis of spectral remote sensing data shows an almost circular alteration zoning in the Mallee Gap Granite, which was emplaced during a late phase of the Mesoproterozoic Williams event. A combination of hyperspectral images, such as white mica, kaolin and MgOH products, were used to map the alteration zoning. The formation of the endoskarn is presumably related to autometasomatism and interaction with fluids released from the country rocks during a late phase of the emplacement. The intrusion of the Mallee Gap Granite has only a local control on the hydrothermal alteration, but high potassic granites of the southern Mount Angelay Granite might have expelled oxidized mineralizing fluids and possibly had a major impact on regional scale alteration. Hyperspectral remote sensing data may be used to estimate the imprint of single igneous bodies on the Mesoproterozoic hydrothermal evolution of the Eastern Fold Belt and are important for the study of ore-forming hydrothermal processes in general.

Satellite-derived mineral mapping and monitoring of weathering, deposition and erosion

The Earth’s surface comprises minerals diagnostic of weathering, deposition and erosion. The first continental-scale mineral maps generated from an imaging satellite with spectral bands designed to measure clays, quartz and other minerals were released in 2012 for Australia. Here we show how these satellite mineral maps improve our understanding of weathering, erosional and depositional processes in the context of changing weather, climate and tectonics. The clay composition map shows how kaolinite has developed over tectonically stable continental crust in response to deep weathering during northwardly migrating tropical conditions from 45 to 10 Ma. The same clay composition map, in combination with one sensitive to water content, enables the discrimination of illite from montmorillonite clays that typically develop in large depositional environments over thin (sinking) continental crust such as the Lake Eyre Basin. Cutting across these clay patterns are sandy deserts that developed <10 Ma and are well mapped using another satellite product sensitive to the particle size of silicate minerals. This product can also be used to measure temporal gains/losses of surface clay caused by periodic wind erosion (dust) and rainfall inundation (flood) events. The accuracy and information content of these satellite mineral maps are validated using published data.

Hydrothermal mineral alteration patterns in the Mount Isa Inlier revealed by airborne hyperspectral data

High-resolution mineral maps derived from hyperspectral imaging (4.5 m pixel) enable the recognition of various types of hydrothermal alteration and the identification of fluid pathways. Airborne hyperspectral images from the Eastern Fold Belt of the Mount Isa Inlier were tested as a new tool for the detection of Fe-oxide Cu–Au (IOCG) related alteration. Four different types of hydrothermal alteration were identified with the hyperspectral mineral maps: (1) Metasomatic 1: white mica mineral maps show the spatial distribution of regional sodic–calcic alteration in metasedimentary successions of the Soldiers Cap Group in the Snake Creek Anticline. (2) Metasomatic 2: alteration zonation is evident from albitised granites assigned to the Williams–Naraku Suite along the Cloncurry Fault. These show characteristic absorption features in the shortwave infrared range (SWIR) which are depicted on the white mica mineral maps (white mica composition, white mica content, white mica crystallinity index). Alteration zonation in gabbros of the Cloncurry District was detected by a combination of MgOH and Fe2+ mineral maps (MgOH content, MgOH composition, amphibole/chlorite and Fe2+ and MgOH) combined with white mica mineral maps (white mica composition and white mica content). (3) Fluid channels 1: major fault zones, such as the Mt Dore fault zone in the Selwyn Corridor, are interpreted as important fluid pathways, where gradual changes in the mineral chemistry are highlighted with mineral maps (e.g. white mica content, white mica composition, white mica crystallinity index). (4) Fluid channels 2: MgOH and Fe2+ mineral maps were used to map breccia pipes in the northern Cloncurry District north of the Saxby Granite (Suicide Ridge). The MgOH and Fe2+ mineral maps were also used to distinguish various mafic rocks from amphibolites, which are host rocks for some of the IOCG deposits in the Eastern Fold Belt (e.g. Mount Elliott), and calcsilicate breccias pipes (e.g. Suicide Ridge).

Short-wavelength infrared spectroscopy of chlorite can be used to determine very low metamorphic grades

This study investigates the use of short-wavelength infrared reflectance spectroscopy for measuring chlorite physicochemical information of value for determination of very low metamorphic grades. The method is termed chlorite spectral index (CSI) and assessed through quantitative comparison with XRD-derived Arkai and Kubler indices (AI and KI, respectively) in a set of 41 pelitic samples. The four CSI measures investigated measure the depth (D) of the absorption features at 1900, 2250 and/or 2350 nm, and are defined as reflectance 2250D/1900D index (CSI(H 2 O)), reflectance (2350D + 2250D × 1.7)/1900D index (CSI(H 2 O) sum ), reflectance 2250D index (CSI(2250)) and reflectance 2350D + 2250D × 1.7 index (CSI sum ). The results demonstrate the applicability of the method: all CSI measures show a good correlation with the AI, and respective correlation factors range between R 2 = 0.840 (for the CSI(H 2 O)) and R 2 = 0.616 (for the CSI sum ). The CSI measures also show a good correlation with the KI with R 2 = 0.704 (for CSI sum ) to 0.769 (CSI(2250)). By applying the correlation equations to define the metamorphic zone boundaries for each of the CSI measures, 64.3–78.6 % of the samples match the metamorphic grade indicated by AI for the epizone, 47.1–70.6 % for the anchizone and 70–80 % for the diagenetic zone. All CSI measures show a good reproducibility, with an intra sample variance ranging between 2.31 % CSI(H 2 O) sum and 3.34 % (CSI(2250)). CSI is not restricted to pelitic rocks and marls, but can also be applied in very low-grade metabasites. In addition to the existing laboratory method, CSI has application potential in remote sensing applications.

Vibrational spectroscopy of calcic amphiboles – applications for exploration and mining

Calcic amphiboles in form of single crystals and in rock samples have been measured using laboratory-based infrared reflectance spectroscopy (IRS) and routine IRS technologies applied in mineral exploration. The composition of amphiboles in selected rock samples was validated with scanning electron microscopy (SEM) and electron microprobe work. Published values for wavenumber frequencies of hydroxyl-related stretching and bending vibrations were compared with the results from our study and both were combined to calculate combinations and overtones of [M1M1M3]-O-H in the short-wave infrared regions of 5000–4080 cm −1 (2000–2450 nm) and 7350 cm −1 (1360 nm) regions, respectively. Detailed comparison of major short-wave infrared absorption features in various calcic amphiboles and talc leads to the conclusion that an absorption feature centred at 2080 nm is diagnostic for talc and can be used to distinguish amphibole from talc. Multiple feature extraction scripts were developed to determine the relative abundance of amphibole and talc, as well as the Mg# of amphiboles in large IRS data sets. Our results show that only the 2390 nm absorption feature in amphibole can be reliably used to determine its abundance and Mg# in mineral assemblages containing other short-wave infrared active minerals. Different mafic and ultramafic lithologies can be inferred from infrared hyperspectral drill core logging and remote sensing datasets, based on the developed scripts.

Mineral system analysis of the Mt Isa–McArthur River region, Northern Australia

The Mt Isa–McArthur region is renowned for a range of commodities and deposit types of world-class proportions. The region is described here in the context of a ‘mineral system,’ through consideration of processes that operate across a range of scales, from geodynamics and crustal architecture, to fluid sources, pathways, drivers and depositional processes. The objective is to improve targeting of Pb–Zn, Cu and Cu–Au deposits. Repeated extension and high heat flow characterise much of the history prior to 1640 Ma. The pre-Barramundi Orogeny (pre-1.87 Ga) metamorphic basement was the substrate on which a volcanic arc developed, focussed along the Kalkadoon-Leichhardt Belt. This is related to an inferred east-directed subduction between 1870 and 1850 Ma. From 1755 to 1640 Ma, three successive volcano-sedimentary basins developed, the Leichhardt, Calvert and Isa Superbasins, in an interpreted distal back-arc environment. The Isan Orogeny, from 1640 to 1490 Ma, overlapped with Isa Superbasin sedimentation, suggesting a transition from back-arc to a foreland basin setting. Most crustal thickening occurred in the Eastern Fold Belt, an area earlier characterised by thinned crust and deep marine environments. This region was deformed into nappe-like structures with high-temperature–low-pressure regional metamorphism and associated granites; the latter are absent from the Western Fold Belt. Metal deposition mainly occurred late in the history, with all known (and preserved) major base metal occurrences either hosted by Isa Superbasin rocks or formed during the Isan Orogeny. Earlier superbasins were potential fluid source regions. Sedimentary formation waters, metamorphic and magmatic fluids were present at prospect scale, while meteoric and possibly mantle sources are also implicated. The spatial distribution of metallogenic associations (i.e. iron oxide–copper–gold, Pb–Zn–Ag, U, Au) across the inlier may result from differences in the geodynamic make-up and evolution of the pre-1.87 Ga tectonic elements. Penetrative faults are interpreted as predominantly steeply dipping and to have acted as pathways for fluids, both in extension and compression. Fluid mixing was a potentially significant ore deposit control. Examples are drawn from the Ernest Henry iron oxide–copper–gold-related hydrothermal breccias in the east and from the Mt Isa Copper deposit in the west. Stress switching during late-stage deformation appears to have triggered a fluid mixing event that led to formation of the major copper deposits.

Regolith-geology mapping with support vector machine: A case study over weathered Ni-bearing peridotites, New Caledonia

Abstract Accurate maps of Earth’s geology, especially its regolith, are required for managing the sustainable exploration and development of mineral resources. This paper shows how airborne imaging hyperspectral data collected over weathered peridotite rocks in vegetated, mountainous terrane in New Caledonia were processed using a combination of methods to generate a regolith-geology map that could be used for more efficiently targeting Ni exploration. The image processing combined two usual methods, which are spectral feature extraction and support vector machine (SVM). This rationale being the spectral features extraction can rapidly reduce data complexity by both targeting only the diagnostic mineral absorptions and masking those pixels complicated by vegetation, cloud and deep shade. SVM is a supervised classification method able to generate an optimal non-linear classifier with these features that generalises well even with limited training data. Key minerals targeted are serpentine, which is considered as an indicator for hydrolysed peridotitic rock, and iron oxy-hydroxides (hematite and goethite), which are considered as diagnostic of laterite development. The final classified regolith map was assessed against interpreted regolith field sites, which yielded approximately 70% similarity for all unit types, as well as against a regolith-geology map interpreted using traditional datasets (not hyperspectral imagery). Importantly, the hyperspectral derived mineral map provided much greater detail enabling a more precise understanding of the regolith-geological architecture where there are exposed soils and rocks.

Remote Sensing-Based Exploration of Structurally-Related Mineralizations around Mount Isa, Queensland, Australia

Hyperspectral imaging is a powerful tool for mineral mapping and increasingly used in poorly-accessible areas. It only requires a limited amount of validation sample points, but can fail to discriminate spectrally-similar features. In this manuscript, we show that we improve the identification of interesting targets by including geomorphological data in the spectral mapping scheme. We jointly use geomorphic and spectral features to locate gossanous ironstone ridges as an indicator for possible Pb-Zn-Ag-mineralization and provide an application around Mount Isa and George Fisher/Hilton mine, Queensland, Australia. We combine hyperspectral HyMap data using mixture tuned matched filtering with topographical indices, such as maximum curvature and the topographical position index. As it is often the case with structurally-controlled mineralization, the amount of training sites is limited, and supervised classification methods cannot be implemented. Therefore, we implement expert knowledge in a decision tree to take advantage of the relationship between mineralization, alteration and structure. Optimized rock sampling and spectral measurements provided data for validation. We are able to map sets of gossanous ridges with a minimum of validation points, not only within the Mount Isa mining area itself, but also outside the commonly-accepted host rocks. The ridges are parallel to north-south trending geomorphological features and probably associated with the Paroo fault zone. Similarities between the ridges were confirmed by field observations, spectral measurements and a qualitative rock sample analysis. We identified new mineralized ridges that we could subsequently attribute to a poorly-known and sub-economic deposit known as the Mount Novit Pb-Zn-deposit.

Multiple methods for regional- to mine-scale targeting, Pataz gold field, northern Peru

Gold production in the Pataz district, northern Peru, is derived from mesothermal veins hosted by the Pataz batholith and basement-hosted epithermal and carbonate–base metal veins. At the regional scale, processing of Advanced Spaceborne Thermal Emission and Reflection Radiometer data can be used to delineate district-scale argillic alteration. One such area extends for tens of kilometres NNW of Vijus in the Maranon Valley. At the southern end of this area, basement-hosted quartz–carbonate–sulfide veins in faults support artisanal gold-mining operations. SEM analyses show that the alteration envelopes around these faults are dominated by illitic clays. These artisanal gold workings highlight the economic potential of the largely unexplored parts of the district-scale argillic alteration zone, further north. At the district scale, paleostress modelling maps areas of low minimum stress during Carboniferous ENE–WSW shortening, based on a new 1:25 000 geological map of the Pataz district. The resulting distribution of low minimum stress is used to predict sites of rock fracture under high fluid pressure, and consequent vein formation. These areas of low minimum stress occupy 11% of the modelled area but contain 50% of the known veins in the Pataz district. Some areas of low minimum stress contain no known veins, and where these are poorly explored or poorly exposed, they are proposed as potential targets for gold exploration. In combination with SEM microanalysis, ASD hyperspectral reflectance analysis of drill core samples shows that visible proximal sericitic alteration around batholith-hosted auriferous veins is predominantly phengitic illite. Automated software interpretation of ASD reflectance spectra using The Spectral Assistant shows that sericite in cryptic alteration distal from auriferous veins varies from mainly illite adjacent to the phengitic illite zone, to more distal muscovite. Reactivation of faults and mineralised vein contacts during the largely Cenozoic Andean orogeny produced chlorite alteration that locally overprints proximal phengitic illite alteration. ASD spectrometry identifies relict phengitic illite in some chloritic alteration zones and thus indicates proximity to mineralised veins at the deposit scale. Elevated pathfinder element concentrations within proximal phengitic illite alteration zones around batholith-hosted veins do not extend more than a few metres beyond the visible alteration envelope. The alkali alteration index [(Rb + Cs)/Th]N is elevated above background levels for up to 15 m beyond the visible sericite alteration zone in one of two holes investigated. In the other hole, both [(Rb + Cs)/Th]N and 3K/Al can be used as a lode-scale vector to gold-bearing veins within broad intersections of visible sericite alteration.

Reflectance spectroscopic characterisation of mineral alteration footprints associated with sediment-hosted gold mineralisation at Mt Olympus (Ashburton Basin, Western Australia)

ABSTRACT Hydrothermal ore deposits are typically characterised by footprints of zoned mineral assemblages that extend far beyond the size of the orebody. Understanding the mineral assemblages and spatial extent of these hydrothermal footprints is crucial for successful exploration, but is commonly hindered by the impact of regolith processes on the Earths surface. Hyperspectral drill core (HyLogger™-3) data were used to characterise alteration mineralogy at the Mt Olympus gold deposit located 35 km southeast of Paraburdoo along the Nanjilgardy Fault within the northern margin of the Ashburton Basin in Western Australia. Mineralogy interpreted from hyperspectral data over the visible to shortwave (400–2500 nm) and thermal (6000–14500 nm) infrared wavelength ranges was validated with X-ray diffraction and geochemical analyses. Spaceborne multispectral (ASTER) and airborne geophysical (airborne electromagnetic, AEM) data were evaluated for mapping mineral footprints at the surface and sub-surface. At the deposit scale, mineral alteration patterns were identified by comparing the most abundant mineral groups detected in the HyLogger data against lithology logging and gold assays. Potential hydrothermal alteration phases included Na/K-alunite, kaolin phases (kaolinite, dickite), pyrophyllite, white mica, chlorite and quartz, representing low-T alteration of earlier greenschist metamorphosed sediments. The respective zoned mineral footprints varied depending on the type of sedimentary host rock. Siltstones were mainly characterised by widespread white-mica alteration with proximal kaolinite alteration or quartz veining. Sandstones showed (1) distal white mica, intermediate dickite, and proximal alunite + kaolinite or (2) widespread white-mica alteration with associated intervals of kaolinite. In both, sandstones and siltstones, chlorite was distal to gold mineralisation. Conglomerates showed distal kaolinite/dickite and proximal white-mica/dickite alteration. Three-dimensional visualisation of the gold distribution and spatially associated alteration patterns around Mt Olympus revealed three distinct categories: (1) several irregular, poddy, SE-plunging zones of >0.5 ppm gold intersected by the Zoe Fault; (2) sulfate alteration proximal to mineralisation, particularly on the northern side of the Mt Olympus open pit; and (3) varying AlIVAlVISiIV–1(Mg,Fe)VI–1 composition of white micas with proximity to gold mineralisation. Chlorite that developed during regional metamorphic or later hydrothermal alteration occurs distal to gold mineralisation. ASTER mineral mapping products, such as the MgOH Group Content used to map chlorite (±white mica) assemblages, showed evidence of correlation to mapped, local structural features and unknown structural or lithological contacts as indicated by inversion modelling of AEM data.