This thesis aims to investigate the spatial drivers of Surface Urban Heat Island Intensity (SUHII) in two European case studies: Wrocław, Poland and Turin, Italy. The objective is to map surface thermal anomalies and to check how far satellite-derived indicators, urban morphology variables and GIS-based modeling can explain their spatial distribution. The workflow was first developed for Wrocław and then applied to Turin following the same main procedure so that the two cities could be compared with consistent data processing steps. Landsat 8/9 Collection 2 Level 2 imagery was processed using Google Earth Engine to derive Land Surface Temperature (LST), SUHI, SUHII, UTFVI, and several surface indicators including NDVI, NDBI, NDWI, NDMI, Albedo, PVI and Emissivity.