ASTER-Based Remote Sensing Image Analysis for Prospection Criteria of Podiform Chromite at the Khoy Ophiolite (NW Iran)
Total Page:16
File Type:pdf, Size:1020Kb
minerals Article ASTER-Based Remote Sensing Image Analysis for Prospection Criteria of Podiform Chromite at the Khoy Ophiolite (NW Iran) Behnam Mehdikhani and Ali Imamalipour * Department of Mining Engineering, Urmia University, Urmia 57561-51818, Iran; [email protected] * Correspondence: [email protected] Abstract: A single chromite deposit occurrence is found in the serpentinized harzburgite unit of the Khoy ophiolite complex in northwest Iran, which is surrounded by dunite envelopes. This area has mountainous features and extremely rugged topography with difficult access, so prospecting for chromite deposits by conventional geological mapping is challenging. Therefore, using remote sensing techniques is very useful and effective, in terms of saving costs and time, to determine the chromite-bearing zones. This study evaluated the discrimination of chromite-bearing mineralized zones within the Khoy ophiolite complex by analyzing the capabilities of ASTER satellite data. Spectral transformation methods such as optimum index factor (OIF), band ratio (BR), spectral angle mapper (SAM), and principal component analysis (PCA) were applied on the ASTER bands for lithological mapping. Many chromitite lenses are scattered in this ophiolite, but only a few have been explored. ASTER bands contain improved spectral characteristics and higher spatial resolution for detecting serpentinized dunite in ophiolitic complexes. In this study, after the correction of ASTER data, many conventional techniques were used. A specialized optimum index factor RGB (8, 6, 3) Citation: Mehdikhani, B.; was developed using ASTER bands to differentiate lithological units. The color composition of band Imamalipour, A. ASTER-Based ratios such as RGB ((4 + 2)/3, (7 + 5)/6, (9 + 7)/8), (4/1, 4/7, 4/5), and (4/3 × 2/3, 3/4, 4/7) produced Remote Sensing Image Analysis for the best results. The integration of information extracted from the image processing algorithms Prospection Criteria of Podiform used in this study mapped most of the lithological units of the Khoy ophiolitic complex and new Chromite at the Khoy Ophiolite (NW Iran). Minerals 2021, 11, 960. prospecting targets for chromite exploration were determined. Furthermore, the results were verified https://doi.org/10.3390/min11090960 by comprehensive fieldwork and previous studies in the study area. The results of this study indicate that the integration of information extracted from the image processing algorithms could be a broadly Academic Editors: Amin Beiranvand applicable tool for chromite prospecting and lithological mapping in mountainous and inaccessible Pour, Omeid Rahmani and regions such as Iranian ophiolitic zones. Mohammad Parsa Keywords: ASTER; chromite; Khoy ophiolite; spectral angle mapper (SAM); band ratio; principal Received: 19 July 2021 component analysis (PCA) Accepted: 19 August 2021 Published: 2 September 2021 Publisher’s Note: MDPI stays neutral 1. Introduction with regard to jurisdictional claims in The mapping of ophiolite sequences has become a research interest of scientists and published maps and institutional affil- iations. exploration geologists in the world because they host economic minerals such as chromium, copper, manganese, gold, nickel, barium, lead, and zinc [1–3]. Ophiolitic ultramafic rocks are the hosts of podiform chromite deposits. Podiform chromite deposits are small magmatic chromite bodies formed in the lower level of an ophiolite complex. Podiform chromite mines have produced 57.4% of the world’s total chromite production [4]. Ophiolite Copyright: © 2021 by the authors. zones in Iran are widespread and are often found in different locations with varying Licensee MDPI, Basel, Switzerland. geologic and tectonic settings. The Khoy ophiolite complex is a part of the Tethyan ophiolite This article is an open access article belt, and it is one of the largest Iranian ophiolite complexes, covering a widespread area in distributed under the terms and conditions of the Creative Commons northwest Iran along the Iran–Turkey border and continuing toward western Turkey [5,6]. 2 Attribution (CC BY) license (https:// Ultramafic rocks, which are often serpentinized, are widespread in 250 km of the Khoy creativecommons.org/licenses/by/ ophiolite [5,6]. The Khoy ophiolite is one of the most promising areas for prospecting 4.0/). chromite deposits because of extensive outcrops of ultramafic rocks. So far, more than Minerals 2021, 11, 960. https://doi.org/10.3390/min11090960 https://www.mdpi.com/journal/minerals Minerals 2021, 11, 960 2 of 17 20 chromite deposits have been identified in this area. These chromite occurrences have lenticular, tubular, and vein-like shapes hosted by serpentinized harzburgite. The chromite deposits in the Khoy ophiolite can be clearly classified into two groups: high-Al chromites (Cr# = 0.38–0.44) from the eastern ophiolite, and high-Cr chromites (#Cr = 0.54–0.72) from the western ophiolite [5,6]. Most Iranian ophiolitic zones are located in mountainous and inaccessible regions. Thus, prospecting for chromite deposits with geological mapping is challenging and time-consuming. Remote sensing analysis plays an important role in the exploration of mineral deposits, as well as in lithological mapping and detection of associated hydrothermal mineraliza- tion, in Iran. The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) is an advanced multispectral satellite imaging system that has created new tools for the mapping of geological structures and detecting certain alteration minerals or assem- blages [2,3,7]. The ASTER sensor launched the TERRA platform in December 1999. The ASTER plat- form travels in a near-circular, sun-synchronous orbit with an inclination of approximately 98.2◦, an altitude of 705 km, and a repeat cycle of 16 days, offering relatively improved spatial, spectral, and temporal resolutions. It is made from three visible and near-infrared spectral bands (VNIR, between 0.52 and 0.86 µm, with 15-m spatial resolution) and infrared, reflecting radiation in six short-wavelength infrared spectral bands (SWIR, between 1.6 and 2.43 µm, with 30 m spatial resolution). Sensor characteristics of the ASTER instruments are shown in Table1[8,9]. Table 1. Sensor characteristics of ASTER instruments [8]. ASTER Sensor Characteristics VNIR SWIR TIR Band01 0.52–0.60 Band04 1.6–1.7 Band10 8.125–8.475 Band02 0.63–0.69 Band05 2.45–2.185 Band11 8.475–8.825 Spectral bandswith range Band03N 0.76–0.86 Band06 2.185–2.225 Band12 8.925–9.275 (µm) Band03B 0.76–0.86 Band07 2.235–2.285 Band13 10.95–10.95 Band08 2.2295–2.365 Backward-looking Band14 10.95–11.65 Band09 2.36–2.43 Spatial resolution (m) 15 30 90 Swathwidth (km) 60 60 60 Due to the great extent of ultramafic rocks, which are the host of chromite deposits in the Khoy ophiolite, the possibility of discovering new chromite deposits is high and more exploration and investigation is needed. Given the extremely rugged topography with difficult access, new exploration methods such as the remote sensing method can be useful for this purpose. The present study evaluates the discrimination of chromite-bearing mineralized zones within the Khoy ophiolite complex by analyzing the capabilities of ASTER satellite data. Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data can easily separate various rock units, the extent of the ultramafic rocks, and it can provide detailed geological maps of the area [8,9]. The extraction of spectral information related to ophiolite mapping can be achieved through image processing techniques such as band ratio (BR) and principal component analysis (PCA) on ASTER bands [8,10]. The color composition of the band ratio (4/1, 4/5, 4/7) is an effective means of determining the lithological ophiolite complexes [11]. Principle component analysis and band ratio methods are very useful for determining the serpentinized dunite that is the host of the chromite veins [7,11–14]. Abdeen used ASTER spectral band ratios RGB color composite of 4/7, 4/1, 2/3, 4/3, and RGB (4/7, 3/4, 2/1) for mapping ophiolitic units, metasediments, volcanoclastic, and granitoids in the southeastern desert of Egypt [15]. Amer used principal component analysis of ASTER data to determine the lithologic units of the ophiolite complexes in Pakistan. In the eastern ophiolites of Egypt, Amer used band ratios of (7 + 9)/8, (5 + 7)/6, (2 + 4)/3, Minerals 2021, 11, x FOR PEER REVIEW 3 of 18 Principle component analysis and band ratio methods are very useful for determin- ing the serpentinized dunite that is the host of the chromite veins [7,11–14]. Abdeen used ASTER spectral band ratios RGB color composite of 4/7, 4/1, 2/3, 4/3, and RGB (4/7, 3/4, Minerals 2021, 11, 960 2/1) for mapping ophiolitic units, metasediments, volcanoclastic, and granitoids in the 3 of 17 southeastern desert of Egypt [15]. Amer used principal component analysis of ASTER data to determine the lithologic units of the ophiolite complexes in Pakistan. In the eastern ophiolites of Egypt, Amer used band ratios of (7 + 9)/8, (5 + 7)/6, (2 + 4)/3, and PCA (4,5,2) andfor the PCA lithological (4,5,2) for mapping the lithological of several mapping units [7]. of severalHashem units and [Pournamda7]. Hashemri andconducted Pournamdari conductedresearch using research ASTER using data ASTERon the Abdasht data on ophiolites the Abdasht in northeastern ophiolites inIran northeastern [12]. Thermal Iran [12]. Thermalinfrared (TIR) infrared bands (TIR) in the bands thermal in therange thermal of spectral range absorption of spectral can absorptionbe used for the can detec- be used for thetion detection of silicate offormations silicate formations [16]. [16]. 2.2. DescriptionDescription of of the the Study Study Area Area The Khoy Khoy ophiolite ophiolite covers covers an anarea area of about of about 3900 3900 km² kmin northwest2 in northwest Iran along Iran the along the Iran–TurkeyIran–Turkey boundary. This This ophiolitic ophiolitic complex complex is is limited limited on on the the west west and and north north by by the the Iran– TurkeyIran–Turkey border border and and on theon the east east and and south south by by a a southeastern-northeastern southeastern-northeastern fault fault (Fig- (Figure1 ).