Linking Dynamic Building Simulation with Long-Term Energy System Planning to Improve Buildings Urban Energy Planning Strategies

The building sector is currently responsible of 40% of global final energy consumption, influencing the broader energy system in terms of new electricity and heat capacity additions, as well as distribution infrastructure reinforcement. Current building energy efficiency potential is largely untapped, especially at the local level where retrofit interventions are typically enforced, neglecting their potential synergies with the entire energy system. To improve the understanding of these potential interactions, this paper proposes a methodology that links dynamic building simulation and energy planning tools at the urban scale. At first, a detailed bottom-up analysis was conducted to estimate the current and post-retrofit energy demand of the building stock. The stock analysis is further linked to a broader energy system simulation model to understand the impact of building renovation on the whole urban energy system in terms of cost, greenhouse gas emission, and primary energy consumption up to 2050. The methodology is suited to analyze the relationship between building energy demand reduction potential and clean energy sources’ deployment to shift buildings away from fossil fuels, the key priority for decarbonizing buildings. The methodology was applied to the case study city of Torino, Italy, highlighting the critical role of coupling proper building retrofit intervention with district-level heat generation strategies, such as modern district heating able to exploit low-grade heat. Being able to simulate both demand and supply future alternatives, the methodology provides a robust reference for municipalities and energy suppliers aiming at promoting efficient energy policies and targeted investments.


Introduction
In the last decade, several actions were deployed to contrast climate change. According to the IEA [1], addressing climate change it is necessary to limit global warming well below 2 • C beyond preindustrial levels in 2050. The European Union, in 2009, decided to plan a reduction of greenhouse gas (GHG) emissions to at least 80% below 1990 levels by 2050. To achieve this goal, the European Climatic Foundation redacted a study analyzing the feasibility of this target [2]. With the new European Green Deal [3], ambitions to reduce CO 2 emission are growing, setting the reduction target for 2030 to 55% and aiming to achieve a net zero emission in 2050. For a future decarbonized energy system, the first step will be improved energy efficiency, especially in building patrimony. Energy consumption of the building sector is 40% of total European energy use, accounting for 36% of associated energy related CO 2 emission [4]. In some areas, space heating accounts for a larger quota of building energy consumption, up to 80% in colder climates [5]. The low energy performance of existing long life management of energy sectors, due to its intuitive interface and simple data requirement, this tool can be used for urban application, with it being more understandable to urban stakeholders compared to other energy planning tools (e.g., IMEAS project [23]).
A complete overview of the overall energy system allows municipalities to handle future urban energy planning policies, developing a central role as mediator between different stakeholders from different sectors, being able to lead the construction of a common pathway for achieving future decarbonization goals.
At regional level, this tool was applied in Hong Kong [24], while at the city level, it was applied by [8], underlining that the tool considers only aggregate heat demand and also aggregate supply, with it not being detailed enough to describe alone an urban energy system. To meet the need of a better demand description in urban areas for reflecting the real status of the building stock and to integrate demand-supply analyses, this paper's goals were threefold: (1) Develop a spatial-based methodology to represent the urban building demand; (2) Link the results of the urban stock analysis with an energy planning tool to analyze both demand and supply decarbonization measures and their economic impact; (3) Apply the integrated tool to develop scenarios to support integrated energy planning strategies of an urban municipality.
While urban building stock studies focused on building renovation scenarios are common in literature [25,26], the implications and synergies of demand measures on the supply side are often neglected in urban analyses [27]. The main innovation of the methodology described in this paper is indeed related to its integrated structure, targeting urban analyses and providing an in depth analysis of the energy consumption of entire building stocks and potential energy refurbishment measures, including an additional analysis about the implication that a change in the demand side could have on the supply sector over the long term.
The integrated analysis was performed over a medium-long term period to help urban energy planners in defining technically feasible and economically affordable future energy efficiency scenarios in terms of CO 2 emissions, primary energy supply, and total annual cost. The integrated demand-supply analysis proposed was applied to a case study, the city of Turin. Future scenarios to achieve energy savings goal include the expansion of the existing district heating network of the city and the introduction of renewable energy technologies in the supply. This paper is structured in four sections: the first one analyzes the case study, the second one explains the methodology, and the third one presents and analyses the results. The applicability of the methodology, the necessary conditions, and its limits are further discussed in the last section.

Case Study
The methodology was applied to the case study city of Turin. Turin is a city situated in northern Italy, with 2706 heating degree days (HDDs) at 18 • C. According to previous literature, heating volume of residential building is almost 139 Mm 3 [28]. Almost half of the heating volume is satisfied by individual natural gas boilers, while the remaining volume is connected to district heating. The district heating (DH) in the city is the largest in Italy, with 550 km of pipes line and a building volume of 60.3 Mm 3 connected to it. The heat production is satisfied for 98% by three different cogenerating thermal plants administrated by IREN Energia [29], with a total of 1200 MW electrical power installed and 740 MW thermal power. The remaining heat need is satisfied by natural gas boilers, with a thermal power installed of 1000 MW. Three thermal storages are distributed in the city, with a total capacity of 430 MW. The expansion of district heating allows for a reduction in energy consumption and the integration of renewable energy sources, without a relevant impact in the total cost, as explained in the introduction. For this reason, the following study was concentrated on the heating volume connected to DH. Real data of heat consumption of district heating were provided by IREN, manager of district heating, for year 2014, with a total annual heat consumption of 1.99 TWh/year (in 2014 district heating supplied heat to 57 Mm 3 ). To evaluate a possible network expansion, another 19.5 Mm 3 were analyzed considering previous literature, as explained in the Section 3.1. The total volume considered in the case study was equal to 76.5 Mm 3 , of which 57 Mm 3 were attributed to the DH (as in the real case) and 19.5 Mm 3 were considered in the base scenario as individual heating, while in future scenarios they will be connected to DH. The assumptions of current and future expansion of the district heating are supported by literature [27,28].

Methodology
The presented methodology was composed of three key phases. The first one consisted of analyzing the energy consumption for space and water heating of identified residential reference building typologies and the potential for energy retrofit interventions, through dynamic simulation (Section 3.1). As a second step, a spatial analysis based on census data allowed researchers to quantify the distribution of reference buildings within the urban area and consequently to estimate the total energy requirement of the residential stock (Section 3.2). In the last step, the whole energy system was analyzed with an energy system simulation tool [30], allowing researchers to extend the analysis to the supply side and to develop multiple scenarios ( Figure 1) (Sections 3.3-3.5).
Smart Cities 2020, 3 FOR PEER REVIEW 4 DH (as in the real case) and 19.5 Mm 3 were considered in the base scenario as individual heating, while in future scenarios they will be connected to DH. The assumptions of current and future expansion of the district heating are supported by literature [27], [28].

Methodology
The presented methodology was composed of three key phases. The first one consisted of analyzing the energy consumption for space and water heating of identified residential reference building typologies and the potential for energy retrofit interventions, through dynamic simulation (Section 3.1). As a second step, a spatial analysis based on census data allowed researchers to quantify the distribution of reference buildings within the urban area and consequently to estimate the total energy requirement of the residential stock (Section 3.2). In the last step, the whole energy system was analyzed with an energy system simulation tool [30], allowing researchers to extend the analysis to the supply side and to develop multiple scenarios ( Figure 1) (Sections 3.3-3.5). During the first phase, the objective was to understand where buildings' energy-saving interventions are a priority (high-energy/high-density areas) and which level of the retrofit intervention can be reasonably pursued. Achieving this type of detail required the use of a dynamic simulation software to analyze the impact of different energy saving measures on heat consumptions. To analyze the built environment, the authors introduced the widely adopted concept of a "reference building" (RB). RBs are ''buildings characterized by and representative of their functionality and geographic location, including in-door and outdoor climate conditions" [6]. Simulating a RB requires a lot of data, especially when referring to dynamic simulation [31], usually difficult to find at an urban scale. The international project Typology Approach for Building Stock Energy Assessment (TABULA), 2009-2012, aimed to define European building typologies [32], providing a robust dataset of reference buildings in European countries, including the ones diffused in North Italy. The data from the TABULA project represented a valid support for this analysis.
The spatial distribution of buildings within the urban area was determined thanks to available census data of the city. This allowed us both to estimate the current consumptions of the building stock, but also their variations under different energy retrofit options. The impact of thermal demand variation on district heat generation plants, as well as new individual solutions, can be further studied During the first phase, the objective was to understand where buildings' energy-saving interventions are a priority (high-energy/high-density areas) and which level of the retrofit intervention can be reasonably pursued. Achieving this type of detail required the use of a dynamic simulation software to analyze the impact of different energy saving measures on heat consumptions. To analyze the built environment, the authors introduced the widely adopted concept of a "reference building" (RB). RBs are "buildings characterized by and representative of their functionality and geographic location, including in-door and outdoor climate conditions" [6]. Simulating a RB requires a lot of data, especially when referring to dynamic simulation [31], usually difficult to find at an urban scale. The international project Typology Approach for Building Stock Energy Assessment (TABULA), 2009-2012, aimed to define European building typologies [32], providing a robust dataset of reference buildings in European countries, including the ones diffused in North Italy. The data from the TABULA project represented a valid support for this analysis.
The spatial distribution of buildings within the urban area was determined thanks to available census data of the city. This allowed us both to estimate the current consumptions of the building stock, but also their variations under different energy retrofit options. The impact of thermal demand variation on district heat generation plants, as well as new individual solutions, can be further studied thanks to the energy planning tool EnergyPLAN [30]. This allows users to analyze how to coordinate the introduction of new heat supply technologies with building energy retrofit, both from a technical and economic perspective. The presented methodology was shaped to represent a resource for urban planners to analyze in depth buildings energy savings and their effects on the urban energy system.

Dynamic Building Simulation and Total Heat Consumption
Turin has 36,158 residential and occupied buildings according to census data, which are subdivided into eight construction classes and four building types ( Figure 2). Smart Cities 2020, 3 FOR PEER REVIEW 5 thanks to the energy planning tool EnergyPLAN [30]. This allows users to analyze how to coordinate the introduction of new heat supply technologies with building energy retrofit, both from a technical and economic perspective. The presented methodology was shaped to represent a resource for urban planners to analyze in depth buildings energy savings and their effects on the urban energy system.

Dynamic Building Simulation and Total Heat Consumption
Turin has 36,158 residential and occupied buildings according to census data, which are subdivided into eight construction classes and four building types ( Figure 2). In Turin, almost 96% of buildings were built before 1980, 40% of them are apartment blocks (ABs) and 28% are multifamily houses (MF) (Figure 3). District heating is connected to buildings with a heated volume higher than 2500 m 3 [33], so mostly to AB and MF. For this reason, in this study the chosen reference buildings were apartment blocks belonging to the three most frequent construction periods: ABs built from 1945 to 1960, ABs from 1961 to 1975, and ABs from 1976 to 1990. For each census section, census data provided information about the number of buildings, in turn separated into four building typologies (AB: apartment block; MF: multi-family; TH: terraced house; SF: single family) and in three age classes (pre-1980, post-1980, post-2005). Also, information about gross heated volume and net surface area were connected to these data. In Turin, almost 96% of buildings were built before 1980, 40% of them are apartment blocks (ABs) and 28% are multifamily houses (MF) (Figure 3). District heating is connected to buildings with a heated volume higher than 2500 m 3 [33], so mostly to AB and MF. For this reason, in this study the chosen reference buildings were apartment blocks belonging to the three most frequent construction periods: ABs built from 1945 to 1960, ABs from 1961 to 1975, and ABs from 1976 to 1990. Smart Cities 2020, 3 FOR PEER REVIEW 5 thanks to the energy planning tool EnergyPLAN [30]. This allows users to analyze how to coordinate the introduction of new heat supply technologies with building energy retrofit, both from a technical and economic perspective. The presented methodology was shaped to represent a resource for urban planners to analyze in depth buildings energy savings and their effects on the urban energy system.

Dynamic Building Simulation and Total Heat Consumption
Turin has 36,158 residential and occupied buildings according to census data, which are subdivided into eight construction classes and four building types ( Figure 2). In Turin, almost 96% of buildings were built before 1980, 40% of them are apartment blocks (ABs) and 28% are multifamily houses (MF) (Figure 3). District heating is connected to buildings with a heated volume higher than 2500 m 3 [33], so mostly to AB and MF. For this reason, in this study the chosen reference buildings were apartment blocks belonging to the three most frequent construction periods: ABs built from 1945 to 1960, ABs from 1961 to 1975, and ABs from 1976 to 1990. For each census section, census data provided information about the number of buildings, in turn separated into four building typologies (AB: apartment block; MF: multi-family; TH: terraced house; SF: single family) and in three age classes (pre-1980, post-1980, post-2005). Also, information about gross heated volume and net surface area were connected to these data. For each census section, census data provided information about the number of buildings, in turn separated into four building typologies (AB: apartment block; MF: multi-family; TH: terraced house; SF: single family) and in three age classes (pre-1980, post-1980, post-2005). Also, information about gross heated volume and net surface area were connected to these data.
With the support of GEOPORTALE of Turin [34] it was possible to consult the map with census sections of the city; in this way it was possible to connect territorial distribution of buildings with information about building typologies and construction periods. Considering, for each construction period, a uniform distribution in the city, it was possible to further subdivide census data, and for each section, the number of buildings comparable to the three reference buildings was calculated. Volumes connected to DH were identified, and 30.7 Mm 3 of them could be represented with RBs, in particular: 14.1 Mm 3 belonged to AB 1945-1960, 10.5 Mm 3 to AB 1961-1975, and 6.1 Mm 3 to AB 1976-1990. This subdivision demonstrates the correct choice of reference buildings, because 54% of volumes connected to DH were associable with analyzed buildings.
To evaluate a possible network expansion in relation with energy savings measures, another 43 Mm 3 were analyzed, considering previous literature and assuming that zones near the river were considered not connectable [33]. Of these, 19.5 Mm 3 connectable were represented by RBs: 9 Mm 3 to AB 1945-1960, 6.5 Mm 3 to AB 1961-1975, and 4 Mm 3 to 1975-1990 (Table 1). These 19.5 Mm 3 , considered as individual heating buildings in the base scenario, will be gradually connected to DH in future scenarios. The TABULA project supported the association of thermal characteristics to every RB. The project presented a matrix with 32 different types of building [35,36], and a schedule for each of them containing the heated volume, the floor area, the number of apartments, the number of floors, the surface to volume ratio (S/V), the construction materials with total U-values, energy systems with their energy efficiency, and heat consumption result from steady-state calculation. Other information, such as building shape or number of windows, were derived by author assumptions and set the basis for the following reference building dynamic simulations.
Stratigraphies were hypothesized for roof, floor, and walls (two types of wall were simulated for each RB). A lower insulation was performed, as requested in TABULA (U ≈ 0.8 W/m 2 K), and material data and thermo-physical properties were taken from the Italian National Unification (UNI) standards UNI 10355 and UNI 10351 [37], while for resistances the reference norm was UNI 6946 [38]. All structures neighboring with nonheated zones, as suggested in UNI 6946, provided a superficial resistance equal to the internal superficial resistance for each side. A building layout was created starting from available data of floor area, gross volume of building, and compactness (S/V). RB simulations were performed by DesignBuilder software, a graphic interface based on EnergyPlus™, a building energy simulation program used to model both energy consumption and water used in buildings [39]. Heat energy requirement was calculated following TABULA indications and procedures reported in UNI part one [40], in this way, simulations results could be compared to TABULA results (Table 2), where the difference in energy requirement between the two models was in the order of 1% for each RBs. a representative occupant behavior pattern. Internal gains as internal light, electrical equipment, and people were aggregated with a single operating schedule, called "internal gain schedule". Hourly profile of the schedule was extrapolated from UNI 11300 part 1 [40], this was done in order to have a reliable hourly profile.
The simulation results showed a heat requirement of 160 kWh/m 2 /y for AB 1945-1960, 136 kWh/m 2 /y for AB 1961-1975, and 66 kWh/m 2 /y for AB 1976-1990. The hourly simulation provided for each RB (Appendix B) has been normalized with gross heated volume of RB, as done in [43], and then multiplied for respective volumes (Table 1), creating the total heat profile 60% of the total DH heat demand.
To identify necessary retrofit interventions to comply with new energy efficiency requirements, the authors referred to the environmental energy appendix of building regulation of Turin [44]. In the Appendix A, two restriction levels to thermal insulation are described: the first level is about prescriptive limits, while the second one concerns incentive requirements (an economic incentive is settled for the buildings that satisfy a higher thermal insulation to further minimize the energy consumption). Heat loads can be reduced using materials with a low thermal transmittance; it is important to define a global strategy for building thermal insulation choosing a correct insulation material and relative thickness (insulation materials chosen related to prospect 2 in UNI:2015 [37], for hygrometric properties the reference was UNI:2008 prospect 4 [45]). It was assumed also that thermal bridges were corrected by retrofit measures.
Three scenarios were created to analyze the energy savings potential in buildings (Table 3): -Scenario WS: substitution of windows (U = 1.20 W/W, g,gl = 0.30) -Scenario SR (standard renovation): all building envelopes are insulated according to transmittance limits in Level 1 [39] a representative occupant behavior pattern. Internal gains as internal light, electrical equipment, and people were aggregated with a single operating schedule, called "internal gain schedule". Hourly profile of the schedule was extrapolated from UNI 11300 part 1 [40], this was done in order to have a reliable hourly profile.
The simulation results showed a heat requirement of 160 kWh/m 2 /y for AB 1945-1960, 136 kWh/m 2 /y for AB 1961-1975, and 66 kWh/m 2 /y for AB 1976-1990. The hourly simulation provided for each RB (Appendix B) has been normalized with gross heated volume of RB, as done in [43], and then multiplied for respective volumes (Table 1), creating the total heat profile 60% of the total DH heat demand.
To identify necessary retrofit interventions to comply with new energy efficiency requirements, the authors referred to the environmental energy appendix of building regulation of Turin [44]. In the Appendix A, two restriction levels to thermal insulation are described: the first level is about prescriptive limits, while the second one concerns incentive requirements (an economic incentive is settled for the buildings that satisfy a higher thermal insulation to further minimize the energy consumption). Heat loads can be reduced using materials with a low thermal transmittance; it is important to define a global strategy for building thermal insulation choosing a correct insulation material and relative thickness (insulation materials chosen related to prospect 2 in UNI:2015 [37], for hygrometric properties the reference was UNI:2008 prospect 4 [45]). It was assumed also that thermal bridges were corrected by retrofit measures.
Three scenarios were created to analyze the energy savings potential in buildings (Table 3): -Scenario WS: substitution of windows (U = 1.20 W/W, g,gl = 0.30) -Scenario SR (standard renovation): all building envelopes are insulated according to transmittance limits in Level 1 [39] a representative occupant behavior pattern. Internal gains as internal light, electrical equipment, and people were aggregated with a single operating schedule, called "internal gain schedule". Hourly profile of the schedule was extrapolated from UNI 11300 part 1 [40], this was done in order to have a reliable hourly profile.
The simulation results showed a heat requirement of 160 kWh/m 2 /y for AB 1945-1960, 136 kWh/m 2 /y for AB 1961-1975, and 66 kWh/m 2 /y for AB 1976-1990. The hourly simulation provided for each RB (Appendix B) has been normalized with gross heated volume of RB, as done in [43], and then multiplied for respective volumes (Table 1), creating the total heat profile 60% of the total DH heat demand.
To identify necessary retrofit interventions to comply with new energy efficiency requirements, the authors referred to the environmental energy appendix of building regulation of Turin [44]. In the Appendix A, two restriction levels to thermal insulation are described: the first level is about prescriptive limits, while the second one concerns incentive requirements (an economic incentive is settled for the buildings that satisfy a higher thermal insulation to further minimize the energy consumption). Heat loads can be reduced using materials with a low thermal transmittance; it is important to define a global strategy for building thermal insulation choosing a correct insulation material and relative thickness (insulation materials chosen related to prospect 2 in UNI:2015 [37], for hygrometric properties the reference was UNI:2008 prospect 4 [45]). It was assumed also that thermal bridges were corrected by retrofit measures.
Three scenarios were created to analyze the energy savings potential in buildings ( These simulations calculated only the space heat requirements, the "useful energy", while information about energy system typologies and the relative energy efficiencies were further inserted as input data in the energy planning tool.
The energy analysis performed under the previous assumptions is called "standard" evaluation, but to take in care occupant behavior, the intermittency of heating and the operating schedules of electrical equipment, it was necessary to relax the assumption made following the normative in previous simulation. In literature, several studies have considered the influence of occupant behavior on heat consumption [41,42]. In this paper this argument is not stressed, since the aim was to use reference building to simulate the entire built environment of the city, therefore taking into account a representative occupant behavior pattern. Internal gains as internal light, electrical equipment, and people were aggregated with a single operating schedule, called "internal gain schedule". Hourly profile of the schedule was extrapolated from UNI 11300 part 1 [40], this was done in order to have a reliable hourly profile.
The hourly simulation provided for each RB (Appendix B) has been normalized with gross heated volume of RB, as done in [43], and then multiplied for respective volumes (Table 1), creating the total heat profile 60% of the total DH heat demand.
To identify necessary retrofit interventions to comply with new energy efficiency requirements, the authors referred to the environmental energy appendix of building regulation of Turin [44]. In the Appendix A, two restriction levels to thermal insulation are described: the first level is about prescriptive limits, while the second one concerns incentive requirements (an economic incentive is settled for the buildings that satisfy a higher thermal insulation to further minimize the energy consumption). Heat loads can be reduced using materials with a low thermal transmittance; it is important to define a global strategy for building thermal insulation choosing a correct insulation material and relative thickness (insulation materials chosen related to prospect 2 in UNI:2015 [37], for hygrometric properties the reference was UNI:2008 prospect 4 [45]). It was assumed also that thermal bridges were corrected by retrofit measures.
Three scenarios were created to analyze the energy savings potential in buildings ( The WS scenario allows the possibility of reaching a 10% of energy savings with respect to base scenario, 80% with the SR, and 88% with AR, considering the application of the energy saving interventions to all the volume connected to DH represented by the RBs (30.7 Mm 3 ).

Modeling District Heat Hourly Profiles
As a second step, the estimated total heat consumption needed to be calibrated to match the urban energy balance allocated to space heating and district heat production. At first, local district heat production was characterized. By analyzing DH real operational data for the years 2012, 2013, and 2014, average heat losses in the network resulted equal to 10% of total heat production. For district heated buildings, the DH quota dedicated to domestic hot water (DHW) was derived by analyzing the summer load; in particular, the average profile of a summer day (in which the heating systems are off according to the Italian normative) was extrapolated from DH real data and chosen as reference for DHW needs and assumed equal for every day of the year. The calculated DHW demand and the losses were subtracted from the total DH demand. By doing this, the remaining DH heat consumption could be compared to the heat profile created through RB simulations ( Figure 4). Smart Cities 2020, 3 FOR PEER REVIEW 8 -Scenario AR (advance renovation): all building envelopes are insulated according to transmittance limits in Level 2 [39] The WS scenario allows the possibility of reaching a 10% of energy savings with respect to base scenario, 80% with the SR, and 88% with AR, considering the application of the energy saving interventions to all the volume connected to DH represented by the RBs (30.7 Mm 3 ).

Modeling District Heat Hourly Profiles
As a second step, the estimated total heat consumption needed to be calibrated to match the urban energy balance allocated to space heating and district heat production. At first, local district heat production was characterized. By analyzing DH real operational data for the years 2012, 2013, and 2014, average heat losses in the network resulted equal to 10% of total heat production. For district heated buildings, the DH quota dedicated to domestic hot water (DHW) was derived by analyzing the summer load; in particular, the average profile of a summer day (in which the heating systems are off according to the Italian normative) was extrapolated from DH real data and chosen as reference for DHW needs and assumed equal for every day of the year. The calculated DHW demand and the losses were subtracted from the total DH demand. By doing this, the remaining DH heat consumption could be compared to the heat profile created through RB simulations ( Figure 4). Although considering standard values comports uncertainties because both internal loads and infiltration rates are largely influenced by the occupant's behavior (they can vary from building to building, but also from floors of the same building), as analyzed in [46], using a dynamic simulation Although considering standard values comports uncertainties because both internal loads and infiltration rates are largely influenced by the occupant's behavior (they can vary from building to building, but also from floors of the same building), as analyzed in [46], using a dynamic simulation tool is the best solution for taking into consideration occupant's behavior. Figure 4 validates the results from the reference building simulation, with the simulated heat profile (blue line) being almost 60% of the total DH profile (red line).
Building volumes not connected to DH were considered heated by individual boilers, and with the same method heat profiles were created. Retrofit actions (WS, SR, AR) were gradually analyzed in the same way, and different heat load profiles were compared to DH real data ( Figure 5). The decrease of heat consumption comporting an increase in heat loss was determined according to Equation (1) [47]. This equation for the estimation of heat losses was previously applied to the same network as in [26,28,48]. y = 0.7321 − 1.221x + 0.6462x 2 if x < 1 y = 0.1773 − 0.014x if x > 1 (1) in which y = loss (%) and x = linear heat density in MWh/m. Considering 500 km of pipelines, linear density and the correlated heat losses were calculated for each building's energy saving scenario. In the year 2017, the linear density was equal to 3.5 MWh/m, and it decreased in each scenario. With a line density lower than 2 MWh/m, heat losses were considered too high and the DH network became not convenient. SR and AR scenarios provided a linear density lower than 2 MWh/m, which means 15% of heat loss-to increase the linear density up to 2 MWh/m, network expansions scenarios were associated with retrofit measures for each scenario. For the WS scenario, a lower network expansion was considered, equal to 2.9 Mm 3 , and for SR and AR a higher expansion was supposed, equal to 19.5 Mm 3 (the remaining volume of all buildings with individual heating). This kept the linear density higher than 3 MWh/m in each scenario, and maximum heat losses were in the order of 13%. The heat losses were added to the heat profile consumption of DH and individual buildings and, together with DHW consumption, the heat load profile as input for thermal plant was created ( Figure 5).

Energy Savings Costs
Costs associated to energy savings measures were derived from "Prezzario Regione Piemonte 2016" [49]. For each RB, 3 different costs analyses were done, one for each type of building retrofit option (WS, SR, AR). Cost included: construction material/technology cost, installation, and manpower costs.
Analyses with EnergyPLAN also required energy prices and investment cost data as input The decrease of heat consumption comporting an increase in heat loss was determined according to Equation (1) [47]. This equation for the estimation of heat losses was previously applied to the same network as in [26,28,48].
in which y = loss (%) and x = linear heat density in MWh/m. Considering 500 km of pipelines, linear density and the correlated heat losses were calculated for each building's energy saving scenario. In the year 2017, the linear density was equal to 3.5 MWh/m, and it decreased in each scenario. With a line density lower than 2 MWh/m, heat losses were considered too high and the DH network became not convenient. SR and AR scenarios provided a linear density lower than 2 MWh/m, which means 15% of heat loss-to increase the linear density up to 2 MWh/m, network expansions scenarios were associated with retrofit measures for each scenario. For the WS scenario, a lower network expansion was considered, equal to 2.9 Mm 3 , and for SR and AR a higher expansion was supposed, equal to 19.5 Mm 3 (the remaining volume of all buildings with individual heating). This kept the linear density higher than 3 MWh/m in each scenario, and maximum heat losses were in the order of 13%. The heat losses were added to the heat profile consumption of DH and individual buildings and, together with DHW consumption, the heat load profile as input for thermal plant was created ( Figure 5).

Energy Savings Costs
Costs associated to energy savings measures were derived from "Prezzario Regione Piemonte 2016" [49]. For each RB, 3 different costs analyses were done, one for each type of building retrofit option (WS, SR, AR). Cost included: construction material/technology cost, installation, and manpower costs.
Analyses with EnergyPLAN also required energy prices and investment cost data as input (investment, operations, fixed and variable maintenance, and manpower costs). Fuels cost, handling costs, and variable operations and maintenance (O&M) costs were also required for each type of technology, together with the eventual CO 2 price and the discount rate (this was assumed equal to 3%). The additional costs, as energy savings measures and network expansion, were inserted in EnergyPLAN in terms of total annual cost (Appendix A).

Global Cost Analysis
The total annual cost provided by EnergyPLAN (according to Equation (2), as reported in [50]) related to a single year in which all investments were supposed to be applied. In this way, it was possible to compare different annual costs provided from different energy efficiency investments, to understand which is the most profitable solution among the considered options.
with C invest = investment costs; i = interest; n = lifetime. Each investment with a different lifetime could therefore be compared as an equivalent annual cost. Total annual costs (M€/yref) were therefore derived by summing A inves with fixed and variable annual costs.
To analyze the global cost of each scenario, according to a more realistic distribution of the investment along the timeline, the methodology in [28] and in [14] was followed: a time period horizon was defined (2017-2050) and divided into 3 period steps (2017-2025, 2025-2035, 2035-2050).
The renovation rate of buildings was set equal to 3% for each year starting from 2018, while DH expansion was planned with a year rate of 4% starting in 2025 (Table 4). Investments in thermal plants were planned according to the end-life of each system based on starting year [51,52], as explained in Table 5.  Building renovation rate and DH expansion rate were used to create a heat distribution demand for DH and individual buildings for each time step. Time periods were simulated considering a linear CO 2 cost and natural gas cost evolution, as explained in Section 3.5.
Variable annual costs and fixed annual costs for the three periods were assumed considering a discount rate of 3%; a present value factor (Ri) was calculated for each mid-year of the rime periods (4, 13, 25.5), and the total discounted cost at 2017 was calculated according to Equation (3).
with C tot = total cost discounted at 2017 level, C inv (RB,i) = investment cost for building retrofit on period i, C inv (dh,i) = investment cost on DH expansion and heat supply on period i, C fix and C var = fixed and variable costs (€/y) on period i provide from EnergyPLAN, t = years in period I, R(i)= present value factor for mid-year of each time step.

Heat Supply Scenario
The existing synergy between district heating expansion and buildings' energy retrofit interventions has been demonstrated in several studies [19,20,28]. The production of heat in the DH system is possible through several options: renewable energies (e.g., solar and geothermal), heat pumps, conventional fossil fuels, or even by exploiting waste heat from industrial processes. Taking into account the current heat generation mix of the Turin DH system (more than 90% heat production based on gas cogeneration plants), this study focused its attention on combined heat and power (CHP) technologies and on the integration of large-scale heat pumps and storage. Two scenarios were analyzed, according to [28], considering already planned power plant installation (biomass CHP) and commercial available capacities. In future scenarios, CHP plants could operate to balance the electricity needs of heat pump. Biomass CHP was chosen considering the huge availability of biomass, given that the city of Turin is close to the mountain communities, though an analysis on local emissions should be explored.
To compare the proposed demand-supply options with a scenario reflecting current trends, a business as usual (BAU) scenario was considered (any change in heat consumption and production). This scenario was used as a reference to be compared with the other simulations in terms of energy savings, total primary energy supply (TPES), CO 2 emissions, and total annual cost for year 2050. In the BAU scenario, costs for maintenance and replacement of current power plants after the end of their technical lifetime and power plant investment costs were considered, while change in fuel prices and CO 2 price were considered following the same assumptions as the ones made in the Stratego project [22], assumed equally suitable for an Italian context. A linear growth was considered for the carbon price in the scenario timeline, starting with the assumption of a carbon price equal to 6 €/tCO 2 in 2017, as reported in [53], and an optimistic carbon price of 25 €/tCO 2 for 2050, following the hypothesis made on Heat Roadmap 2050, in which a price of 20-30 €/tCO 2 is expected. By combining scenarios on retrofit measures with scenarios on heat generation, different integrated scenarios were created up to 2050 (Table 6). Hourly consumption profiles obtained for individual buildings in the previous section were used as input in EnergyPLAN for describing the future evolution of the heat demand, starting from the current DH real profile provided by the local utility assumed as reference for the scenario BAU. Individual heating systems were considered as supplied by natural gas boilers with an average efficiency of 0.85. Each cost was introduced in the base year simulation for 2017, and a technical optimization (optimizing capacity sizing and energy use) was performed.
The final heat demand (district heated buildings plus individual heating) was 2.65 TWh/y, 1.99 TWh/y for DH and 0.66 TWh/y for individual heating (19.5 Mm 3 ). DH demand was supplied by natural gas CHP (1.93 TWh/y) and natural gas boilers (0.06 TWh/y). Total primary energy supply (TPES) was equal to 6.19 TWh/y, disaggregated into 5.41 TWh/y for CHP + boilers and 0.78 TWh/y for individual heating. The primary energy supply was calculated as the fuel consumptions needed from the energy technologies to satisfy the energy demand of the system. The fuel consumption of the plants was calculated by the model considering technologies' efficiencies. CO 2 emissions were 1.250 Mt, and total annual cost (including investment cost, fixed and variable cost, fuel cost, and carbon price) was 295 M€/y. Three different buildings retrofit scenarios with two different levels of DH expansion options were compared for three different heat supply mixes (Table 7). Table 7. Different combinations of building retrofit, network expansion, and heat production mix. * Network expansion: 1. expansion of 2.9 Mm 3 ; 2. expansion of 19.5 Mm 3 . ** HP: heat pump; CHPB: biomass combined heat and power (CHP); HS: heat storage. The scenarios A, B, and C consider the same supply configuration as the BAU scenario, integrating 3 different level of building retrofit. The introduction of a biomass CHP and heat pump (HP) of the MOD heat supply scenario correlated with the building retrofits in scenarios D (windows substitution), E (standard renovation), and F (advanced renovation). The same retrofits were considered for scenarios G, H, and I, together with ADV heat supply scenarios.

Building Retrofit (WS,SR, AR) and Network Expansion (1,2) * Scenarios
As first analysis, the total annual cost of these scenarios was analyzed for the timeline 2017-2050 according to Equation (3).

Results
At first, the achievable CO 2 emission reduction considering buildings retrofit together with DH expansion without changing in heat supply was analyzed (A, B, C). Compared to the BAU scenario, this reduction is equivalent to 2% for A, 15%, for B and 18% for C ( Figure 6). However, demand side measures alone are not enough to achieve 40% emission reductions to 2030, as requested by [54].
(HP) of the MOD heat supply scenario correlated with the building retrofits in scenarios D (windows substitution), E (standard renovation), and F (advanced renovation). The same retrofits were considered for scenarios G, H, and I, together with ADV heat supply scenarios.
As first analysis, the total annual cost of these scenarios was analyzed for the timeline 2017-2050 according to Equation (3).

Results
At first, the achievable CO2 emission reduction considering buildings retrofit together with DH expansion without changing in heat supply was analyzed (A, B, C). Compared to the BAU scenario, this reduction is equivalent to 2% for A, 15%, for B and 18% for C ( Figure 6). However, demand side measures alone are not enough to achieve 40% emission reductions to 2030, as requested by [54]. Assuming a "moderate change" in the heat generation mix (MOD scenario) can allow a higher reduction in CO 2 emissions (42%) and TPES (40%), thanks to the integration of heat pumps (HP) and the biomass CHP, relying on low carbon electricity (2050 carbon intensity equal to 0.66 Mt/y), maintaining a natural gas CHP working for 57% of the time in cogeneration mode for heat production. The annual operating time of CHP higher than 4000 h can assure a positive profit for the investment. The optimization algorithm coordinates electricity production by CHP units to heat pumps' heat generation-CHP operation follows the electricity need of HP and heat storages are used to store CHP heat production excess. This optimization is related to the decision of balancing only the heat consumption, in future studies it might be interesting to study the impact on the electricity market as well. Furthermore, CHP production is primarily satisfied by biomass CHP to reduce the CO 2 emissions. Natural gas consumption and related carbon emissions decrease further with a higher integration of heat pumps in the ADV scenario. Together with the energy savings due to building retrofit measures, the integration of heat pump in MOD scenario guarantees 21% of TPES reduction, while ADV scenario with 250 MW of heat pumps allows a TPES reduction of 49% (Figure 7). Assuming a "moderate change" in the heat generation mix (MOD scenario) can allow a higher reduction in CO2 emissions (42%) and TPES (40%), thanks to the integration of heat pumps (HP) and the biomass CHP, relying on low carbon electricity (2050 carbon intensity equal to 0.66 Mt/y), maintaining a natural gas CHP working for 57% of the time in cogeneration mode for heat production. The annual operating time of CHP higher than 4000 h can assure a positive profit for the investment. The optimization algorithm coordinates electricity production by CHP units to heat pumps' heat generation-CHP operation follows the electricity need of HP and heat storages are used to store CHP heat production excess. This optimization is related to the decision of balancing only the heat consumption, in future studies it might be interesting to study the impact on the electricity market as well. Furthermore, CHP production is primarily satisfied by biomass CHP to reduce the CO2 emissions. Natural gas consumption and related carbon emissions decrease further with a higher integration of heat pumps in the ADV scenario. Together with the energy savings due to building retrofit measures, the integration of heat pump in MOD scenario guarantees 21% of TPES reduction, while ADV scenario with 250 MW of heat pumps allows a TPES reduction of 49% ( Figure  7).  Cost outputs were further analyzed (Figure 8) in order to understand to what extent each intervention influences the global cost. In the BAU scenario, total annual costs would be around 460 M€/y in 2050. Fuel cost is the prevalent cost item, covering 58% of total annual costs for BAU. Cost outputs were further analyzed (Figure 8) in order to understand to what extent each intervention influences the global cost. In the BAU scenario, total annual costs would be around 460 M€/y in 2050. Fuel cost is the prevalent cost item, covering 58% of total annual costs for BAU. Compared to the BAU scenario, each intervention regarding change in heat production (D, E, F, G, H, I) connected with a reduction in heat demand provides a cost decrease. In the scenarios without changes in heat production (A, B, C), fuel costs represent 51% of total costs. In MOD scenarios (D, E, F), this cost voice decreases to 39% and it reaches a minimum value of 27% in ADV scenarios (G, H, I). Lower fuel costs due to fuel shifting in heat production mitigate investment costs for buildings retrofit and for heat generation interventions.
The ADV scenario with standard renovation in buildings, together with DH expansions, provides lower costs with respect to all the other scenarios simulated. Compared to a standard renovation, an advanced renovation on buildings comports a higher reduction in heat consumption, but also higher investment costs. These investment costs are not offset by the savings, highlighting the critical role of defining correct retrofit measures for existing buildings. Meanwhile, considering the very old building stock and the presence of a district heat network, a standard retrofit scenario (H) can provide a better balancing between new investment costs and decreasing costs as a consequence of energy savings measures. Results from energy system simulations show the synergy between buildings' retrofit interventions and changes in district heating, in terms of network expansion and heat production configuration. From a planning perspective, combining buildings' energy savings with DH expansion strategies, together with the integration of cleaner heat generation technologies as heat pumps, allows significant reductions in terms of CO 2 emission, total costs, and primary energy supply. Larger interventions compared to baseline, as the one proposed in the H scenario, could allow reaching reductions of 78% in CO 2 emission (corresponding to a reduction of 970 MtCO 2 /y) and of 72% in TPES, with benefits in terms of annual costs equivalent to 150 M€/y.
To understand which is the overall marginal cost for the energy system decarbonization in the different scenarios, a "decarbonization cost" Cdec (€/tCO 2 ) was calculated, as in Equation (4).
With Cin(i) = investments costs for scenario i, Cin(BAU) = investment costs in the BAU scenario, CO 2 (i)= CO 2 emission in scenario i, CO 2 (BAU)= CO 2 emission in BAU scenario.
The decarbonization cost is the rate between the increasing in investment cost and the CO 2 reduction due to low-carbon measures. It can be used as an index to identify cost-effective decarbonization investments.
From the results, it is possible to affirm that for scenarios with change in the heat Supply (D-I) the decarbonization cost is lower than 100 €/tCO 2 saving, while for scenarios A, B, and C (no change in the heat production mix, but progressive demand reduction thanks to building renovation) the decarbonization cost is almost triple compared to integrated solutions involving the MOD and ADV scenarios (Figures 9 and 10). In scenarios A, B, and C, investment costs are higher than in BAU due to building retrofits, but CO 2 reduction is not enough to justify investments. This consideration is key for building stock retrofits, where the potential for deep retrofit intervention is often limited by technical feasibility or affordability, which needs to be properly defined.  Despite the fact that the cost of decarbonization can already represent a good economic indicator, greater insights can be derived from a more detailed economic analysis that considers the cost distribution over the whole time horizon. As shown by Figure 11, compared with the annual cost analysis performed with EnergyPLAN, a more detailed global cost analysis shows that scenarios with improved heat supply (MOD and ADV) can deliver bigger cost reductions. On the other hand, by applying this methodology, the longer term benefits of advanced building retrofit (C, F, I) can be observed more precisely [55]. This might help to deliver improved energy retrofit considerations.  Despite the fact that the cost of decarbonization can already represent a good economic indicator, greater insights can be derived from a more detailed economic analysis that considers the cost distribution over the whole time horizon. As shown by Figure 11, compared with the annual cost analysis performed with EnergyPLAN, a more detailed global cost analysis shows that scenarios with improved heat supply (MOD and ADV) can deliver bigger cost reductions. On the other hand, by applying this methodology, the longer term benefits of advanced building retrofit (C, F, I) can be observed more precisely [55]. This might help to deliver improved energy retrofit considerations. Despite the fact that the cost of decarbonization can already represent a good economic indicator, greater insights can be derived from a more detailed economic analysis that considers the cost distribution over the whole time horizon. As shown by Figure 11, compared with the annual cost analysis performed with EnergyPLAN, a more detailed global cost analysis shows that scenarios with improved heat supply (MOD and ADV) can deliver bigger cost reductions. On the other hand, by applying this methodology, the longer term benefits of advanced building retrofit (C, F, I) can be observed more precisely [55]. This might help to deliver improved energy retrofit considerations.
indicator, greater insights can be derived from a more detailed economic analysis that considers the cost distribution over the whole time horizon. As shown by Figure 11, compared with the annual cost analysis performed with EnergyPLAN, a more detailed global cost analysis shows that scenarios with improved heat supply (MOD and ADV) can deliver bigger cost reductions. On the other hand, by applying this methodology, the longer term benefits of advanced building retrofit (C, F, I) can be observed more precisely [55]. This might help to deliver improved energy retrofit considerations. A sensitivity analysis variating discount rates might improve the understanding of which retrofit measures may have reasonable payback times-a lower discount rate presents higher influence of future energy costs, while a higher discount rate presents lower influence on future costs and higher influence of investment costs. In Figure 12, total costs for each scenario were analyzed with different discount rates. Comparing BAU scenario with a moderate scenario (for example E), a higher discount rate comports a lower costs difference between BAU and E, this means that a BAU scenario with few energy efficiency investments can be more convenient with respect to a scenario with bigger energy investments (e.g., E,C) if discount rate increases (or equivalently, without proper policies). Targeted tax policies can represent a useful instrument to push the energy efficiency measures and unlock the necessary fiscal resources for driving a clean energy transition.
Smart Cities 2020, 3 FOR PEER REVIEW 16 Figure 11. Total costs discounted at 2017 level for different scenarios, divided into investment costs and fix + var costs.
A sensitivity analysis variating discount rates might improve the understanding of which retrofit measures may have reasonable payback times-a lower discount rate presents higher influence of future energy costs, while a higher discount rate presents lower influence on future costs and higher influence of investment costs. In Figure 12, total costs for each scenario were analyzed with different discount rates. Comparing BAU scenario with a moderate scenario (for example E), a higher discount rate comports a lower costs difference between BAU and E, this means that a BAU scenario with few energy efficiency investments can be more convenient with respect to a scenario with bigger energy investments (e.g., E,C) if discount rate increases (or equivalently, without proper policies). Targeted tax policies can represent a useful instrument to push the energy efficiency measures and unlock the necessary fiscal resources for driving a clean energy transition.

Discussion
Linking energy system analyses to dynamic building simulations to analyze different urban energy patterns connects a data-driven robust building-level analysis with a technology-rich supplyside understanding, allowing researchers catch their interactions. The opportunity to study together the demand side and the supply side allows a better understanding of the decarbonization potential of a city, especially if characterized by an ancient building stock, as in many European cities. A complete overview of the overall energy system allows municipalities to handle future urban energy planning policies, developing a central role as mediator between different stakeholders from different sectors, being able to lead the construction of a common pathway for achieving future decarbonization goals. The application of EnergyPLAN for national-level studies is well known in

Discussion
Linking energy system analyses to dynamic building simulations to analyze different urban energy patterns connects a data-driven robust building-level analysis with a technology-rich supply-side understanding, allowing researchers catch their interactions. The opportunity to study together the demand side and the supply side allows a better understanding of the decarbonization potential of a city, especially if characterized by an ancient building stock, as in many European cities. A complete overview of the overall energy system allows municipalities to handle future urban energy planning policies, developing a central role as mediator between different stakeholders from different sectors, being able to lead the construction of a common pathway for achieving future decarbonization goals. The application of EnergyPLAN for national-level studies is well known in literature, but the proposed methodology allows us to expand its use to an urban level without renouncing the level of detail required from the urban environment. The methodology can be further improved considering the increasing data availability in the residential sector (cadaster, smart meters, energy certifications), allowing the calibration of building simulations with real data.
EnergyPLAN was chosen because it provides several advantages: it's a free download tool; on the web page is possible find documentations, guides, and tutorials to help users; input data are easily to find; and it has a user-friendly interface. On the other hand, as previously identified, EnergyPLAN shows some limits in the analyses of the supply side. The technical simulation to balance the heat demands applied by the tool relies on a predefined hierarchy of technologies (solar thermal, industrial CHP, heat production from waste fuel, heat plant CHP, heat pumps, peak load boilers), limiting the possibility of choosing other technology combinations or different parameters to optimize the energy system. Although, this limitation can be overcome by utilizing the MATLAB Toolbox for EnergyPLAN [56]. Furthermore, it requires aggregated input data for production that do not allow users to analyze each thermal plant singularly. and this means that the introduction of a new thermal plant (e.g., biomass CHP in this case study) cannot be analyzed individually. Usually power plants are projected to satisfy different cities, while an aggregated simulation requires a sort of close control volume, in which entry flows and exit flows are well known. However, these uncertainties are limited for the previous case study because the studied power plants were designed to satisfy mainly the energy demand of Turin. To extend the methodology to other municipalities, a collaboration with the utilities is suggested to increase the representation of the supply side, as done with this paper for the demand side.
As shown by this paper, a discretization of the urban demand through the identification and simulation of reference buildings is a valuable contribution to improving the understanding of building retrofit potential, the reference buildings approach being useful to identify an average behavior for the residential sector. Despite this, in the future, with growing data availability thanks to the spread of digitalization technologies, data-driven approaches might be more suitable for this purpose, overcoming barriers related to RBs applied to large-scale analyses, such as catching shading profiles of closed buildings, internal loads, and infiltration rates, as explained in [43].
Further improvements in the methodology can be related to improving the representation of the influence that the DH network itself has on the heat demand. Delay effect and heat loss, together with mixing effects due to return flows of each building, comport a heat load profile of thermal plants different from that of buildings, especially from the peak load, which is different for magnitude (lower than user peak) and for duration (peak of thermal plant lasts for hours, for buildings lasts for few minutes). For a better demand and supply relationship evaluation, a future development of the methodology could integrate a tool to simulate the network to analyze load profile of thermal plant and to apply regulation efficiency intervention in building (as the anticipation of the thermal request of buildings to reduce morning peak analyzed in [57]).
Additional improvements in future analysis can be done by integrating the effectiveness of buildings savings related to the increase of temperatures (e.g., lower heating degree days might reduce the costs benefits related to energy retrofit, reducing its profitability) [46].

Conclusions
The aim of the proposed approach was to investigate the relation between buildings' energy savings and low-carbon heat generation strategies at an urban scale, to help municipalities to lead the transition towards a decarbonized energy society through a targeted planning of future investment and energy policies. For the achievement of the European goals in terms of CO 2 emission reduction and primary energy consumption reduction, it is important correlate a demand reduction with a low-carbon supply strategy. The correlation of these two aspects is evident in the residential sector, especially in the historical buildings and general old buildings that represent the highest share of the residential sector in European cities. For these buildings, where the installation of new technologies, such as decentralized heat pumps, is difficult (due to either social, economic, and technical reasons), the role of the district heating, thanks to its easy integration with several low carbon technologies, is very relevant to pursue the goals of the energy transition. The proposed methodology allows an in-depth analysis on effective measures applicable for heat saving at the building level, integrated with an analysis of the energy system supply performed with the EnergyPLAN tool. The integration and expansion of DH in the city was analyzed through different hypothesis of future energy system configurations. The results underline the correlation between the demand side and the supply side, showing that a joint effect of both demand reduction at building level and change in the heat supply is needed to achieve the energy savings goals. Scenarios with changes in heat production provide a big reduction on CO 2 and TPES, due to the integration of heat pumps, and on biomass CHP. In particular ADV scenarios, with the integration of 250 MW heat pumps, 560 MW of heat storage, and 106 MW of biomass CHP comport a CO 2 reduction from a minimum of 56% (without considering buildings savings and DH expansion) to a maximum of 76% considering advanced renovation level for building retrofit together with DH expansion. The integration of large-scale heat pumps in district heating provides a strong reduction in natural gas consumption, and a lower dependence on natural gas assures a safe energy system. New investments in district heating and in building stock are supported by the strong decrease of fuel consumption and its relative costs. However, the cost reduction of these scenarios is supported by a presumed increase in the cost of fossil fuel and in the carbon taxes on 2050. The cost-effectiveness of the energy savings investment is strongly related to the trend of future energy policy decisions that can provide a strong push for embrace a low carbon energy system.

Acknowledgments:
The authors would like to thank Roberta Roberto from ENEA (Italian National Agency for New Technologies, Energy, and Sustainable Economic Development) working for the project IMEAS.

Conflicts of Interest:
The authors declare no conflict of interest.  Table A6. Fuel costs 2017-2050 [21].