This study investigated the influences of urban growth and heat waves events on Urban Heat Island in relationship with several biophysical variables in Bucharest metropolitan area of Romania through satellite and in-situ monitoring data. Remote sensing data from Landsat TM/ETM+ and time series NOAA AVHRR and MODIS Terra/Aqua sensors have been used to assess urban land cover– temperature interactions over period between 2000 and 2016 years. Vegetation abundances and percent impervious surfaces were derived by means of linear spectral mixture model, and a method for effectively enhancing impervious surface has been developed to accurately examine the urban growth. The land surface temperature (LST), a key parameter for urban thermal characteristics analysis, was also analyzed in relation with the Normalized Difference Vegetation Index (NDVI) at city level. Based on these parameters, the urban growth, urban heat island effect (UHI) and the relationships of Ts to other biophysical parameters (surface albedo, precipitations, wind intensity and direction) have been analyzed. Results show that in the metropolitan area ratio of impervious surface in Bucharest increased significantly during investigated period, the intensity of urban heat island and heat wave events being most significant. The correlation analyses revealed that, at the pixel-scale, LST possessed a strong positive correlation with percent impervious surfaces and negative correlation with vegetation abundances at the regional scale, respectively. This analysis provided an integrated research scheme and the findings can be very useful for urban ecosystem modeling.