Impact of Climate Change on the Energy and Comfort Performance of nZEB: A Case Study in Italy

Climate change is posing a variety of challenges in the built realm. Among them is the change in future energy consumption and the potential decay of current energy efficient paradigms. Indeed, today’s near-zero Energy buildings (nZEBs) may lose their virtuosity in the near future. The objective of this study is to propose a methodology to evaluate the change in yearly performance between the present situation and future scenarios. Hourly dynamic simulations are performed on a residential nZEB located in Rome, built in compliance with the Italian legislation. We compare the current energy consumption with that expected in 2050, according to the two future projections described in the Fifth Assessment Report (AR5) by the Intergovernmental Panel on Climate Change (IPCC). Implications for thermal comfort are further investigated by assuming no heating and cooling system, and by tracking the free-floating operative temperature. Compared to the current weather conditions, the results reveal an average temperature increase of 3.4 ◦C and 3.9 ◦C under RCP4.5 and RCP8.5 scenarios, estimated through ERA-Interim/UrbClim. This comes at the expense of a 47.8% and 50.3% increase in terms of cooling energy needs, and a 129.5% and 185.8% decrease in terms of heating needs. The annual power consumption experiences an 18% increase under both scenarios due to (i) protracted activation of the air conditioning system and (ii) enhanced peak power requirements. A 6.2% and 5.1% decrease in the hours of adaptive comfort is determined under the RCP4.5 and RCP8.5′s 2050 scenarios out of the concerted action of temperature and solar gains. The results for a newly proposed combined index for long-term comfort assessments reveal a milder future penalty, owing to less pronounced excursions and milder daily temperature swings.


Introduction
CO 2 emissions are causing a prolonged and clear increase in global temperatures [1]. In 2010, the building sector accounted for about 32% of global energy consumption, 19% of CO 2 emissions and 51% of global electricity consumption [1]. If, on one side, buildings and their related activities are responsible for a significant portion of greenhouse gas emissions, on the other side they represent a great opportunity for mitigation and adaptation to climate change effects [2].
In order to limit the temperature increase to 2 • C compared to pre-industrial levels, the Fifth Assessment Report (AR5) by the Intergovernmental Panel on Climate Change (IPCC) considers four different future scenarios (RCP; Representative Concentration Pathways) which show how the climate will likely change by 2100, depending on different levels of counteraction [3]. When bold mitigation strategies are taken into consideration, the greenhouse gas emissions could be halved by 2050 with a maximum temperature increase of 2 • C, while, with a "business-as-usual" approach, the CO 2 in the photovoltaic systems (BIPV) that not only generate electrical energy but also behave like skin for the buildings [28]. These technologies have the potential to become a source of income for the buildings, even without subsidies, due to the increasing efficiency and decreasing costs of PV systems [29]. Smart management through building control and automation systems is also key, as demonstrated by the introduction of the Smart Readiness Indicator with the latest revision of the EPBD in July 2018 [30]. This applies to any technical systems [31] and to ventilation strategies [32]. Innovative lines of research are further looking into refrigerant-free cooling appliances based on caloric materials [33] and year-round passive daytime radiative cooling [34].
Generally speaking, whatever the specific strategy, the results emphasize that a one-fits-all recipe for nearly-zero energy buildings does not exist [26]. Climate, among other factors, calls for the resolution of diversified optimization problems [35,36]. Climate is a spatial and temporal variable. Here, we focus on the temporal variability by challenging established nZEB paradigms in the context of increased global warming.
Concerning the method, the need to reduce global energy consumption and CO 2 emissions has induced the European Committee for Standardization (CEN) to provide an hourly dynamic calculation method that allows buildings' consumptions to be assessed in a more realistic and detailed evaluation, especially during the summer season [37]. This method, described in EN ISO 52016−1:2017 [38], replaces the one described in ISO 13790:2008 [39] by introducing a new methodology to calculate energy needs for heating and cooling, on both an hourly and a monthly basis.
The need for knowledge on future scenarios emerges, specifically about a better understanding of the most effective energy retrofit strategies for existing nZEB buildings in the Mediterranean climate, to guide future legislative amendments, and the identification of which mitigation policies will be most appropriate for new buildings to limit both CO 2 production and global energy consumption.
Therefore, this study aims at assessing the impact of climate change on (i) the heating and cooling consumption of an nZEB multi-family house, located in Rome and designed according to the most recent Italian regulations [40], and (ii) the level of comfort achieved indoors.
The EURO-CORDEX5 [41] models combined with the ERA-Interim/UrbClim model [42,43], used for predicting future scenarios according to the Fifth Assessment Report (AR5) by the Intergovernmental Panel on Climate Change (IPCC), shows that Rome is the place with the highest temperature increase in Italy; therefore, this city has been chosen as a case study for the present study.
To this end, hourly dynamics simulations were performed in TRNSYS, which is a well-established building dynamic simulation software worldwide, capable of fine assessments of both energy and comfort levels [44], thanks to a vast variety of components that can be implemented in different models in order to simulate a wide range of simple to complex systems [45]. Its visual interface, which implements a component-based approach and the possibility of the addition of new mathematical sub-models, motivates its use for building energy simulation (BES) [46].
The rest of the paper is structured as follows: in Section 2, we present methods; in Section 3, we show the results of the simulations and, in Section 4, we draw conclusions and discuss future work.

Methods
In this work, the dynamic simulation software TRNSYS [47] was chosen to assess the thermal behavior of a building located in Rome classified as a nearly-zero energy building (nZEB) according to the Italian regulation enforced on an hourly basis. The focus is on the effects of the expected climate change and urban heat island exacerbation by 2050, both in terms of energy needs and comfort.
To this end, hourly dynamic simulations were performed in TRNSYS according to the following procedure: (i) the production of current and 2050 meteorological input files, (ii) simulations and evaluation of the year-round energy consumption, assuming an infinite power system for heating and cooling, (iii) simulations and evaluation of the free-floating operative temperature assuming no heating/cooling systems, (iv) application of the adaptive comfort [48] theory to rate the quality of the indoor environments and (v) evaluation of the long-term thermal comfort [49].

Climate and Geographical Data
The hourly climate data of the typical year adopted in this study refer to both the current climate conditions and those in 2050. About future scenarios, two Representative Concentration Pathways (RCPs) were analyzed, which refer to climate change projections developed in the AR5 report by IPCC.
The two future projections (RCP8.5 and RCP4.5) compute the expected temperature rise as a function of different greenhouse gas emissions determined by both human activities and different mitigation activities implemented by local policies: • RCP8.5 represents a 'business-as-usual' approach, which considers that atmospheric concentrations of CO 2 will triple or quadruple by 2100 compared to pre-industrial levels, thus increasing the air temperature by about 4 • C; • RCP4.5 contemplates control measure to curtail the greenhouse emissions, assuming a trend reversal (decrease below current levels) by 2070. Consequently, the atmospheric concentrations of CO 2 will be about twice as high as pre-industrial levels by 2100. The increase in air temperature will be capped at 2 • C.
The current climatic data are developed in a data set of temperatures and solar irradiances measured between 2000 and 2009 and 1991 and 2009, respectively, while future scenarios are simulated by Meteonorm 7.3 through a set of EURO-CORDEX5 [41] models combined with the ERA-Interim/UrbClim model [42,43]. The evaluation of the increase in solar radiation is based on the IPCC AR4 A2 models for CPR 8.5 and A1B for CPR 4.5 [51].
In detail, EURO-CORDEX is the European branch of the CORDEX initiative and produces ensemble climate simulations based on multiple dynamical and empirical-statistical downscaling models forced by multiple global climate models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) [41]. The major aims of the CORDEX initiative are to provide a coordinated model evaluation framework, a climate projection framework, and an interface to the applicants of the climate simulations of climate change impact, adaptation, and mitigation studies [52].
The ERA-Interim/UrbClim model allows you to evaluate the intensity of urban heat of any city in the world with a resolution up to 100 m, providing the necessary input data from international satellite-based land cover, vegetation and soil sealing data and meteorological databases [42,43].
In order to perform simulations through TRNSYS, the following input data were considered: (i) latitude, (ii) dry-bulb external temperature, (iii) relative humidity, and (iv) solar global irradiance on the horizontal plane. Data were provided by the WMO weather station #162390 located in Rome/Ciampino. At 41.8 • N (latitude), 12.58 • E (longitude) and 131 m above sea level (altitude), the location features a temperate climate, specifically Csa (Mediterranean hot summer climate), in accordance with the Köppen-Geiger climate classification system.

TRNSYS Model Simulation
Outdoor air temperature, relative humidity and global solar irradiance on the horizontal plane for each of the three different scenarios described in Section 2.1 were provided as inputs through TRNSYS component Type 9 (data reader). Then, these data were transferred to Type 16c, a solar radiation processor which implements the Reindl and Perez 1990 models to obtain beam and diffuse solar irradiances on horizontal and tilted surfaces, respectively. The virtual model of the nZEB case-study was designed through a SketchUp [53] plugin for creating the multi-zone building envelope (TRNSYS3d) [54], and then was provided as input in TRNSYS [47]. Each room of the building was defined by creating a thermal zone for each of them. An interface for the detailed TRNSYS multi-zone building (TRNBuild/Type56) was used for modeling walls, gains and ventilation profiles. In order to comply with the Italian regulation for nZEBs [40], overhangs were integrated to reduce solar gains; therefore, Type 34 was used. An hourly basis zone-by-zone calculation was performed, the results of which in terms of energy needs and indoor operative temperatures were plotted via Type 65. Figure  1 displays the flowchart.

Case-Study
The studied building is a new 3-floor residential nZEB, composed of two apartments per floor. Each apartment has an open space that combines kitchen and living room, a double bedroom, a single bedroom and a bathroom, as shown in Figure 2. Table 1 summarizes the geometrical features of each apartment.  In order to comply with the Italian nZEB requirements the building envelope is made of 2 cm of external plaster, 35 cm of pored bricks, 4 cm of external insulation, and 1 cm of internal render. The 24 cm hollow-core concrete slab of the flat and non-walkable roof is well insulated with 10 cm of expanded polystyrene. The stairwell wall is made of two layers of 12 cm hollow bricks separated by 6 cm of insulation and external plaster, while the partition wall between the two apartments is a cavity wall with two layers of 12 cm of hollow bricks and external plaster. The windows are double glazed with low-emissivity glass. The horizontal overhangs on the south-exposed facades and proper shading devices on all the windows are modeled in order to provide sufficient sunlight control and the reduction of solar gains. The thermo-physical parameters of the main building envelope elements,

Case-Study
The studied building is a new 3-floor residential nZEB, composed of two apartments per floor. Each apartment has an open space that combines kitchen and living room, a double bedroom, a single bedroom and a bathroom, as shown in Figure 2. Table 1 summarizes the geometrical features of each apartment. which in terms of energy needs and indoor operative temperatures were plotted via Type 65. Figure  1 displays the flowchart.

Case-Study
The studied building is a new 3-floor residential nZEB, composed of two apartments per floor. Each apartment has an open space that combines kitchen and living room, a double bedroom, a single bedroom and a bathroom, as shown in Figure 2. Table 1 summarizes the geometrical features of each apartment.  In order to comply with the Italian nZEB requirements the building envelope is made of 2 cm of external plaster, 35 cm of pored bricks, 4 cm of external insulation, and 1 cm of internal render. The 24 cm hollow-core concrete slab of the flat and non-walkable roof is well insulated with 10 cm of expanded polystyrene. The stairwell wall is made of two layers of 12 cm hollow bricks separated by 6 cm of insulation and external plaster, while the partition wall between the two apartments is a cavity wall with two layers of 12 cm of hollow bricks and external plaster. The windows are double glazed with low-emissivity glass. The horizontal overhangs on the south-exposed facades and proper shading devices on all the windows are modeled in order to provide sufficient sunlight control and the reduction of solar gains. The thermo-physical parameters of the main building envelope elements,  In order to comply with the Italian nZEB requirements the building envelope is made of 2 cm of external plaster, 35 cm of pored bricks, 4 cm of external insulation, and 1 cm of internal render. The 24 cm hollow-core concrete slab of the flat and non-walkable roof is well insulated with 10 cm of expanded polystyrene. The stairwell wall is made of two layers of 12 cm hollow bricks separated by 6 cm of insulation and external plaster, while the partition wall between the two apartments is a cavity wall with two layers of 12 cm of hollow bricks and external plaster. The windows are double glazed with low-emissivity glass. The horizontal overhangs on the south-exposed facades and proper shading devices on all the windows are modeled in order to provide sufficient sunlight control and the reduction of solar gains. The thermo-physical parameters of the main building envelope elements, previously described, are given in Table 2, while in Table 3 the nZEB requirements for the Italian territory are described. According to the Italian technical specifications in UNI/TS 11300 [55], the following assumptions were made: usage profiles 24/24 h, internal heat gains rate at 5.72 W/m 2 and ventilation rate at 0.5 h −1 .
Two different models of the building were created: • The first one was equipped with an infinite power system for heating and cooling, to evaluate the energy consumption. In this case, the operative temperature, used to control the system, was set to 20 • C in winter and 26 • C in summer, in accordance with the type of building and categories identified in UNI EN 16798-1-Annex A [48], assuming category II as the reference; • The second one was equipped without any heating and cooling system, to allow free-floating operative temperature for the assessment of thermal comfort. The level of thermal comfort was analyzed through two different approaches. Firstly, we applied the adaptive method, according to which the operative temperatures are compared to the indoor operative temperature ranges defined in UNI EN 16798-1-Annex A [48] for buildings without mechanical cooling systems to identify different comfort levels as a function of the outdoor running mean temperature, calculated as follows: In Equation (1), α is a constant between 0 and 1 (recommended value is 0.8) and θ ed-i is the daily mean outdoor air temperature on the i-th previous day. Secondly, we calculated a newly proposed long-term comfort index (Equation (2)). It outperformed the other 6 types of existing indices (23 total) found in the standards and 5 types of new indices (36 total) for comfort assessments in the long run in a recent comparative study [49]. The correlation with the long-term thermal satisfaction of building occupants was based on continuous thermal comfort measurements and post-occupancy evaluation surveys, but in air-conditioned office buildings. Despite the different casuistry, this combined index proved better at taking into account the triggers behind long-term comfort, namely the pronounced excursions beyond some acceptable temperature ranges and the larger variations in daily temperature than the average experience over time [49]. As such, it is worth investigating its expected variation in future scenarios. This combined, normalized percent index is calculated as follows: index = ( % To outside specified ranges + %To daily range > a threshold 2 ) In this study, the hourly operative temperature range was set between 20 • C and 26 • C, in accordance with the type of building and categories identified in UNI EN 16798-1-Annex A [48], assuming category II as reference, and the daily threshold was set to the 80th percentile of the entire simulated database.

Results
In this section the modeling results are illustrated and discussed. They refer to the two apartments at the second floor of the building described in Section 2.3. Considering the dispersing surface of the roof, the second floor is the most responsive to external climatic variations. For this reason the results obtained for this floor are evaluated and discussed.
In Section 3.1, the results in terms of kWh obtained from the modeling of the building in the presence of a heating/cooling system are analyzed, while in Section 3.2 the results in terms of comfort obtained with the modeling of the building in free-floating are discussed.

Simulations to Assess the Energy Consumption
Based on dynamic hourly simulations, the average annual increases in outdoor temperatures (+3.4 • C under RCP4.5 and +3.9 • C under RCP8.5) were found to enhance the global energy consumption from 5079 kWh × year (in 2020) to 6196 kWh × year (+18.0%), and 6248 kWh × year (+18.7%) by 2050, respectively. The breakdown in heating and cooling needs revealed a divergent trend (see Table 4): the heating energy needs showed a decrease in the order of 129.5% (RCP4.5) and 185.8% (RCP8.5), whereas the cooling energy needs increased by 47.8% (RCP4.5) and 50.3% (RCP8.5). Table 4. Air temperature (T air ), relative humidity (UR), horizontal global solar irradiance (I g,H ) , heating energy needs (ϕ h,tot ), cooling energy needs (ϕ c,tot ), peak power for heating (ϕ h,max ) and peak power for cooling (ϕ c,max ) for the three analyzed scenarios. This trend was reflected in terms of maximum hourly peak power too-that of the heating system dropped from 2.5 kW (2020) to 1.9 kW (RCP8.5), while that of the cooling system rose from 3.0 kW (2020) to 4.0 kW (RCP8.5).

Year
By performing a frequency distribution analysis of hourly energy demand (see Figure 3), an increase of about 6% was shown for the 1-2 kWh (light blue) and 2-3 kWh (blue) cooling ranges for both future scenarios, while a decrease of about 8% was shown for the 1-2 kWh (orange) heating range.
Moreover, for both future scenarios, the 2-3 kWh (red) heating range disappeared and the 3-4 kWh (black) cooling range was introduced due to an increase in outdoor temperature and solar irradiances. Only the heating range 0-1 kWh (yellow) remained unchanged among the different scenarios. By performing a frequency distribution analysis of hourly energy demand (see Figure 3), an increase of about 6% was shown for the 1-2 kWh (light blue) and 2-3 kWh (blue) cooling ranges for both future scenarios, while a decrease of about 8% was shown for the 1-2 kWh (orange) heating range. Moreover, for both future scenarios, the 2-3 kWh (red) heating range disappeared and the 3-4 kWh (black) cooling range was introduced due to an increase in outdoor temperature and solar irradiances. Only the heating range 0-1 kWh (yellow) remained unchanged among the different scenarios.  In addition to the increases in peak power and cooling energy demand, the number of hours over which the cooling system was active increased too-from 3033 h in 2020 to 3983 h (+23.8%) and 4022 h (+24.6%) under RCP4.5 and RCP8.5 scenarios, respectively. Indeed, the length of the cooling season (calculated from the first to the last day the system was on) varied for the different scenarios-174 days in 2020, 224 in 2050 (RCP4.5) and 208 in 2050 (RCP8.5). At the same time, the percentage of hours over which the conditioning system stayed off (T air ≤ 26 • C) was found to be higher for RCP4.5 than RCP8.5 (25.9% and 19.5%, respectively). Therefore, the number of hours the cooling system was switched on for in the two future scenarios was similar-3983 h for RCP4.5 and 4023 h for RCP8.5.
In combination with the increase in temperature, the model ERA-Interim/UrbClim [42] provides also an increase in global solar irradiance on the horizontal plane (I g,H ), and thus on tilted and oriented surfaces. As shown in Figure 4, an evident increase in global solar irradiance was observed for both future scenarios compared to 2020, resulting in an increase in solar gains, especially during the months of January, October, November and December. The increase in solar irradiance in the summer months, although greater than in the winter months, did not lead to an evident increase in internal gains. This is because, in order to design an nZEB building, Italian requirements [40] force the use of mobile shading throughout the summer period.
The average annual percentage of global solar irradiance on the horizontal plane was found to increase (see Table 5) by about 19.7% and 15.5% for the cooling period (May-October) and the heating period (November-April), respectively, for both future scenarios. Figure 4 and Table 5 show that the average monthly solar gains do not increase in summer between 2020 and 2050, due to movable shading devices operating in summertime only. Assuming equal internal gains and ventilation rates, the increase in average summer energy consumption, accounting for 249.3% (+2463.4 kWh) and 293.5% (+2723.5 kWh) under RCP4.5 and RCP8.5, respectively, was ascribable to the corresponding average temperature increase (2.7 • C and 3.2 • C). Differently, the decrease in average winter consumption, accounting for 70.3% (−1350.9 kWh) and 77.0% (−1554.3 kWh) under RCP4.5 and RCP8.5, respectively, was associated with both the average increase in temperature (2.1 • C and 2.5 • C) and the increase in solar radiation.
In combination with the increase in temperature, the model ERA-Interim/UrbClim [42] provides also an increase in global solar irradiance on the horizontal plane (Ig,H), and thus on tilted and oriented surfaces. As shown in Figure 4, an evident increase in global solar irradiance was observed for both future scenarios compared to 2020, resulting in an increase in solar gains, especially during the months of January, October, November and December. The increase in solar irradiance in the summer months, although greater than in the winter months, did not lead to an evident increase in internal gains. This is because, in order to design an nZEB building, Italian requirements [40] force the use of mobile shading throughout the summer period. The average annual percentage of global solar irradiance on the horizontal plane was found to increase (see Table 5) by about 19.7% and 15.5% for the cooling period (May-October) and the heating period (November-April), respectively, for both future scenarios. Figure 4 and Table 5 show that the average monthly solar gains do not increase in summer between 2020 and 2050, due to movable shading devices operating in summertime only. Assuming equal internal gains and ventilation rates, the increase in average summer energy consumption, accounting for 249.3% (+2463.4 kWh) and 293.5% (+2723.5 kWh) under RCP4.5 and RCP8.5, respectively, was ascribable to the corresponding average temperature increase (2.7 °C and 3.2 °C). Differently, the decrease in average winter consumption, accounting for 70.3% (−1350.9 kWh) and 77.0% (−1554.3 kWh) under RCP4.5 and RCP8.5, respectively, was associated with both the average increase in temperature (2.1 °C and 2.5 °C) and the increase in solar radiation.  By performing a non-linear regression analysis between daily heating/cooling needs and outdoor air temperatures (see Figure 5), we observed how future distributions gradually shift towards higher temperatures and consumptions-from a maximum cooling demand of 40.4 kWh/day in 2020 to 64.3 kWh/day under RCP4.5 and 68.0 kWh/day under RCP8.5. Precisely for this reason, the distribution in graphs b and c of Figure 5 loses symmetry, and the density of the points in the left side (heating needs) decreases. Figure 5 also shows that the minima (ϕh/c = 0) of the trendlines for future scenarios (graph b and c) narrow down and get slightly shifted towards higher temperatures (to 16.0 • C and 16.5 • C, respectively), compared to the current scenario (graph a) where 0 ϕh/c conditions occur between 14.5 • C and 16.6 • C. outdoor air temperatures (see Figure 5), we observed how future distributions gradually shift towards higher temperatures and consumptions-from a maximum cooling demand of 40.4 kWh/day in 2020 to 64.3 kWh/day under RCP4.5 and 68.0 kWh/day under RCP8.5. Precisely for this reason, the distribution in graphs b and c of Figure 5 loses symmetry, and the density of the points in the left side (heating needs) decreases. Figure 5 also shows that the minima (φh/c = 0) of the trendlines for future scenarios (graph b and c) narrow down and get slightly shifted towards higher temperatures (to 16.0 °C and 16.5 °C, respectively), compared to the current scenario (graph a) where 0 φh/c conditions occur between 14.5 °C and 16.6 °C.

Simulations to Assess Thermal Comfort
In this section, the two approaches employed for the assessment of indoor thermal comfort are sequentially presented: (i) the adaptive method and (ii) the combined long-term index proposed by Li et al. [49].
The assessment of adaptive comfort can be performed only on buildings used mainly for human occupancy engaged in sedentary activities (e.g., residential buildings and offices), where thermal conditions can be regulated by occupants through accessing operable windows and adapting their behavior. On the other hand, the new combined index considers the long-term thermal satisfaction in response to pronounced excursions and daily variability.
The analysis was performed for two rooms (see Figure 6), one due southeast (living room) and one due northwest (double bedroom), to account for different exposures and thus different solar gains. The evaluation focused on the cooling periods as identified in Section 3.1 (174 days in 2020, 224 in 2050 (RCP4.5) and 208 in 2050 (RCP8.5)); therefore, the assessment of thermal comfort is performed for each climate scenario over different periods, assuming the absence of cooling systems.

Simulations to Assess Thermal Comfort
In this section, the two approaches employed for the assessment of indoor thermal comfort are sequentially presented: (i) the adaptive method and (ii) the combined long-term index proposed by Li et al. [49].
The assessment of adaptive comfort can be performed only on buildings used mainly for human occupancy engaged in sedentary activities (e.g., residential buildings and offices), where thermal conditions can be regulated by occupants through accessing operable windows and adapting their behavior. On the other hand, the new combined index considers the long-term thermal satisfaction in response to pronounced excursions and daily variability.
The analysis was performed for two rooms (see Figure 6), one due southeast (living room) and one due northwest (double bedroom), to account for different exposures and thus different solar gains. The evaluation focused on the cooling periods as identified in Section 3.1 (174 days in 2020, 224 in 2050 (RCP4.5) and 208 in 2050 (RCP8.5)); therefore, the assessment of thermal comfort is performed for each climate scenario over different periods, assuming the absence of cooling systems. Regarding the adaptive method, as shown in Table 6, the living room provided comfort conditions (operative temperature cooling set point of 26 °C) for shorter time periods than the double bedroom, regardless of the scenario. Assuming category II as a reference, the percent decrease between 2020 and 2050 RCP4.5 reached 4.1% and 3.7% for the two considered rooms, against 9.0% for both rooms under RCP8.5. Results for the combined index are shown in Table 7. The double bedroom experienced percent increases in discomfort events of 0.9% and 3.6% between 2020 and 2050 for RCP4.5 and RCP8.5, respectively. In contrast, for the living room the variation stayed below 1% under both the considered scenarios.  Regarding the adaptive method, as shown in Table 6, the living room provided comfort conditions (operative temperature cooling set point of 26 • C) for shorter time periods than the double bedroom, regardless of the scenario. Assuming category II as a reference, the percent decrease between 2020 and 2050 RCP4.5 reached 4.1% and 3.7% for the two considered rooms, against 9.0% for both rooms under RCP8.5. Results for the combined index are shown in Table 7. The double bedroom experienced percent increases in discomfort events of 0.9% and 3.6% between 2020 and 2050 for RCP4.5 and RCP8.5, respectively. In contrast, for the living room the variation stayed below 1% under both the considered scenarios.

Discussion and Conclusions
The global climate is undergoing major upheavals, posing a serious risk of premature obsolescence for the current nZEB paradigms.
In accordance with the experimental study by A. Martinelli et al. [56], wherein it is shown that the Urban Heat Island phenomenon is more accentuated in the areas closer to the city center, the ERA-Interim/CliUrbm predictive model [42,43] used in this study identified Rome as the Italian city that will undergo the greatest temperature increase by 2050. For this reason, the impact of climate change on the energy needs and indoor comfort of an nZEB located in this location was assessed.
Based on the simulations made and detailed in Sections 2 and 3, the following considerations can be made: The respective rises of 3.4 • C and 3.9 • C by 2050 for the RCP4.5 and RCP8.5 scenarios do not lead to decreases in heating needs great enough to meet the sharp increase in cooling needs. Specifically, in fact, compared to the current energy needs, there will be an average annual increase of 1143 kWh (+22.4%). This increase is consistent with the increase obtained by Ciancio et al. [4] (+34.4%) for the building located in Rome, which is not an nZEB, but has similar transmittances to our case study.

a.
Peak electricity demand is especially worrisome since it is usually covered by low-efficiency power plants, yet it is strongly associated with typical nZEB paradigms. In fact, while air conditioning is only a fraction of all building energy uses, it is the primary driver of peak electricity demand [57]. Efficiently curbing the air conditioning needs by targeting a resilient nZEB design will be key in the future. b.
According to the ERA-Interim/UrbClim model, by 2050, not only will temperatures rise, but a 19.7% increase in global solar irradiance on the horizontal plane is also to be expected during the summer months, thus triggering higher solar gains. The increase in the solar irradiance, often underestimated by the models [58,59], is consistent with the study conducted by M. Wild et al. [60], which predicts a decrease in solar radiation (clear-sky condition) for many regions of the world, except for parts of China and Europe. In our case study, the increase in the solar irradiance does not imply an increase in the solar contributions inside the building because the national nZEB regulation foresees the use of mobile shading devices for the whole summer period. These systems appear rather effective, since the increase in solar gains is negligible in summer, but becomes evident in winter (particularly in January, October, November and December). c.
By performing simulations in the absence of cooling systems, 6.2% and 5.1% reductions in the hours of adaptive comfort are determined under the RCP4.5 (2050) and RCP8.5 (2050) scenarios, respectively, out of the concerted actions of temperature and solar gains. The results of the newly proposed combined index for long-term comfort assessments revealed a milder future penalty. The index estimates how the level of occupant adaptation and sensitivity to variation would be affected in the future. It was demonstrated that the comfort implications of the pronounced excursions and large variations in daily temperature will be marginal under both scenarios, with greater influence on northwest rather than southeast oriented thermal zones, likely owing to the effect of the combination of higher temperatures and higher solar irradiation in levelling out the daily swings.
In conclusion, this study adds to the current body of knowledge on the preservation of nZEB performance in the future by performing hourly dynamic simulations on a reference building in Rome, modeled in accordance with the latest legislations. It quantifies potential changes in terms of energy and comfort levels and provides useful recommendations to legislators on building standards, both for the design of new nZEBs, e.g., the presence of solar shading devices, and for the renovation of the existing building stock. Further analysis may target climate dependencies and may include technological variants.