On the Derivation of Input Parameters for Urban Canopy Models from Urban Morphological Datasets

Two techniques are analysed to derive mean street width and mean building width from morphological data of real cities: one based on a two-dimensional simplified morphology, and the other based on a three-dimensional regular simplified morphology. For each simplified morphology (two-dimensional and three-dimensional), the sky-view factors (street-to-sky) are computed and compared with the sky-view factors derived from the real morphology for selected districts of three European and two North American cities. The two-dimensional simplified morphology reproduces the real sky-view factors better than the three-dimensional morphology. Since many urban canopy parameterizations represent the city using simplified morphologies, this can be useful information for the derivation of input parameters for urban canopy parameterizations from real morphological data.     

The Influence Of Urban Canopy Configuration On Urban Albedo

We propose a calculation method for shortwave radiation flux and longwave radiation flux within the urban canopy and investigate the influence of urban canopy configuration on net radiation flux. In the assumed urban configuration, buildings of equal size are arranged in a regular lattice within the urban canopy. The net shortwave radiation flux and longwave radiation flux within the urban canopy were calculated by the photon tracking method based on the Monte Carlo method. The albedo value obtained by this method shows close agreement with experimental data, and the average sky view factor shows almost perfect agreement with the theoretical value. Moreover, we calculated the urban albedo for the urban canopy configuration including roads and building height distribution.%Moreover, we calculated net radiation within the urban canopy in %consideration of roads and building height distribution.We found that the sky view factor of the ground surface is high when building coverage is low, building height is low, open space by roads exists, and building height is non-uniform. Moreover, we found that the albedo value is high when building height is small, open space by roads is wide, and building height is uniform. The albedo value was found to vary in a complicated manner with change in building coverage.     

A Scale-Adaptive Approach for Spatially-Varying Urban Morphology Characterization in Boundary Layer Parametrization Using Multi-Resolution Analysis

Urban morphology characterization is crucial for the parametrization of boundary-layer development over urban areas. One complexity in such a characterization is the three-dimensional variation of the urban canopies and textures, which are customarily reduced to and represented by one-dimensional varying parametrization such as the aerodynamic roughness length $z_{0}$ and zero-plane displacement $d$ . The scope of the paper is to provide novel means for a scale-adaptive spatially-varying parametrization of the boundary layer by addressing this 3-D variation. Specifically, the 3-D variation of urban geometries often poses questions in the multi-scale modelling of air pollution dispersion and other climate or weather-related modelling applications that have not been addressed yet, such as: (a) how we represent urban attributes (parameters) appropriately for the multi-scale nature and multi-resolution basis of weather numerical models, (b) how we quantify the uniqueness of an urban database in the context of modelling urban effects in large-scale weather numerical models, and (c) how we derive the impact and influence of a particular building in pre-specified sub-domain areas of the urban database. We illustrate how multi-resolution analysis (MRA) addresses and answers the afore-mentioned questions by taking as an example the Central Business District of Oklahoma City. The selection of MRA is motivated by its capacity for multi-scale sampling; in the MRA the “urban” signal depicting a city is decomposed into an approximation, a representation at a higher scale, and a detail, the part removed at lower scales to yield the approximation. Different levels of approximations were deduced for the building height $\bar{{H}}$ and planar packing density $\lambda _\mathrm{p}$ . A spatially-varying characterization with a scale-adaptive capacity is obtained for the boundary-layer parameters (aerodynamic roughness length $z_{0}$ and zero-plane displacement $d$ ) using the MRA-deduced results for the building height and the planar packing density with a morphometric model; an attribute that is shown to be of great advantage to multi-scale and multi-resolution numerical weather prediction models.     

Geometric Dependence of the Scalar Transfer Efficiency over Rough Surfaces

We performed a series of wind-tunnel experiments under neutral conditions in order to create a comprehensive database of scalar transfer coefficients for street surfaces using regular block arrays representing an urban environment. The objective is to clarify the geometric dependence of scalar transfer phenomena on rough surfaces. In addition, the datasets we have obtained are necessary to improve the modelling of scalar transfer used for computational simulations of urban environments; further, we can validate the results obtained by numerical simulations. We estimated the scalar transfer coefficients using the salinity method. The various configurations of the block arrays were designed to be similar to those used in a previous experiment to determine the total drag force acting on arrays. Our results are summarized as follows: first, the results for cubical arrays showed that the transfer coefficients for staggered and square layouts varied with the roughness packing density. The results for the staggered layout showed the possibility that the mixing effect of air can be enhanced for the mid-range values of the packing density. Secondly, the transfer coefficients for arrays with blocks of non-uniform heights were smaller than those for arrays with blocks of uniform height under conditions of low packing density; however, as the packing density increased, the opposite tendency was observed. Thirdly, the randomness of rotation angles of the blocks in the array led to increasing values of the transfer coefficients under sparse packing density conditions when compared with those for cubical arrays.     

A New Aerodynamic Parametrization for Real Urban Surfaces

This study conducted large-eddy simulations (LES) of fully developed turbulent flow within and above explicitly resolved buildings in Tokyo and Nagoya, Japan. The more than 100 LES results, each covering a 1,000 $\times $ 1,000 m $^{2}$ area with 2-m resolution, provide a database of the horizontally-averaged turbulent statistics and surface drag corresponding to various urban morphologies. The vertical profiles of horizontally-averaged wind velocity mostly follow a logarithmic law even for districts with high-rise buildings, allowing estimates of aerodynamic parameters such as displacement height and roughness length using the von Karman constant $=$ 0.4. As an alternative derivation of the aerodynamic parameters, a regression of roughness length and variable Karman constant was also attempted, using a displacement height physically determined as the central height of drag action. Although both the regression methods worked, the former gives larger (smaller) values of displacement height (roughness length) by 20–25 % than the latter. The LES database clearly illustrates the essential difference in bulk flow properties between real urban surfaces and simplified arrays. The vertical profiles of horizontally-averaged momentum flux were influenced by the maximum building height and the standard deviation of building height, as well as conventional geometric parameters such as the average building height, frontal area index, and plane area index. On the basis of these investigations, a new aerodynamic parametrization of roughness length and displacement height in terms of the five geometric parameters described above was empirically proposed. The new parametrizations work well for both real urban morphologies and simplified model geometries.     

Surface Albedo in Cities: Case Study in Sapporo and Tokyo,Japan

The surface albedo of two large cities in Japan was measured using a pyranometer mounted on a helicopter to avoid the bidirectional reflectance distribution. The daytime albedo was 0.12 in the cities, which was less than that of a nearby forest (0.16). The albedo was dependent on building structure in the cities; the albedo was lower in areas with more buildings, and decreased as the aspect ratio of street canyons increased. There are two reasons for this dependency: the multiple reflection of radiation in the building canopy, as has been shown in many previous studies, and the sparse vegetation in urban areas. These two factors concurrently determine the albedo in a real city, where the vegetation amount decreases as the plan roof ratio increases.     

Study on Effects of Building Morphology on Urban Boundary Layer Using an Urban Canopy Model

An urban canopy model is incorporated into the Nanjing University Regional Boundary Layer Model. Temperature simulated by the urban canopy model is in better agreement with the observation, especially in the night time, than that simulated by the traditional slab model. The coupled model is used to study the effects of building morphology on urban boundary layer and meteorological environment by changing urban area, building height, and building density.It is found that when the urban area is expanded, the urban boundary layer heat flux, thermal turbulence, and the turbulent momentum flux and kinetic energy all increase or enhance, causing the surface air temperature to rise up. The stability of urban atmospheric stratification is affected to different extent at different times of the day.When the building height goes up, the aerodynamic roughness height, zero plane displacement height of urban area, and ratio of building height to street width all increase. Therefore, the increase in building height results in the decrease of the surface heat flux, urban surface temperature, mean wind speed, and turbulent kinetic energy in daytime. While at night, as more heat storage is released by higher buildings, thermal turbulence is more active and surface heat flux increases, leading to a higher urban temperature.As the building density increases, the aerodynamic roughness height of urban area decreases, and the effect of urban canopy on radiation strengthens. The increase of building density results in the decrease in urban surface heat flux, momentum flux, and air temperature, the increase in mean wind speed, and the weakening of turbulence in the daytime. While at night, the urban temperature increases due to the release of more heat storage.     

Large-Eddy Simulation of Pollutant Removal from a Three-Dimensional Street Canyon

Large-eddy simulations were conducted to investigate the mechanism of pollutant removal from a three-dimensional street canyon. Five block configurations with aspect ratios (building height to length) of 1, 2, 4, 8 and $\infty $ were used to create an urban-like array. A pollutant was released from a ground-level line source at the centre of the target canyon floor. For smaller aspect ratios, the relative contribution of the turbulent mass flux to net mass flux at the roof level, which was spatially averaged along the roof-level ventilation area, was closer to unity, indicating that turbulent motions mainly affected pollutant removal from the top of the canyon. As aspect ratio increased, the relative contribution became smaller, owing to strong upwind motions. However, the relative contribution again reached near unity for the infinite aspect ratio (i.e. a two-dimensional street canyon) because of lowered lateral flow convergence. At least 75 % of total emissions from the three-dimensional street canyon were attributable to turbulent motions. Pollutant removal by turbulent motions was related to the coherent structures of low-momentum fluid above the canyons. Though the coherent structure size of the low-momentum fluid differed, the positions of low-momentum fluid largely corresponded to instantaneous high concentrations of pollutant above the target canyon, irrespective of canyon geometry.     

Numerical Modelling of Urban Heat-Island Intensity

A three-dimensional, non-hydrostatic, high-resolution numerical model was used toanalyse urban heat-island (UHI) intensity in an idealised but realistic configuration.The urban area was 20 km square and lay on flat land at about latitude 50° Nin a maritime climate. In the model the urban area was represented by anomalies ofalbedo, anthropogenic heat flux, emissivity, roughness length, sky-view factor (SVF),surface resistance to evaporation (SRE) and thermal inertia. A control simulationincluded all these factors and the resultant UHI structure, energetics and intensitywere validated against observations. The results also compared favourably withearlier simulations.A series of experiments was conducted in which successively one of the anomaliesthat represented the urban area was omitted from the control simulation so as toprovide the basis for an assessment of its effect. In daytime the individual effectsdue to albedo, anthropogenic heat, emissivity, SVF and thermal inertia ranged from0.2 to 0.8 °C. In common with albedo, anthropogenic heat, emissivity andSVF, the SRE aided the formation of a UHI; it was also the most important factorin increasing its intensity. The roughness length had the opposite effect. At nightemissivity, roughness length, SVF and SRE had effects ranging from 0.3 to0.75 °C, but the largest effect (2 °C) was due to the anthropogenicheat. These results showed a difference in the causes of daytime and nighttime UHIs.In daytime the roughness length and SRE were the most important factors affectingUHI intensity; at night the anthropogenic heat was the most important. The simulationssuggested that the size of the urban area had a minimal effect on UHI intensity.     

A Method for Monitoring the Heat Flux from an Urban District with a Single Infrared Remote Sensor

The proposed methodology relies on the modelling capabilities of the thermo-radiative model Solene to simulate the heat and radiation energy exchanges between an actual urban district and the atmosphere. It is based on the comparison of the simulated upward infrared and sensible heat flux diurnal cycles that may be measured by elevated sensors above the three-dimensional scene, as a function of sensor position: the heat flux is a function of an equivalent surface temperature given by the infrared sensor and an equivalent heat transfer coefficient deduced from Solene simulations with the actual geometry. The method is tested against measurements obtained in the city centre of Toulouse, France during an experimental campaign in 2004–2005. To improve the computation of the heat exchanges between air and building surfaces a new algorithm is first implemented, based on an empirical model of the wind distribution within street canyons. This improvement is assessed by a direct comparison of the simulated brightness surface temperatures of the Toulouse city centre to measurements obtained with an airborne infrared sensor. The optimization of the infrared remote sensor position is finally analyzed as a function of its height above the mean roof level: it allows evaluation of the heat flux from an urban district when the three different classes of surfaces (roofs, walls, grounds) have similar contributions to the infrared flux towards the sensor, and to the heat flux into the atmosphere.     

Estimating Aerodynamic Parameters of Urban-Like Surfaces with Heterogeneous Building Heights

There are many geometrical factors than can influence the aerodynamic parameters of urban surfaces and hence the vertical wind profiles found above. The knowledge of these parameters has applications in numerous fields, such as dispersion modelling, wind loading calculations, and estimating the wind energy resource at urban locations. Using quasi-empirical modelling, we estimate the dependence of the aerodynamic roughness length and zero-plane displacement for idealized urban surfaces, on the two most significant geometrical characteristics; surface area density and building height variability. A validation of the spatially-averaged, logarithmic wind profiles predicted by the model is carried out, via comparisons with available wind-tunnel and numerical data for arrays of square based blocks of uniform and heterogeneous heights. The model predicts two important properties of the aerodynamic parameters of surfaces of heterogeneous heights that have been suggested by experiments. Firstly, the zero-plane displacement of a heterogeneous array can exceed the surface mean building height significantly. Secondly, the characteristic peak in roughness length with respect to surface area density becomes much softer for heterogeneous arrays compared to uniform arrays, since a variation in building height can prevent a skimming flow regime from occurring. Overall the simple model performs well against available experimental data and may offer more accurate estimates of surface aerodynamic parameters for complex urban surfaces compared to models that do not include height variability.     

A Matrix Approach Coupled with Monte Carlo Techniques for Solving the Net Radiative Balance of the Urban Block

A new method is developed for solving the shortwave and longwave net radiative balance of a three-dimensional urban structure, represented by parallelepiped blocks uniformly distributed in each direction. The method is based on a novel approach to determine the shape factors among surfaces, which are estimated by Monte Carlo techniques due to the complex geometry associated with the three-dimensional urban structure. Then, a set of linear equations is solved to quantify the radiative balance, in order to obtain their exact solution, considering all the inter-reflections among surfaces. The comparison between the new and the ray-tracing tracking methods resulted in a Pearson correlation coefficient of 0.996. However, by integrating the linear equations’ exact solution with Monte Carlo techniques, the new method reduces by a factor of 36 the central processing unit (CPU) time used to perform the calculations of the ray-tracing tracking method. The use of the model for a sensitivity study allows us to verify the effective absorptance and emittance increases with the canyon aspect ratio of the urban layout. An urban structure formed by square cross-sectional blocks absorbs more solar radiation than an urban structure formed by rectangular cross-sectional blocks. The approximation of a specific geometry for an equivalent bi-dimensional infinite street can be applied for rectangular cross-sectional blocks, where the width is 11 times or more greater than the depth dimension.     

Influence of Building Areal Density and Roof Shape on the Wind Characteristics Above a Town

Flow characteristics in the lower part of theatmospheric boundary layer developing immediatelyabove building roofs have been studied by physicalmodelling under neutral stratification conditions. Thevertical profiles of velocity, turbulence intensityand Reynolds stress were measured in detail above amodel urban fetch consisting of parallel streetcanyons. Two different street densities and roofshapes were tested. It is found that the influence ofthe buildings on the oncoming wind remains confined towithin three overall building heights above ground.Furthermore, the effect on the wind at roof levelfrom the areal building density is relatively weak, butstrong from the roof shape. Thus, altering roof shapecan have a much more beneficial impact on urban airquality than increasing the spacing betweenbuildings. Moreover, these findings yield a novelmethodology for reliable prediction of urban airquality, by combining numerical mesoscale wind flowmodels with physical street canyon pollutiondispersion models.     

Characteristics of the Drag Coefficient in the Roughness Sublayer over a Complex Urban Surface

The statistics of momentum exchange in the urban roughness sublayer are investigated. The analysis focuses on the characteristics of the dimensionless friction velocity, \({u_{*}}/U\) , which is defined as the square root of the drag coefficient. The turbulence observations were made at a height of 47 m above the ground on the 325-m meteorological tower, which is located in a very inhomogeneous urban area in Beijing. Under neutral conditions, the dependence of the drag coefficient on wind speed varies with wind direction. When the airflow is from the area of densely built-up buildings, the drag coefficient does not vary with wind speed, while when the airflow is from the area covered by vegetation, the drag coefficient appears to decrease with increasing wind speed. Also, the drag coefficient does not vary monotonically with the atmospheric stability. Both increasing stability and increasing instability lead to the decrease of the drag coefficient, implying that the roughness length and zero-plane displacement may vary in urban areas.     

Radiative Exchange in an Urban Street Canyon

The influence of building geometry on the radiation terms ofthe surface energy balance is a principal reason for surfacetemperature differences between rural and urban areas.Methods exist to calculate the radiation balance in an urban area,but their validity across the range of urban geometries andmaterials has not been carefully considered.Here the exchange of diffuse radiation in an urban street canyon isinvestigated using a method incorporating all reflections of radiation.This exact solution is compared to two commonly used approximationsthat retain either no reflections, or just one reflection of radiation.The area-averaged net radiative flux density from the facets of the canyondecreases in magnitude monotonically as the canyon aspect ratio increases.The two approximate solutions possess unphysical differences from thismonotonic decrease for high canyon aspect ratios or low materialemissivities/high material albedos.The errors of the two approximate solutions are small for near blackbodymaterials and small canyon aspect ratios but can be an order ofmagnitude for intermediate material properties and deep street canyons.Urban street canyon models need to consider at least one reflectionof radiation and multiple reflections are desirable for full applicability.     

Investigation of the Flow Structure in Step-Up Street Canyons—Mean Flow and Turbulence Statistics

A step-up street canyon is a characteristic urban element composed of two buildings in which the height of the upwind building ( $H_\mathrm{u}$ ) is less than the height of the downwind building ( $H_\mathrm{d}$ ). Here, the effect of canyon geometry on the flow structure in isolated step-up street canyons is investigated through isothermal wind-tunnel measurements. The measurements were acquired along the vertical symmetry plane of model buildings using two-dimensional particle image velocimetry (PIV) for normal approach flow. The building-height ratios considered were: $H_\mathrm{d}/ H_\mathrm{u} \approx 3$ , and $H_\mathrm{d}/ H_\mathrm{u} \approx 1.67$ . For each building-height ratio, the along-wind lengths (L) of the upwind and downwind buildings, and the street-canyon width (S) were kept constant, with $L \approx S$ . The cross-wind widths (W) of the upwind and downwind buildings were varied uniformly from $W/S \approx 1$ through $W/S \approx 4$ , in increments of $W/S \approx 1$ . The objective of the work was to characterize the changes in the flow structure in step-up canyons as a function of W/S, for fixed L, S, and $H_\mathrm{d}/H_\mathrm{u}$ values. The results indicate that the in-canyon flow structure does not vary significantly for $H_\mathrm{d}/H_\mathrm{u} \approx 3$ for the W/S values considered. Qualitatively, for $H_\mathrm{d}/H_\mathrm{u} \approx 3$ , the upwind building behaves as an obstacle in the upwind cavity of the downwind building. In contrast, the flow patterns observed for the $H_\mathrm{d}/H_\mathrm{u} \approx 1.67$ configurations are unique and counter-intuitive, and depend strongly on building width (W/S). For $W/S \approx 1$ and $W/S \approx 2$ , the effect of lateral flow into the canyon is so prominent that even the mean flow patterns are highly ambiguous. For $W/S \approx 3$ and 4, the flow along the vertical symmetry plane is more shielded from the lateral flow, and hence a stable counter-rotating vortex pair is observed in the canyon. In addition to these qualitative features, a quantitative analysis of the mean flow field and turbulence stress field is presented.     

The Budget of Turbulent Kinetic Energy in the Urban Roughness Sublayer

Full-scale observations from two urban sites in Basel, Switzerland were analysed to identify the magnitude of different processes that create, relocate, and dissipate turbulent kinetic energy (TKE) in the urban atmosphere. Two towers equipped with a profile of six ultrasonic anemometers each sampled the flow in the urban roughness sublayer, i.e. from street canyon base up to roughly 2.5 times the mean building height. This observational study suggests a conceptual division of the urban roughness sublayer into three layers: (1) the layer above the highest roofs, where local buoyancy production and local shear production of TKE are counterbalanced by local viscous dissipation rate and scaled turbulence statistics are close to to surface-layer values; (2) the layer around mean building height with a distinct inflexional mean wind profile, a strong shear and wake production of TKE, a more efficient turbulent exchange of momentum, and a notable export of TKE by transport processes; (3) the lower street canyon with imported TKE by transport processes and negligible local production. Averaged integral velocity variances vary significantly with height in the urban roughness sublayer and reflect the driving processes that create or relocate TKE at a particular height. The observed profiles of the terms of the TKE budget and the velocity variances show many similarities to observations within and above vegetation canopies.     

Scalar Fluxes from Urban Street Canyons Part II: Model

A practical model is developed for the vertical flux of a scalar, such as heat, from an urban street canyon that accounts for variations of the flow and turbulence with canyon geometry. The model gives the magnitude and geometric dependence of the flux from each facet of the urban street canyon, and is shown to agree well with wind-tunnel measurements described in Part I. The geometric dependence of the flux from an urban street canyon is shown to be determined by two physical processes. Firstly, as the height-to-width ratio of the street canyon increases, so does the roughness length and displacement height of the surface. This increase leads to a reduction in the wind speed in the inertial sublayer above the street canyons. Since the speed of the circulations in the street are proportional to this inertial sublayer wind speed, the flux then reduces with the inertial sublayer wind speed. This process is dominant at low height-to-width ratios. Secondly, the character of the circulations within the street canyon also varies as the height-to-width ratio increases. The flow in the street is partitioned into a recirculation region and a ventilated region. When the street canyon has high height-to-width ratios the recirculation region occupies the whole street canyon and the wind speeds within the street are low. This tendency decreases the flux at high height-to-width ratios. These processes tend to reduce the flux density from the individual facets of the street canyon, when compared to the flux density from a horizontal surface of the same material. But the street canyon has an increased total surface area, which means that the total flux from the street canyon is larger than from a horizontal surface. The variations in scalar flux from an urban street canyon with geometry is over a factor of two, which means that the physical mechanisms responsible should be incorporated into energy balance models for urban areas.     

An Urban Surface Exchange Parameterisation for Mesoscale Models

A scheme to represent the impact of urban buildings on airflow in mesoscale atmospheric models is presented. In the scheme, the buildings are not explicitly resolved, but their effects on the grid-averaged variables are parameterised. An urban quarter is characterised by a horizontal building size, a street canyon width and a building density as a function of height. The module computes the impact of the horizontal (roof and canyon floor) and vertical (walls) surfaces on the wind speed, temperature and turbulent kinetic energy. The computation of the shortwave and longwave radiation, needed to compute the temperature of the urban surfaces, takes into account the shadowing and radiation trapping effects induced by the urban canyons. The computation of the turbulent length scales in the TKE equation is also modified to take into account the presence of the buildings.The parameterisation is introduced into a mesoscale model and tested in a bidimensional case of a city over flat terrain. The new parameterisation is shown to be able to reproduce the most important features observed in urban areas better than the traditional approach which is based only on the modification of the roughness length, thereby retaining the Monin–Obukhov similarity theory. The new surface exchange parameterisation is furthermore shown to have a strong impact on the dispersion characteristics of air pollutants in urban areas.