Generalized Scaling of Urban Heat Island Effect and Its Applications for Energy Consumption and Renewable Energy

In previous work from this laboratory, it has been found that the urban heat island intensity (UHI) can be scaled with the urban length scale and the wind speed, through the time-dependent energy balance.The heating of the urban surfaces during the daytime sets the initial temperature, and this overheating is dissipated during the night-time through mean convection motion over the urban surface. This may appear to be in contrast to the classical work by Oke (1973). However, in this work, we show that if the population density is used in converting the population data into urbanized area, then a good agreement with the current theory is found. An additional parameter is the “urban flow parameter,” which depends on the urban building characteristics and affects the horizontal convection of heat due to wind. This scaling can be used to estimate the UHI intensity in any cities and therefore predict the required energy consumption during summer months. In addition, all urbanized surfaces are expected to exhibit this scaling, so that increase in the surface temperature in large energy-consumption or energy-producing facilities (e.g., solar electric or thermal power plants) can be estimated.


Introduction
There have been many studies on the causes and impact of the heat island effect (e.g., [1][2][3][4][5][6][7][8][9]).A wide range of causes have been attributed, including the land surface energy balance (or imbalance) due to urbanized land surface and built structures, anthropogenic heat release, and different atmospheric constituents over the city [10,11].While the significance of the urban heat island effect is universally recognized [12,13], the approaches to its understanding and mitigation are quite varied, ranging from simple "slab" approach [14] to sophisticated models that resolve urban canopy physics [15,16].One of the important aspects is quantifying the magnitude of the urban heat island effect in terms of relevant variables, which again vary widely depending on the authors.A key step in this quantification is "scaling, " that is, developing a functional relationship between the urban heat island effect and causal variables.
In previous work from this laboratory, we have shown that scaling of the urban heat island effect based on timedependent energy balance works quite well in understanding and predicting the UHI dependence on the urban length scale and wind speeds [17].This approach is in contrast to the scaling used by many authors [18][19][20][21][22][23][24], where a steady-state energy balance is used along with horizontal and vertical heat flux terms.It is clear that the urban island heating (as well as cooling) is not a steady-state effect, nor is the UHI intensity determined at some fixed time, as shown in Figure 1 where a typical time evolution of the UHI intensity is plotted.In this work, we show that if the population density is used in converting the population data into urbanized area in Oke's work [9], then a good agreement with the current theory is found.An additional parameter is the "urban flow parameter, " which depends on the urban building characteristics and affects the horizontal convection of heat due to wind.This scaling can be used to estimate the UHI intensity in any cities and therefore predict the required energy consumption during summer months.In addition, all urbanized surfaces are expected to exhibit this scaling, so that increase in the surface temperature in large energy-consumption or energyproducing facilities (e.g., solar electric or thermal power plants) can be estimated.About one-third of the electric 2 Advances in Meteorology power goes to satisfy the residential consumption in United States, and thus understanding of the electricity demand due to UHI increase is a key component in planning for renewable energy production.

Methods of Analyses
A brief summary of the scaling is given below.The starting point is the fundamental energy balance equation [17]: where   is the divergence of the eddy-diffusive heat flux.The air density and specific heat are  and   , respectively, while  is the convection speed, often approximated as the mean wind velocity over the urban area.The spatial gradient in the temperature can be approximated and moved to the right-hand side.In our previous work [17], we considered the time evolution of the urban temperature during night hours, after the heat flux from the surface has led to some "equilibrium" temperature above the urban surface.This argument considers the heat flux term to set some initial urban temperature, prior to subsequent cooling, and thus heat flux term in (1) has been removed and absorbed into the initial temperature difference (  −   )  .For now, we can retain this term, as it can have additive effect on the urban air temperature during the night hours.The temperatures,   and   , represent temperatures in the urban and rural locations, respectively, and (  −  ) has been used as the UHI intensity: If we integrate while retaining only the first term on the right-hand side, we obtain the temporal evolution of the UHI intensity.This assumes that the convection effect is much stronger than the heat flux term: where  = 1 in the limiting case of   = 0 and  < 1 in general, if the diffusive heat flux retarding the night cooling process is taken into account.This exponential decay in the temperature difference is closely reflected by typical data for time evolution at night hours in Figure 1.
On the other hand, if we consider the other extreme case of zero wind speed and retain only the heat flux term, then we obtain From (3), we can see that the wind speed, , has the obvious effect of enhancing the cooling of the heated air over the urban surface.The larger the urban length scale is, , the longer it would take to cool   .The heat flux term, as shown in ( 4), adds to the urban temperature until the heat is completely dissipated.The urban-rural breeze usually results in an enhanced convection motion due to mean wind in urban areas, which dominates the heat transfer process over the heat flux term due to eddy diffusion.For compactness, we retain the solution form of (4) and use a correction factor, , to reflect the facts that (1) the temperature decay will be delayed by the heat flux term and (2) different urban surface morphologies, aerodynamic properties, and atmospheric stability may affect the effectiveness of wind speed in cooling of the urban temperature.Thus, a large heat flux term and structural impedance to convection cooling will both result in smaller value of .We can observe from (3) then that the cooling should undergo an exponential decay as a function of time (Figure 1), and depending on the above factors the final minimum temperature achieved will be the so-called night-time minimum temperature typically reported in meteorological data sets.In addition, the wind speed also should lead to an exponential decrease in UHI.
A large set of data has been used to validate this model, and an example is shown in Figure 2 where the wind speed effect for various urban areas is shown to follow the exponential decrease according to (3).In this work, we have used the data reported by Oke [9], along with the population and urban surface area data dated at the same time period (1968)(1969)(1970)(1971)(1972)(1973)(1974)(1975)(1976)(1977)(1978), in order to offer further proof that the correct scaling for the UHI intensity should be based on time-dependent energy balance and that there is exponential dependence as a function of time, wind speed, and the urban length scale.The following Canadian cities are considered as in Oke [9]: Montreal, St. Hyacinthe, St. Hubert, Chambly, Marieville, St. Basile-Le Grand, St. Cesaire, St. Pie, Ste.Angele-de Monnoir, and Ste.Madeline.Symbols represent data as reported by Oke [9].The upper curve is (3) with  = 1.0, while the lower is for  = 0.25.

Results
Figure 3 shows the UHI intensity as a function of the urban length scale, for all of the cities reported in Oke [9].In Oke's work, the UHI intensity was correlated and plotted as a function of the population ("city size"), but no sound physical basis was given.Other authors have attempted a similar correlation, but different slopes were required for North American, European, Japanese, and Korean cities. Figure 3 shows that the UHI intensity for Canadian cities can be plotted using the current scaling, within the upper ( = 1) and lower ( = 0.25) bounds, using the urban area converted to the urban length scale of the identical time periods.As noted above and in our previous work, the parameter  depends on the urban morphological characteristics.Thus, there is a direct effect of increase in the urban length scale on the UHI intensity.In our previous work [17], we have investigated the growth of two cities, Phoenix and Tucson, based on the satellite image analyses, and corresponding increase in the urban length scale.Of course, the growth of urban areas depends on many economical and sociological parameters, as the expected growth of Phoenix, for example, has been halted in 2008 due to housing meltdown.
Thus far, in this and our previous work [17], we have looked at cities around the world with widely disparate population densities, urban surface morphology, and meteorological conditions.Although the general UHI characteristics are captured by (3), exact functional relationships tend to vary as shown in Figure 3 even for cities within Canada.Thus, incorporating and categorizing the effects of urban surface morphology and meteorological conditions into the parameter  in (3) remain as the next task in this formation for the UHI scaling.
Figures 4, 5, and 6 show the effects of the wind speed on the UHI intensity for the same Canadian cities as in Figure 3.It can be observed that functional relationship between the UHI intensity and the wind speed is accurately captured using (3).Hunt et al. [25] have also looked at the scaling of UHI effects for large urban areas and attempted to link the worsening of UHI effects on global warming.As shown in (3), the initial temperature difference between urban and rural areas does determine the final UHI intensity.However, the more direct factors are the urban length scale and the wind speed, since the global warming would elevate the temperatures in urban and rural areas uniformly.

Conclusions
Effects of the urban length scale and the wind speed on the urban heat island effect can be quantified through timedependent energy balance.The heating of the urban surfaces  during the daytime sets the initial temperature, and this overheating is cooled during the nighttime through mean convection motion over the urban surface, resulting in an exponential decay in the temperature.This overheating is cooled during the nighttime, through mean convection motion over the urban surface, resulting in an exponential decay in the temperature.The solution to the time-dependent energy balance equation reproduces this temporal decay with good accuracy, with the main factors being the length scale of the urban area and the wind speed.The minimum temperature reached at the end of night-time cooling period then corresponds to the UHI intensity, which increases with increasing urban length scale and decreasing wind speed.The wind speed effect is also accurately retraced using this method; however, different correction factors are required for different cities, indicating the effects of the urban surface heat content, structural morphology, and density.Thus, using a small number of readily available data for the urban length scale and the wind speed, the UHI intensity can be described with possible projections for future trends.This approach can be used in planning of energy resources, as well as any large areas of concrete surfaces needed for renewable or other power generation.About one-third of the electric power goes to satisfy the residential consumption in United States, and thus understanding of the electricity demand due to UHI increase is a key component in planning for renewable energy production.

Figure 1 :
Figure 1: Temperature decay as a function of time, in Phoenix, AZ, for years 1960 to 2010, using (3).

Figure 2 :
Figure 2: The effect of the wind speed on the UHI intensity in Seoul, Korea, St. Louis, USA, and Melbourne, Australia.The theoretical lines are obtained using (3).(  −   ) is the UHI intensity defined as the maximum difference in the temperatures in the urban and rural locations.

Figure 3 :
Figure3: UHI intensity as a function of the urban length scale.Symbols represent data as reported by Oke[9].The upper curve is (3) with  = 1.0, while the lower is for  = 0.25.

Figure 4 :
Figure 4: The effect of the wind speed on the UHI intensity in Chambly, Marieville, St. Basile-Le Grand, St. Cesaire, St. Pie, Ste.Angele-de Monnoir, and Ste.Madeline.Symbols represent data as reported by Oke [9], while the theoretical lines are obtained from (3).

Figure 5 :
Figure 5: The effect of the wind speed on the UHI intensity in Chambly, Marieville, and St. Basile-Le Grand.

Figure 6 :
Figure 6: The effects of wind speed on UHI intensity for St. Cesaire, St. Pie, and Ste.Angele-de Monnoir.