The prediction of energy demand from buildings represents a key information for optimizing the energy supply, and for planning the use of renewable energy sources, as well as the adoption of strategies to mitigate undesired urbanization effects, such as higher temperatures associated with urban heat islands (Oke et al., 2017). Indeed, in 2018, buildings accounted for 26% of the total final energy use in the EU, and for about 16% of energy-related greenhouse gas emissions. In particular, 64% of the final energy consumption in the residential sector was due to space heating, with renewables accounting for about a quarter (European Commission, 2020). In this regard, several works highlighted that significant reductions in carbon dioxide (CO 2 ) emissions can be achieved from the optimization of the energy use (Pérez-Lombard et al., 2008). Accordingly, many modeling techniques have been developed to simulate energy consumption for space heating and cooling, both at the building and at the city scale. However, as reported in Swan and Ugursal ( 2009), these models may differ as to the level of detail required for the input data, for the time required for processing, and for the degree of complexity. Models such as the US Department of Energy's Energy-Plus (https://energyplus.net/) or TRNSYS (www.trnsys.com) operate at the building scale (Pappaccogli et al., 2018) and have been extensively applied to assess the profiles of energy demand of specific building types (e.g., Kazas et al., 2015;Moreci et al., 2016). Indeed, such models can reproduce quasi-real building features, such as homogeneous zones inside the building interiors, complex layered walls, glazing and shading systems, as well as solar and occupant passive heat gains. Nevertheless, they generally consider buildings as stand-alone entities, neglecting the interactions with the surrounding neighboring ones. This may result in significant errors in the simulation of the building energy demand, as shown by Han et al. (2017). For example, Pisello et al. ( 2012) report overestimates of the energy demand for space heating up to +22%, if the thermal interactions between surrounding neighboring buildings are neglected. In fact, energy consumption for both space heating and cooling is strongly influenced by the surrounding envi-Abstract This work aims at assessing the ability of the Weather Research and Forecasting model, coupled with a multilayer urban canopy scheme and a building energy model (BEP + BEM), to simulate the building energy consumption in the city of Bolzano (Italy) during wintertime. Estimates of the actual energy consumption, to be adopted as a benchmark for evaluating model results, are obtained from observed total annual values, combined with load curves derived from real-time measurements in various meteorological conditions and for different building types. The validation of numerical results highlights that the modeling set-up adopted is able to well reproduce the spatial distribution of the energy consumption in the urban area and its dependence ...