Assessing the Spatial Mapping of Heat Vulnerability under Urban Heat Island (UHI) Effect in the Dhaka Metropolitan Area

: The urban heat island (UHI) phenomenon gets intensiﬁed in the process of urbanization, which increases the vulnerability of urban dwellers to heatwaves. The UHI-induced vulnerability to heatwaves has increased in Bangladesh during past decades. Thus, this study aims to examine the UHI and vulnerability to heatwaves in the city of Dhaka using a heat vulnerability index (HVI). The HVI is constructed using various demographic, socioeconomic, and environmental risk variables at thana level. Principal component analysis (PCA) was applied to the 26 normalized variables for each of the 41 thanas of Dhaka to prepare the HVI. Result shows that more than 60% of the city is under built-up areas, while vegetation cover and water bodies are in low proportion. Analysis of HVI shows that the very high- and high-risk zones comprise 6 and 11 thanas, while low- and very low-risk zones comprise only 5 and 8 thanas. The correlation of HVI with variables such as exposure (0.62) and sensitivity (0.80) was found to be highly positive, while adaptive capacity had a negative correlation ( − 0.26) with the HVI. Findings of this study can be utilized in the mitigation of UHI phenomenon and maintaining the thermal comfort of Dhaka.


Introduction
During urbanization, pervious lands have been transformed into impermeable surfaces, producing several environmental hazards in large metropolitan areas [1]. Among the several negative repercussions of urbanization, the most widely recorded occurrence in cities is the urban heat island (UHI), which develops because of growing heat-absorbing pavement materials in cities [2][3][4]. UHIs have been linked to regional land-use changes [5][6][7], exacerbated by absorbed heat from greenhouse gas emissions [8]. Furthermore, the UHI phenomenon is caused by increased building heights in cities, which block breezes and significantly restrict ventilation [9,10]. According to studies, the average daily summer temperature in urban areas greater than 500 square kilometers increased by 4.7 degrees Celsius across three annual cycles from 2003 to 2005 [11,12]. Air temperatures in the United States' urban and rural areas have risen by 0.24 and 0.16 degrees Celsius per decade [13]. Similarly, the temperature of cities in south and east Asian countries, such as China, India, Bangladesh, and Pakistan, has risen by 0.23-0.57 degrees Celsius every decade in recent focused on a detailed examination of UHI intensity and spatiotemporal pattern of UHI, finding that the urban center has the highest UHI intensity while the suburb region appears to have the lowest [51,52]. Cross-sectional research in the second category has aimed to study the relationship between UHI and land cover patterns [6,7], land-use type [53], and urban structure [54]. Simulation models, such as weather research and forecasting (WRF) and ENVI-met, have been employed in the third category to predict and minimize future UHI intensity. However, no research on the impact of UHI on human health has been undertaken in the Dhaka Metropolitan Area, although UHI has been associated with poor human health in other Asian cities [55]. Thus, this study aims to use a heat vulnerability index to investigate the vulnerability to heatwaves in the Dhaka urban heat island using principal component analysis (PCA). Although UHI and HVI studies aim to examine the impact of heatwaves and rising temperatures on humans, most of these studies lack socioeconomic and environmental variables. The HVI was calculated in this study using multiple normalized socioeconomic, demographic, and environmental variables. It is believed that by combining the UHI with socioeconomic, demographic, and environmental aspects, more accurate results will be achieved, providing a clearer picture of the sensitivity to heatwaves. The outcome of this study will aid in reducing the vulnerability to heatwaves and sustaining urban thermal comfort.

Description of the Study Area
The Dhaka Metropolitan Area (DMA) was selected as the study area for research ( Figure 1). DMA, which is situated on the eastern bank of the Buriganga River, is, at present, projected to have a population of 21,741,000 and an area of 2161.17 sq. km [56]. The geographical location of DMA is at 23 • 43 0" North latitude and 90 • 24 0" East longitude. The DMA area is covered by a total of 41 thana and the study has been built up with the local level variation of urban heat islands, which is also an ideal example for showing the variations of urban heat vulnerability. Dhaka has been experiencing heatwaves for recent years and this study has been conducted to analyze vulnerability on a macro scale. Urban areas have increased by 20.52% in the DMA from 2000 to 2020 [57]. On the other hand, Land Surface Temperature (LST) increased by 2 • from 1990 to 2011 [58].

Materials
In the present study, to develop a HVI model, various data was used. The authors used a Landsat 8 (OLI) image (date: 18

Materials
In the present study, to develop a HVI model, various data was used. The authors used a Landsat 8 (OLI) image (date: 18 April 2021) from the United States Geological Survey (USGS) website to identify Land Use/Land Cover (LULC) and calculate LST and other indices (Spatial resolution: 30 m). The Bangladesh Bureau of Statistics (BBS) provides demographic and socioeconomic statistics. The Bangladesh Metrological Department provided metrological statistics (BMD). The authors used the National Geographical Data Center (NGDC) data on nighttime lights (NTL).

Preparation of Parameters for HVI Modeling and Their Rationale
In the present study, based on previous literature review, the authors selected twenty-six parameters in terms of physical, land use, land surface temperature, and socio-economic related parameters from different sources. Except for satellite-derived parameters, the authors obtained all parameters directly from different sources in vector format. Then, the authors used the 'vector to raster' tool of ArcGIS software to convert the vector format to a raster format. The authors made all parameters in uniform resolution. The satellite-derived indices (spatial resolution: 30 m) have been used as a baseline for resampling. Through resampling technique, the authors did uniform spatial resolution for all parameters (30 m). The rationale for choosing the parameters is provided in Table 1. The authors followed several steps to extract satellite-derived parameters.

Very young population
Young people <9 years are even more vulnerable to heatwaves than old people, caused by their lower sweating rate and body mass ratio [63].

Female
Women's responses to heat wave sensitivity differs from men due to their higher percentage of body fat (BF%) and physical strength [64]. Disabled people are sensitive to heat exposure due to their dependency on others as well as physical health and fragility [66].

Working age population
Employed people are a huge population working altogether in a city and respond to heat wave sensitivity; they are mostly exposed due to their criteria of work [67].

Household
Because of congested and increased population, the number of households are also increasing to a great extent and are exposed to heat wave sensitivity [68].

Poverty
People below the poverty line lead a miserable life with poor nutrition, poor housing conditions, and obstructive socioeconomic situations that leads to heat wave sensitivity exposure [69].
People/km 2 , Vector Water accessibility Usage of water for different purposes, especially for the tempering process, as well as treatment, refers to heat wave sensitivity [70].

Floating people
Floating people are vulnerable to heat wave sensitivity because they are homeless, they lack water access, acute illness, and electricity access; their very poor living conditions make them sensitive to heatwaves [71].

Kaccha Structure
Assessing Kacha structure of a household can indicate the sensitivity when exposed to heatwaves. Kacha structures are not properly stable, and because of the tin roof, heat affects the people miserably [72].
Number/km 2 , Vector NDBI Normalized Difference Buildup Index is an important indicator for heat wave sensitivity because the increasing build-up area, spatial growth of the buildings, and congested urban lands are simultaneously exposed to sensitivity [73].

USGS, 2021
Built-up Area Built-up area contributes to urban heat sensitivity because of its characterization of the overcrowding nature of buildings, apartments, settlements, and so on [74]. Literate people are aware of the facts and knowledge regarding heatwaves, and they prepare themselves to adapt to the situation smartly [75].

BBS, 2011
Household with electricity Access to electricity accelerates the capacity to resist heatwave disasters by using fans, air conditioners, refrigerators, and other necessary elements [76]. Number/km 2 , Vector Pucca structure Pucca structure contains proper flooring, ceiling, and structural embellishment along with socioeconomic factors, which are expressed as adaptions exposed to heatwaves [77].

Road
Larger and wider roads tend to decrease congestion from roads and land. Thus, this feature proves to be in adaptive capacity [70].
km ratio, Vector The satellite image was found to have 14.47% cloud cover and 14.36% scene cover during the pre-processing stage. The pre-processing stages include removing cloud cover using ENVI 5.3, radiometric calibration, and atmospheric adjustment.
LST was evaluated from the following perspective: retrieval of Landsat 8 reflectance (OLI) [82]. Following the extraction of Landsat 8, the NDBI (Normalized Difference Built-up Index) was computed [83]. The NDVI (Normalized Difference Vegetation Index) and NDWI (Normalized Difference Water Index) were developed using the technique of [84]. Finally, the analysis of nighttime light data was performed using the technique of [85].

Method for LULC Classification and Validation
The LULC classification was performed using the unsupervised classification technique in ERDAS Imagine software. In the first step, the classification was performed in 75 classes with no extensive prior knowledge [86]. In the next step, the 75 classes were shortened to only 5 classes (i.e., waterbody, vegetation, agriculture, bare land, build-up area) by using the recode tool in ERDAS Imagine and getting the information from Google Earth images as well as the ground. Post processing of the satellite image for LULC was performed using a majority function from the statistical filtering tool for the enhancement of the image [87]. Finally, the kappa coefficient was applied for the accuracy assessment using 150 sample points taken from the classified images as well as Google Earth imagery for the similar locations.

Heat Vulnerability Index Modeling
Extreme urban temperatures are dangerous to people's health. Heatstroke, heat fatigue, heat syncope, heat cramps, and even mortality may occur when exposed to extreme temperatures [88,89]. According to one study, heat-related mortality in the New York area might rise by 95% by 2050 if effective risk mitigation measures are not taken [90].
The vulnerability of a population to environmental risks is a function of its exposure, sensitivity, and adaptation capability (Equation (1)). Exposure refers to the intensity and geographical distribution of increased temperatures in the event of heat-related mortality. Extreme heat exposure may vary in time and space because of a changing climate and the urban heat island effect, where some parts of a city are much warmer than others. How effectively a population can deal with increased exposure or the amount to which greater exposure may harm them physically is referred to as sensitivity. For example, people with pre-existing heat-sensitive illnesses, such as diabetes, will be more sensitive to heat than those who do not have a pre-existing medical condition [91]. The climate of a city affects local sensitivity because of acclimatization. For instance, a southern city will be less susceptible to heat than a northern city since these places are more often exposed to greater temperatures because of their close location to the equator [92]. Finally, adaptive capacity refers to a population's ability to actively reduce or deal with increasing individual exposure [93]. Those having access to mechanical air conditioning, for example, have a more adaptive ability than those without, regardless of heat sensitivity. This trio of components of vulnerability has been shown to affect heat-related morbidity and death, and they may forecast heat risk.

Urban heat island vulnerability index f ramework
where, E is the exposure level, S is the sensitivity level and A is the adaptive capacity level in the equation. These variables contain different socioeconomic and index components that quantify the vulnerabilities in broader sense [93].

Modeling Process
The additive relationship of the indicator was discovered by a principal component analysis to obtain the HVI for the DMA (PCA). Various studies [94,95] adopted this strategy at municipal and national levels. The PCA approach was then used to decrease the dimensionality of these variables. Because all of our variables needed to be on the same scale in order to run PCA, the authors used min-max normalization to standardize them all. All the variables reflect the same direction, implying that as the value increased, so did the susceptibility. The Kaiser-Meyer-Olkin test performed PCA (KMO). To reduce the dispersion of loadings and enhance the understanding of principal components (PCs), the varimax (orthogonal) rotation was performed. The Jolliffe [96] criteria were used to decide how many principal components should be preserved, and PCs with eigenvalues larger than 1.0 were chosen. PC scores were created for the retained components. The scores were normalized to give them a value between the standard deviation below and above the mean of 0. Because no data existed to prove the nonlinearity of the variable and vulnerability, it was presumed that they both had equal influence in determining the DMA heat vulnerability index.

Spatial Analysis Using the Statistical Method
The study added another dimension by incorporating two powerful statistical techniques. Cluster and outlier (Anselin Local Moran's I) analysis, as well as a hotspot (Getis-Ord Gi*) analysis, were generated using spatial autocorrelation (Global Moran's I). In the Dhaka Metropolitan Area, the spatial distribution of the urban heat island vulnerability index is boosted using spatial analysis methods.
Spatial autocorrelation determines whether the expressed pattern is clustered, dispersed, or random [97]. It computes the Moran's I Index value, as well as a z-score and p-value, to assess the index's significance. The correlation coefficient limits p-values, which are statistical representations of the area under the curve for the distribution. In this study, z-value > +2.58, p-value < 0.01 was observed using a confidence level of 99%, and a positive Moran's I index was found. The result showed a visualization of clustered patterns in the heat island Dhaka Metropolitan Area.
Cluster and outlier analysis (Anselin Local Moran's I) has been used in this study to locate the clusters in our study area [98]. One significant cluster feature with a positive Moran's I indicates a similar neighboring feature with a high or low attribute value.
However, a negative value of Moran's I indicate dissimilarities with a neighboring feature referred to as an outlier.
Hotspot analysis of the study area has been calculated using Getis-Ord Gi* statistics, which provide a z-score and p-score to identify cluster variation [99]. A high z-score > 1.65 with a low p-value < 0.1 shows the hotspot area. However, the reverse scenario explains an area to be a cold spot. Table 2 shows the distribution of LULC in the study area in 2021. As found by the result, the build-up area was the dominating LULC type with 188.77 km 2 (i.e., 63.07%) in the study area. 15.11% and 10.80% of the areas were found for agriculture and vegetation, respectively. At the same time, 21.01 km 2 and 11.97 km 2 areas were found for water bodies and bare land, respectively. Figure 2 illustrates the spatial distribution of LULC types in the study area. Agriculture and bare land were in the western part of the study area. The accuracy assessment using kappa coefficient showed that the accuracy of the classified LULC map is satisfactory in this study. The accuracy assessment has been done to evaluate the performance of the classification technique as well as to examine whether the LULC classification has been done correctly or not. Studies have pointed out that an accuracy of greater than 80% is satisfactory in the LULC classification [33,40]. The user's and producer's accuracy of the classified LULC map were 90.01% and 87.36%, respectively, while the overall accuracy was 88.17% with the corresponding kappa coefficient of 0.851.

Modeling the Urban Heat Island
The spatial distribution pattern identifies the LST variability in all the thanas of the Dhaka Metropolitan Area (Figure 3). To balance the compatibility of the weather temperature of the heat island, temperatures > 27 • C were distributed as a higher range of temperature, and temperatures of 26 • C were distributed as a lower range of temperature. Analyzing the thana temperatures, the mean LST was identified as 26. the LULC classification has been done correctly or not. Studies have pointed out that an accuracy of greater than 80% is satisfactory in the LULC classification [33,40]. The user's and producer's accuracy of the classified LULC map were 90.01% and 87.36%, respectively, while the overall accuracy was 88.17% with the corresponding kappa coefficient of 0.851.

Modeling the Urban Heat Island
The spatial distribution pattern identifies the LST variability in all the thanas of the Dhaka Metropolitan Area (Figure 3). To balance the compatibility of the weather temperature of the heat island, temperatures > 27 °C were distributed as a higher range of temperature, and temperatures of 26 °C were distributed as a lower range of temperature. Analyzing the thana temperatures, the mean LST was identified as 26.53 °C with a 27.62 °C and 24.74 °C, respectively. The Tejgaon Industrial Area of Thana was identified as having the highest temperature of 27.62 °C. The Tejgaon industrial area was observed close to the core of the Dhaka Metropolitan Area, and the lands are mainly used for industrial and commercial purposes with a high range of built-up areas. However, paltan and syampur thana also had similar temperature patterns of 27.51 °C and 27.47 °C, respectively. On the contrary, the lowest temperature was observed in Uttar Khan Thana, with a value of 24.74.

Statistical Analysis of Variables of Exposure
Exposure variables were observed using the extraction method of principal component analysis (PCA) and the varimax rotation method with Kaiser Normalization. The rotation converged, followed by three iterations, and from the KMO and Bartlett's test, the Kaiser-Meyer-Olkin index was found to be 0.64. Only two principal components were chosen to be significant in the analysis (Table 3) (Table 4).

Statistical Analysis of Variables of Sensitivity
For the sensitivity variables, two methods are combined, identified as the extraction method and the varimax rotation method, to complete the analysis. The KMO index was identified with a value of 0.848. The final result of the rotation has been found with two factors where the eigenvalue > 1.

Modeling Spatial Distribution Pattern of HVI in Dhaka Metropolitan Area
The HVI index was visualized to show the spatial distribution pattern of the heat island. The authors interpreted it using three variables: exposure, sensitivity, and adaptive capacity (Figure 4a-c). The authors obtained a positive correlation between exposure and HVI and sensitivity to HVI. Exposure and sensitivity have coefficient values of 0.62 and 0.80, respectively (Figure 4a,b). However, adaptive capacity had a negative correlation with HVI with a coefficient value of −0.26. According to the index, thana in the Dhaka Metropolitan Area has been categorized into five classes: very low, low, moderate, high, and very high (Figure 4d). The authors observed Lalbagh thana with the highest HVI index value of 1. Shympur, Pallabi, Chak Bazar, Jatrabari, and Kadamtoli thana have also been identified within the very high HVI range (0.70-1.00) with high vulnerability. However, thana, which comprises very low HVI index values, are demarcated within (0.00-0.23). The Lalbagh thana was found to have a very low adaptive capacity because of population concentration, build-up area, number of households, and value of LST (Figure 4d). Again, the Gulshan area has a very high adaptive capacity and is indexed as low HVI because the area has many health institutions, a high number of literate people, adequate electricity, good roads and highways, and pucca structures. A total of 11 thanas among 41 were observed to have a moderate HVI value, allowing a stable environment for some more years. adequate electricity, good roads and highways, and pucca structures. A total of 11 thanas among 41 were observed to have a moderate HVI value, allowing a stable environment for some more years.

Thana-Based Land Use Land Cover Variation of Urban Heat Island Vulnerability Index in Dhaka Metropolitan Area
The LULC of Dhaka metropolitan has been detected with extreme built-up areas compared to other land classes (i.e., water-body, bare land, agriculture, and vegetation).

Thana-Based Land Use Land Cover Variation of Urban Heat Island Vulnerability Index in Dhaka Metropolitan Area
The LULC of Dhaka metropolitan has been detected with extreme built-up areas compared to other land classes (i.e., water-body, bare land, agriculture, and vegetation). The water-body has decreased to 25%, and built-up areas were removed from the chart below. The urban heat island vulnerability index is categorized into five classes: very high, high, moderate, low, and very low. Thus, the LULC classes were calculated according to the thana area for each category. From Figure 5a, the very high range of HVI shows only 0.60% vegetation and 0.41% of waterbodies in Chakbazar, which is highly vulnerable to heatwaves. Chakbazar, Shyampur, Jatrabari, and Kadamtali have a very low value of agriculture area. From Figure 5b, high values of HVI are concentrated in 11 thana areas of the Dhaka Metropolitan Area. The Shahbag thana area was found to have a very high level of vegetation (40.26%), and the Bangshal, Gendaria, Kamrangir char, and Shutrapur areas have the lowest vegetation recorded, facilitating the high score of HVI in these areas. In Figure 5c, the highest vegetation value was recorded in Shah Ali Thana with a 41.52% score. In this category, the land classes of each thana have comparatively similar values. Darus salam, Khilkhet, Ramna, Uttarkhan are calculated with many water bodies. Figure 5d,e show the low and very low categories of HVI where Dhanmondi and Gulshan have been identified with a very high range of vegetation. Thana areas with the highest HVI are vulnerable because they have a very low amount of agriculture and vegetation. Moreover, the very low value of bare land also indicates that the future prediction of increasing buildup areas is on a streak in the Dhaka Metropolitan Area.

Identifying Correlation of NDVI, NDBI, and NDWI with Urban Heat Island Vulnerability Index in Dhaka Metropolitan Area
NDVI, NDBI, and NDWI are three biophysical parameters that were used in this research study ( Figure 6). NDWI shows a negative correlation with HVI, having a score of −0.012, which means HVI has a positive shift if NDWI decreases (Figure 7). However, both NDVI and NDBI have a positive correlation with the heatwave index, having scores of +0.039 and +0.132, respectively (Figure 7). The NDVI range was identified with a scale of 0.07-0.28, and this vegetation index decreases the land surface temperature, which eventually turns the score of the HVI index to low in the Dhaka Metropolitan Area (Figure 6a). NDBI ranges from −0.150-0, and due to increased demand for build-up areas, the Dhaka Metropolitan Area has been congested with structures that raise the HVI score to high levels, leaving the city extremely vulnerable to heatwaves (Figure 6b). Most of the water bodies in the Dhaka Metropolitan Area are congested in specific areas. Figure 6c depicts the high values of waterbodies in the Jatrabari and Shyampur areas, which creates a positive relationship with HVI and accelerates the vulnerability rate. In general, the highly built-up regions with high LST should have a high heat vulnerability index, but the increase of adaptive capacity can reduce the magnitude of heat vulnerability. Therefore, in the present, the places with high buildup and LST also have low HVI because of the presence of high adaptive capacity. For this reason, the present study shows that the NDBI and NDVI had low correlation with HVI. Also, the spatial factors can affect the relationship. If the authors use bivariate local Moran's I, then the authors can get the actual scenarios.

Identifying Spatial Distribution Pattern of Cluster-Outlier Analysis and Hotspot Analysis in Dhaka Metropolitan Area
In this study, Moran's I index was found to have a positive value of 0.250 with a z-score of 2.901 and a p-value of 0.003, which indicates the clusters of the urban heat island to be in clustered patterns. For identification of the clustered locations and patterns where the urban heat island vulnerability has undergone a considerable change, spatial autocorrelation (Global Moran's I), cluster and outlier analysis (Anselin Local Moran's I), and hotspot analysis (Getis-Ord-G*) were introduced (Figure 8). Cluster and outlier analysis has been divided into four patterns: HH (high-high cluster), HL (high-low outlier), LH (low-high outlier), and LL (low-low) clusters. The Dhaka Metropolitan Area was illustrated with only two cluster patterns: the HH cluster and the LL cluster. HH indicates that the thana of Dhaka has higher median values that are significantly higher than their neighboring values, and LL shows vice versa. Chalk bazar, Demra, Genderia, Jatrabari, Kadamtali, Kamramgir Char, and Shyampur are seven thanas that have been identified with a higher cluster value. However, Kalabagan, Ramna, and Tejgaon have been identified as having lower values of clusters. Chalk Bazar, Jatrabari, Kadamtali, and Shyamapur thana had the highest HVI values, which justifies the higher cluster values. Jatrabari (2.634), Chalk Bazar (2.763), and Kamrangir Char (2.683) thana uncovered with the highest z-score values (of Getis-Ord-G* hotspot explaining 99% confidence level and Tejgaon (−2.413) has scored the lowest z-score value. Tejgaon thana has been detected with 95% of the cold spot value, which justifies the lowest HVI value in our study area.  highly built-up regions with high LST should have a high heat vulnerability index, but the increase of adaptive capacity can reduce the magnitude of heat vulnerability. Therefore, in the present, the places with high buildup and LST also have low HVI because of the presence of high adaptive capacity. For this reason, the present study shows that the NDBI and NDVI had low correlation with HVI. Also, the spatial factors can affect the relationship. If the authors use bivariate local Moran's I, then the authors can get the actual scenarios.

Discussion
Previous research [100,101] has highlighted the importance of heat vulnerability studies, particularly in developing nations. Early detection of high heat susceptibility may aid in healthcare and distribution of essential requirements during and after the occurrence. This study's severe heat sensitivity analysis was based on 26 characteristics that might make people more vulnerable to intense heat. The thana-level severe heat sensitivity index that results and the identification of driving elements may help researchers better understand population susceptibility during extreme heat occurrences. Rapid urbanization is happening in the Dhaka metropolitan region. The widespread growth in economic sectors fueled the development of urban areas, such as industry and agriculture. The core parts, southwestern sections, and heavily inhabited areas are the locations with the highest LST values. Due to urban heat island effects, metropolitan regions have greater land surface temperatures than rural areas [102,103]. The land use and land cover of the locations control LST distribution, wherein urban areas with less vegetation have higher land surface temperatures than rural and forest areas with lower land surface temperatures [104,105]. People living in locations with high LST have more significant, unfavorable health consequences than populations living in areas with low LST [103][104][105][106], hence the land surface temperature has been utilized as an extreme heat susceptibility indicator.
Defining high heat vulnerability, social sensitivity, and adaptive capability is crucial. The population density in urban areas is often high. Without good planning, squatters

Discussion
Previous research [100,101] has highlighted the importance of heat vulnerability studies, particularly in developing nations. Early detection of high heat susceptibility may aid in healthcare and distribution of essential requirements during and after the occurrence. This study's severe heat sensitivity analysis was based on 26 characteristics that might make people more vulnerable to intense heat. The thana-level severe heat sensitivity index that results and the identification of driving elements may help researchers better understand population susceptibility during extreme heat occurrences. Rapid urbanization is happening in the Dhaka metropolitan region. The widespread growth in economic sectors fueled the development of urban areas, such as industry and agriculture. The core parts, southwestern sections, and heavily inhabited areas are the locations with the highest LST values. Due to urban heat island effects, metropolitan regions have greater land surface temperatures than rural areas [102,103]. The land use and land cover of the locations control LST distribution, wherein urban areas with less vegetation have higher land surface temperatures than rural and forest areas with lower land surface temperatures [104,105]. People living in locations with high LST have more significant, unfavorable health consequences than populations living in areas with low LST [103][104][105][106], hence the land surface temperature has been utilized as an extreme heat susceptibility indicator.
Defining high heat vulnerability, social sensitivity, and adaptive capability is crucial. The population density in urban areas is often high. Without good planning, squatters may increase, resulting in inconvenient housing and a lack of necessary services [107]. Since 2000, an increase in squatter households has become a significant concern in Dhaka's metropolitan neighborhoods [108]. If the bulk of the vulnerable population lives in specific locations, the situation worsens. Rural locations, in certain situations, are more sensitive than urban areas owing to the large senior population. The senior population in Dhaka is most extensive in the southern and western parts since they have dwelt in the suburban area after working in metropolitan areas for an extended period. The growing cost of living diminishes purchasing power, harming low-income households who struggle to meet basic needs and have reduced access to healthcare services [109].
Critical prerequisites to overcome vulnerability determine an area's adaptation capability [110]. Most of the districts with solid adaptation capability in this research were in the state capital, which houses most governmental agencies' headquarters. The greater one's adaptive capability, the smaller one's vulnerability, and the faster one recover [111]. Considerable vulnerability is expected if the projected exposure and sensitivity are more than the adaptive capability. That explains why locations with great adaptive capacity have a high severe heat risk index.
The authors developed the "heat vulnerability index" and "vulnerability hot spots" in DMA, suggesting that metropolitan areas are more vulnerable to heat than other regions. One of the critical factors that increase a population's susceptibility is urbanization. This situation is quite similar to what has been seen in other regions of the globe [112,113] and in nations with comparable latitudes, such as Lahore, which is seeing a rise in temperature exposure [114]. The highest HVI index value of 1 was found in Lalbagh thana. Within the extremely high HVI range (0.70-1.00), Shympur, Pallabi, Chak Bazar, Jatrabari, and Kadamtoli thana have also been highly vulnerable.
On the other hand, thana is a region with extremely low HVI index values (0.00-0.23). Because of its population density, buildup area, number of homes, and LST value, the Lalbagh thana has a relatively poor adaptation potential. Because the area has many health institutions, a high number of literate people, adequate electricity, good roads and highways, and pucca structures, the Gulshan area has a very high adaptive capacity and is indexed with a low HVI. A total of 11 thanas out of 41 were found to have a moderate HVI value, showing that the environment will remain stable for a few more years.
After the generation of the HVI model with exposure, sensitivity and adaptive capacity, it needs to be identifying the most sensitive parameters for HVI occurrences. Therefore, the correlation has been performed between HVI and HVI conditioning factors. The authors found that the elderly population (r: 0.61), very young or children (r: 0.607), illiterate people (r: 0.585), females (r: 0.57), LST (0.564), working age people (r: 0.539), and poverty (r: 0.51) are variables that are very sensitive to heat vulnerability. It is evident that LST can create an urban heat island effect, which, in turn, causes higher vulnerability to heat. On the other hand, because of chronic medical conditions and health issues, elderly people over the age of 65 are particularly vulnerable to urban heat sensitivity. Young individuals under the age of nine are more sensitive to heatwaves than older adults because of their lower sweating rate and body mass ratio. Women are more sensitive to heatwaves than males due to their greater body fat percentage (BF%) and physical strength. Illiterate people are particularly vulnerable to heatwave vulnerability because they are unaware of the potential dangers of a heatwave. Employed people make up a large portion of the population in a city, and their responses to heatwave sensitivity are varied, and they are mostly exposed due to their work requirements. People living in poverty have poor nutrition, poor housing conditions, and an obstructive socioeconomic situation, all of which contribute to heatwave sensitivity. In addition, kaccha structure (r: 0.217), disabled population (r: 0.34), floating population (r: 0.37), and population density (r: 0.44) are considered as moderately sensitive to heat vulnerability. Huong et al. [71] reported that disabled persons are more vulnerable to heat exposure due to their reliance on others, as well as their physical health and fragility. Furthermore, floating people are susceptible to heatwave sensitivity due to their homelessness, lack of access to water, severe sickness, and lack of access to power, as well as their terrible living conditions [71]. It can be said that the sensitivity of a household's Kacha structure may be determined by assessing its susceptibility to heatwaves. The Kacha constructions are not fully sturdy, and the heat from the tin roof has a negative impact on the residents. In contrast, the authors found that road length (0.11), NDVI (0.06), health institution (−0.046), household with electricity (−0.012), packa structure (0.015) have inverse or very low relationship with HVI, because larger and broader roads can reduce traffic and land congestion. Therefore, this feature is adaptable. On the other hand, increasing the number of health institutions helps treat and prevent heatwave disasters. Pucca buildings include adequate flooring, ceiling, and structural ornamentation, as well as socioeconomic variables that adapt to heatwaves. Electricity facilitates the use of fans, air conditioners, refrigerators, and other cooling devices during a heatwave. Increasing vegetation minimizes urban heat island effects and preserves the environment in an ecofriendly manner. When compared to the average temperature of green space and residential areas, differences in solar radiation heating in urban areas are a primary factor contributing to the increase in air temperature and surface temperature in urban areas, particularly around roads, commercial and industrial areas, and especially around residential areas. Therefore, presence of vegetation can reduce the heat vulnerability effect.
Similar to the outcomes of other worldwide research [91], our index highlighted demographic, socioeconomic, environmental, and healthcare system variables. However, there are significant variations in selecting initial factors that make this index applicable to Bangladesh and developing nation contexts. The HVI thanas with the high and very high values were in the southern region. These thanas, which have a more significant migrant population from rural areas, have been at the bottom of many health, education, economic, and population growth metrics. At the same time, air conditioning has been shown to have the most considerable effect on decreasing heatwave mortality in the United States [115]. It is unlikely to be a viable solution for Bangladesh, particularly for Dhaka, due to the vast numbers of migrants who are designated as highly low-paid employees and do not have the financial means to acquire an air conditioning system.
As a result, appropriate local adaptation measures must be explored. Some of them have been explored in the literature, such as public messaging (radio, television), mobile phone-based text messages, automated phone calls, and amber alerts; other examples include the traditional adaption behaviors of sheltering indoors, adopting comfortable clothes, and diets. These are frequently clear in house design and materials used in building. Shaded windows and underground water storage tanks are simple architectural features that may be beneficial. Similarly, the use of insulator housing materials could be a useful preventative measure. It may be critical to access drinking water and indoor toilets within the housing complex. The authors picked several home amenities as a substitute for income and because of their protective function.

Conclusions
The present study explains the hazardous heat vulnerability of urban heat islands and the spatial distribution of the heat index at the thana level. HVI allows us to introduce the other three most important variables that played a vital role in the analysis, including exposure, sensitivity, and adaptive capacity. Because of the development in built-up regions, reduction in water bodies and vegetation, and other socioeconomic variables, the Dhaka Metropolitan Area represents most incidents of heat events. Six thanas in the Dhaka metropolitan region were discovered to have a very high HVI index, while 11 were discovered to have a high HVI rating. Alarmingly, around 87% of the overall Dhaka Metropolitan Area has been affected by heat events and is recognized as an urban heat island. The spatial distribution of adaptive ability is asymmetric. Mirpur thana has the most significant adaptive capacity value, but it also has the highest exposure value. The highest HVI values were found in Jatrabari, Kadamtali, and Shyampur. Land surface temperature increases and heat events occur because of excessive built-up regions and structural setups. A recent study has linked the rising temperature in the Dhaka City Corporation region to the Matuail dump near Jatrabari. The landfill contributes significantly to methane emissions, a key component of rising land surface temperatures. However, apart from lightening up the issue of increasing vulnerability to heatwaves, the study also has limitations. The satellite images used were taken in the afternoon, which does not explain the night temperature of the urban heat island. Furthermore, the heatwaves were not analyzed in this study using meteorological temperature data and solar radiation which are another limitation of this study. An integrated analysis of meteorological temperature and satellite temperature may provide more accurate results in the analysis of HVI. Thus, it is suggested to integrate the meteorological temperature with the satellite temperature and solar radiation for the analysis of vulnerability to heatwaves in the future studies. Despite the limitations described above, this work has a lot of positive aspects. This is the first study of heat vulnerability in the Dhaka metropolitan region. It is a preliminary screening tool for identifying where to concentrate heat-health and climate adaptation actions. This method might be used in future vulnerability assessments. Moreover, the findings of this study may provide an excellent platform for the UHI mitigation and maintaining the urban thermal comfort in Dhaka.

Data Availability Statement:
The datasets used in the study are available from the corresponding author upon reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.