# Typical and extreme weather datasets for studying the resilience of buildings to climate change and heatwaves

Anaïs Machard, Agnese Salvati, Mamak P. Tootkaboni, Abhishek Gaur, Jiwei Zou, Liangzhu Leon Wang, Fuad Baba, Hua Ge, Facundo Bre, Emmanuel Bozonnet, et al.

► To cite this version:

Anaïs Machard, Agnese Salvati, Mamak P. Tootkaboni, Abhishek Gaur, Jiwei Zou, et al. Typical and extreme weather datasets for studying the resilience of buildings to climate change and heatwaves. Scientific Data, 2024, 11 (1), pp.531. 10.1038/s41597-024-03319-8. hal-04587669

HAL Id: hal-04587669

https://hal.science/hal-04587669v1

Submitted on 22 Jul 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
scientific **data**

<figure><img src="image_2.png" /></figure>

Check for updates

DATA DESCRIPTOR

OPEN

# Typical and extreme weather datasets for studying the resilience of buildings to climate change and heatwaves

Anaïs Machard<sup>1,2✉</sup>, Agnese Salvati<sup>3,4</sup>, Mamak P. Tootkaboni<sup>5</sup>, Abhishek Gaur<sup>6</sup>, Jiwei Zou<sup>7</sup>, Liangzhu Leon Wang<sup>1,7</sup>, Fuad Baba<sup>7,8</sup>, Hua Ge<sup>7</sup>, Facundo Bre<sup>1,9,10</sup>, Emmanuel Bozonnet<sup>1,2</sup>, Vincenzo Corrado<sup>1,5</sup>, Xuan Luo<sup>11</sup>, Ronnen Levinson<sup>11</sup>, Sang Hoon Lee<sup>11</sup>, Tianzhen Hong<sup>1,11</sup>, Marcello Salles Olinger<sup>1,2</sup>, Rayner Maurício e Silva Machado<sup>1,2</sup>, Emeli Lalesca Aparecida da Guarda<sup>1,2</sup>, Rodolfo Kirch Veiga<sup>1,2</sup>, Roberto Lamberts<sup>1,2</sup>, Afshin Afshari<sup>1,3</sup>, Delphine Ramon<sup>1,4</sup>, Hoang Ngoc Dung Ngo<sup>1,5,16</sup>, Abantika Sengupta<sup>1,5</sup>, Hilde Breesch<sup>1,5</sup>, Nicolas Heijmans<sup>1,7</sup>, Jade Deltour<sup>1,7</sup>, Xavier Kuborn<sup>1,7</sup>, Sana Sayadi<sup>1,8</sup>, Bin Qian<sup>1,19,20</sup>, Chen Zhang<sup>1,19</sup>, Ramin Rahif<sup>21</sup>, Shady Attia<sup>1,21</sup>, Philipp Stern<sup>1,22</sup> & Peter Holzer<sup>22</sup>

We present unprecedented datasets of current and future projected weather files for building simulations in 15 major cities distributed across 10 climate zones worldwide. The datasets include ambient air temperature, relative humidity, atmospheric pressure, direct and diffuse solar irradiance, and wind speed at hourly resolution, which are essential climate elements needed to undertake building simulations. The datasets contain typical and extreme weather years in the EnergyPlus weather file (EPW) format and multiyear projections in comma-separated value (CSV) format for three periods: historical (2001–2020), future mid-term (2041–2060), and future long-term (2081–2100). The datasets were generated from projections of one regional climate model, which were bias-corrected using multiyear observational data for each city. The methodology used makes the datasets among the first to incorporate complex changes in the future climate for the frequency, duration, and magnitude of extreme temperatures. These datasets, created within the IEA EBC Annex 80 “Resilient Cooling for Buildings”, are ready to be used for different types of building adaptation and resilience studies to climate change and heatwaves.

<footnote n="1">¹Scientific and Technical Building Research Centre (CSTB), Energy and Environment Department, Grenoble, France. ²Laboratory of Engineering Sciences for the Environment (LaSIE), La Rochelle University, La Rochelle, France. ³Architecture, Energy & Environment (AiE), Polytechnic University of Catalonia, Barcelona, Spain. ⁴Brunel University London, London, UK. ⁵Department of Energy, Politecnico di Torino, Torino, Italy. ⁶National Research Council Canada, Ottawa, Canada. ⁷Centre for Zero Energy Building Studies, Department of Building, Civil and Environmental Engineering, Concordia University, Montreal, Canada. ⁸Faculty of Engineering, British University in Dubai, Dubai, UAE. ⁹Institute of Construction and Building Materials, Technical University of Darmstadt, Darmstadt, 64287, Germany. ¹⁰Centro de Investigación de Métodos Computacionales (CIMEC), UNL/CONICET, Santa Fe, 3000, Argentina. ¹¹Lawrence Berkeley National Laboratory (LBNL), Berkeley, California, USA. ¹²Laboratory for Energy Efficiency in Buildings, Federal University of Santa Catarina, Florianópolis, Brazil. ¹³Fraunhofer Institute for Building Physics IBP, Valley, Germany. ¹⁴Department of Architecture, Faculty of Engineering Science, KU Leuven, Belgium. ¹⁵Building Physics and Sustainable Design, Department of Civil Engineering, KU Leuven, Belgium. ¹⁶Department of Environmental Architecture, Faculty of Architecture and Planning, Hanoi University of Civil Engineering, Ha Noi, Viet Nam. ¹⁷Buildwise, Brussels, Belgium. ¹⁸University of Gävle, 801 76, Gävle, Sweden. ¹⁹Department of the Built Environment, Aalborg University, Aalborg, Denmark. ²⁰China Southwest Architecture Design and Research Institute Corp. Ltd., Chengdu, China. ²¹Sustainable Building Design Lab, Department UEE, Faculty of Applied Sciences, Université de Liège, Liège, Belgium. ²²Institute of Building Research and Innovation (IBR&I), Wien, Austria.
✉e-mail: anais.machard@cstb.fr</footnote>
## Background & Summary

Climate change is among the most significant challenges the global community faces in the 21st century, with direct consequences for the building sector. An increase in the magnitude, frequency, and intensity of natural hazards presents a threat to the structural integrity of the buildings. In contrast, changes in climate characteristics, such as rising temperatures and more frequent extreme heat events, present an unprecedented challenge to building designers to design buildings that can perform efficiently over their durations of use. The performance evaluation of renovated or new buildings should consider not only the current average and extreme climates but also expected future climates and extreme events. To achieve this aim, reliable weather files capturing present, future typical, and extreme weather conditions are necessary to carry out building and resilience strategies studies. To reduce the computational costs associated with running building simulation models over long periods of time, simulations are generally performed over subsets of long-term climate data, typically over one year, referred to as reference meteorological years. Depending on the application, either a typical meteorological year (TMY) or an extreme meteorological year (XMY) is chosen. Many researchers and building practitioners are currently using only future TMYs to assess the impact of climate change on building energy performance because future TMYs are easily accessible and usually built from simplified statistical methods to account for climate change (e.g., the morphing method from Belcher et al.¹). Although morphing offers a quick way to generate weather files, it does not account for complex future changes in climate variables, such as changes in the frequency and duration of extreme heat events. Therefore, the generation of future weather files containing extremes has been an ongoing challenge for the building community in the last decade. A few authors have started to use climate model outputs directly to prepare the building simulation weather files to assemble not only future TMYs but also future extreme weather files such as heatwave events (HWE) or extreme meteorological years (XMYs). For example, Nik² prepared typical and extreme weather files for Stockholm and Geneva. The typical and extreme years were selected solely based on the temperature parameter. These weather files were prepared from raw regional climate model (RCM) data from four different climate models without bias correction. Machard et al.³ prepared typical TMY and future HWE for France using data from four RCM and the Representative Concentration Pathway (RCP) 8.5 at 12.5-km spatial resolution. In Machard⁴, bias-adjustment of the RCM projections was added to the method. The typical years were assembled following ISO EN 15927-4⁵, giving equivalent weight to temperature, humidity, and solar irradiance and secondary weight to wind speed. The heatwaves were selected following the French national heatwave definition, based on daily daytime and nighttime temperatures above specific thresholds validated for France using a CORDEX dataset by Ouzeau⁶. Doutreloup et al.⁷ and Ramon et al.⁸,⁹ used a convection-permitting climate model at 2.8 km resolution driven by the EC-Earth RCM and coupled with the land-surface scheme TERRA_URB. Based on the bias-adjusted data⁹,¹⁰, they prepared TMYs for different locations in Belgium for an RCP 8.5 climate change scenario. They also prepared XMYs, selecting extreme months based on two parameters: temperature and solar irradiance. Gaur et al.¹¹,¹² used the Canadian RCM bias-corrected climate projections to prepare TMYs, typical and extreme moisture reference years, typical downscaled years, and extreme warm and extreme cold years for over 500 locations. Recently, Bass et al.¹³ published future TMYs for 18 cities in the United States based on six climate models and different socioeconomic scenarios, Shared Socioeconomic Pathways (SSP) 5 and RCP 8.5. The TMYs were assembled using data from six climate models to reduce individual model bias.

**Study scope.** Future weather files based on bias-corrected RCM predictions are not easily available to the building scientific community; therefore, a large-scale international collaborative effort was made to curate and produce extreme weather data covering major global cities subject to extreme heat hazards by adopting a standardized procedure. This study prepares building simulation weather files ready to be used by building researchers and practitioners to carry out building energy simulations that are novel in the following respects:

a) they have been prepared to employ a consistent methodology over 15 cities distributed across the globe in different continents and climate types for 10 climate zones worldwide, as defined by the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) 169-2013¹⁴ (Fig. 1);

b) the future weather files are prepared directly from regional climate model simulation results and hence are able to account for complex future changes such as heatwaves in the climate variables projected for each city;

c) the use of a multivariate bias-correction method is employed to correct the bias associated with the regional climate model simulations;

d) the reference typical years and extreme heatwave event files are provided for building energy and overheating applications; and

e) bias-corrected multi-year projections are also made available for additional research and other applications.

These datasets were developed for "Annex 80: Resilient Cooling of Buildings", a research project of the International Energy Agency (IEA) - Energy in Buildings and Communities Programme (EBC)¹⁵, to evaluate the resilience of different passive and active cooling strategies.

They are used within the framework defined in Attia et al.¹⁶ and applied in Rahif et al.¹⁷. These weather files are shared to conduct climate change adaptation studies such as overheating risk assessments or a rise in demand for air conditioning under future typical and extreme weather conditions. The multi-year dataset is provided in comma-separated values (CSV) format so that it can easily be used for adaptation studies in other fields of investigation.

**Selected cities.** The weather datasets have been generated for 15 cities representative of the ten climate zones of ASHRAE classification¹⁴. Cities were selected to include at least one city per zone in climate zones 0 to 6 because climate change is expected to markedly increase cooling demand in these zones¹⁸. Preference was given
<figure><img src="image_3.png" /><figcaption>Fig. 1 15 locations selected and ASHRAE 169–2013 climate classification¹⁴.</figcaption></figure>

<table><caption>Table 1. Cities analyzed and population data⁹⁶.</caption><thead><tr><td>ASHRAE Climate Zone (CZ)</td><td>City</td><td>Population 2022 (M)</td><td>Change % (since 2021)</td><td>Country</td><td>Continent</td></tr></thead><tbody><tr><td>0A</td><td>Singapore</td><td>6.0</td><td>0.80%</td><td>Singapore</td><td>Asia</td></tr><tr><td>0B</td><td>Abu Dhabi</td><td>1.5</td><td>1.86%</td><td>UAE</td><td>Asia</td></tr><tr><td>1A</td><td>Guayaquil</td><td>3.1</td><td>1.62%</td><td>Ecuador</td><td>South America</td></tr><tr><td>2A</td><td>Sao Paulo</td><td>22.4</td><td>0.86%</td><td>Brazil</td><td>South America</td></tr><tr><td>3A</td><td>Buenos Aires</td><td>15.4</td><td>0.74%</td><td>Argentina</td><td>South America</td></tr><tr><td>3A</td><td>Rome</td><td>4.3</td><td>0.47%</td><td>Italy</td><td>Europe</td></tr><tr><td>3B</td><td>Los Angeles</td><td>4.0</td><td>0.05%</td><td>USA</td><td>North America</td></tr><tr><td>4A</td><td>Brussels</td><td>2.1</td><td>0.67%</td><td>Belgium</td><td>Europe</td></tr><tr><td>4A</td><td>Ghent</td><td>0.3</td><td>0.48%</td><td>Belgium</td><td>Europe</td></tr><tr><td>4A</td><td>London</td><td>9.5</td><td>1.22%</td><td>UK</td><td>Europe</td></tr><tr><td>4C</td><td>Vancouver</td><td>2.6</td><td>0.97%</td><td>Canada</td><td>North America</td></tr><tr><td>5A</td><td>Toronto</td><td>6.3</td><td>0.93%</td><td>Canada</td><td>North America</td></tr><tr><td>5A</td><td>Copenhagen</td><td>1.4</td><td>0.85%</td><td>Denmark</td><td>Europe</td></tr><tr><td>6A</td><td>Montreal</td><td>4.3</td><td>0.68%</td><td>Canada</td><td>North America</td></tr><tr><td>6A</td><td>Stockholm</td><td>1.7</td><td>1.36%</td><td>Sweden</td><td>Europe</td></tr></tbody></table>

to cities with high populations and high population growth. Most are in Europe, North America, and Asia due to the limitations of gathering observational data for other locations. However, these are also the continents where the most heatwave events have been recorded in the last decade¹⁹. The cities of interest and population data are presented in Table 1.

## Methods

The flow chart in Fig. 2 illustrates the steps adopted to generate the weather files. In step 1, raw climate data were extracted for the different weather variables that dominantly affect the thermal performance of buildings for historical and two future periods (20 years for each period). In step 2, these raw climate data were bias-corrected using observations of the different weather variables for the specific locations. In step 3, the weather files were assembled from the multiyear bias-adjusted datasets to generate (a) TMYs based on the EN ISO 15927-4 standard⁵ and (b) heatwave years (HWYs), based on the method to detect the heatwaves on a CORDEX dataset proposed by Ouzeau *et al.*⁶, already tested for building performance simulations in²⁰. Our methods are detailed in the following sections.

**Boundary conditions.** The historical and future projected climate simulations needed to prepare the weather files were taken from the Coordinated Regional Downscaling Experiment (CORDEX)²¹,²² results contributed by the scientific community towards the Coupled Model Intercomparison Project 5th Phase (CMIP5²³). The CORDEX climate datasets for CMIP6²⁴ were not available at the time our datasets were being prepared, so they were not considered. The future projections made under the Representative Concentration Pathway (RCP) 8.5 were considered²⁵. RCP 8.5 is the highest baseline scenario in which emissions rise throughout the twenty-first century. In this scenario, the emissions and concentrations of greenhouse gases rise significantly over time, causing a radiative forcing of 8.5 W/m² by the end of the century²⁶. This scenario is the most conservative greenhouse
<figure><img src="image_4.png" /><figcaption>Fig. 2 Methodology used for the weather datasets generation.</figcaption></figure>

gas emission scenario of the Coupled Model Intercomparison Project 5th Phase (CMIP5) which is also in line with the current emission trajectories of greenhouse gases around the globe²⁷ and therefore RCP 8.5 was chosen to evaluate the worst case possible in a resilience and adaptation context.

To select an appropriate climate simulation from the CORDEX database, i.e., containing data for many different General Circulation Model (GCM) and Regional Climate Model (RCM) combinations, we referred to the findings of McSweeney *et al.*²⁸. These authors analyzed all GCMs participating in the CORDEX database, and three reliable GCMs with low, medium, and high global equilibrium climate sensitivity (ECS) were identified as NCC-NORESM (Norwegian Earth System Model, developed by the Norwegian Climate Center), MPI-ESM-LR (Max Planck Institute Earth System Model for the High-Resolution Model), and HadGEM-ES (Hadley Centre Global Environment Model with an Earth-System configuration), respectively. These three GCMs have also been used to conduct coordinated downscaling experiments in CORDEX CORE simulations²⁹. In addition to this, we conducted a review of available CORDEX simulations at the needed temporal frequency (sub-daily) across different CORDEX domains encompassing the different cities we are analyzing. The dry-bulb temperature projections of these three climate models were compared with reference to the evaluation of the climate models report (contribution of Working Group I to the IPCC AR5). Finally, the MPI-ESM-LR (GCM) and REMO (RCM) combination was selected for this work as it was associated with medium global ECS, was found to be the closest to the median temperature of all climate model projections (Fig. 3) and contained simulations in the required temporal frequency (at least 3-hourly or more frequent) for all domains. This selected simulation is henceforth referred to as “MPI-REMO”.

**Downscaled climate simulations.** The selected GCM, MPI-ESM-LR³⁰, is dynamically downscaled by means of an RCM, REMO³¹,³². REMO is a three-dimensional atmosphere model developed at the Max Planck Institute for Meteorology in Hamburg, Germany, and currently maintained at the Climate Service Center Germany (GERICS) in Hamburg. The model is based on the Europa Model, the former NWP model of the German Weather Service. The prognostic variables in REMO are horizontal wind components, surface pressure, air temperature, specific humidity, cloud liquid water, and ice. The physical packages originate from the global circulation model ECHAM4³³, although many updates have been introduced³⁴⁻⁴¹.

The MPI-REMO simulations, summarized in Table 2, were of 12.5 km spatial resolution for the European domain and 25 km resolution for other domains.

RCM files were stored for each weather variable and for one year on the entire domain grid (a domain usually corresponds to an entire continent or parts of a continent) in NETCDF4 format. A Python code provided with this dataset was used to download the different NETCDF4 files, extract the nearest point to each city coordinates, and assemble the different weather variables and years in a single dataset. For each city, the weather variables downloaded are described in Table 3. They include dry-bulb temperature, specific or relative humidity, atmospheric pressure, surface downwelling shortwave irradiance, wind speed, and cloud cover (only for Europe). Additional variables, such as rainfall, wind direction, or longwave irradiance, are also important, but they were not available for all the cities; therefore, they were not used. Data were downloaded for the three time periods referenced in Table 4.

**Bias correction of climate model simulations.** Climate model simulations are known to have bias associated with them because of the coarse spatial resolution at which the global or regional climate simulations
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Selection of the climate model to generate future weather datasets – Position of the temperature projection from HadGEM2-ES, MPI-ESM-LR, and NorESM1-M in comparison with other model climate projections. Modified from: Flato, Gregory, *et al.* 'Evaluation of climate models.' Climate change 2013<sup>69</sup>.</figcaption></figure>

<table><caption>**Table 2.** Climate projections (model, scenario, spatial, and time frequency) used for each location.</caption><thead><tr><th>Continent</th><th>Domain</th><th>Driving model</th><th>Downscaling method</th><th>Socio-economic scenario</th><th>Time frequency</th></tr></thead><tbody><tr><td>Africa</td><td>AFR-22</td><td></td><td></td><td></td><td>3 HOURS</td></tr><tr><td>Asia</td><td>SEA-22</td><td></td><td></td><td></td><td>3 HOURS</td></tr><tr><td>Europe</td><td>EU-11</td><td>MPI-ESM-LR</td><td>REMO 2015</td><td>RCP 8.5</td><td>1 HOUR</td></tr><tr><td>South America</td><td>SAM-22</td><td></td><td></td><td></td><td>3 HOURS</td></tr><tr><td>North America</td><td>NAM-22</td><td></td><td></td><td></td><td>3 HOURS</td></tr></tbody></table>

<table><caption>**Table 3.** Weather variables downloaded from the CORDEX platform. *hurs is required in weather files for building performance simulations but was available only for the EU and SAM domains. For the other domains, the huss and tas variables are used to recalculate the hurs.</caption><thead><tr><th>EUR 11 Domain</th><th>AFR 22, NAM 22, SAM 22, and SEA 22 Domains</th></tr></thead><tbody><tr><td>tas (near-surface air temperature)</td><td>tas (near-surface air temperature)</td></tr><tr><td>hurs (near-surface relative humidity)*</td><td>n/a</td></tr><tr><td></td><td>huss (near-surface specific humidity)*</td></tr><tr><td>ps (surface air pressure)</td><td>ps (surface air pressure)</td></tr><tr><td>rsds (surface downwelling shortwave irradiance)</td><td>rsds (surface downwelling shortwave irradiance)</td></tr><tr><td>clt (total cloud fraction)</td><td>n/a</td></tr><tr><td>sfcWind (near-surface wind speed)</td><td>sfcWind (near-surface wind speed)</td></tr><tr><td>clt (cloud cover)</td><td></td></tr></tbody></table>

are conducted<sup>42</sup>. The biases in the climate simulations, if left uncorrected, have been known to lead to incorrect descriptions of climate-driven hazards, such as floods<sup>43</sup> and wildfires<sup>44</sup>. Many bias-correction methods have been discussed in literature<sup>42</sup>. The complexity of the methods can range from methods correcting simply the mean bias<sup>45</sup> to methods able to perform univariate and multivariate distribution-based corrections<sup>46</sup>. The multivariate bias-correction methods have been found most efficient in correcting bias in the marginal distribution of the climate variables, as well as the inter-relationships between the variables, and have been recommended for accurately describing hazards dependent on multiple climate variables<sup>46</sup>. Therefore, the bias correction of raw climate variables was performed using quantile delta mapping (QDM)<sup>47</sup> and Multivariate Bias Correction with N-dimensional probability density function transform (MBCn)<sup>46</sup> methods. The QDM is a univariate bias-correction method that preserves climate model projected future changes in the quantiles of climate variables while at the same time correcting systematic biases in the quantiles. Climate model data are de-trended and then mapped onto
<table><caption>Table 4. 20-year periods downloaded for each variable from the CORDEX platform.</caption><thead><tr><td>Period Name</td><td>Years</td></tr></thead><tbody><tr><td>Historical – 2010s</td><td>2001–2020</td></tr><tr><td>Mid-term future – 2050s</td><td>2041–2060</td></tr><tr><td>Long-term future – 2090s</td><td>2081–2100</td></tr></tbody></table>

<table><caption>Table 5. Observational data used in the bias-correction step for each city. x* = hourly values estimated from daily values.</caption><thead><tr><th rowspan="2">CZ</th><th rowspan="2">City</th><th rowspan="2">Latitude (°)</th><th rowspan="2">Longitude (°)</th><th rowspan="2">Observational data Period</th><th colspan="6">Data used for bias correction</th></tr><tr><th>tas</th><th>hurs</th><th>rsds</th><th>sfcWind</th><th>ps</th><th>clt</th></tr></thead><tbody><tr><td>0A</td><td>Singapore</td><td>1.37</td><td>103.98</td><td>1996–2015</td><td>x</td><td>x</td><td></td><td>x</td><td>x</td><td></td></tr><tr><td>0B</td><td>Abu Dhabi</td><td>24.42</td><td>54.61</td><td>2008–2012</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td></tr><tr><td>1A</td><td>Guayaquil</td><td>-2.15</td><td>-79.92</td><td>07-2016 – 08-2020</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td><td></td></tr><tr><td>2A</td><td>Sao Paulo</td><td>-23.63</td><td>-46.65</td><td>1986–2005</td><td>x</td><td>x</td><td>x*</td><td>x*</td><td>x</td><td></td></tr><tr><td>3A</td><td>Buenos Aires</td><td>-34.56</td><td>-58.42</td><td>1986–2005</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td></tr><tr><td>3A</td><td>Rome</td><td>41.81</td><td>12.25</td><td>2008–2017</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td></tr><tr><td>3B</td><td>Los Angeles</td><td>33.93</td><td>-118.40</td><td>2000–2019</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td></tr><tr><td>4A</td><td>Brussels</td><td>50.80</td><td>4.36</td><td>2009–2018</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td><td></td></tr><tr><td>4A</td><td>Ghent</td><td>51.05</td><td>3.73</td><td>2009–2020</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td></tr><tr><td>4A</td><td>London</td><td>51.48</td><td>-0.45</td><td>1996–2015</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td></tr><tr><td>4C</td><td>Vancouver</td><td>49.19</td><td>-123.18</td><td>1998–2017</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td></tr><tr><td>5A</td><td>Toronto</td><td>43.67</td><td>-79.40</td><td>1998–2017</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td></tr><tr><td>5A</td><td>Copenhagen</td><td>55.88</td><td>12.41</td><td>2001–2019</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td></tr><tr><td>6A</td><td>Montreal</td><td>45.63</td><td>-73.55</td><td>1998–2017</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td></td></tr><tr><td>6A</td><td>Stockholm</td><td>59.9</td><td>18.03</td><td>1986–2005</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td><td>x</td></tr></tbody></table>

the observations using quantile mapping. After that, future projected bias-corrected datasets are obtained by multiplying/adding to them the climate model projected future relative/additive changes in quantiles. The MBCn method extends the application of the QDM method in a multivariate context. First, individual climate variables are corrected following the QDM method. Thereafter, the dependence structure of climate variables is corrected using an iterative reshuffling process where, in each iteration, climate data are rotated by multiplying them with random orthogonal matrices, QDM is corrected and then re-correlated using inverse random matrices.

While all climate variables were bias corrected using the MBCn method, the QDM method was used to correct global solar irradiance because our analysis shows that the reshuffling of marginally corrected global solar irradiance values, as performed in the MBCn method, breaks the diurnal structure of global solar irradiance. This can subsequently lead to unrealistic values for not only global solar irradiance but also direct and diffused solar irradiance components derived from it. The calibration of MBCn/QDM methods and subsequent prediction of bias-corrected values were performed individually for each month of the year to preserve month-to-month variability in bias-corrected climate data. The methods assume that the bias is the same in the future as in the present. All years with observational data available in different cities were considered for the calibration of bias-correction methods. The length of the observational period and the variables available for each city are reported in Table 5. The observational datasets included hourly values of air temperature (tas), relative humidity (hurs), global horizontal irradiation (rsds), wind speed (sfcWind), atmospheric pressure (ps), and cloud cover (clt). Just for Sao Paulo, for which hourly values could not be found for all weather variables, the hourly values of global horizontal irradiation and wind speed were derived from daily values. Hourly values for irradiation were calculated using the Zhang-Huang solar model⁴⁸. The regression coefficients in the model were calibrated based on the daily values using a least-squares approximation. Hourly values for wind speed were obtained by adjusting the monthly cumulative frequency distributions of historical RCM data to the observational data. Hence, each day of the RCM had its wind speed hourly values multiplied by a factor to match the cumulative frequency of the observational data daily mean values. For Abu Dhabi, data on atmospheric pressure could not be found; a static standard atmospheric pressure was used since the city is located at sea level. Note that observations of solar radiation were not available for Singapore, so its solar irradiance was not bias-corrected when the datasets were prepared.

The coordinates given in Table 5 correspond to the location of the weather station where the observations were used for bias correction for each city. The chosen weather stations are located outside of the cities, usually at airport sites; therefore, the observations and the resulting bias-corrected datasets do not account for urban heat island effects (UHI). We decided not to include urban effects in these datasets for various reasons. First, urban observations are not available for some of the cities analyzed. Secondly, the UHI is not homogeneous across a city, varying significantly depending on the different local climate zones (LCZ). Therefore, it would be necessary to create more than one urban weather file for each city, namely one for each LCZ. Furthermore, it would not be correct to use current urban observations as a reference for future UHI intensities because building
density, vegetation, materials, and anthropogenic heat generation in future cities will probably change, leading to a change in UHI intensity. For all these reasons, even if the datasets refer to cities, they do not include urban effects, like most of the currently available weather datasets for building performance simulations. They can be modeled and added to the datasets in post-processing by using tools and methodologies that are discussed and referenced in the “Usage notes” section.

**Calculating direct and diffuse solar irradiance.** The Boland-Ridley model⁴⁹ was used to calculate the direct and diffuse components of global solar irradiance. This method is a robust and straightforward predictor model that requires few inputs. The Italian National organization for standardization (UNI) has adopted this reliable method to split the global solar irradiance for creating national climatic data (UNI 10349-1:2016)⁵⁰. The model was also validated in a later study⁵¹. The Boland-Ridley model uses a logistic function (sigmoid function) for the diffuse fraction of global solar irradiance on a horizontal surface based on the sky clearness index, which is the ratio of the terrestrial global horizontal solar irradiance to the extraterrestrial horizontal solar irradiance. The extraterrestrial horizontal solar irradiance is calculated from the solar elevation and the extra-atmospheric solar irradiance received on a theoretical surface orthogonal to the sun's rays and at the Earth's mean distance from the sun (depending on the Earth's orbital angle). This fraction includes both the horizontal direct and diffuse solar irradiance components of horizontal solar irradiance. This model is used for the generation of direct-normal solar irradiance⁵², which is required for building energy simulation. It is computed as the ratio of the direct horizontal solar irradiance to the cosine of the solar zenith angle. Calculation of direct-normal solar irradiance can yield unphysical results when the direct-horizontal solar irradiance and the cosine of the solar zenith angle are both small because the sun is low. In this case, a threshold is introduced by applying a physical model⁵³ that considers the Rayleigh optical depth (in the function of the air mass) and the Linke Turbidity (TL)⁵⁴, which accounts for scattering and absorption by both atmospheric aerosols and atmospheric gases.

**Creating typical years from multiyear hourly datasets.** The TMYs were created using the international standard EN ISO 15927-4 – Hygrothermal performance of buildings, Calculation and presentation of climatic data, Part 4: Hourly data for assessing the annual energy use for heating and cooling method⁵⁵. The procedure is applicable for assessing the climate change impact on the long-term mean energy loads of buildings. However, this method based on average values is not suitable for studying extreme meteorological events. TMYs are constructed from 12 representative months (typical months) from multiyear records. Two sets of parameters are considered for selecting the typical months: primary parameters, including dry-bulb air temperature, global solar irradiance, and relative humidity (or air absolute humidity, water vapor pressure, or dew point temperature), and secondary parameters, including wind speed. For each primary climatic parameter, *p*, the daily means, $\bar{p}$, are calculated from all multi-year records of hourly values of *p* (at least ten years). After sorting the $\bar{p}$ values for a specific month, *m*, of all the years in increasing order, the cumulative distribution function is calculated for each parameter and *i*th day as:

$$ \Phi(p, m, i) = \frac{K(i)}{N + 1} \qquad (1) $$

where *K(i)* is the rank order of the *i*th day and *N* is the number of days for a month overall multi-year records. Afterward, the cumulative function is calculated for each year of the multi-year records for a specific month, *m*, and specific year, *y*, according to Eq. 2:

$$ F(p, y, m, i) = \frac{J(i)}{n + 1} \qquad (2) $$

where *J(i)* is the rank order of the *i*th day and *n* is the number of days for the specific month and year. Subsequently, the Finkelstein-Schafer statistic (Fs)⁵⁵ is calculated for all the primary climatic parameters for each calendar month and year of multi-year records. Fs is a goodness-to-fit statistic that proved more potent than conventional alternatives and is calculated as:

$$ F_s(p, y, m = \sum_{i=1}^{n} |F(p, y, m, i) - \Phi(p, m, i)|) \qquad (3) $$

For each calendar month and each year, Fs values are calculated and ranked in increasing order. By calculating the total ranking (the sum of the primary parameter's ranks) for each year, three months with the lowest total ranking are selected for each calendar month. The month with the lowest deviation in wind speed (secondary parameter) is selected as the typical month to be included in the typical year. This method was applied to the 20-year bias-corrected RCM data to generate one TMY for each period. The TMYs were then converted to EnergyPlus weather files (.EPW) for use in building energy simulations. The EnergyPlus auxiliary program “weather converter” tool⁵⁶ was used for this purpose.

**Selecting extreme heatwaves from multi-year datasets.** The method proposed by Ouzeau et al.⁶ was used to select heatwaves from the 20-year periods based on high quantiles of daily temperature distributions. The method was validated for France by comparing heatwave detection on an EURO-CORDEX regional multi-model ensemble with the French SAFRAN thermal indicator, historically used by French authorities for cold spell detection. The adopted method has the advantage of applying to different cities worldwide since it is based on relative thresholds and not absolute thresholds. It detects heatwaves based on three temperature thresholds calculated
<table><caption>**Table 6.** Datasets available for each city and data periods (Historical 2001–2020, Mid-term future 2041–2060, Long-term Future (2081–2100)).</caption><thead><tr><td>Category</td><td>Short description</td><td>Extension</td><td>Link</td></tr></thead><tbody><tr><td>MY - Multiyear dataset</td><td>A file containing hourly values of 20-years bias-corrected climate data</td><td>csv</td><td></td></tr><tr><td>TMY - Typical Meteorological Year</td><td>Weather file to run building performance simulations representative of typical meteorological conditions over 20 years</td><td>epw</td><td>https://www.wdc-climate.de/ui/entry?acronym=WDTF_Annex80_build_v1.0<sup>58</sup></td></tr><tr><td>HWY - Heatwave year (year containing heatwaves)</td><td>Weather file to run building performance simulations including extreme heatwaves (i.e., most severe, longest, or most intense over 20 years)</td><td>epw</td><td></td></tr></tbody></table>

from the historical multiyear period: The 99.5 threshold (99.5 percentile) is used to detect a temperature peak and a potential heatwave. The 97.5 threshold (97.5 percentile) is used to calculate the heatwave duration (days during which the temperature is above the threshold) and severity (degree-days above the threshold). If the temperature goes under this threshold for more than three consecutive days, the heatwave stops. The 95 threshold (95 percentile) is used to end the heatwave drastically if the temperature drops below this threshold. The chosen method was recently demonstrated to be the most effective in detecting and characterizing heat waves for building resilience analysis⁵⁷. The current work builds on the methodology initiated by Machard *et al.*³ to assemble future weather files, including heatwave for building energy and thermal performance simulations from CORDEX climate data. In the proposed approach, each heatwave is characterized by three criteria: intensity (maximum daily mean temperature °C reached during the heatwave), duration (in days), and severity (aggregated temperature above the 97.5 threshold in °C.day). Applying this method, many heatwaves were found during each multiyear period in each city. Since the purpose of the datasets is to carry out building performance resilience assessments, the three most extreme heatwaves were selected, according to these three criteria: the most intense, the most severe, and the longest heatwaves.

## Data Records

The entire datasets (Table 6) produced for this work are organized into three categories:

* Multiyear (MY)

* Typical meteorological year (TMY)

* Heatwave year (HWY)

The datasets are available at the link: https://www.wdc-climate.de/ui/entry?acronym=WDTF_Annex80_build_v1.0⁵⁸.

The first category of files is MY datasets in CSV format. There are three MY files for each city, containing the hourly values of the bias-corrected RCM variables for each 20-year reference period. The variables included in the CSV files are air temperature (tas), near-surface relative humidity (hurs), near-surface specific humidity (huss), surface atmospheric pressure (ps), surface downwelling shortwave irradiance (rsds), and wind speed (sfcWind). Some cities have fewer variables due to missing observational data to perform the bias-correction. Cloud cover (clt) is available for London and Stockholm. The MY file name format is: "climatezone.city.MY_referenceperiod.csv". For instance: "0B_Abu_Dhabi_MY_2081-2100". The MY files were used to create both TMYs and HWYs.

There are three TMYs per city, representing the typical meteorological conditions corresponding to historical (2001–2020), mid-term future (2041–2060), and long-term future (2081–2100) periods. The TMYs are provided in the EnergyPlus weather file (EPW) format. The EPW file details hourly dry bulb air temperature (°C), dew point temperature (°C), relative humidity (%), atmospheric pressure (Pa), global horizontal solar irradiance (Wh/m²), direct normal irradiance (Wh/m²), diffuse horizontal irradiance (Wh/m²), wind speed (m/s), and wind direction (°). For the cities of London and Stockholm, the total sky cover (tenths) is also provided. In TMYs, values for wind direction were extracted from the historical time series of METEONORM⁵⁹ for each city because wind direction is needed to perform building energy simulations but is not available for all CORDEX domains. The EPW files were generated using the EnergyPlus weather converter, auxiliary software of EnergyPlus⁵⁶.

The file name of each TMY has the following format: "climatezone.city.TMY_referenceperiod.epw". For instance, the file "4 A_London_TMY_2041-2060" is the TMY for the city of London, located in the ASHRAE climate zone 4A, for the mid-term future period (2041-2060).

Finally, the HWYs are also provided in EPW format. Each city can have a maximum of nine HWY files, corresponding to the years with the most intense, most severe, and longest heatwaves found in the three reference periods. As the most intense and/or the longest heatwaves are also the most severe in many cases, the total number of HWY files is generally less than nine. The HWY file name format is "climatezone.city.HW_referenceperiod_heatwavetype_year.epw". For instance, the file "6 A_Stockholm_HW_Historical_MostSevere_Longest_2002.epw" contains the most severe and longest heatwave occurring in the historical period, in 2002, in Stockholm (climate zone 6A).

## Technical Validation

For technical validation, the multiyear raw climate outputs, observations, and bias-adjusted datasets were compared and analyzed. The mean values of ambient air temperature, relative humidity, global solar irradiance, and wind speed in the typical years during the historical period were compared to the mean values in the multiyear datasets, showing good agreement in values. The extreme values of ambient air temperature for the heatwave
<table><caption>**Table 7.** Mean temperature, solar irradiance, wind speed, and relative humidity in the cities over the validation time period.</caption><thead><tr><th rowspan="2">CZ</th><th rowspan="2">City</th><th colspan="3">Temperature (°C)</th><th colspan="3">Solar irradiance (W/m²)</th><th colspan="3">Wind speed (m/s)</th><th colspan="3">Relative humidity (%)</th></tr><tr><th>OBS</th><th>RCM (raw)</th><th>RCM (bc)</th><th>OBS</th><th>RCM (raw)</th><th>RCM (bc)</th><th>OBS</th><th>RCM (raw)</th><th>RCM (bc)</th><th>OBS</th><th>RCM (raw)</th><th>RCM (bc)</th></tr></thead><tbody><tr><td>0A</td><td>Singapore</td><td>27.8</td><td>28.1</td><td>27.7</td><td>—</td><td>—</td><td>—</td><td>1.9</td><td>4.7</td><td>2.0</td><td>83.8</td><td>74.1</td><td>83.3</td></tr><tr><td>0B</td><td>Abu Dhabi</td><td>27.6</td><td>29.0</td><td>27.6</td><td>237.9</td><td>246.7</td><td>238.1</td><td>3.2</td><td>4.0</td><td>3.2</td><td>60.0</td><td>55.5</td><td>60.2</td></tr><tr><td>1A</td><td>Guayaquil</td><td>27.0</td><td>27.1</td><td>27.0</td><td>263.1</td><td>218.4</td><td>266.9</td><td>1.7</td><td>1.9</td><td>1.8</td><td>74.6</td><td>77.4</td><td>74.3</td></tr><tr><td>2A</td><td>Sao Paulo</td><td>19.3</td><td>19.4</td><td>19.3</td><td>188.6</td><td>187.9</td><td>188.6</td><td>6.1</td><td>2.73</td><td>6.1</td><td>80.6</td><td>77.1</td><td>80.6</td></tr><tr><td>3A</td><td>Buenos Aires</td><td>18.0</td><td>19.1</td><td>17.7</td><td>191.4</td><td>197.5</td><td>190.0</td><td>4.5</td><td>4.5</td><td>4.4</td><td>72.1</td><td>66.9</td><td>72.3</td></tr><tr><td>3A</td><td>Rome</td><td>16.3</td><td>16.5</td><td>16.3</td><td>187.8</td><td>163.9</td><td>188.0</td><td>3.6</td><td>2.7</td><td>3.6</td><td>72.5</td><td>70.0</td><td>72.4</td></tr><tr><td>3B</td><td>Los Angeles</td><td>16.7</td><td>20.7</td><td>16.7</td><td>215.1</td><td>223.3</td><td>214.7</td><td>1.7</td><td>2.4</td><td>1.7</td><td>72.3</td><td>57.5</td><td>72.4</td></tr><tr><td>4A</td><td>Brussels</td><td>10.8</td><td>11.1</td><td>10.8</td><td>127.2</td><td>109.1</td><td>127.2</td><td>3.6</td><td>3.7</td><td>3.6</td><td>78.5</td><td>82.0</td><td>78.5</td></tr><tr><td>4A</td><td>Ghent</td><td>11.1</td><td>11.3</td><td>11.1</td><td>126.4</td><td>110.3</td><td>126.3</td><td>3.4</td><td>4.1</td><td>3.4</td><td>78.6</td><td>82.8</td><td>78.7</td></tr><tr><td>4A</td><td>London</td><td>11.6</td><td>11.1</td><td>11.8</td><td>118.8</td><td>106.8</td><td>117.5</td><td>4.2</td><td>2.8</td><td>4.0</td><td>75.5</td><td>79.5</td><td>75.7</td></tr><tr><td>4C</td><td>Vancouver</td><td>10.1</td><td>7.8</td><td>10.6</td><td>142.7</td><td>130.5</td><td>153.7</td><td>3.7</td><td>3.3</td><td>4.4</td><td>78.9</td><td>69.2</td><td>74.5</td></tr><tr><td>5A</td><td>Toronto</td><td>9.2</td><td>6.6</td><td>7.8</td><td>159.0</td><td>146.6</td><td>153.7</td><td>4.4</td><td>3.6</td><td>4.4</td><td>69.3</td><td>79.1</td><td>69.2</td></tr><tr><td>5A</td><td>Copenhagen</td><td>8.8</td><td>9.3</td><td>8.8</td><td>118.2</td><td>102.1</td><td>118.2</td><td>3.3</td><td>4.6</td><td>3.3</td><td>82.4</td><td>84.9</td><td>82.4</td></tr><tr><td>6A</td><td>Montreal</td><td>7.7</td><td>4.9</td><td>7.8</td><td>153.8</td><td>134.8</td><td>153.7</td><td>4.4</td><td>3.5</td><td>4.4</td><td>69.2</td><td>83.6</td><td>69.2</td></tr><tr><td>6A</td><td>Stockholm</td><td>6.6</td><td>6.4</td><td>6.6</td><td>116.5</td><td>92.9</td><td>116.6</td><td>3.9</td><td>3.0</td><td>3.9</td><td>79.6</td><td>86.0</td><td>79.6</td></tr></tbody></table>

years were compared to the extreme of the multiyear datasets. An assessment of the future weather files confirms that climate change will increase the mean temperature in all cities. Heatwave frequency, intensity, and duration will also increase in all cities and more drastically in the four hottest cities (Singapore, Abu Dhabi, Guayaquil, and Sao Paulo) analyzed.

**Comparison of raw-output and bias-corrected data.** The validation of the bias-correction step was performed by comparing bias-corrected climate estimates with observations over a validation time-period that varies from city to city depending on the time period of observations available to them. The validation time period is considered the period overlapping between observational and historical time-periods. This allowed us to make use of the entire length of observational data available in different cities for performing validation of bias-correction methods. The validation results show that the QDM/MBCn methods were able to reduce the bias associated with RCM simulations effectively. This can be seen from the results presented in Table 7, in which mean climate statistics from observations, raw RCM, and bias-corrected (bc) RCM are presented for the validation time period. The results show that the projected temperature, solar irradiance, wind speed, and relative humidity from raw RCMs have noticeable bias, which is reduced by the application of the bias-correction step. For instance, RCM over-predicts the mean temperature in Singapore by 0.5 °C, which is effectively eliminated after the bias correction. Table 8 presents the standard deviation of observations (OBS), RCM-raw, and RCM-bs for these four climate variables, which also shows the bias reduction between OBS and bias-corrected RCM data. Not only is the bias correction effective in correcting bias in average climate characteristics over the cities, but it also reduces bias across the whole distribution of climate variables. This is evident from Fig. 4, in which probability density functions (PDFs) of temperature, wind speed, and relative humidity from observations (grey), raw RCM (blue), and bias-corrected RCM (red) datasets are presented for Singapore, London, and Toronto. PDFs of raw RCM are effectively adjusted by the bias-correction procedure to mimic the PDFs of observations. This is true not only for temperature but also for relatively more complex variables such as wind speed, highlighting the effectiveness of the bias-correction step in simulating realistic estimates of a range of climate variables considered in this study.

**Projected changes in weather variables over multi-year (MY) future periods.** The values of mean temperature, solar irradiance, wind speed, and relative humidity over the 2010s, 2050s, and 2090s for all cities are presented in Table 9. In general, between the historical period (the 2010s) and the two future time periods (2050s, 2090s), mean temperatures are projected to increase in all cities located in different climate zones (CZs). In most cities, the increase in MY by the 2050s is about 1 °C, while it will be about 2-3 °C by the 2090s, with the largest increase of 4.2 °C in Abu Dhabi (CZ: 0B – Extremely Hot Dry) and the smallest increase of 1.6 °C in Buenos Aires (CZ: Warm Humid). Mean temperature increases within the same ASHRAE climate zone are consistent: in zone ASHRAE CZ: 4 A - Mixed Humid, the temperature increase in Brussels, Ghent, and London are about 0.8 °C, 0.7 °C, and 0.7 °C between MY-2050s and MY-2010s, and of 2.6 °C, 2.6 °C, and 2.5 °C between MY-2090s and MY-2010s. Global solar irradiance is projected to decrease in the majority of the cities, with the largest decrease of 12.8 W/m² by the 2090s is projected for Stockholm (CZ: Cold Humid), whereas a slight increase of 0.6 W/m² is projected for Abu Dhabi (CZ 0B -: Extremely Hot Dry). Such a reduction in future solar irradiance was also found in other studies⁴⁵,⁶⁰. According to Cutforth and Judiesch⁶¹, this can be the consequence of two factors: (1) higher attenuation of solar irradiance from increased aerosol concentrations and sometimes from increasing cloudiness, and (2) an increase in annual number of precipitation events. These assumptions are coherent since the irradiance is not decreasing in Abu Dhabi, for which cloud cover is very low. However, this trend in decreasing
<table><caption>**Table 8.** Standard deviation of temperature, solar irradiance, wind speed, and relative humidity in the cities over the validation time period.</caption><thead><tr><th rowspan="2">CZ</th><th rowspan="2">City</th><th colspan="3">Temperature (°C)</th><th colspan="3">Solar irradiance (W/m²)</th><th colspan="3">Wind speed (m/s)</th><th colspan="3">Relative humidity (%)</th></tr><tr><th>OBS</th><th>RCM (raw)</th><th>RCM (bc)</th><th>OBS</th><th>RCM (raw)</th><th>RCM (bc)</th><th>OBS</th><th>RCM (raw)</th><th>RCM (bc)</th><th>OBS</th><th>RCM (raw)</th><th>RCM (bc)</th></tr></thead><tbody><tr><td>0A</td><td>Singapore</td><td>2.2</td><td>1.5</td><td>2.1</td><td>—</td><td>—</td><td>—</td><td>1.6</td><td>1.9</td><td>1.6</td><td>9.9</td><td>6.9</td><td>10.0</td></tr><tr><td>0B</td><td>Abu Dhabi</td><td>7.9</td><td>7.7</td><td>7.9</td><td>312.4</td><td>326.9</td><td>313.2</td><td>2.2</td><td>2.1</td><td>2.2</td><td>20.4</td><td>21.4</td><td>20.5</td></tr><tr><td>1A</td><td>Guayaquil</td><td>3.4</td><td>3.3</td><td>3.4</td><td>380.8</td><td>312.1</td><td>382.3</td><td>1.0</td><td>0.9</td><td>1.0</td><td>12.4</td><td>14.4</td><td>12.0</td></tr><tr><td>2A</td><td>Sao Paulo</td><td>4.7</td><td>4.6</td><td>4.7</td><td>256.4</td><td>288.3</td><td>256.4</td><td>2.9</td><td>1.3</td><td>3.0</td><td>16.4</td><td>16.6</td><td>16.4</td></tr><tr><td>3A</td><td>Buenos Aires</td><td>5.6</td><td>5.1</td><td>5.7</td><td>283.5</td><td>297.3</td><td>282.5</td><td>2.4</td><td>1.9</td><td>2.4</td><td>15.2</td><td>16.3</td><td>15.4</td></tr><tr><td>3A</td><td>Rome</td><td>7.1</td><td>7.3</td><td>7.1</td><td>270.4</td><td>255.6</td><td>270.4</td><td>2.2</td><td>1.7</td><td>2.2</td><td>16.7</td><td>17.6</td><td>16.7</td></tr><tr><td>3B</td><td>Los Angeles</td><td>4.4</td><td>7.4</td><td>4.4</td><td>295.9</td><td>310.6</td><td>296.1</td><td>1.0</td><td>1.3</td><td>1.0</td><td>22.3</td><td>22.8</td><td>22.2</td></tr><tr><td>4A</td><td>Brussels</td><td>6.8</td><td>6.8</td><td>6.8</td><td>196.2</td><td>201.0</td><td>196.1</td><td>1.8</td><td>1.8</td><td>1.8</td><td>14.3</td><td>14.1</td><td>14.3</td></tr><tr><td>4A</td><td>Ghent</td><td>6.8</td><td>6.6</td><td>6.8</td><td>201.0</td><td>193.5</td><td>200.6</td><td>1.9</td><td>2.0</td><td>1.9</td><td>15.4</td><td>13.5</td><td>15.2</td></tr><tr><td>4A</td><td>London</td><td>6.1</td><td>6.3</td><td>6.1</td><td>193.1</td><td>195.9</td><td>191.3</td><td>2.2</td><td>1.2</td><td>2.2</td><td>15.9</td><td>14.0</td><td>15.7</td></tr><tr><td>4C</td><td>Vancouver</td><td>5.3</td><td>12.0</td><td>6.0</td><td>231.2</td><td>511.1</td><td>235.4</td><td>2.3</td><td>2.2</td><td>2.5</td><td>13.0</td><td>16.8</td><td>12.7</td></tr><tr><td>5A</td><td>Toronto</td><td>10.9</td><td>9.2</td><td>12.0</td><td>241.4</td><td>222.8</td><td>235.4</td><td>2.7</td><td>1.8</td><td>2.5</td><td>16.2</td><td>14.5</td><td>16.8</td></tr><tr><td>5A</td><td>Copenhagen</td><td>7.2</td><td>6.3</td><td>7.2</td><td>196.0</td><td>191.4</td><td>196.1</td><td>2.1</td><td>2.2</td><td>2.1</td><td>15.4</td><td>11.7</td><td>15.4</td></tr><tr><td>6A</td><td>Montreal</td><td>12.0</td><td>8.7</td><td>12.0</td><td>235.4</td><td>213.4</td><td>253.4</td><td>2.6</td><td>1.8</td><td>2.5</td><td>16.8</td><td>13.4</td><td>16.8</td></tr><tr><td>6A</td><td>Stockholm</td><td>7.9</td><td>7.7</td><td>7.9</td><td>184.7</td><td>178.0</td><td>184.7</td><td>1.7</td><td>1.4</td><td>1.7</td><td>14.6</td><td>11.7</td><td>14.6</td></tr></tbody></table>

global solar irradiance cannot be generalized. It can be due to a coarse representation of rain and cloud events at the model spatial resolution (25 or 50 km depending on the CORDEX domain) and to potential biases for this climate parameter in the selected climate model. In terms of wind speed and relative humidity, a general change is not observed. Most cities have minimal change in wind projections in the future: The largest decrease of 0.4 m/s in wind speed is projected for Buenos Aires (CZ: Warm Humid), whereas the largest increase of 0.3 m/s in wind speed is projected for Sao Paulo (CZ: Hot Humid). Finally, the largest variability in the sign of projected future change is obtained for relative humidity. While the cities of Singapore (CZ: 0 A - Extremely Hot Humid), Guayaquil (CZ: 1 A - Very Hot Humid), Buenos Aires (CZ: 3 A - Warm Humid), Los Angeles (CZ: 3B - Warm Dry) are projected to experience increases in relative humidity of up to 5%, the cities of Sao Paulo (CZ: 2 A - Hot Humid) and Abu Dhabi (CZ: 0B - Extremely Hot Dry) are projected to experience future decreases of up to 4%. Smaller future changes in relative humidity are projected for other cities such as Montreal and Stockholm (CZ: 6 A - Cold Humid) as well as Ghent, Brussels, and London (CZ: 4 A - Mixed Humid and 5 A - Cold Humid).

The change between future 20-year periods (2050s and 2090s) compared to the present period (2010s) is presented for the mean temperature, mean solar irradiance, mean wind speed, and mean relative humidity in Fig. 5.

Table 10 highlights changes at the 99 percentiles of the multi-year distributions. A sharp increase in temperatures is witnessed, especially in the four hottest cities, with changes up to +5.8 °C by the end of the century (i.e., Sao Paulo). For the solar irradiance, wind speed, and relative humidity, similar trends are observed for the mean values.

**Projected changes in weather variables of typical meteorological years (TMY) for building performance simulations.** Table 11 presents the values of mean temperatures, solar irradiance, wind speed, and relative humidity in the three typical meteorological years (TMY) generated from each 20-year dataset. The projected changes in climate variables in the future TMYs are generally consistent with those resulting from the comparison of the 20-year datasets. This means that the TMYs are indeed representative of the climate projections over an interval (i.e., 20 years) and thus suitable for assessing the impact of climate change on building energy loads.

The air temperature is consistently higher in the future weather files for all the cities, with a higher increase in the long-term (2090s) future TMY than in the mid-term (2050s) future TMY. The 2090s-TMY of Abu Dhabi (CZ: 0B Extremely Hot Dry) has the largest increase in temperature of 3.8 °C whereas the TMY of Buenos Aires (CZ: 3 A Warm Humid) has the smallest increase of 1.7 °C for the same period. Many cities are projected to have significantly higher increases in temperature in the long-term than in the mid-term (e.g., Brussels, Ghent, and London). These results are in close agreement with the changes obtained from the 20-year projections. As for the MYs, global solar irradiance will be reduced in the future TMYs of most cities. This is also in agreement with the 20-year projections. The 2090-TMY of Guayaquil (CZ: 1 A Very Hot Humid) has the largest decrease in solar irradiance (23.0 W/m²). The TMYs with slight increases in long-term global solar irradiance are those of Singapore (CZ: 0 A Extremely Hot Humid) and Abu Dhabi (CZ: 0B Extremely Hot Dry). Regarding wind speed, the changes between the 2010s, 2050s, and 2090s weather files are minimal. The 2090s-TMY of Vancouver (CZ: 4 C Mixed Marine) has the largest decrease in mean wind speeds of 0.7 m/s.

Finally, the future TMYs reflect a high variability in the sign of future changes in relative humidity in agreement with the results of the 20-years projections. The cities of Singapore (CZ: 0 A Extremely Hot Humid), Guayaquil (CZ: 1 A Very Hot Humid), Buenos Aires (CZ: 3 A Warm Humid), Los Angeles (CZ: 3B Warm Dry) and Vancouver (CZ: 4 C Mixed Marine) have an absolute increase in relative humidity up to 8% in the 2090-TMYs while Abu Dhabi (CZ: 0B Extremely Hot Dry) has a reduction of relative humidity in the 2090-TMY
Singapore

<figure><img src="image_6.png" /></figure>

London

<figure><img src="image_7.png" /></figure>

Toronto

<figure><img src="image_8.png" /></figure>

Fig. 4 Probability density functions of temperature, wind speed, and relative humidity in Singapore, London, and Toronto from observations (grey), raw RCM (blue), and bias-corrected RCM (red) datasets over the validation time period.

of 2%. The other cities have relatively smaller changes in relative humidity in future TMYs. This variability can be explained by two phenomena. On the one hand, there is general warming, and warmer air can hold more water vapor (air can contain about 7% more moisture for every 1 °C temperature increase according to the Clausius-Clapeyron equation). On the other hand, global warming leads to more evaporation of water and, thus, an increase in specific humidity. Therefore, to keep relative humidity the same, specific humidity must also increase by 7% per °C of warming. However, the oceans are warming more slowly than the land surface, which also means that not enough moisture has evaporated, and relative humidity has, therefore, been reduced.

**Projected changes in heatwaves (HWY) and selected extreme heatwaves for building performance simulations.** Table 12 presents the three thresholds calculated for each city from the 20-year bias-adjusted historical daily temperatures data from 2001 to 2020 for heatwave selection. The relative thresholds are similar for all cities, resulting in different absolute thresholds presented in Table 12. Abu Dhabi is the city with the highest daily mean temperatures. The three European cities in CZ 4 A have equivalent thresholds. For the colder climate zones 5 A and 6 A, Toronto and Montreal in the eastern of Canada have similar thresholds, while European cities Copenhagen and Stockholm also have similar thresholds.

Table 12 also presents the evolution in the number of heatwaves found during each 20-year period. While between 3 and 10 heatwaves are found during the historical period, depending on the cities, a substantial increase in heatwave numbers in the future will be observed in all cities. By 2050, the increase is more pronounced in cities in the four hottest climate zones, followed by cities in North America and then in Europe. Still, in the twenty-year period, every city displays at least one heatwave per summer on average by the mid-century. By the end of the century, the three cities in hot-humid climate zones (Singapore, Guayaquil, and Sao Paulo) showcase an impressive number of heatwaves, beyond a hundred, which would be equivalent to an average of five heatwaves per summer. In these cities, due to the large increase in temperatures, the heatwaves thresholds are exceeded many times during the same summer.

An illustration of the selection of the extreme heatwaves (the most intense, the most severe, and the longest of each period) is made in Fig. 6 for the city of Los Angeles. A bubble represents a heatwave, which size is linked
<table><caption>**Table 9.** 20-year mean temperatures, solar irradiance, wind speed, and relative humidity in the cities over the 2010s, 2050s, and 2090s time periods obtained from multi-year bias-corrected RCM data. Cells with future projected increases (decreases) in climate variables are highlighted in red (green). Grey color means no change. Values in brackets represent the change (absolute value for temperature, solar irradiance, and relative humidity, relative change for wind speed) between the selected term and the 2010s.</caption><thead><tr><th rowspan="2">CZ</th><th rowspan="2">City</th><th colspan="3">Temperature (°C)</th><th colspan="3">Solar irradiance (W/m²)</th><th colspan="3">Wind speed (m/s)</th><th colspan="3">Relative humidity (%)</th></tr><tr><th>2010s</th><th>2050s</th><th>2090s</th><th>2010s</th><th>2050s</th><th>2090s</th><th>2010s</th><th>2050s</th><th>2090s</th><th>2010s</th><th>2050s</th><th>2090s</th></tr></thead><tbody><tr><td>0A</td><td>Singapore</td><td>27.9</td><td>29.1 (1.2)</td><td>30.3 (2.4)</td><td>168.6</td><td>169.0 (0.4)</td><td>166.2 (−2.4)</td><td>2.0</td><td>2.0 (0.0%)</td><td>2.0 (0.0%)</td><td>83.2</td><td>84.7 (1.5)</td><td>84.6 (1.4)</td></tr><tr><td>0B</td><td>Abu Dhabi</td><td>27.7</td><td>29.3 (1.6)</td><td>31.9 (4.2)</td><td>240.3</td><td>241.1 (0.8)</td><td>240.9 (0.6)</td><td>3.2</td><td>3.2 (0.0%)</td><td>3.1 (−3.1%)</td><td>60.0</td><td>58.8 (−1.2)</td><td>57.3 (−2.7)</td></tr><tr><td>1A</td><td>Guayaquil</td><td>26.9</td><td>28.1 (1.2)</td><td>30.2 (3.3)</td><td>263.0</td><td>258.1 (−4.9)</td><td>253.5 (−9.5)</td><td>1.7</td><td>1.7 (0.0%)</td><td>1.7 (0.0%)</td><td>74.9</td><td>75.2 (0.3)</td><td>75.7 (0.8)</td></tr><tr><td>2A</td><td>Sao Paulo</td><td>19.8</td><td>21.3 (1.5)</td><td>23.3 (3.5)</td><td>188.8</td><td>188.5 (−0.3)</td><td>182.8 (−6.0)</td><td>6.1</td><td>6.0 (−1.6%)</td><td>6.0 (−1.6%)</td><td>80.3</td><td>80.0 (−0.3)</td><td>80.3 (0.0)</td></tr><tr><td>3A</td><td>Buenos Aires</td><td>17.9</td><td>18.7 (0.8)</td><td>19.5 (1.6)</td><td>188.4</td><td>184.5 (−3.9)</td><td>178.8 (−9.6)</td><td>4.3</td><td>4.2 (−2.3%)</td><td>3.9 (−9.3%)</td><td>73.5</td><td>75.7 (2.2)</td><td>78.3 (4.8)</td></tr><tr><td>3A</td><td>Rome</td><td>16.1</td><td>17.1 (1.0)</td><td>19.6 (2.5)</td><td>187.7</td><td>183.7 (−4.0)</td><td>185.7 (−2.0)</td><td>3.6</td><td>3.6 (0.0%)</td><td>3.5 (−2.8%)</td><td>72.1</td><td>73.3 (1.2)</td><td>71.4 (−0.7)</td></tr><tr><td>3B</td><td>Los Angeles</td><td>16.7</td><td>17.9 (1.2)</td><td>19.4 (2.7)</td><td>214.7</td><td>211.2 (−3.5)</td><td>205.9 (−8.8)</td><td>1.7</td><td>1.6 (−5.9%)</td><td>1.6 (−5.9%)</td><td>72.4</td><td>74.6 (2.2)</td><td>77.1 (4.7)</td></tr><tr><td>4A</td><td>Brussels</td><td>10.8</td><td>11.6 (0.8)</td><td>13.4 (2.6)</td><td>126.2</td><td>123.3 (−2.9)</td><td>118.3 (−7.9)</td><td>3.6</td><td>3.6 (0.0%)</td><td>3.6 (0.0%)</td><td>78.6</td><td>78.6 (0.0)</td><td>78.7 (0.1)</td></tr><tr><td>4A</td><td>Ghent</td><td>11.0</td><td>11.7 (0.7)</td><td>13.6 (2.6)</td><td>108.1</td><td>105.9 (−2.2)</td><td>101.7 (−6.4)</td><td>4.2</td><td>4.2 (0.0%)</td><td>4.2 (0.0%)</td><td>83.1</td><td>83.1 (0.0)</td><td>83.1 (0.0)</td></tr><tr><td>4A</td><td>London</td><td>12.0</td><td>12.7 (0.7)</td><td>14.5 (2.5)</td><td>118.4</td><td>115.2 (−3.2)</td><td>113.1 (−5.3)</td><td>4.0</td><td>4.0 (0.0%)</td><td>4.0 (0.0%)</td><td>75.1</td><td>75.3 (0.2)</td><td>74.9 (−0.2)</td></tr><tr><td>4C</td><td>Vancouver</td><td>7.8</td><td>9.1 (1.3)</td><td>10.9 (3.1)</td><td>153.8</td><td>149.4 (−4.4)</td><td>142.9 (−10.9)</td><td>4.4</td><td>4.3 (−2.3%)</td><td>4.1 (−6.8%)</td><td>69.2</td><td>69.5 (0.3)</td><td>70.4 (1.2)</td></tr><tr><td>5A</td><td>Toronto</td><td>7.9</td><td>8.9 (1.0)</td><td>11.1 (3.2)</td><td>153.7</td><td>153.4 (−0.3)</td><td>149.6 (−4.1)</td><td>4.4</td><td>4.3 (−2.3%)</td><td>4.2 (−4.6%)</td><td>68.9</td><td>68.7 (−0.2)</td><td>69.2 (0.3)</td></tr><tr><td>5A</td><td>Copenhagen</td><td>8.8</td><td>9.7 (0.9)</td><td>11.3 (2.5)</td><td>117.8</td><td>114.9 (−2.9)</td><td>108.0 (−9.8)</td><td>3.3</td><td>3.2 (−3.0%)</td><td>3.3 (0.0%)</td><td>82.4</td><td>82.6 (0.2)</td><td>83.0 (0.6)</td></tr><tr><td>6A</td><td>Montreal</td><td>7.9</td><td>8.9 (1.0)</td><td>10.9 (3.0)</td><td>154.3</td><td>156.8 (2.5)</td><td>152.0 (−2.3)</td><td>4.4</td><td>4.4 (0.0%)</td><td>4.2 (−4.6%)</td><td>69.0</td><td>68.6 (−0.4)</td><td>69.0 (0.0)</td></tr><tr><td>6A</td><td>Stockholm</td><td>7.7</td><td>8.9 (1.2)</td><td>10.6 (2.9)</td><td>116.4</td><td>110.5 (−5.9)</td><td>103.6 (−12.8)</td><td>3.8</td><td>4.0 (5.3%)</td><td>4.0 (5.3%)</td><td>79.1</td><td>78.9 (−0.2)</td><td>79.1 (0.0)</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Fig. 5** Changes in climatic variables from the 2010s to 2050s and 2090s: (a) absolute change for temperature, (b) relative change in wind speed, (c) absolute change in solar radiation, (d) absolute change in relative humidity.</figcaption></figure>

to its severity. Figure 6a illustrates well the tremendous increase throughout the century and the diversity of heatwaves that are found as well. In comparison with the historical period, during which only very short heatwaves of five days are witnessed, in the mid-term future, longer heatwaves that are both less or more intense
<table><caption>**Table 10.** 20-year 99% temperatures, solar irradiance, wind speed, and relative humidity in the cities over the 2010s, 2050s, and 2090s time periods obtained from multi-year bias-corrected RCM data. Cells with future projected increases (decreases) in climate variables are highlighted in red (green). Grey color means no change.</caption><thead><tr><th rowspan="2">CZ</th><th rowspan="2">City</th><th colspan="3">Temperature (°C)</th><th colspan="3">Solar irradiance (W/m²)</th><th colspan="3">Wind speed (m/s)</th><th colspan="3">Relative humidity (%)</th></tr><tr><th>2010s</th><th>2050s</th><th>2090s</th><th>2010s</th><th>2050s</th><th>2090s</th><th>2010s</th><th>2050s</th><th>2090s</th><th>2010s</th><th>2050s</th><th>2090s</th></tr></thead><tbody><tr><td>0A</td><td>Singapore</td><td>33.0</td><td>34.4 (1.4)</td><td>38.2 (5.2)</td><td>965.1</td><td>960.0 (–5.1)</td><td>949.4 (–15.7)</td><td>6.2</td><td>6.2 (0.0%)</td><td>6.0 (–3.3%)</td><td>99.4</td><td>100 (0.6)</td><td>100 (0.6)</td></tr><tr><td>0B</td><td>Abu Dhabi</td><td>44.3</td><td>46.4 (2.1)</td><td>49.5 (5.2)</td><td>940.3</td><td>934.9 (–5.4)</td><td>927.3 (–13.0)</td><td>9.2</td><td>9.1 (–1.1%)</td><td>9.0 (–2.2%)</td><td>95.1</td><td>94.8 (–0.3)</td><td>94.4 (–0.7)</td></tr><tr><td>1A</td><td>Guayaquil</td><td>34.6</td><td>35.7 (1.1)</td><td>37.8 (3.2)</td><td>1,296.5</td><td>1,281.0 (–15.5)</td><td>1,260.3 (–36.2)</td><td>4.7</td><td>4.7 (0.0%)</td><td>4.5 (–4.4%)</td><td>99.0</td><td>99.0 (0.0)</td><td>99.0 (0.0)</td></tr><tr><td>2A</td><td>Sao Paulo</td><td>31.0</td><td>33.6 (2.6)</td><td>36.8 (5.8)</td><td>895.6</td><td>891.4 (–4.2)</td><td>881.7 (–13.9)</td><td>13.9</td><td>14.1 (1.4%)</td><td>14.5 (4.1%)</td><td>100</td><td>100 (0.0)</td><td>97.7 (–2.3)</td></tr><tr><td>3A</td><td>Buenos Aires</td><td>29.5</td><td>30.0 (0.5)</td><td>30.6 (1.1)</td><td>1,011.0</td><td>1,002.8 (–8.2)</td><td>994.7 (–16.3)</td><td>11.1</td><td>10.9 (–1.8%)</td><td>10.4 (–6.7%)</td><td>100</td><td>100 (0.0)</td><td>100 (0.0)</td></tr><tr><td>3A</td><td>Rome</td><td>30.6</td><td>32.0 (1.4)</td><td>35.9 (5.3)</td><td>906.0</td><td>898.8 (–9.2)</td><td>892.0 (–14)</td><td>10.5</td><td>10.5 (0.0%)</td><td>10.4 (–1.0%)</td><td>100</td><td>100 (0.0)</td><td>100 (0.0)</td></tr><tr><td>3B</td><td>Los Angeles</td><td>28.1</td><td>29.8 (1.7)</td><td>31.0 (2.9)</td><td>956.0</td><td>946.1 (–9.9)</td><td>934.9 (–21.1)</td><td>4.2</td><td>4.2 (0.0%)</td><td>4.2 (0.0%)</td><td>100</td><td>100 (0.0)</td><td>100 (0.0)</td></tr><tr><td>4A</td><td>Brussels</td><td>26.4</td><td>28.0 (1.6)</td><td>29.3 (2.9)</td><td>759.0</td><td>749.7 (–9.3)</td><td>733.1 (–25.9)</td><td>8.7</td><td>9.0 (3.3%)</td><td>9.0 (3.3%)</td><td>98.9</td><td>98.8 (–0.1)</td><td>98.8 (–0.1)</td></tr><tr><td>4A</td><td>Ghent</td><td>26.9</td><td>28.1 (1.2)</td><td>29.5 (2.6)</td><td>779.4</td><td>771.6 (–7.8)</td><td>753.8 (–25.6)</td><td>9.6</td><td>9.8 (2.0%)</td><td>9.9 (3.0%)</td><td>100</td><td>100 (0.0)</td><td>100 (0.0)</td></tr><tr><td>4A</td><td>London</td><td>26.7</td><td>27.8 (1.1)</td><td>29.9 (3.2)</td><td>779.4</td><td>770.6 (–8.8)</td><td>763.6 (–15.8)</td><td>10.5</td><td>10.7 (1.9%)</td><td>10.6 (0.9%)</td><td>98.0</td><td>98.0 (0.0)</td><td>98.0 (0.0)</td></tr><tr><td>4C</td><td>Vancouver</td><td>28.9</td><td>30.4 (1.5)</td><td>31.7 (2.8)</td><td>882.5</td><td>872.7 (–9.8)</td><td>858.6 (–23.9)</td><td>11.8</td><td>11.6 (–1.7%)</td><td>11.5 (–2.6%)</td><td>98.3</td><td>98.1 (–0.2)</td><td>98.1 (–0.2)</td></tr><tr><td>5A</td><td>Toronto</td><td>29.0</td><td>30.7 (1.7)</td><td>32.9 (3.9)</td><td>882.4</td><td>879.6 (–2.8)</td><td>863.6 (–18.8)</td><td>11.9</td><td>11.7 (–1.7%)</td><td>11.6 (–2.6%)</td><td>98.1</td><td>98.2 (0.1)</td><td>98.2 (0.1)</td></tr><tr><td>5A</td><td>Copenhagen</td><td>24.7</td><td>25.7 (1.0)</td><td>26.4 (1.7)</td><td>779.1</td><td>770.9 (–8.2)</td><td>757.8 (–21.3)</td><td>9.2</td><td>9.2 (0.0%)</td><td>9.3 (1.1%)</td><td>100</td><td>100 (0.0)</td><td>100 (0.0)</td></tr><tr><td>6A</td><td>Montreal</td><td>29.0</td><td>30.0 (1.0)</td><td>31.8 (2.8)</td><td>883.6</td><td>881.0 (–2.6)</td><td>872.8 (–10.8)</td><td>11.9</td><td>11.7 (–1.7%)</td><td>11.4 (–4.4%)</td><td>97.8</td><td>97.9 (0.1)</td><td>97.3 (–0.5)</td></tr><tr><td>6A</td><td>Stockholm</td><td>24.2</td><td>24.9 (0.7)</td><td>25.6 (1.4)</td><td>699.0</td><td>696.3 (–2.7)</td><td>682.7 (–16.3)</td><td>8.4</td><td>8.7 (3.4%)</td><td>8.8 (4.5%)</td><td>99.0</td><td>99.0 (0.0)</td><td>98.9 (–0.1)</td></tr></tbody></table>

<table><caption>**Table 11.** Mean temperatures, solar irradiance, wind speed, and relative humidity in the three TMYs weather files generated based on the bias-corrected 20-years datasets for each city.</caption><thead><tr><th rowspan="2">CZ</th><th rowspan="2">City</th><th colspan="3">Temperature (°C)</th><th colspan="3">Solar irradiance (W/m²)</th><th colspan="3">Wind speed (m/s)</th><th colspan="3">Relative humidity (%)</th></tr><tr><th>2010s</th><th>2050s</th><th>2090s</th><th>2010s</th><th>2050s</th><th>2090s</th><th>2010s</th><th>2050s</th><th>2090s</th><th>2010s</th><th>2050s</th><th>2090s</th></tr></thead><tbody><tr><td>0A</td><td>Singapore</td><td>27.9</td><td>29.1</td><td>30.3</td><td>163.7</td><td>167.5</td><td>165.6</td><td>2.1</td><td>2.0</td><td>2.0</td><td>82.7</td><td>84.8</td><td>84.8</td></tr><tr><td>0B</td><td>Abu Dhabi</td><td>27.9</td><td>29.4</td><td>31.7</td><td>233.9</td><td>235.5</td><td>234.7</td><td>3.2</td><td>3.2</td><td>3.1</td><td>59.7</td><td>58.4</td><td>57.7</td></tr><tr><td>1A</td><td>Guayaquil</td><td>27.1</td><td>28.5</td><td>30.4</td><td>255.5</td><td>243.4</td><td>232.5</td><td>1.7</td><td>1.8</td><td>1.6</td><td>73.4</td><td>71.8</td><td>75.8</td></tr><tr><td>2A</td><td>Sao Paulo</td><td>19.8</td><td>21.3</td><td>23.1</td><td>190.9</td><td>193.3</td><td>183.8</td><td>6.0</td><td>6.0</td><td>6.0</td><td>80.3</td><td>80.6</td><td>80.6</td></tr><tr><td>3A</td><td>Buenos Aires</td><td>17.8</td><td>18.9</td><td>19.5</td><td>192.8</td><td>190.1</td><td>182.9</td><td>4.3</td><td>4.1</td><td>3.9</td><td>73.5</td><td>74.9</td><td>78.5</td></tr><tr><td>3A</td><td>Rome</td><td>16.1</td><td>17.2</td><td>19.5</td><td>189.4</td><td>185.6</td><td>187.2</td><td>3.6</td><td>3.6</td><td>3.6</td><td>73.2</td><td>72.5</td><td>70.5</td></tr><tr><td>3B</td><td>Los Angeles</td><td>16.7</td><td>17.8</td><td>19.3</td><td>219.8</td><td>210.1</td><td>206.1</td><td>1.7</td><td>1.6</td><td>1.6</td><td>71.4</td><td>75.5</td><td>79.0</td></tr><tr><td>4A</td><td>Brussels</td><td>11.2</td><td>11.4</td><td>13.5</td><td>124.9</td><td>119.9</td><td>118.8</td><td>3.7</td><td>3.7</td><td>3.7</td><td>79.5</td><td>78.6</td><td>78.4</td></tr><tr><td>4A</td><td>Ghent</td><td>10.9</td><td>11.5</td><td>13.3</td><td>124.5</td><td>119.0</td><td>118.5</td><td>3.4</td><td>3.5</td><td>3.4</td><td>79.0</td><td>79.6</td><td>79.2</td></tr><tr><td>4A</td><td>London</td><td>12.1</td><td>12.7</td><td>14.2</td><td>116.9</td><td>114.4</td><td>110.6</td><td>4.1</td><td>4.0</td><td>4.0</td><td>75.1</td><td>75.2</td><td>76.0</td></tr><tr><td>4C</td><td>Vancouver</td><td>7.7</td><td>9.3</td><td>10.7</td><td>156.3</td><td>147.2</td><td>147.8</td><td>4.7</td><td>4.2</td><td>4.0</td><td>67.6</td><td>71.5</td><td>76.0</td></tr><tr><td>5A</td><td>Toronto</td><td>8.2</td><td>9.3</td><td>11.4</td><td>155.1</td><td>155.4</td><td>149.8</td><td>4.3</td><td>4.4</td><td>4.2</td><td>71.0</td><td>67.4</td><td>69.5</td></tr><tr><td>5A</td><td>Copenhagen</td><td>9.0</td><td>9.7</td><td>11.2</td><td>119.8</td><td>113.1</td><td>107.4</td><td>3.3</td><td>3.2</td><td>3.4</td><td>81.2</td><td>83.3</td><td>83.0</td></tr><tr><td>6A</td><td>Montreal</td><td>7.9</td><td>9.0</td><td>10.8</td><td>153.9</td><td>156.5</td><td>149.0</td><td>4.3</td><td>4.7</td><td>4.4</td><td>69.8</td><td>70.1</td><td>69.9</td></tr><tr><td>6A</td><td>Stockholm</td><td>7.9</td><td>8.9</td><td>10.7</td><td>119.5</td><td>111.9</td><td>108.8</td><td>3.8</td><td>3.9</td><td>4.0</td><td>79.7</td><td>79.7</td><td>78.4</td></tr></tbody></table>

than the most intense heatwave of the historical period are found. By the end of the century, heatwaves are more severe and also longer. In Fig. 6b, the three most extreme heatwaves, the ones that are selected for future periods, are highlighted. During the 2050s:

* the most intense heatwave is 8 days long with an intensity of 30.1 °C and a severity of 14.2 °C.d;

* the most severe heatwave is 21 days long with an intensity of 27.6 °C and a severity of 32.4 °C.d;

* the longest heatwave is 22 days long with an intensity of 25.6 °C and a severity of 17.5 °C.d;

During the 2100s, only two extreme heatwaves are selected:

* the most intense, which is also the most severe: intensity of 32.1 °C, duration of 16 days, and severity of 39.1 °C.d;

* the longest, which is 38 days long with an intensity of 28.7 °C and a severity of 53.6 °C.d

For each city, the three extreme heatwaves (the most intense, most severe, and longest heatwave) are selected. Figure 7 shows the characteristics (intensity, severity, and duration) of the most intense and longest heat waves
<table><caption>**Table 12.** Thresholds used over the historical period 2010s (2001–2020) for heatwave selection and number of heatwaves found per period in each city.</caption><thead><tr><th rowspan="2">CZ</th><th rowspan="2">City</th><th rowspan="2">CORDEX Domain</th><th colspan="3">Threshold to detect heatwaves over 2010s</th><th colspan="3">Number of heatwaves detected</th></tr><tr><th>95 Threshold (°C)</th><th>97.5 Threshold (°C)</th><th>99.5 Threshold (°C)</th><th>2010s</th><th>2050s</th><th>2090s</th></tr></thead><tbody><tr><td>0A</td><td>Singapore</td><td>SEA</td><td>30.4</td><td>30.9</td><td>31.7</td><td>7</td><td>58</td><td>136</td></tr><tr><td>0B</td><td>Abu Dhabi</td><td>SEA</td><td>37.1</td><td>38.1</td><td>39.3</td><td>5</td><td>47</td><td>61</td></tr><tr><td>1A</td><td>Guayaquil</td><td>AFR</td><td>29.3</td><td>29.8</td><td>30.8</td><td>8</td><td>40</td><td>207</td></tr><tr><td>2A</td><td>Sao Paulo</td><td>SAM</td><td>25.3</td><td>26.3</td><td>28.3</td><td>7</td><td>87</td><td>172</td></tr><tr><td>3A</td><td>Buenos Aires</td><td>SAM</td><td>25.3</td><td>26.4</td><td>28.0</td><td>6</td><td>19</td><td>32</td></tr><tr><td>3A</td><td>Rome</td><td>EUR</td><td>25.6</td><td>26.4</td><td>27.9</td><td>7</td><td>21</td><td>36</td></tr><tr><td>3B</td><td>Los Angeles</td><td>NAM</td><td>22.0</td><td>22.8</td><td>24.4</td><td>3</td><td>40</td><td>81</td></tr><tr><td>4A</td><td>London</td><td>EUR</td><td>20.6</td><td>21.9</td><td>24.3</td><td>6</td><td>16</td><td>46</td></tr><tr><td>4A</td><td>Brussels</td><td>EUR</td><td>20.4</td><td>22.0</td><td>24.6</td><td>9</td><td>14</td><td>36</td></tr><tr><td>4A</td><td>Ghent</td><td>EUR</td><td>20.2</td><td>21.8</td><td>25.0</td><td>7</td><td>20</td><td>33</td></tr><tr><td>4C</td><td>Vancouver</td><td>NAM</td><td>24.0</td><td>25.3</td><td>27.5</td><td>8</td><td>23</td><td>54</td></tr><tr><td>5A</td><td>Toronto</td><td>NAM</td><td>23.3</td><td>24.2</td><td>25.7</td><td>4</td><td>39</td><td>85</td></tr><tr><td>5A</td><td>Copenhagen</td><td>EUR</td><td>18.7</td><td>20.1</td><td>22.3</td><td>10</td><td>19</td><td>23</td></tr><tr><td>6A</td><td>Montreal</td><td>NAM</td><td>23.3</td><td>24.1</td><td>25.5</td><td>4</td><td>38</td><td>88</td></tr><tr><td>6A</td><td>Stockholm</td><td>EUR</td><td>18.6</td><td>19.8</td><td>21.7</td><td>9</td><td>25</td><td>27</td></tr></tbody></table>

in each climate zone. Characteristics of the most severe heatwaves are often similar to the longest heatwaves and are not shown here.

The intensity of both extreme heatwaves strongly increases between the three periods and in each climate zone. The increase in intensity of the most intense heatwave by the end of the century is, in each city, between +2°C (European cities in climate zone 4A) and +7°C (Sao Paulo). The intensity of the longest heatwaves is between 0°C and 3.4°C (Vancouver), inferior to the most intense heatwaves in the 2010s, of 0°C and 4.5°C (Los Angeles) inferior in the 2050s, and of 0.2 and 3.8°C (Vancouver) by the 2100s.

The extreme heatwaves' durations strongly increase between the three time periods, especially the one of the longest heatwaves. The increase is more pronounced between 2100s and 2050s than between 2050s and 2010s. By the 2010s, the duration of both the most intense and longest extreme heat waves is generally around one to three weeks, depending on the city. However, by 2050s, the extreme heatwaves last more than a month in Abu Dhabi (41 days), Rome (42 days), and Stockholm (37 days), between 6 and 24 days for the most intense heatwaves in the other cities, between 7 and 49 days for the longest heatwaves in the other cities. By the 2100s, in the five hottest cities (from climate zones 0A, 0B, 1A, 2A, and 3A), the longest and the most intense heatwaves last 3 to 4 months. This high number is found because the temperatures will constantly be above the current thresholds during the hot period of the year. In other parts of the world, the longest heatwave will be between three weeks and 2 months long by the 2100s, except in Buenos Aires. For climate zone 3A, the severity and duration of the heat waves in Rome are more significant than in Buenos Aires. This disparity might be attributed to the heatwave data record, which shows European cities have more exposure to heatwaves¹⁹. As expected, we observe that the durations of extreme intense heatwaves are generally shorter than the longest heatwaves.

**Effect of future TMY and HWY weather files on building performance.** Lee and Levinson⁶² evaluated the effect of cool envelope strategies on heating, ventilation, and air conditioning (HVAC) primary energy use intensity and thermal comfort for a mechanically cooled single-family home in Los Angeles in Fig. 8. They used the future TMYs produced based on the methodology introduced in this paper (named CORDEX 2010, 2050, and 2090) as well as the historical Typical Meteorological Year 3 (TMY3), which spans 1991–2005⁶³. Panel A shows that cooling demand grows over time. They also calculated the thermal sensation scale unit (TSSU) weighted warm discomfort exceedance hours (TSSU-h) to evaluate the Predicted Mean Vote (PMV) based thermal comfort, which is the sum of summer thermal discomfort when PMV exceeds +0.7 according to ISO 1772-2:2018⁶⁴. PMV greater than +0.7 is considered uncomfortably warm during the summer season according to ISO 17772-1:2017 Annex H.1 Category III⁶⁵. Annex H.1 Category III, considered uncomfortably warm during the summer season. Panel B shows that the occupants experience many more TSSU-weighted warm-discomfort exceedance hours in the future because the cooling system is sized based on historical TMY3 weather, which results in many hours during which the cooling system cannot meet future loads. They also show that use of passive strategies such as cool envelope materials, helps decrease these loads. These results emphasize the need to use future TMYs to anticipate an increase in cooling energy use intensity and take necessary action to adapt building design or refurbishment to future climate.

Another example of how these weather files can be used is the work of Sengupta *et al.*⁶⁶ in which they evaluated the overheating of an educational building in Ghent, Belgium, under future weather files, comparing the results with the future TMY and HWY prepared in this paper. Educational buildings in Belgium are not equipped with mechanical air conditioning, and recent heat waves have already posed a threat to occupants' cognitive performance and health conditions. In their paper, they analyzed the thermal resilience of test lecture rooms with open windows at night for natural ventilation to flush heat and equipped with indirect evaporative cooling to cool the air during the daytime. Figure 9 shows the results of unmet degree hours (UDH) for different
<figure><img src="image_10.png" /><figcaption>Fig. 6 Heatwaves in Los Angeles (CZ 3B): (a) All heatwaves detected and (b) extreme heatwaves selection.</figcaption></figure>

weather files: a) TMY and b) HWY (1 A: 2010s intense HW, 1B: 2010s severe and longest HW, 2 A: 2050s intense HW, 2B: 2050s severe HW, 2 C: 2050s longest HW, 3 A: 2100s intense HW, 3B: 2100s severe and longest HW) with and without power outage (PO). The results emphasize that the HWYs present a much larger number of UDH when compared to TMY. The variety of HWY shows that HWY 1B leads to many UDHs due to its length of 28 days, while HWs of the 2100s also predict a very elevated number of UDHs due to the increase in outdoor temperatures (Fig. 7). Additionally, a study by Sengupta *et al.*<sup>67</sup> identifying, quantifying, and comparing different shocks that can increase overheating risk in buildings (e.g., outdoor shocks such as heatwaves and mechanical shocks such as solar shading failure, cooling strategy failure, natural night ventilation failure) proves that heatwaves are by far the most intense shocks for buildings that impact the thermal resilience to overheating. Thus, assessing and improving the buildings' performance against heatwaves are a crucial step to future proof these buildings, emphasizing the robust methodology needed to develop and utilize future weather data and heatwave data to assess and design buildings.

## Usage Notes

The provided typical meteorological years (TMY) and years containing heatwaves (HWY) in both EPW format are ready-to-use weather datasets to perform building performance simulations using Energy Plus, TRNSYS, or any other building energy model. They permit assessment of the thermal performance of buildings under typical and extreme future climate scenarios. Therefore, they help evaluate the efficiency and resilience of building renovation solutions to climate change in different climate zones<sup>68</sup>. In particular, the TMY can be used to analyze changes in building heating and cooling loads under typical future weather conditions. The HWYs allow prediction of building thermal response under extreme heat events, which will be one major issue in the next decades. The multi-year (MY) datasets are also provided in CSV format to allow other authors to test different methods for assembling different types of future typical or extreme weather files for building performance assessments or in other sectors.

The provided datasets were generated based on the bias-corrected climate model MPI-ESM-LR/REMO, whose temperature projections are found to be the closest to the median of all climate model projections<sup>69</sup>. At least two other GCM/RCM model combinations satisfy the required spatial and temporal resolutions in the CORDEX database to generate weather files for building thermal performance analysis. These are the HadGEM2-ES/REMO and the NorESM1-M/REMO. Therefore, the results of this paper can be further expanded by comparing the outputs of all available CORDEX models. This can be used in future work to enrich the datasets. The datasets were generated based on RCP 8.5 climate projections, the worst-case socioeconomic scenario at the time of the IPCC AR5, and the most realistic based on the past and current emissions of greenhouse gases by the global community<sup>70</sup>. This means that they are suitable for applications in studies of system resilience, but they should be used with caution in building retrofits and HVAC system designs to avoid system oversizing or under-sizing.

It is possible to assemble additional weather files for other cities worldwide using other climate models with the same methodology as provided in this paper. The Python code to assemble the datasets from CORDEX climate projections is provided in the section “Code Availability.” Additional weather variables, such as cloud cover, precipitation, and longwave solar irradiance, would be an added value in the datasets. However, these additional variables were not currently available for all cities, neither observation. Indeed, robust climate observations are needed, and in this study, for some cities, only a few years (<10 years) of observations for the bias correction were available, which can affect the final result to an important extent for those cities. Future climate projections available on the CORDEX platform, with the newest SSP scenarios from CMIP6, might allow the include additional climate data in the datasets. In that case, observational data for these specific variables must also be found to correct the model.

Given that the multi-year datasets are provided, they could be used to select heatwaves based on methods different from the one chosen here. The common method used detects heatwaves solely based on the temperature; however, in some hot and humid parts of the world, humidity is known to be an important variable affecting indoor heat stress. The simple method proposed here was validated for several cities in France and allows a standardized approach that fits the purpose of a common method for all cities, climates, building typologies, and other local specificities. Nevertheless, the multiyear datasets allow the use of additional criteria to select the heatwaves with different methods. Beyond a different method, less future extreme heatwaves could also be selected for building design<sup>4</sup>.

As explained in the “boundary conditions” section, the datasets do not incorporate urban effects. In the selected GCM/RCM-REMO model, urban areas are represented as simple impervious surfaces. Recent studies
<figure><img src="image_11.png" /><figcaption>**Fig. 7** Characteristics (intensity, severity, and duration) of the most intense and longest XTRM-HW: (a) intensity of the most intense HW, (b) intensity of the longest HW, (c) severity of the most intense HW, (d) severity of the longest HW, (e) duration of the most intense HW, (f) duration of the longest HW.</figcaption></figure>

have shown that a more detailed urban parametrization allows a better understanding of the regional-urban climate interactions and urban climate effects, such as UHI intensity<sup>71–74</sup>. However, this entails a significant increase in computing power and time, limiting the analysis to shorter time periods. Due to such limitations in modelling urban areas, the RCM REMO model does not accurately simulate climate modifications induced by urban features, such as the urban heat island effect or urban microclimates. Accordingly, bias-correction of the model projections was performed using observational data from weather stations located outside cities. The urban heat island effect and other urban climate modifications can be added to the weather datasets following
Baseline annual HVAC primary energy use intensity

<figure><img src="image_12.png" /></figure>

(A)

Baseline annual PMV-PPD model-based TSSU-weighted warm discomfort exceedance hours

<figure><img src="image_13.png" /></figure>

(B)

<figure><img src="image_14.png" /><figcaption>**Fig. 8** Effect of future TMYs on energy use (A) and summer thermal discomfort (B) in an air-conditioned single-family home in Los Angeles (Lee and Levinson<sup>62</sup>).</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Fig. 9** (a) Impact of future TMYs and (b) Impact of future HWYs on summer thermal discomfort from Sengupta *et al.*<sup>66</sup>.</figcaption></figure>

different methodologies already proposed in building performance simulation studies<sup>75-78</sup>. Most climate models do not explicitly model urban areas and, at best, describe them as rock covers. Nonetheless, the very high resolutions reached now by the regional climate models may justify and require a more realistic parameterization of surface exchanges between urban canopy and atmosphere.

To quantify the potential impact of urbanization on the regional climate and evaluate the benefits of a detailed urban canopy model compared with a simpler approach, a sensitivity study was carried out over France at a 12 km horizontal resolution with the ALADIN-Climate regional model for 1980–2009 time period. Different descriptions of land use and urban modeling were compared, corresponding to an explicit modeling of cities with the urban canopy model TEB, a conventional and simpler approach representing urban areas as rocks, and a vegetated experiment for which cities are replaced by natural covers. A general evaluation of ALADIN-Climate was first done, which showed an overestimation of the incoming solar irradiance but satisfying results in terms of precipitation and near-surface temperatures. The sensitivity analysis then highlighted those urban areas had a significant impact on modeled near-surface temperature. A further analysis of a few large French cities indicated that over the 30 years of simulation, they all induced a warming effect both at daytime and nighttime with values
up to +1.5 °C for the city of Paris. The urban model also led to regional warming extending beyond the boundaries of urban areas. Finally, the comparison to temperature observations available for the Paris area highlighted that the detailed urban canopy model improved the modeling of the urban heat island compared with a simpler approach.

The urban heat island effect could be added to the weather datasets by using offline urban canopy tools like the Urban Weather Generator (UWG)⁷⁹,⁸⁰, the Surface Urban Energy and Water Balance Scheme (SUEWS)⁸¹ module of the Urban Multi-scale Environmental Predictor (UMEP) GIS tool, or other similar urban canopy models⁸². Urban canopy models can also be coupled with mesoscale models such as the Weather Research and Forecasting (WRF) Model⁸³,⁸⁴ or the Global Environmental Multi-scale (GEM) Model⁸⁵ for a better consideration of the urban boundary layer conditions⁸⁶. The UWG is an easy-to-use, computational inexpensive tool that directly outputs urban weather files. However, it assumes that the city's urban fabric is homogeneous and that the city is surrounded by rural areas. This can make its results inaccurate for coastal cities or inhomogeneous urban fabrics⁸⁷,⁸⁸. UWG accuracy may also be limited by the simplified ways in which it calculates latent heat balance flux and urban canyon wind speed⁸⁰. Recently, new stand-alone UCM models have been developed that overcome some of the UWG limitations, such as the Stand-alone Urban Energy/Climate Model (SUECM)⁸⁹. City Fast Fluid Dynamics (CityFFD)⁹⁰ and the Vertical-city Weather Generator (VCWG)⁹¹. Machine learning techniques were also used to interpolate weather data spatially⁹²⁻⁹⁴. Any of these tools can be used to add urban effects as well as the evolution of land use to both the present and future TMYs and heatwave weather files presented in this data paper.

## Code availability

The source codes to generate these datasets from CORDEX climate data can be found at: https://zenodo.org/record/7300024#.ZBbi4XbMI2x⁹⁵.

Received: 18 April 2023; Accepted: 29 April 2024;

Published online: 23 May 2024

## References

1. Belcher, S., Hacker, J. & Powell, D. Constructing design weather data for future climates. *Building Services Engineering Research and Technology* **26**, 49–61 (2005).

2. Nik, V. M. Making energy simulation easier for future climate – Synthesizing typical and extreme weather data sets out of regional climate models (RCMs). *Applied Energy* **177**, 204–226 (2016).

3. Machard, A., Inard, C., Alessandrini, J.-M., Pelé, C. & Ribéron, J. A Methodology for Assembling Future Weather Files Including Heatwaves for Building Thermal Simulations from the European Coordinated Regional Downscaling Experiment (EURO-CORDEX) Climate Data. *Energies* **13**, 3424 (2020).

4. Machard, A. Towards mitigation and adaptation to climate change: Contribution to Building Design. (2021).

5. ISO. ISO. ISO 15927-4:2005 Hygrothermal Performance of Buildings — Calculation and Presentation of Climatic Data — Part 4: Hourly Data for Assessing the Annual Energy Use for Heating and Cooling. (2005).

6. Ouzeau, G., Soubeyroux, J.-M., Schneider, M., Vautard, R. & Planton, S. Heat waves analysis over France in present and future climate: Application of a new method on the EURO-CORDEX ensemble. *Climate Services* **4**, 1–12 (2016).

7. Doutreloup, S. et al. Historical and future weather data for dynamic building simulations in Belgium using the regional climate model MAR: typical and extreme meteorological year and heatwaves. *Earth Syst. Sci. Data* **14**, 3039–3051 (2022).

8. Ramon, D., Allacker, K., van Lipzig, N. P. M., De Troyer, F. & Wouters, H. Future Weather Data for Dynamic Building Energy Simulations: Overview of Available Data and Presentation of Newly Derived Data for Belgium. in *Energy Sustainability in Built and Urban Environments* (eds. Motoasca, E., Agarwal, A. K. & Breesch, H.) 111–138. https://doi.org/10.1007/978-981-13-3284-5_6 (Springer Singapore, Singapore, 2019).

9. Ramon, D., Allacker, K., Trigaux, D., Wouters, H. & Van Lipzig, N. P. M. Dynamic modelling of operational energy use in a building LCA: A case study of a Belgian office building. *Energy and Buildings* **278**, 112634 (2023).

10. Ramon, D., Allacker, K., De Troyer, F., Wouters, H. & Van Lipzig, N. P. M. Future heating and cooling degree days for Belgium under a high-end climate change scenario. *Energy and Buildings* **216**, 109935 (2020).

11. Gaur, A. & Lacasse, M. Climate Data to Support the Adaptation of Buildings to Climate Change in Canada. *Data* **7**, 42 (2022).

12. Zou, J. et al. Assessment of future overheating conditions in Canadian cities using a reference year selection method. *Building and Environment* **218**, 109102 (2022).

13. Bass, B., New, J., Rastogi, D. & Kao, S.-C. *Future Typical Meteorological Year (fTMY) US Weather Files for Building Simulation*. https://zenodo.org/record/6939750 (2022).

14. ASHRAE. ANSI/ASHRAE Standard 169-2013, Climatic Data for Building Design. (ASHRAE, Atlanta, 2013).

15. Holzer, P. & Cooper, P. IEA EBC Annex 80 on Resilient Cooling for Residential and Small Non-Residential Buildings. 14 https://annex80.iea-ebc.org/Data/Sites/10/media/documents/supporting/ebc-annex-80-annex-text-190616.pdf (2019).

16. Attia, S. et al. Framework to Evaluate the Resilience of Different Cooling Technologies. https://doi.org/10.13140/RG.2.2.33998.59208 (2021).

17. Rahif, R. et al. Simulation-based framework to evaluate resistivity of cooling strategies in buildings against overheating impact of climate change. *Building and Environment* **208**, 108599 (2022).

18. Li, D. H. W., Yang, L. & Lam, J. C. Impact of climate change on energy use in the built environment in different climate zones – A review. *Energy* **42**, 103–112 (2012).

19. Centre for Research on the Epidemiology of Disasters - CRED, Ucl, Brussels, Belgium. EM-DAT. https://www.emdat.be/.

20. Machard, A. et al. Climate change influence on buildings dynamic thermal behavior during summer overheating periods: An in-depth sensitivity analysis. *Energy and Buildings* **284**, 112758 (2023).

21. Giorgi, F. et al. The CORDEX-CORE EXP-I Initiative: Description and Highlight Results from the Initial Analysis. *Bulletin of the American Meteorological Society* **103**, E293–E310 (2022).

22. Teichmann, C. et al. Assessing mean climate change signals in the global CORDEX-CORE ensemble. *Clim Dyn* **57**, 1269–1292 (2021).

23. Taylor, K. E., Stouffer, R. J. & Meehl, G. A. An Overview of CMIP5 and the Experiment Design. *Bulletin of the American Meteorological Society* **93**, 485–498 (2012).

24. Eyring, V. et al. Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. *Geosci. Model Dev.* **9**, 1937–1958 (2016).

25. Van Vuuren, D. P. et al. The representative concentration pathways: an overview. *Climatic Change* **109**, 5–31 (2011).
26. Riahi, K. et al. RCP 8.5—A scenario of comparatively high greenhouse gas emissions. *Climatic Change* **109**, 33–57 (2011).

27. Cannon, A. J. & Jeong, D. I. Climate-Resilient Buildings and Core Public Infrastructure 2020: An Assessment of the Impact of Climate Change on Climatic Design Data in Canada. (Environment and Climate Change Canada, Gatineau, QC, CA, 2021).

28. McSweeney, C. F., Jones, R. G., Lee, R. W. & Rowell, D. P. Selecting CMIP5 GCMs for downscaling over multiple regions. *Clim Dyn* **44**, 3237–3260 (2015).

29. World Climate Research Program, C. R. C. D. E. CORDEX CORE Simulations CORDEX Coordinated Output for Regional Evaluations (CORE). https://cordex.org/experiment-guidelines/cortex-core/cortex-core-simulations/.

30. Giorgetta, M. A. et al. Climate and carbon cycle changes from 1850 to 2100 in MPI-ESM simulations for the Coupled Model Intercomparison Project phase 5: Climate Changes in MPI-ESM. *J. Adv. Model. Earth Syst.* **5**, 572–597 (2013).

31. Jacob, D. & Podzun, R. Sensitivity studies with the regional climate model REMO. *Meteorol. Atmos. Phys.* **63**, 119–129 (1997).

32. Jacob, D. A note to the simulation of the annual and inter-annual variability of the water budget over the Baltic Sea drainage basin. *Meteorology and Atmospheric Physics* **77**, 61–73 (2001).

33. Roeckner, E. et al. The atmospheric general circulation model echam-4: model description and simulation of present-day climate. (1996).

34. Hagemann, S. *An Improved Land Surface Parameter Dataset for Global and Regional Climate Models*. (2002).

35. Pfeifer, S. Modeling cold cloud processes with the regional climate model REMO. (Max Planck Institute, Hamburg, 2006).

36. Kotlarski, S. A Subgrid Glacier Parameterisation for Use in Regional Climate Modelling. (Max Planck Institute, Hamburg, 2007).

37. Teichmann, C. Climate and Air Pollution Modelling in South America with Focus on Megacities. (Max Planck Institute, Hamburg, 2010).

38. Semmler, T., Cheng, B., Yang, Y. & Rontu, L. Snow and ice on Bear Lake (Alaska) – sensitivity experiments with two lake ice models. *Tellus A: Dynamic Meteorology and Oceanography* **64**, 17339 (2012).

39. Pietikäinen, J.-P. et al. The regional aerosol-climate model REMO-HAM. *Geosci. Model Dev.* **5**, 1323–1339 (2012).

40. Preuschmann, S. Regional surface albedo characteristics - analysis of albedo data and application to land-cover changes for a regional climate model. (Max Planck Institute, Hamburg, 2012).

41. Wilhelm, C., Rechid, D. & Jacob, D. Interactive coupling of regional atmosphere with biosphere in the new generation regional climate system model REMO-iMOVE. *Geosci. Model Dev.* **7**, 1093–1114 (2014).

42. Maraud, D. Bias Correcting Climate Change Simulations - a Critical Review. *Curr Clim Change Rep* **2**, 211–220 (2016).

43. Faghih, M., Brissette, F., & Sabeti, P. Impact of correcting sub-daily climate model biases for hydrological studies, Hydrol. Earth Syst. Sci., 26, 1545–1563 (2022).

44. Yang, W., Gardelin, M., Olsson, J. & Bosshard, T. Multi-variable bias correction: application of forest fire risk in present and future climate in Sweden. *Nat. Hazards Earth Syst. Sci.* **15**, 2037–2057 (2015).

45. Gaur, A., Lacasse, M. & Armstrong, M. Climate Data to Undertake Hygrothermal and Whole Building Simulations Under Projected Climate Change Influences for 11 Canadian Cities. *Data* **4**, 72 (2019).

46. Cannon, A. J. Multivariate quantile mapping bias correction: an N-dimensional probability density function transform for climate model simulations of multiple variables. *Clim Dyn* **50**, 31–49 (2018).

47. Cannon, A. J., Sobie, S. R. & Murdock, T. Q. Bias Correction of GCM Precipitation by Quantile Mapping: How Well Do Methods Preserve Changes in Quantiles and Extremes? *Journal of Climate* **28**, 6938–6959 (2015).

48. Qingyuan, Z., Huang, J. & Siwei, L. Development of typical year weather data for chinese locations. *ASHRAE Transactions* 2002, v. 108. Pt. 2, held June 2002, Honolulu, Hawaii.

49. Boland, J., Ridley, B. & Brown, B. Models of diffuse solar radiation. *Renewable Energy* **33**, 575–584 (2008).

50. UNI. UNI 10349-1:2016 Heating and Cooling of Buildings-Climatic Data-Part 1: Monthly Means for Evaluation of Energy Need for Space Heating and Cooling and Methods for Splitting Global Solar Irradiance into the Direct and Diffuse Parts and to Calculate the Solar Irradiance on Tilted Planes. (2016).

51. Ridley, B., Boland, J. & Lauret, P. Modelling of diffuse solar fraction with multiple predictors. *Renewable Energy* **35**, 478–483 (2010).

52. Kasten, F. The linke turbidity factor based on improved values of the integral rayleigh optical thickness. *Solar Energy* **56**, 6 (1996).

53. Remund, J., Wald, L., Lefèvre, M., Ranchin, T. & Page, J. Worldwide Linke turbidity information. *ISES Solar World Congress* 14.

54. Remund, J. & Domeisen, D. Aerosol Optical Depth and Linke Turbidity Climatology: Description for Final Report of IEA SHC Task 36. *Meteotest, Bern, Switzerland* (2009).

55. Finkelstein, J. M. & Schafer, R. E. Improved goodness-of-fit tests. *Biometrika* **58**, 6 (1971).

56. U.S. Department of Energy. EnergyPlus™ Version 9.6.0 Documentation/Auxiliary Programs. Chapter 2, 257 (2021).

57. Flores-Larsen, S., Bre, F. & Hongn, M. A performance-based method to detect and characterize heatwaves for building resilience analysis. *Renewable and Sustainable Energy Reviews* **167**, 112795 (2022).

58. Machard, A. et al. IEA EBC Annex 80 'Typical and extreme weather datasets for studying the resilience of buildings to climate change' (Version 1.0). 2024-02-26 https://doi.org/10.26050/WDCC/WDTF_Annex80_build_v1.0 (2024).

59. Meteotest. Meteonorm Software. https://meteonorm.meteotest.ch/en/.

60. Ohunakin, O. S., Adaramola, M. S., Oyewola, O. M., Matthew, O. J. & Fagbenle, R. O. The effect of climate change on solar radiation in Nigeria. *Solar Energy* **116**, 272–286 (2015).

61. Cutforth, H. W. & Judiesch, D. Long-term changes to incoming solar energy on the Canadian Prairie. *Agricultural and Forest Meteorology* **145**, 167–175 (2007).

62. Lee, S. H. & Levinson, R. Cool envelope benefits in future typical weather and heatwave conditions 1 for single-family homes in Los Angeles [accepted paper]. *Sixth International Conference on Countermeasures to Urban Heat Islands, RMIT University, Melbourne, Australia*. https://doi.org/10.20357/B7DK6T (2023).

63. Wilcox, S. & Marion, W. Users Manual for TMY3 Data Sets. Technical Report (NREL/TP-581-43156). National Renewable Energy Lab., (2008).

64. ISO. ISO/TR 17772-2:2018 Energy performance of buildings — Overall energy performance assessment procedures — Part 2: Guideline for using indoor environmental input parameters for the design and assessment of energy performance of buildings. (2018).

65. ISO. ISO 17772-1:2017 Energy performance of buildings — Indoor environmental quality — Part 1: Indoor environmental input parameters for the design and assessment of energy performance of buildings. (2017).

66. Sengupta, A., Breesch, H., Al Assaad, D. & Steeman, M. Evaluation of thermal resilience to overheating for an educational building in future heatwave scenarios. *International Journal of Ventilation* 1–11 https://doi.org/10.1080/14733315.2023.2218424 (2023).

67. Sengupta, A., Al Assaad, D., Bastero, J. B., Steeman, M. & Breesch, H. Impact of heatwaves and system shocks on a nearly zero energy educational building: Is it resilient to overheating? *Building and Environment* **234**, 110152 (2023).

68. IEA EBC Annex 80 - Resilient Cooling of Buildings. Weather Data. https://annex80.iea-ebc.org/weather-data.

69. Flato, G. et al. Evaluation of Climate Models. In: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. (2013).

70. Schwalm, C. R., Glendon, S. & Duffy, P. B. RCP8.5 tracks cumulative CO2 emissions. *Proceedings of the National Academy of Sciences* **117**, 19656–19657 (2020).

71. Daniel, M. et al. Benefits of explicit urban parameterization in regional climate modeling to study climate and city interactions. *Clim Dyn* **52**, 2745–2764 (2019).
72. Langendijk, G. S., Rechid, D., Sieck, K. & Jacob, D. Added value of convection-permitting simulations for understanding future urban humidity extremes: case studies for Berlin and its surroundings. *Weather and Climate Extremes* **33**, 100367 (2021).

73. Trusilova, K. et al. The urban land use in the COSMO-CLM model: a comparison of three parameterizations for Berlin. *metz* **25**, 231–244 (2016).

74. Gao, Y. et al. Implementation of the CORINE land use classification in the regional climate model REMO. **20**, (2020).

75. Taha, H. Site-specific weather files and fine-scale probabilistic microclimate zones for current and future climates and land use. *IBPSA World News* **29–41** (2020).

76. Salvati, A., Palme, M., Chiesa, G. & Kolokotroni, M. Built form, urban climate and building energy modelling: case-studies in Rome and Antofagasta. *Journal of Building Performance Simulation* **13**, 209–225 (2020).

77. Masson, V., Lemonsu, A., Hidalgo, J. & Voogt, J. Urban Climates and Climate Change. *Annu. Rev. Environ. Resour.* **45**, 411–444 (2020).

78. Salvati, A. & Kolokotroni, M. Generating future-urban weather files for building performance simulations: case studies in London. in *Proceedings of the 17th IBPSA Conference*. https://doi.org/10.26868/25222708.2021.30315 (Bruges, Belgium, 2021).

79. Bueno, B., Norford, L., Hidalgo, J. & Pigeon, G. The urban weather generator. *Journal of Building Performance Simulation* **6**, 269–281 (2013).

80. Mao, J., Yang, J. H., Afshari, A. & Norford, L. K. Global sensitivity analysis of an urban microclimate system under uncertainty: Design and case study. *Building and Environment* **124**, 153–170 (2017).

81. Ward, H. C., Kotthaus, S., Järvi, L. & Grimmond, C. S. B. Surface Urban Energy and Water Balance Scheme (SUEWS): Development and evaluation at two UK sites. *Urban Climate* **18**, 1–32 (2016).

82. Leroyer, S., Bélair, S., Spacek, L. & Gultepe, I. Modelling of radiation-based thermal stress indicators for urban numerical weather prediction. *Urban Climate* **25**, 64–81 (2018).

83. Skamarock, W. C. et al. A Description of the Advanced Research WRF Model Version 4. 165.

84. Chen, F. et al. The Integrated WRF/Urban Modeling System: Development, Evaluation, and Applications to Urban Environmental Problems. **40**.

85. Yeh, K.-S. et al. The CMC-MRB Global Environmental Multiscale (GEM) Model. Part III: Nonhydrostatic Formulation. *Mon. Wea. Rev.* **130**, 339–356 (2002).

86. Shu, C. et al. Added value of convection permitting climate modelling in urban overheating assessments. *Building and Environment* **207**, 108415 (2022).

87. Martinez, S. et al. A practical approach to the evaluation of local urban overheating– A coastal city case-study. *Energy and Buildings* **253**, 111522 (2021).

88. Salvati, A., Monti, P., Coch Roura, H. & Cecere, C. Climatic performance of urban textures: Analysis tools for a Mediterranean urban context. *Energy and Buildings* **185**, 162–179 (2019).

89. Afshari, A. & Ramirez, N. Improving the accuracy of simplified urban canopy models for arid regions using site-specific prior information. *Urban Climate* **35**, 100722 (2021).

90. Mortezazadeh, M., Wang, L. L., Albettar, M. & Yang, S. CityFFD – City fast fluid dynamics for urban microclimate simulations on graphics processing units. *Urban Climate* **41**, 101063 (2022).

91. Moradi, M., Krayenhoff, E. S. & Aliabadi, A. A. A comprehensive indoor-outdoor urban climate model with hydrology: The Vertical City Weather Generator (VCWG v2.0.0). *Building and Environment* **207**, 108406 (2022).

92. Gaur, A., Eichenbaum, M. K. & Simonovic, S. P. Analysis and modelling of surface Urban Heat Island in 20 Canadian cities under climate and land-cover change. *Journal of Environmental Management* **206**, 145–157 (2018).

93. Leirvik, T. & Yuan, M. A Machine Learning Technique for Spatial Interpolation of Solar Radiation Observations. *Earth and Space Science* **8**, (2021).

94. Appelhans, T., Mwangomo, E., Hardy, D. R., Hemp, A. & Nauss, T. Evaluating machine learning approaches for the interpolation of monthly air temperature at Mt. Kilimanjaro, Tanzania. *Spatial Statistics* **14**, 91–113 (2015).

95. Machard, A. AMachard/Assembling-future-weather-files-including-heatwaves: v1.0.0 (v1.0.0). Zenodo. https://doi.org/10.5281/zenodo.7300024 (2022).

96. United Nations, Department of Economic and Social Affairs, Population Division. 2022 Revision of World Population Prospects. https://population.un.org/wpp/ (2022).

## Acknowledgements

The authors of the paper wish to thank their colleagues from IEA EBC Annex 80, who provided insight and expertise during discussions on future weather data needed to assess resilient cooling building design. The authors would like to acknowledge the World Climate Research Programme’s Working Group on Regional Climate and the Working Group on Coupled Modelling, the former coordinating body of CORDEX and the responsible panel for CMIP5. We also thank the climate modelling groups (listed in Table 2 of this paper) for producing and making their model output available. We also acknowledge the Earth System Grid Federation infrastructure, an international effort led by the U.S. Department of Energy’s Program for Climate Model Diagnosis and Intercomparison, the European Network for Earth System Modelling, and other partners in the Global Organization for Earth System Science Portals (GO-ESSP). Finally, the authors gratefully thank Alex J. Cannon for providing his code “Multivariate Bias Correction of Climate Model Outputs”. This work was partially supported by the Assistant Secretary for Energy Efficiency and Renewable Energy, Building Technologies Office, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. Dr. Bre performed some of the tasks regarding his contributions to this research during his OE-BUILDINGS action, under grant agreement ID: 101024627. This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement N° 101024627. This research was supported, in part, by Fraunhofer Internal Programs under Grant No. Attract 003-695033. The research was also partially supported by Det Energiteknologisk Udviklingsog Demonstrations Program (EUDP) under grant 64018-0578.

## Author contributions

Anaïs Machard, Agnese Salvati, Mamak P. Tootkaboni, and Abhishek Gaur coordinated the efforts to produce these datasets, designed the methodology, and wrote and edited the draft. All other authors prepared the datasets in the different cities and reviewed and edited the manuscript.

## Competing interests

The authors declare no competing interests.
**Additional information**

**Correspondence** and requests for materials should be addressed to A.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_16.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024