1. Introduction
Ophiolites represent fragments of oceanic crust and upper mantle that were tectonically emplaced onto continents during orogenic events [
1,
2]. Ophiolitic sequences are generally exposed along suture zones and can be used to reconstruct the paleo-geography of orogenic events [
2,
3]. In general, they host various mineral deposits, including podiform chromite deposits (PCD), which are small magmatic bodies unevenly distributed in the ultramafic section of an ophiolite complex. PCDs have been found in the mid-oceanic ridge, off-ridge, and, especially, in suprasubduction tectonic settings. Most of them are found near the mantle transition zone, within ultramafics (mainly dunite) of the ophiolitic sequence, and at the contact of the upper cumulate zones [
4].
While the exploitation of the surficial deposits has a long history, with the development of high-precision and successful methods for mineral exploration, the geophysical techniques had a low success rate and suffered from limitations, such as the small size of chromite pods, the physical similarities with certain rocks, structural features, and presence of iron-rich bodies [
4,
5] Most surficial podiform chromite deposits have been exploited due to their resistance to weathering compared to their serpentinized host, which allows them to stand higher than the surrounding rocks. However, these deposits tend to be found in clusters so that additional orebodies may exist deep below the known deposits. Despite this, there are no proven exploration guidelines for finding associated ore bodies. Since there are rarely geochemical halos around deposits, geochemical approaches have typically been unsuccessful in detecting blind chromite deposits. Furthermore, due to the complexity of the involved geologic factors, no single method of exploration can be universally applied to finding these deposits. Integrated methods, instead, were somewhat successful for chromite exploration; for example, the combination of geological, gravity, magnetic, and electrical surveys [
5].
Satellite remote sensing has proven capabilities in locating exploration targets through the discrimination of hydrothermal alteration zones, structures, and lithologies [
6,
7]. With the advancement of satellite missions and their wide applications for remote reconnaissance of mineralization signs, podiform chromite exploration has benefited from some developed methods [
8]; however, the intrinsic characteristics of podiform chromites still have been a barrier as, for instance, the lack of extensive hydrothermal alteration related to mineralization. Most remote sensing efforts on ophiolites have been concentrated on lithological mapping using conventional methods, including band rationing, Principal Component Analysis (PCA), and spectral transform techniques [
9,
10,
11]. Some authors used machine-learning methods with more potential and accuracy to differentiate chromite host lithological units [
12,
13], particularly when applied to multiple sensors [
14,
15,
16]. However, regarding the low spatial resolution of almost all freely available satellite data compared to the dimension of chromite pods, the direct detection of chromite mineralization is rarely possible.
On the other hand, with the emergence of Unmanned Aerial Vehicles (UAVs) in the mining industry, several applications have been developed to help mineral exploration and geological mapping [
17,
18]. UAV photogrammetry has been proposed to be an efficient tool for the exploration of Cu veins through precise geological mapping [
19], semi-automated rock mass characterization [
20], supervised lithological classification [
21], structural analysis [
22] and semi-automatic lineation extraction [
23,
24]. The ultra-resolution orthomosaics (typically 3–5 cm) acquired from UAV mapping allow the geologists to manually draw the lithological boundaries, hydrothermally altered areas, and structures. However, none of the abovementioned methodologies could produce reliable results without adequate field mapping and analytical and microscopic examinations. Despite the potential of satellite and UAV data synergies [
25] for geological remote sensing applications, there is little research involving both datasets, mainly related to landslide mapping and monitoring [
26] and mineral exploration [
27]. Although UAV-mounted multi-to-hyperspectral sensors provide high-resolution data [
28], the expensive equipment still represents a barrier to developing this groundbreaking exploration method.
The present research, for the first time, aims to build a bridge among UAV photogrammetry, satellite remote sensing, and field data, to constrain better the criteria for finding chromite-bearing mineralized zones. This approach becomes essential for areas having rough topography, such as those analyzed here as a case study (Sabzevar Ophiolite, NE Iran). To achieve this objective, we followed a multi-stage processing method. Firstly, we looked for the most potential lithological units (i.e., serpentinized ultramafic units) hosting chromite mineralization using satellite remote sensing. Then, ultra-resolution RGB orthomosaic acquired from UAV mapping is used to manually differentiate geological units and structures. After verifying results by field checking, the support vector machine was applied to classify lithological units, peridotites with a variable range of serpentinization, and serpentinized dunites. Finally, to evaluate the validation of results, we carried out prospecting and field surveying together with chemical and microscopic analysis. After combining the results, we were able to provide a geological map with more exploration guidelines for prospecting chromite.
2. Geological Background
The present geological configuration of Iran is the result of several geodynamic processes that have operated since the Paleozoic. Among the most important, the closure of the PaleoTethys and Neotethys oceans and the consequent accretion of Gondwana fragments (Cimmerian blocks) over the southwestern flank of the Asia plate has been the main process of orogenesis, continuing today [
29]. Ophiolitic sequences register the timing and modality of oceanic closure, and this is the case of Iranian ones, which belong to the Middle East Tethyan ophiolites. The case study of the present paper regards the western part of Sabzevar Ophiolite (NE Iran), which hosts several known chromite occurrences and mines (
Figure 1). The Sabzevar suture zone is made up of a dismembered ophiolitic complex, extending for more than 400 km along an E-W direction. It is bordered at the north by the major Sangbast–Shandiz strike–slip fault system and at the south by the major Dorouneh sinistral strike–slip fault, which separates Sabzevar from the Lut Block (
Figure 1). The Sabzevar ophiolite would then represent the suture zone of the northern branch of the Neothethys “Sabzevar ocean,” which opened and closed during the Late Cretaceous [
1].
Although highly dismembered and tectonized, the Sabzevar ophiolites show a typical oceanic tectonic melange [
1,
30], representing a complete lithospheric section and containing all the rocks representative of an ophiolitic sequence. In the Sabzevar ophiolitic melange, lithological units and formations can be firstly divided into the pre- (Late Cretaceous–Early Paleocene) and post-ophiolitic sequence (Eocene to Plio-Quaternary arc). Mantle peridotites, cumulate ultramafics, gabbros, diabasic dykes, and volcanic, sedimentary, and metamorphic rocks constitute the ophiolitic sequence [
31]. Mantle peridotites, including lherzolite, harzburgite, dunite, and serpentinite, are the prevalent rock type in the Sabzevar ophiolite. All lithologies in the mantle section are cut by many rodingitized diabase, pegmatitic gabbros, micro-gabbros, and amphibolitic dykes [
31,
32]. There are large masses of podiform chromitite within the mantle sequence [
33], well-developed in the central and western parts of the Sabzevar ophiolite (
Figure 1). In some regions, chromites with both low #Cr and high #Cr have been found, suggesting crystallization of chromites from early MORB-type and then bainitic-type melts [
34]. There are also many ancient and recent mines in the region, but most are inactive due to the small deposit size and or inappropriate exploration work. The crustal section of the Sabzevar ophiolite includes cumulate peridotites, gabbroic rocks, and plagiogranites. Gabbroic rocks are largely found as isolated outcrops and are divided into isotropic to pegmatitic lenses within mantle harzburgites and crustal cumulate layered gabbro/gabbronorite/leucogabbro. The diabasic dykes occur as both isolated intrusive bodies, which cut mantle and plutonic crustal sequences, and sheeted dyke complexes. The extrusive sequence of the Sabzevar ophiolite consists of pillow lavas and sheet flows, intercalated with oceanic fossiliferous pelagic limestone and volcano-sedimentary rocks [
30]. Geochronological data show that the Sabzevar crustal sequence formed over a considerable period in the Late Cretaceous time ([
30,
31,
32] and references therein).
3. Materials and Methods
3.1. Exploration Model
Choosing exploration criteria to develop a model for remote sensing of the PCDs could be ambiguous due to their natural characteristics and the high level of uncertainty. Unlike most mineral deposits with a clear and somewhat zoned hydrothermal alteration halo, these deposits generally do not present such strong absorption features, which are needed as an indicator for use in spectral remote sensing. Moreover, PCD spatial distribution in the ultramafic rocks is irregular, not allowing an exploration pattern recognition. Within this context of difficulties, in the present paper, three assumptions were adopted for the remote sensing study. First, completely covering of chromitite pods by a dunitic layer is a widespread phenomenon, although not all dunites contain mineralization [
4]. We assume that highly serpentinized areas might host chromite pods since these probably have dunitic protolith. Second, the emplacement of a significant share of PCDs at or near (within 1 km) the Moho discontinuity, or transition zone between the overlying cumulate and underlying tectonite zones of ophiolite sequences, which narrows the geologic target [
4,
35,
36]. However, recognition of these zones is not always clear, especially in tectonic mélange zones. The third assumption is that some podiform chromite outcrops have a close spatial relationship with the mafic dykes intruded within serpentinized peridotites, compositionally ranging from gabbro to pyroxenite [
36]. With these limitations in mind, we designed a multi-stage remote sensing method whose details are schematically drawn in
Figure 2.
The whole experiment was divided into six steps, namely: (i) multi-sensor satellite data processing, including investigation of temporal changes over the period of 35 years; (ii) target detection using the area of chromite open pits as reference spectra; (iii) spectral mapping of chromite and associated minerals; (iv) UAV mapping and processing; (v) visual interpretation of various lithological units, field surveying, and sampling; (vi) lithological mapping by a support vector machine algorithm (SVM), accuracy evaluation and assessment.
We found a method for grading the degree of serpentinization in a wide area of the ophiolitic zone using satellite images. Then, with the aim of ultra-resolution orthomosaic obtained from UAV mapping, we verified the results of remote sensing and manually traced the geological units whose extent is beyond the resolution limit of satellite images. During UAV flight missions, geological surveying and sampling were also carried out to combine with the results of remote sensing and UAV mapping. The complementary field prospecting and laboratory studies were performed after preparing the geological map draft. Finally, integrating all information layers, we provided a geological map with indicators for chromite prospecting.
3.2. Satellite Remote Sensing
3.2.1. Sensor Data and Characteristics
Several satellite images from various platforms have been applied, including Aster, Landsat series, and Sentinel-2. Landsat missions have been continuing since 1974, which, as the longest continuous space-based record, provide the opportunity for managing and monitoring Earth’s resources.
The three cloud-free scenes were downloaded from different Landsat mission sensors. These images were taken in 1987 (Thematic mapper, TM), 1999 (Enhanced thematic mapper, ETM+), and 2021 (Landsat8, LC8). The 14-band ASTER Level 1B data acquired on 8 May 2001, was used in this research. The extended 11 bands in the short-wave infrared region (SWIR) and thermal-infrared (TIR) regions enable mineral mapping and surface lithological discrimination.
The Sentinel-2 MultiSpectral Instrument (MSI) at the moment provides the highest resolution among the free satellite images with 13 spectral bands ranging from visible-near infrared (VNIR) to SWIR wavelengths, spatial resolution from 10 to 20 m, and three bands with 60 m spatial resolution for using in the atmospheric correction. The cloud-free image of Sentinel-2B utilized in this study is a Level-1C top-of-atmosphere (TOA) reflectance (100 km × 100 km tile) product which was acquired on 18 September 2021. The sensors’ characteristics and the wavelength regions of all the instruments are given in
Figure 3. All satellite images were obtained from the U.S. Geological Survey Earth explorer website (
https://earthexplorer.usgs.gov/ (accessed on 30 May 2022).
3.2.2. Preprocessing
All satellite images used in this study were georeferenced to the UTM zone 40 North projection using the WGS84 datum. For processing the remote sensing datasets, various software packages, including ENVI version 5.3 (L3Harris Geospatial, Broomfield, CO, USA) and ArcGIS version 10.8 (Esri, Redlands, CA, USA), and PCI Geomatica (version 2015, PCI Geomatics, Markham, ON, Canada), were used.
The reflectance products for Landsat images were downloaded from the mentioned site, and they did not need any further correction. Sentinel-2 data includes radiometric correction along with orthorectification. A layer of VNIR + SWIR bands (bands 2, 3, 4, 5, 6, 7, 8, 8a, 11, and 12) with 10 and 20 m spatial dimensions was generated to obtain a ten-band dataset. After atmospheric correction, all the image bands were resampled to 10 m spatial resolution and layer stacked into one file. Bands 1, 9, and 10 were excluded since their resolution (60 m) and design are beyond the main objective of this study.
The ASTER imagery was supplied in terms of scaled radiance at-sensor data with geometric corrections applied. It was georeferenced in the UTM projection and for the WGS-84 ellipsoid. ASTER data layer stacked of VNIR + SWIR bands with 30 m spatial dimensions was generated. Internal Average Relative Reflectance (IARR) calibration was applied to Crosstalk-corrected ASTER data for atmospheric correction.
3.2.3. Processing Techniques
The conventional techniques, such as some band combinations and band ratio color composites recommended for targeting potential chromite zones [
8,
9,
37], were employed to compare with previously published papers regarding chromite targeting [
10,
13]. Landsat images were used to detect the temporal changes that occurred over nearly five decades (
Appendix A).
Target detection for finding spectral similarities of the whole scene of ASTER and TM with training data (
Appendix A) was performed using SAM and CEM methods. Before that, the MNF Rotation transformations were used to separate noise from the data and reduce the intrinsic dimensionality of data [
38]. The capability of satellite images for lithological mapping and detection of hydrothermally altered areas has been widely accepted, and some methods have been adopted for discriminating lithological units of ophiolites [
13,
16,
39,
40,
41,
42] which were applied as a starting point for targeting chromite hosts.
3.3. UAV Mapping
The main study area, extended 3500 Ha (35 km
2), was selected for detailed topography and geological mapping, according to the exploration license. UAV photogrammetry has been divided into two main phases: (i) fieldwork, including establishing Ground Control Points (GCP) with GNSS (Global navigation satellite system) solutions, planning flight missions, and collecting photos, and (ii) processing, described in the following sections and
Appendix A. Based on the required resolution for distinguishing the dykes as the finest geological units in the area [
43,
44,
45], the ground sampling distance (GSD) was selected to be 6 cm (
Appendix A).
3.3.1. Fieldwork and Image Acquisition
The 70 GCPs were considered to ensure the accuracy requirements and to georectify the dense point cloud. A total of 14 of GCPs were selected to be as check points to double-check the overall accuracy of the DEM derived from the dense point cloud of the whole area. The center of markers was measured using the ground-based Real-Time Kinematic-differential (RTK)-based Integrated GNSS Systems (
Figure 4).
The used UAV was equipped with a GNSS module which determines the coordinates of the center of each taken photo with an accuracy of ±5 cm. After processing GNSS module measurements with a ground-based station using the PPK (post-processing kinematic) method, the images can be geo-tagged with a specialized program that gives an accuracy of better than 5 cm. Therefore, both direct georeferencing (using the on-board GNSS module with in situ ground support) and indirect georeferencing (using GCPs) were applied in this study to satisfy the georeferencing task [
46]. This combined method of georeferencing yields higher accuracy than the other methods.
To obtain the highest possible quality for images, the UAV speed, image sensitivity (ISO), time of flight (sunny days), and other camera parameters with adjusting the manual focus to aperture and shutter were taken into account during final flight planning in the field as described in [
47]. The survey design was optimized in advance to minimize motion blur and to guarantee a 6 cm resolution orthomosaic (Vollgger and Cruden, 2016) [
47].
3.3.2. Processing
After finishing the aerial survey and UAV flight missions, Agisoft metashape professional software (v1.8.4, Agisoft LLC, Saint Petersburg, Russia) was applied for photogrammetric processing. This software uses the popular structure from motion (SfM) procedure [
48] to obtain 3-D reconstruction from photogrammetric imagery. The details of the processing stages are represented in
Appendix A.
3.4. Field and Geological Mapping
Geological surveying, sampling, microscopic studies on thin and polish sections, and whole-rock chemical analysis, including XRF (PW1730, Philips, The Netherlands) and ICP-OES (four-acid digestion, Spectro Arcos, SPECTRO, Kleve, Germany), were carried out to confirm the remote sensing results and outlined lithological units from photogrammetry imagery. A survey was conducted in the study area to verify the spatial distribution of serpentinized areas and other hydrothermally altered and lithological units using both GNSS solution (accuracy of cm) and handheld GPS (average accuracy of 3 m). Moreover, during flight missions, we took the samples and marked the lithological units with estimation for their name and other geological structures. Additionally, numerous photos were taken from alteration zones and lithological units during the field surveys, which were then geo-tagged within a GIS environment to ease the comparison with the UAV-mapping-derived orthomosaic. Thirty hand specimens were collected from highly serpentinized zones, which contained disseminated chromite, alteration zones, ore (Cu, Mn) mineralization, and lithological units for laboratory analysis.
3.5. Machine Learning
Recently, using machine learning methods in geoscience studies has been increasing in various applications, including mapping land cover types and monitoring land usage [
49]. Some work has been carried out in the field of identifying and classifying lithologies of economically valuable minerals [
50,
51]. The Support Vector Machine algorithm (SVM) is one of the widely used machine learning methods for geological applications [
52], which, over the years, has proven to be robust and effective [
53], especially in helping in the identification of rock units (Bentahar and Raji 2021; Shirmard et al., 2022b [
14,
50]).
The SVM method was used to classify lithological units in the study area, particularly those belonging to the ophiolitic sequence, and to discriminate peridotites with variable ranges of serpentinization (
Appendix A). This was carried out after preparing the initial geological map based on the manual interpretation of the orthomosaic obtained from UAV image processing and field checking. The SVM method was applied to ASTER data for classifying lithological units of the Mawat ophiolite complex in northeastern Iraq to help target chromite deposits (Othman and Gloaguen, 2014) [
12]. However, they discriminated ultramafics as a single lithological unit.
4. Results
4.1. Satellite Remote Sensing
The first step of remote sensing started with change detection, by which we can identify and discriminate the exposed areas subject to chromite extraction. Using the change detection, algorithms yielded the delineation of mining areas where the most exploitation activities have been performed in the region. We found that the most prominent activity during the last 35 years occurred on the west side of the surveying area (
Figure 5). The field investigation showed that this area was under mining at the time and was called the Sarvar mine. Additionally, there are some other spots in the ophiolite section which are most probably related to chromite extraction.
Before spectral mapping, the conventional band ratio combinations were also successfully used to highlight the different parts of ophiolitic sections. The ASTER band ratios composite (4/7, 4/1, 2/3 × 4/3;
Figure 6) has been successfully applied to study serpentinized harzburgites [
8,
54]. The application of this band ratio shows the chromite-bearing area of the Sarvar mine as red-pinkish to purple colors, which are related to the serpentinized harzburgites; however, the contact between the serpentinized harzburgites with other ophiolitic rocks is not well discriminated. The ASTER band ratios composite of ((2 + 4)/3, (5 + 7)/6, (7 + 9)/8) image was proposed as the highlighter of the serpentinized harzburgites between other ophiolitic rocks [
9]. However, in the study area, this combination did not work well. Alternatively, we used the band combination of ((7 + 9)/8, 4/3, 5), which, by linear enhancing, provides an image in which the highly serpentinized areas, including the main quarry workings of the chromite mine are discriminated as orange to yellowish-orange color (
Figure 6).
We used the reflectance spectra of the main chromite quarry as the reference spectra at the first stage of target detection using SAM and CEM methods. Two boundaries for the region of interest (ROI) were considered, one of which envelops only the internal parts of the main open-pit. By applying this ROI as target spectra together using extremely conservative thresholds, few pixels were found with maximum similarity. These pixels most probably contain chromite mineralization (
Figure 7). In the second boundary, the wider area of the mine was drawn to define the potential lithology, although we know that the highly serpentinized peridotites and dunites host the chromite mineralization, it is important to find the spectral signature of these rocks in the study area considering possible spectral mixing. The results of target detection with this broader area show the widespread distribution of potential lithological hosts in the study area. The reflectance spectra of both ROI share similarities with the serpentinite, which has a relatively strong absorption feature at 2.3 μm due to the vibrational processes of Mg–OH [
8].
Chromite has distinct absorption features at 0.55, 0.59, 0.69, 0.94, 1.75, and 2.3 μm. The best correlations exist between Cr content and wavelength positions at 0.49, 0.59, and 2 μm in the SWIR absorption bands. Most chromitite bodies, however, have dimensions smaller than the resolution limit of satellite images. Nevertheless, the spectral mapping using standard USGS spectral libraries and the SAM method was carried out for this mineral together with serpentine, talk, antigorite, and lizardite. The results then were compared with the results of target detection using the training classes (chromite mines). We used the “intersect” tool in ArcGIS software to combine the results of target detection layers obtained from different sensors and various methods, which yields the “target overlap” map in
Figure 7. The results of target detection algorithms (SAM and CEM), using reflectance spectra of minerals as a reference, is highly consistent with the results of the same methods by using reflectance spectra of the region of interest (i.e., the quarry of the neighbor chromite mine). This is particularly true for serpentine minerals of the standard spectral library.
4.2. Photogrammetry and Photogeology
The root means square error (RMSE) was measured to test the accuracy of measurements for the ground control points, check points, and camera location (
Table 1). The total RMSE of 8 cm was obtained for fourteen checkpoints that satisfied the required accuracy of less than 10 cm for the orthomosaic image.
The general 3D model and DSM, which were derived from the dense point cloud, are presented in the
supplementary materials Figure S1. Additionally, some magnified parts of the resulting orthomosaic are presented in
Figure 8. The image provides valuable insight into the lithology of the study area. While the orthomosaic has a Red Green Blue (RGB) composite that does not allow conventional spectral mapping, its ultra-resolution provides an image with a rich texture that can be used as the basic data for classification and drawing the geological units. The ultra-spatial resolution of the orthomosaic provided an observation experience with the naked eye at about 0.5 m distances, combined with the ground truth data in the field, achieved satisfactory recognition accuracy. The characteristic texture of orthomosaic, together with the outcrop of rocks with minimum vegetation and debris cover, allows for delineating geological boundaries. Therefore, the primary polygons of rock units were drawn using ArcGIS software. On rare occasions, virtually in situ debris hindered the precise mapping of specific contacts, which then were verified by field prospecting. Moreover, the color differences among the mafic-ultramafic units were too low for a clear distinction. Peridotites, in some places, could be identified through their dark appearance, which their color grade from black–dark green to light green–gray due to serpentinization (
Figure 8). However, the amount of serpentinization is not apparent only based on the orthomosaic image. The serpentinite unit, in some cases, can be distinguished from the less altered peridotites by their brighter color, smoother surface, and fine texture (
Figure 8). In some places, the serpentinites were partly covered by a thin layer of soil having dull-yellow to orange color which is formed during the weathering of Ni-Fe-bearing minerals and iddingitization of olivine (
Figure 8B).
The ultramafic rocks were crosscut by several dykes, which in most cases appeared as fragments due to the effect of faulting or might be just isotropic lenses (
Figure 8). Almost all intruded dykes, regardless of their composition and thickness were detectable and outlined by digitizing through orthomosaic because of their presumed relationship to chromite mineralization. Due to their higher resistance to weathering, these dykes stand higher than the surrounding rocks. Thin dykes with thicknesses smaller than 1 m are common in all ultramafic units, which sometimes appear as bright white, probably because of the alteration or metasomatism. The dykes, according to occurrence and abundance, can be divided into individual dykes (ca. 1–10 m thick) and dyke swarms (thickness up to 1 m), the latter are distributed within highly serpentinized peridotites (
Figure 8A). Due to this character, we separated the mixed lithology of swarm dykes distributed within ultramafics as a distinctive unit which was then partly verified in the field. We assigned its name as MTZ, or transition zone unit of the mantle, since this unit, overlaying the highly serpentinized peridotites, marks a transition zone in lithology. Intrusive bodies as thick bands and stocks, injected into mantle peridotites with discordant dunite lenses, are common in the east-northeast direction of the area (
Figure 8C). This unit was outlined as a distinctive unit called Dcu due to the abundance of intrusions into peridotites with patches of serpentinized dunites (orange to dull-yellow surface color).
A lithological unit in the northeastern corner of the mapped area, which was also the highest altitude part of the study area, had distinctive texture and color, in some parts covered by the stain of Fe-oxides, which made it distinguishable from the neighboring units. According to the regional geological map, this unit is a colored mélange, which was confirmed by field evidence. Apart from this area, two other units can be specified based on the color differences, both in UAV and in satellite images, resemblance the post-ophiolitic subvolcanic rocks. Moreover, all old workings, including small trenches which were for exploring chromite, are detectable in the orthomosaic; hence they were accurately delineated and shown on the geological map.
4.3. Field Work and Laboratory Studies
Geological surveying was carried out in two stages, first during UAV mapping, based on previous satellite remote sensing results, and second after image processing, producing an orthomosaic image and drafting the initial geological map with the help of the SVM method (for ranking serpentinized peridotites).
Field campaigns verified the initial lithological units obtained from the orthomosaic and identified the main lithological units, including mantle sequence rocks, which led to higher precision in outlining the lithological boundaries. Almost all of the lithological units in the study area belonged to the mantle section of the ophiolitic sequence. However, the abundance of gabbro-dioritic dykes and the occurrence of a massive, heavily eroded gabbro-dioritic body points to grading toward the crustal section of the ophiolitic domain. The ultramafic rocks vary from harzburgite to clinopyroxene-bearing harzburgite, and minor lherzolite, with variably sized lenses and dykes of dunite, suffered from variable degrees of serpentinization. Harzburgites consisted of olivine, large orthopyroxene porphyroblasts, and spinel with rare clinopyroxene (
Figure 9A,B). While the lherzolites are less serpentinized (~10%), dunites are strongly serpentinized (~90%), and only a few relicts of olivine and spinel survive. The range of harzburgite serpentinization varied from 20 to 80%. However, some areas exist in which the harzburgites and dunites were completely serpentinized. Lherzolites are most common in higher altitudes and have gradational contacts with harzburgites with similar mineralogy and higher content of clinopyroxene.
Microscopic investigation of thin and polished sections highlighted the main lithological units and their degree of serpentinization and mineralization. Together with the ICP-OES and XRF chemical analysis (
supplementary material Table S1), it can be concluded that there are two types of chromite mineralization in the area, high-Cr, low-Al, and low-Cr, high-Al; however, the Al contents rise up to 15% in the high-Al category [
55]. The concentration of chromite range from 5% in the form of dissemination into harzburgite to 50% in nodules. The texture of chromites varied from disseminated, leopard, and shliern to nodular. Chromitites were hosted by a thin envelope of highly serpentinized dunite to serpentinized harzburgite without dunite traces. These chromite types can be seen in a short distance and have a spatial relationship with the occurrence (high Al type) or absence of mafic dykes. The occurrence of magnesite and hydromagnesite minerals around some chromite pods reflects the impact of faulting and circulation of hydrothermal fluids and probably the local redistribution of minerals, including chromite. Apart from chromite, other mineralization types, including Fe, Cu, and Mn, were also observed hosted by serpentinized peridotites and mafic dykes which highlights the importance of accurate discrimination of these units. In the mineralized areas, the effects of leaching and hydrothermal alteration are almost visible locally.
During the microscopic study, caution has been taken to discriminate peridotites according to their degree of serpentinization. With this in mind, four types of peridotite were recognized, including weakly serpentinized harzburgite to lherzolite (mainly harzburgite), moderately serpentinized, highly serpentinized, and serpentinite (
Figure 9). Due to the higher content of olivine in dunite and clinopyroxene-poor harzburgite, we presume that lithological units with the highest serpentinite have dunitic or clinopyroxene–poor harzburgitic protolith because of more vulnerability of olivine and orthopyroxene to serpentinization in comparison with clinopyroxene. Additionally, serpentinites show two surficial appearances and colors, which might be related to their weathering pattern and protolith. These patterns are also partly observable in orthomosaic images (
Figure 8). The first pattern of serpentinites is represented by the gray to dull-white color and gradual contact with moderate to highly serpentinized harzburgites (
Figure 8D and
Figure 9E). The protolith of this serpentinite type was probably clinopyroxene-poor harzburgite. The second type, typically seen as lenses, dykes, and individual bodies distributed within serpentinite type1, had an orange to dull-yellow surface color (
Figure 8B), probably due to the release of Iron from the ferromagnesian minerals and iddingitization of olivine. The thin bright color crust was slightly enriched in Ni, according to chemical analysis. This serpentinite (marked as serpentinite2 in the SVM classification scheme,
Figure 10 and
Figure 11) likely had dunitic protolith and should be considered the most promising unit hosting chromite mineralization.
The dykes which crosscut the ultramafic units (including both depleted harzburgites/dunites) based on their mineralogy vary from diabase, gabbro to norite, and wehrlite to pyroxenite. However, more felsic lithological names can also be assigned to the sampled dykes, including hornblende micro-diorite to hornblende micromonzodiorite, some of which are post-ophiolitic intrusions. The highly serpentinized harzburgite to serpentinite, within which dyke swarms with various thicknesses intruded, resemble the cumulate section of the mantle transition zone. Additionally, other areas with dyke and banded mafic intrusions exist in which the serpentinite lenses and dykes (probably concordant dunitic protolith) occur with a prominent iron oxide (limonitic) crust. These intrusions range from crustal cumulate gabbro/gabbronorite/ferrogabbro (with minor diorite) with local layering associated with minor ultramafic cumulates. The microdioritic dykes intruded into cumulate gabbros which show clear evidence of magma mingling. Small (4–5 mm, rarely 2–3 cm) clinopyroxenite-wehrlitic veins are common within all ultramafic units of mantle sequence, including both types of serpentinites (
Figure 8).
A subvolcanic body with microdiorite mineralogy was also observed, which regarding the alteration halo, was considered a post-ophiolitic magmatic unit. The colored mélange, with its characteristic colors and iron oxide stains on the outer surface of rocks, can be easily recognized both in RGB orthomosaic and in hand samples. Moreover, a volcanic unit characterized by spilitic lavas and dykes is seen in the eastern part of the area, which, together with small blocks of pelagic limestone, represents a part of the extrusive section of the ophiolite sequence. The boundary of limestone is distinctive in the orthomosaic image; however, the volcanic unit does not show an apparent boundary regarding its texture and color and is verified by field evidence.
The occurrence of depleted harzburgite, dunite, wehrlitic dykes/sills, and chromitite marks a Moho Transition Zone in the region, as suggested by previous works [
31,
56]. Tectonically overlying the transition zone (marked by thrust fault) is a melange zone containing lenses of cumulate rock, highly metamorphosed ultramafics (up to greenschist facies), pillow lavas, and blocks of pelagic limestone. These rocks occur as blocks in a highly serpentinized, ultramafic matrix. Individual lithologies are irregularly distributed, and all the rocks are highly deformed.
4.4. SVM
The second stage of satellite remote sensing based on the spectral properties in the VNIR to SWIR region of the electromagnetic spectrum was carried out after drafting the initial geological map. This aimed to classify rock units based on the spectral properties of the representative lithology. Because of rugged mountain topography with high altitude variation (up to 850 m) at short distances and steep valleys, it is impossible to conduct field surveying with specified narrow traverses as recommended for geological mapping and prospecting at a large scale. On the other hand, from the orthomosaic alone, the distinction among the peridotites with different ranges of serpentinization based on their color and texture was not possible. Therefore, by combining the results of drawing geological units from the orthomosaic, field surveying, sampling, and also the first results of satellite remote sensing, the 12 areas, representative of 12 presumed lithological units, were separated and traced to use as training data in SVM method, of which five areas are related to peridotites and serpentinized derivative products, as described in the previous sections. The difference among these five units is related to the intensity of reflectance spectra which itself is a function of the severity of serpentinization (
Figure 10). To select the training data related to peridotite, we only traced areas with at least one thin section from a representative sample of the area, to ensure the high confidence level. One additionally used training data which also may belong to serpentinized peridotite classes is the location of an open pit of chromite mining (the ROI used for SAM method).
All the presumed lithological units were successfully reproduced for the entire study area. However, individual dykes are usually thinner than the spatial resolution limits of satellite images. Additionally, the training data in the area of dykes swarm enclosed in serpentinized peridotites (Dcu unit) display reflectance spectra that are affected by reflectance spectra of both serpentinite and gabbroic dykes (
Figure 8C). The classified lithological map after post-processing (
Figure 11) is similar to the initially digitized lithological map of the study area, while the latter has more details. Therefore, the lithological boundaries were retraced to be in accordance with the new results of the SVM method. The peridotites without the effect of serpentinization are rare in the study area, so the few separated pixels with these characters were discriminated after the classification of the ASTER image. The moderately serpentinized peridotites and peridotites with the slightest affection by serpentinization may have clinopyroxene-bearing harzburgite to lherzolite protolith.
The rationale behind the successful classification is revealed in the average reflectance spectra of each training area (
Figure 10). The weak to moderately serpentinized peridotites have average reflectance spectra with absorption features around 0.65 μm to 0.7 μm (correspondent with band2 of ASTER), attributable to ferrous Iron in pyroxene and olivine minerals. All of the serpentinized peridotite and serpentinite areas have absorptions around 0.8 (3rd band position of ASTER), which might be the effect of the presence of the ferrous Iron on the weathered surface of rocks or abundance of the olivine, hypersthene, enstatite, augite, and diopside minerals [
8]. However, the chromite open-pit together serpentinites did not show evident absorption around 0.65 μm to 0.7 μm, probably due to higher Mg content than Fe in serpentinized dunites or the content of chromite [
57]. The higher reflectance in band1 of ASTER than band3 of ASTER was also the result of higher content of Mg than Fe
+3 in ultramafic rocks.
The most striking resemblance among peridotite-related classes was a strong absorption feature around 2.3 μm due to the vibrational processes of Mg\OH, i.e., the occurrence of serpentine) [
8]. The absorption features around 2.35 μm are also responsible for C\Obonds, such as magnesite and Fe
+2 content in olivine [
8,
11] The abundance of serpentine (antigorite-lizardite) mineral, which also has two absorption features at 1.4 and 2.35 μm, defines the intensity of the absorption [
9]. Therefore, the serpentinites with higher content of serpentine mineral together with chromite open-pit showed higher intensity of absorption feature in the 2.35 μm (band 8 of ASTER). In addition, the spectra of serpentinites are closer than other rock types to the current chromite open-pit spectrum. In band4 of ASTER (around the at 1.6 to 1.7 μm) peridotite related classes displayed different behavior; the weakly serpentinized one showed an absorption feature while the serpentinite2 (equivalent to serpentinite grading to talcitization with Ni-iron rich stain) showed the reflectance. Two reasons for showing reflectance might be the occurrence of thin gabbroic dykes or a dunitic lense containing chromite [
12] which have high reflectance in both wavelengths of 1650 nm and 2150 nm (bands 4 and 5 of ASTER reflectance data, respectively). Since both serpentinite2 and cumulate area classes show reflectance in these wavelengths, the two assumed reasons should be verified in the field. The reason for showing absorption around 1.6 to 1.7 μm is due to the presence of hydroxyl molecules and the hydration effects of water molecules in the altered minerals or ferrous Iron. The spectral plot of colored mélange shows substantial similarity to the hematite spectral plot, which correspond with the iron-oxide alteration seen on the orthomosaic image and field evidence. Instead, the mélange2 area, while having a similar spectral plot with mélange, showed smoother absorption features due to the effect of spectral mixing.
All of the classes related to gabbro and subvolcanic bodies show higher reflectance than ultramafics, both in SWIR and visible bands, which also appeared in their brighter colors in orthomosaic image. The area characterized as a subvolcanic body did not have strong absorption features but shares spectral similarity with dioritic rocks, which have absorption at around 2.1 μm and 2.2 μm. The gabbro-dioritic area shows an absorption feature in band 2 ASTER and high reflectance at bands 4, 6, 8, and 9, similar to the spectral feature of amphibole-rich rock, as well as the presence of pyroxene minerals.
Two approaches were considered for testing the accuracy of classified maps, the first of which was by calculating the confusion matrix and Kappa Coefficient for the results. The overall accuracy for Sentinel-2 and ASTER datasets is 83.9% and 86.6%, respectively. The Kappa Coefficient was also assessed at 0.81 for Sentinel-2 and 0.84 for the ASTER dataset (
Table 2). Although the total pattern of discriminated units by the two satellite platforms is the same, the visual inspection showed that Sentinel-2 is not capable of distinguishing among the variously serpentinized units fully. This is because of the subtle changes in spectral absorption features within 1.9 to 2.3, while the Sentinel-2 image has only one band in this SWIR region (
Figure 3).
The second approach, which we believe is more accurate, is the comparison with the other results of this study, including field evidence, the orthomosaic, the initial geological map, and the second stage of field prospecting (see
Section 5). With this criterion, the ASTER image has more capability to discriminate mafic-ultramafic lithological units and the degree of serpentinization due to its higher number of bands in the SWIR region (
Figure 3). The Sentinel-2 combined two classes distinguished based on the training data for weakly and moderately serpentinized peridotite and classified them as a single unit. On the other hand, the Sentinel-2 accurately revealed some small size areas of serpentinite units, according to other information layers, which are not discriminated in ASTER image due to its higher spatial resolution.
Since the focus was on the ophiolite sequence, especially the mantle section, which is the dominant lithology in the study area, we did not try to define the training data for the classification of the volcano-sedimentary units younger than ophiolites including quaternary deposits and debris, which are weakly developed in the study area, but widespread at the outside of study area. So, there may be some misclassification due to spectral mixing in these areas, which are not the subject of this study. Moreover, we can precisely digitize such boundaries using the ultra-resolution orthomosaic.
5. Discussion
According to the regional geological map of the area, the ultramafic section of the ophiolitic complex has been divided into three main categories: harzburgite-lherzolite, dunite, and serpentinized counterparts. However, we note that the degree of serpentinization and the protolith of serpentinites are two main factors for prospecting chromite pods which were neglected in the regional scale maps, which is a routine procedure in these geological scale mapping. Some other mentioned lithological units do not have accurate boundaries. Moreover, the diversity of dykes in terms of composition, mineralogy, and pattern of distribution was not reflected in previous geological maps. An additional geological survey was conducted to finalize the geology map. The sampling locations and one of the promising areas were positioned on the updated geological map after lithological classification by SVM to facilitate the field checking and decrease the cost of the prospecting plan. In comparison to the regional geological map (
Figure 1), except for the existence of mafic and ultramafic sections of the ophiolite sequence, all other parameters have entirely changed, including the number, size, position, and distributions of lithological units (
Figure 12), which was mainly caused by the difference in cartographic scale.
Accompanying chromitite lenses by mafic (gabbro to pyroxenite) dykes, dunitic pods, and cumulate layers [
35,
36,
58,
59] has been reported, which might be related to the spatial preference of chromite mineralization. Most of the chromite pods were formed at or in the vicinity (within 1 km) of the petrological Moho discontinuity, or MTZ, where there are abundant olivine-rich (mainly dunitic) residual lenses, discontinuous bands or lenses of wherlite, troctolite, pyroxenite and gabbros, and overlying cumulate units [
4,
35]. The identification of MTZ is not an easy task, especially in highly tectonized dismembered ophiolites and due to its highly varied thickness (from a few meters to several hundred meters) [
60,
61]. However, we had this chance to trace the petrological moho, including transition zone unit (MTZ) with serpentinites, dunitic horizons, and abundant dyke swarms (
Figure 8A,
Figure 11 and
Figure 12); to cumulate unit (Dcu) with abundant cumulate dyke and banded gabbro and dunitic lenses (highly altered,
Figure 8C,
Figure 11 and
Figure 12); and to outline these units in the final geological map (
Figure 12), due to gathering valuable high spatial and radiometric resolution data. Outlining the thin mafic dykes and small bodies is a task that could be quickly carried out using ultra-resolution orthomosaic images and field evidence (
Figure 8), while they are absent in the low to medium-resolution satellite images investigated in this study. The microscopic studies can also be served to define the mineralogy and corresponding name of the dykes.
On the other hand, the main factor for prospecting chromite pods is the necessity to find dunite, troctolite, and more serpentinized harzburgite and dunites [
4,
8,
35,
62]. In the study area, unaltered dunites have not been seen except on rare occasions. One reason might be the increased vulnerability of olivine (as the main component of dunite) to serpentinization, compared with harzburgite and lherzolite. It is widely accepted that podiform chromites are products of melt–rock interaction in which the rocks are mantle peridotite (mainly depleted harzburgites), and melts are usually arc-related [
34,
58,
59,
63,
64]. This interaction leads to the creation of replacive dunites enveloping chromitites, especially in supra-subduction zone ophiolites [
59]. Additionally, due to this process, the host rock harzburgites show a gradual change in their composition, mineralogy, and texture as the modal abundance of pyroxene decreases near the dunite channels and the modal olivine increases [
55,
59,
63]. The cpx-poor, high-Cr harzburgites have irregular and continuous boundaries with the chromitite–dunite zone and the surrounded cpx-bearing, low-Cr harzburgites; however, such subtle distribution has been proposed to be used as a vector to locate the podiform chromitites [
65]. Another effect of this zoning is that the extent of serpentinization and iddingitaztion of olivines increase toward chromitites, which in turn results in contrasting yellow to dull yellow-white colors. These changes in color are marked partly by the orthomosaic image (
Figure 8), but the increasing degree of serpentinization is more evident in processed (
Figure 6 and
Figure 7) and classified multispectral images (
Figure 10 and
Figure 11).
The serpentinization process affects spectral absorption features and roughness of rocks since serpentine is usually readily weathered, and consequently, the surface of serpentinized rocks is smoother than their neighboring lithological units. The surface textures can be compared through visual interpretation of orthomosoaic while extracting the spectral features need multi and hyperspectral images. Different physical processes, including weathered surfaces, might produce spectral features in the VNIR spectral region. However, most mineral absorption features are located in the SWIR wavelengths [
8,
11,
57]. Combining the satellite images with those acquired from UAV photogrammetry may resolve the problem of classifying lithological units confined within the peridotite to serpentinite domain. Machine learning methods, such as SVM, could help to classify complex lithological units, including those seen in the study area, based on spectral similarities [
12]. Nevertheless, due to the spectral mixing effect and fuzziness of spectral features in low spatial resolution of satellite images [
15] we assigned more weight to visual inspection and interpretation to finalize the lithological mapping. Therefore, we prioritized our knowledge of the study area and field evidence in a few cases where orthomosaic and SVM results did not overlap. On rare occasions, the chromitites were seen in the crustal section of ophiolites [
59]. The crustal section in the surveying area includes a gabbro-dioritic body, an area covered by volcanic rocks (mainly spilitic lavas) and gabbroic dykes, among which the intrusive and volcanic bodies have higher reflectance in visible bands and different spectral absorption feature, which make them discriminable in both satellite images and orthomosaic images. In contrast, in the case of dykes, especially thin dyke swarms with thicknesses smaller than 1 m, the only option, except field mapping, is delineation from the orthomosaic. While some researchers have tested the capability of using machine learning methods for the automatic classification of lithological units based on photogrammetry-obtained orthomosaics [
21,
66], the extent of the area and diversity of lithological units matter. They used the ML method in small areas with a significant difference among the units, two features that were not the case in our study area.
Benefiting from the synergy of optical remote sensing and UAV imagery, the resulting detailed map and lithological boundaries were acquired in spending much shorter time and with higher accuracy compared to traditional field-based techniques. Moreover, we were able to identify subtle lithological/compositional differences that would not have been identified in the field, as the high-resolution orthomosaic allowed rapid viewing of the outcrop at multiple scales and from different perspectives and scales [
44]. Nevertheless, field-based geological mapping still was the backbone of the whole exploration procedure described in this study for verifying remote sensing data.
The tectonic structures have to be outlined for their effect on the spatial distribution of podiform chromite deposits [
67,
68]. In some cases, the newly arranged deposits are associated with shear zones formed in or around zones of weakness in preexisting chromitite bodies [
4]. These new chromite occurrences and original chromitite bodies are affected by displacement and deformation through faulting and tectonic processes. Therefore, finding tectonic colored mélange zones can help to avoid prospecting in areas where the size and position of chromite deposits have been substantially reduced and changed due to intense faulting and folding.
Integrating the lithological map classified through machine learning, visual interpretation, and field checking yielded a more detailed geological map from which the promising rock units for prospecting chromite were delineated. (
Figure 12). A combination of layers of this map with results of target detection and spectral mapping (
Figure 7) was used to find the most promising areas for chromite prospecting, one of which is highlighted in
Figure 11. Some areas were then verified by field checking, which revealed the validity of the methods after finding the chromite outcrops (
Figure 13) or overlapping with the old outcrops. The new outcrops are in a highly altered serpentinite with a thin covering soil rich in Iron and nickel, reminiscing the dunitic protolith. These new chromite-bearing outcrops linearly extend over 100 m at a narrow distance from a thick mafic dyke. The chromites outcrop as nodules or lumps within the host rock, and the chromium content is high (Cr
2O
3 ~51.78%,
supplementary material). Chromite textures vary from massive to more nodular, disseminated types and rarely appear as banded schlieren.
6. Conclusions
Detailed lithological mapping and chromite ore exploration in the mountainous study area with altitude variation reaching 850 m and narrow, steep gorges is costly and time-consuming using the traditional exploration methods. The differentiating ultramafics in ophiolitic assemblages in terms of their serpentinization degree, while crucial for finding chromite deposits, is not practical in the field or by traditional geological mapping. To our knowledge, all remote sensing studies on ophiolite domains have been concentrated on separating serpentinized peridotites as mantle section domains from other crustal sections. This is the first effort to look closer into the mantle section and transitional zone rock units in an ophiolitic sequence, differentiate them and discriminate ultramafics based on the degree of serpentinization, which led to classify these rocks in our study area. Moreover, the identification of the mantle transition zone, including the ultramafic cumulate unit and overlying cumulate gabbro, is a probable task using the above-mentioned techniques, which is an efficient way for prospecting dunite–chromite lenses which always concentrate in and around this zone. To aim for this, we took advantage of the ultra-spatial resolution of UAV images with poor spectral information and the high radiometric and spectral resolution of the satellite images with poor spatial resolution. Combining spectral mapping methods such as SAM and CEM methods, RGB orthomosaic interpretations obtained from UAV photogrammetry, and initial field surveying, we prepared the primary draft of the lithological map and the promising areas for prospecting chromite. Then, training data from this map were extracted for supervised classification using the machine-learning methods (SVM method). Applying the SAM method on ASTER data was more successful than other medium-resolution satellite images due to the spectral signature of different mafic-ultramafic rocks and their serpentinized equivalents in the SWIR region. This statement is also true for the SVM method because other platforms, such as Sentinel-2 and Landsat series, have only one band at the 1.9 to 2.3 interval of the electromagnetic spectrum.
The results show that integrating satellite remote sensing and UAV photogrammetry is an efficient way to identify different rock units, including five units of peridotite with variable ranges of serpentinization, together with two units ascribed to cumulate and transition zone layers. UAV-based photogrammetry generates models that are useful for multi-scale analysis not easily achieved using traditional field techniques. The subtle changes in texture, topography, and color of the distributed lithologic units on the Earth’s surface could help to discriminate dykes and other rock units, which is not achievable by satellite images. However, verification of remote sensing results is crucial following each stage by field checking, prospecting, and sampling.
The results of this study are promising in identifying the potential zones of chromite ore mineralization within the ophiolite region. There are many ancient and recent chromite mines in the region, but most are inactive due to the small deposit size. Applying this method could help to place the most potential areas for chromite mineralization and discover more outcrops which might help to locate clusters of chromite pods.