This article aims to aid in exploring and forecasting hydrothermal minerals and gold deposits in Central Saudi Arabia (SA), with a focus on structural contexts. Remote sensing (RS) and satellite-based Bouguer gravity (SBG) data were integrated in order to create a mineral prediction map for the researched location. Data from the Landsat Operational Land Imager (OLI) and Shuttle Radar Topography Mission (SRTM) were transformed and enhanced using a variety of approaches. The delineation of hydrothermal alteration zones (HAZs) and highlighting of structural discontinuities in the OLI data were made possible using band ratios and oriented principal component analysis (PCA). Additionally, the underlying structural features were successfully exposed by processing the SBG using a variety of edge detection techniques, like the analytical signal (AS), total horizontal derivative (THD), tilt angle (TA), horizontal tilt angle (TDX), theta map (TM), horizontal derivative of the tilt derivative (HD_TDR), horizontal gradient of the tilt angle (HGTA), tilt angle of the analytical signal (TAAS), and soft sign function (SF). As a result, more prominent lineaments were found in the NW–SE, NNW–SSE, NE–SW, and NNE–SSW directions than in the N–S and E–W directions. The GIS incorporated surface/subsurface geological structure density maps with zones of hydrothermal alteration. It was found that the lineaments derived from the analysis of the RS and SBG data were more in line with the HAZs, which demonstrated the common connection between alteration zones and deep lineaments. The findings revealed a mineral prediction map with extremely low to extremely high probabilities. Overall, combining RS and SBG data effectively identified probable mineralization sites associated with hydrothermal processes and made it easier to create this study’s final predictive mineralization map.

The present research’s main objective was to apply thorough exploration approaches that combine remote sensing data with geochemical sampling and analysis to predict and identify potential chromitite locations in a complex geological site, particularly in rugged mountainous terrain, and differentiate the ultramafic massif containing chromitite orebodies from other lithologies. The ultramafic massif forming the mantle section of the Kırdağ ophiolite, located within the Erzurum–Kars Ophiolite Zone and emerging in the east of Oltu district (Erzurum, NE Turkey), was selected as the study area. Optimum index factor (OIF), false-color composite (FCC), decorrelation stretch (DS), band rationing (BR), minimum noise fraction (MNF), and principal and independent component analyses (PCA-ICA) were performed to differentiate the lithological features and identify the chromitite host formations. The petrography, mineral chemistry, and whole-rock geochemical properties of the harzburgites, which are the host rocks of chromitites in the research area, were evaluated to verify and confirm the remote sensing results. In addition, detailed petrographic properties of the pyroxenite and chromitite samples are presented. The results support the existence of potential chromitite formations in the mantle section of the Kırdağ ophiolite. Our remote sensing results also demonstrate the successful detection of the spectral anomalies of this ultramafic massif. The mineral and whole-rock geochemical features provide clear evidence of petrological processes, such as partial melting and melt–peridotite interactions during the harzburgite formation. The chromian spinels’ Cr#, Mg#, Fe³⁺, Al₂O₃, and TiO₂ concentrations indicate that the harzburgite formed in a fore-arc environment. The Al₂O₃ content and Mg# of the pyroxenes and the whole-rock Al₂O₃/MgO ratio and V contents of the harzburgite are also compatible with these processes. Consequently, the combined approaches demonstrated clear advantages over conventional chromitite exploration techniques, decreasing the overall costs and supporting the occurrence of chromite production at the site.