The numerical weather prediction model HARMONIE-AROME and a multiple linear regression model (referred to in this article as the TURCLIM model after the local climate observation network) were used to model surface air temperature for 25–31 July 2018 in the City of Turku, Finland, to study their performance in urban areas and surrounding rural areas. The 0200 LT (local standard time) temperatures modeled by the HARMONIE-AROME and TURCLIM models were compared to each other and against the observed temperatures to find the model best suited for modeling the urban heat island effect and other spatial temperature variabilities during heatwaves. Observed temperatures were collected from 74 sites, representing both rural and urban environments. Both models were able to reproduce the spatial nighttime temperature variation. However, HARMONIE-AROME modeled temperatures were systematically warmer than the observed temperatures in stable conditions. Spatial differences between the models were mostly related to the physiographic characteristics: for the urban areas, HARMONIE-AROME modeled on average 1.4°C higher temperatures than the TURCLIM model, while for other land-cover types, the average difference was 0.51°C at maximum. The TURCLIM model performed well when the explanatory variables were able to incorporate enough information on the surrounding physiography. Respectively, systematic cold or warm bias occurred in the areas in which the thermophysically relevant physiography was lacking or was only partly captured by the model.