GIS, geochemical and geophysical based identification of geothermal energy resources in sedimentary basins: A case study of Zindapir Anticlinorium zone, Eastern Sulaiman Fold and Thrust Belt, Pakistan

Aqeel Rasheed¹Ali Wahid^1*Hammad Tariq Janjuhah^2*Fahad Hameed¹Syed Haroon Ali³Malik Talha Riaz⁴George Kontakiotis⁵Assimina Antonarakou⁵Evangelia Besiou⁵Ahmed Shafique⁶

• ¹Institute of Geology, University of Azad Jammu and Kashmir, 13100, Muzaffarabad, Pakistan
• ²INTI International University, Nilai,, Nagri Sembilan, Malaysia
• ³Department of Earth Sciences, University of Sargodha, 40100,, Sargodha, Pakistan
• ⁴University of Silesia, Institute of Earth Science, Natural Sciences,, Silesia, Poland
• ⁵Department of Historical Geology and Paleontology, Faculty of Geology and Geoenvironment, National and Kapodistrian University of Athens, Panepistimiopolis, 15784 , Greece;, Athens, Greece
• ⁶Pakistan Institute of Nuclear Science and Technology, Islamabad, Pakistan

Globally, fold and thrust belts are considered promising zones for geothermal energy exploration due to their complex structural geometries. These zones provide favorable conditions that contribute to the trapping and circulation of geothermal resources. However, their inherent structural complexity poses significant challenges in delineating geothermal reservoirs. This research focuses on analyzing the geothermal reservoir potential in the Zindapir Anticlinorium Zone, which is part of the Eastern Sulaiman fold and thrust belt in Pakistan. The objectives of this research were achieved by utilizing an integrated approach involving remote sensing, geochemical, gravity, magnetic, seismic, and petrophysical analyses. The remote sensing data indicate that regional and local tectonic activity, along with the uplift of Eocene and Paleocene rocks, play a significant role in the development of geothermal reservoirs. Additionally, higher land surface temperature anomalies, active seismicity, and high-gradient streams passing through the hot springs are key indicators of geothermal resource potential. The higher ionic concentration of Na+ in geochemical analysis indicates deep circulation of water that affects thermal activities. Temperature estimations suggest that the Gulki hot spring lies in a zone characterized by immature thermal water, whereas the hot springs of Jaffarabad and Zindapir originate from subsurface zones with partially equilibrated water. Interpretations of gravity and magnetic data suggest that the anticlines of the Zindapir Anticlinorium Zone have less sedimentary thickness and favorable heat flow conditions. Moreover, seismic and petrophysical analyses indicate that the selected reservoir formations (Pab sandstone, Lower Goru Formation and Chiltan Limestone) have potential to act as geothermal reservoirs. The Comprehensive Evaluation Index (CEI), which incorporates temperature, porosity, and permeability parameters, suggests that the Dhodhak, Affiband, Rodho, and Zindapir anticlines are ideal locations for the extraction of economic geothermal reservoirs in the Zindapir Anticlinorium Zone. These anticlines have less sedimentary thickness and compressional tectonic deformation, which allows upward transfer of heat from deeper sources. It is suggested that the workflow of this integrated approach can be implemented in any sedimentary basin to identify the potential and economic locations of geothermal reservoirs.