Urban Heat Island Simulations in Guangzhou, China, Using the Coupled Wrf/ucm Model with a Land Use Map Extracted from Remote Sensing Data

The Weather Research and Forecasting (WRF) model coupled with an Urban Canopy Model (UCM) was used for studying urban environmental issues. Because land use data employed in the WRF model do not agree with the current situation around Guangzhou, China, the performance of WRF/UCM with new land-use data extracted from Remote Sensing (RS) data was evaluated in early August 2012. Results from simulations reveal that experiments with the extracted data are capable of reasonable reproductions of the majority of the observed temporal characteristics of the 2-m temperature, and can capture the characteristics of Urban Heat Island (UHI). The " UCM_12 " simulation, which employed the extracted land-use data with the WRF/UCM model, provided the best reproduction of the 2-m temperature data evolution and the smallest minimum absolute average error when compared with the other two experiments without coupled UCM. The contributions of various factors to the UHI effect were analyzed by comparing the energy equilibrium processes of " UCM_12 " in urban and suburban areas. Analysis revealed that energy equilibrium processes with new land use data can explain the diurnal character of the UHI intensity variation. Furthermore, land use data extracted from RS can be used to simulate the UHI.


Introduction
Urbanization is a population shift from rural to urban areas, and also covers the ways in which society adapts to this change.It predominantly results in land use and land cover changes and increased building density, both horizontally and vertically, in urban areas.Land use and land cover changes associated with urbanization can have a considerable impact on urban climate [1].One phenomenon is a consistent rise in temperature in the urban atmosphere caused by the land use changes, and the accompanying effects on the physical processes governing energy, momentum, and matter exchanges between the land surface and the atmosphere [2,3].Guangzhou, the capital and largest city of Guangdong Province and the third largest city in China, is undergoing rapid and widespread urbanization and a deterioration of the urban environment caused by Urban Heat Island (UHI).
The Weather Research and Forecasting (WRF) model, coupled with an Urban Canopy Model (UCM), was developed as a community tool for studying urban environmental issues [4,5].The coupled WRF/UCM modeling system has been applied to many regions, such as Houston, USA [6], Nanjing, China [7], and New York, USA [8], and its performance has been evaluated against observations.Many studies have performed simulations using land use data from different years with the single-layer urban canopy model in the Noah land surface model to demonstrate that urban development and the accompanying land use changes can make a significant contribution to extreme heat events [9][10][11][12].However, the default data (U.S. Geological Survey (USGS) land use) employed in the WRF model are based on the NationalOceanic and Atmospheric Administration (NOAA) 1-km Advanced Very High Resolution Radiometer (AVHRR) data obtained from 1992 to 1993.Owing to the rapid urban expansion in the past decade in China, the USGS data are considered outdated.Numerical simulations in Guangzhou [9] and Chengdu [13] with the new land use data reproduced better 2-m temperature evolution data with a smaller minimum absolute average error compared with results obtained using the default USGS land use data in the WRF model.Versions of the WRF model after version 3.1 provide an alternative land use dataset based on the Moderate Resolution Imaging Spectroradiometer (MODIS) 2001 satellite products.Implementation of MODIS data results in better performance of the coupled WRF/UCM modeling system, although the urban area in MODIS also falls short of reality because of fast urbanization [14].The diurnal variations of UHI intensity and the spatial distribution of the UHI effect in Beijing have been reproduced well by the WRF/UCM with MODIS data [15].The simulation results with MODIS data improved predictions of the accumulated rainfall when compared with the simulation performed using USGS data [16].Numerical simulations of Taiwan showed that results obtained using MODIS land use data are in better agreement with the observed data than those obtained using USGS data, although MODIS land use data overestimated the urban area [17].When compared with observational data, numerical simulations in Guangzhou [18] showed improved accuracy using MODIS data compared with USGS data, but poorer accuracy than that obtained using land use data extracted from the Landsat-7 remote sensing dataset.The effect of urbanization on the weather and climate in Hangzhou were investigated using the coupled WRF/UCM model with updated land use data, and the research showed that updated land use data can reproduce the local climate accurately, and that urban land use has a significant impact on the simulated UHI effect [19].
In this study, the land-use data were extracted from the remote sensing dataset of Landsat-7 using a previous research [20] method and then up-to-date extracted land-use data classified as urban land cover were divided into three urban subcategories by satellite-measured night time light data and the normalized difference vegetation index dataset.The mesoscale coupled WRF/UCM model with different land-use data are used to simulate the formation of high-temperature synoptic conditions in the area of Guangzhou in early August 2012.The aim of this research is to demonstrate that, when using the extracted and up-to-date urban land use data from a remote sensing dataset, the WRF/UCM modeling system provides a more accurate simulation of urban temperatures and the UHI effect in Guangzhou.Furthermore, the study aims to simulate the formation of high-temperature synoptic conditions in the area of Guangzhou, and to investigate the properties of coupled WRF/UCM model and their effects on the UHI simulations.

Study Area
The study area is located in Guangzhou, south China, with a population of 12.78 million as of the 2010 census.Located in the south-central portion of Guangdong, Guangzhou spans between 112 ˝57 1 and 114 ˝03 1 E longitudes, and 22 ˝26 1 and 23 ˝56 1 N latitudes.Located just south of the Tropic of Cancer, Guangzhou has a humid subtropical climate that is influenced by the East Asian monsoon.Summers are wet with high temperatures, and high heat index.The hottest period in summer usually lasts from 11 July to 20 August, with an average period of 41 days.The mean annual temperature in the city is 22.3 ˝C, with monthly daily averages ranging from 12.3 ˝C in January to 28.4 ˝C in July.
Guangzhou is expanding at a very high speed.Statistical data from the China City Statistical Yearbook show that the built-up area in Guangzhou city has increased more than 4 times since the year 2000, which is much faster than the average expansion pace in China.Such tremendous urban expansion accompanied with increasing built-up areas and intensifying human activities causes remarkable modifications of the underlying surface properties, thereby significantly enhancing the UHI effect [17].Previous studies show that the UHI effect noticeably contributes to regional warming in Guangzhou area [21].This study area includes the seven main districts of Guangzhou, all of which are undergoing rapid urban expansion.

Model Configuration
The WRF version 3.2 model, coupled with a sophisticated single-layer UCM, was used in this study [5].Three one-way nested domains with 72 ˆ72, 101 ˆ101, and 121 ˆ121 grid points, and grid spacings of 25, 5, and 1 km, respectively, were used (Figure 1).The coordinate of the domain's center is 22 ˝61 57 11 N by 113 ˝31 1 8 11 E. The largest domain (Domain 1) covers most of southern China and the second domain covers the whole Guangdong area.The main districts of Guangzhou are covered by Domain 3. The integration starts at 0:00 GMT on 20 July 2012 for a period of one month.Four days from 1 August to 4 August were selected to evaluate performance of the WRF/UCM with different land-use data.year 2000, which is much faster than the average expansion pace in China.Such tremendous urban expansion accompanied with increasing built-up areas and intensifying human activities causes remarkable modifications of the underlying surface properties, thereby significantly enhancing the UHI effect [17].Previous studies show that the UHI effect noticeably contributes to regional warming in Guangzhou area [21].This study area includes the seven main districts of Guangzhou, all of which are undergoing rapid urban expansion.

Model Configuration
The WRF version 3.2 model, coupled with a sophisticated single-layer UCM, was used in this study [5].Three one-way nested domains with 72 × 72, 101 × 101, and 121 × 121 grid points, and grid spacings of 25, 5, and 1 km, respectively, were used (Figure 1).The coordinate of the domain's center is 22°6′57''N by 113°31′8''E.The largest domain (Domain 1) covers most of southern China and the second domain covers the whole Guangdong area.The main districts of Guangzhou are covered by Domain 3. The integration starts at 0:00 GMT on 20 July 2012 for a period of one month.Four days from 1 August to 4 August were selected to evaluate performance of the WRF/UCM with different land-use data.The parameterization schemes used in our simulation are listed in Table 1, including long-and short-wave radiation processes, planetary boundary layer processes, land surface processes, microphysical processes, etc.In the WRF model, a single-layer UCM was implemented in the NOAH land surface model to account for the thermal and dynamic effects of urban buildings, including the trapping of radiation within the urban canopy.The UCM standard values for heat capacity, conductivity, albedo, and emissivity roughness length for heat and momentum of roof, road, and wall surfaces were employed [4,22].Initial and lateral boundary conditions were taken from the National Centers for Environmental Prediction (NCEP) Global Forecast System Final Analyses (horizontal resolution of 1° × 1°) with 6-h intervals.The parameterization schemes used in our simulation are listed in Table 1, including longand short-wave radiation processes, planetary boundary layer processes, land surface processes, microphysical processes, etc.In the WRF model, a single-layer UCM was implemented in the NOAH land surface model to account for the thermal and dynamic effects of urban buildings, including the trapping of radiation within the urban canopy.The UCM standard values for heat capacity, conductivity, albedo, and emissivity roughness length for heat and momentum of roof, road, and wall surfaces were employed [4,22].Initial and lateral boundary conditions were taken from the National Centers for Environmental Prediction (NCEP) Global Forecast System Final Analyses (horizontal resolution of 1 ˝ˆ1 ˝) with 6-h intervals.

Updated Land Use Data
To derive and classify each land use class from the Landsat-7 ETM+ images, a supervised classification with the maximum likelihood method, which is the same method as the previous research [20], was only applied to the Domain 3. Supervised classification is a procedure most often used for quantitative analysis of remote sensing image data.It is based on using suitable algorithms to label the pixels in an image as representing particular ground cover classes [23].The extracted land use classes are essentially consistent with the ones defined from MODIS land use classes of forest, grassland, wetland, cropland, urban, barren land, and water.These land use data are almost the same as the previous research result, here named RS_12 [20].There is only one urban land category in RS_12 and thus these data cannot reproduce urban effects on local climate caused by heterogeneity in an urban area.A human settlement index (HSI) [24], based on satellite-measured nighttime light imagery and the normalized difference vegetation index (NDVI), was used to divide the extracted urban land cover into three urban subcategories: commercial/industrial/transportation (CIT), high-intensity residential (HIR), and low-intensity residential (LIR).
HSI data were obtained from a combination of nighttime light imagery and the normalized difference vegetation index (NDVI) [24], as expressed in Equation (1).

HSI "
p1 ´NDVI max q `OLS nor 1 ´OLS nor `NDVI max `NDVI max ˆOLS nor (1) where OLS nor is the normalized value of the 2012 Defense Meteorological Satellite Program/Operational Linescan System (DMSP/OLS) image [25,26], expressed by Equation (2).A threshold of 23 was selected to minimize the effects of blooming, a phenomenon observed in DMSP/OLS nighttime imagery.OLS max is the maximum value in the 2012 DMSP/OLS nighttime light image.
OLS nor " pOLS ´23q pOLS max ´23q (2) To separate human settlements from bare land effectively, and to remove the effect of cloud contamination, 2012 multi-temporal SPOT NDVI images at a resolution of 1 km ˆ1 km were used to generate a new NDVI composite index, as shown in Equation (3) [27].
HSI values were used to classify urban land in Guangzhou into the following subcategories: (1) pixels with HSI values ě80th percentile, represented as CIT areas; (2) pixels with HSI values between the 30th and the 80th percentiles, represented as HIR areas; and (3) pixels with HSI values ď30th percentile, represented as LIR areas.The extracted urban land use data were classified into the three urban land cover categories to represent up-to-date urbanization conditions in Guangzhou.These Landsat-7 extracted and then enhanced land use data, named UCM_12, are shown in Figure 2a and the default Modis land use data are shown in Figure 2b.The area in Figure 2 indicates the square in Domain 3.

Numerical Experiments
To understand how well the WRF/UCM simulates the UHI effect in urban environments, three numerical experiments were designed to run with or without the UCM model, and using either the MODIS data or the extracted land use data (Table 2).Sensitivity experiments were conducted as follows: (1) MODIS experiment, integrated solely by the WRF model not coupled with UCM and based on the MODIS land use data, the detailed land use map is shown in Figure 2b; (2) RS_12 experiment, also integrated without UCM and based on the Landsat-7 extracted land use data for 2012; and (3) an experiment based on the enhancement of the extracted land use data from Landsat-7 with HSI and coupled with the UCM model.The physics schemes remained the same for all runs.

Observational Data
The simulation lasts for one month, as obtaining long time data to evaluate performance of the WRF/UCM and the land use data is better.Only four days of climate data from 1 August 0:00 (GMT) to 4 August 24:00 (GMT) 2012 were acquired from the Guangzhou Meteorological Bureau.The first two days were sunny and high-temperature weather, whereas rainy and cloudy weather was observed on the last two days.The observational data were gathered from four weather stations (one in a suburban area, and three in urban areas) in Domain 3. The suburban weather station (suburban point in Figure 1) is located on an elevation in the mountainous and hilly terrain far away from the downtown of Guangzhou, such that there are no tall buildings and medium-sized enterprises within 3 km, and the station is not shaded.The other three weather stations (Points 1-3 in Figure 1) are located in the urban area.The area around Point 1 experiences rapid urbanization, and thus its representation greatly differs between the MODIS and extracted land use map.Points 2 and 3 are located in the old down town.Points 1 and 2 were representative points for the model evaluation.Hourly temperature and wind velocity data were collected to investigate the UHI characteristics and the performance of the WRF/UCM model with the extracted land use data.

Air Temperature at 2-m Height
To investigate the performance of the WRF/UCM model with the extracted land use data, we compared the results of the three experiments with the 2-m temperature data observed at the

Observational Data
The simulation lasts for one month, as obtaining long time data to evaluate performance of the WRF/UCM and the land use data is better.Only four days of climate data from 1 August 0:00 (GMT) to 4 August 24:00 (GMT) 2012 were acquired from the Guangzhou Meteorological Bureau.The first two days were sunny and high-temperature weather, whereas rainy and cloudy weather was observed on the last two days.The observational data were gathered from four weather stations (one in a suburban area, and three in urban areas) in Domain 3. The suburban weather station (suburban point in Figure 1) is located on an elevation in the mountainous and hilly terrain far away from the downtown of Guangzhou, such that there are no tall buildings and medium-sized enterprises within 3 km, and the station is not shaded.The other three weather stations (Points 1-3 in Figure 1) are located in the urban area.The area around Point 1 experiences rapid urbanization, and thus its representation greatly differs between the MODIS and extracted land use map.Points 2 and 3 are located in the old down town.Points 1 and 2 were representative points for the model evaluation.Hourly temperature and wind velocity data were collected to investigate the UHI characteristics and the performance of the WRF/UCM model with the extracted land use data.

Air Temperature at 2-m Height
To investigate the performance of the WRF/UCM model with the extracted land use data, we compared the results of the three experiments with the 2-m temperature data observed at the automatic weather stations.Figure 3a presents the observed and modeled air temperature at 2-m height lasting 20 days from 20 July to 9 August in the suburban area, and the four-day comparison from 1 August to 4 August was analyzed in the urban areas (Figure 3b).
From Figure 3a, it can be found that the modeled daily temperature cycles are generally consistent with the observational data on sunny and hot days, but vary in magnitude on rainy days.The same characteristic can be found in Figure 3b: all models performed better on 1 and 2 August than on 3 and 4 August because of the rainy and cloudy weather seen on the latter two days.All of the experimental schemes showed the diurnal variation of temperature in urban and suburban areas, though with varying magnitudes.The air temperature from the MODIS experiment was lower than the observed data during the period 0:00-12:00 at the suburban station, but the RS_12 and UCM_12 experiments gave lower temperatures from 06:00-12:00.In Figure 3b, all the modeled air temperatures were higher than the observed data, and there was a good fit for the maximum air temperature between the three experiments and the observed data for 1 and 2 August.
Sustainability 2016, 8, 628 6 of 15 automatic weather stations.Figure 3a presents the observed and modeled air temperature at 2-m height lasting 20 days from 20 July to 9 August in the suburban area, and the four-day comparison from 1 August to 4 August was analyzed in the urban areas (Figure 3b).
From Figure 3a, it can be found that the modeled daily temperature cycles are generally consistent with the observational data on sunny and hot days, but vary in magnitude on rainy days.The same characteristic can be found in Figure 3b: all models performed better on 1 and 2 August than on 3 and 4 August because of the rainy and cloudy weather seen on the latter two days.All of the experimental schemes showed the diurnal variation of temperature in urban and suburban areas, though with varying magnitudes.The air temperature from the MODIS experiment was lower than the observed data during the period 0:00-12:00 at the suburban station, but the RS_12 and UCM_12 experiments gave lower temperatures from 06:00-12:00.In Figure 3b, all the modeled air temperatures were higher than the observed data, and there was a good fit for the maximum air temperature between the three experiments and the observed data for 1 and 2 August.The simulated and observed daily mean (T2M), maximum (T2MAX), and minimum (T2MIN) 2-m air temperatures are shown in Table 3.In general, the UCM_12 experiment efficiently simulated the daily variations in near-surface air temperature, and the RS_12 results were better than the MODIS on the near-surface air temperature simulation results.The simulated and observed daily mean (T2M), maximum (T2MAX), and minimum (T2MIN) 2-m air temperatures are shown in Table 3.In general, the UCM_12 experiment efficiently simulated the daily variations in near-surface air temperature, and the RS_12 results were better than the MODIS on the near-surface air temperature simulation results.Willmott suggested using the root mean square error (RMSE), systematic (RMSE S ) and unsystematic (RMSE U ) errors, and an index of agreement d to evaluate the simulation performance [28].The index of agreement d can be interpreted as a measure of how error-free a model predicts a variable.The RMSE summarizes the magnitude of the average difference between observation and prediction.Systematic errors, described by RMSE S , result from causes which occur consistently.RMSE U gives information about unsystematic errors, which consist of a number of small effects are positive and some are negative in terms of affecting the final output value.The good model therefore has a systematic difference of zero while the unsystematic difference should approach the RMSE.The d value of 1.0 indicates perfect agreement between observation and prediction.The d can be calculated by Equation (4) [28,29]. where O is the average of the observational data.The index of agreement d describes how accurately a model predicts a variable.When d value is equal to 1.0, it indicates perfect prediction for the variable.
The difference measures for the four comparison points 2-m temperature are summarized in Table 4.The RMSE for the three simulations at the four weather stations ranged from 1.41 to 2.86 ˝C.Additionally, the systematic errors (RMSE S ) were comparatively small and the unsystematic errors (RMSE U ) approached the RMSE in UCM_12, indicating that the UCM_12 experiment fits the criteria of the systematic error.The d values for all of the simulations at the four stations ranged from 0.74 to 0.95, suggesting that all three experiments provided a reasonable approximation of the observed data.Comparing the d values for the three simulations for each point, the value at UCM_12 is highest, followed by the RS_12 d value, which is in turn larger than the MODIS d value.This suggests that the extracted land use data used in WRF/UCM model provide the best reproduction of temperature variations among the three simulations in both suburban and urban areas.Additionally, the results show that the extracted land use data provided a better simulation than the MODIS data.

Urban Heat Island Intensity (UHII)
By comparing the results of the three numerical experiments, their ability to simulate the urban heat island effect can be examined.Figure 4 presents the evolution of the urban heat island intensity (UHII) as simulated by the three experiments.The intensity of the heat island is calculated as the difference of temperature between urban Point 1 and the suburban point at 2-m height.As shown by the observations, the UHII gradually increases after midday, reaching a maximum at night, with a maximum temperature difference of more than 5 ˝C.The UHII then gradually decreases in the morning, even becoming negative at noon.
the other experiments at other times.
The maximum UHII and the average local UHII reflect the average difference over the whole day, and were selected to analyze the UHI effect.
Modeled and observed values of maximum UHII at Point 1 and Point 2 are compared in Figure 5a,b, respectively.At Point 1, the RS_12 and UCM_12 results have better correspondence to the observations than the MODIS results.The UCM_12 experiment did not show a noticeably better performance among the simulations despite the implementation of the coupled WRF/UCM model.At Point 2, although the UCM_12 experiment showed a higher maximum UHII, there is still a large deviation from the observations.On 3 August, none of the experiments agreed with the observations, predominantly because of the rainy weather.The average local UHII were calculate by Equation ( 5), where, T u stands for the hourly 2-m height temperature in the urban area, T r stands for the hourly 2-m height temperature in the suburban area, t 0 stands for initial time of the simulation, and t 1 stands for end time of the simulation.All experiments plausibly reproduce the characteristics of UHII evolution, except for on 3 August, when the negative UHII is weakly reproduced and its overall variation is underestimated.UCM_12 produces a better reproduction of the UHII at peak values, but a poorer reproduction than the other experiments at other times.
The maximum UHII and the average local UHII reflect the average difference over the whole day, and were selected to analyze the UHI effect.
Modeled and observed values of maximum UHII at Point 1 and Point 2 are compared in Figure 5a,b, respectively.At Point 1, the RS_12 and UCM_12 results have better correspondence to the observations than the MODIS results.The UCM_12 experiment did not show a noticeably better performance among the simulations despite the implementation of the coupled WRF/UCM model.At Point 2, although the UCM_12 experiment showed a higher maximum UHII, there is still a large deviation from the observations.On 3 August, none of the experiments agreed with the observations, predominantly because of the rainy weather.By comparing the results of the three numerical experiments, their ability to simulate the urban heat island effect can be examined.Figure 4 presents the evolution of the urban heat island intensity (UHII) as simulated by the three experiments.The intensity of the heat island is calculated as the difference of temperature between urban Point 1 and the suburban point at 2-m height.As shown by the observations, the UHII gradually increases after midday, reaching a maximum at night, with a maximum temperature difference of more than 5 °C .The UHII then gradually decreases in the morning, even becoming negative at noon.
All experiments plausibly reproduce the characteristics of UHII evolution, except for on 3 August, when the negative UHII is weakly reproduced and its overall variation is underestimated.UCM_12 produces a better reproduction of the UHII at peak values, but a poorer reproduction than the other experiments at other times.
The maximum UHII and the average local UHII reflect the average difference over the whole day, and were selected to analyze the UHI effect.
Modeled and observed values of maximum UHII at Point 1 and Point 2 are compared in Figure 5a,b, respectively.At Point 1, the RS_12 and UCM_12 results have better correspondence to the observations than the MODIS results.The UCM_12 experiment did not show a noticeably better performance among the simulations despite the implementation of the coupled WRF/UCM model.At Point 2, although the UCM_12 experiment showed a higher maximum UHII, there is still a large deviation from the observations.On 3 August, none of the experiments agreed with the observations, predominantly because of the rainy weather.The average local UHII were calculate by Equation ( 5), where, T u stands for the hourly 2-m height temperature in the urban area, T r stands for the hourly 2-m height temperature in the suburban area, t 0 stands for initial time of the simulation, and t 1 stands for end time of the simulation.The average local UHII were calculate by Equation ( 5), where, T u stands for the hourly 2-m height temperature in the urban area, T r stands for the hourly 2-m height temperature in the suburban area, t 0 stands for initial time of the simulation, and t 1 stands for end time of the simulation.The average local UHII for all experiments and the observation data at Point 1 and Point 2 are shown in Figure 6a,b.The daytime was constrained between 8:00 and 20:00 and the rest was considered as nighttime.All experiments have a better agreement with the observations and represent larger magnitudes of average local UHII at Point 2 than at Point 1.All experiments underestimate both day-and nighttime UHII.In urban areas, the UCM_12 experiment produces poor statistical scores for the average daytime local UHII, but better ones for the average nighttime local UHII.The UCM_12 experiment produces the highest average nighttime local UHII, although the value is still lower than the observed one.The average local UHII for all experiments and the observation data at Point 1 and Point 2 are shown in Figure 6a,b.The daytime was constrained between 8:00 and 20:00 and the rest was considered as nighttime.All experiments have a better agreement with the observations and represent larger magnitudes of average local UHII at Point 2 than at Point 1.All experiments underestimate both day-and nighttime UHII.In urban areas, the UCM_12 experiment produces poor statistical scores for the average daytime local UHII, but better ones for the average nighttime local UHII.The UCM_12 experiment produces the highest average nighttime local UHII, although the value is still lower than the observed one.Although the results of the coupled WRF/UCM model do not accurately reproduce the observations and reveal consistently underestimated UHII, they feasibly reflect the diurnal variation of the UHII.Additionally, the coupled WRF/UCM model exhibits better performance with the extracted land use data than with the MODIS land use data.

Wind Velocity at 10-m Height
Regarding wind velocity, all of the simulation results were close to the observation data in the daytime.The observed and simulated wind velocity at 10-m height at the suburban point and Point 1 are shown in Figure 7a,b, respectively.In general, wind velocity simulated in all three experiments is overestimated, and the distribution trends of all the experiments are similar when comparing the simulations and observations.For the urban area, both modeled and observed wind velocity is lower in the nighttime, which corresponds to the changes in UHI intensity.It is believed that the strong thermal stability of the surface layer is the primary reason for nighttime enhancement of UHI.It should also be noted that surface wind speed is generally reduced in a stable surface layer.All experiments have large magnitudes compared with the observations.Even using the coupled WRF/UCM model, the influence of the buildings on the wind environment cannot be reproduced.Although the results of the coupled WRF/UCM model do not accurately reproduce the observations and reveal consistently underestimated UHII, they feasibly reflect the diurnal variation of the UHII.Additionally, the coupled WRF/UCM model exhibits better performance with the extracted land use data than with the MODIS land use data.

Wind Velocity at 10-m Height
Regarding wind velocity, all of the simulation results were close to the observation data in the daytime.The observed and simulated wind velocity at 10-m height at the suburban point and Point 1 are shown in Figure 7a,b, respectively.In general, wind velocity simulated in all three experiments is overestimated, and the distribution trends of all the experiments are similar when comparing the simulations and observations.For the urban area, both modeled and observed wind velocity is lower in the nighttime, which corresponds to the changes in UHI intensity.It is believed that the strong thermal stability of the surface layer is the primary reason for nighttime enhancement of UHI.It should also be noted that surface wind speed is generally reduced in a stable surface layer.All experiments have large magnitudes compared with the observations.Even using the coupled WRF/UCM model, the influence of the buildings on the wind environment cannot be reproduced.The average local UHII for all experiments and the observation data at Point 1 and Point 2 are shown in Figure 6a,b.The daytime was constrained between 8:00 and 20:00 and the rest was considered as nighttime.All experiments have a better agreement with the observations and represent larger magnitudes of average local UHII at Point 2 than at Point 1.All experiments underestimate both day-and nighttime UHII.In urban areas, the UCM_12 experiment produces poor statistical scores for the average daytime local UHII, but better ones for the average nighttime local UHII.The UCM_12 experiment produces the highest average nighttime local UHII, although the value is still lower than the observed one.Although the results of the coupled WRF/UCM model do not accurately reproduce the observations and reveal consistently underestimated UHII, they feasibly reflect the diurnal variation of the UHII.Additionally, the coupled WRF/UCM model exhibits better performance with the extracted land use data than with the MODIS land use data.

Wind Velocity at 10-m Height
Regarding wind velocity, all of the simulation results were close to the observation data in the daytime.The observed and simulated wind velocity at 10-m height at the suburban point and Point 1 are shown in Figure 7a,b, respectively.In general, wind velocity simulated in all three experiments is overestimated, and the distribution trends of all the experiments are similar when comparing the simulations and observations.For the urban area, both modeled and observed wind velocity is lower in the nighttime, which corresponds to the changes in UHI intensity.It is believed that the strong thermal stability of the surface layer is the primary reason for nighttime enhancement of UHI.It should also be noted that surface wind speed is generally reduced in a stable surface layer.All experiments have large magnitudes compared with the observations.Even using the coupled WRF/UCM model, the influence of the buildings on the wind environment cannot be reproduced.The difference measures of the modeled wind velocity at 10-m height of four comparison points are summarized in Table 5.The d values for all experiments at the four stations ranged from 0.34 to 0.70, indicating low consistency of simulations with the observed data.At the same time, the systematic errors ( RMSE S ) were comparatively small and the unsystematic errors ( RMSE U ) approached the RMSE in UCM_12, indicating better performance of coupled WRF/UCM with extracted land use data compared to other cases, although the improvement is relatively small.The difference measures of the modeled wind velocity at 10-m height of four comparison points are summarized in Table 5.The d values for all experiments at the four stations ranged from 0.34 to 0.70, indicating low consistency of simulations with the observed data.At the same time, the systematic errors (RMSE S ) were comparatively small and the unsystematic errors (RMSE U ) approached the RMSE in UCM_12, indicating better performance of coupled WRF/UCM with extracted land use data compared to other cases, although the improvement is relatively small.

Energy Balance Data
The UHI is mainly determined by the surface energy balance [30], and the contributing factors during land surface energy equilibrium processes on the UHI can be represented in models.To focus on the formation of the UHI on 1 August and to investigate the performance of the UCM in simulating the UHI, the energy equilibrium processes between urban and suburban areas, and between experiments UCM_12 and RS_12, were compared and analyzed.Figure 8 shows the temporal variability of parameters relevant to the land surface energy equilibrium in the urban area during the UCM_12 experiment from 1 August 00:00 (GMT) to 23:00 (GMT).The absorption of solar short-wave radiation is predominantly an energy input, whereas the outputs are latent heat flux and ground heat flux.Long-wave radiation shows weak daily variation (approximately 450 W/m 2 ).Ground heat flux and upward heat flux reach their maximum values of 250 and 340 W/m 2 , respectively, in the daytime around noon.The latent heat flux only reaches a maximum of 100 W/m 2 owing to evaporation, which is absent at night.At night, the ground heat flux acts as an input, becoming a major source of energy, while the upward heat flux, though much smaller than during the daytime, remains positive as it continues to transfer heat to the atmosphere.The characteristics of these energy equilibrium processes suggest a positive performance when simulating the UHI using the WRF/UCM with the land use data extracted from remote sensing data.
Figure 9 presents the differences in the energy equilibrium between Point 1 and the suburban point simulated in the UCM_12 experiment.Urban areas absorb more short-wave radiation than suburban areas, although the difference is small.The diurnal latent heat flux was significantly different between the urban and suburban areas.During the day, the latent heat flux was reduced in urban areas because of the limited vegetation and less evaporation, while it remained small at night.During the day, the upward heat flux was significantly higher in the urban area than in the suburban areas, close to 200 W/m 2 .After sunset, the difference in the ground heat flux turns from negative to positive, while that of upward heat flux remained positive.The positive heat fluxes, though small, resulted in a higher temperature at night in urban areas than in suburban areas.The WRF model simulation of UHI with the land use data extracted from remote sensing data can reveal the causes of the UHI.  Figure 10 shows the differences between the individual terms of energy equilibrium between the UCM_12 and RS_12 experiments.The improvement of the simulation performance by the UCM model can be illustrated from the difference.The urban latent heat flux is reduced in the daytime when coupled UCM the model in UCM_12, the reduction reaches a maximum of 100 W/m 2 .In addition, the UCM_12 simulation suggested a larger diurnally upward heat flux and ground heat flux than the RS_12 simulation.This suggests that the urban canopy model can be used to improve nighttime UHI simulations.

Summary and Conclusions
The thermal environment over Guangzhou from late July to early August 2012 was investigated using the coupled WRF/UCM model.Simulation experiments were performed using both the urban land use data extracted from remote sensing datasets and the default MODIS data.

Summary and Conclusions
The thermal environment over Guangzhou from late July to early August 2012 was investigated using the coupled WRF/UCM model.Simulation experiments were performed using both the urban land use data extracted from remote sensing datasets and the default MODIS data.
(1) Using the new land use data, the simulated 20 days of daily temperature cycles at the suburban weather station and four days of daily temperature cycles at urban point showed an encouraging agreement with observations, especially on hot sunny days.UCM_12 simulated maximum diurnal temperatures closer to the observed values than the simulations that excluded the UCM model.The coupled WRF/UCM with the new land use data improved the simulation performance.(2) Compared with the maximum UHII and the average local UHII, all of the simulated results reproduced the diurnal characteristics of UHI intensity.Both RS_12 and UCM_12 simulations performed better than the default geographic model, although they did not perfectly replicate the observations.(3) The modeled wind velocity results were higher than observations.The new land use data produced similar results as the default land use data.For the urban area, all wind velocity simulation results and observations were lower, and the UHII remains high during the nighttime.(4) The UCM_12 experiment successfully reproduced the differences in the energy equilibrium both in the urban and suburban areas.In the urban area, most of the input energy is used in sensible heating, which resulting the higher temperature, can be simulated by the WRF/UCM.The latent heat flux in the urban area is lower than in suburban area because lack of evapotranspiration of water vapor in the urban area.For the nocturnal situation, relatively high minimum temperature is reproduced in the urban owing to the sustained upward ground heat flux heats the atmosphere in the urban area.
As changing the land use data can have a significant impact on the WRF/UCM simulation result, this study confirmed that the WRF/UCM model with land use data extracted from a remote sensing dataset can be an effective tool for urban thermal environment analyses.

Figure 1 .
Figure 1.The three nested domains used for the numerical simulation and the location of comparison points.

Figure 1 .
Figure 1.The three nested domains used for the numerical simulation and the location of comparison points.

Landsat- 7
with HSI and coupled with the UCM model.The physics schemes remained the same for all runs.

Figure 2 .
Figure 2. Types of land use data in Domain 3: (a) UCM_12, enhancement of the extracted land use data from Landsat-7 with HIS; and (b) Modis, default land use data in WRF.

Figure 2 .
Figure 2. Types of land use data in Domain 3: (a) UCM_12, enhancement of the extracted land use data from Landsat-7 with HIS; and (b) Modis, default land use data in WRF.

Figure 3 .
Figure 3.Comparison of the air temperature at 2-m height simulated in the three experiments with the observed temperature (OBS): (a) Suburban station; and (b) Point 1.

Figure 3 .
Figure 3.Comparison of the air temperature at 2-m height simulated in the three experiments with the observed temperature (OBS): (a) Suburban station; and (b) Point 1.

Figure 4 .Figure 5 .
Figure 4. Heat island intensity of the different experiments at Point 1.

Figure 4 .
Figure 4. Heat island intensity of the different experiments at Point 1.

Figure 4 .Figure 5 .
Figure 4. Heat island intensity of the different experiments at Point 1.

Figure 7 .
Figure 7.Comparison of the wind velocity at 10-m height simulated in the three experiments with observations: (a) Suburban station; and (b) Point 1.

1 -Figure 7 .
Figure 7.Comparison of the wind velocity at 10-m height simulated in the three experiments with observations: (a) Suburban station; and (b) Point 1.

Sustainability 2016, 8 , 628 11 of 15 Figure 8 .Figure 8 .
Figure 8. Temporal evolution of parameters relevant to the land surface energy equilibrium in the urban area during the UCM_12 experiment on 1 August.The absorption of solar short-wave radiation is predominantly an energy input, whereas the outputs are latent heat flux and ground heat flux.Long-wave radiation shows weak daily variation (approximately 450 W/m 2 ).Ground heat flux and upward heat flux reach their maximum values of 250 and 340 W/m 2 , respectively, in the daytime around noon.The latent heat flux only reaches a maximum of 100 W/m 2 owing to evaporation, which is absent at night.At night, the ground heat

Figure 9 .Figure 9 .
Figure 9. Differences in the energy equilibrium between Point 1 and the suburban point for the UCM_12 simulation.

Figure 9 .
Figure 9. Differences in the energy equilibrium between Point 1 and the suburban point for the UCM_12 simulation.

Figure 10 .
Figure 10.Differences in the energy equilibrium between the UCM_12 and RS_12 simulation.

( 1 )Figure 10 .
Figure 10.Differences in the energy equilibrium between the UCM_12 and RS_12 simulation.

Table 3 .
Comparison of the simulated and observed 2-m temperatures (T2).
*: The average results over 3 urban points.

Table 3 .
Comparison of the simulated and observed 2-m temperatures (T2).
*: The average results over 3 urban points.

Table 4 .
Quantitative evaluation of three cases' performance on temperature compared with observed data.

Table 5 .
Quantitative evaluation of three cases' performance on wind velocity compared with observed data.