Global land cover dynamics alter energy, water, and greenhouse gas exchange between land and atmosphere, affecting local to global weather and climate change. Although reforestation can provide localized cooling, ongoing land use land cover (LULC) shifts are expected to exacerbate urban heat island impacts. In this study, we monitored spatiotemporal changes in green cover in response to land use transformation associated with the Khyber Pakhtunkhwa (KPK) provincial government’s Billion Tree Tsunami Project (BTTP) and the Ravi Urban Development Plan (RUDP) initiated by the provincial government of Punjab, both in Pakistan. The land change modeler (LCM) was used to assess the land cover changes and transformations between 2000 and 2020 across Punjab and KPK. Furthermore, a curve fit linear regression model (CFLRM) and sensitivity analysis were employed to analyze the impacts of land cover dynamics on land surface temperature (LST) and carbon emissions (CE). Results indicated a significant increase in green fraction of +5.35% under the BTTP, achieved by utilizing the bare land with an effective transition of 4375.87 km². However, across the Punjab province, an alarming reduction in green fraction cover by −1.77% and increase in artificial surfaces by +1.26% was noted. A significant decrease in mean monthly LST by −4.3 °C was noted in response to the BTTP policy, while an increase of 5.3 °C was observed associated with the RUDP. A substantial increase in LST by 0.17 °C was observed associated with transformation of vegetation to artificial surfaces. An effective decrease in LST by −0.21 °C was observed over the opposite transition. Furthermore, sensitivity analysis suggested that LST fluctuations are affecting the % of CO₂ emission. The current findings can assist policymakers in revisiting their policies to promote ecological conservation and sustainability in urban planning.

Spatial land surface heterogeneities are widespread at various scales and represent a great challenge to leaf area index (LAI) retrievals and product validations. In this paper, considering the mixed water and vegetation pixels prevalent at moderate and low resolutions, we propose a methodological framework for conducting global comparisons of heterogeneous land surfaces based on criterion setting and a global search of high-resolution data. We construct a global network, Heterogeneous Surface Network aiming Water and Vegetation Mixture (HESNet-WV), comprised of three vegetation types: grassland, evergreen broadleaf forests (EBFs), and evergreen needle forests (ENFs). Validation is performed using the MCD15A3H Global 500-m/4-day and GLASS 500-m/8-day LAI products. As the water area fraction (WAF), LAI values and LAI inversion errors increase in the MODIS and GLASS products, the GLASS product errors (relative LAI error (RELAI): 94.43%, bias: 0.858) are lower than the MODIS product errors (RELAI: 124.05%, bias: 1.209). The result indicates that the proposed framework can be applied to evaluate the accuracy of LAI values in mixed water-vegetation pixels in different global LAI products.