Building Thermal and Energy Performance of Subtropical Terraced Houses under Future Climate Uncertainty

: Due to global temperature increases, terraced house (TH) residents face a threat to their health due to poor indoor thermal environments. As buildings are constructed by low-income residents without professional guidance, this study aims to investigate the indoor thermal comfort and energy resilience of THs under the future climate and determine the optimal passive design strategies for construction and retroﬁtting. By exploring the effects of building envelope structures, adjusting the window-to-wall ratio (WWR) and designing shading devices, EnergyPlus version 22.0 was used to optimize the thermal environment and cooling load of THs throughout their life cycle under future climate uncertainties. Unimproved THs will experience overheating for nearly 90% of the hours in a year and the cooling load will exceed 60,000 kWh by 2100 under the Representative Concentration Pathways (RCP) 8.5 scenario. In contrast, optimization and improvements resulted in a 17.3% reduction in indoor cooling load by increasing shading devices and the WWR, and using building envelope structures with moderate thermal insulation. This study can guide TH design and renovation, signiﬁcantly reducing indoor cooling load and enabling residents to better use active cooling to combat future overheating environments.


Background
The first part of the Intergovernmental Panel on Climate Change's (IPCC) Sixth Assessment Report (AR6) mentions that global temperature will rise by from 1.5 • C to 2 • C in the coming decades [1]. Global warming is projected to bring about more severe heatwave events and urban heat island effects, along with their synergistic effects [2]. Furthermore, the probability of cardiovascular, respiratory, and cerebrovascular patient mortality will increase in high-temperature environments [3,4]. Therefore, thermal conditions and energy demand play crucial roles in the survival of residential buildings in climate change, especially under future scenarios with more intense climate changes, which necessitates greater attention to assessments of building performance.
Improving the thermal comfort and building energy demands of residential buildings is particularly important as they serve as crucial living spaces for humans. Residential buildings lacking sufficient resilience will pose higher heat risks, more extreme indoor temperatures, and longer durations of overheating in higher-outdoor-temperature environments [5,6]. Furthermore, humidity can significantly impact indoor thermal comfort and energy demand [7]. These factors can result in further health issues and even fatalities for occupants [8]. Excessive indoor temperatures can lead to heatstroke and heat exhaustion [9], as well as cardiovascular, kidney, and respiratory diseases [10][11][12], and may even have psychological impacts that contribute to incidents of suicide [13]. The energy demands of residential buildings also increases with rising temperatures, and reports indicate that the surge in building energy demands during high-temperature extreme-weather periods poses significant challenges and risks to the electricity sector [14]. The additional energy consumed during heatwaves also contributes to increased carbon emissions [15]. Therefore, improving the resilience of residential buildings to future climate change is of great significance in reducing the probability of heat risks, mitigating indoor overheating, safeguarding the physical and mental well-being of occupants, alleviating stress on urban power systems, and reducing carbon emissions.
In China, the rural population are allowed to build their own houses on homestead. Terraced houses (THs) are designed and built by the building owners without technical guidance, and are common outside urban areas. Improving the energy performance and thermal comfort of naturally ventilated THs is important to enhance the living quality and ensure better health for residents, especially in hot-humid areas under climate change. Due to a lack of relevant architectural expertise among TH owners and their tight budget, the energy performance of poor insulated THs is low. On the other hand, despite prolonged exposure to overheating conditions, rural residents with a high tolerance on high temperatures often choose to reduce their energy expenses by refraining from or minimizing the use of air conditioning [16]. However, the two objectives of energy performance and thermal comfort under natural ventilation are somehow conflicting [17][18][19]. It is necessary to offer a quantitative assessment of the building performance of THs in the context of global warming to help residents to make wise choices regarding building retrofitting.

Literature Review
Climate change trends are expected to raise atmospheric temperatures, leading to various challenges for residential buildings in overheated environments. Researchers use Global Climate Models (GCMs) to obtain monthly data of various meteorological parameters and morph these data using statistical methods to generate hourly data. Sailor et al. [20] simulated that indoor temperatures without air conditioning would exceed 40 • C in two warm US cities (Chicago, Illinois, and Houston, Texas) in 2050. Baba et al. [21] simulated a school in Montreal, Canada, discovering that, without mitigation measures, classroom overheating hours would increase to 333 and 437 h by 2044 and 2090, respectively. Baglivo et al. [22] found that, from 2020 to 2080, the hours worked in temperatures higher than 26 • C in apartments in Lecce, Italy, had increased by around 16% while the external temperature increased by 11%. Andric et al. [23] further found that the hours with temperature exceeding 38 • C in Qatar would increase by 17% and 21% for 2050 and 2080, respectively. Muñoz González et al. [24] simulated a scenario in southern Europe, with discomfort indoor projected to increase by 20-30% in warm months in 2050. Ciancio et al. [25] simulated several European cities, such as Clermont-Ferrand, Plovdiv, Granada, and Salamanca, finding that temperature would increase between 1 • C and 5 • C in July and August in 2050 and between 2 • C and 9 • C in 2080. Pérez-Andreu et al. [26] found that the indoor temperatures of residential buildings in Valencia, Spain, along the Mediterranean coast are projected to increase by 3.60 • C and 5.33 • C under the CNRM-CM, and MPI-ESM-LR circulation scenarios, respectively, from 2096 to 2100. Flores-Larsen et al. [27] found that, from 2020 to 2080, temperatures in Santa Rosa, Mendoza, Córdoba, and Orlán would rise by 2.2 • C, 2.7 • C, 2.8 • C, and 3.8 • C, respectively (annual average). Overall, under future climate change, levels of overheating, maximum temperatures, and overheating risks in residential buildings will all increase, and the challenges faced by residents deserve attention.
As indoor thermal comfort is expected to decrease under future climate change, researchers have studied energy demands to mitigate this discomfort. Velashjerdi Farahani et al. [28] studied Helsinki in southern Finland and found that cooling load was projected to increase by 47% under the future average climate in 2050, and up to 128% under extreme climate conditions. Zou et al. [29] found that, under the SSP3-7.0 scenario, residential energy use in humid regions of China would increase by 44.5% to 59.1% from 2075 to 2099. Hosseini et al. [30] simulated the southeastern environment of Sweden, finding that the peak cooling load load of typical years from 2040 to 2069 increased by 210% compared to 2010-2039, and by 290% from 2070 to 2099. Bell et al. [31] studied Köppen-Geiger climate zones and found that, under extreme weather conditions in 2050, the peak cooling demand increased by 24-35%, and unsatisfied cooling hours increased by 41-189%. Jafarpur et al. [32] predicted that, in 2056-2075, Quebec City, Toronto, and Vancouver in Canada are expected to increase their cooling load by 34.6% (7.9 kW/m 2 ), 32.2% (8.6 kW/m 2 ), and 27.8% (6.5 kW/m 2 ), respectively, using EnergyPlus version 22.0. Moazami et al. [33] studied a district in Geneva, finding that peak electricity demand increased by 4.0%, 7.6%, and 16.8% for 2011-2040, 2041-2070, and 2071-2100, respectively, under future extreme climate conditions. In the simulation of Seville, southwestern Spain, Sánchez-García et al. [34] found that the cooling load demand increased from 14.6% in 2020 to 57.5% in 2080. Xiang Li et al. [35] found that Switzerland's service industry's spatial cooling load demand would increase by 400%, 500%, and 600% by 2050 under Representative Concentration Pathway (RCP) 2.6, RCP 4.5, and RCP 8.5, respectively. Ramon et al. [36] found that, in Belgium, the number of cooling days increased from 167 days to 401 days, an increase of 2.4 times between 1976-2004 and 2070-2098. Spinoni et al. [37] simulated the last 10 years of the 21st century in Europe and found that the annual number of cooling days under RCP 4.5 increased by approximately 45% on average, while, under RCP 8.5, the number increased by 170%. In conclusion, under gradually increasing temperatures in the future, residents will consume more energy to cool down overheating indoor environments, which will result in higher energy costs, indicating that low-income groups will face various health risks they cannot afford, and there will be a significant challenge for cities' energy load.
Facing the dual challenge of the increasing risk of building overheating and energy demand in the context of future temperature rises, some researchers are attempting to save energy by improving the thermal performance of building envelopes; however, this may lead to disadvantages in the thermal environment under natural ventilation. Zhang et al. [38] simulated the environment in Shanghai and concluded that adding an insulation board could increase the indoor winter temperature by 7.8% and decrease the summer temperature by 7.5%, resulting in a certain insulation and heat-preservation effect. Moreles et al. [39] improved the thermal performance of building envelopes using PCM walls with thicknesses of 6 mm and 18 mm, and found that the maximum energysaving effect in air-conditioning environments increased by 4% and 29%, respectively. Alshuraiaan et al. [40] discovered that using phase change materials (PCMs) to improve building thermal performance can reduce heat flow into the building by more than 50%. Wu et al. [41] found that using multi-layer assemblies of phase change materials and bio-based concrete to improve building thermal performance in summer scenarios could reduce the maximum heat load by 8.2%, temperature fluctuations by 46.3%, and local vapor pressure fluctuations by 43.7%, which could enhance indoor thermal comfort. Aslani et al. [42] found that buildings with DSW R-50 walls in the building envelope consume 75% less energy than conventional buildings. Tarek et al. [43] found that replacing traditional bricks with FS geo-polymeric bricks in the city of New Aswan, Egypt, reduced the indoor cooling demand from 32.41103 kWh to 30.64103 kWh, a reduction of about 5.7% compared to the use of traditional bricks in summer peak periods. Additionally, other geo-polymeric bricks can reduce building energy demand by from 5.7% to 14.9%. Dora et al. [44] improved the thermal performance of buildings by using large encapsulated PCM materials, which reduced the total peak heat flux by 27.32%, lowered the cooling load by 38.76%, and saved 28.31 rupees/day (about 0.40 USD/day) in electricity costs. Improving the thermal performance of the building envelope has a certain energy-saving effect, but also has some impact on the thermal environment. Mourid et al. [45] found that the use of PCM to improve building thermal performance can lead to a temperature increase of about 6 • C indoors at night and significantly increase the surface thermal quality of the room. Al-Yasiri et al. [46] found, in their experiment in Iraq, that the average indoor air temperature in PCM rooms increased slightly compared to another identical non-PCM room, regardless of natural ventilation. Yang et al. [47] optimized energy demand by using genetic algorithms to develop improvement strategies. Yu et al. [48] found that improving the building envelope's insulation could reduce heat transfer from the outdoor environment to the indoor environment, which is advantageous for energy conservation during the day but disadvantageous for heat dissipation at night. Al-Yasiri et al. [49] found that using PCM to improve building thermal performance in hot climates results in a higher operating temperature in PCM rooms than in non-PCM rooms due to the PCM absorbing heat during the day and releasing heat uncontrollably to the indoor environment in the evening. Zou et al. [50] effectively improved indoor thermal comfort by improving the passive survivability of four indigenous buildings in humid and hot areas under future climate scenarios.
Therefore, when considering upgrading building resilience, researchers need to conduct more comprehensive studies to balance the impact of these improvement measures on energy efficiency and the thermal environment, and find the best solution that can meet the needs of different users after a comprehensive assessment of both factors. However, for low-income populations, there is no definitive conclusion regarding the priorities of improving the thermal environment under natural ventilation and optimizing energy demand. Therefore, this study aims to explore the best passive design strategies for improving indoor thermal comfort under natural ventilation and reducing energy demand under air-conditioning, finding a balance between the two factors under future climates using sensitivity analyses, and providing references for the construction and renovation of thermal habitats.

Aim of This Study
To our knowledge, it is expected that indoor temperatures around the world will increase under future climate change trends, leading to a reduction in indoor thermal comfort [30]. This increase in temperatures will also cause a rise in building energy demand, particularly due to the increased use of cooling equipment [51]. Additionally, there are limitations to the use of enclosure structures to improve the indoor thermal environment under natural ventilation [52]. Previous research on optimizing indoor thermal environment and energy demand in Guangzhou under future climates is scarce. This article aims to study the following: (1) changes in outdoor temperature in the Guangzhou region under future climate conditions; (2) evaluation of indoor thermal comfort and energy demand of THs in Guangzhou under 24 different future weather conditions; (3) improvement potential of energy performance and thermal comfort for THs under future climate by retrofitting.

Methodology
This article presents a research process illustrated in Figure 1. The building performance simulation engine EnergyPlus [53], a popular choice for building performance calculations, was used for the building performance calculations. The Ladybug tools [54], a visualization programming interface, were utilized for modeling. Regarding climate data, TMY data and meteorological data for the years 2030-2100 under three different RCP scenarios were employed as inputs. The TMY dataset was jointly developed by the China Meteorological Administration Meteorological Information Center and Tsinghua University, spanning from 1994 to 2003 [55].
To assess the degree of overheating both outdoors and indoors, the Apparent Warmth Degree (AWD) and Indoor Overheating Hours (IOH) metrics were referenced. Considering that the research targets contained a region with hot summers and mild winters, the primary energy demand of local buildings was focused on cooling load, while heating load was not considered. As this study aims to optimize and retrofit TH buildings as a whole, exploring the most suitable retrofitting strategies applicable to the entire TH building cluster, it does not consider phased retrofitting. The final results are displayed through Python, presenting the average IOH and average cooling load for each room. The schedule of occupancy, lighting, and equipment is illustrated in Figure 2. The Simulation settings are shown in Table 1.  As shown in Figure 3, the study site is located in the hot-summer and warm-winter region of southern China [56], with dense distribution in the urban-rural transition zone. The region is located in the subtropics and is influenced by warm and humid air currents, resulting in hot and humid summers with abundant rainfall. On the other hand, winters are relatively mild, with rare occurrences of severe cold weather and widespread freezing and snow. As this article aims to investigate the risk of overheating, the study period is limited to the summer season from April to October in Guangzhou. As shown in Figure 4, the outdoor dry-bulb temperatures are generally higher than 26 • C, with an average temperature of 27.22 • C and maximum temperatures exceeding this. Representative Concentration Pathways (RCP) are jointly formulated by the international climate science community to assess and simulate future climate change scenarios. They provide predicted concentrations of greenhouse gases in the atmosphere under different emission scenarios and serve as input parameters for climate models. RCP2.6 represents a relatively low greenhouse gas emission scenario, requiring significant global greenhouse gas reduction measures in the mid-term to keep greenhouse gas concentrations at lower levels. RCP4.5 represents a relatively moderate greenhouse gas emission reduction scenario, predicting a moderate increase in greenhouse gas concentrations, with weaker emission reduction targets. RCP8.5 represents a relatively high greenhouse gas emission scenario, predicting a rapid increase in greenhouse gas concentrations, with very limited emission reduction measures [57]. This study explores the climate change scenarios for the three RCP environments from 2030 to 2100, considering a total of 24 potential future climates. The relevant meteorological data were obtained from the Meteonorm database.
In the 24 future climates from 2030 to 2100, the outdoor dry-bulb temperatures generally increase. In scenarios with a stronger greenhouse effect, the outdoor dry-bulb temperatures are higher, and the duration with temperatures exceeding 26 • C will be longer, indicating a greater risk of overheating and less comfortable indoor thermal environment, as shown in Figure 3. From 2030 to 2100, the average temperature from April to October is projected to increase by 2.05% in the RCP4.5 scenario and by 5.79% in the RCP8.5 scenario compared to the RCP2.6 scenario. The temperature in the period from April to October is predicted to increase at a rate of approximately 0.2% per decade in the RCP2.6 scenario, 0.86% per decade in the RCP4.5 scenario, and 1.94% per decade in the RCP8.5 scenario. The maximum temperature is projected to reach 38.4 • C in the RCP2.6 scenario, 39.8 • C in the RCP4.5 scenario, and 42.1 • C in the RCP8.5 scenario, as shown in Figure 5. THs were located in the southern part of China, characterized by hot summers and cold winters. The cooling load dominates in the summer, while heating load is generally not considered in the winter. Local residents mainly use two cooling modes indoors: natural ventilation to mitigate indoor overheating and air-conditioned with air conditioning to alleviate indoor overheating. Buildings with good insulation can save cooling energy, but this may also reduce the heat dissipation from the indoors to outdoors through natural ventilation, leading to higher overheating risks. Therefore, this study discusses IOH under natural ventilation and cooling load under an air-conditioned layout separately. For natural ventilation, we considered an all-day window opening schedule, with 50% of the window glass area being openable and 100% of window glass height being openable. For the air-conditioned area, the windows are closed.

Ambient Warmness Degree
Apparent Warmth Degree (AWD), proposed by Hamdy et al. [58], is a metric used to evaluate the severity of the outdoor temperature. It calculates the mean cooling degree hours based on a baseline temperature of 26 • C during the period when the outdoor temperature is not lower than 26 • C [59].
where T a,i is the outdoor dry bulb temperature, T b is the average temperature during the summer in Guangzhou (26 • C), N is the number of hours that T a,i is greater than 26 • C, and t i is set as 1 h for evaluation, as shown in Equation (1).

Indoor Overheating Hours (IOH)
Indoor Overheating Hours (IOH) is a metric used to evaluate non-comfortable indoor environment caused by heat. Operative temperature is a comprehensive parameter that measures the combined effect of indoor air temperature and surface temperature on human comfort. It takes into account the heat radiation experienced by the human body and its perception of the air temperature [60]. Therefore, IOH calculates the annual hours when the indoor operative temperature exceeds the comfort temperature of the region, as shown in Equations (2) and (3).
where z is the occupied zone counter, i represents the occupied working time, N occ(z) is the total working time of occupied zone z during the statistical period, t denotes time step (1 h), OT i,z represents the operative temperature of zone z at time step i, and TL com f ,i,z denotes the maximum comfortable temperature limit of zone z at time step i. According to the research conducted by Zhang et al. [61], the acceptable range of indoor temperatures for people in hot summer and warm winter regions of China is 18.0-28.5 • C. Therefore, TL com f ,i,z is set at a maximum value of 28.5 • C.

Sensitivity Analysis
Sensitivity analyses have been utilized to assess building energy demand and peak demand [61]. The aim of this study is to investigate the impact of various component parameters on building performance indicators, with IOH under natural ventilation and cooling load under air-conditioning being the primary building performance indicators. Researchers conducted a sensitivity analysis of the influence of sunshade length and window-wall ratios (WWR) for four building orientations, as well as the thermal insulation capabilities (U-values) of windows, roofs, and walls, on the building performance indicators. The sensitivity coefficient formula for the sensitivity analysis is as follows: where SC represents sensitivity coefficient, ∆OP represents the difference in energy demand or overheating hours before and after the strategy change, and ∆IP represents the difference in building performance before and after the strategy change. OP represents the average energy demand difference or overheating hours before and after the strategy change, and IP represents the average building performance before and after the strategy change. This study selected the strategy change that resulted in the greatest difference in energy demand or overheating hours to investigate the factors that have the greatest impact on indoor overheating hours and cooling load, as shown in Equation (4).

Building Geometry
As shown in Figure 6, THs distributed along the streets have four different orientation types. TH is a type of housing built by low-income individuals to meet their own housing needs without professional guidance. THs are 17.5 m long and 10 m wide, consisting of a ground floor with a shop, kitchen, dining room and bathroom, a second floor with a living room, bathroom and two bedrooms, a third floor with a bathroom and four bedrooms, and a fourth floor with a bathroom, a terrace, and two bedrooms. They also feature a central courtyard and a double staircase that runs through all four floors, as shown in Figure 7. The original TH design strategies are as shown in Table 2.

Severity of Outdoor Climate
As can be seen from Figure 8, in accordance with the assessment criterion for ambient temperature, the degree of outdoor overheating is expected to gradually climb from 2030 to 2100. For the RCP2.6 condition, it will rise at a rate of 1.23% per decade and peak at 5.09% by 2100. Under RCP4.5, the corresponding increase will be 3.02% per decade and a staggering 20.94% by 2100. If placed under the RCP8.5 condition, the degree of outdoor overheating will mount at a rate of 8.01% every decade, reaching an astounding 72.47% by the end of the century.

Sensitivity Coefficients of Various Components under Four Climates
From Figure 9, it can be seen that, under TMY conditions, the WWR has the greatest impact on the indoor overheating risk, with an average sensitivity coefficient of 0.078, while the overheating sensitivity coefficients of other design parameters are less than 0.03. In the cooling load sensitivity analysis, the sensitivity coefficients of WWR and wall U-value are significantly greater than those of other design parameters, being 0.095 and 0.087, respectively.
As can be seen from Figure 10, under the RCP2.6 scenario from 2030 to 2100, the WWR remains the factor with the greatest impact on indoor overheating risk at an average sensitivity coefficient of 0.056, while the roof U-value has the lowest sensitivity. In the cooling load sensitivity analysis, wall U-value and WWR showed a greater sensitivity, with coefficients of 0.112 and 0.086, respectively.  From Figure 11, it can be seen that under the RCP4.5 scenario from 2030 to 2100, the WWR is still the factor with the greatest impact on indoor overheating risk, at an average sensitivity coefficient of 0.048. In the cooling load sensitivity analysis, wall U-value has the greatest sensitivity coefficient at an average of 0.126, followed by the WWR at 0.096.
Under the RCP8.5 scenario from 2030 to 2100, the WWR remains the factor with the greatest impact on indoor overheating risk at an average sensitivity coefficient of 0.012. In the cooling load sensitivity analysis, wall U-value has the greatest sensitivity coefficient, at an average of 0.146, followed by the WWR at 0.089, as shown in Figure 12.  The study found that, in the TMY, RCP2.6, RCP4.5, and RCP8.5 scenarios in 2100, the WWR has the greatest impact on indoor overheating risk with a sensitivity coefficient ranging from 147.86% to more than 253.00% greater than other components in the same scenario. Meanwhile, the roof has the smallest impact on overheating risk. Additionally, the WWR and wall U-value have a much greater impact on indoor cooling load compared to other factors. In the TMY scenario, the impact of the WWR and wall U-value is almost the same, but as the RCP increases, the impact of wall U-value on cooling load becomes more significant than that of the WWR. In the RCP8.5 scenario, wall U-value has the greatest impact on cooling load, 163.2% greater than the WWR. Therefore, the WWR affects ventilation efficiency and the amount of solar radiation entering the room, while the wall U-value controls the external heat radiation entering and internal heat radiation exiting the room.

Thermal Discomfort and Energy Demand
To investigate the impact of enclosure structure performance on indoor thermal environment and cooling load, four different performance parameters were set for each component, including walls, roofs, windows, WWR, and shading, for four different orientations. The overheating and cooling load in the indoor environment from 2030 to 2100 were studied under TMY, RCP2.6, RCP4.5, and RCP8.5 scenarios.
Under the TMY scenario, the best improvement in cooling load was observed when the shading length was 0.2 m (Shading 2.0), WWR was 20% (WWR 0.2), the window U-value was 1, the roof U-value was 0.25, and the wall U-value was 1. Moreover, a longer shading length and smaller WWR, as well as lower changes in U-values for windows, roof, and wall, led to lower cooling loads in the indoor environment. As for the risk of overheating indoors, only the increase in shading length and gradual decrease in changing trend of roof U-value led to fewer overheating hours, while the opposite was observed for changes in WWR, window, and wall improvement strategies, as shown in Figure 13. Under the RCP2.6 scenario, the lowest probability of indoor overheating risk was observed when the modification strategies for wall, roof, window, WWR, and shading were 1.5, 0.25, 2, 0.5, and 2, respectively. On the other hand, the lowest cooling load in the building was achieved when the modification strategies were 1, 1.25, 2, 0.5, and 2, respectively, as shown in Figure 14. Under the RCP4.5 scenario, the lowest probability of indoor overheating risk was observed when the modification strategies for wall, roof, window, WWR, and shading were 1.5, 0.25, 2, 0.5, and 2, respectively. On the other hand, the lowest cooling load in the building was achieved when the modification strategies were 1, 0.25, 2, 0.5, and 2, respectively, as shown in Figure 15. Under the RCP8.5 scenario, the lowest probability of indoor overheating risk was observed when the modification strategies for wall, roof, window, WWR, and shading were 1.5, 0.25, 2, 0.2, and 2, respectively. On the other hand, the lowest cooling load in the building was achieved when the modification strategies were 1, 1.5, 2, 0.5, and 2, respectively, as shown in Figure 16. As shown in Figure 17, under the RCP2.6 scenario, the median of improved indoor overheating hours ranges from 57.7% to 58.6%, with a maximum of 62.51%. Under the RCP4.5 scenario, the median of overheating hours ranges from 65.5% to 66.2%, with a maximum of 76.43%. Meanwhile, under the RCP8.5 scenario, the median of overheating hours ranges from 75.2% to 75.8%, with a maximum of 90.08%. As observed in Figure 18, under the RCP2.6 scenario, the median improved indoor cooling load ranges from 36,969 kWh to 38,594 kWh, with a maximum of 42,784 kWh. Under the RCP4.5 scenario, the median cooling load ranges from 41,335 kWh to 43,141 kWh, with a maximum of 52,003 kWh. Meanwhile, under the RCP8.5 scenario, the median cooling load ranges from 48,214 kWh to 49,910 kWh, with a maximum of 70,751 kWh. Additionally, the cooling load is lower for buildings facing north or south, but higher for those facing east or west. The study found that under the RCP2.6, RCP4.5, RCP8.5, and TMY scenarios, buildings with a north-south orientation had a lower probability of indoor overheating and lower cooling loads compared to those with an east-west orientation. In scenarios with a north orientation, buildings had the lowest probability of indoor overheating and cooling loads, while those with a west orientation had the highest. To improve the situation, reducing solar radiation into the building and lowering the U-value of the building envelope to enhance the building's heat dissipation were effective in reducing the cooling load. Strategies such as increasing shading length, WWR and wall U-value were effective in reducing the probability of indoor overheating risk by reducing solar radiation and increasing natural ventilation while reducing the inflow of outdoor heat radiation.
Through simulating various improvement strategies, we selected the optimal renovation strategies for residents with different orientations and needs. Residents primarily concerned with improving the thermal environment under natural ventilation can refer to the design strategies in the IOH section, while those primarily concerned with reducing cooling load with air-conditioning can refer to the design strategies in the cooling section, as shown in Table 3. Incorporating the TMY, RCP2.6, RCP4.5, and RCP8.5 scenarios, the integration of Strategy 1, which corresponds to the south orientation, resulted in a cooling load reduction of from 12.81% to 21.11% when compared to the original strategy. Similarly, Strategy 2, aligned with the west orientation, led to cooling load reductions of from 17.3% to 23.56%. Strategy 3, representing the north orientation, yielded cooling load reductions ranging from 9.71% to 18.17%, while Strategy 4, embracing the east orientation, exhibited cooling load reductions of from 4.1% to 16.7%. This information is illustrated in Figure 19.

Discussions
The trend of rising temperatures in future is inevitable. For buildings located in humid and hot regions, the challenges of future energy efficiency and maintaining indoor comfort will be more pronounced. Terraced houses, as a predominant residential building typology in the study area, have a closely intertwined relationship between their architectural performance, local energy-saving initiatives, and residents' quality of life and well-being. In comparison to historical typical climates, Guangzhou's future temperatures are projected to experience a significant increase. The extent of the rise in cooling load for THs is contingent upon the deterioration of thermal conditions under natural ventilation, in relation to various climates. Thus, it becomes essential to explore the potential for enhancing TH building performance through retrofitting under future climates, thereby offering insights to building owners to allow for informed decisions to be made regarding architectural renovations.
In general, the various design parameters of TH, the WWR and wall U-value have the most significant impact on indoor IOH and cooling load. By adjusting these design parameters, the optimization of indoor thermal environment and cooling load can be achieved. It is worth noting that window opening size directly affects the ventilation rate and solar radiation entering the interior. Under natural ventilation, although a larger WWR value results in more solar heat gain, it removes indoor heat more efficiently by promoting air movement. As a result, this could mitigate indoor overheating to some extent [62]. On the other side, in the air-conditioned case, increased solar radiation results in a higher cooling load [63], as seen in the results presented in Section 3.3, where a larger WWR leads to higher IOH but lower cooling load, and vice versa with reduced WWR.
The uncertainty of future climate can impact the effectiveness of architectural renovations. Under natural ventilation, with future temperature increases, the impact of TH components on indoor thermal environment diminishes, and improving building envelope insulation performance alone is insufficient to enhance indoor thermal comfort. In the RCP2.6 scenario, it is possible to mitigate indoor overheating by adjusting the WWR, wall U-value, shading device length, roof U-value, and window U-value. In air-conditioned environments, measures such as reducing WWR, decreasing the envelope's U-value, and increasing the length of shading devices significantly reduce cooling loads. However, with worsening climatic conditions under the RCP4.5 and RCP8.5 scenarios, it becomes challenging to alleviate indoor overheating through design parameter modifications of building components under natural ventilation. The results indicate that any optimized design strategies maintain the IOH at around 75%. Due to the rising outdoor temperatures, even higher ventilation rates fail to effectively cool the indoors. For TH residents with high heat tolerance, the implementation of building performance upgrade and retrofit may not necessarily be beneficial, both ecologically and economically, if the frequency of air conditioning usage does not significantly increase with rising temperatures.
Conversely, significant reductions in cooling load can be achieved through retrofitting, aligning with findings from some previous studies. In an air-conditioned environment, reducing the WWR, decreasing the envelope's U-value, and increasing the length of shading devices all contribute to reducing external heat radiation. Therefore, in an environment where outdoor temperatures are gradually increasing, employing optimized design strategies can still significantly decrease indoor cooling loads. This finding aligns with relevant research conducted in other regions. Albayyaa et al. [64] achieved a 35% energy saving in buildings through retrofitting walls and floors. El-Darwish et al. [65] found that simple retrofitting strategies such as shading, window glazing, airtightness, and insulation could lead to an average energy consumption reduction of 33%. Yılmaz's research on Turkish buildings revealed that altering design strategies had a significant positive impact on energy efficiency [66]. Considering the stock of TH buildings, the overall improvement in energy efficiency for this building type is of great significance in achieving carbon emission targets, especially in scenarios with substantial temperature increases.
Given the private nature of residential properties, whether to conduct renovations and how to do so largely depend on the users' economic conditions and the expected benefits. The high initial investment in renovation is one of the main obstacles preventing residents from implementing building renovations [47]. In fact, adopting a reasonable retrofit scheme has long-term economic benefits (often with a 20-or 30-year cycle), characterized by an acceptable payback period and positive net present value [55,[67][68][69]. Additionally, initial investments in renovations can be covered through bank loans. However, most homeowners are not aware of the economic benefits of building energy-saving retrofits due to a lack of knowledge in the field. As long as a thorough calculation ensures the economic benefits of the building energy-saving retrofit plan and s explained to homeowners, their enthusiasm for renovating their homes can be motivated. Clearly, the government and professionals need to coduct a lot of work in this regard.

Conclusions
The present study examines the impact of design strategies on the natural ventilation, indoor thermal environment and energy demand of THs under scenarios of increasingly severe climate change in the coming decades. Sensitivity analysis was used to investigate the impact of individual components on the overheating hours and cooling load. Multiple improvement strategies were simulated using EnergyPlus software, v. 23.1.0, to analyze indoor thermal environment and energy performance. The optimal set of strategies was selected for each of the four building orientations.
Multiple studies indicate that climate conditions will worsen from 2030 to 2100, with outdoor temperatures increasing at a rate of 0.2% to 1.94% per decade. This will lead to severe urban heat island effects and heatwaves, with the synergistic effects between the two becoming stronger. The rising outdoor temperatures result in more adverse indoor environments and increasing risks to the psychological and physiological health of residents. The burden of indoor cooling also significantly increases, posing enormous challenges to urban power systems. Some low-income populations reside in poorly constructed houses without professional guidance, leading to an inadequate thermal performance and low energy efficiency, which results in a higher risk of overheating and an increased cooling energy demand. Therefore, improving the thermal comfort environment in THs is of great significance for the health and energy burden of residents.
Environmental warmth is used to evaluate outdoor temperature under future climate conditions. The temperature change varies significantly across different climate scenarios. Between 2030 and 2100, outdoor temperature is expected to rise dramatically. Under the RCP2.6 scenario, the temperature is projected to increase by 5.09% by 2100 relative to 2030. Under the RCP4.5 scenario, the increase will be 20.94%, and under RCP8.5, the phenomenon will be even more pronounced, with a growth of 72.47% by 2100 compared to 2030. However, the trend of rising outdoor temperature is evident in all climate scenarios.
This study simulates the impact of various building components of THs on indoor thermal environment and energy demand under future climate conditions. It investigated the sensitivity of the U-values of walls, roofs, windows, as well as the length of shading devices and WWR. The results indicate that the WWR has the greatest impact on the overheating hours, while the WWR and wall U-value have the greatest impact on the cooling load. Moreover, a larger WWR can provide better natural ventilation to cool the indoor environment, while longer shading devices have a better shading effect, reducing the solar radiation entering the room and effectively reducing the overheating hours. However, larger U-values of walls, roofs, and windows are not conducive to the dissipation of indoor heat, leading to an increase in overheating hours. Changes in the length of shading devices in the east-west direction and control of WWR in the north-south direction have some impact on the overheating hours. The optimization effect of mitigation measures on overheating hours decreases as the climate becomes more severe. Under the worst-case RCP8.5 scenario, the best optimization effect on overheating hours is only 1.29%, while it is still effective in reducing indoor cooling load, with the best optimization effect of up to 7.62%. To find the optimal combination of mitigation strategies, this study selects one set of improvement combinations for each of the four building orientations. The result indicates that the combination of THs facing south can reduce cooling load by from 12.81% to 21.11% compared with the original strategy under the TMY, RCP2.6, RCP4.5, and RCP8.5 scenarios, while the combination for the west-facing THs reduced cooling load by 17.3% to 23.56%, the combination for the north-facing THs reduced the cooling load by 9.71% to 18.17%, and the combination for the east-facing THs reduced the cooling load by 4.1% to 16.7%.
Mitigating the indoor thermal environment under natural ventilation requires a larger WWR, but also introduces more solar radiation, thereby increasing the cooling load. In contrast, mitigating the indoor thermal environment in air-conditioned THs requires active cooling measures to reduce solar radiation and lower the cooling load. It is noteworthy that THs face challenges in alleviating the IOH under future climate changes through natural ventilation alone and must rely on active cooling measures. This study provides guiding suggestions for residents with cooling load demands by exploring improvement strategies for TH building performance. The aim is to help residents better cope with future IOH, reduce energy expenses, and alleviate the burden on urban electricity consumption.
The limitations of this study are that, although future climate changes are considered, the calculations only account for limited possibilities, treating each possibility's probability as equal. The actual uncertainty of climate change is much more complex than the settings used in this study. One approach to address this limitation is to use data from multiple Global Climate Models (GCMs) to more comprehensively consider climate change factors and obtain probability distributions regarding building performance in the future, providing more accurate assessment results. Additionally, the study calculates building performance using fixed building operation schedules without considering residents' complex indoor behaviors. The results cannot offer guidance to families with significantly different indoor behavioral patterns. Residents' behavior patterns are highly diverse, and developing a stochastic indoor behavior model requires extensive research. Addressing the mentioned limitations will be the focus of the research team's future work.
Author Contributions: Conceptualization, writing-review and editing, D.X.; formal analysis, writing-original draft preparation, W.X.; software, writing-review and editing, J.G.; resources, supervision, conceptualization, methodology, funding acquisition, Y.Z.; formal analysis, Z.W.; software, Y.F. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement:
Publicly available meteorological datasets were used in this study. This data can be found here: https://meteonorm.com/en/.