Florent Renard

Lucille Alonso

Magalie Técher
email: magalie.techer@etu.univ-montp3.fr

Hassan Ait Haddou

Rahim Aguejdad

Integrating Urban Planning Policies

Keywords: vulnerability, urban heat islands, urban planning, exposure, sensitivity, overheating, mitigation measures, France

The energy balance differences between urban and rural areas give rise to a representative phenomenon of urban climate change called the urban heat island (UHI). The UHI phenomenon, which is exacerbated by global warming effects and urban expansion, exposes urban areas and populations to thermal stress and overheating. Therefore, measuring the intensity of UHI vulnerability is needed to mitigate its effects on health, comfort, energy consumption and infrastructure. Additionally, urban decision makers are looking for methodologies to identify the most vulnerable areas in urban planning. This research aims to investigate areas prone to urban overheating in Montpellier Mediterranean metropolitan area as a representative case study of the French urban planning context. The objective is to develop a methodology for assessing UHI vulnerability based on the interaction of exposure and sensitivity to urban planning policies. Accordingly, a set of fifteen indicators is developed to visualize urban planning policies, exposure and sensitivity to urban overheating. A statistical-based method is applied to determine the individual weights based on the indicators' variance, which is calculated using principal component analysis (PCA). The spatial and urban planning distribution of vulnerability is then analyzed. The results indicate that the areas with the highest urban planning challenges are most vulnerable to UHI. In fact, the results reveal that these areas have urban planning characteristics likely to aggravate urban overheating. Therefore, this research proposes a methodology for integrating mitigation solutions into urban planning: increase in permeable and green areas; improvement of urban ventilation and heat exchange and decrease in anthropogenic heat generation. It demonstrates that combining urban planning policies with exposure and sensitivity is a new perspective for vulnerability studies. The findings of this study can be utilized by decision makers to better visualize areas susceptible to UHI and plan appropriate mitigation strategies.

Introduction

The concentration of the world's population and human activities within urban environments has an impact on the urban microclimate [START_REF] Lowry | Empirical Estimation of Urban Effects on Climate: A Problem Analysis[END_REF]. Urban areas are subject to urban overheating, characterized by higher temperatures compared to surrounding rural areas [START_REF] Oke | The energetic basis of the urban heat island[END_REF]. This phenomenon, described since the 19th century [START_REF] Howard | The Climate of London, Volume 1: Deduced from Meteorological Observations[END_REF] and known as the urban heat island (UHI), is a consequence of the modification of the environment by human occupation and urban planning [START_REF] Yin | Effects of urban form on the urban heat island effect based on spatial regression model[END_REF][START_REF] Lin | Effects of urban planning indicators on urban heat island: A case study of pocket parks in high-rise high-density environment[END_REF]. Indeed, this phenomenon is reinforced by a set of urban parameters such as urban morphology which disrupts air circulation, leading to the trapping of heat and pollutants [START_REF] Liao | The effect of spatial heterogeneity in urban morphology on surface urban heat islands[END_REF]. They are composed of few vegetated spaces, which limit the cooling power of plants through evapotranspiration [START_REF] Tan | Urban tree design approaches for mitigating daytime urban heat island effects in a high-density urban environment[END_REF]. The materials used in construction and public spaces disrupt the absorption and reflection of solar radiation, resulting in the trapping during the day and the release at night of the accumulated heat [START_REF] Jandaghian | Increasing urban albedo to reduce heat-related mortality in Toronto and Montreal, Canada[END_REF]. Finally, the presence of human activities and the energy used in urban areas are at the origin of the generation of anthropogenic heat in the ambient air [START_REF] Shahmohamadi | The Impact of Anthropogenic Heat on Formation of Urban Heat Island and Energy Consumption Balance[END_REF]. In recent decades, more and more territories are confronted with the UHI phenomenon. According to Santamouris [START_REF] Santamouris | Recent progress on urban overheating and heat island research. Integrated assessment of the energy, environmental, vulnerability and health impact. Synergies with the global climate change[END_REF][START_REF] Santamouris | Innovating to Zero the Building Sector in Europe: Minimising the Energy Consumption, Eradication of the Energy Poverty and Mitigating the Local Climate Change[END_REF], this phenomenon is well documented in more than 400 cities in Europe and around the world, and the intensity of overheating is currently close to 10 • C. This increase in urban temperatures is not without consequences. Indeed, combined with global warming and urban densification, UHI's effects are particularly visible in cities which are confronted by heat waves, leading to health risks and discomfort for the most fragile inhabitants and over-consumption of energy linked to cooling [START_REF] Buchin | Evaluation of the health-risk reduction potential of countermeasures to urban heat islands[END_REF][START_REF] Fouillet | Excess Mortality Related to the August 2003 Heat Wave in France[END_REF].

Within the scientific community, these concerns have led to improved knowledge in urban climatology of the UHI phenomenon [START_REF] Ren | The Application of Urban Climatic Mapping to the Urban Planning of High-Density Cities: The Case of Kaohsiung, Taiwan[END_REF]. Intensive research has, therefore, been carried out to develop and test measures to combat the effects of UHI at the micro and macro level [START_REF] Rosenfeld | Mitigation of urban heat island: Material, utility, programs, updates[END_REF][START_REF] Tsoka | Urban warming and cities' microclimates: Investigation methods and mitigation strategies-A review[END_REF][START_REF] Szali Ńska | Local Urban Risk Assessment of Dry and Hot Hazards for Planning Mitigation Measures[END_REF]. One of the major advances in this field has highlighted solutions linked to urban forms and adaptation, such as increasing green spaces and the installation of reflective or green materials and evaporative systems. Hulley in 2012 [START_REF] Hulley | The Urban Heat Island Effect: Causes and Potential Solutions[END_REF], Akbari in 2016 [START_REF] Akbari | Three decades of urban heat island and mitigation technologies research[END_REF] and Santamouris in 2017 [START_REF] Santamouris | Passive and Active Cooling for the Outdoor Built Environment-Analysis and Assessment of the Cooling Potential of Mitigation Technologies Using Performance Data from 220 Large Scale Projects[END_REF] carried out state-of-the-art studies of mitigation solutions and their cooling potential. More recently, Bernard and al. [START_REF] Bernard | Rafraîchissement des Villes: Solutions Existantes et Pistes de Recherche[END_REF] proposed a classification of solutions into three main categories, namely: green solutions which are based on the systemic services of nature (i.e., water and vegetation) which enables to cool urban temperatures; gray solutions based on engineering systems using water, materials, urban form or the regulation of energy consumption and soft solutions based on political strategies and uses.

Interest in the UHI phenomenon has also enabled the development of research aimed at identifying and locating areas most prone to overheating. Numerous studies have been conducted to determine the vulnerability to the effects of UHI of several international cities and metropolitan areas [START_REF] Gonzalez-Trevizo | Research Trends on Environmental, Energy and Vulnerability Impacts of Urban Heat Islands: An Overview[END_REF]. Vulnerability is described by Turner as the susceptibility to experience and resist harm due to exposure to hazard, disturbance or stress [START_REF] Turner | A framework for vulnerability analysis in sustainability science[END_REF]. Thus, studies of vulnerability to the effects of UHI combine the scientific knowledge acquired in several fields such as urban climatology and social sciences, which impact and exacerbate the vulnerability of a territory [START_REF] Cheng | Approaches for Identifying Heat-Vulnerable Populations and Locations: A Systematic Review[END_REF]. On the one hand, by considering exposure to urban overheating, these studies report parameters that intensify vulnerability [START_REF] Apreda | A Climate Vulnerability and Impact Assessment Model for Complex Urban Systems[END_REF]. On the other hand, sensitivity and adaptive capacity refer to the population and reflect their fragility and degree of resistance to overheating [START_REF] Alonso | A Comparative Study of the Physiological and Socio-Economic Vulnerabilities to Heat Waves of the Population of the Metropolis of Lyon (France) in a Climate Change Context[END_REF][START_REF] Sabrin | Quantifying Environmental and Social Vulnerability: Role of Urban Heat Island and Air Quality, a Case Study of Camden, NJ[END_REF]. These studies use various indicators, which can be classified into two groups. The first group includes the socioeconomic characteristics that allow to identify populations based on the conditions that heighten their risk exposure such as age, gender, health conditions and financial resources [START_REF] Cheng | Approaches for Identifying Heat-Vulnerable Populations and Locations: A Systematic Review[END_REF][START_REF] Alonso | A Comparative Study of the Physiological and Socio-Economic Vulnerabilities to Heat Waves of the Population of the Metropolis of Lyon (France) in a Climate Change Context[END_REF][START_REF] Sabrin | Quantifying Environmental and Social Vulnerability: Role of Urban Heat Island and Air Quality, a Case Study of Camden, NJ[END_REF]. The second group concerns the environmental characteristics which include climatic, urban, architectural and landscape conditions of risk exposure [START_REF] Cheng | Approaches for Identifying Heat-Vulnerable Populations and Locations: A Systematic Review[END_REF][START_REF] Apreda | A Climate Vulnerability and Impact Assessment Model for Complex Urban Systems[END_REF].

The urban microclimate and its impacts on populations are at the heart of reflections for urban decision makers and constitute a challenge for the sustainable development of their territory. However, until recently, within urban planning, the consequences of town planning and urbanization on the urban climate are less integrated [START_REF] Mills | Progress toward sustainable settlements: A role for urban climatology[END_REF]. Indeed, the prioritization of these questions within the territories is still recent, as a result, the integration of knowledge acquired in urban climatology into employable methodologies and solutions is complex for urban decision makers [START_REF] Ng | The Urban Climatic Map for Sustainable Urban Planning[END_REF]. However, there are advances in trying to better integrate this knowledge into urban planning [START_REF] Hurlimann | Urban Planning Policy Must Do More to Integrate Climate Change Adaptation and Mitigation Actions[END_REF][START_REF] Bitan | The methodology of applied climatology in planning and building[END_REF][START_REF]Climate Change 2014: Impacts, Adaptation and Vulnerability. Part B: Regional Aspects; Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[END_REF]. This concerns the integration, within studies of vulnerability to UHI effects, of parameters making it possible to identify the effects of urban form on UHIs to create a link between urban planning and exposure to urban overheating. Weber et al. [START_REF] Weber | Policy-Relevant Indicators for Mapping the Vulnerability of Urban Populations to Extreme Heat Events: A Case Study of Philadelphia[END_REF] carried out an analysis of Philadelphia to propose a grid of reproducible indicators in the United States that can be integrated into urban planning. In France, urban planners are wondering the scale of their territory about strategies and mitigation measures to adopt to reduce urban overheating vulnerability. This is the case in the Ile-de-France region [START_REF]Adapter l'Ile de France à la Chaleur Urbaine: Identifier les Zones à Effet D'îlot de Chaleur Urbain (ICU) et Etablir Leur Degré de Vulnérabilité Afin de Mieux Anticiper[END_REF], which carried out a UHI vulnerability study based on Locale Climate Zone classification [START_REF] Stewart | Local Climate Zones for Urban Temperature Studies[END_REF] and used cross-referencing with socioeconomic data. The results of this work enabled stakeholders to define adaptation solutions and to integrate new tools into the urban planning practice. The French Environment Agency, ADEME, has also released a report of French territorial applications for urban overheating diagnostics [START_REF] Ademe | Diagnostic de la Surchauffe Urbaine-Méthodes et Applications Territoriales[END_REF]. This report enables local authorities to integrate this issue by proposing replicable methodologies according to the needs of the territory. These different studies and methodologies are promising and should be considered as a major asset for urban planning to help reduce UHI effects and determine vulnerable areas, in the context of global warming and urban densification. There is a need to integrate this issue into the urban planning process and urban projects [START_REF] Mills | Progress toward sustainable settlements: A role for urban climatology[END_REF][START_REF] Ng | The Urban Climatic Map for Sustainable Urban Planning[END_REF][START_REF] Hurlimann | Urban Planning Policy Must Do More to Integrate Climate Change Adaptation and Mitigation Actions[END_REF][START_REF] Bitan | The methodology of applied climatology in planning and building[END_REF][START_REF]Climate Change 2014: Impacts, Adaptation and Vulnerability. Part B: Regional Aspects; Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[END_REF]. In the long term, approaches for urban climate remain limited and have little impact, while there is a need to transfer knowledge from urban climatology into urban planning.

Based on several vulnerability studies of the effects of UHIs, this research attempts to develop a methodology for integrating urban climatology and urban planning for better consideration of UHI within urban environments in a French context. Our hypothesis is that the integration of urban planning rules into the calculation of UHI effects will provide decision makers with a vulnerability index that can be interpreted and utilized in the mitigation of UHI and maintaining the thermal comfort of urban dwellers. Actually, the use of the indicators derived from the regulations will provide planners with data which are adapted to their practices and on which they can act to reduce the effects of the UHI.

Accordingly, the objective of this work is to carry out a spatial mapping of UHI vulnerability at the territory level using urban planning policies, exposure to UHIs and sensitivity of populations to urban overheating. The combination of these three components permits us to identify the most vulnerable sectors of the territory which make the link between urban planning and knowledge of urban climatology.

To operate the articulation between climatology science and urban planning policies, the article proceeds in two stages. On the one hand, the components and their indicators are described, as well as their known behaviors on UHI vulnerability from the scientific literature. On the other hand, the calculation of the vulnerability is carried out from the superposition of all the indicators and components, leading to the generation of a global map of day and night vulnerability. We will be dealing with demonstrating that, from the location of vulnerable areas and knowledge of their characteristics, urban decision makers are owning the means to mitigate the effects of UHI within urban planning documents. Specifically, the analysis of the urban planning component will let them define appropriate adaptation solutions that can be integrated into urban planning plans and documents. This paper is organized into five sections. The first section analyzes the results of the existing research studies that begin to make a link between UHI vulnerability and urban planning policies. The second section presents the study area and materials. In this section, we expose our method development; whereas, the impact of each indicator on UHI is discussed according to the scientific literature review. The third section highlights the main results including an investigation of the regulatory characteristics of the different zones, land uses and urban planning documents. These results are then discussed in the fourth section which, also, outlines the different regulatory solutions that are proposed in order to mitigate urban heat island effects. The main conclusions and future work are addressed in the last section.

Materials and Methods

The study methodology consists of two stages. First, fifteen indicators are integrated as input data. Urban planning policies and the impacts of exposure and sensitivity to overheating are all considered, and their effects compiled to generate a UHI vulnerability map of the territory (Figure 1). The second step is to analyze and highlight vulnerable areas characteristics. Finally, the scientific literature, rich in solutions for mitigating UHI, enables us to develop recommendations for French urban planning.

characteristics. Finally, the scientific literature, rich in solutions for mitigating UHI, enables us to develop recommendations for French urban planning.

Study Area

The territory of Montpellier Méditerranée metropolitan area is in the southeast of France; it is bordered by the Mediterranean Sea (Figure 2). With a land area of approximately 49,000 ha, it contains 31 municipalities which hosted 481,276 inhabitants in 2018 [START_REF]Populations Légales 2018-Recensement de la Population-Regions, Départements, Arrondissements, Cantons et Communes[END_REF]. According to the Köppen-Geiger classification [START_REF] Peel | Updated World Map of the Köppen-Geiger Climate Classification[END_REF], the climate is typically Mediterranean with hot, dry summers, with a mean annual temperature of 20.9 °C in 2015 and soft winters with a mean annual temperature of 11.3 °C in 2015 [START_REF] Audric | Montpellier Méditerranée Métropole: Diagnostic de Territoire[END_REF]. This territory has experienced urban sprawl in recent years caused by its attractiveness, with an annual population growth of 1.85% from 2010 to 2015, i.e., nearly 8000 new inhabitants per year [START_REF] Weber | Policy-Relevant Indicators for Mapping the Vulnerability of Urban Populations to Extreme Heat Events: A Case Study of Philadelphia[END_REF]. As a result, UHI affects urban areas of the territory, causing discomfort and health problems. Indeed, according to Météo France, in Figure 3 the annual mean temperature of the territory and number of heat waves day have been increasing for several decades [START_REF]Plan Climat Air Énergie Territorial-Diagnostic Énergie-Climat Territorial[END_REF]. This work focuses on an urban vulnerability assessment of the metropolitan area by using various urban indicators to elaborate urban planning measures.

Study Area

The territory of Montpellier Méditerranée metropolitan area is in the southeast of France; it is bordered by the Mediterranean Sea (Figure 2). With a land area of approximately 49,000 ha, it contains 31 municipalities which hosted 481,276 inhabitants in 2018 [START_REF]Populations Légales 2018-Recensement de la Population-Regions, Départements, Arrondissements, Cantons et Communes[END_REF]. According to the Köppen-Geiger classification [START_REF] Peel | Updated World Map of the Köppen-Geiger Climate Classification[END_REF], the climate is typically Mediterranean with hot, dry summers, with a mean annual temperature of 20.9 • C in 2015 and soft winters with a mean annual temperature of 11.3 • C in 2015 [START_REF] Audric | Montpellier Méditerranée Métropole: Diagnostic de Territoire[END_REF].

Study Area

The territory of Montpellier Méditerranée metropolitan area is in the southeast of France; it is bordered by the Mediterranean Sea (Figure 2). With a land area of approximately 49,000 ha, it contains 31 municipalities which hosted 481,276 inhabitants in 2018 [START_REF]Populations Légales 2018-Recensement de la Population-Regions, Départements, Arrondissements, Cantons et Communes[END_REF]. According to the Köppen-Geiger classification [START_REF] Peel | Updated World Map of the Köppen-Geiger Climate Classification[END_REF], the climate is typically Mediterranean with hot, dry summers, with a mean annual temperature of 20.9 °C in 2015 and soft winters with a mean annual temperature of 11.3 °C in 2015 [START_REF] Audric | Montpellier Méditerranée Métropole: Diagnostic de Territoire[END_REF]. This territory has experienced urban sprawl in recent years caused by its attractiveness, with an annual population growth of 1.85% from 2010 to 2015, i.e., nearly 8000 new inhabitants per year [START_REF] Weber | Policy-Relevant Indicators for Mapping the Vulnerability of Urban Populations to Extreme Heat Events: A Case Study of Philadelphia[END_REF]. As a result, UHI affects urban areas of the territory, causing discomfort and health problems. Indeed, according to Météo France, in Figure 3 the annual mean temperature of the territory and number of heat waves day have been increasing for several decades [START_REF]Plan Climat Air Énergie Territorial-Diagnostic Énergie-Climat Territorial[END_REF]. This work focuses on an urban vulnerability assessment of the metropolitan area by using various urban indicators to elaborate urban planning measures. This territory has experienced urban sprawl in recent years caused by its attractiveness, with an annual population growth of 1.85% from 2010 to 2015, i.e., nearly 8000 new inhabitants per year [START_REF] Weber | Policy-Relevant Indicators for Mapping the Vulnerability of Urban Populations to Extreme Heat Events: A Case Study of Philadelphia[END_REF]. As a result, UHI affects urban areas of the territory, causing discomfort and health problems. Indeed, according to Météo France, in Figure 3 the annual mean temperature of the territory and number of heat waves day have been increasing for several decades [START_REF]Plan Climat Air Énergie Territorial-Diagnostic Énergie-Climat Territorial[END_REF]. This work focuses on an urban vulnerability assessment of the metropolitan area by using various urban indicators to elaborate urban planning measures.

Data Acquisition

A geographic information system (GIS) is used as a data processing and generation tool (QGIS software, version 3.22). Its use in this work may allow later use of output data produced [START_REF] Mineau | L'apport des SIG en urbanisme (The contribution of GIS to town planning)[END_REF]. The data are acquired from the following climatic, urban, environmental, demographic and socio-economic databases. Table 1 summarizes the data sources and spatial resolution. Data relating to urban planning are taken from Local Urban Plans, which are urban planning documents operating at the municipal or inter-municipal scale. Urban planning regulations of each municipality in the Montpellier metropolitan area have been integrated into the GIS. Land surface temperatures were calculated from satellite observation data and remote sensing. Daytime data are derived from the LANDSAT-8 satellite image (22 July 2019) and the daytime land surface temperature is calculated using algorithms developed by Avdan and Jovanovska [START_REF] Avdan | Algorithm for Automated Mapping of Land Surface Temperature using LANDSAT 8 Satellite Data[END_REF]. Nighttime data are acquired by the Advanced Space-Borne Thermal Emission and Reflection Radiometer (ASTER) satellite image (6 July 2015) and the nighttime land surface temperature is calculated thanks to the single-channel algorithm presented by Jimenez-Munoz and Sobrino [START_REF] Jimenez-Munoz | A Single-Channel Algorithm for Land-Surface Temperature Retrieval from ASTER Data[END_REF]. Demographic and socio-economic data of the metropolis of Montpellier come from the National Institute of Statistics and Economic Studies (INSEE). These are census data collected in 2015 [START_REF]Documentation-Données Carroyées FILOSOFI[END_REF]. The data provided are represented by a mesh that can have a spatial resolution varying from 200 m to 1 km. Data relating to the categorization of land uses come from the open data of the Montpellier Mediterranean metropolis, which provides the public with updated territory data, and from aerial view interpretation (2019-2020). Characterization of building age was extracted from land files produced by CEREMA in 2020 [START_REF]Cerema Données Foncières Retraitées par le Cerema-DGALN DGFIP[END_REF].

Data Acquisition

A geographic information system (GIS) is used as a data processing and generation tool (QGIS software, version 3.22). Its use in this work may allow later use of output data produced [START_REF] Mineau | L'apport des SIG en urbanisme (The contribution of GIS to town planning)[END_REF]. The data are acquired from the following climatic, urban, environmental, demographic and socio-economic databases. Table 1 summarizes the data sources and spatial resolution. Data relating to urban planning are taken from Local Urban Plans, which are urban planning documents operating at the municipal or inter-municipal scale. Urban planning regulations of each municipality in the Montpellier metropolitan area have been integrated into the GIS. Land surface temperatures were calculated from satellite observation data and remote sensing. Daytime data are derived from the LANDSAT-8 satellite image (22 July 2019) and the daytime land surface temperature is calculated using algorithms developed by Avdan and Jovanovska [START_REF] Avdan | Algorithm for Automated Mapping of Land Surface Temperature using LANDSAT 8 Satellite Data[END_REF]. Nighttime data are acquired by the Advanced Space-Borne Thermal Emission and Reflection Radiometer (ASTER) satellite image (6 July 2015) and the nighttime land surface temperature is calculated thanks to the single-channel algorithm presented by Jimenez-Munoz and Sobrino [START_REF] Jimenez-Munoz | A Single-Channel Algorithm for Land-Surface Temperature Retrieval from ASTER Data[END_REF]. Demographic and socio-economic data of the metropolis of Montpellier come from the National Institute of Statistics and Economic Studies (INSEE). These are census data collected in 2015 [START_REF]Documentation-Données Carroyées FILOSOFI[END_REF]. The data provided are represented by a mesh that can have a spatial resolution varying from 200 m to 1 km. Data relating to the categorization of land uses come from the open data of the Montpellier Mediterranean metropolis, which provides the public with updated territory data, and from aerial view interpretation (2019-2020). Characterization of building age was extracted from land files produced by CEREMA in 2020 [START_REF]Cerema Données Foncières Retraitées par le Cerema-DGALN DGFIP[END_REF].

Identification of Areas Subject to Urban Heat Island Effect

From data obtained previously which constitute the input layers of each component, indicators of urban planning policies, exposure and sensitivity to the risks of urban overheating are developed. Collected data are unified, cut and divided in pixels of 200 m × 200 m size. Table 1 presents the processes that allow us to assign a value to each cell of the study area based on the various data collected upstream. The socio-economic indicators, namely population density, sensitive population and low-income households were available on INSEE database at 200 m spatial resolution. There was no adaptation of these data. For indicators with higher spatial resolution, a cell-based clipping was performed and the value (class) of the spatial entity representing the largest area was assigned to the entire cell. A class was assigned to a cell as soon as this class was predominant over the others. Land surface temperatures were calculated by averaging the land surface temperature of 30 m spatial resolutions of each entity contained in each cell. Finally, the urban regulation was assigned from the superposition of the zone where each rule applies with each cell. A rule value is assigned to the cell if it overlaps a zone with the same rule value.

This identification provides users with a simple but detailed cutting of a similar size to the district and block scale. Additionally, this mesh allows better representation of the results and subsequent analysis by urban decision makers. Due to the large number of indicators that enter the calculation of urban vulnerability, two methods are becoming widespread for assigning weights in a relevant way: explicit weights and statistical methods [START_REF] Cheng | Approaches for Identifying Heat-Vulnerable Populations and Locations: A Systematic Review[END_REF]. For this study, the method of explicit weights was used, assigning an equal weight to each of the components and indicators involved in the calculation of urban vulnerability. For this study, a statistical and sensitive analysis method for multivariate outputs-was used based on principal component analysis (PCA). A PCA-based approach is commonly applied in vulnerability studies because they allow correlated variables to be converted into uncorrelated variables and reduce the number of dimensions to be considered. The method implemented by Wolf et al. [START_REF] Wolf | Integrated Assessment of Vulnerability to Heat Stress in Urban Areas[END_REF] for the development of a heat wave vulnerability index for London (UK) was applied [START_REF] Wolf | The development of a heat wave vulnerability index for London, United Kingdom[END_REF]. It was also implemented by Alonso [START_REF] Alonso | Intérêt de la Modélisation de la Température de l'air Associé à la Nécessité de la Caractérisation des Vulnérabilités Territoriales pour une[END_REF] in the calculation of socio-economic vulnerability to high heat in Lyon (France) and Tokyo (Japan).

A principal component analysis with Varimax rotation was performed. Thus, PCA allows to retain the three vulnerability's components: urban planning policies, urban overheating exposure and urban overheating sensitivity, which are retained because their eigenvalues are greater than 1. The components are then weighted by the value of the variance calculated by the PCA. Each of the indicators included in the components is weighted by the value of correlation coefficient obtained by the Varimax rotation. This method was chosen because all the indicators do not have the same behavior on the vulnerability's component. Finally, each component is integrated into the following equation (Equation ( 1)) to calculate the UHI vulnerability (V):

V = 0.26 P + 0.39 E + 0.33 S (1)
where P corresponds to urban planning policies; E is the urban overheating exposure and S is the urban overheating sensitivity.

Component 1: Urban Planning Policies

Local Urban Plans are the urban planning documents which define territory development [START_REF]PLUI-Un Outil pour L'avenir des Territoires Guide[END_REF]. These documents are of capital importance because they determine urban morphology and aesthetics. By orientations taken during their elaboration, urban planning documents have an impact on socio-economic development but also on climatic behavior of territories. Local Urban Plans are made up of an urban planning regulation, which defines for each plot of the territory architectural, urban and landscape rules [START_REF] Ecriture | [END_REF]. Plots are cadastral land located in the same municipality, subject to a set of urban regulation rules. Regulation defines the form and appearance of each construction on the plots. Consisting of 16 articles (Table 2), it is a document that can be opposed to any construction project or extension. As part of this study, five articles of urban planning regulation are considered to define the challenges represented by urban planning policies (P) according to the following equation (Equation ( 2)):

P = 0.43CES + 0.64Vg + 0.54H + 0.20A + 0.20I ( 2
)
where CES is the building footprint coefficient (%); Vg is the proportion of green and planted spaces (%); H is the building height (m); A corresponds to the installation of air conditioning systems and I represents the construction's location to parcel boundaries. These indicators (Table 3) are known throughout the scientific literature to have impacts on the UHI intensification (Figure 4). On the one hand, the increase in planted spaces and vegetation within the urban environment has a positive impact because it leads to a decrease in urban temperatures. It is negatively correlated with the surface urban heat island (SUHI) [START_REF] Wolf | The development of a heat wave vulnerability index for London, United Kingdom[END_REF]. Indeed, a 10% increase in green spaces leads to a reduction of 1 • C in land surface temperature [START_REF] Alonso | Intérêt de la Modélisation de la Température de l'air Associé à la Nécessité de la Caractérisation des Vulnérabilités Territoriales pour une[END_REF] and a tree cover of at least 16% achieves a 1 • C drop in air temperature [START_REF]PLUI-Un Outil pour L'avenir des Territoires Guide[END_REF]. According to the literature, the lower the percentage of green/planted space, the higher the risk of overheating [START_REF] Suzán-Azpiri | The role of urban vegetation in temperature and heat island effects in Querétaro city, Mexico[END_REF][START_REF] Li | The Use of Reflective and Permeable Pavements as a Potential Practice for Heat Island Mitigation and Stormwater Management[END_REF][START_REF] Marando | Urban Heat Island Mitigation by Green Infrastructure in European Functional Urban Areas[END_REF][START_REF] Abreu-Harbich | Effect of Tree Planting Design and Tree Species on Human Thermal Comfort in the Tropics[END_REF][START_REF] Dimoudi | Vegetation in the Urban Environment: Microclimatic Analysis and Benefits[END_REF]. The proportion of planted space is described in Article 13 of the Local Urbanism Plan and is determined by a coefficient ranging from 0 to 1 for each plot, from the lowest proportion to the highest. On the other hand, the other indicators retained seem to have a negative impact on the intensification of the urban heat island. Thus, a study revealed that densely built environments, characterized by their high proportion of impermeable spaces, undergo greater urban overheating [START_REF] Escourrou | Climat et la Ville[END_REF][START_REF] Arnfield | Two Decades of Urban Climate Research: A Review of Turbulence, Exchanges of Energy and Water, and the Urban Heat Island[END_REF]. Indeed, it also concludes that the increase in built On the other hand, the other indicators retained seem to have a negative impact on the intensification of the urban heat island. Thus, a study revealed that densely built environments, characterized by their high proportion of impermeable spaces, undergo greater urban overheating [START_REF] Escourrou | Climat et la Ville[END_REF][START_REF] Arnfield | Two Decades of Urban Climate Research: A Review of Turbulence, Exchanges of Energy and Water, and the Urban Heat Island[END_REF]. Indeed, it also concludes that the increase in built density and urban height lead to an intensification of UHI [START_REF] Li | On the Influence of Density and Morphology on the Urban Heat Island Intensity[END_REF][START_REF] Goncalves | The Environmental Performance of Tall Buildings[END_REF]. Within Local Urban Plans, the built area is described in Article 9 while the height of the constructions appears in Article 10 of the Local Urban Plans. Determined by a coefficient ranging from 0 to 1, the built area and the maximum height of constructions are indicators of urban density and disturb sunshine, shade and heat exchange. Similarly, the location of buildings on the plot also has a negative impact on urban heat islands by altering air circulation [START_REF] Nakamura | Temperature and Stability Conditions in an East-West Oriented Urban Canyon[END_REF]. Indeed, the implantation can lead to a change in the speed, turbulence and direction of the winds, and, therefore, lead to the trapping of heat and increasing the UHI's effects [START_REF] Athamena | Proposal of indicators dedicated to the analysis of contribution of urban projects to urban Heat Island[END_REF]. This indicator is defined in Articles 6 and 7 of the Local Urban Plans, which determines, concerning the location of the construction, either an alignment or a withdrawal from the parcel boundaries. Finally, the last indicator considered in this work accounts for the contribution of anthropogenic heat to the intensification of UHIs. An example of this is the installation of air conditioners on the buildings, resulting in cooling the interior and a generation of heat which heats the outside air and the unequipped surrounding buildings [START_REF] De Munck | Modélisation de la Végétation Urbaine et des Stratégies D'adaptation au Changement Climatique pour L'amélioration du Confort Climatique et de la Demande Energétique en Ville[END_REF]. A study has shown that the urban temperature is 1 • C warmer near buildings equipped with air conditioners, thus increasing the effects of the UHI [START_REF] Salamanca | Anthropogenic Heating of the Urban Environment Due to Air Conditioning[END_REF]. For this work, all the areas whose urban regulations, within Article 11, authorize the installation of air conditioners on the main facades are more likely to be vulnerable to urban overheating than those which are prohibited from installing them.

Component 2: Exposure to Urban Overheating

This component considers the exposure of the territory (E) to urban overheating. It is interested in characterizing the extent of the UHI phenomenon and the probabilities of aggravation of the phenomenon linked to the particularities of the territory. For this, more specifically, this work determines the SUHI, impacted by its topography, proximity to natural areas and the presence of population and activities according to the following equation (Equation (3)):

E = 0.22T + 0.79D + 0.83U + 0.73∆T + 0.90Pg + 0.90Pw (3)
where T is the topography (m); D is the population density (number of people per cell); U corresponds to land uses; ∆T is the SUHI intensity ( • C); Pg is the proximity to green areas and Pw, the proximity to water areas (m). Among the indicators selected to characterize exposure to UHI effects (Table 4), topography and proximity to natural areas show a positive impact by reducing UHI intensity (Figure 5). On the one hand, when the altitude increases, the air temperature decreases under the effects of changes in atmospheric pressure. It is estimated that the average air temperature decreases by 0.6 • C every 100 m of elevation [START_REF] Ren | The Application of Urban Climatic Mapping to the Urban Planning of High-Density Cities: The Case of Kaohsiung, Taiwan[END_REF]. The topography of the metropolitan area of Montpellier is obtained from the DTM (Digital Terrain Model) and is relatively flat, except for certain areas at the limits of the territory. It, therefore, appears that the impact of the topography on the temperature will be minimal on the majority of the study area, given that the 0-150 m class is largely in the majority. On the other hand, it has been shown that proximity to natural spaces allows an aircooling effect, thus reducing UHI intensity. It appears that green spaces allow, by providing shade, the reduction in the heating of materials [START_REF] Suzán-Azpiri | The role of urban vegetation in temperature and heat island effects in Querétaro city, Mexico[END_REF], and by evapotranspiration, the cooling of the air, especially at night [START_REF] Doick | The Role of One Large Greenspace in Mitigating London's Nocturnal Urban Heat Island[END_REF]. Similarly, water spaces, especially the larger ones, can absorb heat, cool the atmosphere and, thus, reduce overheating [START_REF] Revaud | Thermoregulatory effect of green spaces and wetlands in Paris[END_REF]. Like On the other hand, it has been shown that proximity to natural spaces allows an aircooling effect, thus reducing UHI intensity. It appears that green spaces allow, by providing shade, the reduction in the heating of materials [START_REF] Suzán-Azpiri | The role of urban vegetation in temperature and heat island effects in Querétaro city, Mexico[END_REF], and by evapotranspiration, the cooling of the air, especially at night [START_REF] Doick | The Role of One Large Greenspace in Mitigating London's Nocturnal Urban Heat Island[END_REF]. Similarly, water spaces, especially the larger ones, can absorb heat, cool the atmosphere and, thus, reduce overheating [START_REF] Revaud | Thermoregulatory effect of green spaces and wetlands in Paris[END_REF]. Like Yiannakou and Salata [START_REF] Yiannakou | Adaptation to climate change through spatial planning in compact urban areas: A case study in the city of Thessaloniki[END_REF] concerning Park Cool Island's effect, which defines an area within which the cooling effect of a natural space is perceived, the indicators selected represent the proximity of urban spaces to natural spaces. This is to highlight the proximity or distance of built cells to green or blue cells and the possibility of benefiting from their cooling effects. For green spaces, it is considered that beyond 300 m distance for spaces of more than 20 ha, the cooling effect no longer has an impact. Concerning the aquatic spaces, the metropolis is crossed by two major rivers, and is bordered by the Mediterranean Sea and by ponds, so beyond 200 m of distance with these spaces, the effects of freshness are no longer perceived. Cells far from these natural spaces are, therefore, likely to be less cooled and, furthermore, to suffer more from urban overheating.

The UHI phenomenon presents a temperature difference between the urban environment and the rural environment. It is a function of three types of phenomena: surface UHI, the canopy layer and the boundary layer [START_REF] Bahi | Review of Methods for Retrieving Urban Heat Islands[END_REF][START_REF] Veena | Urban heat island studies: Current status in india and comparison with the international studies[END_REF]. Thus, the intensity of the surface urban heat island (∆T) is a good indicator for studying the UHI phenomenon in a territory (Figure 6). Calculated from obtaining the ground surface temperatures (LST) for day and night [START_REF] Mineau | L'apport des SIG en urbanisme (The contribution of GIS to town planning)[END_REF][START_REF] Avdan | Algorithm for Automated Mapping of Land Surface Temperature using LANDSAT 8 Satellite Data[END_REF]. For this work, it is the thermal gradient which is the difference at time t between average temperature of natural areas N listed in the Local Urban Plans and average temperatures recorded for each cell. Natural areas (N) considered characterize restricted or prohibited urbanization to protect the green areas of the territory. These are natural areas, as notified in urban planning documents, which circle urban areas. It excluded urban green spots whose cells are classified as urban or to be urbanized. The following equation (Equation ( 4)) was used to calculate the thermal gradient:

∆T (cell n) = T (cell n) -T (N) (4)
where T (cell n) is the average land surface temperature of cell n ( • C) and T (N) is the average land surface temperature of natural areas (N) listed in urban planning regulations, which serves as a reference ( Yiannakou and Salata [START_REF] Yiannakou | Adaptation to climate change through spatial planning in compact urban areas: A case study in the city of Thessaloniki[END_REF] concerning Park Cool Island's effect, which defines an area within which the cooling effect of a natural space is perceived, the indicators selected represent the proximity of urban spaces to natural spaces. This is to highlight the proximity or distance of built cells to green or blue cells and the possibility of benefiting from their cooling effects. For green spaces, it is considered that beyond 300 m distance for spaces of more than 20 ha, the cooling effect no longer has an impact. Concerning the aquatic spaces, the metropolis is crossed by two major rivers, and is bordered by the Mediterranean Sea and by ponds, so beyond 200 m of distance with these spaces, the effects of freshness are no longer perceived. Cells far from these natural spaces are, therefore, likely to be less cooled and, furthermore, to suffer more from urban overheating.

The UHI phenomenon presents a temperature difference between the urban environment and the rural environment. It is a function of three types of phenomena: surface UHI, the canopy layer and the boundary layer [START_REF] Bahi | Review of Methods for Retrieving Urban Heat Islands[END_REF][START_REF] Veena | Urban heat island studies: Current status in india and comparison with the international studies[END_REF]. Thus, the intensity of the surface urban heat island (ΔT) is a good indicator for studying the UHI phenomenon in a territory (Figure 6). Calculated from obtaining the ground surface temperatures (LST) for day and night [START_REF] Mineau | L'apport des SIG en urbanisme (The contribution of GIS to town planning)[END_REF][START_REF] Avdan | Algorithm for Automated Mapping of Land Surface Temperature using LANDSAT 8 Satellite Data[END_REF]. For this work, it is the thermal gradient which is the difference at time t between average temperature of natural areas N listed in the Local Urban Plans and average temperatures recorded for each cell. Natural areas (N) considered characterize restricted or prohibited urbanization to protect the green areas of the territory. These are natural areas, as notified in urban planning documents, which circle urban areas. It excluded urban green spots whose cells are classified as urban or to be urbanized. The following equation (Equation ( 4)) was used to calculate the thermal gradient:

ΔT (cell n) = T (cell n) -T(N) ( 4
)
where T (cell n) is the average land surface temperature of cell n (°C) and T (N) is the average land surface temperature of natural areas (N) listed in urban planning regulations, which serves as a reference (°C). One of the main causes of the UHI formation is the modification of the urban environment by human presence. Therefore, population density and land use are suitable indicators for studying UHIs (Figure 6). On the one hand, the more densely populated an area is, the higher the generation of anthropogenic heat is, which increases population exposure to urban overheating [START_REF] Cheng | Approaches for Identifying Heat-Vulnerable Populations and Locations: A Systematic Review[END_REF][START_REF] Apreda | A Climate Vulnerability and Impact Assessment Model for Complex Urban Systems[END_REF][START_REF] Alonso | A Comparative Study of the Physiological and Socio-Economic Vulnerabilities to Heat Waves of the Population of the Metropolis of Lyon (France) in a Climate Change Context[END_REF][START_REF] Sabrin | Quantifying Environmental and Social Vulnerability: Role of Urban Heat Island and Air Quality, a Case Study of Camden, NJ[END_REF][START_REF] Mills | Progress toward sustainable settlements: A role for urban climatology[END_REF][START_REF] Ng | The Urban Climatic Map for Sustainable Urban Planning[END_REF][START_REF]Adapter l'Ile de France à la Chaleur Urbaine: Identifier les Zones à Effet D'îlot de Chaleur Urbain (ICU) et Etablir Leur Degré de Vulnérabilité Afin de Mieux Anticiper[END_REF]. On the other hand, the uses of the ground influence the UHI, because of forms, density and energy consumption characteristics which are distinct and dependent on the moment of the day [START_REF] Ren | The Application of Urban Climatic Mapping to the Urban Planning of High-Density Cities: The Case of Kaohsiung, Taiwan[END_REF][START_REF] Zhang | Quantifying the Relationship between Land Use Features and Intra-Surface Urban Heat Island Effect: Study on Downtown Shanghai[END_REF][START_REF] Giridharan | Daytime urban heat island effect in high-rise and high-density residential developments in Hong Kong[END_REF]. During the day, occupied uses such as employment areas will be more exposed because they generate heat, while at night, the most exposed uses will be more residential.

Component 3: Sensitivity to Urban Overheating

This component, which is systematically integrated into urban vulnerability studies, is particularly interested in humans and the constructions present on the territory. It allows us to analyze the fragility of populations, their habitat and their resources in coping with heat waves and urban overheating. It results from an analysis of indicators of sensitivity to the effects of climate change cited in the scientific literature [START_REF] Cheng | Approaches for Identifying Heat-Vulnerable Populations and Locations: A Systematic Review[END_REF][START_REF] Apreda | A Climate Vulnerability and Impact Assessment Model for Complex Urban Systems[END_REF][START_REF] Alonso | A Comparative Study of the Physiological and Socio-Economic Vulnerabilities to Heat Waves of the Population of the Metropolis of Lyon (France) in a Climate Change Context[END_REF][START_REF] Sabrin | Quantifying Environmental and Social Vulnerability: Role of Urban Heat Island and Air Quality, a Case Study of Camden, NJ[END_REF][START_REF] Weber | Policy-Relevant Indicators for Mapping the Vulnerability of Urban Populations to Extreme Heat Events: A Case Study of Philadelphia[END_REF][START_REF]Adapter l'Ile de France à la Chaleur Urbaine: Identifier les Zones à Effet D'îlot de Chaleur Urbain (ICU) et Etablir Leur Degré de Vulnérabilité Afin de Mieux Anticiper[END_REF] and the data available in the case of the metropolis of Montpellier. The sensitivity component (S) is calculated from the following equation (Equation ( 5)):

S = 0.66PA + 0.16R + 0.40HC + 0.73CA ( 5
)
where PA is the proportion of sensitive populations by age (%); R is the proportion of lowincome households (%); HC is the presence of health and care facilities and CA corresponds to average building age (year).

The sensitivity indicators (Table 5 and Figure 7) used in this study have an impact on the increase in urban vulnerability to urban heat islands. Indeed, the proportion of fragile population, the presence of retirement and health establishments and household income are indicators of the presence of a population vulnerable due to their state of health, their age or their resources. On the one hand, the age of the population reflects the fragile persons that require more assistance and/or often find themselves isolated during heat waves [START_REF] Fouillet | Excess Mortality Related to the August 2003 Heat Wave in France[END_REF][START_REF] Li | Heat Waves and Morbidity: Current Knowledge and Further Direction-A Comprehensive Literature Review[END_REF]. It appears that the youngest and oldest populations show greater difficulty in regulating their metabolism to evacuate their body heat and, therefore, suffer more from the effects of urban overheating. The location of retirement homes and health establishments reflect the concentration of frail and elderly populations at one point. On the other hand, the economic income of the populations is considered because low economic resources impact their vulnerability by reducing their possibilities of adapting as well as possible to face urban overheating. Indeed, many studies have shown a difficulty in developing housing (ventilation, insulation) for the most economically fragile populations [START_REF] Alonso | A Comparative Study of the Physiological and Socio-Economic Vulnerabilities to Heat Waves of the Population of the Metropolis of Lyon (France) in a Climate Change Context[END_REF][START_REF] Alizadeh | Increasing Heat-Stress Inequality in a Warming Climate[END_REF]. On the one hand, the age of the population reflects the fragile persons that require more assistance and/or often find themselves isolated during heat waves [START_REF] Fouillet | Excess Mortality Related to the August 2003 Heat Wave in France[END_REF][START_REF] Li | Heat Waves and Morbidity: Current Knowledge and Further Direction-A Comprehensive Literature Review[END_REF]. It appears that the youngest and oldest populations show greater difficulty in regulating their metabolism to evacuate their body heat and, therefore, suffer more from the effects of urban overheating. The location of retirement homes and health establishments reflect the concentration of frail and elderly populations at one point. On the other hand, the economic income of the populations is considered because low economic resources impact their vulnerability by reducing their possibilities of adapting as well as possible to face urban overheating. Indeed, many studies have shown a difficulty in developing housing (ventilation, insulation) for the most economically fragile populations [START_REF] Alonso | A Comparative Study of the Physiological and Socio-Economic Vulnerabilities to Heat Waves of the Population of the Metropolis of Lyon (France) in a Climate Change Context[END_REF][START_REF] Alizadeh | Increasing Heat-Stress Inequality in a Warming Climate[END_REF].

Finally, the age of buildings (Figure 7) is also included in the study because it can also exacerbate urban vulnerability to urban heat island effects. Indeed, the construction period of the buildings is an indicator of their potential discomfort. Indeed, in France, studies have revealed that buildings built before the first thermal regulations of 1974 are likely to have insulation problems making them thermal sieves [START_REF] Cheng | Approaches for Identifying Heat-Vulnerable Populations and Locations: A Systematic Review[END_REF][START_REF]Adapter l'Ile de France à la Chaleur Urbaine: Identifier les Zones à Effet D'îlot de Chaleur Urbain (ICU) et Etablir Leur Degré de Vulnérabilité Afin de Mieux Anticiper[END_REF]. This lack of insulation can lead to difficulties in adapting populations and to energy overconsumption linked to heating in winter and air conditioning in summer.

Results

This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn. Finally, the age of buildings (Figure 7) is also included in the study because it can also exacerbate urban vulnerability to urban heat island effects. Indeed, the construction period of the buildings is an indicator of their potential discomfort. Indeed, in France, studies have revealed that buildings built before the first thermal regulations of 1974 are likely to have insulation problems making them thermal sieves [START_REF] Cheng | Approaches for Identifying Heat-Vulnerable Populations and Locations: A Systematic Review[END_REF][START_REF]Adapter l'Ile de France à la Chaleur Urbaine: Identifier les Zones à Effet D'îlot de Chaleur Urbain (ICU) et Etablir Leur Degré de Vulnérabilité Afin de Mieux Anticiper[END_REF]. This lack of insulation can lead to difficulties in adapting populations and to energy overconsumption linked to heating in winter and air conditioning in summer.

Vulnerability Assessment to UHI Effects

Results

This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.

Vulnerability Assessment to UHI Effects

The vulnerability to UHI effects is analyzed with regards to the selected drivers set. These fifteen input parameters are related to urban planning policies, exposure and sensitivity to urban overheating. The final map of urban vulnerability to UHI effects is generated at the Montpellier metropolitan area scale based on a geographic cross-check of the different scores assigned to each cell. The daytime, nighttime and global vulnerability maps reflect the climatic issues and challenges facing the Montpellier metropolitan area (Figure 8). The analysis of the spatial variations of the vulnerability intensity during the day, at night and globally is carried out according to the number of cells, the associated surface area and the most affected areas of the territory.

The vulnerability to UHI effects is analyzed with regards to the selected drivers set. These fifteen input parameters are related to urban planning policies, exposure and sensitivity to urban overheating. The final map of urban vulnerability to UHI effects is generated at the Montpellier metropolitan area scale based on a geographic cross-check of the different scores assigned to each cell. The daytime, nighttime and global vulnerability maps reflect the climatic issues and challenges facing the Montpellier metropolitan area (Figure 8). The analysis of the spatial variations of the vulnerability intensity during the day, at night and globally is carried out according to the number of cells, the associated surface area and the most affected areas of the territory. During the day, the low-vulnerability class is dominant for half of the study area. In fact, it accounts for 57.9% of the total number of cells and 56.6% of the total cell surface During the day, the low-vulnerability class is dominant for half of the study area. In fact, it accounts for 57.9% of the total number of cells and 56.6% of the total cell surface area. The results also show that 42.1% of cells present a medium level of vulnerability and account for 43.4% of the cell surface area. Actually, this proportion corresponds to more than 4800 cells and 19,000 ha of the territory. Conversely, the proportion of high level of vulnerability does not exceed 12% of cells and 12.4% of surfaces, which corresponds to 1367 cells and 5456 ha of the study area compared to the high nighttime vulnerability surfaces which account for 11.9%. Such a difference is explained by the surface UHI which is derived based on land surface temperature that generates higher thermal gradients during the day. In addition to that, the zones of high daytime vulnerability are found to be more located around urban centers.

At night, the low-vulnerability class is also dominant. It represents 60.5% of the number of cells and 59.2% of the surface area. Actually, the number of cells affected by low vulnerability at night increased by 2.6% compared to the daytime ones. Therefore, the proportion of cells associated with medium-to-high vulnerability levels decreased compared to the daytime score. Both medium-to-high vulnerability levels together account for 39.5% of the number of cells and 40.8% of the surface area, which corresponds to 4500 cells and nearly 18,000 ha.

Globally, the same observation can be made. Very low-to-low vulnerability levels are dominant. They represent 59.1% of the cells and 57.8% of the total surface area. The surface affected by medium-to-very high vulnerability corresponds to 48.8% of the cells and 42.2% of the total surface area. Those levels concern 4700 cells and 18,500 ha.

The spatial patterns of the vulnerability intensity exhibit significant variations among sectors and between the city center and the surrounding suburban areas. The results show that 66.2% of the core of the metropolitan area is exposed to medium and high levels of vulnerability, with 67.5% during the day and 64.8% at night. Actually, 10,500 ha of the total surface area are affected by urban overheating during the day. In contrast, the first and second suburban sectors exhibit a lower level of vulnerability, probably due to lower density of buildings, higher vegetation coverage and abundant water bodies. In fact, 34.2% and 25.7% of these areas are, respectively, concerned by medium-to-high levels of vulnerability globally. Furthermore, Montpellier municipality is the most affected sector followed by "Vallée du Lez". In fact, 87.5% of Montpellier municipality is vulnerable and 44.6% is exposed to a high level of vulnerability. The results indicate that 63% of the "Vallée du Lez" sector is exposed to medium and high levels of vulnerability, where 17% is affected by high-level intensity. The sectors "Cadoule et Bérange" and "Plaine Ouest" are the least affected areas with, respectively, 27.7% and 28.8% of their surface classified with mediumto-very high vulnerability levels. Particularly, their surface affected by high-to-very high vulnerability concerned, respectively, 5.8% and 4.7%.

Vulnerability Assessment Concerning "Urban Planning Policies" Component

Vulnerability Assessment of Global Component

The vulnerability intensity is assessed according to the level of urban planning challenges by quantifying the area of the associated surface. The results reveal that areas, which present the highest level of urban planning challenges, represent one third of the total cells affected by the highest vulnerability score. Actually, these areas account for 35.2% and 33.7% of the total surface affected by the highest daytime and nighttime vulnerability intensity, respectively. These areas, which amount to more than 3800 cells, cover 15,000 ha. Furthermore, areas with medium-to-high planning challenges characterized by a high level of vulnerability account for 12.1% of the total surface during the day, 11.6% at night and 11.8% globally. These affected surfaces correspond, respectively, to 5308 ha, 5109 ha and 5168 ha. Additionally, areas with medium-to-high planning challenges characterized by a medium level of vulnerability represent 23.2% of the surface during the day, 22% at night and 22.8 globally. Their total area reaches more than 10,000 ha. Moreover, areas with medium-to-high planning challenges represent 19.9% and 21.5% of the total zones concerned by low daytime and nighttime vulnerability, respectively. These areas cover more than 8700 ha during the day and 9400 ha at night.

Vulnerability Assessment of "Urban Planning Policies" Indicators

The results of the vulnerability assessment based on urban planning regulation indicators are presented in Figure 9. Cells with a medium-to-high building footprint coefficient refer to a footprint greater than or equal to 0.35 (classes 2 and 3). Actually, their area accounts for 9.18% of the metropolitan area of Montpellier. The results show that 96.6% of these affected cells show a medium-high vulnerability score, which corresponds to 4000 ha. In addition to that, 36.9% of these medium-to-high building footprint cells are found to be highly vulnerable during the day compared to 34.8% at night.

The results of the vulnerability assessment based on urban planning regulation indicators are presented in Figure 9. Cells with a medium-to-high building footprint coefficient refer to a footprint greater than or equal to 0.35 (classes 2 and 3). Actually, their area accounts for 9.18% of the metropolitan area of Montpellier. The results show that 96.6% of these affected cells show a medium-high vulnerability score, which corresponds to 4000 ha. In addition to that, 36.9% of these medium-to-high building footprint cells are found to be highly vulnerable during the day compared to 34.8% at night. Moreover, cells with low-to-medium percentage of planted space refer to a percentage less than or equal to 65% (classes 2 and 3). Actually, these areas account for 24.39% of the total surface of the territory. The results indicate that 88.1%, 85.9% and 87.60% of the affected cells have a medium-to-high vulnerability score during the day, at night and globally, respectively. In fact, 37.9% of these cells show a high score of vulnerability for the daytime compared to 37% for the nighttime.

Cells with a construction height of more than 10 m (classes 3, 4 and 5) represent 13.84% of the study area. However, cells associated with buildings less than 10 m in height (classes 1 and 2) account for 57.46% of the total area. In addition to that, cells concerned with high heights are more confronted with a medium-to-high level of vulnerability than cells concerned by low heights. Actually, these cells refer to 82.9% and 47.5% of the total affected cells, respectively. Such a difference is more pronounced when the analysis is Moreover, cells with low-to-medium percentage of planted space refer to a percentage less than or equal to 65% (classes 2 and 3). Actually, these areas account for 24.39% of the total surface of the territory. The results indicate that 88.1%, 85.9% and 87.60% of the affected cells have a medium-to-high vulnerability score during the day, at night and globally, respectively. In fact, 37.9% of these cells show a high score of vulnerability for the daytime compared to 37% for the nighttime.

Cells with a construction height of more than 10 m (classes 3, 4 and 5) represent 13.84% of the study area. However, cells associated with buildings less than 10 m in height (classes 1 and 2) account for 57.46% of the total area. In addition to that, cells concerned with high heights are more confronted with a medium-to-high level of vulnerability than cells concerned by low heights. Actually, these cells refer to 82.9% and 47.5% of the total affected cells, respectively. Such a difference is more pronounced when the analysis is based on the high level of vulnerability. In fact, the results indicate that cells with more than 10 m in height account for an average of 36% of the total affected cells compared to cells less than 10 m which represent 12.2%. Moreover, cells concerned by the minimum alignment restriction on one of the parcel boundaries (classes 3 and 4) constitute 35.99% of the total metropolitan area. An average of 51.4% of this proportion is affected by a medium-to-high level of vulnerability. Indeed, 16.5% of these cells face a high level of vulnerability, which represents 2200 ha. Additionally, cells concerned by the regulation which authorizes the installation of air conditioners on the main facade (class 2) account for 22.69% of the territory. The results show that 39.6% of these cells are affected by medium-high vulnerability, whereas 10.3% are exposed to high vulnerability.

Vulnerability Analysis of Local Urban Plans

Vulnerability Assessment of Urban Planning Zoning

The UHI vulnerability is analyzed according to the zoning of Local Urban Plans (Figure 10). Local Urban Plans are composed, as well as urban regulation, of urban zoning on which the regulation is applied. Zoning plans are cartographic documents representing a sectoral division of each territory [START_REF] Ecriture | [END_REF]. They are composed of four types of zones which have specific land uses: urbanized zones (U); areas to be urbanized (AU); natural areas (N) and agricultural areas (A). Actually, U and AU zones account for 32% of the entire territory compared to A and N zones, which represent 67.3%. The results exhibit that 93.6% (daytime), 90.9% (nighttime) and 92.6% (global) of these zones, which refer to areas where urbanization is likely to occur, are affected by medium-to-high vulnerability. Their total area accounts for more than 3500 cells and 13,000 ha. However, 22.5%, 19.3% and 21% of medium-to-high vulnerability cells during daytime, nighttime and globally, respectively, are located in agricultural areas. Moreover, only 14.6% (daytime), 12.4% (nighttime) and 13.5% (global) of natural areas are affected by medium-to-high vulnerability levels. The total amount of these affected cells is 5200 ha during the day, 4400 ha at night and 4800 ha globally.

Vulnerability Assessment of Authorized Land Use

The results of the vulnerability assessment according to land use type are shown in Figure 11. Actually, authorized land use within each zone is regulated by the Local Urban Plan through the two articles, 1 and 2 [START_REF] Ecriture | [END_REF]. These articles define the type of permitted use: economic activities, housing, mixed zones, services and equipment, agricultural, natural, without construction or prohibited. The analysis of their vulnerability highlights that the

Vulnerability Assessment of Authorized Land Use

The results of the vulnerability assessment according to land use type are shown in Figure 11. Actually, authorized land use within each zone is regulated by the Local Urban Plan through the two articles, 1 and 2 [START_REF] Ecriture | [END_REF]. These articles define the type of permitted use: economic activities, housing, mixed zones, services and equipment, agricultural, natural, without construction or prohibited. The analysis of their vulnerability highlights that the areas that can accommodate housing (housing + mixed areas) amount approximately to 24% of the territory's cells. The results indicate that 96.2% (daytime), 95.2% (nighttime) and 95.9% (global) of these areas are affected by medium-to-high vulnerability. Their total area accounts for more than 10,000 ha of the territory. Areas that can accommodate activities and employment represent 20.3% of the territory. Indeed, the results show that 93.9% (daytime), 90.3% (nighttime) and 92.4% (global) of these areas are affected by medium-to-high vulnerability, which corresponds to more than 8500 ha. In addition to that, Figure 11 shows that slight variation in medium and high vulnerability between daytime and nighttime is observed. and 95.9% (global) of these areas are affected by medium-to-high vulnerability. Their total area accounts for more than 10,000 ha of the territory. Areas that can accommodate activities and employment represent 20.3% of the territory. Indeed, the results show that 93.9% (daytime), 90.3% (nighttime) and 92.4% (global) of these areas are affected by medium-to-high vulnerability, which corresponds to more than 8500 ha. In addition to that, Figure 11 shows that slight variation in medium and high vulnerability between daytime and nighttime is observed.

Vulnerability Assessment according to the Date of Approval of Local Urban Plans

The Local Urban Plan (PLU) must be approved before it can be adopted. In fact, it has to be subject to a local public enquiry after the preparation and consultation processes. The publication date refers to the date from which the Local Urban Plan document is applied. Once adopted, the PLU has a legal and binding force. The Local Urban Plans of Montpellier metropolitan area were approved from 2013 to 2020 (Figure 12). In fact, the majority of these 31 urban planning documents were approved in 2018 and 2019, with a rate of 6 and 7 documents per year, respectively. However, the approved documents in 2013 and 2019 covered the largest cumulative surface area which accounted for 41% of the

Vulnerability Assessment According to the Date of Approval of Local Urban Plans

The Local Urban Plan (PLU) must be approved before it can be adopted. In fact, it has to be subject to a local public enquiry after the preparation and consultation processes. The publication date refers to the date from which the Local Urban Plan document is applied. Once adopted, the PLU has a legal and binding force. The Local Urban Plans of Montpellier metropolitan area were approved from 2013 to 2020 (Figure 12). In fact, the majority of these 31 urban planning documents were approved in 2018 and 2019, with a rate of 6 and 7 documents per year, respectively. However, the approved documents in 2013 and 2019 covered the largest cumulative surface area which accounted for 41% of the total metropolitan area. The correlation between the vulnerability gradient and the different approval dates of the 31 urban planning documents is carried out. The results reveal that the territories, whose planning documents were approved in 2013, 2019 and 2020, are composed of approximately half of the cells affected by medium-to-high vulnerability with 48.7% (daytime), 46% (nighttime) and 47.4% (global).

Moreover, medium-to-high vulnerability, which is associated with the approved urban planning documents in 2017, accounted for 58.8% of the total affected cells during the day, 55.5% at night and 57.1% globally. In 2020, it reached more than 83.6% and 81.1%, respectively during the day and the night.

Discussion

Mitigation and Adaptation Measures

The results of the vulnerability analysis enable the localization and the quantification of the areas that are most vulnerable to UHI effects. Additionally, this research adopted an innovative approach which permits the integration of urban planning indicators in the analysis of climatic behavior. These indicators are derived from the regulation rules that are specifically used in the urban planning of Montpellier Metropolitan. Therefore, this study is the first known attempt of integrating urban planning in mapping and analyzing the vulnerability intensity and patterns at the Montpellier Metropolitan scale.

A set of 15 indicators is derived from regulatory, climatic, environmental and socioeconomic databases. These input drivers are organized into urban planning policy, The correlation between the vulnerability gradient and the different approval dates of the 31 urban planning documents is carried out. The results reveal that the territories, whose planning documents were approved in 2013, 2019 and 2020, are composed of approximately half of the cells affected by medium-to-high vulnerability with 48.7% (daytime), 46% (nighttime) and 47.4% (global).

Moreover, medium-to-high vulnerability, which is associated with the approved urban planning documents in 2017, accounted for 58.8% of the total affected cells during the day, 55.5% at night and 57.1% globally. In 2020, it reached more than 83.6% and 81.1%, respectively during the day and the night.

Discussion

Mitigation and Adaptation Measures

The results of the vulnerability analysis enable the localization and the quantification of the areas that are most vulnerable to UHI effects. Additionally, this research adopted an innovative approach which permits the integration of urban planning indicators in the analysis of climatic behavior. These indicators are derived from the regulation rules that are specifically used in the urban planning of Montpellier Metropolitan. Therefore, this study is the first known attempt of integrating urban planning in mapping and analyzing the vulnerability intensity and patterns at the Montpellier Metropolitan scale.

A set of 15 indicators is derived from regulatory, climatic, environmental and socioeconomic databases. These input drivers are organized into urban planning policy, exposure and sensitivity to urban overheating components. These three components are combined in order to generate a map of urban vulnerability at the territory scale. The results indicate that areas most vulnerable to UHI have various characteristics which impact urban overheating in several ways. Actually, these affected areas present high urban planning challenges which consist of the planning regulations that impose few planted areas, most construction, a medium-to-large height (over 10 m), minimal alignment and planning permission for installing air conditioning on the main facade. These results are consistent with those of other published studies that deal with the impact of urban form parameters on UHI [START_REF] Yin | Effects of urban form on the urban heat island effect based on spatial regression model[END_REF][START_REF] Lin | Effects of urban planning indicators on urban heat island: A case study of pocket parks in high-rise high-density environment[END_REF][START_REF] Liao | The effect of spatial heterogeneity in urban morphology on surface urban heat islands[END_REF][START_REF] Tan | Urban tree design approaches for mitigating daytime urban heat island effects in a high-density urban environment[END_REF][START_REF] Marando | Urban Heat Island Mitigation by Green Infrastructure in European Functional Urban Areas[END_REF][START_REF] Abreu-Harbich | Effect of Tree Planting Design and Tree Species on Human Thermal Comfort in the Tropics[END_REF][START_REF] Dimoudi | Vegetation in the Urban Environment: Microclimatic Analysis and Benefits[END_REF][START_REF] Li | On the Influence of Density and Morphology on the Urban Heat Island Intensity[END_REF][START_REF] Goncalves | The Environmental Performance of Tall Buildings[END_REF][START_REF] De Munck | Modélisation de la Végétation Urbaine et des Stratégies D'adaptation au Changement Climatique pour L'amélioration du Confort Climatique et de la Demande Energétique en Ville[END_REF][START_REF] Salamanca | Anthropogenic Heating of the Urban Environment Due to Air Conditioning[END_REF][START_REF] Yiannakou | Adaptation to climate change through spatial planning in compact urban areas: A case study in the city of Thessaloniki[END_REF][START_REF] Agathangelidis | Integrating Urban Form, Function, and Energy Fluxes in a Heat Exposure Indicator in View of Intra-Urban Heat Island Assessment and Climate Change Adaptation[END_REF][START_REF] Erdem | An Analysis of Urban Form Factors Driving Urban Heat Island: The Case of Izmir[END_REF]. They permit us to define measures and to assist planners in drawing up urban planning documents. Therefore, a series of mitigation and adaptation solutions is suggested based on the findings of the current study and the literature which has grown over the past few years (Table 6). The literature review and the current study reveal that increasing permeable and natural spaces in vulnerable areas can mitigate the UHI phenomenon [START_REF] Rosenfeld | Mitigation of urban heat island: Material, utility, programs, updates[END_REF][START_REF] Tsoka | Urban warming and cities' microclimates: Investigation methods and mitigation strategies-A review[END_REF][START_REF] Szali Ńska | Local Urban Risk Assessment of Dry and Hot Hazards for Planning Mitigation Measures[END_REF][START_REF] Hulley | The Urban Heat Island Effect: Causes and Potential Solutions[END_REF][START_REF] Akbari | Three decades of urban heat island and mitigation technologies research[END_REF][START_REF] Santamouris | Passive and Active Cooling for the Outdoor Built Environment-Analysis and Assessment of the Cooling Potential of Mitigation Technologies Using Performance Data from 220 Large Scale Projects[END_REF][START_REF] Bernard | Rafraîchissement des Villes: Solutions Existantes et Pistes de Recherche[END_REF][START_REF] Yiannakou | Adaptation to climate change through spatial planning in compact urban areas: A case study in the city of Thessaloniki[END_REF]. In fact, trees and green areas refresh urban air and improve its quality. In French urban planning, increasing the percentage of minimal planted spaces within urban planning regulations could bring more freshness in the summer heat. Additionally, green roofs and green facades are also beneficial for improving the thermal comfort during the day [START_REF] De Munck | Modélisation de la Végétation Urbaine et des Stratégies D'adaptation au Changement Climatique pour L'amélioration du Confort Climatique et de la Demande Energétique en Ville[END_REF]. Accordingly, urban planners could focus on more authorization of this material within urban planning regulations. For instance, the ALUR law of 2014 encourages urban decision makers to set a biotope coefficient in urban planning regulation. The biotope coefficient consists of introducing or maintaining eco-friendly and non-waterproofed surfaces in the plot (floor, wall, roof, etc.) [START_REF] Ecriture | [END_REF].

Moreover, the current study suggests that further research should study relationships between urban form, the coefficient of building footprint and the objectives of decision makers in terms of urban morphology and density. The implicit solution, which consists of going up in height, cannot be applied everywhere. Therefore, an optimal solution, which permits to restrict the building footprint while respecting the density and morphology objectives, must be found. This solution consists of developing a compromise between all the density parameters, including the height and the building footprint

Solutions to Improve Urban Ventilation and Heat Exchanges

The urban form characterizes all the elements that compose the urban environment, such as roads, plots, buildings and vacant spaces. The urban form has an impact on urban overheating because it affects the circulation of the winds, the trapping of heat and the creation of shading [START_REF] Agathangelidis | Integrating Urban Form, Function, and Energy Fluxes in a Heat Exposure Indicator in View of Intra-Urban Heat Island Assessment and Climate Change Adaptation[END_REF][START_REF] Erdem | An Analysis of Urban Form Factors Driving Urban Heat Island: The Case of Izmir[END_REF]. Open islands of controlled density (height, prospect, etc.) allow better ventilation and promote urban cooling [START_REF] Nakamura | Temperature and Stability Conditions in an East-West Oriented Urban Canyon[END_REF][START_REF] Athamena | Proposal of indicators dedicated to the analysis of contribution of urban projects to urban Heat Island[END_REF]. Urban planners should focus on the regulation of the buildings' height and the road width in order to optimize the wind circulation and the contribution of solar radiation. Moreover, it is known that increasing the albedo of materials is a cooling solution [START_REF] Rosenfeld | Mitigation of urban heat island: Material, utility, programs, updates[END_REF][START_REF] Tsoka | Urban warming and cities' microclimates: Investigation methods and mitigation strategies-A review[END_REF][START_REF] Szali Ńska | Local Urban Risk Assessment of Dry and Hot Hazards for Planning Mitigation Measures[END_REF][START_REF] Hulley | The Urban Heat Island Effect: Causes and Potential Solutions[END_REF][START_REF] Akbari | Three decades of urban heat island and mitigation technologies research[END_REF][START_REF] Santamouris | Passive and Active Cooling for the Outdoor Built Environment-Analysis and Assessment of the Cooling Potential of Mitigation Technologies Using Performance Data from 220 Large Scale Projects[END_REF][START_REF] Bernard | Rafraîchissement des Villes: Solutions Existantes et Pistes de Recherche[END_REF]. In fact, it permits to reduce the absorption of solar radiation from surfaces and, therefore, helps to fight against their heating [START_REF] Abreu-Harbich | Effect of Tree Planting Design and Tree Species on Human Thermal Comfort in the Tropics[END_REF]. However, this solution must be handled with care because of the risk of increasing energy consumption during winter in temperate climates or creating dazzling for pedestrians. Consequently, urban decision makers should impose materials with characteristics that cool the outdoor temperature without creating discomfort.

Solutions to Reduce Anthropogenic Heat Generation

This last solution is based on strategies that refer to the uses present in urban environment [START_REF] Akbari | Three decades of urban heat island and mitigation technologies research[END_REF][START_REF] Santamouris | Passive and Active Cooling for the Outdoor Built Environment-Analysis and Assessment of the Cooling Potential of Mitigation Technologies Using Performance Data from 220 Large Scale Projects[END_REF][START_REF] Bernard | Rafraîchissement des Villes: Solutions Existantes et Pistes de Recherche[END_REF]. Actually, human activities are responsible for the heat generation, particularly through the use of air conditioning which contributes to the increase in outside temperatures [START_REF] Salamanca | Anthropogenic Heating of the Urban Environment Due to Air Conditioning[END_REF]. Accordingly, one of the solutions that help reduce energy consumption consists of encouraging bioclimatic buildings by supporting construction projects. In addition to that, inhabitants should have the ability to create conditions for adapting to heat waves actively. The challenge in urban planning is to continue this effort by restricting the use of air conditioning in buildings.

Research Highlights and Limitations

This research reveals that areas mostly vulnerable are located within the sectors which are associated with the most urban planning challenges, and the most factors of exposure and sensitivity to urban overheating. These results are in accordance with several case studies carried out in France [START_REF] Alonso | A Comparative Study of the Physiological and Socio-Economic Vulnerabilities to Heat Waves of the Population of the Metropolis of Lyon (France) in a Climate Change Context[END_REF][START_REF]Adapter l'Ile de France à la Chaleur Urbaine: Identifier les Zones à Effet D'îlot de Chaleur Urbain (ICU) et Etablir Leur Degré de Vulnérabilité Afin de Mieux Anticiper[END_REF][START_REF] Ademe | Diagnostic de la Surchauffe Urbaine-Méthodes et Applications Territoriales[END_REF] and in other countries [START_REF] Cheng | Approaches for Identifying Heat-Vulnerable Populations and Locations: A Systematic Review[END_REF][START_REF] Apreda | A Climate Vulnerability and Impact Assessment Model for Complex Urban Systems[END_REF][START_REF] Sabrin | Quantifying Environmental and Social Vulnerability: Role of Urban Heat Island and Air Quality, a Case Study of Camden, NJ[END_REF][START_REF] Weber | Policy-Relevant Indicators for Mapping the Vulnerability of Urban Populations to Extreme Heat Events: A Case Study of Philadelphia[END_REF][START_REF] Sabrin | Developing Vulnerability Index to Quantify Urban Heat Islands Effects Coupled with Air Pollution: A Case Study of Camden, NJ[END_REF][START_REF] Leal Filho | Coping with the Impacts of Urban Heat Islands. A Literature Based Study on Understanding Urban Heat Vulnerability and the Need for Resilience in Cities in a Global Climate 927 Change Context[END_REF][START_REF] Bhattacharjee | Assessment of Different Methodologies for Mapping Urban Heat Vulnerability for Milan, Italy[END_REF]. In fact, previous studies have demonstrated that the higher the concentration of factors, the higher the vulnerability and the damage caused by hazards [START_REF] Yang | natural hazards and risk: Spatial assessment of urban heatwave vulnerability of coupling adaptability based on BPNN model: A case study of Xiamen City, China. Geomat[END_REF].

This research demonstrates that integrating urban planning regulation in assessing vulnerability to overheating is promising to better control UHI impacts and climate change adaptation. The results will allow drawing the first solutions to fight against the main causes of a UHI: decrease in permeable and natural spaces, disruption of ventilation and radiative exchanges and overheating due to human activity. It may provide considerable support to policy makers, designers and city planners to address the urban heat issue for maintaining urban thermal comfort, especially in the French urban context. Actually, it provides them with practical solutions and measures to identify vulnerable areas, to prepare strategies and where to concentrate mitigation and adaptation actions. In addition to that, taking into consideration urban planning policies through the integration of regulatory indicators enables them to target the morphological and architectural characteristics of each specific urban zone. Thus, specific and adapted regulatory solutions can be proposed locally to limit the formation and exacerbation of UHI impacts in future urban planning documents.

The use of a GIS tool leads urban planners to perform geospatial data analysis, mapping and visualization of vulnerability and different indicators. Nevertheless, urban planners are still looking for suitable tools that can permit them to integrate knowledge from urban climatology in urban planning.

However, a number of limitations needs to be considered. First, a limitation of this study is that the number of the selected indicators was limited. Particularly, only five indicators are derived from the local urban planning policy. Secondly, the current investigation is limited by data availability. Actually, certain factors are not examined due to the lack of information, even though they have been already addressed in previous studies [START_REF] Weber | Policy-Relevant Indicators for Mapping the Vulnerability of Urban Populations to Extreme Heat Events: A Case Study of Philadelphia[END_REF][START_REF]Adapter l'Ile de France à la Chaleur Urbaine: Identifier les Zones à Effet D'îlot de Chaleur Urbain (ICU) et Etablir Leur Degré de Vulnérabilité Afin de Mieux Anticiper[END_REF][START_REF] Sandholz | Rethinking Urban Heat Stress: Assessing Risk and Adaptation Options across Socioeconomic Groups in Bonn, Germany[END_REF]. Therefore, the results suggest that future investigations on the UHI effect should consider the climatic and meteorological conditions. For example, further climatic data collection related to air temperatures or humidity parameters is required to exactly determine the exposure to UHI intensity. Thirdly, this work is limited by the fact that the calculation of the UHI vulnerability developed in his study uses a statistical approach by a principal component analysis to weight by the variance the impact of the various components. This approach has been applied many times in urban overheating studies and allows us to make the results more reliable. However, it can affect the importance of certain variables, which increase the UHI vulnerability as the SUHI intensity is characterized by the thermal gradient. Consequently, a search will have to be undertaken in order to answer the questions concerned with the weight attributed to all indicators. Future developments should compare vulnerability maps from several calculation methods (equal explicit weighting, weighting by several groups of experts, etc.). Finally, the most important limitation lies in the division of the study area. The cells that represent the basic unit of this work measure 200 m × 200 m, and by their size do not allow the integration of phenomena occurring at the microclimatic scale in the vulnerability analysis. Indeed, phenomena such as cool areas offered by natural spaces of small sizes in urban sectors are not taken into account precisely. Therefore, future developments of this work could improve the spatial resolution by the use of a finer division of the territory and the integration of more precise indicators.

Despite the limitations described above, this research has several positive aspects. In fact, this work is one of the few studies that is interested in characterizing UHI vulnerability and its connection with urban planning regulation in France. This research, which will serve as a base for future studies, has several practical applications. The methods used for the vulnerability analysis may be applied to other cities elsewhere in France, as well as in the world. In fact, the application of the current research methodology to other French territories, which are controlled by Local Urban Plan, is possible thanks to the integration of data from national public databases (IGN, INSEE, etc.) and urban planning documents. Indeed, one of the key strengths of this study is its use of remote sensing and geographic information systems which are commonly used in urban planning [START_REF] Mineau | L'apport des SIG en urbanisme (The contribution of GIS to town planning)[END_REF]. In addition to that, the spatial distribution of vulnerability could be improved by increasing the number and nature of input data.

Conclusions

The purpose of this study is to contribute to the existing knowledge and practices that attempt to cope with the problem of increasing urban heat vulnerability of cities. By selecting the Mediterranean city of Montpellier as the case study of the French urban plan-ning context, this research aims to develop a methodology for assessing UHI vulnerability based on the integration of urban planning policy. Actually, the vulnerability index is developed by including various variables and input data related to urban planning regulation, exposure and sensitivity to urban overheating based on an extensive review of the existing literature. Fifteen indicators are weighted by a statistical approach and combined to calculate the global map of UHI vulnerability, which is ranked from very low to very high intensity levels. Indeed, regulation articles from the Local Urban Plans of Montpellier metropolitan area are integrated.

The analysis of the spatial patterns of the vulnerability intensity reveals that the most vulnerable zones are located close to the urban centers of the metropolitan area. In fact, the amount of these zones accounts for 42.2% of the total surface area of the territory. The results also indicate that these affected zones are mainly composed of residential and employment areas. In addition to that, approximately 34.6% of the most vulnerable areas present a medium-to-high level of urban planning challenges. They are also associated to high footprint coefficient, low planted space, alignments at the limits and authorization of installing air conditioning.

At the same time, this research reveals many questions in need of further investigation. In fact, a comprehensive analysis should be conducted by considering more input data associated with LULC databases, transportation, climate and socio-economic dimensions. However, the availability and spatial resolution of data, especially climatic, social vulnerability and economic information at the whole territory level, are often a serious constraint. Additionally, more data from urban regulations could also be incorporated. For example, regulation articles that refer to approved or alternative building materials, which vary in terms of energy efficiency and environmental impact, could be considered.

Figure 1 .

 1 Figure 1. Conceptual framework of UHI vulnerability assessment scheme.

Figure 2 .

 2 Figure 2. (a) Location of Montpellier Méditerranée Metropolitan area in France; (b) 31 municipalities of Montpellier Méditerranée Metropolitan area.

Figure 1 .

 1 Figure 1. Conceptual framework of UHI vulnerability assessment scheme.

Sustainability 2023 ,

 2023 15, x FOR PEER REVIEW 4 of 27 characteristics. Finally, the scientific literature, rich in solutions for mitigating UHI, enables us to develop recommendations for French urban planning.

Figure 1 .

 1 Figure 1. Conceptual framework of UHI vulnerability assessment scheme.

Figure 2 .

 2 Figure 2. (a) Location of Montpellier Méditerranée Metropolitan area in France; (b) 31 municipalities of Montpellier Méditerranée Metropolitan area.

Figure 2 .

 2 Figure 2. (a) Location of Montpellier Méditerranée Metropolitan area in France; (b) 31 municipalities of Montpellier Méditerranée Metropolitan area.

Figure 3 .

 3 Figure 3. Annual mean temperature and number of heat wave days (temperature 5 °C higher than normal for at least 5 consecutive days) from 1980 to 2018.

Figure 3 .

 3 Figure 3. Annual mean temperature and number of heat wave days (temperature 5 • C higher than normal for at least 5 consecutive days) from 1980 to 2018.

Figure 4 .

 4 Figure 4. (a) Building footprint coefficient; (b) proportion of green spaces; (c) building height; (d) implantation of air-conditioning systems; (e) construction's location to plot boundaries.

Figure 4 .

 4 Figure 4. (a) Building footprint coefficient; (b) proportion of green spaces; (c) building height; (d) implantation of air-conditioning systems; (e) construction's location to plot boundaries.

Figure 5 .

 5 Figure 5. (a) Topographical map; (b) proximity to green spaces; (c) proximity to water spaces.

Figure 5 .

 5 Figure 5. (a) Topographical map; (b) proximity to green spaces; (c) proximity to water spaces.

Figure 6 .

 6 Figure 6. (a) SUHI intensity (daytime); (b) SUHI intensity (nighttime); (c) population density; (d) land use classification.

Sustainability 2023 , 27 Figure 7 .

 2023277 Figure 7. (a) Proportion of sensitive population by age; (b) proportion of low-income households; (c) presence of health and care facilities; (d) building age.

Figure 7 .

 7 Figure 7. (a) Proportion of sensitive population by age; (b) proportion of low-income households; (c) presence of health and care facilities; (d) building age.

Figure 8 .

 8 Figure 8. (a) Map of daytime UHI vulnerability; (b) map of nighttime UHI vulnerability; (c) map of global UHI vulnerability.

Figure 8 .

 8 Figure 8. (a) Map of daytime UHI vulnerability; (b) map of nighttime UHI vulnerability; (c) map of global UHI vulnerability.

Figure 9 .

 9 Figure 9. (a) Distribution of urban vulnerability levels of cells affected by medium-to-high classes of "urban planning policies" indicators to total study area; (b) Distribution of urban vulnerability levels of cells affected by medium-to-high classes of "urban planning policies" indicators to the total area of affected cells.

Figure 9 .

 9 Figure 9. (a) Distribution of urban vulnerability levels of cells affected by medium-to-high classes of "urban planning policies" indicators to total study area; (b) Distribution of urban vulnerability levels of cells affected by medium-to-high classes of "urban planning policies" indicators to the total area of affected cells.

Sustainability 2023 , 27 Figure 10 .

 20232710 Figure 10. (a) Distribution of urban vulnerability levels of urban planning zoning to total study area; (b) distribution of urban vulnerability levels of urban planning zoning to the total area of affected cells; (c) urban planning zoning of Montpellier Metropolitan area representing the 4 regulatory land uses.

Figure 10 .

 10 Figure 10. (a) Distribution of urban vulnerability levels of urban planning zoning to total study area; (b) distribution of urban vulnerability levels of urban planning zoning to the total area of affected cells; (c) urban planning zoning of Montpellier Metropolitan area representing the 4 regulatory land uses.

Figure 11 .

 11 Figure 11. (a) Distribution of urban vulnerability levels of authorized land use to total study area; (b) distribution of urban vulnerability levels of authorized land use to the total area of affected cells.

Figure 11 .

 11 Figure 11. (a) Distribution of urban vulnerability levels of authorized land use to total study area; (b) distribution of urban vulnerability levels of authorized land use to the total area of affected cells.

Sustainability 2023 , 27 Figure 12 .

 20232712 Figure 12. (a) Distribution of urban vulnerability levels of the date of approval of Local Urban Plans to total study area; (b) distribution of urban vulnerability levels of the date of approval of Local Urban Plans to total area of affected cells.

Figure 12 .

 12 Figure 12. (a) Distribution of urban vulnerability levels of the date of approval of Local Urban Plans to total study area; (b) distribution of urban vulnerability levels of the date of approval of Local Urban Plans to total area of affected cells.

Table 1 .

 1 Source, spatial resolution, years of collection and method to assign a value to the cells.

	Data

Source Spatial Resolution Years of Collection Method to Assign a Value to the Cells

	Socio-economic	INSEE	200 m	2015	No adaptation
	Land uses	Montpellier open data, aerial view interpretation	50 m	2019-2021	The value of the largest entity contained in the cell is
	Building age	Montpellier open data, CEREMA	50 m	2020	assigned as the value for the cell
	Topography	Digital Elevation Model, IGN	30 m	2016
	Land surface Temperature (LST)	LANDSAT-8 ASTER	30 m	2019 2015	Average LST of each entity contained in the cell
	Urban regulations	Local Urban Plans	50 m	2020	A same value is assigned to all the cells overlapping an area which share the same value of an urban rule
	Health and care facilities	Montpellier open data, aerial view interpretation	50 m	2019-2021	A value of presence is assigned to the cells if it contained an health and care facilities

Table 1 .

 1 Source, spatial resolution, years of collection and method to assign a value to the cells.

	Data

Source Spatial Resolution Years of Collection Method to Assign a Value to The Cells

	Socio-economic	INSEE	200 m	2015	No adaptation
	Land uses	Montpellier open data, aerial view interpretation	50 m	2019-2021
	Building age	Montpellier open data, CEREMA	50 m	2020	The value of the largest entity contained in the cell is assigned as the value for the cell
	Topography	Digital Elevation Model, IGN	30 m	2016
	Land surface Temperature (LST)	LANDSAT-8 ASTER	30 m	2019 2015	Average LST of each entity contained in the cell
					A same value is assigned to all the
	Urban regulations	Local Urban Plans	50 m	2020	cells overlapping an area which share the same value of
					an urban rule
	Health and care facilities	Montpellier open data, aerial view interpretation	50 m	2019-2021

A value of presence is assigned to the cells if it contained an health and care facilities

Table 2 .

 2 Articles of urban regulation of Local Urban Plan.

	Article	Title	Urban Form Elements
	1	Forbidden or allowed land uses	Building
	2	Land uses subjected to particular rules	Building
	3	Access and roads	Road and public area
	4	Technical networks	Road and public area
	5	Minimum parcel area	Parcel
	6	Building location related to road and public area	Road and public area
	7	Building location related to parcel borders	Parcel
	8	Building location related to other buildings	Building
	9	Building footprint	Parcel
	10	Building maximal height	Building
	11	Exterior appearance	Building
	12	Parking	Unbuilt space
	13	Unbuilt space, vegetation and protected woods	Unbuilt space
	14	Floor area ratio	Parcel
	15	Energetic and environmental performance	Building
	16	Communication network and infrastructure	Road and public area

Table 3 .

 3 Urban planning policies indicators and values for each cell.

	Indicators	Definition	Article	Class	Categories	Effects on UHI Vulnerability	Weight Value	References
	Green spaces	Percentage of green areas on the plot	13	1 2 3	65-100% 35-65% 4-35%	Decrease	0.64	[51-55]

Table 3 .

 3 Cont.

	Indicators	Definition	Article	Class	Categories	Effects on UHI Vulnerability	Weight Value	References
	Building Footprint	Percentage of building footprint on the plot	9	1 2 3		5-35% 35-65% 65-100%		Increase	0.43	[56-58]
				1		2-5 m
	Building Height Sustainability 2023, 15, x FOR PEER REVIEW Building maximal height	10	4 2 3		30-50 m 5-10 m 10-30 m		Increase	0.54	[56,57,59] 8 of 27
				5		50-70 m
			5	1	Withdrawal on all 50-70 m			0.20
	Construction's location Air conditioning Construction's location Air conditioning	Construction's implantation on the plot Implantation of systems Construction's implantation on the 6 and 7 plot Implantation of systems 11	6 and 7 11 1 2 3 4 1 2	2 3 4 1 Withdrawal on all 1 withdrawal min. 1 alignment min. Alignment on all Prohibited Increase 1 withdrawal min. 1 alignment min. Alignment on all 2 Authorized Prohibited Increase Authorized	Increase Increase	0.20 0.20	0.20 [60,61] [62,63]	[60,61] [62,63]

Table 4 .

 4 Urban overheating exposure indicators and values for each cell.

	Indicators	Definition	Class	Categories	Effects on UHI Vulnerability	Weight Value	References
			1	>300 m
	Topography	Altitude	2	150-300 m	Decrease	0.22	[14]
			3	<150 m

Table 4 .

 4 Cont.

	Indicators	Definition	Class	Categories	Effects on UHI Vulnerability		Weight Value	References
	Proximity to natural areas	Distance	1 2	Near natural areas Away from natural areas		Decrease		0.90	[51-55,64-66]
			1		<0 • C
			2	0-1 • C
	SUHI intensity (day and night)	Surface Urban Heat Island Intensity	3 4	1-2 • C 2-5 • C		Increase		0.73	[42,43,67,68]
			5	5-10 • C
			6	>10 • C
		1 Sustainability 2023, 15, x FOR PEER REVIEW	0-48				10 of 27
			2	48-126
	Population density	Number of inhabitants	3 4 5	126-231 231-384 1 384-627 2	0-48 48-126	Increase		0.79	[24-29,33]
	Land uses (daytime) Land uses (nighttime)	Population density Majority uses Land uses (daytime) Land uses	6 7 1 Number of inhabitants 2 3 1 2 Majority uses Uses not occupied at day 627-1231 1231-2439 3 126-231 4 231-384 5 384-627 Unbuilt areas 6 627-1231 7 1231-2439 Uses occupied at day Other uses 1 Unbuilt areas 2 Uses not occupied at day Increase 3 Uses occupied at day Uses occupied at night 1 Other uses	Increase Increase	0.83	0.79 0.83	[24-29,33] [14,69,70] [14,69,70]
		(nighttime)		2	Uses occupied at night

Table 5 .

 5 Urban overheating sensitivity indicators and values for each cell.

	Indicators	Definition	Class	Categories	Effects on UHI Vulnerability	Weight Value	References
	Sensitive population	Proportion of children under 10 years and adults over 65 years/ total inhabitants	1 2 3 4 5	0-23% 23-31% 31-38% 38-48% 48-97%	Increase	0.66	[13,27,34,71]
	Health and care facilities	Implantation	0 1	Absence Presence	Increase	0.16	[13,27,34,71]
	Low-income households	Percentage on low-income households/number of households	1 2 3 4 5	0-16% 16-32% 32-48% 48-64% 64-80%	Increase	0.40	[26,72]
	Building age	Average age	0	>1974	Increase	0.73	[24,34]

Table 6 .

 6 Mitigation and adaptation solutions from French urban regulation.

	Article	Title	Decrease Waterproof Areas	Increase Natural Areas	Improve Urban Ventilation	Improve Heat Exchanges	Reduce Anthropogenic Heat Generation
	1	Forbidden or allowed land uses	•	•			•
	2	Land uses subjected to particular rules	•	•			•
	3	Access and roads	•		•
	4	Technical networks					•
	5	Minimum parcel area	•
		Building location
	6	related to road and			•	•
		public area
		Building location
	7	related to			•	•
		parcel borders
		Building location
	8	related to			•	•
		other buildings
	9	Building footprint	•		•
	10	Building maximal height			•	•
	11	Exterior appearance		•		•	•
	12	Parking	•	•		•
	13	Free space, vegetation and protected woods	•	•		•
	14	Floor area ratio	•
		Energetic and
	15	environmental		•		•	•
		performance
		Communication
	16	network and
		infrastructure

• 4.1.1. Solutions to Increase Permeable and Natural Spaces

Acknowledgments:

The authors gratefully acknowledge the support provided by the funders of this research which has enabled and encouraged the production of this article and thank the anonymous reviewers for their helpful suggestions.

Funding: This research was supported by ADEME, the Ecological Transition Agency and the French Ministry of Culture.

Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.

Data Availability Statement: Publicly available datasets were analyzed in this study. This data can be found here: https://data.montpellier3m.fr/; https://www.insee.fr; https://datafoncier.cerema.fr; https://www.usgs.gov; https://search.earthdata.nasa.gov and Local Urban Plans on request from the municipalities of Montpellier metropolitan area.

Conflicts of Interest:

The authors declare no conflict of interest.