Forest cover change alters land–atmosphere exchanges of energy, water, and carbon, thereby influencing local and regional climate. This study assessed climatic patterns associated with afforestation in Khyber Pakhtunkhwa, Pakistan, from 2003 to 2023 using remote sensing data and Google Earth Engine. Land surface temperature (LST) was treated as the primary response variable, while evapotranspiration (ET) was analyzed as a secondary response variable. Air temperature; precipitation; vegetation indices, including the normalized difference vegetation index (NDVI) and the enhanced vegetation index (EVI); and elevation were used as supporting variables to interpret the broader climatic and biophysical responses of afforestation. MODIS land-cover, LST, ET, and vegetation-index products, together with climate research unit (CRU) climate data and ALOS-PALSAR DEM, were used to evaluate spatiotemporal trends and variable relationships. The results showed that mean LST increased by 0.520 ± 0.070 °C across KP during 2003–2023; however, areas classified as forest gain showed a localized cooling pattern of 0.490 ± 0.050 °C during the 2013–2023 forest-cover transition assessment window. Afforested areas also exhibited increased ET, whereas forest-loss areas showed reduced ET and higher LST. Specifically, ET increased by 0.013 ± 0.002 mm/8-day in afforested areas, whereas forest-loss areas showed a decline of 0.005 ± 0.001 mm/8-day. CRU-derived regional air temperature showed an increasing tendency of 0.310 ± 0.050 °C, whereas precipitation showed only a weak and statistically non-significant regional tendency; therefore, precipitation was used only as background climatic context. The NDVI and the EVI were negatively correlated with daytime LST, and elevation showed a strong negative relationship with LST. Overall, the findings indicate that forest-cover gain was associated with localized surface cooling patterns and improved vegetation–climate regulation indicators in the study area.
Khan et al. (2026) studied this question.