Development of a Building ‐ Scale Meteorological Prediction System Including a Realistic Surface Heating

: Microscale urban meteorological models have been widely used in interpreting atmospheric flow and thermal discomfort in urban environments, but most previous studies examined the urban flow and thermal environments for an idealized urban morphology with imposing neutral or homogeneous thermal forcing. This study has developed a new building ‐ scale meteorological prediction system that extends the ability to predict microscale meteorological fields in real urban environments. A computational fluid dynamics (CFD) model has been developed based on the non ‐ hydrostatic incompressible Reynolds ‐ averaged Navier ‐ S tokes (RANS) equations with a standard k ‐ ε turbulence model, and the microscale urban surface energy (MUSE) model was coupled with the CFD model to provide realistic surface thermal boundary conditions in real urban environments. It is driven by the large scale wind and temperature fields predicted by the Korean operational weather prediction model. The validation results of the new building ‐ scale meteorological prediction system were presented against wind tunnel data and field measurements, showing its ability to predict in ‐ canyon flows and thermal environments in association with spatiotemporal variations of surface temperatures in real urban environments. The effects of realistic surface heating on pedestrian level wind and thermal environments have been investigated through sensitivity simulations of different surface heating conditions in the highly built ‐ up urban area. The results implied that the inclusion of surface thermal forcing is important in interpreting urban flow and thermal environment of the urban area, highlighting a realistic urban surface heating that should be considered in predicting building ‐ scale meteorology over real urban environments. this study is capable of predicting the atmospheric wind and temperature structures under neutral and unstable stability conditions.


Introduction
Rapid urbanization, which is projected that two-thirds of the world's population would live in cities by 2050 [1], causes various environmental problems in association with unique micrometeorological states in urban areas (e.g., [2,3]). Urban surfaces are mostly characterized by buildings, roads, and permeable natural surfaces, and various physical processes such as radiative trapping, turbulence heat exchange, sub-surface thermal conduction occur, which results in unique urban micrometeorology. Thus, understanding of the local micrometeorology is important to minimize various environmental problems that have occurred in urban areas. Many wind tunnel experiments and field measurements have been conducted to understand the urban meteorological environments (e.g., [4][5][6][7][8]). Wind tunnel experiments have advantages in collecting high-density spatial data and in interpreting microscale flow and thermal structures because wind tunnel data are obtained from well-designed experimental setups in terms of ideal/real cities' morphological structures and meteorological driving forcing (e.g., [4,6]). However, they have limitations in that the multiscale meteorological influences in the urban boundary layer (UBL) are not considered. Field measurements can collect the meteorological data that reflect the UBL influences in the real urban area of interest, which compensates the wind tunnel experiments to some extent. However, the spatial density of the measurement data is generally not enough to resolve the microscale meteorological states of an urban area (e.g., [7,8]).
Microscale meteorological models have been widely used in investigating atmospheric flow, dispersion, and thermal comfort conditions in urban environments, compensating for the limitations in the wind tunnel experiments and field measurements (e.g., [9][10][11][12][13][14][15][16][17]). Kim and Baik [9] found that urban morphology is an important factor in determining in-canyon flow and turbulence in an ideal two-dimensional canyon morphology using a computational fluid dynamics (CFD) model. Park et al. [12] examined the effects of homogeneous wall surface heating on in-canyon flow applying a largeeddy simulation (LES) model for a single canyon circulation simulation, which showed the wall surface heating can alter the in-canyon vortex circulation in shape and strength. Recently, microscale meteorological applications to real urban environments are increasing by virtue of the increase of available computer resources and the enhancement of the geographic information system (GIS) data (e.g., [18][19][20]). Gronemeier et al. [18] and Wang et al. [20] examined air ventilation under unstable atmospheric boundary layer in the highly built-up urban area of Hong Kong by imposing a homogeneous surface heating within the parallelized large-eddy simulation model (PALM). The studies reported that homogeneous surface heating led to an increase in pedestrian-level velocity ratio and wind directional change. Meanwhile, Nazarian and Kleissl [15] investigated the effects of in-canyon heterogeneous surface heating on atmospheric flow over an idealized cube array configuration, for which an urban surface energy balance model [21] was used to predict spatial and temporal variations of the surface temperatures. The results showed that heterogeneous surface heating modified the mean atmospheric flows and turbulent transfer efficiency particularly under strong unstable atmospheric conditions.
The previous studies have examined characteristic features of urban flow and thermal environments for an idealized urban morphology with imposing heterogeneous surface heating or for a real urban area with imposing neutral or homogeneous thermal forcing. This study aims to develop a new building-scale meteorological prediction system in order to extend the ability to predict microscale meteorological fields with realistic surface heating in real urban environments. In doing so, heterogeneous surface temperatures of the urban facets of a highly built-up real urban area have been considered using the microscale urban surface energy (MUSE) model. It also investigates the effects of heterogeneous thermal forcing on pedestrian-level wind and thermal environments in the highly built-up urban area. The remainder of the manuscript is structured as follows: Section 2 describes the newly developed building-scale meteorological prediction system and Section 3 presents validation results of the building-scale meteorological prediction system against wind tunnel data and field measurements. Section 4 discusses the effects of realistic surface heating on pedestrian-level wind and thermal environments in the highly built-up urban area, and summary and conclusions are followed in Section 5.

The Governing Equations
The CFD model has been developed based on the dynamic equations and numerical methods in Baik et al. [10]. The governing equations consist of the Reynolds-averaged conservation equations of mass, momentum, and heat under the Boussinesq approximation as follows [10,22,23]: .
Here, is the Reynolds-averaged wind velocity component in the th space coordinate of , is the air density, * is the pressure perturbation from the reference state, is the Reynoldsaveraged air temperature; and are the perturbation fields of th wind velocity component and temperature, respectively; and are are the air density and temperature of the reference state, respectively; and are the kinematic viscosity and molecular thermal diffusivity of the air, respectively; is the Kronecker delta and is the gravitational acceleration; is the diabatic source term; is the time. and represent the Reynolds stresses and turbulent heat fluxes, respectively, which are parameterized using the K-theory as follows [24]: where and are the eddy diffusivity of momentum and heat, respectively, and is the turbulent kinetic energy. The eddy diffusivity of momentum is expressed as a function of the turbulent kinetic energy ( ) and its dissipation rate ( ) as follows: , (6) where is the empirical constant (=0.09). The eddy diffusivity of heat can be expressed by the turbulent Prandtl number ( ) (=0.7) as follows: .

Numeric Representation
The governing equations are fully implicitly discretized on the Arakawa-C staggered grid system using a finite volume method and a power-law interpolation scheme [10], which are solved numerically with the semi-implicit method for a pressure-linked equation (SIMPLE) algorithm [25,26]. The power-law scheme, which is based on the analytical solution of an advection-diffusion equation and is known as an accurate and efficient scheme for CFD [25,26], was used to interpolate physical quantities at the grid cell faces. The SIMPLE algorithm iteratively calculates converged divergent-free solution of the wind velocity components ( ) and pressure perturbation ( * ) at every time step. It has a merit that a stable numerical solution can be obtained due to the fully implicit discretization. Each grid is shaped with a uniform rectangular grid on the Cartesian coordinate. The scalar variables of temperature, pressure, turbulent kinetic energy and dissipation rate are located at the center of the grid cell and the wind velocity components are at the center of the grid cell faces. In the staggered grid system, it is advantageous that the wind velocity components are forced by the pressure difference between the two adjacent cell faces defined at each grid without numerical approximation. It is also an important advantage that the three-dimensional scalar transport across the faces of grid cells can be exactly conserved and calculated without interpolation of the wind velocity components. In the CFD model, urban form structures with varying heights of buildings, trees, and topography are represented by 3-dimensional rectangular grid cells stacked in multiple vertical layers. It can easily represent complex morphological structures in real urban environments. The adjacent atmospheric grid cells surrounding the surface obstacle grid form the bottom boundaries of the domain.
The free-slip conditions are applied at the upper boundary, while the wall function is applied for the bottom boundary conditions based on Launder and Spalding [27] and Versteeg and Malalasekera [26].
The lateral boundary conditions can be selected from a constant inflow, zero-gradient, and periodic conditions depending on an experimental setup.

The Microscale Urban Surface Energy (MUSE) Model
Various urban physical processes occur in association with complex urban surface conditions (e.g., [28]). The MUSE model has been developed to support microscale meteorological models of CFD and LES in simulating atmospheric flow and thermal environments in real urban environments [29], which was designed to represent urban surfaces compatible with 3-dimensional rectangular grid cells the same as in microscale meteorological models so that it can be easily coupled with the microscale meteorological models. It includes explicit parameterizations of shortwave and longwave radiative transfer processes, turbulence momentum and heat exchanges, and sub-surface thermal conduction based on surface energy balance. An embedded shadow model predicts time-varying shadows by considering the detailed distribution of buildings in real urban environments, and the 3dimensional view factors, which are calculated from the analytically-based numerical method [17] that is used to take the shortwave and longwave radiation trapping effects within urban canyons into account. The MUSE model calculates surface temperatures and surface energy balance fluxes at each grid cell using the meteorological forcing of wind velocity components, air temperature, specific humidity, pressure, air density, radiation, and precipitation. In this study, the performance of the MUSE model was presented and validated against the field measurement obtained from a highly built-up urban area in South Korea. A detailed description of the MUSE model will be given in Lee and Lee [29].

The Building-Scale Meteorological Prediction System
The new building-scale meteorological prediction system has been developed with coupling the CFD model and the MUSE model to predict microscale atmospheric flow and thermal environments for real urban areas ( Figure 1). The large scale wind and temperature profiles predicted by the Local Data Assimilation and Prediction System (LDAPS), which is a Korean operational local weather forecast system, are used as a meteorological driving forcing of the CFD and MUSE models. The LDAPS routinely produces +48 h hourly local weather forecasts including the wind velocity components and temperature at a horizontal resolution of 1.5 km and 70 vertical layers every 6 h intervals over South Korea. The meteorological profiles of the nearest LDAPS grid to the CFD domain are interpolated to the vertical grid of the CFD model for use in lateral boundary conditions. Overall, the three models were coupled in a one-way interacting approach. The MUSE model calculates the surface temperatures at each grid cell using the predicted meteorological forcing variables by the LDAPS model, which is provided as the bottom thermal boundary condition of the CFD model. The CFD model predicts microscale atmospheric flow and thermal environments using the realistic surface heating in real urban environments by the MUSE model and the meteorological wind and temperature profiles by the LDAPS. In an offline mode, the meteorological forcing of the MUSE model can be obtained from the CFD model and/or field measurements.
The profiles of turbulent kinetic energy ( ) and its dissipation rate ( ), which are also required in the CFD model as lateral boundary conditions, were calculated as follows: , where * is the friction velocity, is the wind speed at the top of the CFD domain, z is the roughness length, δ is the boundary layer height, is the von Karman constant (=0.4), and C is the empirical constant (=0.09) [10].

Validation of the CFD Model against Wind Tunnel Data
The CFD model of the new building-scale meteorological prediction system has been validated against the wind tunnel data of Uehara et al. [6] to investigate the performance of the CFD model in predicting atmospheric profiles of urban flow and air temperature. Uehara et al. [6] examined the effects of atmospheric stability on atmospheric flow in urban street canyons through the wind tunnel experiments with varying atmospheric stability, which measured atmospheric flow and temperature fields within and above the street canyon. The urban geometry was configured using building blocks with a size of 0.1 × 0.1 × 0.1 m 3 , which were arranged homogeneously with a spacing of 0.1 m in the x-direction and 0.05 m in the y-direction, respectively. Thus, the urban geometry parameters of H/W (building height/canyon width) and H/L (building height/canyon length) are 1.0 and 0.5, respectively. The urban geometry corresponds to the skimming flow regime [30]. The inflow turbulent boundary layer formed 7H in depth with a roughness length of 0.0033 m and a zero-plane displacement height of 0.035 m before impinging the array of the building blocks. The thermal forcing in the wind tunnel experiments was assigned with the air temperature of 20 °C and the surface temperature of 79 °C, which results in an unstable atmospheric condition with the bulk Richardson number of −0.21. The measurements were conducted between the fifth and sixth rows of the building blocks and the vertical sections of wind speed and temperature were taken at the center of the street canyon.
The domain of the CFD model was set identical to the wind tunnel experiment (Figure 2). The model domain has 13.3 m (x) × 0.8 m (y) × 0.6 m (z). The building blocks were resolved with a grid spacing of 0.01 m. The geometric and dynamic similarities were taken into consideration. The inflow flow condition was assigned using a power-law profile of . , where the reference wind speed ( ) was set to 3 m s −1 . With the constant wind profile, the turbulent boundary layer was developed for a patch distance of 12 m with the same surface roughness conditions of the wind tunnel experiment, and then the well-developed turbulent flow impinged the rows of building blocks. A zero-gradient condition was applied for the outflow and lateral boundary conditions. The atmospheric temperature and the ground surface temperature were set identical to the wind tunnel experiment, configuring the atmospheric stability of Ri = −0.21.  Figure 3 compares the measured and simulated mean wind and temperature profiles at the center of the street canyon. The mean wind profiles were normalized by the streamwise wind velocity at 2H and the temperature profiles were expressed by the temperature difference from the temperature at 2H normalized by the temperature difference between the ground (Tg) and the temperature at 2H (T2H). It is observed that the approaching mechanical forcing formed a single vortex structure within the canyon and the sharp gradient just above the building block in the streamwise direction. The center height of the vortex is located at approximately 0.6H. In the neutral simulation, the observed primary vortex structure and the sharp vertical gradient were reasonably captured because the observed characteristics are primarily mechanically driven. The simulation with the unstable atmospheric stability compares against the wind tunnel data better than the neutral simulation, enhancing the in-canyon vortex intensity and the vertical gradient ( Figure 3a). The improvement in the streamwise wind velocity is more clearly identified near the ground, which might be influential in predicting the pedestrian-level wind environments. It is also shown that the surface heating formed the characteristic temperature profiles that have a strong gradient near the ground (z/H < 0.1) and the well-mixed structure within the canyon (0.1 < z/H < 1.0). The temperatures gradually decrease above the building block (z/H > 1.0). The observed temperature profiles within and above the street canyon were well simulated by the CFD model (Figure 3b). This result indicates that the CFD model developed in this study is capable of predicting the atmospheric wind and temperature structures under neutral and unstable stability conditions. . Figure 3. Comparison of the simulated and measured (a) wind and (b) temperature profiles at the center of a street canyon.

Validation of the MUSE Model against Field Measurements in a High-Rise Commercial Area
The MUSE model of the new building-scale meteorological prediction system has been validated against field measurements collected in a highly built-up real urban area. Many aspects associated with heterogeneous surface temperatures in real urban environments were considered for validation and meteorological simulation of the building-scale meteorological prediction system. A high-rise built-up area in Seoul, South Korea (37.572° E, 126.978° N) was selected for this study. The area can be characterized by a high dense commercial/residential area with high and low-rise buildings (Figure 4a). Various meteorological variables were continuously measured by the 3 dimensional sonic anemometer, thermometer, and infrared thermal imagery camera deployed at a 7-m meteorological tower on a building rooftop (71 m above the ground level) [8]. The simulation domain of the MUSE model was configured by 120 × 120 × 27 grid mesh with a horizontal and vertical spacing of 5 m and 3 m, respectively, covering the area of 600 (x) × 600 (y)× 81 (z) m 3 surrounding the measurement tower ( Figure 4b). The surface grid representation was used identically with the CFD model except that the vertical grid was stretched up to 300 m in meteorological simulations. The thermal properties of road and building surfaces were assigned with asphalt and concrete, respectively [31]. The simulation was conducted for a clear day of 7 September 2015 when weak easterly winds were dominant throughout the day. Figure 5 shows the diurnal variations of the radiative and meteorological forcing variables on the day. The downward shortwave and longwave radiation were obtained from the automatic surface observation station operated by the Korea Meteorological Administration (KMA), which is located about 1 km west. The downward longwave radiation was calculated from the observed atmospheric temperature and cloudiness at the site [31,32]. The atmospheric temperature, pressure, and wind were obtained from the tower measurements ( Figure  5b). The vertical profiles of atmospheric temperature were assigned homogeneously, while the vertical wind profiles were modified using the wind fields of the neutral CFD simulation.    Figure 7 shows the spatial distributions of the simulated surface temperatures at the wall and ground surfaces and the frequency distribution over the urban area during the day. The model simulated the heterogeneous distribution of the surface temperatures at different times of the day through the consideration of the surface energy balance of each urban facet. At 09 LST, the spatial variation of the wall ground surface temperatures was relatively small, approximately 20 °C, within the model domain ( Figure 7a). As the downward shortwave radiation increased, the surface temperature of each urban facet gradually increased up to approximately 50 °C (Figure 7b-d). In addition, the spatial heterogeneity in the surface temperatures also significantly increased over 20 °C during the daytime (Figure 7e). The large spatial heterogeneity of the surface temperatures was primarily due to the difference in the incident direct shortwave radiation on each urban facet. This result showed that the realistic representation of the radiative energy exchanges was primarily important in determining the urban surface temperatures and forming thermal heterogeneity in the real urban environment.   (P2 and P4), and a north-facing wall (P3) were obtained from the infrared thermal imagery measurement at the meteorological tower site (Figure 8a). The roof surface temperatures, which were obtained from the infrared thermometer at the rooftop of the meteorological measurement tower site, varied 18-42 °C during the day with a very large diurnal range of 24 °C. The nocturnal surface temperatures were well simulated by the MUSE model, but the daytime surface temperatures were slightly underestimated (Figure 8b). The model well simulated the observed characteristic diurnal variations with a reasonable diurnal range at east-facing and north-facing wall surfaces, but it showed systematic cold biases at all the wall surfaces throughout the day (Figure 8c). The diurnal variation of the ground surface temperatures was well compared with the roof surface temperatures showing a slightly reduced diurnal range of 21 °C. The model simulated the observed diurnal variation of the road surface temperatures, but it slightly overestimated the observed diurnal range (Figure 8d). The model-measurement discrepancies might attribute to various factors such as thermal properties of the urban surfaces, uncertainties in the meteorological forcing and measurement, but it is difficult to separate the errors due to complex interactions of the urban physical processes in the real simulation. Despite the discrepancies, overall the validation results showed that the model represented well the spatial and temporal variations of the urban surface temperatures.

Validation of the Building-Scale Meteorological Prediction System against Field Measurement
The building-scale meteorological prediction system has been applied over the real urban area for the day of 7 September 2015, in which the surface temperatures predicted by the MUSE model were incorporated as a thermal bottom boundary conditions and the wind and temperature profiles operationally predicted by the LDAPS were used as a large scale meteorological forcing. The inflow profiles of turbulent kinetic energy and its dissipation rate were calculated with the roughness length and the boundary layer height of 0.1 m and 1000 m, respectively. The building-scale meteorological prediction system produced hourly microscale wind and temperature fields that are influenced by the realistic surface heating and the large scale meteorological forcing in the highly built-up urban area. Figure 9 compares the atmospheric temperature and wind fields simulated by the buildingscale meteorological prediction system against the field measurements at the meteorological tower site. The operational forecast by the LDAPS was also compared to evaluate the potential performance of the new microscale model in simulating atmospheric temperature and wind fields over the urban area. The atmospheric temperatures measured at 78 m a.g.l. ranged from 16 °C at 06 LST to 26 °C at 17 LST, which were simulated better by the building-scale meteorological prediction system than the LDAPS especially during the daytime (Figure 9a). The improvement was largely attributed to the realistic representation of the urban surface temperatures shown in Figure 7. The low predicted atmospheric temperature at night might be largely attributed to the discrepancy in the LDAPS prediction. The measured wind speeds and directions ranged 1-3 m s −1 and 60-90° (easterly) throughout the day, respectively (Figure 9b,c). The diurnal change in the measured wind direction was reasonably well simulated by both the models (Figure 9c). The LDAPS overestimated the measured wind speeds throughout the day, whereas the building-scale meteorological prediction system showed a better agreement with the measurement (Figure 9b). The better agreement of the wind speed in the building-scale meteorological prediction system might be attributed to the explicit representation of surface momentum drag of the real urban geometries, in which the pressure drag by buildings was predominant compared to the Reynolds stresses. Overall, the validation results showed that the new building-scale meteorological prediction model had a better capability than the LDAPS in simulating microscale meteorological fields in the highly built-up urban areas.

The Effects of Realistic Surface Heating on Pedestrian-Level Wind and Temperature Fields
Two sensitivity simulations with no surface heating and homogeneous surface heating have been further conducted for the built-up urban area using the building-scale meteorological prediction system to investigate the effects of the realistic surface heating on pedestrian level wind and thermal environments. The former simulation considered no surface heating at all urban facets within the domain and the latter simulation assigned the domain averaged surface temperature identically to all the urban facets at every hour (Figure 7e). All three simulations were conducted for 7 September 2015 with the same configuration except for the bottom thermal boundary conditions. The wind velocity ratio (V ), which is commonly used in the ventilation studies (e.g., [33]), was used in the analysis of pedestrian-level wind environment, which is defined by Here, is the wind speed at the pedestrian level and is the wind velocity just above the atmospheric boundary layer. In this study, and were used by the wind velocity at the first vertical grid level of 1.5 m and at 300 m above the ground level, respectively. The analyses of the realistic surface heating effects were investigated for the daytime period with enhanced thermal heterogeneity in the urban area. Figure 10 compares the spatial distributions of the pedestrian-level wind field and wind velocity ratio in the highly built-up urban area simulated with the different surface heating conditions at 10 LST and 14 LST. At 10 LST, the large scale inflow wind direction and speed were approximately 60° and 4 m s −1 (Figure 9), respectively, and the surface temperatures in the urban facets of the domain ranged 14.9-29.7 °C and 22.6 °C on average (Figure 7e). The simulated pedestrian-level wind fields showed very complicated variations depending on the thermal forcing as well as the buildings' shape and their spatial arrangement (Figure 10a-c). It is found that that the spatial distributions in wind direction show better consistency among the three simulations than in wind speed, which is mainly attributed to the formation of channeling flows between the buildings and wake flows behind buildings. In addition, the thermal forcing also significantly modified the wind fields, especially in their directions. Wang and Ng [33] reported the enhanced variance of wind direction in unstable atmospheric conditions in their PALM simulation with homogeneous surface heat fluxes in Hong Kong. These characteristic urban flows were also reasonably simulated by the PALM model in the high-rise building district in Hong Kong [33,34]. At 14 LST, the large scale inflow slightly changed counter-clockwise with the wind direction and speed of approximately 30° (northeasterly) and 3 m s −1 (Figure 9), respectively, and the surface temperatures and their spatial heterogeneity were significantly increased ranging 19.6-42.5 °C and 29.7 °C on average (Figure 7e). The slight change in large scale inflow direction and speed, compared to 10 LST, led to large differences both in the pedestrian-level wind direction and speed as a consequence of the changes in the characteristic urban flows (Figure 10d-f). The surface thermal forcing also led to large changes from the neutral simulation in the pedestrian-level wind direction and speed, which is more clearly identified in 14 LST than 10 LST following the increased thermal forcing. Figure 11 shows the differences in the simulated pedestrian-level wind velocity ratio in the highly built-up urban area among the sensitivity simulations. The heterogeneous surface heating changed the spatial distribution of the wind velocity ratio within the built-up area ranging by −0.17-0.07 in 10 LST and by −0.31-0.17 in 14 LST. The air ventilation within the building blocks was overall enhanced on the domain through the inclusion of the realistic surface heating, but the changes are significantly dependent on specific locations of the area and the intensity of the surface thermal forcing (Figure 11a,c). When compared to the homogeneous heating simulation, the spatial heterogeneity in urban surface temperatures changed the wind velocity ratio within the built-up area ranging by −0.04-0.04 in 10 LST and by −0.16-0.05 in 14 LST. The relatively large difference in the wind velocity ratio at 14 LST is attributed to the strong spatial heterogeneity in the urban surface temperatures. Overall, these results indicate that it is important to resolve the realistic prediction of the urban surface temperatures in predicting wind fields and assessing ventilation at a specific location in real urban areas.   Figure 12 compares the simulated wind velocity ratio values in the highly built-up urban area from 10-16 LST. The wind velocity ratio was binned by 0.05 and the occurrence frequency of each bin over the urban area was compared among the three simulations. The surface heating decreased the simulated wind velocity ratio values in the lower V classes less than 0.5, while it increased the values in the higher V classes (Figure 12a). This is a similar result found in Wang and Ng [33] which showed the increased wind velocity ratio values in the higher V classes greater than 0.2 over a dense urban district in Hong Kong. This study shows that the simulation with heterogeneous surface temperatures gives slightly lower thermal influences on the wind velocity ratio than the homogeneous surface heating, which is also shown in the statistical distributions (Figure 12b). The wind velocity ratio in the building district ranged 0.03-0.76 in the neutral simulation, 0.03-0.90 in the homogeneous simulation, and 0.03-0.85 in the heterogeneous simulation with a median wind velocity ratio of 0.31 in the neutral simulation, 0.33 in the homogeneous simulation, and 0.32 in the heterogeneous simulation.  Figure 13 compares the spatial distributions of the pedestrian-level air temperature in the highly built-up urban area simulated with the homogeneous and heterogeneous surface heating conditions. The simulated pedestrian-level air temperatures ranged 19.6-21.8 °C at 10 LST and 24.0-26.1 °C at 14 LST, showing spatial variations of approximately 2 °C within the urban district (Figure 13b,e). The spatial variations in pedestrian-level air temperature are not resolved in the operational mesoscale weather forecast model of the LDAPS. Meanwhile, the simulated pedestrian-level air temperatures with heterogeneous surface heating were similar in statistical distribution within the urban district (Figure 13c,f). However, the spatial distribution of the simulated temperatures was quite different between the homogeneous and heterogeneous surface heating simulations, ranging −0.41-0.30 °C at 10 LST and −0.81-0.46 °C at 14 LST, respectively ( Figure 14). Large differences in the pedestrian-level air temperature were found near the wall surfaces of the buildings, which indicates that accurate representation of the urban surface temperatures is important.

The Thermal Effects on Pedestrian-Level Thermal Environment
Overall, these results inferred that reasonable representation of both the urban mechanical and thermal forcing is a prerequisite in accurately predicting the pedestrian-level wind and thermal environments of the real urban area. The new building-scale meteorological prediction system well captured the characteristic urban flows (e.g., channeling flows between buildings and wake flows behind buildings) and the spatial variations of the pedestrian-level wind and air temperature fields in the highly built-up real urban area by virtue of the inclusion of the realistic surface heating by the MUSE, the explicit representation of mechanical forcing of buildings, and the large-scale meteorological forcing of the LDAPS into the CFD model.

Summary and Conclusions
Microscale urban meteorological models have been widely used in interpreting urban flow and thermal discomfort mostly under idealized urban morphological and physical environments. In this study, a new building-scale meteorological prediction system has been developed to extend the ability to predict microscale meteorological fields in real urban environments. A CFD model has been developed following Baik et al. [10], which solves the non-hydrostatic incompressible RANS equations with a standard k-ε turbulence model. The prognostic variables were discretized on the Arakawa-C grid system by a finite volume method in which the interpolation process was calculated by the power-law scheme. The new building-scale meteorological prediction system has been developed through a one-way coupling of the CFD model and the MUSE model, which is driven by the large scale wind and temperature fields predicted by the LDAPS/KMA. The MUSE model includes the urban physical parameterizations of building shadow effects, in-canyon direct and diffuse shortwave radiative transfer processes, turbulent sensible heat exchange processes between artificial surfaces and their surrounding atmosphere, and sub-surface thermal conduction processes, consequently being capable of predicting urban surface temperatures for real urban environments. The MUSE model provides realistic surface heating conditions for the CFD simulations in real urban environments.
The new building-scale meteorological prediction system was validated in its ability of canyon flow and thermal structure prediction against the wind tunnel data and field measurements obtained in a highly built-up urban area in Seoul, Korea. The results showed that the CFD model well captured the characteristic features of the measured in-canyon vertical temperature structure, showing the importance of surface heating. Additionally, the surface temperatures of different urban facets simulated by the MUSE model compared well against the IR imagery surface temperatures in the highly built-up urban area, highlighting the importance of the solar position and building morphology in urban surface temperature prediction. Overall, the validations showed its ability to predict in-canyon flows and thermal environments in association with spatial and temporal variations of urban surface temperatures in the real urban environment.
The effects of realistic surface heating on pedestrian level wind and thermal environments have been further investigated through three sensitivity simulations of no heating, homogeneous heating, and heterogeneous heating conditions in the highly built-up urban area for a clear day of 7 September 2015. The results showed that either of homogeneous or heterogeneous surface heating increased domain-mean pedestrian level wind speeds and their spatial heterogeneity compared to the neutral simulation, resulting in the higher frequency of high V classes and the lower frequency of low V classes. The change in mean air temperature was minor because the thermal forcing of the homogeneous heating case was assigned identically with a domain mean surface temperature. Meanwhile, a similar statistical frequency distribution was found both in the pedestrian level velocity ratio and the air temperature at the homogeneous and heterogeneous surface heating cases, even though the wind speed at each grid location has changed in time and space according to the thermal boundary conditions. These results imply that the inclusion of surface heating, as a forcing of microscale meteorological models, is important in interpreting urban flow and thermal comfort of a certain urban area in a statistical sense, for which homogeneous heating becomes a good assumption. However, it is identified that realistic urban surface heating should be considered in predicting building-scale meteorology over real urban environments. Most previous studies examined the urban flow, dispersion, and thermal discomfort of an urban area imposing neutral or homogeneous thermal forcing for an idealized urban morphology, this study highlights the importance of imposing realistic urban surface heating in microscale urban meteorological simulations. In this sense, the new building-scale meteorological prediction system developed in this study can be a useful tool for the evaluation and prediction of urban flow and thermal conditions for real urban environments.
More research is needed to evaluate and improve the building-scale meteorological prediction system. Intensive field measurements are necessary to evaluate the model's capabilities faithfully in the real urban environment. In addition, a two-way coupling of the CFD model and the MUSE model will be able to provide a more realistic modeling strategy.