Seasonal Variation of Urban Heat Island and Its Impact on Air-Quality Using SAFAR Observations at Delhi , India

This paper discussed the urban heat island (UHI) intensity and local air quality by using observational data of project of the System of Air Quality Forecasting and Research (SAFAR) over Delhi during the month of May and December 2013. It is found that UHI magnitudes ~2.2 ̊C and ~1.5 ̊C are formed at the evening traffic hours during May and December respectively. Also, intensity of UHI < 0 ̊C over daytime is referred as Urban Cool Island (UCI) during May and December. The diurnal PM2.5 concentration shows a bimodal pattern with peaks at morning and evening traffic hours during May and December. The planetary boundary layer height (PBLH) values show higher in magnitude during the daytime and lower in magnitude during the night-time. Whereas, the Ground Heat Flux values are lower during the daytime and higher during the night-time. The wind speed shows lower values during the UHI and higher magnitudes during the UCI formation hours. Concentration of PM2.5 and wind speed shows a strong negative correlation during May (r = −0.56, p = 0.002) and December (r = −0.57, p = 0.001) at C V Raman (CVR) site, however, high values in the concentration of PM2.5 during the low wind speed favour the condition for the formation of UCI. The regression analysis indicated that PM2.5 plays a significant role in the daytime cooling and nighttime warming over the urban areas during the low wind speed condition.


Introduction
The urban heat island (UHI) effect is simply defined as the temperature difference between an urban and the surrounding rural regions.This impacts to re-lease in additional heat from the urban environment which is a major contribution to the urban heat island [1].The UHI effect appears more frequently under calm and weak wind conditions with a cloudless sky at night and that its intensity is affected by daytime and nocturnal circulations [2] [3].Many studies have mentioned that the association between the UHI phenomenon and air pollution with the favourable conditions such as high temperature, low wind speed, low relative-humidity and cloudless, which are important for UHI development, easily triggers air pollution episodes [4] [5].
The basic heat transfer and energy conservation processes, such as conduction, convection and radiation play their characteristic roles in this heat exchange.The structures on ground level, such as walls and roof facets, irrigated gardens, non-irrigated green spaces, lawns and paved areas among others, capture solar radiation to different extents.These natural and man-made structures continuously absorb and store this radiation in the form of heat energy from sunrise till late afternoon.Afterwards, the sun starts setting a side and the environment starts cooling down.The heat energy stored in structures is then released to the environment.The method and quantity of heat released by the urban structures, however, depend on other controllable factors such as the sky view factor and building material.In a typical urban area, massive construction material is placed within a very small space that captures high intensity of solar radiation.The ability of heat release by long-wave radiation in cities is low due to decreased sky view which results in high heat storage in building structures.It is also believed that air pollutants, in particular aerosols that are abundant over polluted urban areas, can absorb and re-radiate long wave radiation and inhibit the corresponding radiative surface cooling producing a pseudo-greenhouse effect, which is responsible for causing UHI [6].PM 2.5 , which refers to particulate matter (PM) in air that is less than 2.5 μm in aerodynamic diameter, is a key pollutant affecting radiation balance and a main factor in decreasing air quality [7].Urban landscape features have obvious effects on particulate matter, related closely with PM 2.5 concentration [8].PM 2.5 mainly emitted from our daily activities, such as vehicle exhaust, marine aerosols, coal and fuel combustion, burning of agricultural wastes, paved road dust, and secondary sulfates, etc., [9].
Dispersion and transport of lower atmospheric pollutant depend largely on the local planetary boundary layer height (PBLH) structure.Thus, PBLH has been used as a key in weather, climate, and air quality models to determine turbulence mixing, vertical diffusion, convective transport, and atmospheric pollutant deposition [10] [11].As this is one of the critical physical parameters in atmospheric environmental evaluation, studying PBLH's association with ground level PM 2.5 concentration is very important for air pollution and environment quality monitoring.
The UHI has been the focus in many studies [12] [13] [14] [15] while a significant number of these studies document that cities often remain cooler than the countryside from the early morning until the early afternoon during fair weather and low wind speed conditions.Some earlier studies qualitatively suggest that the UHI may originate from shadow effects in the urban canyon [16], the daytime energy storage in the urban, attenuate the net radiation due to aerosols [17], or the difference in land cover and the available soil moisture altering the surface energy balance.

Observational Datasets
In this study we mainly focus on the hourly Meteorological data of air tempera- in Figure 1(a).The UHI intensity is calculated by the simple temperature gradient between the urban areas (T u ) to its surrounding rural areas (T r ) (i.e., T UHI = T u − T r ).

Model Datasets
The The main physics schemes tested and used in the model are microphysics scheme from [18] shortwave radiation scheme from [19] longwave radiation scheme [20] from Rapid Radiative transfer model, Monin-Obukhov surface layer scheme [21], Yonsei University boundary layer scheme and the Noah land surface model.While no cumulus convection scheme is used for the inner domains i.e. (2 & 3) of the model, the Kain Fritch cumulus convection scheme is applied for the outer domain.
A simulation of the particular intricacies of such an environment is heavily influenced by its depiction of the planetary boundary layer (PBL)-that portion of the lower troposphere directly affected by the earth's surface via tropospheresurface exchanges of heat, moisture, and momentum on sub hourly time scales [22] [23].Exchanges of moisture, heat, and momentum occur within the PBL through mixing associated with turbulent eddies.These eddies influence the way in which lower-tropospheric thermodynamic and kinematic structures evolve.
Such eddies operate on spatiotemporal scales that cannot be explicitly represented on grid scales and time steps employed in most Mesoscale models.As such, their effects are expressed in these models via the use of PBL parameterization schemes, whose theoretical development is outlined in multiple sources addressing the subject [22] [23] [24], Here we use YSU PBL scheme the main advantage of this scheme is More accurately simulates deeper vertical mixing in buoyancy-driven PBLs with shallower mixing in strong-wind regimes [25].during the month of December 2013, that UHI intensity shows a peak late in the evening around 20:00 hrs with a magnitude of ~1.5˚C, thereafter it decreases till the sunrise in the morning.During the daytime the magnitude in the intensity of UHI continuously increases with the morning UCI to the night time UHI effect.The PM 2.5 concentrations over the CVR site shows a bimodal pattern in both the month May and December with the morning peak around 10:00 hrs and the second peak in the evening peak around 20:00 hrs, the peaks in the PM 2.5 concentration hours are related to the morning and evening traffic hours [26].The peaks in the concentration of PM 2.5 during the traffic hours are due to the enhanced anthropogenic activity during the traffic hours [27].Whereas the concentration of PM 2.5 over IITMD site during May 2013 starts increasing after the sunrise hours and peaks around 10:00 hrs and remains almost constant till the evening around 17:00 hrs and then it starts decreasing and remains lower level during the night time as in Figure 2(a) and diurnal variation of PM 2.5 concentration at the IITMD location follows the same trend with a lesser in the concentrations compared to the CVR site during December 2013 as in Figure 3(a).The PBLH is mainly affected by temperature, ground heat flux and wind speed among others [28].Here the PBLH remain almost at constant height after sunset hours to the sunrise hours and starts increasing thereafter and peaks to a maximum during the morning traffic hours and remains constant till the solar insolation is maximum and there after starts decreasing till sunset hours.The PBLH follows the same pattern at the IITMD site with a little lower value at the night time and higher during the day during both months May and December 2013.The UHI formation during the evening traffic hours is due to increase in the concentration of PM 2.5 with decrease in the height of the PBLH and higher in the magnitude of the ground heat flux and lower wind speed.Before sunrise the concentration in the PM 2.5 are lower in the magnitude and with the decrease in the PBLH and lower in the ground heat flux values and lower in the wind speed, which leads to the formation of the Urban Cool Island (UCI).The formation of UHI is mainly favoured with the low wind speed and high concentration of PM 2.5 with a positive correlation between the UHI and PM 2.5 concentrations.Whereas, a negative correlation between the wind speed and concentrations of PM 2.5 [31] during both May and December.The results from the correlation analysis in the present study suggest "Urban Cool Island" formation over Delhi during daytime in the months of May and December under low wind conditions.Thus it is found that during the daytime the extinction of incoming solar radiation by the aerosols are more dominating than the absorption of outgoing long wave radiation.Thus results suggest the daytime negative forcing by the aerosols during May and December under low wind speed and clear sky conditions [32] [33].

Conclusions
We performed analysis of site-specific air quality and weather data from the project SAFAR over Delhi during May and December 2013.UHI of ~2.2˚C and ~0.5˚C observed over nighttime with a peak around 20:00 hrs during May and December respectively.The PM 2.5 concentration shows a bimodal pattern with first peak at 10:00 hrs and the second peak at 20:00 hrs.The peaks in the diurnal concentration of PM Also, concentration of PM 2.5 and wind speed shows a strong negative correlation at CVR (r = −0.56,p = 0.002; r = −0.57,p = 0.001) and IITMD (r = −0.51,p = 0.002; r = −0.56,p = 0.001).The results from the regression analysis of wind speed and PM 2.5 indicated the significant role played by the PM 2.5 levels in the daytime cooling and nighttime warming over the urban areas which moderated the diurnal UHI during the low wind speed condition.

Figure 1 .
Figure 1.Land use land cover pattern of Delhi and its surroundings, (a) where star mark indicates the station from the SAFAR network and solid line indicated the Delhi boundary, and images from the Google earth for the location of (b) CVR and (c) IITMD and its surroundings.
Mesoscale model used in this study is the Advanced Research WRF (ARW) model version 3.7.Our model simulation covers the PBLH and ground heat flux for the months of May and December 2013.The simulations were run at a three interactive domains with spatial resolution of 9 km, 3 km and 1 km which covers part of North India (22˚N -34˚N -71˚E -84˚E), second Domain covers complete Delhi and outer boundaries (25˚N -32˚N -73˚E -80˚E) region and Third domain covers Delhi (27˚N -30˚N -75˚E -78˚E).The vertical grid contained 32 full sigma levels from the surface up to 50-hPa.Approximately eight of these levels were below 1 km, thereby providing a fine vertical resolution within the planetary boundary layer.The terrain and land use data for the entire 3 domain were taken from 30 s data available from the United States Geological Survey (USGS).The National Centre for Environmental Prediction (NCEP) Final Analysis data (FNL) available at 6 hourly intervals on a 1˚ × 1˚ resolution are used as initial and boundary conditions.The model uses two-way nested boundary conditions and they are updated every 6 hours.

Figure 3 .
Figure 3.Diurnal cycle of UHI in comparison with the (a) PM 2.5 concentration, (b) PBLH, (c) Ground heat flux and (d) Wind speed at CVR and IITMD sites during December 2013.

Figure 2 (
b) and Figure 3(b) show diurnal cycle of UHI and PBLH during the month of May and December 2013.

Figure 2 (
c) and Figure 3(c) show the diurnal cycle in the UHI and the Ground heat flux during both months May and December 2013.Here the Ground heat flux value remains higher from sunset to the sunrise hours and starts decreasing there after the sunrise hour and attains a least value during 09:00 hrs and there after it starts increasing and attains a maximum during 17:00 hrs in May month and in December month it attains the least values around 06:00 hrs and there after it starts increasing and attains a maximum during 11:00 hrs.Whereas at the IITMD site the ground heat flux values follows the same trend as CVR site values with lesser in magnitude in the morning hours and a little slight higher in magnitude during the night hours during both the months May and December 2013.The strong decrease in the ground heat flux values between the nighttime and daytime suggest that more thermal energy absorption by the urban surfaces [29].The diurnal variation in the wind speed during May and December 2013 over the CVR and IITMD sites are seen in Figure 2(d) and Figure 3(d).Here the wind speed shows peaks during 06:00 hrs and 08:00hrs after the sunrise and then start decreasing and attains a least values and remain constant from 13:00 hrs to 22:00 hrs in the mid-night during May 2013 as in Figure 2(d).But the wind speed shows a peak during 08:00 hrs in the morning and 15:00 hrs in the evening, whereas during the night hours the winds are low as compares to that of the morning hours.The winds are lower in the magnitude during the same evening traffic hours when the UHI intensity is higher around 20:00 hrs during December 2013 as in Figure 3(d).The wind speed at the IITMD site follows the same trend with a lesser in the magnitude during both months May and December 2013.Changes in the concentration of PM 2.5 , PBLH and Ground heat flux are mainly caused due to the changes in the magnitude of the wind speed [30].

Figures 4 (Figure 4 .
Figures 4(a)-(d) show regression analysis to explore the correlation between the UHI and its two possible dependent variables during May 2013.The first variable is the concentration of PM 2.5 , here Figure 4(a) shows the correlation of UHI with the first variable exists a strong positive correlation at both CVR (r =

Figures 5 (Figure 5 .
Figures 5(a)-(d) show the regression analysis to explore the correlation between the UHI and its two possible dependent variables during December 2015.The first variable is the concentration of PM 2.5 , here Figure 5(a) shows that the correlation of UHI with the first variable is positive at both CVR (r = 0.62, p = 0.001) and IITMD (r = 0.40, p = 0.003) sites.The second variable, wind speed from Figure 5(b) shows a negative correlation of UHI at the CVR (r= −0.40, p = 2.5 correspond to morning and evening traffic hours during May and December.The evening peak in the concentration of PM 2.5 , ground heat flux values, low PBLH and wind speed leads to the formation of UHI.On the other hand, during daytime, peak in the concentration of PM 2.5 during the morning traffic hours, high PBLH, low intensity of Ground heat flux and less wind speed leads to a peak in the intensity of UHI < 0˚C leading to the formation of UCI during May and December.During May and December, the concentration of PM 2.5 shows a positive correlation with UHI intensity at CVR site (r = 0.61, p = 0.001; r = 0.57, p = 0.001), and IITMD site (r = 0.65, p = 0.001l; r = 0.40, p = 0.003) and negative correlation of wind speed correlation with UHI intensity at the CVR (r = −0.40,p = 0.001; r = −0.40,p = 0.002) and IITMD (r = −0.41,p = 0.001; r = −0.27,p = 0.005).