“…To understand the complex impacts of urban environment on the spatiotemporal variability of canopy temperature, we incorporated the major known environmental factors in the GAMs model analysis, namely, greenspaces (green area fraction, GAF), bluespaces (distance to waterbodies, DTW), and surrounding building morphological characteristics (building height, BH). GAMs are a flexible approach allowing the quantification of non-linear controls of confounding factors on spatially and diurnally 31 , 44 . Due to the common coexistence of various greenery, e.g., trees, shrubs, and grassland, in the urban landscape and their similar diurnal temperature patterns (Fig.…”