Evaluation of Urban Heat Island (UHI) Using Satellite Images in Densely Populated Cities of South Asia

: Rapid Urbanization, and other anthropogenic activities, have ampliﬁed the change in land-use transition from green space to heat emission in built-up areas globally. As a result, there has been an increase in the land surface temperature (LST) causing the Urban Heat Island (UHI) effect, particularly in large cities. The UHI effect poses a serious risk to human health and well-being, magniﬁed in large developing cities with limited resources to cope with such issues. This study focuses on understanding the UHI effect in Kathmandu Valley (KV), Delhi, and Dhaka, three growing cities in South Asia. The UHI effect was evaluated by analyzing the UHI intensity of the city with respect to the surroundings. We found that the central urban area, of all three cities, experienced more heat zones compared to the peri-urban areas. The estimated average surface temperature ranged from 21.1 ◦ C in March 2014 to 32.0 ◦ C in June 2015 in KV, while Delhi and Dhaka experienced surface temperature variation from 29.7 ◦ C in June 2017 to 40.2 ◦ C in June 2019 and 23.6 ◦ C in March 2017 to 33.2 ◦ C in March 2014, respectively. Based on magnitude and variation of LST, highly built-up central KV showed heat island characteristics. In both Delhi and Dhaka, the western regions showed the UHI effect. Overall, this study ﬁnds that the UHI zones are more concentrated near the urban business centers with high population density. The results suggest that most areas in these cities have a rising LST trend and are on the verge of being UHI regions. Therefore, it is essential that further detailed assessment is conducted to understand and abate the impact of the temperature variations.


Introduction
The frequency and magnitude of extreme weather events, such as heatwaves, are expected to rise with an increase in air temperature [1]. Such events are exacerbated when coupled with Urban Heat Island (UHI) effect. UHI is the phenomenon where urban air temperatures are higher than the surrounding rural areas [2]. Several factors, such as an increase in anthropogenic heat flux's emission [3], change in urban geometry, and population density [4], and change in land-use and land cover LULC [5], results in the UHI phenomenon. With the rapid increase in urbanization, the green land cover is replaced by impervious land surfaces, such as concrete buildings and bituminous roads [6,7]. Change in land cover properties alters the thermal properties, surface radiation, and humidity of the urban area [8], leading to the UHI effect. Evaluation of UHI in urbanized and populated cities is crucial to analyze the change in surface albedo, emissivity, and evapotranspiration [6]. The UHI phenomenon has been widely studied [9][10][11] since its first observation by

Landsat-8 Data
We used Landsat 8 satellite images, OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) 15 to 30-m multispectral data from Landsat-8 C1 Level-1, downloaded from United States Geological Survey (USGS). Thermal band, Band 10, is provided as the atmospheric brightness temperature in Kelvin (K), and the multispectral-bands of Landsat-8 OLI is provided as surface reflectance. Band 10 is available at 100m horizontal resolution. Since the UHI effect is weakened by cloud cover [36], the satellite imagery with cloud coverage <10% were selected (cloud-free). Cloud cover in the entire study area is not similar. Data was selected to maintain temporal uniformity over the study period (Table 1). This data selection should not be hindered as our study focused on LST variations and UHI over modified land use rather than the absolute value of LST. A comprehensive outline of the necessary metadata is outlined in Table 1. Data required were downloaded from Earth Explorer [37] and masked to the study area. Then, the LST map for each study area was generated for different data acquisition dates as shown in Table 1. The extracted images for most of the years were obtained for summer (March-August) so that the high surface temperatures could be captured. Since we filtered the cloud cover to be less than 10%, the downloaded satellite image could not coincide for the same day in each year and most of the images were acquired for March. The scaling-factors used for determining LST and UHI were extracted from metadata (K 1 = 774.4483, K 2 = 1321.0789, M L = 0.0003342, A L = 0.1, please refer to Section 2.3.1 for details).

Population Distribution
The annual population distribution in the KV, Delhi and Dhaka is shown in Table 2 [38]. The increase in average rates of the population for KV, Delhi and Dhaka were found to be 3.97%, 3.25%, and 3.61% respectively, the highest being in KV. The rate of population change increased until 2014, then decreased until 2017, and again increased in the year 2018 for KV whilst remained constant throughout the study period for Delhi and Dhaka except in the year 2019. The rate of population change is expected to drop to 3.49% for KV, 3.03% for Delhi and 3.56% for Dhaka in the year 2020. Table 2 shows the increasing population trend for all study areas which could eventually impact the land use pattern and the surface temperature. The rapid increase in the urban population in two decades has led to the rise in the rate of urbanization too in these cities.

Observed Air Temperature
The average observed air temperature acquired [39] for all the study areas is shown in Table 3. The observed temperature for the missing satellite data acquisition months in 2015 and 2018 for Delhi were obtained on 28 July 2015 and 20 June 2018 respectively. Likewise, the observed air temperature for Dhaka in 2013 was obtained on the day of 20 March 2013. Meanwhile, the missing temperature data on the satellite data acquisition day were obtained by averaging the temperature from the previous day and the following day of the missing date.

Study Area
South Asian (SA) nations, namely: Afghanistan, Pakistan, India, Nepal, Maldives, Bhutan, Sri Lanka and Bangladesh combined have more than 18 billion population [40] where 31% resides in the urban area [41]. The capitals of SA nations are economically rich, urbanized, and highly populated. Delhi (30.2 million) and Dhaka (21.0 million) are the top 10 most populated capital cities in SA, while KV is the most populated and emerging urbanized city in Nepal. The rising population has induced urbanization in all these capitals. The rate of urbanization is increasing rapidly with the increase in population to 130 million in just a span of 10 years (2001 to 2011) and is expected to reach 250 million by 2025 [42] in SA. A brief description of each study area is discussed below under each sub-sections. The section deals with the physiography and climate of each study area.

Kathmandu Valley
KV is one of the biggest cities, in terms of population and economic development, in Nepal comprising an area of approximately 664 km 2 ( Figure 1). The bowl-shaped valley is located at 85 • 11 E 27 • 32 N inscribing the major cities Kathmandu, Lalitpur and Bhaktapur [43,44]. KV lies in the warm temperate zone [45] with a fair climate varying from 2.4 • C to 37 • C in the period of 1981-2010 [46]. Geographically, the central lower part of the valley is situated at an elevation of 1425 m above mean sea level MSL and is surrounded by four mountain ranges namely Shivapuri, Phulchowki, Nagarjun, and Chandragiri Hills. KV is the most developed city in Nepal with a greater portion occupied by the built-up area. On the contrary, KV has some open spaces namely Tundikhel, Tribhuvan park, Sankha park, etc., which are expected to serve as green space thereby reducing the rising daytime temperature.

Delhi
Delhi located in Northern India, at 77 • 14 E 28 • 36 N is bordered by Uttar Pradesh in the East and by the states of Haryana on the West, North and South ( Figure 1). It covers an area of 1484 km 2 , of which 783 km 2 is designated as rural, and 700 km 2 as urban therefore making it the largest city in terms of area in India. It has a length of 51.9 km and a width of 48.48 km. Delhi has a dry subtropical and semi-arid type of climate [47]. The annual average temperature is 25 • C with monthly means varying from 13 • C to 32 • C. The warm season begins in early April and peaks in late May or early June with an average temperature of about 38 • C but occasional heat waves may result in high temperatures of over 45 • C on some days. Per capita availability of green space in Delhi is about 20 m 2 whereas that of open space is 30 m 2 [48].

Dhaka
Dhaka is located in central Bangladesh at 90 • 22 E 23 • 42 N (Figure 1) along the eastern banks of the Buriganga River. The city lies on the lower reaches of the Ganges Delta and covers a total area of 306.38 km 2 . Dhaka has a tropical savanna climate. The city has a distinct monsoonal season, with an annual average temperature of 26 • C and monthly means varying between 19 • C in January and 29 • C in May [49]. There are many parks within Dhaka city, including Ramna Park, Suhrawardy Udyan, Shishu Park, National Botanical Garden, Baldha Garden, Chandrima Uddan, Gulshan Park, and Dhaka Zoo. There are lakes within the city, such as Crescent lake, Dhanmondi lake, Baridhara-Gulshan lake, Banani lake, Uttara lake, and Hatirjheel-Begunbari lake. These parks and lakes act as green space and help in minimizing the rising surface temperature. Figure 2 shows the flowchart of methodologies adopted in this study. A step-wise description is discussed below:

Methodology
Step-Wise Methodologies to Estimate UHI Conversion of Spectral Radiance to Top of Atmospheric Brightness Temperature Landsat 8 satellite imagery obtained was processed using 32-bit floating-point calculations. These values were then converted to 16-bit integer values in the finished level 1 product. Conversion to spectral radiance was done using the radiance scaling factors provided in the metadata file [37] using Equation (1).
The calculated spectral radiance was then converted to brightness temperature which was often determined as the effective temperature under unit emissivity. Top of Atmosphere (TOA) Brightness Temperature was obtained from spectral radiance in degree Celsius using Equation (2).
where, L λ = TOA spectral reflectance (watts/(m 2 * sr * µm)), M L = Band specific multiplicative rescaling factor, Q cal = Quantized and calibrated standard product pixel values (DN), A L = Band specific additive rescaling factor, Urban Heat Island (UHI) Step 1 Step 2 Step 3 Step 4 Landsat-8 Bands Landsat Surface Temperature Generation OLI and TIRS band data were converted to radiance using the radiance scaling factors provided in the metadata file using Equation (1) and then to TOA brightness temperature using Equation (2). After acquiring brightness temperature, LST for Landsat-8 was estimated using Equations (3)- (7). Parameters such as NDV I, proportion of vegetation (P v ), and ground emissivity (e) were required to estimate the LST. Each parameters NDV I, P v and e were calculated using Equations (4), (6), and (7) respectively. NDV I was calculated using Near InfraRed Band (B 5 ) and Red Band (B 4 ) of Landsat imagery. Normalized Difference Built-up Index (NDBI) was calculated using 6th (Middle InfraRed) and 5th (Near InfraRed) of Landsat 8 bands using Equation (5). where, , NDV I = Normalized Difference Vegetation Index calculated using Equation (4), NDBI = Normalized Difference Built-up Index calculated using Equation (5), P v = Proportion of vegetation calculated using Equation (6), e = Ground Emissivity calculated from Equation (7), λ and C 2 are constants having values of 10.8 and 14,388 respectively.

Urban Heat Island
After obtaining LST, Equations (8) and (9) were used for determining UHI used by Kaplan et al. [34]. Here, µ is the mean LST value for the study area, and σ is the standard deviation of the LST. UH

Kathmandu Valley
The estimated spatial distribution of LST for KV is shown in Figures 3 and 4 for the dates reported in Table 1. We observed an increase in surface temperature up to more than 30 • C for 2015 and 2018. In the year 2015, LST was estimated for June so most of the region in KV had higher temperature compared to other years. Kathmandu and Lalitpur Metropolitan City (referred to as Kathmandu and Lalitpur) showed higher surface temperature compared to other nearby municipalities (Madhyapur Thimi, Kirtipur etc.) and rural municipalities (Konjyosom) (refer Figure 1 for the locations of administrative units). It may be noted that in June 2015, Konjyosom rural municipality exhibited higher temperature than in other years.
The peri-urban municipalities like Nagarjun, and the higher region of Budhanilkantha, and Godawari had the temperatures ranging from 15 • C to 25 • C most of the years. These regions also had surface temperature higher than 27 • C during 1 June 2015, 22 April 2018, and 8 May 2018. This might be due to the rise in air temperature as a result of reflectance in the urban and barren land. Kathmandu and Lalitpur experienced higher temperature most of the summer time. However, the boundaries of KV such as Chandragiri, Dakshinkali, Nagarjun, Gokarneshwor, Budhanilkantha, and Konjyosom experienced temperatures less than 10 • C.
The second and third columns of Figures 3 and 4 show the temporal and spatial distribution of NDVI and NDBI for KV. The figures depicted that vegetative area (NDVI > 0) was higher in the peri-urban area compared to central KV which had NDVI < 0, negative NDVI for March 2013 and 2014. NDVI values for central KV were negative which support the presence of non-vegetative areas usually barren land or built-up area. It was observed, as expected, that the built-up area increased with the increase in population. Further, the result of NDBI confirmed that the built-up area (NDBI < 0) was concentrated in the central KV, mainly in Kathmandu, compared to peri-urban areas of KV such as Godawari and Nagarjun municipalities. Central KV, including Kathmandu and Lalitpur, had lower values of NDVI and higher values of NDBI which states that it KV was less vegetative compared to other regions. NDVI ranged from −0.  [50,51] which eventually reduces the soil moisture.

NDVI
The estimated average surface temperature ranged from 20 • C in 2016 to 32 • C in 2015 ( Table 4). The higher discrepancy in observed and estimated temperature occurred in the year 2017 and 2019 resulting in higher differences of 5.04 and 5.3 • C, respectively. While the discrepancies in other years were found to be comparatively less as shown in Table 4.   Table 5 showed that the maximum difference in the temperature was 9.3 • C in the year 2017, whereas the least bias of 2.3 • C found in the year 2019.
LST in Delhi was found to be inconsistent intra-annually. During the year 2013, June was considered for the evaluation purpose and the mean surface temperature was found to be 32. The result of NDVI portrayed that vegetative area was less compared to non-vegetative areas. The non-vegetative area increased with the increase in LST. The vegetative area was more in south Delhi compared to north Delhi which holds more built-up towns and industrial areas. However, in, 2016, NDVI map showed that Delhi has a more vegetative area compared to non-vegetative which might be the result of lesser air temperature (observed 33.5 • C) and higher relative humidity (RH) compared to 2015 and 2017 for the same month. RH was found to be 21% in the data acquisition day in 2019.  Similarly, 40 • C observed temperature and 30% RH, was observed in 2013. NDVI values ranged from −0.09 in 2013 to 0.52 in 2016 (August). This suggested that Delhi has a lesser vegetative area compared to the non-vegetative area. A similar pattern was observed from NDBI. Most regions in Delhi have NDBI value ranging from −0.3 to 0.51, the distribution shows higher percentage of the built-up area than the non-built-up area.
The NDBI values suggest that more areas in Delhi radiate the incoming sunlight thereby increasing the temperature of the study area.

Dhaka
The spatial and temporal distribution of LST, NDVI, and NDBI were plotted for Dhaka district as shown in Figures 7 and 8. The spatial distribution of LST over Dhaka showed the temperature distribution in March 2014 was higher than 35 • C, a few areas have a temperature less than 35 • C. This was validated by the observed temperature as shown in Table 6. The bias in observed air temperature and estimated surface temperature was found to be −1.23 • C in the year 2014. The maximum bias of +5.5 • C was found in 2018.
The spatial distribution of NDVI showed that most of the western districts in Dhaka such as Dhamari, Savar, Nawabgang and Dohar are vegetative compared to eastern districts such as Tejgaon and Keraniganj. It further showed that during the years 2016 and 2017, the non-vegetative area has increased in the eastern districts. The similar pattern was observed from the NDBI too. The eastern districts are highly built-up in comparison with the western. The NDBI values ranged from −0.10 in the year 2019 to +0.55 in the year 2018. Negative values of NDBI showed a non-built-up region while the positive value reflected the built-up area of the study area.

Relationship between LST, NDVI and NDBI
The relationship between LST, NDVI and NDBI is necessary to understand the UHI phenomenon in urban cities. Increase in the LST values is governed by the climatic (rainfall, air temperature, humidity) and non-climatic (land use and land cover, aerosols, air particles in the atmosphere, built-up area) variables. Spatial average values of LST, NDVI, and NDBI were used to understand the relationship between them.

Kathmandu Valley
The relationship between LST and NDVI showed a mixed relation trend for KV. The negative linear trend exists for temperature less than 25 • C and the parabolic relationship for higher values with p-value = 0.08 (p-value > 0.05) and R 2 = 0.05. The statistical performance of p-value shows that there is no statistically significant relation between mean LST and mean NDVI. The non-vegetative area increased with increment in the surface temperature (Figure 9). This implies that likely the increase in non-vegetative areas is directly associated with the change in surface temperature and built-up area [26]. For the surface temperature ranging from 22.5 • C to 30.0 • C, the NDVI increased. This indicates an increase in greenness together with an increase in surface temperature. This variation might be attributed to rainfall days before the measurement was taken.

Delhi
The relationship between mean LST and mean NDVI was established for the study period 2013-2019 as shown in Figure 10. The statistical performance of the relationship among mean LST and mean NDVI shows a negative linear trend with R 2 value of 0.79 and p-value = 0.043 (p-value < 0.05). This shows a statistically significant relationship between mean LST and mean NDVI. Since Delhi is often viewed as a single administrative unit, we have used only 5 temporal data points to develop the relationship. The Figure 10 shows the well-distributed relationship between mean LST and NDVI (R 2 = 0.79). The increase in mean LST was concomitant in the decrease in vegetative area for Delhi. Thus with the increase in the surface temperature, the number of UHI regions are likely to increase. This situation might exacerbate the risk of heat-stroke in Delhi.

Dhaka
Relationships between mean LST and mean NDVI were studied for Dhaka district and found that NDVI value increased slightly with the increase in LST. The statistical result shows very low correlation exists between mean LST and mean NDVI (R 2 ∼ 0). The p-value is equal to 0.75 (p > 0.05) is not statistically significant and cannot provide enough evidence for rejecting the null hypothesis of similarity. During data acquisition day in 2014-2019, most of the administrative units had mean LST ranging from 24.0 • C to 28.0 • C with the variation of NDVI ranging from 0.0 to 0.35 ( Figure 11).

UHI for the Study Area
UHI was calculated based on mean and standard deviation (SD) of LST. Green color shows the non-UHIs while red color shows the UHI regions for each of the study areas. The threshold value of UHI ranged from 21.6 • C to 24.2 • C for KV, from 26.2 • C to 41.3 • C for Delhi and from 24.6 • C to 34.9 • C for Dhaka (Figures 12-14 respectively). The variation in the UHIs was governed by the LST each land cover type had in the study area. Further, the availability of green space (as seen from NDVI) and clustering of the built-up area (as seen from NDBI) have a significant impact on the UHI for each study area.

Kathmandu Valley
The spatial and temporal distribution of UHI for the KV was developed for each year using the surface temperature for the same month ( Figure 12). Distribution of UHI showed that for all the years, central KV experienced higher temperature and thus may be considered as heat islands which were further supported by LST, NDVI and NDBI in earlier sections. Compared to southern regions, the Northern region is highly populated (as can be seen from NDBI maps) and thus these regions experienced higher temperatures. On 26 March 2013, the central and Northern region experienced temperatures higher than 28.3 • C. The boundary administrative such as Budhanilkantha, Gorkarneshwor, Changunarayan, Godavari, Chandragiri, etc. had comparatively lesser UHI effects throughout the study period. The result of UHI analysis also show an increasing trend in UHI especially in the peri-urban areas. Spatial analysis revealed that the UHI zones were more concentrated in the central and northern regions of the KV. UHI at the different administrative units of KV reflected the increasing trend of UHI zones which might be the impact of increasing population and intensified urbanization. Developing the new greener space might help in reducing the impact of increasing UHI in the densely populated urban cities like KV [52].  Figure 13 showed the spatio-temporal distribution of UHI for Delhi. The minimum threshold of UHI for Delhi was found to be 26.2 • C in the year 2016 on the data acquisition date. The maximum threshold of UHI was estimated to be 41.3 • C for the year 2019. With the higher values of mean LST and the least value of SD, most of the areas in Delhi are estimated to behave non-UHI. Since the UHI threshold is less and the temperature range is small, more UHI areas were observed in Delhi in 2016 compared to other high UHI threshold years. This point to a limitation of UHI threshold for cities like Delhi. Residential or mixed regions have a higher potential risk of UHI in Delhi as portrayed by Mohan et al. [53]. Researchers [53] further stressed that UHI in summer is expected to increase and become more dominant in the densely urbanized built-up areas.

Dhaka
The minimum and maximum threshold of UHI for Dhaka were found to be 24.6 • C and 34.9 • C respectively ( Figure 14). The result showed that western districts such as Dohar, Nawabganj, Dhamrai behave as non-UHI while the eastern districts such as Tejgaon, Savar, and Kernaiganj behave as UHI zones. The North-West region of Dhaka was less impacted by the rising surface temperature and is thus identified as non-UHI regions throughout the temporal study period. Eastern parts of Tejgaon and Keraniganj were identified as UHI regions for each year is associated with a higher density of impervious structures such as concrete buildings and paved roads in the region. The potential impact of UHI is increased by the rising surface temperature in the Dhaka city from 28.5 • C in 2002 to 40.1 • C in 2014 [54].

Discussions
Satellite imageries provide an appropriate platform to evaluate LST and UHI at any spatial and temporal scale. The primary objective of the study is to evaluate the UHIs in densely populated cities of South Asia namely, Kathmandu, Delhi, and Dhaka using satellite imageries. Also, we examined the LST, NDVI, and NDBI for the same to observe the state of surface temperature, wetness, and dryness of the land and Built-up intensity for the temporal period of 2013-2019.

Relationship between LST, NDVI and NDBI
Visual analysis of the LST showed that June 2015 was the most hottest month during the temporal study period in KV. In 2014 and 2016, the month of March was found to have less warm days in the high elevated regions where the LST is <10 • C (Figures 3 and 4). In the case of Delhi, high LST was observed in 2019 followed by 2013 and 2014. Similar trends of the LST variation were observed in Dhaka. NDVI values for KV and Dhaka were found to increase at the rate of 0.007 and 0.004 respectively while a negative trend was observed in Delhi ( Figure 15). The negative trend in Delhi might be the result of the presence of a relatively higher amount of the open spaces as compared to KV and Dhaka. The results depict an increasing trend of NDBI values for Delhi; almost no trend for Dhaka and decreasing trend for KV ( Figure 15). The spatio-temporal variation in the NDVI and NDBI values might be the consequences of the local climatic conditions [55] of the study area. The spatio-temporal variation over the different topographical regions has been well established in the study domains. Climates of the study area have an important role in governing LST, NDVI, NDBI, and UHI as computed in this study. The large orographic differences over a short latitude change could be responsible for lesser LST and higher NDVI values in the KV [45,56,57]. The differences in spatio-temporal LST retrieval in the study domain might have been affected by the biophysical effects, evapotranspiration, and albedo that are eventually influenced by precipitation and local climate [58]. Densely populated zones in the study areas are found to have higher LST values compared to surrounding areas. The higher values of LST in the central zone of the study area is likely due to the densely built-up area and paved roads [26,59]. Pan et al. [60] found that the LST values at built-up areas are higher than 40 • C in the humid subtropical climate. The variation in LST was also attenuated by the change in elevation. Increase in the LST values was also concomitant with increased population in all the cities of South Asia [61].   [62] and increased population [63] in each study area were found to have higher UHI values compared to surrounding areas. The increased built-up areas and the paved roads might be the driving factors that alter the spatio-temporal alteration of UHI zones. Further, the reduced open spaces (green areas) and current development works such as the construction of roads, buildings might have aggravated the increase in the UHI zones. Growth and development activities increases the impervious surfaces, resulting in reduced evapotranspiration and lesser soil moisture [64], which ultimately have a direct impact on the LST of the urban areas. The spatio-temporal variation in UHI values across the study area has been impacted by increased NDBI index [34]. UHI magnitudes increased across the regions with increased NDBI and decreased NDVI. El Niño might also have contributed to the wider variability of UHI effects in the study regions [65].

Impact of UHIs and Mitigation Strategies to Minimize Rising Surface Temperature
The UHIs has diverse impacts on different elements of the society such as energy consumption, human health, biodiversity, agriculture, water availability and others [66]. With the increase in the urban population (Table 2) and urbanization, the intensity and frequency of the heatwaves have increased. Further, an increase in the frequency and magnitude of hotter days attributed by the rise in LST and UHI zones have a direct impact on the energy consumption. The urban population tends to consume more electricity to make themselves more comfortable against the increasing heat. In recent years, the trend of electricity consumption in SA region is increasing [67,68]. The increase in electricity consumption might be the cumulative impact of the increasing population, rising urbanization, growing wealth, and climate extremities. Further, the increased climate extremes (such as heat strokes, UHI, increased LST) has also impacted the health of the people in SA. The comfort level induced by the climatic extremes in the health of the residents of any city is measured in terms of discomfort index (DI) and physiological equivalent temperature (PET) index. The previous research in SA and West Bengal (India) showed that the area with shades due to high rise buildings and trees have comfortable conditions compared to the one with no shades [66,69]. The increase in surface temperature has increased the number of heat strokes in urban cities of SA [70]. The number of patients suffering from heat stroke was found to be comparatively higher in the urban centers than in the peri-urban areas of SA. This supports the idea that increasing UHI and LST has increased the risk of heatwave globally and regionally [71,72]. The increased surface temperature has a significant impact on diarrheal disease and heat stress in Bangladesh [73]. The policy makers and planners of each study region now need to focus on the proper mitigation and adaptation strategies to cope against the rising LST and increasing UHIs. This study shows that more planning and perhaps enforcement are required to reduce the impact of the rising surface temperature. Increase in the green land area and afforestation activities can reduce the impact of excess heating from the solar radiation. This also enhances the asthetics of the city. Few mitigation strategies have been considered by the local government in Kathmandu valley such as cleaning of the Bagmati river corridor and increasing the number of new recreational parks. In India, Niti Aayog has proposed to ban the diesel vehicle and sell the electric vehicle by 2030 to reduce the air pollution. Such major mitigation measures are necessary for reduction of UHI too for the region to be better prepared for climate extremities.

Limitations of the Research
The current study has limitations in the spatial and temporal domains. Temporally, we only focused on the summer days to assess the summer surface temperature and subsequent heat island. Spatially we focused on the highly urbanized and rapidly rising cities at SA. The study focused on three major cities, however several smaller cities in SA may face similar problems. A coordinated effort is needed to understand the regional LST and UHI in detail. Also, understanding of antecedent conditions coupled with on ground sensing may be beneficial for future work.

Conclusions
This study evaluated the UHI using the Landsat-8 satellite images in three densely populated cities of South Asia namely, Kathmandu Valley, Delhi, and Dhaka. These are the growing cities in terms of economic development, urbanization and population rise. The spatial and temporal variations of LST, NDVI, NDBI, and UHI were analyzed in these three cities to assess the impact of urbanization on the surface temperatures. An increase in the impervious areas, such as concrete buildings and paved roads reduces the recharge capacity of the soils thereby reducing the soil moisture. This leads to an increase in barren and non-vegetative lands in these urban cities. The analysis of LST and UHI in these urban cities demonstrate the importance of urban planning to mitigate the effects of future climate. The analysis was focused from 2013-2019 when the cloud cover is less than 10%. The following conclusions are drawn from the research:

1.
Results of LST showed the surface temperature is more in the Kathmandu and Lalitpur Metropolitan City while the regions that are situated at the boundaries of KV experienced LST less than 10 • C, below that of KV. Similarly, for Delhi higher LST is observed in the western region of Delhi and the eastern region of Dhaka. It can be inferred that the zones which are densely populated experience higher LST.

2.
Like LST, NDVI of the study area shows more vegetative regions in peri-urban areas and less in the central KV. NDVI result of Delhi shows a lesser vegetative area than the non-vegetative area. In the case of Dhaka, the eastern district (Tejgaon) has less vegetative area compared to the western region. 3.
NDBI shows the concentration of built-up areas in most regions of the KV (NDBI > 0). Regions with NDBI < 0 are concentrated in the peripheries of KV. Similarly, NDBI result exhibits that the built-up area is concentrated in the western region for Delhi and the eastern region for Dhaka.

4.
The results of this study imply that the spatial distribution of LST magnitude and UHI zones are greater in Delhi and Dhaka compared to KV. However, the core center of the KV has a higher rate of LST magnitude and UHI effects are increasing faster annually.

5.
The results of the research provide insights into urban microclimates and changes in the environment that may be used for drafting the city planning legislation to mitigate the rising LST.

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

Abbreviations
The following abbreviations are used in this manuscript: