Land Surface Temperature Differences within Local Climate Zones, Based on Two Central European Cities

The main factors influencing the spatiotemporal variability of urban climate are quite widely recognized, including, for example, the thermal properties of materials used for surfaces and buildings, the mass, height and layout of the buildings themselves and patterns of land use. However, the roles played by particular factors vary from city to city with respect to differences in geographical location, overall size, number of inhabitants and more. In urban climatology, the concept of “local climate zones” (LCZs) has emerged over the past decade to address this heterogeneity. In this contribution, a new GIS-based method is used for LCZ delimitation in Prague and Brno, the two largest cities in the Czech Republic, while land surface temperatures (LSTs) derived from LANDSAT and ASTER satellite data are employed for exploring the extent to which LCZ classes discriminate with respect to LSTs. It has been suggested that correctly-delineated LCZs should demonstrate the features typical of LST variability, and thus, typical surface temperatures should differ significantly among most LCZs. Zones representing heavy industry (LCZ 10), dense low-rise buildings (LCZ 3) and compact mid-rise buildings (LCZ 2) were identified as the warmest in both cities, while bodies of water (LCZ G) and densely-forested areas (LCZ A) made up the coolest zones. ANOVA and subsequent multiple comparison tests demonstrated that significant temperature differences between the various LCZs prevail. The results of testing were similar for both study areas (89.3% and 91.7% significant LST differences for Brno and Prague, respectively). LSTs computed from LANDSAT differentiated better between LCZs, compared with ASTER. LCZ 8 (large low-rise buildings), LCZ 10 (heavy industry) and LCZ D (low plants) are well-differentiated zones in terms of their surface temperatures. In contrast, LCZ 2 (compact mid-rise), LCZ 4 (open high-rise) and LCZ 9 (sparsely built-up) are less distinguishable in both areas analyzed. Factors such as seasonality and thermal anisotropy remain a challenge for future research into LST differences.


Introduction
The main factors influencing the spatiotemporal variability of urban climate are quite widely recognized, including, for example, the thermal properties of materials used for surfaces and buildings, the mass, height and layout of the buildings themselves and patterns of land use (e.g., [1,2]). However, urban climate research has long been limited by unsatisfactory quantification of these factors and non-standardized descriptions of them [3]. To address these problems in urban climatology, the concept of local climate zones (LCZs) has emerged over the past decade [4]. Stewart and Oke [4] define local climate zones as "regions of uniform land cover, surface structure, construction material and human activity that span hundreds of meters to several kilometers on a horizontal scale". The thermal climate of each LCZ is determined by its surface character: structure (building and tree height/density), cover (permeability), fabric (albedo, thermal admittance) and metabolism (waste heat from transportation and space heating/cooling). Unique combinations of these properties provide a distinctive thermal regime for each LCZ, particularly the characteristic temperature regime at screen height, best observed in clear weather conditions and areas of simple relief [4,5].
The concept of the LCZ emerged out of the need to improve the documentation of atmospheric heat island observations, but its relation to spatial and temporal temperature variability has recently become a key research question [4,6,7]. Most recent studies investigating urban temperature fields utilize the LCZ classification and largely support the correspondence of LCZs with air temperature fields in cities and their surroundings (e.g., [6,[8][9][10]). However, a number of authors working with LCZs have pointed out that the influence of thermal, radiative, metabolic and surface-cover properties on the formation of local climate may significantly vary with respect to the geographical location of the zone, the size of the city, the position within the city (distance from the city center, neighboring zones) and relief [11][12][13]. Simultaneously, Bechtel and Daneke [14] and Lelovics et al. [9] created the first LCZ mapping methods and applied the LCZ concept to secondary uses. With such an extension of LCZ beyond its primary function, certain new methodological problems appeared: the size of a spatial mapping unit, the method used for generalization, the temporal variability of the physical properties of the environment, the objectification and standardization of the classification procedure, and others [8,9,14]. The extension of LCZ to regional mapping significantly widens the field of LCZ applications. For instance, Alexander et al. [15] revealed distinctive heat fluxes for some LCZs and concluded that LCZ classification provides a useful sampling framework for the derivation of the land-cover fractions needed to run surface-energy balance models. Zuvela-Aloise et al. [16] and Geletič et al. [17] used LCZs as input to MUKLIMO_3, a non-hydrostatic micro-scale air temperature model. Further, Skarbit et al. [18,19] worked on the relations between LCZ and LST. The authors show that the different LCZ classes have different surface temperature characteristics as for their typical values and their frequency distribution. The latter will be addressed in this study in more detail because the statistical significance of the differences between the typical LST of climate zones will be analyzed.
From the nature of LCZ classification, standard or special-purpose near-surface air temperature measurements may be best suited to analyzing temperature contrasts among LCZs [6]. However, such measurements, even supplemented by measurements from a dense network of automatic stations, often prove insufficient to describe and understand fully the spatial variability of air temperature in urban areas [17,20]. These measurements are insufficient especially in terms of their ability to characterize temperature spatial distribution in detail, as the result may be dependent on the number of "point" measurements, the method of interpolation and more. Thus LSTs derived from various airborne or satellite remote sensing systems may be a suitable alternative, as their spatial coverage is complete. Moreover, spatial and temporal resolution is constantly improving [21][22][23]. In spite of the fact that surface temperature modulates the air temperature of the boundary layer atmosphere [21], relationships between the two types of temperature are quite complicated in urban environments [24][25][26]. Apparent urban LSTs are modified by numerous static (e.g., land cover type, urban morphology) and dynamic (e.g., geometry of thermal imagery acquisition, solar elevation) factors. These factors not only modify the energy balance of urban areas [1,4], but they also generate effects, such as thermal anisotropy [27]. Nevertheless, the use of land surface temperature for the analysis of temperature spatial variability and the estimation of urban heat island intensity remains challenging.
From all of this, it becomes clear that the LST characteristics typical of LCZs (if such definitive factors indeed exist) should be characterized in more detail in order to understand the influence of local climate zones on local climate formation better. Case studies using remote sensing data to investigate land surface temperatures (LSTs) from various cities may therefore contribute significantly to the study of LCZs. The particular aim of this study is to analyze relations between surface temperatures and local climate zones, taking two cities in the Czech Republic as examples. A new GIS-based method is used for LCZ delimitation. Compared to the widely-used method of Bechtel and Daneke [14], this method does not use thermal imagery, which is an important aspect as the paper explores the extent to which LCZ classes are delineated with respect to land surface temperatures.

Study Area
This study focuses on Prague and Brno, the two largest cities in the Czech Republic. The capital Prague, with around 1.3 million inhabitants, as well as Brno, with fewer than 0.4 million inhabitants, may be considered as medium-sized cities typical of the central European region (Figure 1, Table 1). From all of this, it becomes clear that the LST characteristics typical of LCZs (if such definitive factors indeed exist) should be characterized in more detail in order to understand the influence of local climate zones on local climate formation better. Case studies using remote sensing data to investigate land surface temperatures (LSTs) from various cities may therefore contribute significantly to the study of LCZs. The particular aim of this study is to analyze relations between surface temperatures and local climate zones, taking two cities in the Czech Republic as examples. A new GIS-based method is used for LCZ delimitation. Compared to the widely-used method of Bechtel and Daneke [14], this method does not use thermal imagery, which is an important aspect as the paper explores the extent to which LCZ classes are delineated with respect to land surface temperatures.

Study Area
This study focuses on Prague and Brno, the two largest cities in the Czech Republic. The capital Prague, with around 1.3 million inhabitants, as well as Brno, with fewer than 0.4 million inhabitants, may be considered as medium-sized cities typical of the central European region (Figure 1, Table 1).  In terms of urban morphology and relief, the cities are quite similar. Brno and Prague are both located in hilly landscapes. Altitudes vary from under 200 m to over 500 m in Brno and under 200 m to 400 m in Prague (mean elevations are 259 m for Brno and 288 m for Prague). Both cities have grown in concentric zones, reflecting the historical stages of their development. Building patterns show this clearly, and neither city lost significant proportions of their architectural heritage to the Second World War (clearly-defined historical center, residential buildings old and new, industrial areas, housing estates, modern shopping centers, malls and allotments). Because of this, neighborhoods with irregular, as well as regular (grid and radial grid) street plans are present in both cities. Differences between the cities are especially the presence of a large river (the Vltava) in Prague and a higher proportion of open vegetated (low plant) spaces in Brno. The landscapes beyond the city's boundaries  In terms of urban morphology and relief, the cities are quite similar. Brno and Prague are both located in hilly landscapes. Altitudes vary from under 200 m to over 500 m in Brno and under 200 m to 400 m in Prague (mean elevations are 259 m for Brno and 288 m for Prague). Both cities have grown in concentric zones, reflecting the historical stages of their development. Building patterns show this clearly, and neither city lost significant proportions of their architectural heritage to the Second World War (clearly-defined historical center, residential buildings old and new, industrial areas, housing estates, modern shopping centers, malls and allotments). Because of this, neighborhoods with irregular, as well as regular (grid and radial grid) street plans are present in both cities. Differences between the cities are especially the presence of a large river (the Vltava) in Prague and a higher proportion of open vegetated (low plant) spaces in Brno. The landscapes beyond the city's boundaries are of a predominantly agricultural character, with agricultural areas represented by low plants (cultivated plants like maize, barley, wheat or rapeseed; growing period beginning in April and harvest in August; fields not irrigated). Patches of orchards (LCZ B) and mixed forests (LCZ A) complete the landscape mosaics of the hinterlands (Brno also has large forested areas north of the city).
Mean annual temperature stands at 9.1 • C for Brno and 8.4 • C for Prague (1961-2000 reference period). The annual temperature regime is characterized by a single maximum (the warmest month being July with mean temperatures of 19.3 • C in Brno and 18.2 • C in Prague), with a single minimum (the coldest month is January with a mean temperature −1.9 • C in Brno and −1.4 • C in Prague). Mean annual precipitation totals are quite similar in both (505 mm for Brno and 526 mm for Prague) [28,29].

Local Climate Zones
A GIS-based approach was taken for the delimitation of LCZs [30]. This method had already been developed and tested in the area of Brno and its surroundings. It is based on the measurable physical properties of the environment and a clearly-defined decision-making algorithm. The algorithm derives from the basic physical parameters defined by Stewart and Oke [4]: building surface fraction (BSF), pervious surface fraction (PSF), impervious surface fraction (ISF) and height of roughness elements (HRE). These were supplemented by number of buildings (NoB). For the classification process, the study area was divided into a regular grid with cells of 100 m × 100 m, each cell exhibiting these physical parameters.
In the first step, cells were divided by BSF into two basic classes: building types and natural land cover types. In the second step, the cells in the first step categorized as building types (BSF ≥ 10%) were classified into particular LCZs on the basis of the smallest deviation from the optimum interval, as defined by Geletič and Lehnert [30] (against the background of values suggested by Stewart and Oke [4]; see Table 2) in terms of physical parameters (BSF, PSF, ISF and HRE). Industrial zones (LCZs 8 and 10) were refined using NoB. Cells that had been categorized as land-cover types in the first step were automatically classified into particular LCZs in the second step by means of the vector land use geodatabase [30]. When all cells had been assigned to an appropriate LCZ, the LCZ areas were delineated using a majority filter to smooth the results of classification. Independent testing of the decision-making algorithm for defining the percentage coverage for individual LCZs demonstrated close agreement with areas defined on the basis of expert knowledge (e.g., there was 90% agreement for Brno, not shown here). New LCZ maps for Brno and Prague were created for this particular study.
Altogether, 16 thermal images, eight for each city, four from LANDSAT-8 and four from ASTER, were used for the comparison with LCZs. Basic scene information is summarized in Table 3. For the interpretation of further analyses carried out in this study, it is important to emphasize that all thermal images were recorded in the morning daylight hours (except 19 July). Two algorithms were used for the estimation of LST: split-window and multispectral. The split-window technique uses two thermal infrared (TIR) bands, typically located in the atmospheric window between 10 µm and 12 µm [31]. In this study, the algorithm applied was based on the estimation of top-of-atmosphere spectral radiance and at-satellite brightness temperature. LANDSAT-8 provides two thermal bands, Band 10 and Band 11 [32]. However, LANDSAT-8 does not provide data for surface emissivity calculation. Because of this, a land-surface emissivity (LSE) algorithm was used, estimating emissivity from the normalized difference vegetation index (NDVI) [33]. The scheme of LANDSAT-8 scene processing is summarized in Figure 2. The thermal infrared bands on ASTER are located between 8 μm and 12 μm [34]. Two different approaches to the calculation of LST can be taken. The first exploits Level-1B data, which are produced by applying radiometric calibration and geometric correction coefficients. The multispectral algorithm for ASTER multiplies thermal bands by a scale factor, then atmospheric corrections are applied. The normalized emissivity (NEM) algorithm was used for LST calculation [34]. A schematic representation of this method appears in Figure 3. The ASTER product AST_08 (L2 surface kinetic temperature) may also be downloaded, but this option has only been available since 1 April 2016.  The thermal infrared bands on ASTER are located between 8 µm and 12 µm [34]. Two different approaches to the calculation of LST can be taken. The first exploits Level-1B data, which are produced by applying radiometric calibration and geometric correction coefficients. The multispectral algorithm for ASTER multiplies thermal bands by a scale factor, then atmospheric corrections are applied. The normalized emissivity (NEM) algorithm was used for LST calculation [34]. A schematic representation of this method appears in Figure 3. The ASTER product AST_08 (L2 surface kinetic temperature) may also be downloaded, but this option has only been available since 1 April 2016. The thermal infrared bands on ASTER are located between 8 μm and 12 μm [34]. Two different approaches to the calculation of LST can be taken. The first exploits Level-1B data, which are produced by applying radiometric calibration and geometric correction coefficients. The multispectral algorithm for ASTER multiplies thermal bands by a scale factor, then atmospheric corrections are applied. The normalized emissivity (NEM) algorithm was used for LST calculation [34]. A schematic representation of this method appears in Figure 3. The ASTER product AST_08 (L2 surface kinetic temperature) may also be downloaded, but this option has only been available since 1 April 2016.

Comparison of Land-Surface Temperatures in Local Climate Zones
The analysis herein is based on the assumption that individual LCZs should demonstrate certain features typical of a given LST regime. LST fields were overlaid with LCZs, and typical LSTs were calculated for each zone. Differences between mean LCZ temperatures were evaluated by one-way analysis of variance (ANOVA). Before applying ANOVA, a careful assessment of assumptions for the application of this method was made. This comprised, in particular, testing for normality using the Kolmogorov-Smirnov tests, Q-Q plots and comparisons of LCZ temperature variability. When the ANOVA F-test indicated statistically-significant differences in LST, the Tukey HSD test was employed to control for the effect of multiple comparisons [35] and to reveal which LCZs differentiate significantly as to their mean LST. The analysis was performed separately for all 16 temperature fields constructed (two study areas with eight temperature fields each). Finally, the score counting the number of tests indicating significant temperature differences was calculated for each LCZ. Interpretation of this score is straightforward: LCZs with higher scores are better differentiated from other zones in terms of their typical LSTs. The other way round, LCZs with lower scores are less distinguishable from other zones in terms of their typical LSTs.

Local Climate Zones
The LCZ distributions exhibit concentric patterns, especially for Prague. The historical centers of Brno and Prague form compact areas of LCZ 2 (Figures 4 and 5), clearly surrounded by LCZ 5, the most frequent type of development in both cities. Towards the outer boundaries of compact development, there is an increase in areas classified as LCZ 6. Together with LCZ 9, LCZ 6 is also typical of the villages around the cities. An increased occurrence of LCZ B is typical of the boundary between compact development and its surrounding landscape. LCZ B could, however, be found even near the centers of both cities (parks). LCZs 8, 10 and E appear especially on the outer part of the cities; these zones are frequently present in projections of development (city) into the surrounding countryside. Nevertheless, there are also the sites of LCZ 8, 10 and E in the core areas of both Brno ( Figure 4) and Prague ( Figure 5). LCZ D dominates in Prague's surroundings, as it does in the landscape south of Brno, but to the north of Brno lie large areas of LCZ A. Also in the latter area is a large body of water, the Brno town reservoir, which creates a 259-ha area of LCZ G. Similarly, a line of LCZ G marks the course of the River Vltava in Prague.

Comparison of Land-Surface Temperatures in Local Climate Zones
The analysis herein is based on the assumption that individual LCZs should demonstrate certain features typical of a given LST regime. LST fields were overlaid with LCZs, and typical LSTs were calculated for each zone. Differences between mean LCZ temperatures were evaluated by one-way analysis of variance (ANOVA). Before applying ANOVA, a careful assessment of assumptions for the application of this method was made. This comprised, in particular, testing for normality using the Kolmogorov-Smirnov tests, Q-Q plots and comparisons of LCZ temperature variability. When the ANOVA F-test indicated statistically-significant differences in LST, the Tukey HSD test was employed to control for the effect of multiple comparisons [35] and to reveal which LCZs differentiate significantly as to their mean LST. The analysis was performed separately for all 16 temperature fields constructed (two study areas with eight temperature fields each). Finally, the score counting the number of tests indicating significant temperature differences was calculated for each LCZ. Interpretation of this score is straightforward: LCZs with higher scores are better differentiated from other zones in terms of their typical LSTs. The other way round, LCZs with lower scores are less distinguishable from other zones in terms of their typical LSTs.

Local Climate Zones
The LCZ distributions exhibit concentric patterns, especially for Prague. The historical centers of Brno and Prague form compact areas of LCZ 2 (Figures 4 and 5), clearly surrounded by LCZ 5, the most frequent type of development in both cities. Towards the outer boundaries of compact development, there is an increase in areas classified as LCZ 6. Together with LCZ 9, LCZ 6 is also typical of the villages around the cities. An increased occurrence of LCZ B is typical of the boundary between compact development and its surrounding landscape. LCZ B could, however, be found even near the centers of both cities (parks). LCZs 8, 10 and E appear especially on the outer part of the cities; these zones are frequently present in projections of development (city) into the surrounding countryside. Nevertheless, there are also the sites of LCZ 8, 10 and E in the core areas of both Brno ( Figure 4) and Prague ( Figure 5). LCZ D dominates in Prague's surroundings, as it does in the landscape south of Brno, but to the north of Brno lie large areas of LCZ A. Also in the latter area is a large body of water, the Brno town reservoir, which creates a 259-ha area of LCZ G. Similarly, a line of LCZ G marks the course of the River Vltava in Prague.

Land Surface Temperature
Brno and Prague show similar spatial distribution of surface temperatures. Industrial zones and the city centers with historical buildings were always warmer compared to their surroundings (Figures 4 and 5). Away from the densely built-up areas, the warm parts of the cities also occurred near large patches of relatively flat, impervious surfaces (car parks, railway marshalling yards, stations, etc.). Hotspots were very often associated with large commercial and distribution areas and the international trade-fair grounds (which are quite extensive) regardless of whether they lay in the core of the city or in the countryside. On the other hand, bodies of water (the Brno reservoir and Prague's River Vltava) and forested areas formed the coldest localities (see Figures 4 and 5).

Local Climate Zones and Land Surface Temperatures
The spatial distributions of LCZs and LST presented above indicate that there is a certain degree of correspondence. Further analyses were therefore carried out to investigate their relations.
Box-plots summarizing typical LSTs for individual LCZs show relatively consistent results for the two cities (Figures 6 and 7). However, the LSTs for LCZs derived from individual scenes cannot be further compared directly as they represent different days within a year, a period spanning 20 March-28 September (Table 3). Particular LSTs arose not only out of static factors (e.g., land cover), but were influenced by dynamic elements, such as the synoptic situation on a given day. Therefore, LSTs of LCZs were compared in terms of rank order within a set of eight analyzed scenes for each city. The mean LSTs of LCZs were sorted in descending order for each scene, and mean rankings were calculated from all scenes individually.
The majority highest mean LST emerged for LCZ 10 (heavy industry), which was the warmest for all eight scenes in Prague and for seven scenes in Brno. LCZ 3 (dense mix of low-rise buildings) and LCZ 2 (compact mid-rise buildings) were the second-and third-warmest zones in Brno and Prague. Next, higher LSTs were associated with LCZ 8 (large low-rise) in both study areas. The lowest ranking and, thus, the coolest LCZ was related to bodies of water (LCZ G) followed by LCZ A, areas with dense tree cover. Lower mean LSTs also occurred in other zones where pervious surfaces prevailed: LCZ B (scattered trees) and LCZ C (mixture of shrubs, woody trees and agricultural areas).
The box-plots in Figures 6 and 7 present not only typical (median) LSTs for LCZs, but also permit the evaluation of their variability and the occurrence of outliers. Whereas LCZ temperature variations are more or less comparable and boxes and whiskers symmetrical, indicating no significant departures from normality, outliers typically occur in certain zones; they are analyzed further in the

Land Surface Temperature
Brno and Prague show similar spatial distribution of surface temperatures. Industrial zones and the city centers with historical buildings were always warmer compared to their surroundings (Figures 4 and 5). Away from the densely built-up areas, the warm parts of the cities also occurred near large patches of relatively flat, impervious surfaces (car parks, railway marshalling yards, stations, etc.). Hotspots were very often associated with large commercial and distribution areas and the international trade-fair grounds (which are quite extensive) regardless of whether they lay in the core of the city or in the countryside. On the other hand, bodies of water (the Brno reservoir and Prague's River Vltava) and forested areas formed the coldest localities (see Figures 4 and 5).

Local Climate Zones and Land Surface Temperatures
The spatial distributions of LCZs and LST presented above indicate that there is a certain degree of correspondence. Further analyses were therefore carried out to investigate their relations.
Box-plots summarizing typical LSTs for individual LCZs show relatively consistent results for the two cities (Figures 6 and 7). However, the LSTs for LCZs derived from individual scenes cannot be further compared directly as they represent different days within a year, a period spanning 20 March-28 September (Table 3). Particular LSTs arose not only out of static factors (e.g., land cover), but were influenced by dynamic elements, such as the synoptic situation on a given day. Therefore, LSTs of LCZs were compared in terms of rank order within a set of eight analyzed scenes for each city. The mean LSTs of LCZs were sorted in descending order for each scene, and mean rankings were calculated from all scenes individually.
The majority highest mean LST emerged for LCZ 10 (heavy industry), which was the warmest for all eight scenes in Prague and for seven scenes in Brno. LCZ 3 (dense mix of low-rise buildings) and LCZ 2 (compact mid-rise buildings) were the second-and third-warmest zones in Brno and Prague. Next, higher LSTs were associated with LCZ 8 (large low-rise) in both study areas. The lowest ranking and, thus, the coolest LCZ was related to bodies of water (LCZ G) followed by LCZ A, areas with dense tree cover. Lower mean LSTs also occurred in other zones where pervious surfaces prevailed: LCZ B (scattered trees) and LCZ C (mixture of shrubs, woody trees and agricultural areas).
The box-plots in Figures 6 and 7 present not only typical (median) LSTs for LCZs, but also permit the evaluation of their variability and the occurrence of outliers. Whereas LCZ temperature variations are more or less comparable and boxes and whiskers symmetrical, indicating no significant departures from normality, outliers typically occur in certain zones; they are analyzed further in the Discussion section below. The normality of the LST distribution in each LCZ was further tested more formally, since it is an important assumption for subsequent analysis of variance. A Gaussian curve was fitted to the LST frequency distribution for each zone and differences addressed by the Kolmogorov-Smirnov test. Figure 8 presents an example for a single date. Significant departures from normality occurred for several zones in Brno (LCZ 10 and G) and one in Prague (LCZ 10). Discussion section below. The normality of the LST distribution in each LCZ was further tested more formally, since it is an important assumption for subsequent analysis of variance. A Gaussian curve was fitted to the LST frequency distribution for each zone and differences addressed by the Kolmogorov-Smirnov test. Figure 8 presents an example for a single date. Significant departures from normality occurred for several zones in Brno (LCZ 10 and G) and one in Prague (LCZ 10).    Table 3 for the scene ID; for the explanation, see Figure 6.  Table 3 for the scene ID; for the explanation, see Figure 6. The results of the normality tests and comparable variations in LCZ land-surface temperatures indicated that ANOVA may be used for testing whether the differences in LCZ mean temperatures outlined above are significant or not.
Indeed, the results of the one-way ANOVA F-test (p < 0.001) determined that there exist significant differences among the means of LCZ land surface temperatures. Moreover, these results stand for the two study areas and for all 16 dates analyzed (scenes). Subsequently, Tukey HSD tests revealed which pairs of LCZs significantly differentiate in terms of their mean LSTs (Figures 9 and  10).
Interpretation of all tests for both cities and analyzed scenes is provided in a simple binary format in Figures 9 and 10. If the result of the test indicates no significant difference in average LSTs between two zones in question (p-value higher than 0.05), there is a blue point in the corresponding cell of the table. Blue points represent "negative" result in our analysis ("misses"), that is no difference in average LSTs for a given pair of LCZs. In contrast, if the result of the test shows significant difference (p-value lower than 0.05), there is an empty space in the corresponding cell of the table. Significant differences in average LSTs represent "positive" results in this analysis and may be denoted as "hits".
Specifically, for both study areas, LCZ pairs with statistically-significant differences in mean LCZ temperatures ("hits") prevail, while pairs with no temperature differences ("misses") are less frequent, and such LCZ pairs are different in the individual scenes analyzed. This result strongly supports the method used for LCZ delimitation in both cities. That is, LCZs show typical surface temperatures that are largely quite different from the surface temperatures of other zones. The robustness of this finding is further confirmed by the fact that the relative number of "hits" is comparable for both cities (89.3% for Brno and 91.7% for Prague). Slightly larger differences may be found between the LSTs calculated from the two satellite data sources. Whereas ASTER-based LSTs The results of the normality tests and comparable variations in LCZ land-surface temperatures indicated that ANOVA may be used for testing whether the differences in LCZ mean temperatures outlined above are significant or not.
Indeed, the results of the one-way ANOVA F-test (p < 0.001) determined that there exist significant differences among the means of LCZ land surface temperatures. Moreover, these results stand for the two study areas and for all 16 dates analyzed (scenes). Subsequently, Tukey HSD tests revealed which pairs of LCZs significantly differentiate in terms of their mean LSTs (Figures 9 and 10).
Interpretation of all tests for both cities and analyzed scenes is provided in a simple binary format in Figures 9 and 10. If the result of the test indicates no significant difference in average LSTs between two zones in question (p-value higher than 0.05), there is a blue point in the corresponding cell of the table. Blue points represent "negative" result in our analysis ("misses"), that is no difference in average LSTs for a given pair of LCZs. In contrast, if the result of the test shows significant difference (p-value lower than 0.05), there is an empty space in the corresponding cell of the table. Significant differences in average LSTs represent "positive" results in this analysis and may be denoted as "hits".
Specifically, for both study areas, LCZ pairs with statistically-significant differences in mean LCZ temperatures ("hits") prevail, while pairs with no temperature differences ("misses") are less frequent, and such LCZ pairs are different in the individual scenes analyzed. This result strongly supports the method used for LCZ delimitation in both cities. That is, LCZs show typical surface temperatures that are largely quite different from the surface temperatures of other zones. The robustness of this finding is further confirmed by the fact that the relative number of "hits" is comparable for both cities (89.3% for Brno and 91.7% for Prague). Slightly larger differences may be found between the LSTs calculated from the two satellite data sources. Whereas ASTER-based LSTs indicate 89.2% hits, the LANDSAT-8-based LSTs show 91.9% hits. This implies that LSTs computed from LANDSAT-8 data better differentiate among various LCZs compared to those from ASTER.
A closer look at multiple comparison results facilitates the identification of zones that are well differentiated in terms of mean surface temperatures ( Figure 11). The highest number of "hits" (significant temperature differences) was registered for LCZ 8 (large low-rise), LCZ 10 (heavy industry) and also for LCZ D (low plants). Such results confirm what is generally known of the spatial LST distribution in the areas analyzed. LCZs 8 and 10 consist of areas where impervious surfaces with lower proportions of vegetation cover prevail and where less surface moisture is available, features that are responsible for the higher surface temperatures of these zones. In both Czech cities, industrial areas are well delimited from neighborhoods that consist largely of LCZs with lower surface temperatures. In spite of the fact that LCZ D occurs in rural areas in which low plants prevail, typical land surface temperatures in this zone are markedly higher compared to their surroundings, due to the radiation response of soil (Figures 4 and 5); see the Discussion section for details. LCZs with the lowest number of "hits" (those less distinguishable in terms of typical surface temperatures) may be identified in the same way. They stand for LCZ 2 (compact mid-rise), LCZ 4 (open high-rise) and LCZ 9 (sparsely built-up) in both areas analyzed. While LCZs 4 and 9 are characterized by somewhat lower LSTs, LCZ  A closer look at multiple comparison results facilitates the identification of zones that are well differentiated in terms of mean surface temperatures ( Figure 11). The highest number of "hits" (significant temperature differences) was registered for LCZ 8 (large low-rise), LCZ 10 (heavy industry) and also for LCZ D (low plants). Such results confirm what is generally known of the spatial LST distribution in the areas analyzed. LCZs 8 and 10 consist of areas where impervious surfaces with lower proportions of vegetation cover prevail and where less surface moisture is available, features that are responsible for the higher surface temperatures of these zones. In both Czech cities, industrial areas are well delimited from neighborhoods that consist largely of LCZs with lower surface temperatures. In spite of the fact that LCZ D occurs in rural areas in which low plants prevail, typical land surface temperatures in this zone are markedly higher compared to their surroundings, due to the radiation response of soil (Figures 4 and 5); see the Discussion section for details. LCZs with the lowest number of "hits" (those less distinguishable in terms of typical surface temperatures) may be identified in the same way. They stand for LCZ 2 (compact mid-rise), LCZ 4 (open high-rise) and LCZ 9 (sparsely built-up) in both areas analyzed. While LCZs 4 and 9 are characterized by somewhat lower LSTs, LCZ 2 represents zones with the highest mean surface temperatures; see   Table 3 for the scene ID; blue points indicate pairs of LCZs for which LSTs are not significantly different (p > 0.05), while empty cells denote LCZs pairs for which average LSTs are significantly different.

Discussion
This study is based on satellite (ASTER and LANDSAT-8) observations; therefore, bird's eye-view surface temperatures from nadir (more precisely, nadir and very small off-nadir angles) are analyzed. The results here cannot therefore be equalized with respect to the general thermal responses of any particular LCZ, as vertical surfaces are neglected [20]; the results of this study thus provide a partial, but coherent, insight into the nature of surface temperatures in LCZs. Some noteworthy particularities and queries are addressed previously.
It was mentioned above that the highest LSTs were observed in LCZ 10 (heavy industry). However, when interpreting the thermal behavior of LCZ 10, it must be taken into account that LCZ 10 delineation was based on the Geletič and Lehnert [30] method, which employed a higher percentage of ISF and BSF for delineation of LCZ 10 than that originally suggested by Stewart and Oke [4]. Nevertheless, the higher land surface temperatures of industrial areas in daytime have previously been described by many authors, e.g., [28,36,37]. Unexpectedly, the LSTs of LCZ 8 (large low-rise) were in most cases considerably lower than the LSTs of LCZ 10. However, LCZ 8 was characterized by high statistical variability in LST. This may correspond to the different thermal properties implicit in the large roof surfaces that occur quite often in LCZ 8 [38]. LCZ A (dense trees) was identified as a lower LST zone, but a large number of positive outliers are evident (i.e., highest LST). This is a reaction to the relatively high number of patches, mostly consisting of clear-cuts of areas that are quite often below the spatial resolution of the data used for LCZ classification. On the other hand, the high LST variability of LCZ G (water) should be considered as realistic, since temperature differences reflect the different characters of individual bodies of water (depth, rate of flow, etc.). The LST of LCZ D also exhibited a relatively wide interquartile range with frequent outliers, especially for the upper part of the distribution. Considerable temperature differences between LCZ D and other zones are especially clearly expressed at the beginning and the end of the growing season, periods characterized by lower proportions of vegetation and the stronger radiation response of bare soil (lower albedo, reduced latent heat flux). This corresponds with the findings of [20,39,40] who described different spatial patterns of LST for spring and autumn compared to summer due to seasonal changes in land cover. Therefore, when working with LCZs, seasonal changes in land cover must be taken into consideration.
There were statistically-significant differences in mean LSTs between most LCZs, but LCZs 2, 4 and 9 were recognized as zones less distinguishable from other LCZs. A feature common to these three zones is the higher heterogeneity of land cover typical of central European cities. For instance, LCZ 2 might exhibit a range of building patterns (e.g., street canyons, courtyards) of various orientations to incoming solar radiation. This may significantly influence radiation response regarding both real surface temperature differences and the magnitude of effective anisotropy [14,27]. Moreover, the LSTs of LCZs 2 and 3 are similar because the percentages of PSF, ISF and BSF are also similar and the same types of materials for infrastructure, roofing, etc., are used (at least in central European conditions). One particular factor that differentiates these two zones is the height of the buildings; however, considering the nadir (or very close to nadir) view position of the sensors, the effect of vertical surfaces on apparent surface temperature is neglected, and simultaneously, dense development significantly reduces increase in shaded area (in the sense of direct solar radiation) with building height [41]. This analysis has indicated certain differences between the two sensor data sources (ASTER and LANDSAT-8). There may be several reasons for this. Each span leads to different LST retrieval algorithms to be used for LST calculation from ASTER and LANDSAT, arising out of the varying viewing geometry (IFOV (instantaneous field of view)) of the two satellite systems. In addition to this, land surface temperatures may be biased by several other factors arising from the nature of the remote sensing method. A higher proportion of vertically-oriented surfaces in urban areas compared to rural environments causes uneven solar heating of those surfaces and induces a thermal anisotropy effect [41]. Different solar heating of vertical surfaces arising out of variations in viewing geometry renders the interpretation of land surface temperatures problematic [27], and the effect of thermal anisotropy may influence the results of this study. According to Krayenhoff and Voogt [27], thermal anisotropy in an urban environment depends strongly on solar elevation and irradiance, and it increases in tall, closely-spaced buildings. Further, compact and high-rise zones generate greater anisotropy than open low-rise zones, while the regularity of street orientation increases anisotropy. As follows from the description of the study areas (see Section 2) and also from the results of LCZ mapping (Section 3), typical values of several of the above parameters should suffer less from the thermal anisotropy effect. Very high and closely-spaced buildings are quite rare in Prague and Brno. The average building height in the LCZ 2 zone of Prague is 19.9 m (typically between five and seven stories) and 20.9 m in Brno (typically 5-6 stories). Compact mid-rise zones (LCZ 2) constitute 1.32% and 0.31% of the Prague and Brno areas, respectively. Finally, streets are not regularly oriented in the two cities in the fashion typical of some larger cities elsewhere, e.g., in North America.
Moreover, one may assume that the negative effect of thermal anisotropy is partly suppressed in this study because eight different thermal images acquired from May-September were used. Different dates of data acquisition have slightly different viewing geometry and, thus, partly "average" the negative effect of thermal anisotropy. In spite of this, not only the bias due to thermal anisotropy, but also other aspects of LST (e.g., atmospheric corrections, LST~air temperature relationships) should be further clarified. However, these problems are beyond the scope of the current article.
It is important to stress that the same sets of warmest (coldest) LCZs were identified by both data sources. Last, but not least, an important precondition for LST comparison between different zones is mutual independence of LCZ classification and the surface temperature field. The GIS-based method used for LCZ delimitation in this study met this requirement. This is an important methodological aspect of this study with respect to widely-used imagery-based methods [14,42,43], including thermal bands for LCZ delimitation.

Conclusions
This paper compared the spatial distribution of land surface temperatures with local climate zones for the cities of Brno and Prague in the Czech Republic, employing data provided by two satellites. The cities studied are characterized by a similar spatial distribution of local climate zones previously constructed on the basis of a GIS method. This method utilizes data completely independent of those used for LST derivation, a feature that permitted the use of LSTs for unbiased investigation of the extent to which LCZs correspond with surface temperatures.
The key findings of this study may be summarized as follows: (i) regardless of the city studied or the satellite data employed, LSTs show typical surface temperatures that differ significantly between zones (in 89.3% and 91.7% of all tests applied for Brno and Prague, respectively); (ii) the warmest (LCZ 10, heavy industry; LCZ 3, dense low-rise buildings; LCZ 2, compact mid-rise buildings) and the coldest (LCZ G, water bodies; LCZ A, areas with dense trees) zones were identified in the two study areas; (iii) LCZ 8 (large low-rise), LCZ 10 (heavy industry) and LCZ D (low plants) were well differentiated in terms of surface temperatures; LCZ 2 (compact mid-rise), LCZ 4 (open high-rise) and LCZ 9 (sparsely built-up) constitute less clearly distinguishable zones in both cities.
It may be concluded that our findings generally support the concept of LCZs and the GIS-based method used for their delimitation for two typical central European cities of different sizes. Although the concept of LCZ was developed for the classification of air temperature measurements, the results of this paper confirm that individual LCZ prove also characteristic features of surface temperatures. Such conclusions were indicated in other studies [18,19,44]. In combination with the use of the GIS-based method, this LCZ-LST comparison may be considered the most important finding of this this contribution. Certain questions, such as seasonality in LST differences and thermal anisotropy (complete surface temperature differences), remain open to future research.