The remaining 18 studies were performed in different districts of Pakistan. The study by Abbasi et al.
[56] in the forest divisions of the districts of Sukkur and Shikarpur (1997 km
2 area), Sindh, used Landsat MSS and TM satellite images to generate landcover maps of riverine forests for 1979, 1992, 1998, 2000, 2006, and 2009. Their study examined a considerable decline, 85%, of riverine forests from 1979–2009. Forest cover in 1979 was observed as 22.67% of the total study area, while in 2009, it was 5.97% (119.2 km
2). The total forest cover of both districts was calculated as 12 km
2 using the GFW data set (2019). The study by Iqbal and Khan
[57] in the Muzaffarabad district of AJK (740 km
2 area) used Landsat TM satellite data to map subtropical, evergreen, dry, broad-leaved forests, subtropical Chir pine forests, and temperate broad-leaved forests for 1998 and 2009 using the pixel-based supervised image classification method. The classification results revealed that, from 1998 to 2009, over a period of about 11 years, forest cover and low vegetation decreased at the annual rate of 2.70% and 2.60%, respectively. This study reported a forest cover of 155.5 km
2 for 2009 compared to 724 km
2 calculated through the GFW data set (2019). Batool et al.
[58] used Landsat ETM data to map forest cover in the Thak valley (Diamir District), KP. Using supervised image classification, the forest change assessment in the study area (213 km
2) was performed for 1989, 1999, and 2009, where forest decreased from 85.83 km
2 in 1989 to 34.4 km
2 in 2009. The total forest area in the entire district was calculated as 570 km
2 using the GFW data set for 2019. Among all the reviewed papers in this review, only one study, conducted by Baig et al.
[59], used active remotely sensed data. Using a calibration technique, they mapped an irrigated forest plantation (47 km
2 area) in the Sahiwal district, Punjab, through SAR ALOS-2 PALSAR (the Phased Array type L-band Synthetic Aperture Radar). Along with this, they used WorldView-3 imagery to verify their results. They mapped various forest tree species, including Shisham (
Dalbergia sissoo), Sufeda (
Eucalyptus camaldulensis), Toot or Mulberry (
Morus alba), and Simal (
Bombax ceiba) in pure and mixed forms with naturally grown Mesquite (
Prosopis juliflora). The total forest area in the entire district was calculated as 3 km
2 using the GFW data set for 2019. Younis and Ammar
[60] mapped the forest area in Besham watershed (which lies within Buner and Mansehra districts), KP, using Landsat 5 TM satellite data from 2000–2010 with an area of 6812 km
2. This study reported a net loss of 17% (3.88 km
2) from 2000–2010 in the forest area. This study reported forest cover of 1892.6 km
2 for 2010 compared to 1849 km
2 forest cover in both of these districts (Buner and Mansehra) through the GFW data set for 2019. The following study, led by Rashid and Iqbal
[61], used Landsat 4–5 TM, Landsat 7 ETM, and Landsat 8 OLI to map conifer and pine forests along the Karakoram Highway (N-35), through KP, to GB (4,200 km
2 area). The study used imageries of 1990, 2000, 2010, and 2016 and reported a 26% decrease in the forest cover during this period. The total forest cover of KP and GB combined was 11,042 km
2, as per the GFW data set for 2019. The following study, by Khan et al.
[62], used Landsat 5, Landsat 7, and Landsat 8 images for 1996, 2003, and 2016 to map the urban forest in Peshawar district (1257 km
2 area), Pakistan. This study used the supervised image classification method to report vegetation (including forest) increase from 1996 to 2003 and then decrease from 2003 to 2016. An overall vegetation decrease of 21.90% was described between 2003 to 2016. The total forest area in the entire district was calculated as 9 km
2 using the GFW data set for 2019. The study by Urooj and Ahmad
[63] incorporated Landsat 5 TM, Landsat 7 ETM+, and Landsat 8 OLI to map the dry, subtropical thistle and scrub forests in and around the surroundings of Mangla Dam (3053 km
2 area) in the Mirpur district of AJK. Using the supervised image classification method, they mapped the land cover of 1992, 2002, and 2013 and reported that deforestation occurred with a net change of 547.45 km
2 from 1992 to 2013. The total forest area in the entire district was calculated as 63 km
2 using the GFW data set for 2019. The study by Khan et al.
[64] used Landsat TM, Landsat ETM+, and Landsat OLI satellite data sets to map temperate forests in Sudhnuti district (471 km
2 area), AJK. They extracted vegetation fractions (forest, non-forest area) using Linear Spectral Mixture Analysis (LSMA), a supervised image classification approach, for 1989, 1993, 1999, 2005, 2010, 2015, and 2018. As a result, they reported that the forest area declined between 1989–1993 and 1993–1999; it increased between 1999–2005, 2005–2010, 2010–2015, and 2015–2018. This study reported approximately 203 km
2 of forest cover for 2018 compared to 297 km
2 for 2019, as calculated through the GFW data set. The subsequent study by Mahmoudi et al.
[65], used Moderate Resolution Imaging Spectroradiometer (MODIS) data to map regional mixed forests of Balochistan and the eastern Iran region (769,824 km
2 area). Using MODIS land cover type product (MCD12Q1) on a yearly basis, this study reported that the forest area increased between 2001–2013. The total forest area in the entire province of Balochistan was calculated as 18 km
2 using the GFW data set for 2019. The following research, by Ali et al.
[66], used SPOT-5 satellite imagery to map Subalpine, dry temperate, moist temperate, oak, subtropical broad-leaved, subtropical pine, and dry tropical thorn forests of KP (11,336 km
2 area) using OBIA. The total forest area in the entire province of KP was calculated as 10,123 km
2 using the GFW data set for 2019. The following study by Saddique et al.
[67] used Landsat TM, Landsat ETM+, and Landsat OLI to map the evergreen and deciduous forests in the river Jhelum basin (Mangal Dam watershed), AJK (33,397 km
2 area). Using a machine-learning algorithm (RF) in supervised image classification, they reported that the forest area showed a positive difference of 2806.87 km
2 between 2001–2009 and 2009–2018. This study reported a forest cover of 12,118 km
2 for the entire watershed for 2018. The total forest area in the entire AJK region was calculated as 4061 km
2 using the GFW data set for 2019. The following study, conducted by Hussain et al.
[68], used Landsat TM and Landsat OLI satellite imagery in the urban forest (3650 km
2 area) of Multan district, Punjab. They mapped the study area for two cropping seasons of
Rabi and
Kharif. They reported that the forest area gained a negative change of –41 km
2 from 1988–2017 during the
Rabi season and –48.6 km
2 for the same period during the
Kharif season. This study reported a forest cover area of 53 km
2 compared to 59 km
2 extracted through the GFW data set for 2019. Another research work by Hussain et al.
[69] for mapping urban forest in Lodhran district, Punjab, used Landsat 4–5 TM, Landsat 7 ETM+, and Landsat 8 OLI satellite images for 1977, 1987, 1997, 2007, and 2017. Using a maximum likelihood classifier in the supervised image classification, they reported that the vegetation (forest) area increased from 87.9% (of the total study area) in 1987 to 90.8% in 2017. The total forest area in the entire district was calculated as only 1 km
2 using the GFW data set for 2019. That the subsequent study used the supervised image classification method was published by Khan et al.
[70], who used Landsat 7 and Landsat 8 satellite images to map the Himalayan moist temperate and sub-alpine temperate forests’ AGB of Battagram (1507 km
2 area), KP, for 2000 and 2015. Their study indicated a –16.88% loss in the forests of the study area with an annual deforestation rate of 2.51%. This study reported a forest cover area of 450.8 km
2 in 2015 compared to 462 km
2 extracted through the GFW data set for 2019. Ali and Nayyar
[71] used Landsat 8 OLI satellite imagery to assess mangrove forests located in the Karachi region’s coastal belt (2030 km
2 area) for 2017 using the unsupervised image classification method. They used pixel-based spectral indexes, including NDVI, Normalized Difference Moisture Index (NDMI), Ratio Vegetation Index (RVI), EVI, Combined Mangrove Recognition Index (CMRI), and Soil Adjusted Vegetation Index (SAVI) to delineate mangrove and non-mangrove land covers. This study reported a mangrove forest cover area of 228.6 km
2 in 2017 compared to 34 km
2 extracted through the GMW data set for 2019. The difference in the areas is because mangroves do not have an exact boundary. The paper by Haq et al.
[72] used Landsat 3 MSS, Landsat 7 ETM+, and Sentinel 2A satellite images to map dense Deodar and Pine forest,
Juniperus Communis (juniper),
Pinus Wallichiana (blue pine),
Abies Webbiana (silver fir),
Aesculus Indica (bankhor), Pinus species,
Cedrus Deodara (deodar),
Abies Pindrow (palunder),
Pinus Gerardiana (chalghoza), and
Juglans Regia (walnut) in the Palas valley in Kohistan district (7492 km
2 area), KP. Using the maximum likelihood classification technique in the supervised image classification method, this study reported that the forest cover declined between 1980–2017, with a change of −12.23%. This study reported a forest cover area of 198.6 km
2 in 2017. The forest cover of the entire Kohistan district, extracted through the GFW data set, was 1304 km
2. The last reviewed study, performed by Zafar et al.
[73], used MODIS data to map the coniferous, broad-leaved, and mixed forests in GB (68,601 km
2 area) using a pixel-based supervised image classification method. This study reported that a substantial increase in the forest was observed between 2008–2017 yearly. This study reported an overall forest cover in the study area as 0.081% of the total area, while the GFW data set showed 919 km
2 of forest cover in the entire GB for 2019.