Land surface temperature (LST) in urban areas is an important environmental variable considered a reliable indicator of the urban heat island (UHI) phenomenon. LST is affected by various factors such as solar irradiance, cloudiness, wind or urban morphology. Traditionally, LST is observed and recorded by thermal remote sensors. For example, thermal satellite sensors are very popular for assessing the UHI effect on a global scale such as MODIS, Sentinel 3, ASTER, Landsat 7 ETM+, or Landsat 8 TIRS. However, these sensors provide rather low spatial (60 m to 1000 m) and temporal resolutions (several hours to days) of satellite observations that limit the accurate estimation of LST in urban areas for local studies and specific time periods (Mushore et al., 2017), (Hu and Wendel, 2019). Airborne or terrestrial remote sensing can be viewed as another option for capturing higher spatial resolution of thermal data but it is not feasible to be used for large urban areas with increased periodicity. However, the increasing availability of the high-resolution geospatial data and adequate modeling techniques provide an alternative approach to high-resolution estimation of LST in urban areas. Several studies showed the potential of geographic information system (GIS) tools, digital surface models (DSM) and 3-D city models for the estimation of solar radiation in urban areas (e.g., Hofierka and Kaňuk, 2009; Hofierka and Zlocha, 2012; Freitas et al., 2015; Biljecki et al., 2015). Solar irradiance is a key factor affecting LST during daylight periods, especially under clear sky situations. Nevertheless, LST assessment requires a physical model combining surface-atmosphere interactions and energy fluxes between the atmosphere and the ground. Properties of urban materials, in particular, solar reflectance, thermal emissivity, and heat capacity influence the LST and subsequently the development of UHI, as they determine how the Sun’s radiation energy is reflected, emitted, and absorbed (Hofierka et al., 2020b; Kolečanský et al., 2021). It is clear, that the problem complexity requires a comprehensive GIS-based approach. Our solution is based on open-source solar radiation tools available in GRASS GIS, a 3D city modeling and spatially distributed data representing thermal properties of urban surfaces and meteorological conditions (Hofierka et al., 2020a, 2020b; Kolečanský et al., 2021) . The proposed LST model is calculated using the methodology implemented in GRASS GIS as a LST module written using a script (shellscripts, Python). In these scripts, the r.sun and v.sun solar radiation models in GRASS GIS were used to calculate the effective solar irradiance for selected time horizons during the day . The solar irradiance calculation accounts for attenuation of beam solar irradiance by clouds estimated by field measurements. The proposed LST model also accounts for a heat storage in urban structures depending on their thermal properties and geometric configuration. The 2D LST model uses the output of the r.sun solar radiation model and a DSM representing urban surfaces and the 3D LST model uses the output of the v.sun solar radiation model and a vector-based 3D city model. Computed LST values for selected urban surfaces were validated using field measurements of LST in 10 locations within the study area with acceptable accuracy. The proposed approach has the advantage of providing high spatial detail coupled with the flexibility of GIS to evaluate various geometrical and land surface properties for any daytime horizon. The methodology can be used for evaluation of proposed UHI mitigation measures such as increasing albedo of urban surfaces or expanding green areas including green roofs and trees.