The increase in built surfaces constitutes the main reason for the formation of the Urban Heat Island (UHI), that is a metropolitan area significantly warmer than its surrounding rural areas. The urban heat islands and other urban-induced climate feedbacks may amplify heat stress and urban flooding under climate change and therefore to predict them correctly has become essential. Currently in the COSMO model, cities are represented by natural land surfaces with an increased surface roughness length and a reduced vegetation cover, but this approach is unable to correctly reproduce the UHI effect. By increasing the model resolution, a representation of the main physical processes that characterize the urban local meteorology should be addressed, in order to better forecast temperature, moisture and precipitation in urban environments. Within the COSMO Consortium a bulk parameterization scheme (TERRA_URB or TU) has been developed. It parametrizes the effects of buildings, streets and other man-made impervious surfaces on energy, moist and momentum exchanges between the surface and atmosphere, and additionally accounts for the anthropogenic heat flux as a heat source from the surface to the atmosphere. TU implements an impervious water-storage parameterization, and the Semi-empirical Urban canopy parametrization (SURY) that translates 3D urban canopy into bulk parameters. This paper presents evaluation results of the TU scheme in high-resolution simulations with a recent COSMO model version for selected European cities, namely Turin, Naples and Moscow. The key conclusion of the work is that the TU scheme in the COSMO model reasonably reproduces UHI effect and improves air temperature forecasts for all the investigated urban areas, despite each city has very different morphological characteristics. Our results highlight potential benefits of a new turbulence scheme and the representation of skin-layer temperature (for vegetation) in the model performance. Our model framework provides perspectives for enhancing urban climate modelling, although further investigations in improving model parametrizations, calibration and the use of more realistic urban canopy parameters are needed.