In recent years, the Quetta Valley and surrounding areas have experienced unprecedented levels of subsidence, which has been attributed mainly to groundwater withdrawal. However, this region is also tectonically active and is home to several regional strike-slip faults, including the north–south striking left-lateral Chaman Fault System. Several large earthquakes have occurred recently in this area, including one deadly Mw 6.4 earthquake that struck on 28 October 2008. This study integrated Interferometric Synthetic Aperture Radar (InSAR) results with GPS, gravity, seismic reflection profiles, and earthquake centroid-moment-tensor (CMT) data to identify the impact of tectonic and anthropogenic processes on subsidence and earthquake patterns in this region. To detect and map the spatial-temporal features of the processes that led to the surface deformation, this study used two Synthetic Aperture Radar (SAR) time series, i.e., 15 Phased Array L-band Synthetic Aperture Radar (PALSAR) images acquired by an Advanced Land Observing Satellite (ALOS) from 2006–2011 and 40 Environmental Satellite (ENVISAT) Advanced Synthetic Aperture Radar (ASAR) images spanning 2003–2010. A Small Baseline Subset (SBAS) technique was used to investigate surface deformation. Five seismic lines totaling ~60 km, acquired in 2003, were used to map the blind thrust faults beneath a Quaternary alluvium layer. The median filtered SBAS-InSAR average velocity profile supports groundwater withdrawal as the dominant source of subsidence, with some contribution from tectonic subsidence in the Quetta Valley. Results of SBAS-InSAR multi-temporal analysis provide a better explanation for the pre-, co-, and post-seismic displacement pattern caused by the 2008 earthquake swarms across two strike-slip faults.

Spatial mapping of potential geothermal areas is an effective tool for a preliminary investigation and the development of a clean and renewable energy source around the globe. Specific locations within the Earth’s crust display some manifestations of sub-surface geothermal occurrences, such as hot springs, a volcanic plug, mud volcanoes, and hydrothermal alterations, that need to be investigated further. The present area of investigations also reveals some of these manifestations. However, no attempt was made to examine the prospectivity of this terrain using the efficient GIS-based multicriteria evaluation (MCE) within the scope of the Analytic hierarchy process (AHP). The integration of remote sensing, Geographic information system (GIS), and other geophysical methods (Magnetic and gravity) was performed to map the promising geothermal areas. Multiple input data sets such as aero-magnetic, aero-gravity, aero-radiometric, digital elevation model (DEM), geological map, and Landsat-8 Operational Land Imager (OLI) data were selected, processed, and use to generate five thematic layers, which include heat flow, temperature gradients, integrated lineaments, residual gravity, and lithology maps. The five thematic layers were standardized and synthesized into a geothermal prospectivity map. The respective ranks and weight of the thematic layers and their classes were assigned based on expert opinion and knowledge of the local geology. This research aims to apply an efficient method to evaluate the factors influencing the geothermal energy prospects, identify and map prospective geothermal regions, and, finally, create a geothermal prospectivity model.

The Suleiman Fold-Thrust Belt represents an active deformational front at the western margin of the Indian plate and has been a locus of major earthquakes. This study focuses on the western part of the Suleiman Fold-Thrust Belt that comprises two parallel NW–SE oriented faults: Harnai Fault and Karahi Fault. These faults have known thrust components; however, there remains uncertainty about the lateral component of motion. This work presents the new observation of surface deformation using the Small Baseline Subset (SBAS), Interferometric Synthetic Aperture Radar (InSAR) technique on Sentinel-1A datasets to decompose displacement into the vertical and horizontal components employing ascending and descending track geometries. The subsurface structural geometry of this area was assessed using 2D seismic and well data. In addition, geomorphic indices were calculated to assess the relative tectonic activity of the area. InSAR results show that the Karahi Fault has a ~15 mm right-lateral movement for descending and ~10 mm/for ascending path geometries. The Harnai Fault does not show any lateral movement. Seismic data are in agreement with the InSAR results suggesting that the Harnai Fault is a blind thrust. This work indicates that the block between these two faults displays a clockwise rotation that creates the “bookshelf model”.