Assessing the Effects of Land-Use Types in Surface Urban Heat Islands for Developing Comfortable Living in Hanoi City

Hanoi City of Vietnam changes quickly, especially after its state implemented its Master Plan 2030 for the city’s sustainable development in 2011. Then, a number of environmental issues are brought up in response to the master plan’s implementation. Among the issues, the Urban Heat Island (UHI) effect that tends to cause negative impacts on people’s heath becomes one major problem for exploitation to seek for mitigation solutions. In this paper, we investigate the land surface thermal signatures among different land-use types in Hanoi. The surface UHI (SUHI) that characterizes the consequences of the UHI effect is also studied and quantified. Note that our SUHI is defined as the magnitude of temperature differentials between any two land-use types (a more general way than that typically proposed in the literature), including urban and suburban. Relationships between main land-use types in terms of composition, percentage coverage, surface temperature, and SUHI in inner Hanoi in the recent two years 2016 and 2017, were proposed and examined. High correlations were found between the percentage coverage of the land-use types and the land surface temperature (LST). Then, a regression model for estimating the intensity of SUHI from the Landsat 8 imagery was derived, through analyzing the correlation between land-use composition and LST for the year 2017. The model was validated successfully for the prediction of the SUHI for another hot day in 2016. For example, the transformation of a chosen area of 161 ha (1.61 km2) from vegetation to built-up between two years, 2016 and 2017, can result in enhanced thermal contrast by 3.3 ◦C. The function of the vegetation to lower the LST in a hot environment is evident. The results of this study suggest that the newly developed model provides an opportunity for urban planners and designers to develop measures for adjusting the LST, and for mitigating the consequent effects of UHIs by managing the land use composition and percentage coverage of the individual land-use type.


Introduction
In the context of climate change and global warming, it is important to monitor the signatures of urban heat islands (UHIs) and to understand their impacts on ecosystems and human health.UHI arises from the phenomenon of relatively higher temperature in the urban center over its surrounding rural environment.The phenomenon behind UHI has been studied for a long time.It was first described by Luke Howard in the 1810s [1].From 1964 to 1968, Bornstein [2] used a helicopter to study the UHI of New York City, and determined the effect of UHI in both vertical and horizontal directions.The results display a maximum intensity of UHI near the ground surface, and a decrease to zero at a height of 300 m [2].Ackerman [3] studied the diurnal and seasonal variations of UHI in Chicago, recording an increase-averaging 1.85 • C-in temperature inside the city most of the time.Since the 1990s, 3D models have been developed to examine the effects of UHI in Tokyo, Japan, using satellite and land survey data [4].In Nagoya, Japan, seasonal changes pertaining to UHI were analyzed using Landsat and ASTER images taken during the day, as well as at night, which were modeled to determine whether or not the heat fluxes are natural or artificial [5].In Washington, the surface temperature of the city center was found to be higher than the surrounding vegetative areas by up to 10 • C [6].These studies indicated that vegetation cover plays a key role in minimizing the UHI effect.The minimization of the UHI effect tends to be beneficial to the community because it may result in the enhancement of dangerous natural phenomena, in addition to its impacts on ecosystems and human health.For example, it has been reported that the UHI effect may alter the precipitation [7], characteristics of cloud-to-ground lightning activities through increased aerosols [8,9] and their enhancement [10,11] in response to urbanization, and modify the environmental and regional climate by reshaping the boundary layer and land-sea circulation [12].
Many other studies also demonstrated correlations between UHI and land-use composition in a city.Weng et al. [13] reported a correlation between the surface temperature and vegetation in Indianapolis, USA.Chen et al. [14] examined the relationship between UHI and land-use change, in certain cities of the Guangdong Province in Southeast China, using Landsat images from 1990 to 2000.A similar study on the relationship between UHI, land-use change, and population density was also conducted in Nagpur, India [15].Other studies relevant to UHI effects were also conducted [16][17][18][19][20].These studies provide similar conclusions that urban temperatures are highly correlated with land-use composition (water, vegetation, built-up, among others) in the cities.Due to easy access and wall-to-wall continuous coverage, LST derived from thermal infrared remote sensors are one of the most commonly used indicators for surface UHI (SUHI) analysis [21][22][23][24].In this study, the SUHI is defined as the magnitude of the temperature differentials between any two land-use types, a more general way than that which is typically adopted in the literature.From the physical point of view, LST and air temperature are different entities, while strong correlations were found between them by many researchers in the literature [21,25,26].
More recently, Deilami et al. [27] provided a systematic and overarching review of different spatiotemporal factors that affect the UHI effect.It is indicated that the UHI effect can be considered as a critical factor contributing to heat-related mortalities, and unpredictable climatic changes.Lai et al. [28] were concerned with the quality control of the satellite data for investigating the SUHI.They used eighty-six major cities across mainland China, and analyzed SUHI intensity (SUHII) derived from Moderate Resolution Imaging Spectroradiometer (MODIS) LST data.Their findings suggested the need to be extremely cautious when using LST product-based SUHIIs to interpret SUHIs.Li et al. [29] presented a new method to quantify the SUHI.They were concerned with the effective evaluation of potential heat risk.A new approach was proposed to quantify the SUHII by using the relationship between MODIS LST and impervious surface areas (ISA).The calculated SUHII shows high values in summer and during the day than in winter and at night.Despite a great effort being devoted to study the UHI and SUHI in the major cities of the world, fewer investigations have been conducted to solve the associated problems for developing countries such as Southeast Asian countries, among which Vietnam implemented a Master Plan 2030 for its capital, Hanoi, in 2011.How will the Master Plan impact Hanoi's overall environmental conditions?Nam et al. [30] evaluated the influence of UHI under the Hanoi Master Plan 2030 on the energy consumption for space cooling in residential buildings.They found that the increments in built-up areas were larger than those in existing built-up areas, and that the cooling load in an apartment is approximately half of that in a detached house, which itself approximately half of that in a row house.It was also observed that although sensible cooling loads increased with the increase in outdoor temperature, the latent cooling loads decreased due to the decrease in absolute humidity and the increase in air temperature.Trihamdani et al. [31] assessed the UHI effects in the city under the present land use conditions, as well as those conditions proposed by the Hanoi Master Plan 2030 through numerical simulation, using Weather Research and Forecasting (WRF).They found that the peak air temperature in the built-up areas (approximately 1 • C higher at the maximum) was not significantly modified, but high temperature areas, with temperatures of 40-41 • C, would expand widely over the new built-up areas.They also stressed that the number of hotspots increased further when the strategic green spaces in the master plan were not taken into account.
Based on the literature survey, state-of-the-art of correlations between UHI or SUHI and land-use composition are being widely analyzed in the recent years.In contrast, the study of UHI or SUHI with respect to Hanoi City Master Plan 2030 is rather limited; for example, about spatial energy consumption [30], assessment of UHI effects based on WRF simulations [31].Then, what would be the most updated status of land use change with respect to the city Master Plan?What would be the impact on the thermal signatures?What possible measures can the city take to mitigate UHI?Therefore, the objectives of this research are to (i) assess the land use changes in Hanoi, (ii) assess the quantitative relationships between the composition of the main land-use types and SUHI in Hanoi, (iii) analyze the effects of land-use composition on SUHI on an extremely hot day, (iv) derive a regression model for the prediction of SUHI, and (v) suggest the measures applicable for minimizing the SUHI impacts on human health, due to increased urban temperature.

Study Area
Inner Hanoi was chosen as our study site.It is situated in the capital city of Vietnam, encompassing an area of 160 km 2 (yellow polygon in Figure 1c).The studied location mainly comprises built-up areas, vegetation covers, and water bodies.The geographical location of the study area from a global view, to the national scale, to the regional scale is illustrated in Figure 1.The figure also depicts the boundary of Hanoi City before and after 2008, a composited Landsat 8 satellite image acquired on 4 June 2017, and locations of the meteorological stations Lang (yellow) and Ha Dong (orange) were used to acquire air temperature data.
In terms of population, Hanoi is inhabited by more than 7.7 million people (whole Hanoi as shown in Figure 1b) and it has a high population density (12,340 people/km 2 ) in urban districts in 2015.As for climate, it is located in a tropical belt so that its summer lasts from May to August with a hot and humid atmosphere, abundant rainfall, and an average temperature of 29 • C. On average, July is the hottest month.Its annual average rainfall is 1680 mm, and the average temperature 23.6 • C.

Air Temperature Data
The air temperature (°C) of the monthly hottest days observed at two meteorological stations (Lang and Ha Dong as shown in Figure 1) in Hanoi in 2016 and 2017 are listed in Table 1.Therefore, the days with the hottest air temperature are selected and considered as having the most negative effects on human health that are associated with UHI effect.From Table 1, it shows that the monthly hottest days are 2 June and 4 June in 2016 and 2017, respectively.

Acquisition and Pre-Processing of Satellite Data
We acquired the Landsat 8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) images for the same day (4 June 2017; scene ID: LC08_L1TP_127045_20170604_20170615_01_T1), when Hanoi experienced the greatest heat wave in recent 40 years.On that day, the air temperature in Hanoi reached up to 42 °C.Since Landsat 8 images were not available for 2 June 2016, the hottest day of 2016, we acquired the Landsat 8 OLI and TIRS images for the second hottest day of 2016, June 1 (scene ID: LC08_L1TP_127045_20160601_20170324_01_T1).For simplicity, 1 June is still named as the hottest day in 2016 hereafter, in the paper.On that day, the air temperature in Hanoi reached up to 39 °C.The Landsat 8 image of 2016 are later used to validate the models that are

Air Temperature Data
The air temperature ( • C) of the monthly hottest days observed at two meteorological stations (Lang and Ha Dong as shown in Figure 1) in Hanoi in 2016 and 2017 are listed in Table 1.Therefore, the days with the hottest air temperature are selected and considered as having the most negative effects on human health that are associated with UHI effect.From Table 1, it shows that the monthly hottest days are 2 June and 4 June in 2016 and 2017, respectively.

Acquisition and Pre-Processing of Satellite Data
We acquired the Landsat 8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) images for the same day (4 June 2017; scene ID: LC08_L1TP_127045_20170604_20170615_01_T1), when Hanoi experienced the greatest heat wave in recent 40 years.On that day, the air temperature in Hanoi reached up to 42 • C. Since Landsat 8 images were not available for 2 June 2016, the hottest day of 2016, we acquired the Landsat 8 OLI and TIRS images for the second hottest day of 1 June 2016 (scene ID: LC08_L1TP_127045_20160601_20170324_01_T1).For simplicity, 1 June is still named as the hottest day in 2016 hereafter, in the paper.On that day, the air temperature in Hanoi reached up to 39 • C. The Landsat 8 image of 2016 are later used to validate the models that are developed from the Landsat 8 images of 2017.The OLI data were converted into Top-Of-Atmosphere (TOA) spectral reflectance, applying the rescaling coefficients that are available in the metadata file.The Landsat 8 images utilized for the study area are cloud-free.

Land-Use Mapping and Validation
We adopted the ISODATA-based (Iterative Self-Organized Data analysis) unsupervised classification approach for mapping the three major land-use types, built-up area, water body, and vegetation cover, which are prevalent in inner Hanoi City.We carried out tasseled cap transformation for the Landsat 8 OLI images, and calculated three tasseled cap indices, Greenness, Wetness, and Brightness, following the methodology and transformation coefficients (Table 2) provided by Baig et al. [32].
These three tasseled cap indices were used for ISODATA-based clustering and unsupervised mapping of the land-use types for the years 2016 and 2017.Four main land-use types, including built-up, vegetation, water, and others (others mean that the pixels (regions) are mixed with built-up, vegetation, and water) were mapped.For validation, we made a common map (changed and un-changed areas) between 1 June 2016 and 4 July 2017, including four classes: built-up, vegetation, water, and others & changes.The changed areas are not stable pixels (and also needed to be checked).We chose 100 geo-location points belonging to the un-changed area of each class (built-up, Vegetation, and water) plus the others & changes areas, to validate the classification results.A total of 300 points (for three target classes) were selected.The random points were checked through visual interpretation of Google Earth images and Google Map, acquired in June 2016 and August 2017, respectively.The distribution of the chosen points is displayed in Figure 2.   Table 2. Tasseled cap transformation coefficients for Landsat 8 at-satellite reflectance of six bands: Blue, Green, red, near-infrared (NIR), shortwave infrared (SWIR 1), and shortwave infrared (SWIR 2) (adopted from Baig et al. [32]).According to the United States Geological Survey (USGS), Thermal infrared sensor (TIRS) bands of Landsat 8 have been affected by stray light from far out-of-field since its launch in 2013.A new stray light correction algorithm (SLCA) has been implemented into the USGS ground system since February 2017 and applied to reprocess historical Landsat 8 images.After SLCA implementation, Wang and Lentilucci [33] had a study to compare Landsat 8 TIRS Stray Light Correction with Multi-Sensor Measurements.It was concluded that the maximum difference in a temperature varies from 0.5% to 0.7% only.García-Santos et al. [34] validated the SLCA implementation using in situ LST measurements and three different LST estimation method algorithms (radiative transfer equation (RTE), single-channel algorithms (SCA), and split-window algorithms (SWA)) were applied for 21 scenes of Landsat 8 images.The in situ measured site is composed of different types of land covers, such as buildings, asphalt roads, vegetation regions, and so on.The study concluded that the SWA shows the best result for LST calculation from the Landsat 8 image with the lowest root mean square error (RMSE) (within 1.6-2 K).Therefore, the SWA was selected to calculate the LST for this study.

Tasseled Cap Indexes
LST maps were generated by using the thermal infrared bands 10 and 11, which are available in the Landsat 8 TIRS images (1 June 2016 and 4 June 2017).We applied the SWA (Equation ( 1)) adopted from the literature [35,36] for the generation of the LST maps: where T i and T j are the at-sensor brightness temperatures at the thermal infrared bands i and j (in Kelvins), respectively, ε is the mean emissivity, ε = 0.5(ε i + ε j ); ∆ε is the emissivity difference, ∆ε = (ε i − ε j ); w is the total atmospheric water vapor content (in g•cm −2 ); and c0 to c6 are SW coefficients to be determined from simulated data.T i and T j were calculated on the basis of the following formula (Equation ( 2)): where T is at-sensor brightness temperatures; L λ is TOA spectral radiance in W/(m 2 ster µm); K 1 and K 2 are the pre-launch calibration constants (from metadata file of Landsat 8 image).
The TOA spectral radiance (L λ ) (in Equation ( 2)) was calculated from the radiance rescaling factors provided in the metadata file, applying the following formula (Equation (3)): where M L is the band-specific multiplicative rescaling factor; A L is the band-specific additive rescaling factor; Q cal are the quantized and calibrated standard product pixel digital numbers (DN).The land surface emissivity (ε) was estimated from the Landsat 8 imagery using the Normalized Difference Vegetation Index (NDVI) threshold method [37].The total atmospheric water vapor content coefficient was obtained from NASA's Atmospheric Correction Parameter Calculator (http: //atmcorr.gsfc.nasa.gov/)[38].At the location of Hanoi (N21.00,E105.83), the total atmospheric water vapor was 5.17 g•cm −2 for 4 June 2017, and 5.31 g•cm −2 for 1 June 2016.The coefficients c0-c6 in Equation (1) were determined from the simulated data provided by Jimenez-Munoz et al. [36].It was found that the mean error of the LST was less than 1.5 K.The original LST calculated from the Equation ( 1) is in Kelvin degrees.The LST was converted to Celsius degrees for regression calculation.

Preparation of Statistical Datasets
The 30 m resolution land-use and LST maps of the years 2016 and 2017 produced in the research were used to prepare statistical datasets for inner Hanoi City.We applied the moving window method with varying window sizes, from 120 × 120 m (4 × 4 pixels) to 570 × 570 m (19 × 19 pixels), for the preparation of the datasets.Note that the LST data have a native pixel size of 100 m (resampled by USGS at 30 m), and the window averaging worsens the LST resolution and the ability to capture details in a heterogeneous area.The aim of the moving window (5 × 5 pixels) throughout the study area is to gather a sufficient number of samples (as large as 7770) for the statistical analysis, and selection of window size for further analysis.A "window" with a size and shape of interest is moved over the data with a moving distance that is equal to the window's width.For each window, the percentage coverage of the land-use types, urban built-up (%U), vegetation cover (%V), and water body (%W), were calculated.We also calculated the mean land surface temperature (µLST) for each window.The statistical datasets of 2016 were used for the validation of the model derived with 2017 datasets.

Regression Analysis, Modeling, and Validation
The statistical datasets constituting land-use composition (%U, %V, and %W) and µLST with varying window sizes (from 120 × 120 m to 570 × 570 m) of 2017 were used to analyze the relationship between land-use composition and LST by performing a regression analysis.Based on these relationships, a multivariate regression model was derived for the prediction of LST from the land-use composition data.According to the European Green Capital report [39], the number of houses located at a distance of more than 300 m from a 0.5 ha adjacent green (or larger) is considered as a basis for the evaluation of a green city.If this number is large, the green score of the city will be reduced.It means that the smallest urban area used for evaluation is around 28 ha (=3.14 × 300 m 2 ).In addition, the Ministry of Construction of Vietnam issued a Circular No. 10/2008/TT-BXD to guide the assessment and recognition of model new urban centers on 22 April 2008.In the Circular, the first requirement is that the urban area must be 50 ha or more.We assumed that a window size of 510 × 510 m, 17 × 17 pixels of Landsat 8 image, close to 25 ha (a half of the minimum requirement of the Circular No. 10) can be considered representative enough as a suitable unit for urban land management and planning.The performance of the newly derived regression model with a 510 × 510 m window size was assessed for the prediction of LST with different window sizes.The model was validated with the statistical datasets of the hottest day of 2016.Implications of the derived regression model for urban planning and design were discussed.The outline of the research flowchart is shown in Figure 3.
The model was validated with the statistical datasets of the hottest day of 2016.Implications of the derived regression model for urban planning and design were discussed.The outline of the research flowchart is shown in Figure 3.

Land Use Maps, Land Use Changes, and Validation Results
The produced 30 m resolution land use maps of Hanoi for the years 2016 and 2017 are displayed in Figures 4a and 4b, respectively.Some areas with substantial land use changes are visually chosen and labeled as A, B, C, E, and F in Figure 4c,d, with coverage percentage statistics for the three major types of land use of concern in the years 2016 and 2017, respectively.The three major land-use types include built-up areas, vegetation covers, and water bodies.It is found that the changes in land-use coverage percentage are 4.4%, −4.4%, −0.8%, and 0.8% for built-up, vegetation, water, and others, respectively, from years 2016 to 2017.The 4.4 % area of inner Hanoi City is equivalent to 704 ha (7.04 km 2 ).Obviously, the increase in built-up area is mainly contributed by the decrease in vegetation coverage by 4.4%.However, it does not mean that builtup area is completely transformed from vegetation cover, since there are other minor land use covers.Even though we only consider the land-use changes in a one year interval from 2016 to 2017, the fast change of land use in inner Hanoi City is obvious by comparing the two consecutive years of images.
After the quantified areas of land use changes, qualitative land use changes may be easily observed by comparing the Google Earth images, which are conveniently available online, as well as field survey photos, as shown in Figure 5. Field survey photos were taken on 14 October 2018.The transformation of vegetation cover into built-up areas is easily detected by bare eyes.Figure 5a shows Google Earth image of Inner Hanoi City with areas B and E for the demonstration of land use changes and (red) dots with field survey photos as the ground truth for reference.Figures 5b and 5c are Google Earth images of area B that were acquired in June 2016 and August 2017, respectively, with four field survey photos showing the current in-situ land use status. Figure 5d

Land Use Maps, Land Use Changes, and Validation Results
The produced 30 m resolution land use maps of Hanoi for the years 2016 and 2017 are displayed in Figure 4a,b, respectively.Some areas with substantial land use changes are visually chosen and labeled as A, B, C, E, and F in Figure 4c,d, with coverage percentage statistics for the three major types of land use of concern in the years 2016 and 2017, respectively.The three major land-use types include built-up areas, vegetation covers, and water bodies.It is found that the changes in land-use coverage percentage are 4.4%, −4.4%, −0.8%, and 0.8% for built-up, vegetation, water, and others, respectively, from years 2016 to 2017.The 4.4 % area of inner Hanoi City is equivalent to 704 ha (7.04 km 2 ).Obviously, the increase in built-up area is mainly contributed by the decrease in vegetation coverage by 4.4%.However, it does not mean that built-up area is completely transformed from vegetation cover, since there are other minor land use covers.Even though we only consider the land-use changes in a one year interval from 2016 to 2017, the fast change of land use in inner Hanoi City is obvious by comparing the two consecutive years of images.
After the quantified areas of land use changes, qualitative land use changes may be easily observed by comparing the Google Earth images, which are conveniently available online, as well as field survey photos, as shown in Figure 5. Field survey photos were taken on 14 October 2018.The transformation of vegetation cover into built-up areas is easily detected by bare eyes.Figure 5a shows Google Earth image of Inner Hanoi City with areas B and E for the demonstration of land use changes and (red) dots with field survey photos as the ground truth for reference.Figure 5b,c    The accuracy assessment of LULC classification was performed to quantitatively assess how effective the pixels were sampled into the correct LULC classes.The accuracy assessment of pixel selection was on areas that could be clearly identified on both Landsat high-resolution images, and on Google Earth and Google Map.A total of 300 points (locations) for both years 2016 and 2017 were created in classification images of the study area, as presented in Section 2.4.The study had an overall classification accuracy of 92%, a kappa coefficient (K) of 0.88, an overall accuracy of 93%, and kappa coefficient (K) of 0.9, 0 for 2016 and 2017, respectively.The confusion matrices of LULC classification maps derived from Landsat 8 images in 2016 and 2017 are shown in Tables 3 and 4, respectively.The kappa coefficient is rated as being substantial, and hence, the classified image was found to be appropriate for further analysis.The accuracy assessment of LULC classification was performed to quantitatively assess how effective the pixels were sampled into the correct LULC classes.The accuracy assessment of pixel selection was on areas that could be clearly identified on both Landsat high-resolution images, and on Google Earth and Google Map.A total of 300 points (locations) for both years 2016 and 2017 were created in classification images of the study area, as presented in Section 2.4.The study had an overall classification accuracy of 92%, a kappa coefficient (K) of 0.88, an overall accuracy of 93%, and kappa coefficient (K) of 0.9,0 for 2016 and 2017, respectively.The confusion matrices of LULC classification maps derived from Landsat 8 images in 2016 and 2017 are shown in Tables 3 and 4, respectively.The kappa coefficient is rated as being substantial, and hence, the classified image was found to be appropriate for further analysis.Overall accuracy: 93% Kappa coefficient: 0.90 Table 5 shows the land-use classification and its changes from years 2016 to 2017 in inner Hanoi City.It reveals a big increase for built-up areas by 8.9% in a year, but a decrease for vegetation, water, and others by 10.9%, 7.0%, and 4.0%, respectively.For assessing the quantity of transformations among land-use types from the years 2016 to 2017, a transformation matrix is determined, as shown in Table 6.Two key points are observed from the table.First, increased built-up area 841 ha is primarily changed from vegetation (618 ha), followed by others (137 ha) and water (86 ha).Second, the transformations of built-up to vegetation, water, and others are 0.33% (27/7937), 0.27% (22/7937), and 0.92% (74/7937), respectively.It is evident that, under the Master Plan 2030, the inner Hanoi City has been changing quickly with time, since the implementation of the Plan in 2011.This can be easily justified by the dramatic land cover changes in a year, by comparing the land-use maps of years 2016 and 2017.Under such a fast-changing land-use situation, any suggested measures for mitigation of UHI for the inner Hanoi City are urgently required.June 2017, respectively.That is, the average LST was slightly higher for the hottest days in 2016 than 2017 by 0.8 • C. Note that even with an overall warmer thermal environment for the whole of inner Hanoi City in 2016 than 2017, LSTs were lower in the areas with significant transformation of land-use from vegetation cover to built-up, labeled as B, E, and F, in 2016 than 2017 by 0.9, 0.8, and 2.5 • C, respectively.The number of pixels and mean LST for the three chosen areas B, E, and F on the two hottest days in 2016 and 2017 and the difference in LST between the two hottest days are shown in Table 7. Results indicate that the transformation of land-use from vegetation to built-up has enhanced the contrast in thermal signatures, i.e., LSTs, by 1.7, 1.6, and 3.3 • C in the three fast-changing land-use regions B, E, and F, respectively, in one year between years 2016 and 2017.The contrast confirms the cooling effect of vegetation cover on the SUHI intensity that cannot be overstressed in the hot urban cites.As for the signature of the traditionally defined SUHI, it is analyzed here by assuming the boundary of Hanoi's downtown, that consists of seven districts, including Tay Ho, Hoan Kiem, Ba Dinh, Dong Da, Hai Ba Trung, Cau Giay, and Thanh Xuan Districts, as a mask to divide the urban  As for the signature of the traditionally defined SUHI, it is analyzed here by assuming the boundary of Hanoi's downtown, that consists of seven districts, including Tay Ho, Hoan Kiem, Ba Dinh, Dong Da, Hai Ba Trung, Cau Giay, and Thanh Xuan Districts, as a mask to divide the urban (pink color) and suburban areas (green color), as shown in Figure 7.The statistics of urban and suburban LSTs are given in Table 8.The STDs in Table 8 represent 8 indicates the mean LST of the urban area was higher by 0.53 °C than that in the suburban area in 2016.In contrast, in 2017, the LST of the urban area was lower than that in the suburban area, by 0.21 °C.This was due to the massive land use change in the suburban area, possibly resulting from the implementation of Hanoi Master Plan 2030, where many buildings are being constructed.That is, areas with vegetation and water surfaces in the suburban area are significantly reduced.Table 9 shows that the built-up area is increased by 4.28 km 2 and vegetation decreased by 4.01 km 2 in the suburban belt in only a one year interval from 2016 to 2017.Consequently, fast urbanization in suburban area results in slightly higher LSTs in the suburban Table 8 indicates the mean LST of the urban area was higher by 0.53 • C than that in the suburban area in 2016.In contrast, in 2017, the LST of the urban area was lower than that in the suburban area, by 0.21 • C.This was due to the massive land use change in the suburban area, possibly resulting from the implementation of Hanoi Master Plan 2030, where many buildings are being constructed.That is, areas with vegetation and water surfaces in the suburban area are significantly reduced.Table 9 shows that the built-up area is increased by 4.28 km 2 and vegetation decreased by 4.01 km 2 in the suburban belt in only a one year interval from 2016 to 2017.Consequently, fast urbanization in suburban area results in slightly higher LSTs in the suburban area than in urban area in 2017.This finding is also in line with the comparative results of the LST in the chosen areas B, E, and F, which are located in the suburban belt to demonstrate a negative SUHI in response to urban expansion in inner Hanoi City.

Individual Land Use Coverage versus LST
Results from the regression analysis show that there exist high correlations between the percentage coverage of the individual land-use types (%U, %V, and %W) and the mean LSTs (µLSTs) for all window sizes considered in the research.However, in our study area, we detect that smaller window sizes (less than 300 × 300 m) do not properly represent a heterogeneous mixture of the land-use composition.
In most cases, one or two components of the land-use composition (%W, %U, and %V) are absent.Therefore, we present the analysis results pertaining to window sizes of larger than 300 × 300 m.On the other hand, the correlations between the percentage coverage of the land use types (%U, %V, and %W) and µLSTs are much stronger for large than for small window sizes.
Figure 8 demonstrates the relationships between the percentage coverage of the land-use types (%U, %V, and %W) and µLSTs in the case of 510 × 510 m window size on the hottest day of 2017.The regression coefficients are found to be the highest for the water coverage (%W) with R 2 = 0.70, followed by urban built-up cover (%U), with R 2 = 0.67, and then vegetation cover (%V), with R 2 = 0.43.Note that water (%W) and vegetation (%V) coverages display negative correlations with µLST, as they play a known cooling effect on the SUHI intensity, whereas urban built-up coverage (%V) is positively correlated with µLST to enhance the UHI phenomena.

Individual Land Use Coverage versus LST
Results from the regression analysis show that there exist high correlations between the percentage coverage of the individual land-use types (%U, %V, and %W) and the mean LSTs (µLSTs) for all window sizes considered in the research.However, in our study area, we detect that smaller window sizes (less than 300 × 300 m) do not properly represent a heterogeneous mixture of the land-use composition.In most cases, one or two components of the land-use composition (%W, %U, and %V) are absent.Therefore, we present the analysis results pertaining to window sizes of larger than 300 × 300 m.On the other hand, the correlations between the percentage coverage of the land use types (%U, %V, and %W) and µLSTs are much stronger for large than for small window sizes.
Figure 8 demonstrates the relationships between the percentage coverage of the land-use types (%U, %V, and %W) and µLSTs in the case of 510 × 510 m window size on the hottest day of 2017.The regression coefficients are found to be the highest for the water coverage (%W) with R 2 = 0.70, followed by urban built-up cover (%U), with R 2 = 0.67, and then vegetation cover (%V), with R 2 = 0.43.Note that water (%W) and vegetation (%V) coverages display negative correlations with µLST, as they play a known cooling effect on the SUHI intensity, whereas urban built-up coverage (%V) is positively correlated with µLST to enhance the UHI phenomena.

Land-Use Composition-Driven Prediction of LST
A multivariate regression model between land-use composition (percentage coverage of urban built-up area (%U), vegetation (%V), and water (%W)) and LST on the hottest day of 2017 (4 June) is derived.For this purpose, we utilized a 510 × 510 m window size, close to 25 ha, assumed to be a suitable unit for urban land management and planning.The resultant regression function is given in Equation ( 4

Land-Use Composition-Driven Prediction of LST
A multivariate regression model between land-use composition (percentage coverage of urban built-up area (%U), vegetation (%V), and water (%W)) and LST on the hottest day of 2017 (4 June) is derived.For this purpose, we utilized a 510 × 510 m window size, close to 25 ha, assumed to be a suitable unit for urban land management and planning.The resultant regression function is given in Equation ( 4 The performance of the above derived regression model (Equation ( 4)) for a 25 ha land management unit is evaluated by comparing its predicted LST for different window sizes with that derived from Landsat-8 data (LST LS8 ).It was found that the regression function gave reasonable results, with R 2 ~0.9 for all cases (window sizes from 300 × 300 m to 570 × 570 m).Almost 90% of the variation in LST can be predicted by the model for window sizes of 330 × 330 m and 510 × 510 m (Figure 9).Since the R 2 alone may not be sufficient to evaluate the performance and accuracy of a prediction model, the RMS error (RMSE) and the ratio between RMSE and the standard deviation (STD) of the Landsat 8 LST map after windowing and averaging, i.e., RMSE/STD_LST LS8 , are also reported.Note that if the ratio between RMSE/STD_LST LS8 is lower than 0.5, the model can be considered reliable.It is found that the RMSE and RMSE/STD_LST LS8 for the window sizes of 510 × 510 m and 330 × 330 m are 1.77 • C, 0.38, and 1.83 • C, 0.4, respectively.Such statistics justify that the regression model that is derived by using the 510 × 510 m window size can be applied for the prediction of LST for the other window sizes as well.
Remote Sens. 2018, 10, x FOR PEER REVIEW 15 of 20 results, with R 2 ~ 0.9 for all cases (window sizes from 300 × 300 m to 570 × 570 m).Almost 90% of the variation in LST can be predicted by the model for window sizes of 330 × 330 m and 510 × 510 m (Figure 9).Since the R 2 alone may not be sufficient to evaluate the performance and accuracy of a prediction model, the RMS error (RMSE) and the ratio between RMSE and the standard deviation (STD) of the Landsat 8 LST map after windowing and averaging, i.e., RMSE/STD_LSTLS8, are also reported.Note that if the ratio between RMSE/STD_LSTLS8 is lower than 0.5, the model can be considered reliable.It is found that the RMSE and RMSE/STD_LSTLS8 for the window sizes of 510 × 510 m and 330 × 330 m are 1.77 °C, 0.38, and 1.83 °C, 0.4, respectively.Such statistics justify that the regression model that is derived by using the 510 × 510 m window size can be applied for the prediction of LST for the other window sizes as well.

Validation of the Land-Use Driven Model
We test the derived regression model (Equation 4) for its predicted LSTs on another hottest day (1 June 2016) by using the statistical dataset of 2016. Figure 10

Validation of the Land-Use Driven Model
We test the derived regression model (Equation ( 4)) for its predicted LSTs on another hottest day (1 June 2016) by using the statistical dataset of 2016. Figure 10

Discussion
Hanoi is one of the hottest cities in Asia.It is highly vulnerable to the detrimental consequences of heat islands on urban ecology and human health.On 4 June 2017, Hanoi experienced a heat wave to the hottest level in recent 40 years, reaching 42 °C [40].According to the study by Nguyen et al. [41], Vietnam's average temperature has increased at a rate of 0.26 °C per decade since the 1970s, which is approximately twice the rate of global warming over the same period of time.Over the last century, Hanoi has been experiencing a drastic increase in population, high-speed urbanization, and transformation of natural forests into urban built-up areas.All of these socioeconomic and biophysical changes have had a critical impact on the urban microclimate.As a result, severe heat waves have become common during the summer months in recent years.UHI is mainly caused by the modification of land surfaces and the concentration of the population [42].The process of land modification generally increases the use of materials that retain heat (concrete surfaces) and sacrifices air-conditioning elements such as trees and water surfaces and moreover, waste heat is generated by energy usage [43].In tropical cities, where there are a great number of very hot days, UHI causes the temperature to rise dramatically, negatively affecting the quality of the living environment in these cities.
By assuming that the unit of urban land management is around 25 ha (510 × 510 m), on a very hot day such as 4 June 2017, the LST (°C) can be predicted with respect to land-use composition using the regression model derived in this research, as shown in Table 10.

Discussion
Hanoi is one of the hottest cities in Asia.It is highly vulnerable to the detrimental consequences of heat islands on urban ecology and human health.On 4 June 2017, Hanoi experienced a heat wave to the hottest level in recent 40 years, reaching 42 • C [40].According to the study by Nguyen et al. [41], Vietnam's average temperature has increased at a rate of 0.26 • C per decade since the 1970s, which is approximately twice the rate of global warming over the same period of time.Over the last century, Hanoi has been experiencing a drastic increase in population, high-speed urbanization, and transformation of natural forests into urban built-up areas.All of these socioeconomic and biophysical changes have had a critical impact on the urban microclimate.As a result, severe heat waves have become common during the summer months in recent years.UHI is mainly caused by the modification of land surfaces and the concentration of the population [42].The process of land modification generally increases the use of materials that retain heat (concrete surfaces) and sacrifices air-conditioning elements such as trees and water surfaces and moreover, waste heat is generated by energy usage [43].In tropical cities, where there are a great number of very hot days, UHI causes the temperature to rise dramatically, negatively affecting the quality of the living environment in these cities.
By assuming that the unit of urban land management is around 25 ha (510 × 510 m), on a very hot day such as 4 June 2017, the LST ( • C) can be predicted with respect to land-use composition using the regression model derived in this research, as shown in Table 10.A substantial impact of land-use composition on LST and UHI is apparent from Tables 8 and 10.This analysis shows that the LST at each management unit (e.g., 25 ha) can be decreased from 45.23 • C (100% urban built-up cover) to 37.91 • C (100% vegetative cover), and then to 28.24 • C (100% water cover).The urban planner and designer indeed can neither change the LST directly to mitigate the UHI effects, and nor would it be possible to convert all urban built-up coverage in the cities to water bodies or vegetative areas.However, even a 20% conversion of urban built-up areas into vegetative (10%) and water bodies (10%) could reduce LST by 2.43 • C. The results of this research present an opportunity for urban planners and designers-the LST and associated effects of SUHI can be adjusted by managing the land-use composition and percentage coverage of the individual land-use types (%U, %V, and %W) in each urban land management unit.Reduction in the greenery has been described as a major cause of the rising temperature in cities by many researchers.It has been suggested that urban greening activities are the efficient nature-friendly solutions for mitigating UHI effects [44][45][46].Such suggestions are also justified by the observed negative correlations between LST and NDVI, resulting in increased eco-environmental vulnerability [47][48][49][50][51].

Conclusion
We analyzed the correlations between land-use composition, the percentage coverage of three major land use types (%U, %V, and %W), and LST for different window sizes or urban land management units.The land use maps showed that inner Hanoi experienced substantial changes in land-use in a one-year interval from years 2016 to 2017, likely influenced by the City's Master Plan 2030.High correlations were observed between the percentage coverage of each land-use type (i.e., %U, %V, and %W) and LST.They are coherent with the findings of the previous studies, while considering different window sizes on the hottest day of 2017 (4 June).The vegetation and water coverages, acting as cooling effect, exhibited negative correlations with LST.In contrast, the urban coverage was proportionate to the LST.Consequently, a multivariate regression model was derived by taking into account the correlations between land-use composition and LST.We successfully validated our model for the prediction of LST on another hot day (1 June 2016).The validated model in this research provides an opportunity for urban planners and designers to suggest measures for adjusting the LST and the associated effects of SUHI by managing the land use composition and percentage coverage of the individual land-use types (%U, %V, and %W).This model demonstrates that a substantial decrease in LST can be achieved by increasing the proportion of vegetation and water for improved comfortable and sustainable living in the city.Despite this, we have provided some interesting findings, and we did not deeply analyze the traditionally-defined SUHI.In addition, impacts of currently rapid urbanization on SUHI and its consequence in suburban area of inner Hanoi can be further investigated.

Figure 1 .
Figure 1.Location of the study area: (a) Location of the study area in Vietnam; (b) detail of the study area with the boundary of Hanoi City before and after 2008; (c) inner Hanoi city (yellow polygon) on the composited Landsat 8 satellite image acquired on 4 June 2017, and locations of the meteorological stations Lang (yellow) and Ha Dong (orange) used to acquire air temperature data.

Figure 1 .
Figure 1.Location of the study area: (a) Location of the study area in Vietnam; (b) detail of the study area with the boundary of Hanoi City before and after 2008; (c) inner Hanoi city (yellow polygon) on the composited Landsat 8 satellite image acquired on 4 June 2017, and locations of the meteorological stations Lang (yellow) and Ha Dong (orange) used to acquire air temperature data.

Figure 2 .
Figure 2. Distribution of the ground truth points on a common area map chosen in the research.

Figure 2 .
Figure 2. Distribution of the ground truth points on a common area map chosen in the research.

Figure 3 .
Figure 3. Outline of the research flowchart.
and 5e are Google Earth images of area E, acquired in June 2016 and August 2017, respectively, with two field survey photos confirming the current in situ land use status of either built-up or ongoing construction circumstances.

Figure 3 .
Figure 3. Outline of the research flowchart.
are Google Earth images of area B that were acquired in June 2016 and August 2017, respectively, with four field survey photos showing the current in-situ land use status. Figure 5d,e are Google Earth images of area E, acquired in June 2016 and August 2017, respectively, with two field survey photos confirming the current in situ land use status of either built-up or ongoing construction circumstances.

Figure 4 .
Figure 4. Land-use maps of inner Hanoi city for (a) 1 June 2016, and (b) 4 June 2017.Areas with substantial land-use changes are visually chosen and labeled as A, B, C, E, and F in (c) and (d) with coverage percentage statistics for the three major land-use types of concern and others in the years 2016 and 2017, respectively.

Figure 4 .
Figure 4. Land-use maps of inner Hanoi city for (a) 1 June 2016, and (b) 4 June 2017.Areas with substantial land-use changes are visually chosen and labeled as A, B, C, E, and F in (c) and (d) with coverage percentage statistics for the three major land-use types of concern and others in the years 2016 and 2017, respectively.

Figure 4 . 20 Figure 5 .
Figure 4. Land-use maps of inner Hanoi city for (a) 1 June 2016, and (b) 4 June 2017.Areas with substantial land-use changes are visually chosen and labeled as A, B, C, E, and F in (c) and (d) with coverage percentage statistics for the three major land-use types of concern and others in the years 2016 and 2017, respectively.

Figure 5 .
Figure 5. Field survey photos taken on 14 October 2018 are overlaid onto Google Earth images.(a) Google Earth image of Inner Hanoi City with areas B and E for the demonstration of land use changes and (red) dots with field survey photos as the ground truth for reference; (b,c) are Google Earth images of area B acquired in June 2016 and August 2017, respectively, with four field survey photos confirming the current built-up status; (d,e) are Google Earth images of area E acquired in June 2016 and August 2017, respectively, with two field survey photos confirming the current built-up status and on-going construction circumstances.

Figure 6
Figure 6 shows LST maps of the inner Hanoi City for (a) 1 June 2016 and (b) 4 June 2017, and their corresponding histogram distributions in (c) and (d), respectively.Three areas with substantial changes in LST are randomly chosen and labeled as B, E, and F in both (a) and (b) for a further interpretation of SUHI intensity's variation.It is found that the average LSTs are 40.9• C and 40.1 • C for 1 June 2016 and 4 June 2017, respectively.That is, the average LST was slightly higher for the hottest days in 2016 than

20 Figure 6 .
Figure 6.Land surface temperature (LST) maps of inner Hanoi City for (a) 1 June 2016 and (b) 4 June 2017; and their corresponding histogram distributions in (c) and (d), respectively.Three areas with substantial changes in LST between two hottest days are randomly chosen and labeled as B, E, and F in both (a,b).The pixel size is 30 × 30 m.

Figure 6 .Table 7 .
Figure 6.Land surface temperature (LST) maps of inner Hanoi City for (a) 1 June 2016 and (b) 4 June 2017; and their corresponding histogram distributions in (c,d), respectively.Three areas with substantial changes in LST between two hottest days are randomly chosen and labeled as B, E, and F in both (a,b).The pixel size is 30 × 30 m. Table 7. Number of pixels (1 pixel = 30 × 30 m) and mean LST for the three randomly chosen areas B, E, and F on the two hottest days in 2016 and 2017, and the difference in LST between the two hottest days.

20 Figure 7 .
Figure 7.The map shows downtown of inner Hanoi (pink color) and suburban (green color).

Figure 7 .Table 8 .
Figure 7.The map shows downtown of inner Hanoi (pink color) and suburban (green color).Table 8.LSTs of urban and suburban areas of inner Hanoi.STD stands for standard deviation.Areas Mean LST ( • C)/STD 2016 Mean LST ( • C)/STD 2017

Figure 9 .
Figure 9. Land-use composition-driven predictions of LST for the hottest day of 2017 (4 June) versus LST observed by Landsat-8 for (a) 510 × 510 m and (b) 330 × 330 m window sizes.
shows the land-use composition-driven predictions of LST for the hottest day of 2016 (1 June) versus the LST observed by Landsat-8 for (a) 510 × 510 m and (b) 330 × 330 m window sizes.High correlations of 0.897 and 0.869 for 510 × 510 m and 330 × 330 m window sizes, respectively, were obtained.RMSEs for the window sizes 510 × 510 m and 330 × 330 m were 1.71 °C and 1.94 °C, respectively.The ratios RMSE/STD_LSTLS8 are also low (below 0.5), at 0.44 and 0.5 for the window sizes 510 × 510 m and 330 × 330 m, respectively.They indicate that the regression model reasonably predicts the LST for the hottest day of 2016.That is, the developed model can be used to retrieve LST for the needs of land-use management and planning with high reliability.

Figure 9 .
Figure 9. Land-use composition-driven predictions of LST for the hottest day of 2017 (4 June) versus LST observed by Landsat-8 for (a) 510 × 510 m and (b) 330 × 330 m window sizes.

Figure 10 .
Figure 10.Land-use composition-driven predictions of LST for the hottest day of 2016 (1 June) versus LST observed by Landsat-8 for (a) 510 × 510 m and (b) 330 × 330 m window sizes.

Figure 10 .
Figure 10.Land-use composition-driven predictions of LST for the hottest day of 2016 (1 June) versus LST observed by Landsat-8 for (a) 510 × 510 m and (b) 330 × 330 m window sizes.

Table 1 .
Maximum monthly air temperatures (°C) in Hanoi in 2016 and 2017.

Table 1 .
Maximum monthly air temperatures ( • C) in Hanoi in 2016 and 2017.

Table 3 .
Confusion matrix of land-use map classified from the Landsat 8 image in 2016.

Table 3 .
Confusion matrix of land-use map classified from the Landsat 8 image in 2016.

Table 4 .
Confusion matrix of land-use map classified from the Landsat 8 image in 2017.

Table 5 .
Land-use changes between years 2016 and 2017 in inner Hanoi.

Table 6 .
Transformation among land use types from years 2016 and 2017 in inner Hanoi.

Table 7 .
Number of pixels (1 pixel = 30 × 30 m) and mean LST for the three randomly chosen areas B, E, and F on the two hottest days in 2016 and 2017, and the difference in LST between the two hottest days.

Table 8 .
LSTs of urban and suburban areas of inner Hanoi.STD stands for standard deviation.

Table 9 .
Land use/land cover (LULC) change in urban and suburban areas of Hanoi in 2016-2017.

Table 10 .
Land-use composition-driven prediction of land surface temperature.

Table 10 .
Land-use composition-driven prediction of land surface temperature.