Interpreting the Shortwave Infrared & Thermal Infrared Regions of Remote Sensed Electromagnetic Spectrum with Application for Mineral-Deposits Exploration ()
1. Introduction
Since the appearance of ASTER in Dec. 1999, the use of ASTER data in mineral prospecting has increased in recent years because of its relatively low cost, broad coverage, and unique integral bands highly sensitive to alteration minerals (i.e. minerals known to surround target minerals). The ASTER data set contains visible, shortwave infrared and thermal bands. Because the spectral bands of ASTER are sensitive to absorption of some minerals. The proper preprocessing and combination of these bands can produce relative mineral alteration distributions such as illite, alunite, kaolinite, carbonates and so on. Kanlinowski and Oliver [1] offered a suite of standard ASTER alteration enhancements that can be used to broadly delineate these alteration mineral zones or assemblages.
Since then the remote sensing anomaly technology has been spread throughout the country. By the same time our research effort was concentrated on the application of ASTER data for resolving the problems, left with ETM, on the classification of remote sensing anomalies, on the anomaly extraction of new mineral types. We have developed a simplified algorithm for temperature and emission separation using de-correlation stretch function, proposed at first by Gilespie, et al. in 1998 [2]. Up to now we’ve got anomalies of mineral 12 types [3]. Here the anomaly extraction of mineral deposits will be explored. In comparison with the works by others in this field, our new achievements comprise the extraction of anomalies for Ni-Cu deposit, gypsum, and borax.
2. Spectroscopy for Remote Sensing in Geological Application
2.1. Visible-Reflected Infrared Remote Sensing
Spectroscopy for RS may be defined as the study of the interaction of earth matter with electromagnetic radiation. The regions of the RS EMR for geological work with goal in detection and exploration of mineral resources are visible (0.35 - 0.70 μm), near-infrared and SWIR (0.7 - 2.5 μm), TIR (8 - 13 μm) (see Table 1). Within these ranges, intrinsic interactions are caused by electronic and/or vibrational processes.
For geologic materials, electronic processes provide information about the presence and nature of a particular iron and its environment, about adjacentions, or about defects in specific crystal locations. Thus, information concerning bulk composition is of an indirect nature.
Vibrational processes are governed by three characteristic properties of the material: the chemical composi- tion, the geometry and equilibrium positions of the constituent atoms, and potential field of the interatomic
Table 1. Discrimination of deposit types by multi-channel r.s.
forces representative of inter-atomic bonds and angles. Consequently the information available is directly related to bulk composition.
Rocks are basically assemblages of minerals. So much of our referred spectra presented by USGS and ASU are mostly for pure minerals since the spectrum of a rock is typically some composition of the spectra of its constituent minerals.
So the inverse problem is usually quite complicated. Nevertheless we have succeeded in extracting RS anomalies of 12 different types of mineral deposits (Table 1).
2.2. Thermal Infrared Remote Sensing
Thermal infrared spectral range of 8 - 13 μm, the thermal infrared spectra in the ASTER data include information of temperature and emissivity. In order to geological applications, it’s necessary to get the thermal infrared spectral emissivity. The method “discrimination method of extremes on the spectral shape” was applied based on de-correlation stretch (or transformation) in this paper. The information of carbonate, sulphate, and ultrabasic rocks were extracted.
3. Processing Techniques of ASTER Images
The ASTER data (all of the Visual-near Infrared, the Short Wave IR and the Thermal IR) were used for the extraction of R.S. anomaly for the hydroxyl minerals in hydrothermal alteration zone of metallic ore deposit , the carbonaceous (marble, smithsonite, dolomite, calcite) and sulfates (gypsum, celestine) sediments.
3.1. De-Interfered Anomalous Principal Component Thresholding Technique
The early use of PCA for extraction of TM anomaly was made by Crosta and Moore [4] and Loughlin [5]. Then we expanded its use for ASTER data [6]-[9]. An advanced system for extraction of hydrothermal alteration assemblages from ETM+ and Aster data, “De-interfered Anomalous Principal Component Thresholding Technique”, is illustrated in Figure 1. The RM, SAM and PCA are compared by the analysis of their principles and the modeling test. The pre- and post-processing have been done for limiting the interference of water, ice, cloud, snow, vegetation, relief shadow, cloud shadow, salt-lick.
Figure 1.De-interfered anomalous principal component thresholding technique.
3.2. Discrimination Method of Extremes on the Spectral Shape
The ASTER data (all of the Visual-near Infrared, the Short Wave IR and the Thermal IR) were used for the extraction of R.S. anomaly for the carbonaceous (marble, smithsonite, dolomite, calcite) and sulfates (gypsum, celestine) sediments. The main processing techniques used are:
・ The De-correlation Stretch (or transformation);
・ The Band-Selective Principle Component Analysis (PCA);
・ The Combinatorial(combinational) Band Ratio Analysis;
・ The Relative Band Depth (RBD);
・ The Piece-Wise Thresholding.
The De-correlation stretch is especially effective on ASTER-TIR images in which channels are highly correlated by temperature. De-correlation is used for separation of temperature and thermal emission (see Figure 2).
4. The Discovery of Luodong Ni Deposit
Using these techniques we have inversed quite a lot science of ETM + and ASTER data and gained high quality result. For example we take the discovery of Luodong Ni deposit, the determination of gypsum, and the determination of borax .
4.1. The Discovery of Luodong Ni Deposit
Figure 3 shows the complex information for discovery of Luodong deposit in Ruoqiang county of Xinjiang. For discrimination of the ETM+ anomaly characteristic, the SAM was used in interpretation. The spectrum data of Poshi known small Ni deposit was taken as a reference spectrum. The inversion result was verified by field work, PIMA determination and electro-probing.
4.2. The Determination of Gypsum
Owing to the severe shortage of potash fertilizer in China, great attention has been paid to make full use of achievements of geophysics and remote sensing in exploration of potash salt in China during the past more than 50 years. Because of the lacking of diagnostic spectral property for potash salt, we tried to use the gypsum for indirect path for this problem.
Figure 3. Complex information for discovery of Luodong deposit. (A) RS image; (B) RS anomaly image; (C) Spectra; (D) Photo of rock; (E) Simplified geological map; (F) Microscopic photo of olivine; (G) Microscopic photo of gabbro; (H) PIMA result.
The specific spectral characteristic for gypsum shown in Figure 4 is displaying one example of determination of gypsum by processing of ASTER-TIR data using flow-chart in Figure 2. This interpretation method was used for Qiangtang Basin, where the condition is quite severe for field geologic work. The average elevation is >5000 m; the air is very thin; it is extremely cold. The basin is almost a restricted zone. So the application of remote sensing for geological work in this kind of area is especially important.
Figure 4. (A) ASTER anomaly image; (B) ASTER image (7, 3, 1); S1. Photo, taken of point s
1 in
(A), showing gyosum; S2. Photo of point S
2 in
(A), showing gypsum covered by mud.
4.3. The Determination of Borax
The usage of borax on Qinghai-Tibet Plateau for medical, glazier and other purposes takes special role. The exploitation of borax is one of the important component parts of Tibetan industry all along. BanKogCo was taken as a known borax deposit for the application of ASTER data. BanKogCo belongs to the moderately mineralized carbonate type. It is located
10 km
in the east from SilingCo or 120 km in the north-east from the famous NamCo as the biggest saline lake in Tibet. BankogCo consists of 3 small lakes. The sedimentary sub-environment in BankogCo district is shown in Figure 5 [10].
The Combinatorial Band-Ratio Analysis of the SWIR was studied for the borax anomaly and to avoid the false anomaly from carbonates and argillaceous sediments. Besides BanKogCo the borax anomaly appears at the east of NamkaCo. This region is an undiscovered area. Since the anomaly is even stronger than that for BanKogCo, it is worth carrying on exploration after some ground truth investigations (Figure 6).
Figure 5. Sedimentary sub-environment in BanKogCo district [10]. Perennial brackish lake; 2-Perennial saline-lake; 3-Ephemeral lake; 4-Salt pan; 5-Salt mud-flat; 6-Spring pond-marsh; 7-Sand flat; 8-Alluvial fan; 9-Bench and strand; 10-Bench pad and platform bench; 11-Bed rock.
Figure 6. Results of processing of ASTER image for BanKogCo.
Acknowledgements
This research was jointly supported by National Natural Science Foundation of China (No. 41102205), Special Funds Projects for Basic Scientific Research Business Expenses of Mineral Resources Research Institutes in Chinese Academy of Geological Sciences (Nos. K1208, K1213). We thank greatly Prof. Yang Jianmin for Field work in discovery of Luodong Ni Deposit.
NOTES
*Corresponding author.