Effects of anthropogenic heat due to air-conditioning systems on an extreme high temperature event in Hong Kong

Anthropogenic heat flux is the heat generated by human activities in the urban canopy layer, which is considered the main contributor to the urban heat island (UHI). The UHI can in turn increase the use and energy consumption of air-conditioning systems. In this study, two effective methods for water-cooling air-conditioning systems in non-domestic areas, including the direct cooling system and central piped cooling towers (CPCTs), are physically based, parameterized, and implemented in a weather research and forecasting model at the city scale of Hong Kong. An extreme high temperature event (June 23–28, 2016) in the urban areas was examined, and we assessed the effects on the surface thermal environment, the interaction of sea–land breeze circulation and urban heat island circulation, boundary layer dynamics, and a possible reduction of energy consumption. The results showed that both water-cooled air-conditioning systems could reduce the 2 m air temperature by around 0.5 °C–0.8 °C during the daytime, and around 1.5 °C around 7:00–8:00 pm when the planetary boundary layer (PBL) height was confined to a few hundred meters. The CPCT contributed around 80%–90% latent heat flux and significantly increased the water vapor mixing ratio in the atmosphere by around 0.29 g kg−1 on average. The implementation of the two alternative air-conditioning systems could modify the heat and momentum of turbulence, which inhibited the evolution of the PBL height (a reduction of 100–150 m), reduced the vertical mixing, presented lower horizontal wind speed and buoyant production of turbulent kinetic energy, and reduced the strength of sea breeze and UHI circulation, which in turn affected the removal of air pollutants. Moreover, the two alternative air-conditioning systems could significantly reduce the energy consumption by around 30% during extreme high temperature events. The results of this study suggest potential UHI mitigation strategies and can be extended to other megacities to enable them to be more resilient to UHI effects.


Introduction
The substitution of natural surfaces with impervious urban structures because of rapid ongoing urbanization and human activities has resulted in changes in the urban surface energy and water balance, and has influenced the local atmospheric boundary layer structure and urban climate (Oke 1976, 1988, Bonacquisti et al 2006. Almost all the energy consumption for human activities can eventually transform into anthropogenic heat and be released into the atmosphere within the Earth's land-atmosphere system (Taha 1997, Flanner 2009, Xie et al 2016. Anthropogenic heat fluxes can increase turbulence fluxes, including both sensible and the latent heat fluxes (Oke 1988). Moreover, anthropogenic heat fluxes influence surface meteorological conditions, particularly air temperature and water vapor pressure; the vertical motion of urban air flow, which in turn affects urban heat island circulation (UHIC) and dynamics and thermal dynamics of the boundary layer; and air pollutants within the urban canopy layer (Oke 1987, 1988, Ichinose et al 1999, Kikegawa et al 2003, Sailor and Lu 2004, Feng et al 2012, Salamanca et al 2014, Xie et al 2016, Ma et al 2017.
Buildings are found to be the major source of anthropogenic heat fluxes; residential and commercial buildings account for 40% of total energy consumption (EIA 2012, Quah andRoth 2012). In the global-scale urban consumption of energy model, heat release from buildings is the largest contributor (89%-96%) of heat emission globally (Allen et al 2011). The energy consumption of buildings is likely to increase (Mansur et al 2008). Studies on the effects of air-conditioning systems on the urban climate were only found for the cities of Houston, Wuhan, Tokyo, Madrid, and Paris (Ohashi et al 2007, Wen and Lian 2009, Salamanca et al 2011, de Munck et al 2013. These studies focused mainly on the effect of air-conditioning systems on the outdoor air temperature by using models with different levels of sophistication, did not directly parameterize air-conditioning systems into urban canopy models and neglected the effects of the main physically-based processes on the atmosphere. Wen and Lian (2009) revealed that the domestic use of air conditioners increases the mean air temperature by 0.2 • C or 2.56 • C under neutral and inversion conditions, respectively. Further, de Munck et al (2013) found that the temperature increase reached 0.5 • C and 2 • C due to the current and double waste heat release into the atmosphere, respectively. However, the simplified model failed to capture the effects of sensible and latent heat fluxes within the urban canopy layer. Estimation of anthropogenic sensible and latent heat fluxes is difficult, and indirect measurement, empirical estimation, and numerical models are mostly used for this (Grimmond 1992, Kikegawa et al 2003, Sailor and Lu 2004, Pigeon et al 2007. The city of Hong Kong is within a humid subtropical climate according to Köppen classification. Electricity consumption in non-domestic buildings accounted for 61.8% of total consumption in June 2016 (www.censtatd.gov.hk). Statistical analysis shows that air-conditioning accounts for around 30% of total electricity consumption, of which 68% is in nondomestic premises (Ho et al 2007). The Hong Kong government encourages the adoption of water-cooled air-conditioning systems, including direct cooling systems (DCS) and central piped cooling towers (CPCT), which seems to be an effective method of territory-wide energy improvement and reduction of electricity consumption and greenhouse gas emissions from 2020 onward (Ho et al 2007). The effects of the altered air-conditioning systems on the urban thermal and dynamic environment at the city scale are rarely studied but essential. A central component of our study is whether the implementation of two alternative air-conditioning systems can result in an improvement of urban resilience to UHI effects.
Extreme high temperature events and abnormally hot periods are among the most significant climatic stressors for public health, ecosystems, economies, societies, energy consumption, UHI threat, and more (De Bono et al 2004, Sharma et al 2016, Zampieri et al 2016. June 2016 had the second highest temperature in Hong Kong since 1884 (www.weather.gov.hk/) and was the warmest June on record (SCMP, www.scmp.com/tech/scienceresearch/article/1992366/june-2016-warmest-monthrecord-globally). Therefore, such situations could represent the worst conditions for cooling-energy consumption.
To help overcome the aforementioned gaps in research on the influence of water-cooled airconditioning systems on climate, in this study, we assess the potential effects of DCS and CPCT systems on an urban scale. Changes in urban air temperature, moisture level, energy fluxes, urban boundary layer structure, local circulations, and potential energy consumption improvement are assessed. Although the results of this study are location dependent, the effectiveness of alternative systems is expected to make other urbanized areas more resilient.

Overview of the model
The mesoscale meteorological model has been used to study various meteorological phenomena. The advanced research weather research and forecasting (WRF) model (ARW version 3.6.1) is a non-hydrostatic, compressible model with a mass coordinate system (Skamarock et al 2008). The model was set up with four nested domains with 79 × 79, 79 × 79, 118 × 118 and 142 × 136 grid points possessing spatial resolutions of 13.5, 4.5, 1.5 and 0.5 km, respectively (figure 1(a)). The innermost domain covered Hong Kong entirely (figure 1) and provided sufficient information on UHI effects and sea-land breeze circulation (SLBC). The vertical grid contains 51 terrains following full sigma levels from the ground up to 50 hPa. The planetary boundary layer (PBL) scheme of Bougeault and Lacarrere (1989) was used. The other physical schemes used in this study were the singlemoment three-class microphysics scheme (Hong et al 2004), Dudhia shortwave radiation (Dudhia 1989), RRTM long-wave radiation (Mlawer et al 1997), the Noah land surface model (Chen and Dudhia 2001), and Kain-Fritsch (Kain and Fritsch 1990) cumulus schemes for the two outer domains. A 120 h simulation (from 0000 UTC June 23 to 0000 UTC June 28, 2016) was conducted with the initial and boundary conditions from the National Centers for Environmental Prediction (NCEP) Global Forecast System Final Analyses data (FNL) at 1 • at 6 h intervals and daily updated sea surface temperatures from NCEP Marine Modeling and Analysis Branch data at 0.5 • .
Hong Kong has a large spatial and vertical land surface heterogeneity ( figure 1(b)). To better represent the dynamics and thermal interacts in the urban canopy layers, multilayer building effect parameterization (BEP, Martilli et al 2002) and the building energy model (BEM, Salamanca et al 2010) were chosen. The multilayer canopy model, BEP/BEM, with a revised drag coefficient, based on plan area fraction, p , has been tested over Kowloon Peninsula in Hong Kong, and has shown reasonably good results for capturing spatial air temperature distribution, diurnal wind rotation, and stagnation wind phenomena in the lower level convergence zone due to the combination of urban heat island circulation and sea-land breeze circulation (Wang et al 2017). The latest gridded building information and land cover and land use data (figure 1(c)) were used for the entire Hong Kong area for better representation of spatial heterogeneity. The default airconditioning system in BEM is the wall-type scenario (shown in figure 2) in all urban areas. In this study the wall-type air-conditioning system is assumed to be used only in residential areas and central air-cooled air systems in non-domestic areas according to the real situation in Hong Kong (shown in figure 2). To study the effect of air-conditioning systems in non-domestic areas on the urban climate in the city of Hong Kong, several configurations or scenarios of air-conditioning systems (shown in figure 2) in non-domestic areas were simulated, as summarized in table 1.

Air-conditioning systems formulation
The new cooling tower scheme (CPCT case) coupled with the building energy model (BEM) is illustrated. In the counter-flow cooling tower system, the gas phase (cool air) flows upwards and the liquid phase (water) downwards (figure 2) and a large interface exists between the two phases (Milosavljevic and Heikkilä 2001). Heat and mass transfer takes place in the cooling tower system, and the water droplets carry out evaporative cooling with the dry and cool air from the atmosphere (Gutiérrez et al 2015). After the process, the water temperature decreases, but the air increases both the temperature and humidity and is released into the atmosphere; further, make-up water is needed for the cooling tower system to maintain the process.
The heat exchange effectiveness, CPCT , between the refrigerant and the water in the condenser is the heat that is generated inside a building at each floor, H , which is calculated via the BEM model, and the heat that is pumped into the cooling tower system, H CPCT . Assume that no heat is lost in the pipe between the condenser and the cooling tower. Then, the heat exchange effectiveness, CPCT , can be defined as follows:   where m represents the total number of floors in the building.
If we assume that the water temperature difference between the inlet and the outlet of the condenser system is constant, then the water volume flow rate is as follows: where denotes the water density, C represents the water heat capacity, and ΔT stands for the water temperature difference between the inlet and the outlet of the condenser and the cooling tower.
In the counter-flow cooling tower system, the air volume flow rate, q , is as follows: where v denotes the ratio of the water mass flow rate and the air flow rate in the contour-flow cooling tower system, which is determined by the inter-facial contact area and the heat-transfer coefficient, and represents the air density. Here, assume that the ratio v is 0.7 . Then, the mass flow rate of the air in the cooling tower system is defined as follows: = . (4) The inlet air in the cooling tower is assumed to be from the atmosphere at the appropriate building height level; thus, the enthalpy of the outlet air, h ao , from the cooling tower system can be defined as follows:  where h ai denotes the enthalpy of the inlet air and can be defined as follows: where C denotes the specific heat of air at constant pressure, C pw represents the specific heat of water at constant pressure, L indicates the latent heat of vaporization, and T ai and d ai refer to the air temperature and the water vapor mixing ratio of the inlet air, respectively. If the outlet air is saturated, then the water vapor mixing ratio of the outlet air from the cooling tower is as follows: = 6.11 × 10 7.5 × ao 237.3+ ao where e denotes the water vapor pressure calculated using the Magnus (1844) empirical formula. Direct cooling systems (DCS) use sea water to cool the condenser in the buildings (Ho et al 2007). If we assume that the sea water volume flow rate, q sw , is constant, then the temperature difference between the inlet and outlet of the system is as follows: where sw denotes the sea water density, C sw represents the sea water heat capacity. For the DCS systems in this study, the estimation of the temperature rises of the sea water during one day is around 1 × 10 −5 • C, which is ignored in this study.

Air-conditioning heat partition in BEM
BEM starts from the assumption that the internal target temperature and humidity levels remain within the comfort range  and then estimates the cooling load inside the buildings. The default air-conditioning system in BEM is the wall-type scenario (shown in figure 2), and is assumed to be used in residential areas in this study. The BEM model calculates the sensible heat fluxes, out , and latent heat fluxes, LE out , for each floor, i, into the atmosphere required to maintain standard indoor comfort conditions . The details of the sensible and latent heat flux calculation can be found in Salamanca et al (2010). Then, the sensible and latent heat fluxes generated at each floor, i, and released into the atmosphere through the air-conditioning systems  can be expressed as follows: where COP denotes the coefficient of performance for the air-conditioning systems, which is assumed to be constant for all wall-type systems. The Baseline case (figure 2) assumes that the systems are settled at the rooftop level. Then, the sensible and latent heat fluxes generated at each floor, i, are summed up and released into the atmosphere only at the rooftop level through the air-conditioning systems (followed by Krpo et al 2010); these fluxes can be expressed as follows: The sensible or latent heat fluxes at the rooftop are given by the sum of turbulent fluxes H t ur or LE tur , from the roof surfaces, which is computed using the Louis formulation (Martilli et al 2002) and Base or LE Base , followed by the method proposed by Krpo et al (2010).
For the DCS case (figure 2), the sensible heat fluxes generated at each floor, i, are summed up and go back to the sea water through DCS systems, which is assumed to be located in the equipment unit in the building (ASHRAE Handbook 2008); this can be expressed as follows: For the CPCT case (figure 2), the sensible, CPCT , and latent heat fluxes, CPCT , to the atmosphere at the rooftop level are shown below: (17) Note that for the residential areas, wall-type systems are kept unchanged for all the cases.

Model evaluation
To validate the performance of the Baseline model, the simulation results are compared with surface level observations. The model could reasonably reproduce the physical characteristics of the temperature and wind field distribution (shown in figures 3, 4 and S1 available at stacks.iop.org/ERL/13/034015/mmedia). Figure 3 shows that there are some convergence zones (Tong et al 2005) formed at 2:00 pm on 25 June in different areas in Hong Kong. The temperature distribution is basically captured, although the model underestimated the air temperature at some stations, like SEK on 24 (shown in figure S1) and 25 June, which may due to the uncertainty of soil conditions and also the 1 • NCEP FNL data chosen as the boundary condition (Takane and Kusaka 2011). The Baseline model could well capture the high temperature area in the convergence zone in Kowloon Peninsula at 3:00 pm on 24 June (figure S1), which is dominated by a westerly sea breeze and easterly wind and 2:00 pm on 25 June (figure 3), which is dominated by a southeastern wind. The observed and simulated 2 m air temperature for all urban and rural sites are compared (shown in figure 4). Note that because the horizontal spatial resolution is limited to 0.5 km, the stations along the coast (CCH, WGL, PEN, and SLW, as shown in figure 1(d)) are excluded because they may represent sea surface characteristics due to the land cover and land use classification. The Baseline model could capture the spatial distribution of 2 m air temperature well. The agreement between the Baseline model and observations was set in accordance with the criteria value of mean bias error (MBE), root-mean-square error (RMSE), and index of agreement (IOA) for 2 m air temperature for 27 available stations, which are 0.224, 1.242 and 0.90 respectively. The MBE value between the Baseline model and observations for the urban areas for 2 m air temperature is −0.09 • C. The Baseline model could capture the wind rotation from southwestern to southeastern on 24 June well, particularly in the central area of the Kowloon Peninsula, as shown in figure  5. The model underestimated the wind field when the extreme high-temperature event ended on 28 June, as shown in figure 5. In the Baseline case, the anthropogenic heat releases at the rooftop level and can affect the street level by vertical diffusion (Martilli et al 2002, which is not significant in a high-rise city. In general, the Baseline case could capture the temperature distribution and wind rotation in the urban areas well, which can fulfill the goal of this study.

Effects of air-conditioning systems on air temperature, humidity, and surface energy balance
To evaluate the contribution of air-conditioning systems on temperature and humidity levels, the 2 m air temperature and the water vapor mixing ratio are chosen as the criteria for comparison. Figure S2 shows the 2 m air temperature distribution at 2:00 pm on June 25, 2016, i.e. the time when the maximum air temperature was recorded, and the temperature difference between the Baseline case and different air-conditioning    Further, the large depth evolution of the PBL height during the daytime inhibited the effects of air-conditioning on the urban air temperature as compared to the relatively shallow PBL height and stable atmospheric conditions during the night. If the anthropogenic heat were excluded in the urban areas, the temperature difference could have reached around 4 • C in the commercial areas in Kowloon Peninsula, which was also verified by Wang et al (2017). As the non-domestic buildings are mostly concentrated in Kowloon Peninsula and the northern part of Hong Kong Island, an area, marked by the red rectangle in figure 6, is chosen for further investigation. The averaged 2 m air temperature, water vapor mixing ratio, sensible heat fluxes, and latent heat fluxes are compared among the different scenarios ( figure 7). The use of water-cooled air systems can reduce the 2 m air temperature by around 0.5 • C-0.8 • C during the daytime as shown in figure 6(a), and the DCS case presents around 0.1 • C-0.15 • C lower air temperature as compared to the CPCT case. The temperature differences between the Baseline and the alternative systems reach 1.5 • C at 7:00-8:00 pm, as the impact of air-conditioning systems is more significant during the night-time (Fan and Sailor 2005). The use of cooling tower systems in commercial areas can significantly increase the water vapor mixing ratio in the atmosphere by around 0.29 g kg −1 on average. However, an increase in the water mixing ratio in the air will result in an increase in the apparent temperature and the heat index level, which will in turn increase energy consumption to remove the humidity in the supply air (Steadman 1984, Gutiérrez et al 2015. The cooling tower diminishes the sensible heat by around  80%-90% and transforms it into latent heat flux (Milosavljevic and Heikkilä 2001) during the afternoon with the maximum value delayed for around 1-2 as compared to the maximum sensible heat fluxes for the Baseline case, as shown in figure 7(b). As sensible heat fluxes are the sum of H tur and Base , the differences between the Baseline case and the CPCT and DCS cases with respect to the sensible heat fluxes indicate the heat fluxes from the air-conditioning systems ( figure 7(b)). The sensible heat fluxes in the CPCT case are slightly higher than those in the DCS case, while in the CPCT case, the latent heat fluxes contribute almost 90% of the energy emitted by the air-conditioning systems ( figure 7(b)). Ichinose et al (1999) reported that the anthropogenic heat emission in central Tokyo exceeds 400 W m −2 in the daytime, and the maximum value could be 1590 W m −2 . The anthropogenic heat flux in Hong Kong using satellite data was studied by Wong et al (2015), and the results show that commercial areas emit the largest anthropogenic heat fluxes around 500-600 W m −2 compared with other land-use types. The use of different airconditioning systems affects the surface energy balance, and air-cooled air systems (Baseline case) could significantly increase sensible heat fluxes through the air-conditioning systems and enhanced turbulence fluxes.

Effects on vertical structures
The effects of city-scale air-conditioning systems on vertical structures were investigated. The use of different air-conditioning systems could not only modify the surface energy balance but also extended to the areas above the ground level, which affects the development of the PBL height (Fan and Sailor 2005, Chen et al 2009, Sailor 2011, Gutiérrez et al 2015, as shown in figure S3. The PBL height presents diurnal variation, with the peak value of around 1400 m to 1600 m occurring at around 13:00-14:00 pm and around 500-600 m during the night (shown in figure S3). The implementation of the two alternative air-conditioning systems modifies the heat emitted and generation of turbulence, which induces a 100-150 m reduction in the PBL height. A reduction of the PBL height could cause potential air pollution issues (Han et al 2009, Gutiérrez et al 2015. The Baseline case has a relatively high urban canopy temperature, as shown in figures 8(a1)-(a4) on the vertical potential temperature, which leads to thermal discomfort for pedestrians and high building energy consumption (Nikolopoulou et al 2001. The vertical profile of the potential temperature (figure 8) shows a higher value below 300 m in the Baseline case than in the other cases at the same height, which illustrates that anthropogenic heat from the air-conditioning systems enhances vertical mixing (Salamanca et al 2014). Differences in the potential temperature between the different cases are reduced with an increase in height (figures 8(a1)-(a4)), which shows the effect of anthropogenic heat fluxes on the potential temperature distribution at the lower level (Chen et al 2009).   The momentum fluxes for the horizontal surfaces are determined by the bulk Richardson number based on Louis formulation, which is a function of roof roughness length, wind speed and the air temperature (Martilli et al 2002). The differences in the potential temperature affect the wind speed, and the Baseline case has the largest horizontal wind velocity (figures 8(b1)-(b4)). Figures 8(b1)-(b4) further illustrate that the wind direction is not affected by the different schemes of air-conditioning systems.
The selection of an air-conditioning system affects the sensible heat fluxes released into the atmosphere, as discussed in section 2.2, which has a direct effect on the buoyant production of turbulent kinetic energy (TKE) (Martilli et al 2002). The results show that the buoyant production of TKE has the largest magnitude for the Baseline case and the smallest for the no-heat case (figures 8(c1)-(c2)). The largest differences in TKE were found 100 m above the ground level, above the mean building height of the region. Furthermore, figure 9 shows that the differences in the buoyant TKE are extended to a height of around 1.3 km above the ground level in the commercial areas, which has no influence in the non-domestic areas. The decreases in TKE within the urban boundary layer lead to potential air pollution issues, particularly pollutants near the ground level (Chen et al 2009, Salamanca et

Effects on SLBC and UHIC
The effect of different air-conditioning systems on local wind circulation, particularly UHIC and SLBC, is investigated in the urban areas (vertical East-West cross-section, transect A-B, shown in figure 9(a)) in Kowloon Peninsula. At 06:00 UTC on 24 June 2016, a convergence zone was observed in the urban areas of Kowloon Peninsula (Tong et al 2005), with sea breeze penetration on the left side and a downslope wind on the right side. The urban areas in the convergence zone area are characterized by a low water vapor mixing ratio and enhanced UHIC of up to around 1.5 km (as shown in figures 10(a) and (d)). On the right side of the convergence zone, the UHIC is stronger than on the left side with the assistance of the airflow from across the mountain (as shown in figure 10(a)). Across the mountain, the water vapor mixing ratio is also reduced ( figure  10(a)). To further investigate the effects of different air-conditioning scenarios on the local circulation, the difference in vertical wind fields is examined. The CPCT case has a higher water vapor mixing ratio near the ground level in the urban areas (on the left side of the convergence zone, shown in figure 10(e)), which is in agreement with Gutiérrez et al (2015). Figures 10(b) and (c) show that the Baseline case has the largest sensible heat fluxes and can enhance the downslope wind on the mountain side merged into the UHI circulation. The DCS case has less sensible and latent heat fluxes, which leads to a large vertical velocity difference, as shown in figure 10(c). Both cases ( At 06:00 UTC on 25 June 2016, the background wind was southeasterly, and the urban areas in the cross-section were dominated by a strong eastern wind and no sea breeze penetration ( figure S4(a)). The urban areas were characterized by a low water vapor mixing ratio as well. The largest vertical wind velocity was near the western edge of the cross-section, and the UHIC moved west in comparison to that shown in figure 10(a). The Baseline case, with the largest sensible heat fluxes within the urban canopy layer, had a strong UHIC. The changes in air-conditioning systems could not only affect the vertical distribution of the air temperature and the water vapor mixing ratio, but also the secondary local circulation, particularly the sea breeze circulation and the UHI circulation.

Discussion
The Baseline case has a higher air temperature within the urban canopy layer than the CPCT and DCS cases, which implies that the alternative systems can enhance pedestrian thermal comfort from this aspect. However, the CPCT and DCS cases show a decrease in sensible heat fluxes (figure 7), horizontal and vertical wind speed (figure 8), buoyant production of TKE Figure 11. Comparison of air-conditioning energy consumption in the non-domestic areas in the specific areas shown in figure 4 for the Baseline case, CPCT case, DCS case, and no-heat case (upper), and the energy saving ratio for the CPCT and DCS cases as compared to that in the Baseline case (below).
(figures 8 and 9), and PBL height (figure S3). The airconditioning systems also affect secondary circulation, SLBC, and UHIC. The changes may cause air stagnation and air pollution issues near the ground level, particularly in the convergence areas in Kowloon Peninsula (Tong et al 2005, Chen et al 2009, Xie et al 2016. Xie et al (2016) clearly showed that the changes in the meteorological conditions due to the anthropogenic heat fluxes can significantly affect the spatial and vertical distributions of air pollutants, particularly in the big cities in China. The role of sea land breeze circulation on the coastal and inland areas has been investigated (Kitada 1987, Melas et al 1995, Clappier et al 2000. An increase in the pollutant level usually combines with a relatively weak sea breeze and unfavorable conditions for pollutant dispersion (Papanastasiou and Melas 2009). Therefore, the hypothesis for air quality issues needs to be confirmed in the future.
The energy consumption using different airconditioning systems is investigated, and in this study, the COP for the Wall_type, Baseline, CPCT, and DCS cases are set as 2.0, 2.5, 3.5, and 3.5, respectively, on the basis of the findings of other studies (Ashie et al 1999, Ho et al 2007, Huang et al 2013. Figure 11 shows the averaged air-conditioning fluxes in the non-domestic areas (as shown in figure 6) and the energy saving with the alternative air-conditioning systems referring to the Baseline case. The results clearly show that in the CPCT and DCS cases, the energy consumption of the air-conditioning systems was effectively reduced, with 30% energy savings in comparison to that in the Baseline case. Note that here, the cost of changing the systems and maintenance is not considered. In making the decision to choose alternative systems, one must consider the competing benefits and disadvantages.

Conclusions
This study is a step forward in investigating the effects of two alternative air-conditioning systems (DCS and CPCT) on urban climate according to the real-life conditions of Hong Kong. The two air-conditioning systems are physically based, parameterized, implemented in the WRF model coupled with the BEP-BEM urban canopy model, and assessed during an extreme high temperature event at the city scale of Hong Kong. Note that the performance of the alternative airconditioning systems under different climatic scenarios is beyond the scope of this study but will be studied in the future.
The results show that both water-cooled air systems could reduce the 2 m air temperature by around 0.5 • C-0.8 • C during the daytime, and the temperature differences between the Baseline and the two alternative systems could reach 1.5 • C at 7:00 pm-8:00 pm when the PBL height was confined to 100 m. The sensible heat fluxes in the CPCT case were slightly higher than in the DCS case, while in the CPCT case, the latent heat fluxes contributed almost 90% of the energy emission from the air-conditioning systems. The cooling tower (CPCT case) diminished around 80%-90% of the sensible heat fluxes and transformed them into latent heat fluxes, which could significantly increase the water vapor mixing ratio in the atmosphere by around 0.29 g kg −1 on average.
Environmental and meteorological effects are not confined to the rooftop where the air-conditioning systems are located, but can extend to the surface and alter the vertical distribution of various parameters within the urban boundary layer. The implementation of the two alternative air-conditioning systems modified the heat and turbulence, which affected the evolution of the PBL height (with a reduction in height of 100-150 m), and also reduced vertical mixing. Moreover, the CPCT and DCS cases exhibited lower horizontal wind speed and buoyant production of TKE, which might cause air pollution issues, particularly near the ground level. Hence, the hypothesis on the air quality issue needs to be confirmed in the future. The alternative airconditioning systems could affect the secondary local circulation, which reduces the strength of the SLBC and UHIC in the urban areas.
Certain aspects of energy consumption were illustrated, and the CPCT and DCS cases could effectively reduce the air-conditioning energy consumption by 30% as compared to the Baseline case. The non-extreme high temperature period (June 8-9, 2016) showed that water-cooled air systems could reduce the 2 m air temperature by around 0.6 • C-1.0 • C in the commercial areas, which was higher than in the June 23-28, 2016, period, and had similar effects on the boundary dynamics. The new results from this study are expected to help define UHI mitigation strategies. However, installation and maintenance costs are not considered. In making the decision to choose alternative systems, one must consider the competing benefits and disadvantages.