Development of an operational modelling system for urban heat islands : an application to Athens , Greece

Introduction Conclusions References


Introduction
Urban expansion results to significant land cover changes.As cities sprawl into the natural habitat, roads, buildings and other structures, necessary to accommodate growing populations, replace vegetation and natural surfaces.Such land cover changes significantly influence the local climate and weather.This has been extensively documented in the scientific literature, so that climatologists consider cities to have unique local climates, pretty much like forests, lakes, valleys and coastlines.
The urban heat island (UHI) effect, describing the excess warmth of the urban atmosphere compared to the rural surroundings (Voogt and Oke, 2003), is one prominent nomenon that is significant in itself since it occurs in that layer of the atmosphere where almost all daily human activities take place.Research has shown that heat islands have considerable implications for air quality (Rosenfeld et al., 1998;Sarrat et al., 2006;Davies et al., 2007), energy demand for cooling purposes (Santamouris et al., 2001;Konopacki and Akbari, 2002), human health (Conti et al., 2005;Haines et al., 2006), and regional/local atmospheric circulations (Lin et al., 2008;Miao et al., 2009).Further, although the UHI in itself does not influence global temperatures (Houghton et al., 2001), it does have an effect on local temperatures used for assessing climate change (Van Wevenberg et al., 2008).
According to the pioneering study of Oke (1982), heat islands develop as a result of the differences in the surface energy balance (SEB) between urban and rural areas.These differences are due to the combination of the thermal, radiative, aerodynamic and moisture properties of the urban fabric that significantly differ from those of the natural landscape (Oke, 1987).Briefly, urban areas store more heat during the day than do rural areas, due to the generally lower surface albedo (Taha, 1997;Christen and Vogt, 2004) and the surface enlargement provided by the canyon-like geometry (Christen and Vogt, 2004).Following sunset, the rural surroundings begin cooling quickly due to the generally unobstructed sky view.Meanwhile, the cooling rate of the urban environment is significantly reduced because of the decreased sky view factor (SVF) and the increased release of storage heat from the urban surfaces (Grimmond and Oke, 1995).As a result, urban and rural temperatures begin to diverge, generating the UHI, until after sunrise the daily solar radiation cycle begins and the urban/rural areas start warming quickly.Therefore, the UHI is primarily a nocturnal phenomenon resulting from the diverging urban/rural cooling rates (Oke, 1982(Oke, , 1987)).Numerical modelling is considered to be an appropriate tool for studying the UHI effect, showing a wide area of successful applications (Rizwan et al., 2009).In particular, meso-scale atmospheric models have been extensively used during the past decade for investigating heat islands.This has been made possible through the introduction of a new generation of urbanization schemes to meteorological models such as the Fifth-Introduction

Conclusions References
Tables Figures
Properly designed modelling systems can be used as a powerful tool for detailed and continuous monitoring of the urban thermal environment.The implementation of such systems for investigating the UHI could provide valuable information to civil protection and urban planning authorities, but also to the general public.Up to nowadays, this has been well documented in the scientific literature and numerous heat island modelling studies have been conducted in cities as diverse as Athens, Greece (Giannaros et al., 2013), Goteborg, Sweden (Miao et al., 2007), Oklahoma City, USA.(Liu et al., 2006), and Taipei, Taiwan (Lin et al., 2008).To the authors' knowledge, however, modelling of the UHI has yet to be realised in the operational context.The current paper presents the development of an urban modelling system, designed to support operational real-time weather forecasting activities.The modelling system exploits the capabilities of a state-of-the-art numerical weather prediction (NWP) model together with techniques for improving the representation of urban areas.The philosophy upon which it is built is the provision of reliable, high-resolution forecasting products, related to the heat island effect and human thermal comfort.In this context, several new developments and modifications were carried out, aiming to improve the quality of the final products without increasing computational demands.The presented modelling system is practically designed to support applications over any urban area of the European continent.This study focuses on the operational implementation of Introduction

Conclusions References
Tables Figures

Back Close
Full the system for the city of Athens, Greece, during a 1.5-month period in the summer of 2010.

Description of the study area
Athens (37 • 58 N, 23 • 43 E) is the capital and largest city of Greece.It is located in a small peninsula situated in the south-eastern end of the Greek mainland (Fig. 1a).
The greater Athens area (GAA) covers approximately 450 km 2 and the urban zone sprawls across a basin that is often referred to as the Attica basin.The city is surrounded by fairly high mountains to the west (Mount Aigaleo, 469 m), north (Mount Parnitha, 1413 m), north-east (Mount Penteli, 1109 m) and east (Mount Hymettus, 1026 m).
Athens is also influenced by the sea to the south (Saronic Gulf).The industrial zone of the city is located in the western part of the Attica basin (Thriassion Plain), while the Mesogeia Plain dominates in the south-east of the peninsula (Fig. 1b).
Athens enjoys a typical Mediterranean climate with hot, dry summers and cold, mild winters.The average daily maximum temperature in the summer months of July and August approximates 32 • C, often exceeding 37 • C when a heat wave occurs.The city also experiences a strong UHI effect as reported in several previous studies (e.g.Katsoulis and Theoharatos, 1995;Philandras et al., 1999;Santamouris et al., 2001;Keramitsoglou et al., 2011).
3 Overview of the modelling system The initialization can be carried out using data from a larger scale atmospheric model, such as the Global Forecasting System (GFS) operated by the National Centre for Environmental Predictions (NCEP).The primary output of the UHI-MS consists of three-day high-resolution forecasts for near-surface air temperature.
A detailed description of the various components, input/output data, and procedures of the UHI-MS is given in the following sections.

The WRF/Noah modelling system
As seen in Fig. 2, the meteorological driver of the UHI-MS is the WRF model, version 3.2 (Skamarock et al., 2008), coupled with the Noah LSM (Chen and Dudhia, 2001).The key advantage of this architecture lies in that it allows for a computationally efficient representation of the urban thermal environment.In particular, the coupled WRF/Noah modelling system includes a bulk urban parameterization that can be used for representing zero-order effects of urban surfaces (Liu et al., 2006).Within this simple urban scheme, urban areas are treated as homogenous surfaces whose radiative (albedo, emissivity), thermal (heat capacity, thermal conductivity) and other (roughness length, green vegetation fraction) properties are properly modified to account for the urban effects.This approach is considered to be suitable for operational forecasting applications, as highlighted in the study of Liu et al. (2006).

Representation of urban land use
The employment of the urbanization scheme of Liu et al. (2006) requires that land use be represented as accurately as possible.This is because the land surface parameters are specified in the WRF/Noah as a function of the land use category.
In the standard WRF configuration, land use can be represented using either the 24-category and 30-arc-sec spatial resolution US Geological Survey (USGS) global land cover characteristics dataset or the 20-category and 30-arc-sec spatial resolution International Geosphere-Biosphere (IGBP) global land cover dataset.Despite their rel-Introduction

Conclusions References
Tables Figures

Back Close
Full atively high horizontal resolution, which approximates 1 km, both datasets are considered to be inadequate for urban modeling purposes.This is mainly because they include a single category for representing urban areas (Skamarock et al., 2008).Further, the USGS dataset is considered to be outdated since it is derived from 1 km AVHRR (Advanced Very High Resolution Radiometer) data spanning from April 1992 through May 1993.Nevertheless, both datasets do provide an option for detailed urban land use representation, as long as such data are available.
To enable the detailed urban land use representation in the WRF/Noah modeling system, the 250 m spatial resolution European Environment Agency (EEA) CORINE (Coordination of Information on the Environment) land cover (CLC) dataset, version 12/2009, was employed.As shown in Table 1, all urban elements of the CLC dataset were first extracted and remapped to the corresponding IGBP land use categories.The remapped 250 m resolution CLC urban land use elements were then overlaid with the 30-arc-sec resolution IGBP land use grid.Finally, each IGBP grid cell was assigned to the maximum percentage CLC urban land use category, as illustrated in Fig. 3.
At a second stage, the detailed urban land use information was exploited for specifying new values for the critical parameters of surface albedo and emissivity, and roughness length.These values were defined per CLC urban land use category using satellite remote sensing data collected in the frame of the "Urban Heat Islands and Urban Thermography" project (hereafter "the UHI project"), funded by the European Space Agency (ESA).The description of the methodology applied for deriving these parameters is beyond the scope of this paper, but it can be found in Viel et al. (2011).The satellite-retrieved surface albedo and emissivity, and roughness length values were corresponded to the IGBP urban land use categories by applying the mapping scheme presented in Table 1.For clarity, Table 2 highlights the values that were specified for the city of Athens, Greece.Introduction

Conclusions References
Tables Figures

Back Close
Full

The downscaling mask
Air temperature variations due to the UHI effect can be significant on the scale of a few hundreds of meters.Therefore, the numerical simulation of heat islands requires in principle the adoption of a high horizontal grid resolution (e.g.≤ 1 km).However, this is not always possible because of limitations in the availability of computational resources.This is particularly true in operational forecasting systems, as the one presented herein, where forecasts should be always delivered in time.Within this context, downscaling techniques can provide an alternative solution for producing high-resolution forecasts without increasing computational demands.
In the present study, a downscaling mask was developed exploiting the artificial neural network (ANN) technology.The principle idea was to design and train an ANN to predict air temperature at high spatial resolution.For this, the Tiberius Data Mining software (version 6.1.9,http://www.tiberius.biz/),a predictive modeling tool allowing for the design and implementation of ANNs, was selected for setting up a high-resolution temperature predictive neural network.The feed-forward architecture was chosen for designing the ANN, based on literature review and a preliminary trial and error procedure for various types of neural networks.In this case, all of the data information flows in one direction.The neurons of one layer are connected to the neurons of the following layer using adjustable weights and there is no feedback (Fig. 4).The initial weights of the neural network were randomly specified.
As seen in Fig. 4, the designed ANN consists of one hidden layer with seven neurons, followed by one output layer with a single neuron.The back propagation learning algorithm was used for training the network, while the tangent sigmoid function was used for transferring data from one layer to the other.The input layer consists of three neurons, responsible for feeding the ANN with the input data.ture data, bilinearly interpolated onto the locations of the measurement sites that were used during the training of the ANN.It should be noted that prior to the selection of the above proxy data, several other combinations of input data were tested without inducing significant improvements on the network's performance.The downscaling neural network was trained and tested for the city of Athens, Greece, during a one-month period in the summer of 2009.The training/testing methodology and the corresponding results are described in detail in Giannaros (2013).
In brief, this procedure resulted to a statistical, non-linear downscaling equation of the following form: where AT(i c ) and LU(i c ) are the air temperature and land use category of the i c -th CLC grid point, respectively, MAT is the model simulated air temperature, bilinearly interpolated onto the location of the i c -th CLC grid point, and h is the time of the day, transformed using the cosine function.Equation ( 1) is applied as a post-processing step after the implementation of the WRF/Noah modelling system (Fig. 2), resulting to an air temperature field whose spatial resolution equals that of the CLC grid (i.e.250 m).
4 Numerical simulations and evaluation methods

Modelling system configuration
The UHI-MS was implemented operationally during a 1.5-month period in the summer of 2010, spanning from 13 July to 27 August 2010.Three 1-way nested domains with horizontal grid resolutions of 30 km (d01; mesh size of 199 × 175), 10 km (d02; mesh size of 214 × 175), and 2 km (d03; mesh size of 75 × 75) were defined (Fig. 1c).The outermost 30 km domain covers most of the European continent and is used for simulating large-scale meteorological conditions.The inner two domains (d02 and d03) are used for simulating meso-scale and local-scale meteorological patterns, of which Introduction

Conclusions References
Tables Figures

Back Close
Full the innermost domain (d03) focuses on the study area.All modeling domains have 33 layers in the vertical dimension.The lowest layer is found at approximately 10 m above ground level, while the model top is defined at 100 hPa with radiative boundary conditions for all domains.
The WRF single-moment six-class scheme (WSM6) was used for the parameterization of microphysical processes (Hong et al., 2004).Long-wave radiation processes were parameterized using the Eta Geophysical Fluid Dynamics Laboratory (GFDL) scheme (Schwarzkopf and Fels, 1991) while the MM5 scheme (Dudhia, 1989) was adopted for short-wave radiation.The MM5 similarity scheme (Zhang and Anthes, 1982) was selected for parameterizing the surface layer.Planetary boundary layer (PBL) processes were handled using the Yonsei University (YSU) scheme (Hong et al., 2006) and the Kain-Fritsch scheme (Kain, 2004) was used for parameterizing cumulus convection.Land surface processes were parameterized by coupling the WRF model with the Noah LSM (also refer to Sect.3).
The numerical simulations were initialized on a daily basis using the operational 0.5 • × 0.5 • spatial resolution and 6 h temporal resolution NCEP/GFS data.The lateral boundary conditions for the outermost modeling domain (d01) were obtained by linearly interpolating the 6-hourly NCEP/GFS data, while for the innermost domains (d02 and d03) the lateral boundary conditions were formed through interpolation from the parent domain (d01).Soil moisture and temperature data were also initialized from the NCEP/GFS data.
For the previously defined 1.5-month period, 84 h numerical simulations were conducted operationally, on a daily basis, providing model output at 1 h intervals.The first 12 h of each simulation were discarded as the warm-up period of the model.The evaluation of the modeling system was conducted using the first 24 h of the remaining 72 h of each daily simulation.Introduction

Conclusions References
Tables Figures

Back Close
Full

Evaluation methods
The performance of the UHI-MS was evaluated using in situ data collected from a network of eight weather stations installed and operated in the GAA by the Hydrological Observatory of Athens (HOA) of the National and Technical University of Athens (NTUA).The locations of the stations are shown in Fig. 1b.Table 3 summarizes the main characteristics of the measurement sites.
The conducted evaluation focuses on the near-surface air temperature, defined at 2 m above ground level (a.g.l.) and thermal comfort conditions; the latter being described by two indices.For this, the downscaled 250 m resolution model results were first interpolated onto the locations of the measurement sites using bilinear interpolation.The interpolated forecast data were then compared with the observations.The statistical measures that were computed to quantitatively evaluate the performance of the modeling system include: 1. the mean bias error (MBE), and, 2. the root mean squared error (RMSE), (3) where letters "f " and "o" denote modeled and observed data, respectively.Introduction

Conclusions References
Tables Figures

Back Close
Full

Near-surface air temperature
Figure 5 presents a summary of the overall performance of the modeling system during the two-month operational implementation.It can be seen that the UHI-MS tends to overestimate near-surface air temperature, showing a warm bias on approximately 60 % of the examined cases.Nevertheless, model errors are found to rarely exceed 3 • C, while the majority of them (∼ 80 %) lies in between −2 • C and 2 • C.This results to a low MBE of 0.45 • C and a RMSE that is lower than 2 • C (Fig. 5a).Moreover, there appear to be no significant differences in the modeling system's performance over the different land cover types, as seen in Fig. 5b.The highest MBE values are found for the DUF and URB sites, although not exceeding 1 • C (0.91 • C and 0.81 • C, respectively).On the other hand, the computed RMSE values show little variation, ranging from 1.43 • C (IOC) to 1.99 • C (DUF).Putting this information together, it can be claimed that, in overall, the UHI-MS performs well in simulating the near-surface air temperature field over the study area.
Figure 6 shows the time series of observed and modeled daily mean, maximum and minimum near-surface air temperature, averaged over all measurement sites.The UHI-MS is found to simulate well the day-to-day variations of temperature during the entire study period.In particular, it captures successfully the prolonged period of high temperatures, from 13 to 20 August 2010 (observed and modeled daily mean temperatures exceeding 30 • C), as well as the abrupt cooling of the study area on 21 August 2010 (Fig. 6a).It is also evident that maximum temperatures are generally underestimated by approximately 1.0-1.5 • C (Fig. 6b), whereas minimum temperatures are systematically overestimated by about 1.5-2 • C (Fig. 6c).
The UHI-MS is also found to perform well on the average diurnal course (Fig. 7).The excess warmth of the urbanized sites (Fig. 7a-c) relative to the rural ones (Fig. 7d), indicating the existence of the UHI, is generally simulated well.The timing of observed and modeled temperature minima and maxima also coincide well.The modeling sys-Introduction

Conclusions References
Tables Figures

Back Close
Full tem captures successfully the magnitude of the average daily maximum temperature, showing absolute errors that do not exceed 1 • C. On the other hand, the errors of the modeled average daily minimum temperature are generally larger, exceeding 2 • C in the case of DUF sites (Fig. 7a).In general, it can be seen that the UHI-MS performs better during the daytime (07:00-20:00 LT, LT=UTC+3), than during the night (21:00-06:00 LT).Indeed, daytime model errors vary, on average, from −0.46 • C (IOC, Fig. 7b) to 0.56 • C (URB, Fig. 7c).During the night, the respective values lie in between 0.08 • C (RUR, Fig. 7d) and 1.98 • C (DUF, Fig. 7a), indicating that nighttime temperatures are overestimated.Looking at Fig. 7, one can finally notice that the modeling system performs remarkably better over the RUR sites (Fig. 7d) than over the urbanized land cover types (Fig. 7a-c).

Thermal comfort
When high temperatures are combined with low wind speeds and high humidity, the population may be put under severe thermal stress.The deterioration of human thermal comfort conditions in urban areas is a well-known and documented problem, potentially associated with the UHI effect (e.g.Giannaros and Melas, 2012;Poupkou et al., 2011).Accurate forecasting of urban thermal comfort conditions is thus critical for limiting the adverse impacts of heat waves and heat islands on human health.Such forecasts could be used for both informing the public of a city and developing successful shortterm mitigation strategies.
In the present study, the ability of the UHI-MS to simulate/forecast thermal comfort conditions over the study area is evaluated using two biometeorological indices: (a) the discomfort index (DI), and (b) the approximated wet-bulb globe temperature (AWBGT).DI, originally proposed by Thom (1959), is a simple and widely used index of thermal comfort.It can be computed using Eq. ( 3 AWBGT is an index that can be used for evaluating outdoor thermal stress conditions.It can be computed using Eq. ( 4) (Steeneveld et al., 2011): where e is the water vapor pressure (hPa).Contrary to DI, the threshold values of AWBGT depend largely on personal clothing and level of activity.Nevertheless, the computed AWBGT values can be interpreted using a more general classification scheme as proposed by Steeneveld et al. (2011).AWBGT values lower than 27.7 • C (AWBGT < 27.7 • C) represent absence of heat stress conditions.For 27.7 • C < AWBGT<32.2 • C heat stress increases, while for AWBGT > 32.2 • C great heat stress danger occurs.
Observed and modeled hourly DI and AWBGT values were computed applying Eqs. ( 3) and (4), respectively.The water vapor pressure data (e), needed for the calculation of AWBGT, were derived using the air temperature and relative humidity data.
First, the saturation vapor pressure was estimated, applying the Goff-Gratch equation (WMO, 2000); the vapor pressure was subsequently computed using the observed and modeled relative humidity.
Figure 8 presents the time series of observed and modeled DI and AWBGT, averaged over all measurement sites.Clearly, the UHI-MS reproduces successfully the phase of the diurnal cycle of both indices, as well as the day-to-day variations.For instance, it can be seen that it captures the low DI and AWBGT values during the period from 27 to 31 July 2010.The increased values during the period from 17 to 20 August Introduction

Conclusions References
Tables Figures

Back Close
Full 2010 are also well reproduced.However, it can be noted that the modeling system performs better in simulating maximum DI and AWBGT values, but overestimates the minimum values of both indices.This can be attributed to the overestimation of nighttime minimum temperatures, as previously discussed (Sect.5.1).This overestimation is found to be larger in the case of AWBGT (Fig. 8b), resulting to larger MBE and RMSE values than for DI (Fig. 8a).Tables 4 and 5 present the frequency distributions of the observed and modeled hourly DI and AWBGT values, respectively, within the specified ranges, as a function of the land cover category.It is clear that the UHI-MS underestimates the occurrences of the lowest DI (Table 4) and AWBGT (Table 5) values, but simulates adequately well the highest ones.As seen in Table 4, the occurrences of DI values in the third class ((24, 27]) are systematically overestimated; the overestimation being greater for the DUF and URB sites.On the other hand, the frequencies of the fourth DI class ((27, 29]) are, more or less, successfully reproduced over all land cover categories.The model-observations agreement as regards the variation of the occurrences of DI ranges among the different land cover categories is also found to be good.In particular, DUF and URB sites exhibit the highest observed and modeled frequencies for DI > 24 • C, while RUR sites show the highest observed and modeled frequencies for As for AWBGT (Table 5), the modeling system shows good skill in simulating values in the third class (> 32.

Conclusions
The development of a heat island modeling system, which could be used in the context of operational real-time weather forecasting applications, is presented in the current paper.The performance of the modeling system is evaluated for a 1.5-month period in the summer of 2010, when it was implemented operationally, on a daily basis, for the city of Athens, Greece.The evaluation procedure focuses on the near-surface air temperature, which is the variable that is directly influenced by the UHI effect.In addition, the ability of the modeling system in terms of simulating thermal comfort conditions is evaluated, using two appropriate indices.
Comparisons made with observations revealed an overall satisfactory performance of the modeling system.The following points summarize the findings of this study: 1.The modeling system simulates well the near-air temperature field over the study area, showing MBE and RMSE values of 0.45 • C and 1.75 • C, respectively.Among the different land cover types, model performance metrics show little variation.
2. Daily maximum temperatures are more accurately modeled, exhibiting a cold bias that does not exceed 1.5 • C. On the other hand, the cold bias of daily minimum temperatures is found to be larger, ranging from 1.5 • C to 2.0 • C.
3. On the average diurnal course, the modeling system is found to perform better during the day than during the night.
4. The two thermal comfort indices are simulated adequately well.However, the modeling system appears to be more able to capture the highest values of both indices, rather than the lowest ones.In summary, our results indicate that given the necessary modifications and adaptations, meso-scale meteorological models constitute a promising tool for the operational simulation/forecasting of the urban thermal climate and thermal comfort conditions.Nevertheless, there is still room for improvement of the presented modeling system.For this, the model performance should be further evaluated, focusing explicitly on the heat island effect and its spatial and temporal variations.This remains to be carried out in the future, together with testing the modeling system over another urban area.Introduction

Conclusions References
Tables Figures

Back Close
Full  Full  Full  Full and well-documented environmental problem associated with urbanization.It is a phe-Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Figure 2
Figure2presents the flow chart of the modelling system (hereafter "the UHI-MS") that was developed for the operational simulation of the UHI effect.The modelling system consists of three major parts: (1) the meteorological WRF model, (2) the Noah land surface model (LSM), and (3) a statistical downscaling mask, developed to increase spatial resolution.The input of the UHI-MS involves static terrestrial data, such as topography and land use/cover, and land surface data, such as albedo and emissivity.
Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | The latter include: (1) the time of the day, transformed using the cosine function, to represent the diurnal cycle of air temperature, (2) the high-resolution CLC land use data to represent the effect of land surface type on air temperature, and (3) the WRF/Noah simulated air tempera-Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | ): DI = T a − 0.55 • (1 − 0.01 • RH) • (T a − 14.5Discussion Paper | Discussion Paper | Discussion Paper | where T a is the air temperature ( • C) and RH the relative humidity (%).DI values below or equal to 21 • C (DI ≤ 21 • C) correspond to absence of thermal discomfort.When 21 • C < I ≤ 24 • C less than half the population is expected to feel discomfort, while when 24 • C < DI ≤ 27 • C the percentage of the population feeling discomfort rises to 50 %.The majority of the population is anticipated to feel discomfort when 27 • C < DI ≤ 29 • C, and the entire population is feeling discomfort when 29 • C < DI ≤ 32 • C. DI values exceeding 32 • C (DI > 32 • C) correspond to sanitary emergency conditions.
Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 2 • C).The occurrences of AWBGT values in the second class ([27.7,32.2]) are generally overestimated; the overestimation being again larger for the DUF sites.Nevertheless, the UHI-MS manages to reproduce adequately well the differences between the examined land cover categories.Most occurrences of AW-BGT < 27.7 • C are successfully modeled for the RUR sites, contrary to DUF and URB sites that both exhibit most occurrences of observed and modeled AWBGT > 32.2 • C. Discussion Paper | Discussion Paper | Discussion Paper |

5.
The observed and modeled frequencies of occurrences for the various thermal comfort and heat stress classes agree reasonably well.This is particularly clear in the case of the classes representing more thermally deteriorated conditions.Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Figure 1 . 6 Fig. 1 .Figure 2 .Fig. 2 .Figure 6 .Fig. 6 .
Figure 1.Topography of (a) Greece, with identification of the city of Athens, and (b) GAA, 3 with identification of the locations of the measurement sites.(c) Configuration of the three 4 one-way nested domains used for the WRF simulations.5 6 Fig. 1.Topography of (a) Greece, with identification of the city of Athens, and (b) GAA, with identification of the locations of the measurement sites.(c) Configuration of the three one-way nested domains used for the WRF simulations.

Table 1 .
Summary of the mapping scheme between the CLC urban land use elements and IGBP (WRF) urban land use categories.

Table 4 .
Frequency distributions (%) of observed and modeled hourly DI values within the specified ranges, grouped by land use category.

Table 5 .
Same as Table 4, but for AWBGT.