Urban Heat Island and Its Interaction with Heatwaves: A Review of Studies on Mesoscale

: With rapid urbanization, population growth and anthropogenic activities, an increasing number of major cities across the globe are facing severe urban heat islands (UHI). UHI can cause complex impacts on the urban environment and human health, and it may bring more severe effects under heatwave (HW) conditions. In this paper, a holistic review is conducted to articulate the ﬁndings of the synergies between UHI and HW and corresponding mitigation measures proposed by the research community. It is worth pointing out that most studies show that urban areas are more vulnerable than rural areas during HWs, but the opposite is also observed in some studies. Changes in urban energy budget and major drivers are discussed and compared to explain such discrepancies. Recent studies also indicate that increasing albedo, vegetation fraction and irrigation can lower the urban temperature during HWs. Research gaps in this topic necessitate more studies concerning vulnerable cities in developing countries. Moreover, multidisciplinary studies considering factors such as UHI, HW, human comfort, pollution dispersion and the efﬁcacy of mitigation measures should be conducted to provide more accurate and explicit guidance to urban planners and policymakers.


Introduction
Due to increasing populations, cities all over the world have been developing both horizontally and vertically to accommodate their increasing needs [1]. Such urbanization has contributed to urban climate change due to the replacement of natural elements (e.g., water, vegetation) with artificial heat-absorbing features (e.g., concrete, asphalt) that have a relatively low albedo, high solar absorption and high heat capacity [1][2][3][4][5][6][7]. These impervious materials can also reduce evapotranspiration, thus reducing latent heat flux and altering the urban thermal environment [4,[7][8][9][10][11][12]. In terms of urban morphology development in the vertical direction, the increasing height of buildings can not only trap the outgoing radiation due to in-canyon reflections but also reduce turbulence transfer and overall ventilation in dense urban environments as a result of increased surface roughness [3,4,7]. Moreover, the released anthropogenic heat from human activities, such as traffic and air conditioning systems, also increases heat accumulation in urban areas [4,6,7]. As a result, cities often present a higher ambient temperature than their surrounding suburban and rural areas, especially at night, and this phenomenon is widely known as urban heat island (UHI) [4,13,14].
UHI has enormous direct and indirect impacts on urban regions and dwellers, such as causing a considerable rise in energy consumption for cooling purposes [4,[14][15][16] and increasing heat-related mortality and morbidity, especially for vulnerable groups such as senior and low-income citizens [4,9,[17][18][19][20]. These impacts are reported to be intensified during heatwave (HW) events, a period of excessively hot days [21][22][23]. Recent research 2 of 26 also shows that the UHI effect can increase the frequency and duration of extreme heat events as it can extend the duration of hot conditions and reduce the time window for relief from high heat exposures [7]. As reported in [24], the UHI effect can lead to a 3.7-day increase in the number of annual HW days in the Yangtze River Delta region in China. More frequent, longer and stronger extreme heat events are expected to occur in the future [6,7,25]. A study conducted by Lhotka et al. [26] predicted that the frequency of HWs would be nearly twice as high in the period 2020-2049 compared to the historical period 1970-1999, and they are projected to occur 3-4 times per year on average by the end of the 21st century. On the other hand, HWs may strongly affect the spatial and temporal variability of UHI by altering sensible and latent heat transfer and changing wind conditions [4,27,28]. While many researchers [27][28][29] have reported synergistic interactions between UHI and HW, other researchers [25,30] found no synergistic effect. Nevertheless, the interaction between UHI and HW is complex and remains a subject of debate.
In order to mitigate the negative effects of UHI, especially during HWs, numerous strategies have been proposed, such as using permeable pavements [31][32][33] and reflective materials [34,35], increasing vegetation cover [36][37][38][39][40][41][42] and so on. Susca et al. [43] evaluated the effectiveness of using vegetation for UHI mitigation in New York City and reported a positive impact since vegetation effectively reduced the building cooling/heating load and the consequent building energy consumption. Although many strategies have been proven to be effective in mitigating or suppressing UHI, it remains unclear whether they are still effective during HW events. For example, considering that grass is sensitive to extreme heat events, its impact on UHI during HW events is still an open question.
Santamouris [4] has highlighted the need for a more integrated and holistic analysis of UHI combining climate issues in cities such as HWs and global climate change in future studies. Due to the severe impact and increasing intensity of UHI and HW, it is necessary to look at the existing studies related to both UHI and HW. This review aims to articulate existing understandings of UHI under HW conditions by answering the following key questions:

•
How does HW affect UHI? • What are the main drivers controlling the synergies between UHI and HW? • What can be done to mitigate UHI under HW conditions?
A comprehensive literature search is conducted to find all peer-reviewed articles related to UHI and HW using three major databases, namely ScienceDirect (access date: 28 April 2021), Scopus (access date: 28 April 2021) and Web of Science (access date: 1 May 2021). The selection process of the literature and the number of articles reviewed are shown in Figure 1. First, different combinations of keywords (urban heat island, heat wave and heatwave) are applied to form advanced query strings for different databases. Then, the review is conducted on the articles that meet the following criteria: published in English; peer-reviewed journals; specifically studying UHII under HW conditions concerning synergistic effect and/or mitigation measures; considering the comparison between HW and non-HW days when studying the synergistic effect; focusing on mesoscale (brown circle in Figure 1). Finally, 64 articles are selected for in-depth review, as is shown in the green circle in Figure 1. Among these articles, 37 of them focus on the synergistic effect; 23 of them focus on mitigation measures; and 4 of them cover both topics.
The remainder of this paper is organized as follows: Section 2 presents the current status of research related to UHI during HWs. Different classifications of UHIs and HWs, interactions between UHI and HW and mitigation solutions are discussed in this section. The research gaps and future research suggestions are discussed in Section 3, and observations from this review are summarized in Section 4.

Current Status of Research
The 64 selected studies are further classified and analyzed. Overall, the number of publications studying UHI under HW conditions has increased over the last decade, as shown in Figure 2, which indicates increasing research attention in this topic. Among these studies, most of them are conducted for Chinese and American cities, as shown in Figure 3, and many of the remaining studies focus on European cities.

Methods for Retrieving UHI
In this section, an overview of the methods adopted in the literature for retrieving UHI is provided. Three major methods have been adopted in this regard, including insitu measurement, remote sensing and numerical modeling [44]. The reviewed articles are classified by research methods in Figure 4. As is shown in this figure, numerical simulation is the most widely adopted method, followed by in situ measurements and remote sensing, while some studies adopt two or more methods. More details can be found in Supplementary materials Table S1. In the subsections to follow, the applications of each of these methods will be briefly discussed.

In Situ Measurement
The in situ measurement approach/method has certain advantages, such as allowing direct, instantaneous and continuous measurements of air temperature, wind speed, humidity and other variables, but it also has some limitations, such as the inability to provide a high spatial resolution of measured variables due to the limited number of weather stations [44]. Ngarambe et al. [45] emphasized the importance of the appropriate selection of urban and rural stations to avoid inaccurate UHI estimations. They suggested that the selected rural area should be located in a relatively flat terrain surrounded by natural properties, without impervious surfaces such as densely constructed buildings. Such a criterion of urban and rural station selection is also adopted in other studies. Jiang et al. [46] and Chew, Liu, Li and Norford [25] considered an urban station in the central urban area with more than 70% impervious surfaces and less than 20% vegetation cover within 1 km 2 of the station. It should be noted that stations in large urban parks have limited representativeness of urban stations [46] since they may underestimate the actual urban heat island intensity (UHII) in cities' built-up zones [47]. In contrast, rural stations were selected with 65% or higher vegetation or cropland cover and 30% or lower impervious surfaces [25,46]. Other criteria, such as population [48] and satellite night-time brightness index [49], have also been applied in some studies. Generally, the selection of weather stations is important when applying the in situ measurement method, and UHII should rather be estimated by fewer representative sites than more unrepresentative sites [50], especially for studying UHI under HWs, as the performance of different land types may change under extreme temperature conditions.

Remote Sensing
In contrast to the in situ measurement, thermal infrared remote sensing from satellites allows global monitoring of geophysical parameters such as land surface temperature, and thus is widely used for studying surface urban heat island [51,52]. A few studies [29,53] also used ground-based or aircraft-based remote sensors to observe the urban boundary  [54][55][56][57][58][59]. These datasets have different spectral, temporal and spatial resolutions, typical ranges of which are shown in Table 1. It should be noted that satellite-based sensors cannot continuously record information over the diurnal cycle due to having an orbital tour around the Earth [44]. In addition, cloud cover may affect dataset quality, and images showing more than 10% cloud cover are usually excluded from analysis [56]. Therefore, although this method can provide continuous spatial data for studying UHI under HWs, it cannot monitor continuous temporal variations in UHI under HWs.  [62], ALARO-0 [63], UrbClim [64], the non-hydrostatic mesoscale atmospheric model (Meso-NH) [65,66], MUKLIMO_3 [67,68] and the Community Earth System Model (CESM) [69], have been developed to simulate, analyze and predict the urban thermal environment and flows. The modeling approach offers a solution to overcome the limitations of in situ measurement and remote sensing as it provides information at high spatial and temporal resolutions and can also test hypothetical scenarios and assess potential future impacts. Promising results have been obtained using the modeling approach in the past few decades [44,70,71]. The details of the various numerical models are summarized in Table 2. These models all comply with the same fundamental mass, momentum and energy conservation laws but differ in hydrostatic assumptions, formulation of equations and available options in parameterization schemes [71]. Among them, WRF, developed by the National Centre for Atmospheric Research, is the most popular one, which can capture temporal-spatial variations in regional climate (temperature, wind, shortwave and longwave radiation, rainfall, surface heat fluxes and other environmental variables) caused by urbanization and has been successfully applied to simulate many UHI and HW events [9,[72][73][74][75].
To precisely resolve the surface energy balance in the atmospheric model, land surface models (LSMs) are usually coupled to provide parameterizations for sub-grid-scale atmospheric processes [71]. These models vary in the level of complexity from simplified bulk representations (such as Noah LSM [76], which treats urban areas as homogeneous and flat surfaces with properly modified radiative, thermal and aerodynamic properties [71,75,77]) to single-and multi-layer urban canopy models (UCMs) that take into account in more detail urban morphology characteristics, building material parameters and anthropogenic heat release [71,78]. Single-layer UCM (SLUCM) [79] considers two-dimensional and symmetrical street canyons with an infinite length. It resolves mass and energy exchanges between the lowest level of the atmosphere and the canopy's roof, wall and road surfaces [6,71,78]. Building Effect Parameterization (BEP) [80], which is a multi-layer UCM (MLUCM), accounts for the three-dimensional nature of urban surfaces. It treats several layers within the urban canopy and simulates the heat, momentum and turbulent kinetic energy exchange across multiple atmospheric levels [37,71,78,81,82]. In WRF simulations, BEP is usually coupled with the Building Energy Model (BEM) [83] to improve the calculation of energy exchange between building interiors and the outdoor atmosphere [78,84]. The BEM, for each floor, considers the conductive/sensible heat transport through walls, Sustainability 2021, 13, 10923 6 of 26 roofs and floors; radiation exchange through windows; longwave radiation exchange between indoor surfaces; production of heat due to occupants and equipment; and air conditioning, ventilation and heating [72,85,86]. It has been demonstrated that BEP-BEM provides more accurate predictions of sensible heat flux than BEP [87]. Princeton UCM (PUCM) [88] is also adopted in WRF simulations. Similar to other UCMs, it considers three main facets-roofs, walls and roads. However, PUCM also includes sub-facet-level heterogeneity for urban fluxes as each facet in the PUCM can be divided into multiple sub-facets with different materials [89,90]. For example, a roof surface can be a combination of a green roof and a cool roof. Moreover, PUCM adopts more realistic representations of hydrological processes, such as in-canyon vegetated soils and water storage capacity for impervious materials [91]. The Town Energy Balance Model (TEB) [92] is usually coupled with a Meso-NH atmospheric model. It considers a simplified two-dimensional building layout, including roofs, walls and roads, and models the thermal exchanges between urban surfaces and the atmosphere. In addition, a thermal comfort model and a building thermal model are included in TEB, which allow the simulation of the universal thermal climate index, estimation of human body thermal comfort and calculation of energy consumption for heating and cooling [44]. Some studies [93,94] combine TEB with the Interaction Soil Biosphere Atmosphere (ISBA) [95] to better represent urban green areas and simulate the interactions between build-up and natural cover. In general, choosing appropriate models to accurately study the interactions between UHI and HW is important, and simulation results should be validated against real measurement data wherever possible. In situ measure [110] WRF-Chem SLUCM, Noah LSM 3000 In situ measure [82] WRF-Chem BEP, Noah LSM 3000 In situ measure [111] WRF-Chem MLUCM, BEM, Noah LSM 1000 In situ measure [112] ALARO-0 TEB 1000 In situ measure  [47] UrbClim N/A 1000 In situ measure [93] Meso-NH TEB, ISBA 1000 In situ measure [94] Meso-NH TEB, ISBA 500 N/A [70] MUKLIMO_3 N/A 100 In situ measure [113] CESM N/A 99,000 × 138,750 Remote sensing

Classifications of Urban Heat Island
There are mainly four types of UHI, as illustrated in Figure 5: subsurface urban heat island (UHI sub ), surface urban heat island (UHI surf ), canopy layer urban heat island (UHI UCL ) and boundary layer urban heat island (UHI UBL ) [114]. Their definitions, processes, data collection methods and impacts are summarized in Table 3. Oke [114] pointed out that the types of UHI measured or modeled should be identified first because each corresponds to different processes and thus requires different measurement techniques. Figure 5 shows that among the reviewed papers, only two types of UHI have been studied, i.e., UHI surf and UHI UCL . Most of the papers have considered the UHI UCL, i.e., between the surface and roof level, such as at 2 m height [25,29,48,73,103,115], since a high air temperature at this level is detrimental to human comfort and wellbeing [57]. Ward, Lauf, Kleinschmit and Endlicher [57] mentioned that surface temperature is the only multi-temporal and spatial information that is widely available for large comparative analyses, while air temperature measurements have many constraints, such as cost, location settings and scale. Nevertheless, studying different types of UHI and their interactions with HW deserves more attention.

Features, Causes and Definitions of Heatwaves
A common feature of all HW events is a persistent high-pressure synoptic system (anticyclone) [116,117]. Other features, such as soil desiccation [48,89], deeper atmospheric boundary layer [89], lower wind speed [45,118,119] and cloudiness [118], are also observed, while some studies report contradictory observations, such as decreased atmospheric boundary layer and abundant soil moisture [29].
Typical HWs occur due to a stubborn jet stream flow, which splits the upper-level atmospheric wind, so that cooler air from the poleward side cannot mix with the hotter air at the equatorial side and warm air starts to build up in the blocked region [116]. Other HWs occur at lower latitudes, such as in Australia [116], where high-pressure systems remain stationary until a low-pressure system shifts them, and they are usually associated with the movement of atmospheric Rossby wave trains [116]. It is reported that anthropogenic activity may cause increased intensity, frequency and duration of HWs, and these three characteristics are predicted to increase remarkably in the future [113,116]. There is no standard and universal definition of HW as different regions have different climatic conditions. Generally, HWs differ in three aspects: temperature metric, intensity (temperature threshold) and duration [120]. In terms of the temperature metric, HWs may be characterized by air temperature, apparent temperature [121] and heat index (HI, see the equation below Table 4). Maximum air temperature is widely used for characterizing HWs in the literature [28,34,53,58,90,99,103,113,118,122,123]. Another temperature indicator is the mean temperature, which is adopted in [30,110]. Combined maximum and minimum air temperatures are also used to characterize HWs [48,112,124]. Other combinations, such as the maximum and mean temperatures, are adopted in certain studies [25]. Compared to the temperature-based indicators, the apparent temperature and HI are reported to be better indicators of HWs because they consider both temperature and humidity, which may provide a better estimation of heat stress on humans, especially in tropical cities [25,120,121,125,126]. High humidity has an adverse effect on the human body's efficiency in discharging heat through perspiration, which in turn results in an increase in the body's internal temperature [29]. Table 3. Classification of UHI [114,127].

Type
Processes Impacts Pyrgou, Hadjinicolaou and Santamouris [119] compared both air temperature and HI indicators under HW events, which are characterized by increased air temperature, lower wind speed and increased absolute humidity. They concluded that using HI to define HWs may predict the intolerable weather conditions in the Eastern Mediterranean region more accurately since high humidity can exacerbate the outdoor thermal comfort in this region. As a result, HWs characterized by HI occur more frequently and last for a longer period of time. An, Dou, González-Cruz, Bornstein, Miao and Li [29] also highlighted the importance of HI since the variation in the daily maximum building energy demand correlates better with the change in the HI than with temperature.
Both absolute and relative thresholds are adopted to define a HW [29,45,120,128]. The latter considers the percentile (90th, 95th, 98th, etc.) over a given period rather than the absolute values (e.g., 32 • C, 35 • C, 37 • C, etc.) to account for the climatic adaptability factor [128]. To a certain extent, these two types of definitions are correlated. Another criterion defining a HW is temporal persistence, i.e., a specified number of consecutive days when the temperature or HI surpasses a certain threshold [45,48]. There is no consensus on the minimum number of days to be considered as a HW [128]. All the HW definitions considering metric, threshold and duration from the reviewed literature are summarized in Table 4. It is clear that most of the studies use 3 days as a criterion to define the minimum duration.

Urban Heat Island Intensity
UHI intensity (UHII) is generally defined as the temperature differences between urban and rural areas [29,48], as shown in Equation (1) [54,57,59,122]: where T urban is the average temperature for all urban stations (in measurements) (T urban = 1 n ∑ n i=1 T i ) or the average of temperatures in grid cells (in simulations) across an urban area, and T rural is the average temperature of reference rural stations (T rural = 1 m ∑ m j=1 T j ) or the average of temperatures in grid cells across a rural area. Certain numerical studies adopt an urban increment method [25,74,130] and calculate the UHII as the difference between a control case and a non-urban case, as shown in Equation (2).
where the control case is the original case of urban areas, and the non-urban case is a hypothetical case with all urban areas replaced by rural areas. A previous study [131] has demonstrated that both methods provide similar diurnal variation of the UHII. However, the latter method can provide the spatial distribution of UHII. Heat magnitude (HM), which refers to the added heat load during HWs to the UHII during average summer or non-HW days [57], is adopted to represent the change in UHII during a HW. HM is calculated using Equation (3): where UH I I HW is the intensity of UHI during HW days and UH I I Re f is the intensity of UHI during a reference period (non-HW days or summer background period). The synergistic interaction between UHI and HW, or, more specifically, whether the joint effect of the two phenomena is greater than the sum of their separate effects [119] or not, is still an open question. A number of studies have explored this question, and their results are discussed below and also summarized in Figure 6 (note that the study in [113] is not included in this figure due to the lack of location information). More details including the location, climate type [132], UHI type, time period, method and HM can be found in Table S1 (in Supplementary Materials).

North America Australia
Europe Asia Most of the studies find that UHII correlates positively with HWs, i.e., UHI is more intense during HW periods than non-HW periods [27,29,47,53,59,74,108,115,117,122,123,129]. This phenomenon is particularly common in major European cities and Chinese metropolitan cities. The HM usually ranges from 0 to 5 K. Among these studies, some [29,59,74,117,122,129] find that the synergistic interaction between UHI and HW is more pronounced in the nighttime, and the UHII can be 2-3 times higher than that over non-HW periods.
In contrast, some studies report that the UHII is not amplified during the HW periods [25,30,128,133]. The diurnal variations of the UHII during the HW period are also similar to those during normal days [25], with the peak UHII recorded at midnight. A few studies [48,56,124,128] even report reduced UHII during HW events. For instance, Kumar and Mishra [56] investigated 89 cities in India and reported that the surface UHII was reduced in the majority of urban areas (63% at daytime and 74% at nighttime) during HWs, the causes of which are discussed in Section 2.4.3.

Energy Balance Models
The origin of a UHI is the change in energy balance due to urbanization [72,114]. As suggested in [114,134], a conceptual control volume, integrating the total urban ecosystem from the top of the roughness sublayer down to a certain depth in the substrate (the vertical heat flux by conduction across the bottom is negligible over the period of consideration), can be defined to study the energy balance. The surface energy balance for such a control volume is illustrated in Figure 7a and described by the following equation [114,[134][135][136]: where Q * is the net all-wave radiative heat flux, Q F is the anthropogenic heat flux released within the volume due to human activities such as traffic, buildings and industry activities, Q H is turbulent sensible heat flux, Q E is latent heat flux, ∆Q S is stored heat flux referring to all energy storage mechanisms within the volume, such as in air, soil, building fabric and trees, and ∆Q A is the net horizontal advective heat flux through the sides of the volume. For rural areas, the counterpart of ∆Q S is the ground heat flux Q G , which represents the sensible heat by conduction to the soil for rural surface energy balance (refer to Figure 7b). where Q * is the net all-wave radiative heat flux, Q F is the anthropogenic heat flux released within the volume due to human activities such as traffic, buildings and industry activities, Q H is turbulent sensible heat flux, Q E is latent heat flux, ∆Q S is stored heat flux referring to all energy storage mechanisms within the volume, such as in air, soil, building fabric and trees, and ∆Q A is the net horizontal advective heat flux through the sides of the volume. For rural areas, the counterpart of ∆Q S is the ground heat flux Q G , which represents the sensible heat by conduction to the soil for rural surface energy balance (refer to Figure 7b).
tical heat flux by conduction across the bottom is negligible over the period of consideration), can be defined to study the energy balance. The surface energy balance for such a control volume is illustrated in Figure 7a and described by the following equation [114,[134][135][136]: where * is the net all-wave radiative heat flux, is the anthropogenic heat flux released within the volume due to human activities such as traffic, buildings and industry activities, is turbulent sensible heat flux, is latent heat flux, ∆ is stored heat flux referring to all energy storage mechanisms within the volume, such as in air, soil, building fabric and trees, and ∆ is the net horizontal advective heat flux through the sides of the volume. For rural areas, the counterpart of ∆ is the ground heat flux , which represents the sensible heat by conduction to the soil for rural surface energy balance (refer to Figure 7b).  [114,137]).
The net all-wave radiative heat flux can be calculated as [114]: where ↓ is incoming shortwave radiation, ↑ is outgoing (scattered and reflected) shortwave radiation, ↓ is incoming longwave radiation and ↑ is outgoing longwave radiation. As shown in Figure 7c,d, the main difference between daytime and nighttime surface radiation is with the shortwave radiation. The shortwave radiation comes from the sun, The net all-wave radiative heat flux can be calculated as [114]: where K ↓ is incoming shortwave radiation, K ↑ is outgoing (scattered and reflected) shortwave radiation, L ↓ is incoming longwave radiation and L ↑ is outgoing longwave radiation. As shown in Figure 7c,d, the main difference between daytime and nighttime surface radiation is with the shortwave radiation. The shortwave radiation comes from the sun, and the received amount at a specific area depends on various factors, including the latitude, longitude, time period, the ability of the atmosphere to transmit and absorb shortwave and the ability of air, clouds and surrounding surfaces to scatter and reflect radiation [114]. It is reported that the incoming solar radiation increases during HWs due to less cloudiness [25]. Other studies observe a decrease in K ↓ in cities with significant photochemical smog or industrial aerosol [135]. The outgoing shortwave radiation K ↑ is primarily determined by surface reflectance, characterized by the albedo value. The incoming longwave radiation L ↓ is the radiation emitted within and from the Earth-Atmosphere System [114]. The value of incoming longwave radiation depends on the ability of the atmosphere and surrounding surfaces to emit and reflect [114]. It is generally agreed that L ↓ will increase in urban areas due to the higher air temperature (UHI) or enhanced atmospheric emissivity resulting from particulate and gaseous pollutants, but it may be offset by the effect of the higher surface temperature on L ↑ [135].
One of the main differences between urban and rural areas is the anthropogenic heat flux. Q F is the heat that is mainly converted from chemical or electrical energy and released into the atmosphere due to human activities in urban areas [114]. It includes waste heat from transport-related fuel combustion, building heating, cooling and ventilation, The net all-wave radiative heat flux can be calculated as [114]: where K ↓ is incoming shortwave radiation, K ↑ is outgoing (scattered and reflected) shortwave radiation, L ↓ is incoming longwave radiation and L ↑ is outgoing longwave radiation. As shown in Figure 7c,d, the main difference between daytime and nighttime surface radiation is with the shortwave radiation. The shortwave radiation comes from the sun, and the received amount at a specific area depends on various factors, including the latitude, longitude, time period, the ability of the atmosphere to transmit and absorb shortwave and the ability of air, clouds and surrounding surfaces to scatter and reflect radiation [114]. It is reported that the incoming solar radiation increases during HWs due to less cloudiness [25]. Other studies observe a decrease in K ↓ in cities with significant photochemical smog or industrial aerosol [135]. The outgoing shortwave radiation K ↑ is primarily determined by surface reflectance, characterized by the albedo value. The incoming longwave radiation L ↓ is the radiation emitted within and from the Earth-Atmosphere System [114]. The value of incoming longwave radiation depends on the ability of the atmosphere and surrounding surfaces to emit and reflect [114]. It is generally agreed that L ↓ will increase in urban areas due to the higher air temperature (UHI) or enhanced atmospheric emissivity resulting from particulate and gaseous pollutants, but it may be offset by the effect of the higher surface temperature on L ↑ [135].
One of the main differences between urban and rural areas is the anthropogenic heat flux. Q F is the heat that is mainly converted from chemical or electrical energy and released into the atmosphere due to human activities in urban areas [114]. It includes waste heat from transport-related fuel combustion, building heating, cooling and ventilation, industrial manufacturing, human metabolism, etc. [72,136]. All of them can be categorized into three sectors: transport, buildings and human metabolism [29,114]. It is reported that Q F from the building sector accounts for the major proportion of the total Q F [72]. However, Q F is difficult to accurately estimate since it is affected by urban morphology, weather conditions, varying human activities and technology adoption, etc. [72].
The turbulent sensible heat flux Q H and latent heat flux Q E are mostly investigated and considered as primary causes of urban and rural temperature variations [28,128]. The Bowen ratio, β, which is defined as Q H /Q E , is an important indicator adopted in UHI studies. β > 1 indicates that more heat is converted into sensible form, which warms the lower atmosphere, and β < 1 indicates that latent heat flux is the dominant one, which cools the lower atmosphere and brings humidity to the environment [114]. Latent heat flux is usually higher in rural areas than in urban areas because of abundant vegetation [128].
The storage heat flux ∆Q S is usually higher in urban areas than that in rural areas due to the higher heat storage capacity of building/pavement materials compared to vegetation and soil [10,128]. In urban areas, heat is usually stored during daytime and mostly released as sensible heat flux at nighttime since urban areas are mostly covered with impervious surfaces [25]. Therefore, UHII is expected to be enhanced during HWs [34,90,[138][139][140][141].
The advective heat flux ∆Q A is the net energy carried through the control volume by the wind-driven transport of sensible and latent heat [114]. Some studies attempt not to consider ∆Q A by assuming a relatively homogeneous urban surface [114,136]. However, the horizontal properties are usually different in urban areas, and the thermal gradient between urban and rural areas can result in secondary circulation [89]. In addition, specific morphologies such as hills/mountains and lakes/sea can significantly influence the (cooling) air flow towards the city, and thus influence urban temperatures. The effects of the ventilation corridors, such as rivers and large boulevard streets, are also remarkable. The study conducted by Ching et al. [142] revealed that the horizontal heat fluxes may dominate vertical fluxes when there is a strong horizontal temperature gradient.

Variation in Energy Budgets and Corresponding Drivers
As mentioned above, the relative warming of the urban environment is related to the energy balance [122]. A number of drivers controlling the energy budget change can be distinguished: atmospheric stability, wind speed, air pressure variation, clouds, evapotranspiration, surface albedo, exhaust gases from traffic, pollution, soil moisture, etc. [48,54,81,103,108,122,143]. The relationships between the energy budget and possible drivers are illustrated in Figure 8.
Theoretical models have been developed to quantify the relationships between different drivers. For instance, a biophysical factorization method, which decomposes the urbanization-induced changes in the surface air temperature (air temperature near ground) into multiple components, was proposed in [73,144]. Chew, Liu, Li and Norford [25] pointed out that both momentum and thermal roughness lengths and the surface resistance are important contributing factors and have developed correlations to quantify their effects accordingly based on Monin-Obukhov similarity theory [145]. Evapotranspiration is considered one of the significant drivers in lowering temperature [119,146], and dew formation is another important driver since the evaporation of night dew during the day may bring about some cooling effects [119,128]. This phenomenon occurs when the surface temperature is less than or equal to the dew point temperature, which correlates with the atmospheric pressure and humidity [119]. Efficient nocturnal radiation cooling occurring in clear sky and light wind conditions is also reported to be an important factor for dew formation [119,147].
The synergistic interactions between UHI and HW are attributed to the urban-rural contrasting responses to different drivers and the consequent changes in the energy budgets [29,73]. The findings of the studies investigating different energy budget changes and relevant drivers are summarized in Table S2. Three distinct scenarios have been reported in the literature, which are described below.

• UHII increases during HWs
Most studies have observed the synergies between UHI and HW. The net all-wave radiation is generally higher under HWs than that under normal conditions due to the increase in incoming shortwave radiation in clear sky or lower water vapor conditions [28,29,46,148]. Incoming longwave radiation was also found to increase by Li, Sun, Liu, Yang, Wang and Gao [28]. Jiang, Lee, Wang and Wang [46] found that the daily total solar radiation during HWs was approximately 1.5-times higher than that in non-HWs in all three cities studied. In [28,115], the differences between urban and rural areas were investigated and it was found that the difference in the net all-wave radiation becomes larger during HWs than that under normal conditions, suggesting that urban areas receive more radiative energy input than rural/suburban areas. In contrast, He, Wang, Feng, Yan, Miao, Zhang and Xia [73] reported that the difference in the surface net radiation between HW and regular episodes is negligibly small. They argued that the downward shortwave radiation during daytime is affected more by urbanization on both HW days and regular days.
Enhanced anthropogenic heat flux during HWs is reported to be one of the key contributors to the synergistic effect at both daytime and nighttime [29,73,74,115]. This is primarily attributed to the increased energy consumption by air conditioning systems in urban areas during HWs [73,115]. This increased energy input will be further released into the atmosphere, thus also altering the global urban temperature.
The increase in urban-rural difference in latent heat flux is suggested to be another important cause leading to the increase in UHII during HWs [74,90]. In urban areas, urbanization leads to a greater amount of input energy being transformed into sensible heat and reduced latent heat due to limited vegetation cover and soil moisture [73,108,115]. This effect becomes more significant during HW events, as observed in [73]. In rural areas, latent heat flux by evapotranspiration increases due to abundant vegetation and available water [28,118]. Oláh [149] also pointed out that vegetation can utilize incoming solar radiation for photosynthesis, and thus reduce the heat transformed by radiation. During HWs, the higher atmospheric temperature could enhance the capacity of air to absorb water vapor and the potential evapotranspiration mechanisms [45,119]. Pyrgou, Hadjinicolaou and Santamouris [119] observed higher dew formation during HWs and reported that the dew deposited on the upper top layer of soil may suppress higher temperatures in rural sites by soil evaporation. Ngarambe, Nganyiyimana, Kim, Santamouris and Yun [45] and Ao, Wang, Zhi, Gu, Yang and Li [115] reported that the urban-rural difference in humidity becomes larger during HWs. Studies conducted by Founda and Santamouris [123] found that the higher water content in rural areas is attributed to sea breeze circulation; they also highlighted the significant effect of background wind field on the synergistic interactions.
Most studies [29,45,54,103,108,115,119,122,129] observe a lower wind speed during HWs than non-HW periods, where the reduced wind speed increases atmospheric stability, prolongs stagnant conditions and weakens horizontal advective cooling. Ngarambe, Nganyiyimana, Kim, Santamouris and Yun [45] found higher UHII during low wind speed conditions than during medium and high wind speed conditions. They also found that the synergies between UHI and HW are more intense in densely built areas, which is attributed to the presence of tall buildings hindering natural wind flow, resulting in reduced heat removal. Prominent wind directions are observed to change during HWs in the studies conducted by Jiang, Lee, Wang and Wang [46] and Ramamurthy and Bou-Zeid [89]. Jiang, Lee, Wang and Wang [46] observed wind originating from neighboring hot cities in the southwest instead of the cool sea breeze from the southeast during HWs, which leads to a significant increase in the daytime UHII. Reduced moderation of UHI by the sea breeze due to the change in wind direction was also reported by Ramamurthy and Bou-Zeid [89].
A significant increase in heat storage in urban areas during HWs was observed by Li, Sun, Liu, Yang, Wang and Gao [28], the additional heat gain of which was released during night, resulting in more significant nocturnal synergistic interactions [45,74,103,123]. Ramamurthy and Bou-Zeid [89] revealed that the variations in heat storage among different cities are related to the size of the cities. They found that New York City, where the UHI is amplified during HWs, has around 3.5 GW heat stored, compared to 20 MW for Philadelphia.

• UHII decreases during HWs
In contrast to the studies described above, some studies find that the intensity of UHI decreases during HWs, and this is mainly caused by reduced latent heat flux in rural areas due to soil moisture depletion, which results in a decrease in evapotranspiration [48,56,128]. Urban areas are mostly covered by impervious surfaces, and thus the evapotranspiration is less likely to have an impact on the urban temperature [48]. Richard, Pohl, Rega, Pergaud, Thevenin, Emery, Dudek, Vairet, Zito and Chateau-Smith [48] studied six HW periods and found that little or no rainfall was recorded during five of them. Such low rainfall leads to a decrease in soil wetness and consequent evapotranspiration. Kumar and Mishra [56] revealed that the soil moisture and vegetation cover decline after harvesting since most of the rural areas are located in agriculture-dominated regions. Khan, Paolini, Santamouris and Caccetta [128] found that the fast depletion of moisture and latent heat flux in rural areas during daytime is caused by the advection of dry and warm wind blowing from the desert, while the urban area is affected by the advection of humid and cool air from sea at the same time. Therefore, the advective heat flux is recognized as a dominant factor in addition to the latent heat flux, especially in coastal cities [128].

• UHII remains unchanged during HWs
Although most studies observe a lower wind speed during HWs, which results in amplified UHII, Chew, Liu, Li and Norford [25] found no significant deviation in the wind speed and diurnal profile compared to pre-HW and post-HW days. Since no change in wind speed is observed, the urban-rural soil moisture deficit is not significantly altered, and the heat storage is not amplified. Consequently, the UHII remains unaltered during a HW episode, as reported by Chew, Liu, Li and Norford [25]. However, a limitation of their numerical modeling is that they assume the same anthropogenic heat on HW days as on normal days, which may affect the validity of their results. Ramamurthy and Bou-Zeid [89] found that UHII remains unchanged in small cities, and they attribute this to the low thermal capacity and small soil moisture gradient across the short urban-rural transect of such cities.

Mitigation of Urban Heat Island during Heatwaves
To alleviate extreme urban temperatures, various mitigation strategies and their effectiveness for different regions have been studied. The strategies discussed in this section include increasing vegetation coverage (e.g., green roof and park), adopting high albedo surfaces (e.g., modified radiative properties of walls, roofs and roads) and adding water features (e.g., irrigation), etc. [70,96,105,150]. Since the occurrence of HWs may become more frequent, understanding the impact of these mitigation solutions in diverse cities and extreme weather conditions is important. All the reviewed studies in mesoscale related to mitigation solutions of UHI under HW conditions are summarized in Figure 9 and Table S3, and the various strategies are described below.

Concerns over Heatwaves
Light Pollution Wilting and Thermal Stress Thermal Stress, Availability & Cost Figure 9. Summary of studies related to UHI mitigation measures under HW conditions (ET. is abbreviation for evapotranspiration).

Increasing Albedo
Albedo refers to the ratio of the reflected solar radiation from a surface to the incoming solar radiation on that surface [77,100]. Lower albedo leads to reduced reflection and increased absorption of solar radiation, and consequently to a higher surface net radiation and temperature. In many cities, the roofs, walls and roads are dark-colored, which increases net radiation and leads to higher temperatures due to their low solar radiation reflectance [106]. The mean albedo of urban surfaces is typically around 0.2, which means that approximately 20% of the incoming solar radiation is reflected [98]. For example, asphalt, which is the main material for streets, has albedo values between 0.12 and 0.24 [100]. To mitigate UHI, deliberately increasing the albedo of urban surfaces, such as having cool/white roofs, high reflective walls and roads, has drawn widespread attention.
Materials with higher albedo usually have a lighter color compared with ordinary materials. They can reflect more of the incoming solar radiation in urban areas and thus reduce the net radiation Q * [84,91,96,105]. Ultimately, the sensible heat flux Q H , which heats up the urban atmosphere, is reduced, and the urban ambient temperature is potentially lower. The study conducted by Ma, Shao and Song [106] revealed that around 183.3 W/m 2 net solar radiation and 125.8 W/m 2 sensible heat flux can be reduced when the roof albedo is increased from 0.15 to 0.8. Imran, Kala, Ng and Muthukumaran [6] found that cool roofs can reduce net radiation by up to 100, 160 and 180 W/m 2 as well as sensible heat flux by up to 100, 170 and 220 W/m 2 by changing the albedo values from 0.2 to 0.5, 0.7 and 0.85, respectively. It should be noted that, with reduced surface temperature, the outgoing longwave radiation is also reduced. However, the net radiation still decreases since the increase in the reflected shortwave radiation outpaces the reduction in outgoing longwave radiation [91]. Previous studies have also reported that increasing albedo can reduce building energy consumption for cooling due to the reduced urban ambient temperature and the reduced heat storage ∆Q S in roofs, and can consequently reduce the anthropogenic heat emission (e.g., heat pumped into the exterior by air conditioners) in urban areas [6,91,98,106] if solar shading is handled properly; otherwise, as indicated in [151], increasing albedo may lead to an increase in solar heat gain through windows and consequently an increase in the building cooling load.
Existing studies suggest that increasing albedo is an effective way of reducing the urban heat island during summer extreme heat events [94,107,111], although it may not completely offset the HW effect [100,106]. HWs generally have more cloud-free conditions than non-HW periods, and thus allow more incoming solar radiation. Cities are likely to receive more effective cooling outcomes under such conditions with increased albedo [96,98]. Chen, Zhang, Zhu, Wang and Skorokhod [110] found that the UHII is decreased by 28% on normal days and by 29.9% on HW days during summer when albedo increases from 0.2 to 0.85. Typical values of the albedo ranging from 0.12 to 0.88 are reported to have positive relationships with temperature reduction. Geletič, Lehnert and Jurek [70] reported that the daily average air temperature in both Prague and Brno is reduced by 0.2 K with an albedo increase of 0.25, while the temperature difference can be up to 0.5 K in Prague and 0.4 K in Brno when albedo is increased by 0.5. Imran, Kala, Ng and Muthukumaran [6] found that the maximum roof surface temperature can be reduced by up to 3, 4 and 5 K in urban areas with albedos of 0.5, 0.7 and 0.85, respectively. Zhou and Shepherd [105] noted that albedo has a moderate effect, with around a 2 K reduction in the 2 m air temperature when it is increased from 0.15 to 0.45, while the albedo has no effect on lowering the temperature at the 2 m height in Atlanta when it is increased from 0.15 to 0.3.
Previous studies have reported that high albedo has more distinct effects on reducing surface temperature than reducing air temperature because the surface is directly affected by the distinct reduction in the net solar radiation [106]. The decrease in the surface temperature can be up to 6.2 K, which is much more than the decrease in the air temperature (1.5 K), as reported by Ma, Shao and Song [106]. The same conclusions are also found in [6,91]. Some studies also compared the daytime and nighttime effects of such mitigation strategies and reported that increasing albedo is most effective during daytime as its cooling mechanism occurs via the reflection of sunlight [84,96,98,111]. At nighttime, altering albedo no longer plays an important role in the surface energy balance due to the absence of incoming solar radiation [91]. Nevertheless, it is interesting to observe that nighttime temperature is slightly lower with higher albedo according to [77,85,100,106,107], due to reduced urban heat storage during the day [102].
Air speed is also reported to be reduced by increasing urban albedo during HWs [6,82,100,106,109]. This is mainly because the lower surface temperature caused by increasing albedo leads to a reduced buoyancy effect and in turn weakens turbulence [82]. It is found that the turbulent kinetic energy is reduced with increasing albedo [6]. Meanwhile, it is also observed that the planetary boundary layer depth is reduced by 300 m over urban areas [6]. It is interesting to note that some studies report a small temperature increase in discrete areas, and they speculate that it is likely due to the changes in vertical and horizontal mixing associated with the increase in albedo [106].

Increasing Vegetation Coverage
Vegetation changes the surface energy balance by increasing the latent heat flux via evapotranspiration relative to the sensible heat flux for a given net radiation [6,9,57,91,[152][153][154][155]. The changes in surface energy balance lead to a decrease in the surface temperature and the near-surface air temperature [34,102]. Vegetation also reduces the amount of solar radiation transmitted to underlying surfaces by providing shade and thus cools the soil and, as a result, also the air [105,156]. Therefore, increasing vegetation is regarded as an effective mitigation strategy to reduce UHI.
Strategies for increasing vegetation coverage in urban areas include applying green roofs, increasing the number of trees and grassy areas and adding urban parks, etc. In general, it is found that increasing vegetation is effective in reducing the temperature during HW days [6,91,99,104,105]. Imran, Kala, Ng and Muthukumaran [6] reported that the daily sensible heat flux can be reduced by up to 150 W/m 2 by increasing the green roof fraction from 0 to 90%, and the roof surface UHI and near-surface UHI can be reduced by 3.8 K and 1.15 K, respectively. Sun, Grimmond and Ni [104] revealed that the temperature reduction has a strong linear relationship with the green roof coverage fraction. A decrease of 0.16 K and 0.22 K in the mean and maximum 2 m air temperature can be achieved by a 10% increase in the green roof coverage fraction. Imran, Kala, Ng and Muthukumaran [101] investigated three different combinations of vegetation types: mixed forest, mixed forest and grasslands and mixed shrublands and grasslands. They found that all these methods have a positive impact on reducing surface and air temperatures. The mixed forest and grasslands method is slightly more effective in mitigating UHI than the other two methods. The cooling effect of increasing vegetation coverage during the night is observed in [82,91,101,104,105]. As revealed by Imran, Kala, Ng and Muthukumaran [101], this phenomenon can be explained by the reduction in heat storage during the daytime since UHI is driven by the release of heat stored during the night.
In contrast, Alavipanah et al. [157] noted that the effectiveness of vegetation is less pronounced during HWs compared to other periods. They explained that it is mainly because vegetation dries out during HWs due to water shortage. The importance of soil moisture in the performance of vegetation is also highlighted in the literature [57,91,158]. Li, Bou-Zeid and Oppenheimer [91] revealed that the cooling effect provided by green roofs may be eliminated when the soil moisture is very close to its wilting point.

Irrigation
As mentioned above, soil moisture plays an important role in evapotranspiration, which is a combination of evaporation and transpiration [9,159,160]. The rate of evapotranspiration is limited by the water availability, and it affects the cooling effect by vegetation [9,161]. This highlights the importance of irrigation.
The effectiveness of irrigation strategies during HWs has been studied by Gao, Santamouris and Feng [9] and Daniel, Lemonsu and Viguié [93]. Both studies show that cooling by irrigation is more effective during HW periods compared with non-HW periods. The daily average latent heat flux is significantly increased and the urban air temperature is reduced with irrigation [9]. Daniel, Lemonsu and Viguié [93] observed that the efficacy of vegetation watering based on the 2 m air temperature slightly increases with increasing duration and intensity of HWs. Apart from the ambient temperature, the study conducted by Gao, Santamouris and Feng [9] revealed that the cooling effect of irrigation is also sensitive to other environmental variables, such as the soil moisture level before irrigation, urban fraction (fraction of impervious surface in an urban area) and sea breeze. It is found that the effectiveness of irrigation is enhanced when soil is drier before irrigation, but it decreases with a higher urban fraction. It is also found that irrigation is less effective in coastal regions affected by the sea breeze.
It should be noted that the increase in humidity by irrigation may have an adverse effect on reducing the thermal stress in hot and humid regions and in turn on outdoor human thermal comfort [162]. Daniel, Lemonsu and Viguié [93] found that the performance of irrigation in mitigating heat stress is slightly reduced under very severe HW conditions.

Combined Strategies
Some studies have investigated the combination of different mitigation strategies. For instance, Masson, Lion, Peter, Pigeon, Buyck and Brun [94] studied the effects of landscape change (forest, local farming and water) coupled with reflective surface coverage and found that the UHI in Paris can be reduced by 2-3 K. Zhou and Shepherd [105] revealed that the cooling effect of combined albedo change and green space is not as significant as individual strategies during HW days. The combination of five different urban hydrological processes, including anthropogenic latent heat, urban oasis effect, urban irrigation, evaporation from paved surfaces and green roofs, was studied by Wang, Liu, Miao, Wu, Zhang and Qiao [97]. They found that the urban hydrological processes can effectively mitigate UHI surf and UHI UBL , and the effectiveness is more distinct with reduced density of urbanization. However, they also pointed out that the reduced vertical airflow speed due to the mitigation method may inhibit moisture transport to the higher atmosphere and in turn enhance the wetting effect in the lower PBL, which may partly counteract its cooling effect on human comfort.

Addressing the Important Questions
The primary goal of this review is to articulate existing understandings of UHI under HW conditions. Although understanding the synergies between UHI and HW is still an open research issue, some general findings can be summarized, which provide insight into the three key questions outlined in the Introduction.

•
How does HW affect UHI?
HWs usually result in higher temperatures in both urban and rural areas. The intensity of UHI depends on the background climate and the different responses to HWs between urban and rural sites. Most studies note intensified UHII during HWs, while other studies find no apparent synergies between them or even a reduced UHII in HWs.

•
What are the main drivers controlling the synergies between UHI and HW?
The synergies between UHI and HW are influenced by various factors, such as the morphology, population, climate background and many more [29,45]. For the daytime, the main changes in the energy budgets identified in most studies are the differences in latent heat flux and net all-wave radiation between urban and rural areas, and the relevant drivers include vegetation fraction, evapotranspiration, soil moisture (precipitation) and albedo, etc. For the nighttime, the main changes lie in the differences in the heat storage fluxes (to be released at night) between urban and rural areas, and the relevant drivers include city size and albedo, etc. Increased anthropogenic heat flux due to greater usage of air-conditioning systems and changes in wind conditions can also significantly amplify both daytime and nighttime UHII.

•
What can be done to mitigate UHI under HW conditions?
Several strategies have been proposed to mitigate UHI under heatwave conditions. First, increasing surface albedo can significantly reduce the net solar radiation on surfaces and thus lower the urban ambient temperature. It should be noted that high albedo values above 0.8 are only applicable for clean coated surfaces and may cause visual discomfort (glare), and thus intermediate albedo values around 0.5 are more practical [106]. Second, increasing vegetation coverage in urban areas may be effective, but city planners should consider irrigation availability. In addition, the increased humidity due to irrigation may reduce human comfort in hot and humid regions. Finally, the combination of different strategies may reduce UHII, but their effectiveness needs to be further investigated.

•
Time duration: Some studies, such as [25,108], only consider a short time period or even a single heatwave event, which may not be sufficiently representative of longterm behavior. However, the evolution of climate change is a long-term and complex process, and HWs are reported to become more intense and last for longer durations in the future [57]. Therefore, future studies exploring the synergies between UHI and HW and mitigation strategies should consider multi-year datasets or even future weather datasets. • Study area: The synergistic effect between UHI and HW has been quantified mostly for metropolitan cities, but few studies have considered the different responses of urban and non-urban sites to HWs. Future studies should focus more on developing countries and the comparison between urban and rural sites. More detailed classifications such as low-intensity urban areas, high-intensity urban areas, forest, grassland and croplands, may be considered and compared. Moreover, although individual cities in different regions are worth investigating due to their unique geographies and climates [163], further studies may pay more attention to the comparisons among cities of different sizes or with different climates to reveal the general trend.

•
Research method: In terms of numerical simulation, more complex UCMs such as MLUCMs may be preferred over the SLUCM for better prediction [25,84,164,165], but the latter is computationally inexpensive. In addition, high-resolution downscaling methods [166] may be considered in further studies. In terms of onsite observations, unified and dense weather stations can be established to improve spatial resolution. The data obtained through onsite measurements can be used to validate the numerical simulations. It is also advisable that a hybrid method incorporating both numerical simulations and onsite observations may be adopted alongside machine learning methods to construct accurate prediction models [167]. • Data collection: Some studies, such as [25,89,105], do not consider the anthropogenic heat flux or its variations during HWs. Future studies should include this as it plays an important role in the interactions between UHI and HW. • UHI types: Despite a great deal of studies on UHIs during HWs, few studies have considered the differences in the thermal conditions, such as the temperature between urban and rural areas in the entire urban boundary layer. In addition, the relationships between the different UHI types also need to be investigated as studying the connections between air temperature and surface temperature may enable us to better understand the heat transfer in the atmospheric boundary layer and propose effective mitigation methods. • Efficacy of mitigation strategies: When considering vegetation as a strategy, factors such as the efficiency of water use for different irrigation methods must be considered. Moreover, the choice of optimal vegetation type for different urban morphologies should be studied for cities in different climates. The sensitivity to soil moisture is worth investigating to provide more information about variations in latent heat flux. Mitigation by highly reflective surfaces may be effective in summer but may be undesirable in winter. The efficacy of these strategies for the mitigation of the HW effect and their impact in various seasons and in future climates should be investigated. More studies on the efficiency of combined mitigation strategies involving various physical mechanisms and processes should be conducted. In addition, studies on the efficiency of lightweight building materials may be conducted as they have a smaller heat capacity, which may reduce heat release at night. It is also worth noting that the efficacy of mitigation strategies may need to be assessed with reference to the 'absolute' urban thermal climate in addition to UHII due to the large diurnal variation in UHII and other factors influencing thermal sensation, as suggested by Martilli et al. [168]. • Health, economy and environmental considerations: Future work may consider other environmental and socio-economic impacts of the mitigation strategies, such as their impact on air quality, storm water, energy use, health, thermal comfort, visual comfort (albedo change), social aspects as wellbeing, costs and so on.

•
Other considerations: A holistic approach considering changes in both climate and urbanization (horizontal and vertical spatial extent, building size, building material) may be adopted in future studies. Ramamurthy and Bou-Zeid [89] pointed out that the soil state and properties in urban areas are different from natural soils, which leads to variability in their thermal and hydrological characteristics. Accordingly, more research effort may be devoted to this topic.

Concluding Remark
HWs are expected to occur more and more frequently. Exploring the synergies between UHI and HW is necessary as it helps to identify the key drivers that may intensify urban overheating and thus may be helpful for the development of adequate mitigation strategies. The present review provides a comprehensive overview of studies on the interactions between UHI and HW and the mitigation methods during HW days. It is found that, although a variety of research topics have been covered in the open literature, a number of research gaps, such as overlooking the long-term effect of heatwave events and failure to consider different UHI types and their connections, etc., still exist and deserve more research attention. It should be noted that this review only considers papers published in English and included in the three main databases (i.e., ScienceDirect, Scopus and Web of Science). There may be other studies related to this topic but not included in this review due to the adopted selection criteria. Nevertheless, the majority of the studies on this topic have been covered in this study.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/su131910923/s1, Table S1: Summary of the reviewed articles related to studying UHI in HW, Table S2: Summary of studies investigating energy budget changes and relevant drivers, Table S3: Summary of the reviewed studies related to mitigation solutions of UHI in HWs.