How Climate Change A ﬀ ects the Building Energy Consumptions Due to Cooling, Heating, and Electricity Demands of Italian Residential Sector

: Climate change a ﬀ ects the buildings’ performance, signiﬁcantly inﬂuencing energy consumption, as well as the indoor thermal comfort. As a consequence, the growing outdoor environmental temperatures entail a slight reduction in heating consumption and an increase in cooling consumption, with di ﬀ erent overall e ﬀ ects depending on the latitudes. This document focuses attention on the Italian residential sector, considering the current and reduced meteorological data, in anticipation of future climate scenarios. According to a sample of 419 buildings, referring to the climatic conditions of Milan, Florence, Rome, and Naples, the heating and cooling needs are calculated by a simpliﬁed dynamic model, in current and future conditions. The e ﬀ ects of the simplest climate adaptation measure, represented by the introduction of new air conditioners, have been also evaluated. The simulations results show an important reduction in complex energy consumption (Milan − 6%, Florence − 22%, Rome − 25%, Naples − 30%), due to the greater incidence of heating demand in the Italian context. However, the increase in air conditioning electrical consumption over the hot season (Milan + 11%, Florence + 20%, Rome + 19%, Naples + 16%) can play a critical role for the electrical system; for that reason, the introduction of photovoltaic arrays as a compensatory measure have been analysed.


Introduction
Mitigation and adaptation to climate change are key challenges of the 21st century. For this reason the European Union has undertaken numerous initiatives intended to contain its effects in the medium and long term [1][2][3], highlighting the role of energy saving and RES energy production. Within the Paris Agreement document [4] the role of stakeholders in coping with the climate change is strongly recognised, highlighting also the necessary participation of cities, subnational authorities, civil society, and the private sector as well; all these entities are invited to increase their efforts for supporting several initiatives so as to shrink the pollutants emissions along with the vulnerability to the adverse effects of climate change. To do so, regional and international cooperation have to be promoted as much as possible.
Indeed, the global warming is strongly affecting our environment, enhancing both the frequency and the intensity of the unexpected extreme weather events. In accordance with the IPCC [5] the anthropic activities have already induced an average temperature increase of 1 • C more, compared

Materials and Methods
The future projections of building stock energy consumption, in general, do not provide for an explicit distinction between existing and new buildings; taking into account the buildings average life span, it is clear that the current building stock will have to face the medium-term impacts of climate change in the future.
For that reason, this document analyses the implications of climate change on very large and varied sample buildings, both in terms of size and construction features. The sample consists of 419 dwellings and has been created by the collaboration of the students of the Faculty of Architecture of Sapienza University of Rome (Italy).
Useful information for energy characterization has been collected by the online questionnaire, from September 2018. The survey takes into account the following building parameters: (i) The envelope peculiarities in terms of geographic location, surfaces orientation, envelope components U-value, shading devices and ventilation rate; (ii) plants typology (i.e., heating system, cooling system, domestic hot water (DHW) systems, PV modules, solar collectors); (iii) the most common electric devices along with their typical use time-scheduling (i.e., ovens, cooking planes, refrigerators, washing machines, cleaning and ironing, lighting systems, audio/video, personal care, and other equipment); (iv) the occupants number in terms of times and attendance; finally, the actual electricity and natural gas consumptions data hailing from bills.
The same building sample was used in other recent works of the authors so as to identify what are the flexible loads [39], to analyse the energy retrofitting effects on the potential of flexibility on the grid [40], and to evaluate the impact of building automation control systems on energy performance [41]; Section 3.1 provides a more extensive description of the sample characteristics.
The sample dwellings energy performance can be simulated through the same data collection questionnaire, created in an Excel tool by implementing macros and functions written in VBA (Visual Basic for Applications).
Analysing the literature on that topic, it was observed that generally the energy performance of buildings, aimed at the future projection of energy consumption was carried out using dynamic simulation. In many papers [6,9,16] the simulation was performed using leading software (e.g., TRNSYS, EnergyPlus); in other projects the simulation was performed by calculation procedures developed in-house, [7,[42][43][44] validating those outcomes by a direct comparison with the aforementioned tools.
In that work, the authors used their own calculation code for carrying out their analysis. That simulation tool is based on the single thermal zone modelling and it is able to dynamically analyse the building characteristics and performance. More in detail, the heat balance method (HBM) together with a solver algorithm, conduction finite differences-based (CondFD) [45], have been implemented for figuring out the calculations. The self-developed code allows to get a higher flexibility, as well as the easy way to implement additional tools and calculation options if needed.
To examine the sample buildings behaviour in different climatic conditions, it has been assumed as it was located in four different cities, such as Milan, Florence, Rome, Naples ( Figure 1).
In the first phase, the buildings energy simulation has been performed considering the standard climate data related to the reference cities [46] in order to provide a preliminary energy characterization.
Energies 2020, 13, 410 5 of 24 the building characteristics and performance. More in detail, the heat balance method (HBM) together with a solver algorithm, conduction finite differences-based (CondFD) [45], have been implemented for figuring out the calculations. The self-developed code allows to get a higher flexibility, as well as the easy way to implement additional tools and calculation options if needed.
To examine the sample buildings behaviour in different climatic conditions, it has been assumed as it was located in four different cities, such as Milan, Florence, Rome, Naples ( Figure 1).  In the second phase, the simulation has been carried out by increasing uniformly the average monthly temperature for all months (+1.0 • C; +2.0 • C), in order to evaluate the effects of these variations on heating and cooling consumptions; those temperature increases have been chosen accounting for the average temperature registered in the last five years in the four sample cities. It is important to highlight that all of simulations consider only the outdoor temperature changes neglecting the relative humidity modifications. These latter, in some cases, can lead to passive solar gains lessening due to the potential fog formation. Anyway, that issue has to be considered depending on the location latitude. Figure 2 shows the average monthly temperature related to the reference cities, comparing the simulated values with the monthly average values registered in the time span of 2014-2019 [47].
Although the Italian standards are very recent (2016) [46] in the last five years, for all cities, the outdoor temperatures, on average, were higher than the reported ones (i.e., Milan: +0.5 • C; Florence: +0.2 • C; Rome: +1.2 • C; Naples: +1.0 • C). Moreover, what is relevant pertain to the summer average values (June, July, August) of Milan (+0.7 • C), Rome (+1.8 • C), and Naples (+1.2 • C), while for Florence there are no differences referring to the standard values (indeed it being the smallest of the examined cities and it is likely less affected by the heat island effects).
In the third phase, the effects (on cooling consumption) of the simplest and the most widespread adaptation measure, aimed at accomplishing the summer comfort, has been considered. To do so, the installation of small air-to-air cooling systems have been assumed (energy label A++). In detail, the introduction of further new air conditioners (i.e., +1; +2) have been fixed, in addition to the existing ones; furthermore, an upper limit has been assumed for the total number of air conditioners, having considered one unit for each 20 m 2 of floor surface. It is worth noticing how in this work the small retrofitting interventions have been considered. Such measures consist of technical options entailing a capital expenditure less than 2000 €. For that reason, other refurbishment interventions on older buildings have not been accounted for (i.e., reduction of U-value of building envelope and roof). Moreover, for the same economic reason, the new cooling system installation along with a ventilation one for the whole flat have been neglected and, consequently, they have not been analysed. Table 1 shows a summary of the simulated scenarios. important to highlight that all of simulations consider only the outdoor temperature changes neglecting the relative humidity modifications. These latter, in some cases, can lead to passive solar gains lessening due to the potential fog formation. Anyway, that issue has to be considered depending on the location latitude. Figure 2 shows the average monthly temperature related to the reference cities, comparing the simulated values with the monthly average values registered in the time span of 2014-2019 [47].    The fourth phase, provides, in a proactive way, a potential solution to the problems hailing from the average outdoor temperature increase; that solution essentially consists of the installation of new photovoltaic systems to compensate the greater electricity demand over the hot season. The required gross surface on the roof for PV modules installation, including the distance to avoid the mutual shading, have been computed using the average indicator equal to 19.8 m 2 /kW [48,49]. Thus, the yearly Energies 2020, 13, 410 7 of 24 PV capability have been calculated using the online tool PVGIS developed by a previous European Commission program [50].

Dwellings Description and General Analysis on Consumptions Typology
The sample buildings consists of 419 dwellings with uneven characteristics in terms of construction year, size, and occupation rate [39].
Referring to the construction year classification, the sample buildings composition is basically equal to the Italian building stock one. Thus, it has been decided to subdivide all those periods after 1991 differently from what was reported by the National Institute of Statistics [51,52]; indeed, in Italy starting from 1991 several regulations dealing with energy saving in the residential sector were promoted subsequently (see Figure 3).  All of the dwelling's characteristics registration has been made by the use of a specific questionnaire survey dedicated to nonexpert users as well (see Appendix A, Table A1). The typical building technical parameters, such as U-value, thermal inertia, etc., have been deduced from the construction age, the climate zone, and the potential energy refurbishment interventions. Moreover, the solar gains have been evaluated taking into account the walls and roofs colour (i.e., very light colour, light colour, medium colour, dark colour, very dark colour), as well as the shading degree in terms of time periods over the day.
Additionally, Figure 3 shows how the larger part (64.1%) of the sample buildings dates back to 1976, since just in that year the Italian Government issued the first version of a technical specification for improving the buildings energy saving; only 17 flats (4.1%) were built more recently, i.e., after 2005 (Table 2). More than one half of sample buildings (55.1%) underwent light refurbishment interventions, which most frequently consisted of windows substitution (47.0%).  All of the dwelling's characteristics registration has been made by the use of a specific questionnaire survey dedicated to nonexpert users as well (see Appendix A, Table A1). The typical building technical parameters, such as U-value, thermal inertia, etc., have been deduced from the construction age, the climate zone, and the potential energy refurbishment interventions. Moreover, the solar gains have been evaluated taking into account the walls and roofs colour (i.e., very light colour, light colour, medium colour, dark colour, very dark colour), as well as the shading degree in terms of time periods over the day. Additionally, Figure 3 shows how the larger part (64.1%) of the sample buildings dates back to 1976, since just in that year the Italian Government issued the first version of a technical specification for improving the buildings energy saving; only 17 flats (4.1%) were built more recently, i.e., after 2005 ( Table 2). More than one half of sample buildings (55.1%) underwent light refurbishment interventions, which most frequently consisted of windows substitution (47.0%). From data analysis, five building dimensional classes, in terms of floor surface, have been identified, so that Figure 4a outlines their frequency within the sample. Specifically, small dwellings are characterised by a floor surface lower than 50 m 2 , the small-medium class refers to 50-85 m 2 , the medium one to 85-115 m 2 , while the medium-large class ranges in 115-145 m 2 . Finally, all those dwellings larger than 145 m 2 have been considered the upper class. By averaging data it emerges how the most common class is the medium one since the average floor surface is equal to 112.4 m 2 . Moreover, 38.8% of building sample are occupied by four people, while 26.0% by three people, and 20.6% by two people, as reported in Figure 4b. That entails an average occupants' number equal to 3.2.  The heating system along with the domestic hot water production device has been installed in all dwellings sample. The main typology consists of autonomous boilers (73.3%) which are mostly fuelled with natural gas (98.8%). With regards to domestic hot water production systems, natural gas is still the main feeding fuel (85.4%), since those devices are very often integrated within the aforementioned boilers (see Figure 5a). Then, the heat emission terminals consists typically of radiators (95.9%) while low temperature devices, such as fan-coils and radiant floors, have been registered only in a few cases, representing only 1.5% and 2.7%, respectively.  The heating system along with the domestic hot water production device has been installed in all dwellings sample. The main typology consists of autonomous boilers (73.3%) which are mostly fuelled with natural gas (98.8%). With regards to domestic hot water production systems, natural gas is still the main feeding fuel (85.4%), since those devices are very often integrated within the aforementioned boilers (see Figure 5a). Then, the heat emission terminals consists typically of radiators (95.9%) while low temperature devices, such as fan-coils and radiant floors, have been registered only in a few cases, representing only 1.5% and 2.7%, respectively. all dwellings sample. The main typology consists of autonomous boilers (73.3%) which are mostly fuelled with natural gas (98.8%). With regards to domestic hot water production systems, natural gas is still the main feeding fuel (85.4%), since those devices are very often integrated within the aforementioned boilers (see Figure 5a). Then, the heat emission terminals consists typically of radiators (95.9%) while low temperature devices, such as fan-coils and radiant floors, have been registered only in a few cases, representing only 1.5% and 2.7%, respectively. In the end, fixed chillers for inner space cooling have been installed only in 207 dwellings showing that the deployment of those devices is limited to 49.4%. Specifically, only few rooms are served by split units and they are characterised also by high values in energy label, since they have been installed recently, as reported in Figure 5c. Indeed, more than 65% of them have an efficiency class equal to or greater than class A+. Thereafter, only in 77 dwellings (18.4%) there are portable equipment such as fans or dehumidifiers.  Figure 5. Equipment: (a) Heating systems typology; (b) chilled rooms; (c) energy label of cooling system.
In the end, fixed chillers for inner space cooling have been installed only in 207 dwellings showing that the deployment of those devices is limited to 49.4%. Specifically, only few rooms are served by split units and they are characterised also by high values in energy label, since they have been installed recently, as reported in Figure 5c. Indeed, more than 65% of them have an efficiency class equal to or greater than class A+. Thereafter, only in 77 dwellings (18.4%) there are portable equipment such as fans or dehumidifiers.
During the simulation process, the reference value for the indoor comfort temperature in the summertime has been fixed equal to 26 • C. Beyond that threshold and depending on the occupants' number, the chilling units will be switched-on. Notwithstanding, due to climate change effects, discomfort conditions will occur in a huge number of Italian dwellings, since cooling devices are not widespread as stated above. For that reason, the installation of further new air conditioners have been hypothesised, in the present analysis, as the easiest adaptation measure to overcome those features.
The spreading out of both heating and cooling systems in the sample buildings is similar to that was reported in the ISTAT survey on household energy consumption [53]; that document also reports the diffusion of system types on the Italian territory according to the area (North-West, North-East, Middle Regions, South); in regards to cooling systems, the document reported an average diffusion of 70.8%, with weak variations at territorial level (North 66.5%; South 67.8%; Middle Regions 76.0%), apparently not influenced by climatic conditions.
The energy characterization reported in [39] shows that energy consumption for the inner space heating is on average equal to 43.5% of the total; energy consumption for cooling, in dwellings where the service is present, is on average 3.6%; the other services have a total incidence of 52.9% characterised by relative shares equal to 14.1% for DHW, 12.4% for cooking purpose, 5.6% for washing machines, and 4.1% for computer/Internet refrigeration 3.8%; thereafter the cleaning and ironing use represents only 3.5% while the care person lighting, audio/video, and the other equipment are equal to 2.9%, 2.7%, 1.8%, and 1.9%, respectively. For more information on the building sample characteristics and its energy characterization see [39][40][41].
It is important to highlight that the calculation tool validation has been carried out comparing the simulation results with the real consumption data hailing from the energy bills [39]. Furthermore, the predictive model is characterised by high values of correlation coefficient R 2 , i.e., 0.8993 for the electricity needs and 0.7716 for the natural gas volumes.

Energy Characterization of Dwellings in the Four Considered Cities
Considering the standard climatic data of the four cities, simulations have been performed to identify the heating consumption, cooling, and other services; in the following charts, the buildings are grouped by dimensional classes, in order to verify any differences depending on the size; the results are expressed in the form of primary energy and are normalised by the floor surface.
With regards to heating consumption ( Figure 6), for the reference cities and accounting for the different size, no significant modifications are observed, except for a slightly lower specific consumption for larger flats.
the predictive model is characterised by high values of correlation coefficient R 2 , i.e., 0.8993 for the electricity needs and 0.7716 for the natural gas volumes.

Energy Characterization of Dwellings in the Four Considered Cities
Considering the standard climatic data of the four cities, simulations have been performed to identify the heating consumption, cooling, and other services; in the following charts, the buildings are grouped by dimensional classes, in order to verify any differences depending on the size; the results are expressed in the form of primary energy and are normalised by the floor surface.
With regards to heating consumption ( Figure 6), for the reference cities and accounting for the different size, no significant modifications are observed, except for a slightly lower specific consumption for larger flats.
Nevertheless, a noncontinuous trend has been registered. Furthermore, as it was expected, specific heating consumption is strongly influenced by the degree days; Milan has the highest consumption, on average equal to 105.0 kWh/ym 2 ; whilst Florence (70.5 kWh/ym 2 ), Rome (64.1 kWh/ym 2 ), Naples (54.5 kWh/ym 2 ). Similarly, with reference to cooling consumption (Figure 7), having considered the different size, no significant variations are observed, except for a slightly lower specific consumption for larger flats. In that case too, a trend that is not continuous has been found, except for Milan. However, a weaker dependence of energy consumption on the cities geographical location can be noticed (Milan 1.3 kWh/ym 2 ; 1.8 kWh/ym 2 ; Rome 1.8 kWh/ym 2 ; Naples 1.3 kWh/ym 2 ). Nevertheless, a noncontinuous trend has been registered. Furthermore, as it was expected, specific heating consumption is strongly influenced by the degree days; Milan has the highest consumption, on average equal to 105.0 kWh/ym 2 ; whilst Florence (70.5 kWh/ym 2 ), Rome (64.1 kWh/ym 2 ), Naples (54.5 kWh/ym 2 ).
Similarly, with reference to cooling consumption (Figure 7), having considered the different size, no significant variations are observed, except for a slightly lower specific consumption for larger flats. In that case too, a trend that is not continuous has been found, except for Milan. However, a weaker dependence of energy consumption on the cities geographical location can be noticed (Milan 1.3 kWh/ym 2 ; 1.8 kWh/ym 2 ; Rome 1.8 kWh/ym 2 ; Naples 1.3 kWh/ym 2 ). Similarly, with reference to cooling consumption (Figure 7), having considered the different size, no significant variations are observed, except for a slightly lower specific consumption for larger flats. In that case too, a trend that is not continuous has been found, except for Milan. However, a weaker dependence of energy consumption on the cities geographical location can be noticed (Milan 1.3 kWh/ym 2 ; 1.8 kWh/ym 2 ; Rome 1.8 kWh/ym 2 ; Naples 1.3 kWh/ym 2 ). Referring to the charts depicted in Figure 7, it is worth noticing that only 49.4% of the buildings sample are equipped with cooling systems. Moreover, the values shown there include the dwellings without a system, where consumption is necessarily equal to zero.
For the other energy consumptions (Figure 8), a decreasing trend can be noted, starting from the smaller flats (120.0 kWh/ym 2 ) to the larger ones (55.0 kWh/ym 2 ); from a geographical point of view Referring to the charts depicted in Figure 7, it is worth noticing that only 49.4% of the buildings sample are equipped with cooling systems. Moreover, the values shown there include the dwellings without a system, where consumption is necessarily equal to zero.
For the other energy consumptions (Figure 8), a decreasing trend can be noted, starting from the smaller flats (120.0 kWh/ym 2 ) to the larger ones (55.0 kWh/ym 2 ); from a geographical point of view there are no relevant variations. It is important to point out that the other energy uses, in this model, are not influenced by the climatic conditions and therefore, in this study, they have been considered constant for all the examined scenarios.  According to the previous statements, the impact of heating and cooling needs on the dwelling energy consumption are summarised in Figure 9; this incidence varies according to the location and the house size; the heating consumption in Milan has a variable incidence ranging between 49.4% and 61.0%; in Florence that span is limited to 40.2% and 51.1%; in Rome it is comprised between 37.9% and 48.9%; finally, in Naples, the range is equal to 34.2-45.6%. The incidence of cooling energy consumption is much lower (Milan 0.4% ÷ 1.0%; Florence 1.0 % ÷ 1.4%; Rome 1.0% ÷ 1.4%; Naples 0.8% ÷ 1.1%). According to the previous statements, the impact of heating and cooling needs on the dwelling energy consumption are summarised in Figure 9; this incidence varies according to the location and the house size; the heating consumption in Milan has a variable incidence ranging between 49.4% and Energies 2020, 13, 410 12 of 24 61.0%; in Florence that span is limited to 40.2% and 51.1%; in Rome it is comprised between 37.9% and 48.9%; finally, in Naples, the range is equal to 34.2-45.6%. The incidence of cooling energy consumption is much lower (Milan 0.4% ÷ 1.0%; Florence 1.0 % ÷ 1.4%; Rome 1.0% ÷ 1.4%; Naples 0.8% ÷ 1.1%). According to the previous statements, the impact of heating and cooling needs on the dwelling energy consumption are summarised in Figure 9; this incidence varies according to the location and the house size; the heating consumption in Milan has a variable incidence ranging between 49.4% and 61.0%; in Florence that span is limited to 40.2% and 51.1%; in Rome it is comprised between 37.9% and 48.9%; finally, in Naples, the range is equal to 34.2-45.6%. The incidence of cooling energy consumption is much lower (Milan 0.4% ÷ 1.0%; Florence 1.0 % ÷ 1.4%; Rome 1.0% ÷ 1.4%; Naples 0.8% ÷ 1.1%).

Changes in Heating and Cooling Energy Consumption Due to Temperature Enhancement
This section illustrates the simulations results to identify variations associated to heating and cooling energy consumption caused by the outdoor temperature, once average increases equal to 1 °C (scenario #1.0) and 2 °C (scenario #2.0) are assumed.
For Milan, referring to the heating consumption ( Figure 10a) in scenario #1.0 there is a reduction of 11.0%, while in scenario #2.0 the reduction is equal to 21.8%; the difference related to the dwelling size are minimal. Conversely, considering the cooling consumption (Figure 10b)  Other Heating Cooling Figure 9. Specific primary energy consumption for heating, cooling, and other.

Changes in Heating and Cooling Energy Consumption Due to Temperature Enhancement
This section illustrates the simulations results to identify variations associated to heating and cooling energy consumption caused by the outdoor temperature, once average increases equal to 1 • C (scenario #1.0) and 2 • C (scenario #2.0) are assumed.
For Milan, referring to the heating consumption ( Figure 10a) in scenario #1.0 there is a reduction of 11.0%, while in scenario #2.0 the reduction is equal to 21.8%; the difference related to the dwelling size are minimal. Conversely, considering the cooling consumption (Figure 10b  Finally, referring to Naples, the heating consumption (Figure 13a) lessens up to 16.6% in scenario #1.0, while in scenario #2.0 that reduction is equal to 33.1%; the dwelling floor surface affects weakly the differences. With regards to the cooling consumption (Figure 13b   Finally, referring to Naples, the heating consumption (Figure 13a) lessens up to 16.6% in scenario #1.0, while in scenario #2.0 that reduction is equal to 33.1%; the dwelling floor surface affects weakly the differences. With regards to the cooling consumption (Figure 13b  Finally, referring to Naples, the heating consumption (Figure 13a) lessens up to 16.6% in scenario #1.0, while in scenario #2.0 that reduction is equal to 33.1%; the dwelling floor surface affects weakly the differences. With regards to the cooling consumption (Figure 13b), there is an increase of 45.2% in scenario #1.0, while 101.3% is accomplished in scenario #2.0; also in this case the floor surface values lead to limited variations: The increase is larger for small houses (#1.0 + 48.9%; #2.0 + 108.5%) and smaller for large houses (#1.0 + 33.9%; #2.0 + 76.0%). As it was expected, the simulation results show a reduction in heating consumption and an increase in cooling consumption. Basically, it is possible to state that heating consumption tends to decrease, while cooling consumption grows up, and they are not strongly influenced by floor surface. Indeed, the percentage reductions in the cold season are strictly related to the location and are smaller where the degree days are greater. More generally, the percentage reductions are very similar for all case studies except for Naples, where they are higher. As it was expected, the simulation results show a reduction in heating consumption and an increase in cooling consumption. Basically, it is possible to state that heating consumption tends to decrease, while cooling consumption grows up, and they are not strongly influenced by floor surface. Indeed, the percentage reductions in the cold season are strictly related to the location and are smaller where the degree days are greater. More generally, the percentage reductions are very similar for all case studies except for Naples, where they are higher.

Change in Cooling Energy Consumption with the Addition of Air Conditioners
This section illustrates the simulations results to identify the enhancement of cooling energy need, when outdoor environmental conditions change. Six different scenarios have been built once average temperature increases equal to 1 and 2 • C are implemented along with an additional air conditioner (scenarios #0.1, #1.1, #2.1) or two air conditioners (scenarios #0.2, #1.2, #2.2); the variations in cooling consumption, with the existing equipment, (scenarios #0.0, #1.0, #2.0), have been already examined in the previous section and are reported here for convenience of reading (Figures 13a, 14a, 15a and 16a).   For Milan (Figure 14), the main findings are listed below: • For small houses in scenario #0.0 the consumption for cooling is equal to 0.9 kWh/ym 2 ; in scenario #2.1 they are 5.9 kWh/ym 2 ; in scenario #2.2 they are 6.0 kWh/ym 2 ; percentage variations are equal to 556% and 567%, respectively; • For large dwellings in scenario #0.0 the consumption for cooling is equal to 1.4 kWh/ym 2 ; in scenario #2.1 they are 3.7 kWh/ym 2 ; in scenario #2.2 they are equal to 4.1 kWh/ym 2 ; the percentage variations rise up to 164% and 193%, respectively; • On average in scenario #0.0 the consumption for cooling is equal to 1.3 kWh/ym 2 ; in scenario #2.1 they are equal to 4.3 kWh/ym 2 ; in scenario #2.2 they are 4.7 kWh/ym 2 ; the relative variations are 231% and 262%, respectively.
For Florence case study ( Figure 15): • For small dwellings in scenario #0.0 the consumption for cooling is 1.9 kWh/ym 2 ; in scenario #2.1 they are equal to 9.2 kWh/ym 2 ; in scenario #2.2 they are equal to 9.3 kWh/ym 2 ; percentage variations are equal to 384% and 389%, respectively; • For large dwellings in scenario #0.0 the consumption for cooling is equal to 1.7 kWh/ym 2 ; in scenario #2.1 they are 4.4 kWh/ym 2 ; in scenario #2.2 they are 4.9 kWh/ym 2 ; percentage variations rise up to 159% and 188%, respectively; • On average in scenario #0.0 the consumption for cooling is equal to 1.8 kWh/ym 2 ; in scenario #2.1 they are 5.8 kWh/ym 2 ; in scenario #2.2 they are equal to 6.4 kWh/ym 2 ; relative variations are 222% and 256%.
For Rome case study ( Figure 16): • For small dwellings in scenario #0.0 the consumption for cooling is 1.9 kWh/ym 2 ; in scenario #2.1 they are 9.0 kWh/ym 2 ; in scenario #2.2 they are 9.1 kWh/ym 2 ; percentage variations are equal to 374% and 379%, respectively; • For large dwellings in scenario #0.0 the consumption for cooling is 1.6 kWh/ym 2 ; in scenario #2.1 they are 4.2 kWh/ym 2 ; in scenario #2.2 they are 4.7 kWh/ym 2 ; percentage variations rise up to 163% and 194%, respectively; • On average in scenario #0.0 the consumption for cooling is equal to 1.8 kWh/ym 2 ; in scenario #2.1 they are equal to 5.6 kWh/ym 2 ; in scenario #2.2 they are 6.2 kWh/ym 2 ; relative variations are 211% and 244%, respectively.
For Naples case study ( Figure 17): • For small houses in scenario #0.0 the consumption for cooling is 1.4 kWh/ym 2 ; in scenario #2.1 they are equal to 7.7 kWh/ym 2 ; in scenario #2.2 they are 7.8 kWh/ym 2 ; percentage variations are equal to 450% and 457%, respectively; • For large dwellings in scenario #0.0 the consumption for cooling is 1.2 kWh/ym 2 ; in scenario #2.1 they are equal to 3.5 kWh/ym 2 ; in scenario #2.2 they are equal to 3.9 kWh/ym 2 ; percentage variations rise up to 192% and 225%; • On average in scenario #0.0 the consumption for cooling is equal to 1.3 kWh/ym 2 ; in scenario #2.1 they are 4.7 kWh/ym 2 ; in scenario #2.2 they are equal to 5.2 kWh/ym 2 ; relative variations are 262% and 300%. In all case studies and scenarios, the simulations outcomes show important variations in cooling consumption owing to temperature changes and to the installation of further air conditioner; this latter, entails a moderate increase in energy consumption, with significant effects only for large houses.
In percentage terms, when temperature profiles shift up of 2 °C and two additional air conditioners are considered, variation values exceed 350% for small houses and 150% for the large ones; on average such variations range between 211% and 300%. In absolute terms, the energy consumption associated to the cooling needs still remains low, showing also a modest impact on the overall consumption.

Change in Overall Energy Consumption for Simulated Scenarios
In this section a comprehensive overview of simulations results has been presented and discussed; specifically, Table 3 shows the total primary energy consumption, by summing the In all case studies and scenarios, the simulations outcomes show important variations in cooling consumption owing to temperature changes and to the installation of further air conditioner; this latter, entails a moderate increase in energy consumption, with significant effects only for large houses.
In percentage terms, when temperature profiles shift up of 2 • C and two additional air conditioners are considered, variation values exceed 350% for small houses and 150% for the large ones; on average such variations range between 211% and 300%. In absolute terms, the energy consumption associated to the cooling needs still remains low, showing also a modest impact on the overall consumption.

Change in Overall Energy Consumption for Simulated Scenarios
In this section a comprehensive overview of simulations results has been presented and discussed; specifically, Table 3 shows the total primary energy consumption, by summing the primary energy consumption for heating, cooling, and for other uses; for each simulated scenario, the relative change compared to the reference scenario #0.0 is reported as well. In previous sections it has been demonstrated that an increase in the outdoor temperature leads to the lessening of consumption for heating purposes, while for the cooling ones it tends to enhance. Notwithstanding, the increase in energy consumption for cooling is greater if the addition of more air conditioners is considered. It is important to highlight how the incidence of heating consumption is higher, and therefore, adding the two variations, it is possible to get better energy gains globally. That circumstance occurs for all scenarios and for all case studies; only scenarios #0.1 and #0.2, referred to the Milan case, are the exceptions. Substantially, for the major part of simulated scenarios, reduction values beyond 20% can be accomplished anyhow.
On the basis of a yearly energy balance, it can be stated that, in the residential sector, the outdoor temperature shift-up favours the overall energy consumption diminishing. As a matter of fact, that feature is mainly caused by the conspicuous reduction in heating consumption. More specifically, referring to the Italian environmental context, it is possible to forecast a marked decrease in local polluting emissions, since the largest part of the heating systems is based on gas-fired technologies; that occurrence certainly represents an important benefit for the air quality, especially in the urban areas.

Monthly Analysis of Cooling Consumption: Criticalities and Foreseeable Solutions
The analysis over one-year period, however, does not allow to highlight properly all those issues that the growing energy consumption for cooling purposes implies on the electricity system in the summertime.
For the typical Italian climatic conditions, the cooling season lasts no longer than four months, i.e., June, July, August, and September; in these months, particularly in July and August, the electricity consumption growth for chilling the inner spaces is concentrated.
Having considered all dwellings and referring also just to the summer months, Table 4 shows accurately the relative variation of total electricity consumption. Those numerical values come out by comparing all the simulated scenarios with the reference scenario #0.0. From data it emerges how July and August are the most critical months for the electrical grid. The outdoor temperature shift-up along with the installation of new air conditioners can lead to an increase in electricity consumption of over 20%. Furthermore, it is well-known how the month of July is the time period characterised by the greatest electricity off-takes [54]. Therefore, a further increase could affect negatively the grid stability or the capacity availability, especially in large cities.
A viable option to overcome those criticalities is to adopt largely self-production, as well as self-consumption of electrical energy, for instance, hailing from photovoltaic systems installed on building roofs. Once well-oriented PV arrays (in terms of slope and azimuth, focused on the maximum yearly production) are assumed, and they are made of crystalline silicon, by the annual capability it is possible to deduce the required receiving surfaces (using the average surface equal to 19.8 m 2 /kW). Those surfaces have to match the additional electricity needs which have been calculated before. Table 5 shows the results of such estimates indicating the relative surface increase when the flat surface of each apartment is considered as reference. Referring to the most critical month for the electricity system (July), the estimates indicate that the surface of new photovoltaic systems (to meet the additional demand for electricity) can exceed 4.5% of the flat surface related to the single apartment.
Those data represent average values which are usable only for preliminary evaluations; yet, in the urban context, they are very often approximated since: i.
There are typically multi-storey buildings in cities; consequently, the usable roof surface must be divided by the various floors; ii.
Not all roofs can be used for PV installation (surfaces already used, architectural constraints, poor irradiation conditions); iii. The modern construction strategies lead to build new apartments characterised by a not wide floor surface (i.e., small < 50 m 2 ; small-medium 50-85 m 2 ), as well as by high specific consumptions for air conditioning, it being even double compared to the larger ones (see .
Whether some of the adverse circumstances listed above occurred at the same time, the buildings' roof surface will be likely insufficient to figure out the criticalities pointed out in this article.
For those reasons, the potential solution to the outdoor temperature shift-up must be identified on a broad territorial scale, involving a dwellings number as large as possible. Thereafter, an inclusion process, taking into account what has just been presented and discussed, it should be developed.
In the end, depending on the territorial context where dwellings are located (e.g., city centre, outskirts, countryside), taking into account the presence of cultural heritage and listed buildings, as well as the prevalent construction typology (e.g., multi-storey building, terraced buildings, single dwellings, etc.), the best territorial scale will be defined (small group of buildings, neighbourhood, urban cell, whole city). To do so, the most appropriate subject able to coordinate the interventions, such as cooperatives, energy communities and municipalities, has to be involved.

Conclusions
In this work, the effects of an outdoor temperature shift up due to climate change, on the energy consumption of the Italian residential sector have been assessed. Using a statistical sample of 419 dwellings, the demand for heating and cooling has been calculated for different scenarios by a simplified dynamic simulation model. Those scenarios consider an increase in the outdoor temperature together with the simpler measure of climate adaptation, represented by the introduction of new air conditioners. Thus, four different climatic conditions have been assumed in order to calculate the sample buildings energy balance. The environmental data refer to Milan, Florence, Rome, and Naples.
The simulation results indicate a general reduction in energy consumption for all scenarios; the reduction in heating demand is much higher than the increase in cooling demand. On the basis of a one-year period the achievable primary energy savings can be even higher than 20%.
That result is undoubtedly positive since it is caused by the particular climatic conditions and equipment of Italian residences (heating systems are installed everywhere, and they serve the whole house; cooling systems are present only in part of the sample buildings and only in a few rooms). In terms of pollutants emissions, a general reduction can be accomplished since the local emissions are strictly related to natural gas-based heating systems.
Nevertheless, it emerges how the increase in cooling demand (due to the introduction of new air conditioners) could generate technical issues in summer months, when the greatest electricity demand occurs. As a matter of fact, the simulations indicate an increase in electricity demand more than 20% during July and August, for many scenarios.
For that reason, the PV systems introduction has been analysed as a compensatory measure, calculating the required receiving surface, as well as the gross surface for the installation. Referring to July, the estimates indicate that in many cases that area can exceed 4.5% of the flat surface related to the single house. Additionally, the single building roof might be not enough, in all those cases which are particularly unfavourable, such as multi-storey buildings, not usable roof surfaces and small houses.
This latter occurrence suggests that the broad territorial scale approach has to be adopted for figuring out all issues related to the outdoor temperature shift up. Thereafter, the most suitable dwellings for the interventions have to be identified along with all of those subjects able to coordinate effectively the activities. That choice has to be made among cooperatives of users, energy communities, or municipalities. Acknowledgments: This work is a part of a wider research activity dealing with: "Study of an aggregator model for a smart district". The project has been carried out in cooperation with ENEA-DTE-SEN-SCC (Italian National Agency for New Technologies, Energy and Sustainable Economic Development-Department of Energy Technologies-Smart Energy Division) and CITERA (Sapienza University of Rome-Interdepartmental Research Center for Territory, Construction, Restoration and Environment). The aforementioned institutions are gratefully acknowledged by the authors for their support and funding.

Conflicts of Interest:
The authors declare no conflict of interest.