•Quantify the interaction between urban ecological land and environmental factors.•Urban ecological land show spatiotemporal coupling with surface urban heat island.•The average explanatory power of patch dominance and richness was 19.95% and 16.03%.•Interaction of urban ecological land with topography rapidly increased in 2015–2020.•Interaction between urban ecological land and anthropogenic factors was dominant. The surface urban heat island effect (SUHI) that occurs during rapid urbanization increases the health risks associated with high temperatures. Urban ecological land (UEL) has been shown to play an important role in improving urban heat stress, however, the impact of UEL interactions with the natural-anthropogenic environment on SUHI at the urban agglomeration-scale is less explored. In this study, the Google Earth Engine and GeoDetector were applied to characterize the spatiotemporal patterns of UEL and SUHI in the Guangdong-Hong Kong-Macao Greater Bay Area from 2000 to 2020 by extracting major built-up urban areas and quantifying the impacts of UEL and its interactions with the natural-anthropogenic factors on SUHI. The results show that the evolution of the UEL landscape structure exhibits clear spatiotemporal coupling with SUHI. Specifically, the UEL underwent a dispersion and degradation process in 2000–2015 and a convergence and restoration process in 2015–2020, the SUHI correspondingly transitioned from intensification and continuity to mitigation and contraction. The UEL landscape structure showed a notable impact on the SUHI reduction, and the dominance and richness of the patches explained an average of 19.95% and 16.03% of the SUHI, respectively. Moreover, the interaction between UEL and land urbanization rate and anthropogenic heat release had a dominant effect on SUHI, but this effect significantly declined from 2015 to 2020. With the implementation of ecological restoration projects, the interaction of UEL with topography rapidly increased and the SUHI gradually dominated by the joint interaction of UEL and natural-anthropogenic factors. A synthesis of the varying effects of several factors showed that the dynamic relationship between the development stages of the urban agglomeration’s regional system and SUHI may conform to the Environmental Kuznets Curve. SUHI reduction strategies should therefore comprehensively optimize the rational allocation of UEL landscape structures and natural-human elements to promote the well-being of residents.