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.     

Large-eddy Simulation of Turbulent Flow across a Forest Edge. Part II: Momentum and Turbulent Kinetic Energy Budgets

Momentum and turbulent kinetic energy (TKE) budgets across a forest edge have been investigated using large-eddy simulation (LES). Edge effects are observed in the rapid variation of a number of budget terms across this vegetation transition. The enhanced drag force at the forest edge is largely balanced by the pressure gradient force and by streamwise advection of upstream momentum, while vertical turbulent diffusion is relatively insignificant. For variance and TKE budgets, the most important processes at the forest edge are production due to the convergence (or divergence) of the mean flow, streamwise advection, pressure diffusion and enhanced dissipation by canopy drag. Turbulent diffusion, pressure redistribution and vertical shear production, which are characteristic processes in homogeneous canopy flow, are less important at the forest transition. We demonstrate that, in the equilibrated canopy flow, a substantial amount of TKE produced in the streamwise direction by the vertical shear of the mean flow is redistributed in the vertical direction by pressure fluctuations. This redistribution process occurs in the upper canopy layers. Part of the TKE in the vertical velocity component is transferred by turbulent and pressure diffusion to the lower canopy levels, where pressure redistribution takes place again and feeds TKE back to the streamwise direction. In this TKE cycle, the primary source terms are vertical shear production for streamwise velocity variance and pressure redistribution for vertical velocity variance. The evolution of these primary source terms downwind of the forest edge largely controls the adjustment rates of velocity variances.     

Simulating Australian Urban Climate in a Mesoscale Atmospheric Numerical Model

We develop an urban canopy scheme coupled to a mesoscale atmospheric numerical model and evaluate the simulated climate of an Australian city. The urban canopy scheme is based on the Town Energy Budget approach, but is modified to efficiently represent the predominately suburban component of Australian cities in regional climate simulations. Energy conservation is improved by adding a simple model of air-conditioning to prevent the urban parametrization acting as an energy sink during the Australian summer. In-canyon vegetation for suburban areas is represented by a big-leaf model, but with a largely reduced set of prognostic variables compared to previous approaches. Although we have used a recirculation/venting based parametrization of in-canyon turbulent heat fluxes that employs two canyon wall energy budgets, we avoid using a fixed canyon orientation by averaging the canyon fluxes after integrating over 180° of possible canyon orientations. The urban canopy scheme is evaluated by simulating the climate for Melbourne, Australia after coupling it to The Air Pollution Model. The combined system was found to predict a realistic climatology of air temperatures and winds when compared with observations from Environmental Protection Authority monitoring stations. The model also produced a plausible partitioning of the urban energy budget when compared to urban flux-tower studies. Overall, the urban canyon parametrization appears to have reasonable potential for studying present and predicting changes in future Australian urban climates in regional climate simulations.     

Turbulent Transfer Between Street Canyons and the Overlying Atmospheric Boundary Layer

The turbulent exchange of momentum between a two-dimensional cavity and the overlying boundary layer has been studied experimentally, using hot-wire anemometry and particle image velocimetry (PIV). Conditions within the boundary layer were varied by changing the width of the canyons upstream of the test canyon, whilst maintaining the square geometry of the test canyon. The results show that turbulent transfer is due to the coupling between the instabilities generated in the shear layer above the canyons and the turbulent structures in the oncoming boundary layer. As a result, there is no single, unique velocity scale that correctly characterizes all the processes involved in the turbulent exchange of momentum across the boundary layer. Similarly, there is no single velocity scale that can characterize the different properties of the turbulent flow within the canyon, which depends strongly on the way in which turbulence from the outer flow is entrained into the cavity and carried round by the mean flow. The results from this study will be useful in developing simple parametrizations for momentum exchange in the urban canopy, in situations where the street geometry consists principally of relatively long, uniform streets arranged in grid-like patterns; they are unlikely to be applicable to sparse geometries composed of isolated three-dimensional obstacles.     

Large-Eddy Simulation of Flow and Pollutant Transport in Urban Street Canyons with Ground Heating

Our study employed large-eddy simulation (LES) based on a one-equation subgrid-scale model to investigate the flow field and pollutant dispersion characteristics inside urban street canyons. Unstable thermal stratification was produced by heating the ground of the street canyon. Using the Boussinesq approximation, thermal buoyancy forces were taken into account in both the Navier–Stokes equations and the transport equation for subgrid-scale turbulent kinetic energy (TKE). The LESs were validated against experimental data obtained in wind-tunnel studies before the model was applied to study the detailed turbulence, temperature, and pollutant dispersion characteristics in the street canyon of aspect ratio 1. The effects of different Richardson numbers (Ri) were investigated. The ground heating significantly enhanced mean flow, turbulence, and pollutant flux inside the street canyon, but weakened the shear at the roof level. The mean flow was observed to be no longer isolated from the free stream and fresh air could be entrained into the street canyon at the roof-level leeward corner. Weighed against higher temperature, the ground heating facilitated pollutant removal from the street canyon.     

Boundary layer structure over an inhomogeneous surface: Simulation with a non-hydrostatic mesoscale model

A three-dimensional, non-hydrostatic mesoscale model is used to study boundary-layer structure over an area characterized by the city of Copenhagen, the Øresund strait, and adjacent coastal farmland. Simulations are compared with data obtained on June 5, 1984 during the Øresund experiment.Under moderately strong wind conditions, a stable internal boundary layer (IBL) developed over the Øresund strait during the day. Near-surface winds decelerate over water due to diminished vertical momentum transfer.The turbulent kinetic energy field closely reflects the surface roughness distribution due to the imposed relatively strong wind forcing. TKE budgets over water, farmland and a city area are discussed by inspection of vertical profiles of the individual terms. The buoyancy term is used to indicate IBL heights because it changes sign at the boundary between different stability regimes. Measured and simulated dissipation rates show a decrease in the transition zone as the air travels over water and an abrupt increase when the IBL over a downwind city area is intersected. The top of the stable IBL is characterized by a minimum in the vertical TKE profile.     

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.     

Intercomparison of three urban climate models

An intercomparison of the surface energy budgets from three urban climate models was made to assess the comparability of results, and to evaluate the surface energy fluxes from each model. The three models selected spanned the continuum of approaches currently employed in the treatment of the effects of urban geometry. The first model was an urban canopy-layer model which explicitly examined urban canyon geometry. The second model treated the city as a warm, rough, moist plate but included greatly simplified parameterizations of urban geometry. Neither model included a dynamic link to the urban boundary-layer. The third model was a one-dimensional urban boundary-layer model which utilized a simple warm, rough, moist plate approach but included a dynamic coupling of the urban surface layer to the urban boundary-layer.Results showed considerable disagreement between the three models in regards to the individual energy fluxes. Average rankings of the energy fluxes in terms of comparability from high-to-low similarity were: (1) solar radiation, (2) sensible heat flux, (3) conduction, (4) latent heat flux, (5) longwave re-radiation, and (6) longwave radiation input. In general, the urban canopy-layer model provided more realistic results, although each model demonstrated strong and weak points. Results indicate that current urban boundary-layer models may produce surface energy budgets with lower sensible heat fluxes and substantially higher latent heat fluxes than is supported by field evidence from the literature.     

Relative influence of vegetation on urban energy budgets and surface temperatures

Two deterministic models were combined: one for canopy leaf energy budgets and one for street canyon energy budgets. The effects of street parks and roof gardens in contrast to non-vegetated city blocks were examined by the use of four typical urban morphologies, which were exposed latitudinally to summer and winter simulations. A variety of increases and decreases in shortwave radiation, net radiation, sensible heat flux, and system reradiation resulted. These changes appear to represent the generalized limits of the possible responses to the addition of vegetation to non-vegetated city blocks.     

Large-Eddy Simulation of Turbulent Organized Structures within and above Explicitly Resolved Cube Arrays

Large-eddy simulations have been performed for fully developed turbulent flow within and above explicitly resolved simple cube arrays. The results from our model, hereafter LES-CITY, are shown to agree with laboratory experiments. We investigated the systematic influence of cube density on turbulent flow characteristics by performing numerical experiments for cube areal densities from 0 to 44%. The following results were obtained: (1) The dispersive momentum flux was quite large within the canopy layer due to a mean stream re-circulation, whereas it was smaller above the canopy. The spatial variation of temporally averaged momentum in the roughness sub-layer was 20% or less of the total kinematic surface drag. (2) The temporally and spatially-averaged flow structure confirmed the existence of conventionally described canyon flow regimes; isolated, interfacial, and wake. However, the intermittency of the canyon flow for all cube densities was quite large and the stream patterns were never persistent. (3) Turbulent organized structures (TOS) similar to those observed in turbulent surface-layer flows were simulated, which are characterized by longitudinally-elongated low speed streaks and the corresponding shorter streamwise vortices. The streaks in sparse and dense canopy flows were likely to be aligned to the street line and to the roof lines, respectively. Such heterogeneity of TOS partially accounts for the large spatial variation of momentum flux. (4) In contrast to the mixing layer analogy of vegetation flows, the TOS and the resulting turbulent statistics of urban flow above the canopy resembled those in surface layers. The recirculation within the canopy significantly influenced the turbulent statistical properties.     

Surface heating in relation to air temperature,wind and turbulence in an urban street canyon

Wind and temperature measurements from within and above a deep urban canyon (height/width = 2.1) were used to examine the thermal structure of air within the canyon, exchange of heat with the overlying atmosphere, and the possible impacts of surface heating on within-canyon air flow. Measurements were made over a range of seasons and primarily analysed for sunny days. This allowed the study of temperature differences between opposing canyon walls and between wall and air of more than 15°C in summer. The wall temperature patterns follow those of incoming solar radiation loading with a secondary daytime effect from the longwave exchange between the walls. In winter, the canyon walls receive little direct solar radiation, and temperature differences are largely due to anthropogenic heating of the building interiors. Cool air from aloft and heated air from canyon walls is shown to circulate within the canyon under cross-canyon flow. Roofs and some portions of walls heat up rapidly on clear days and have a large influence on heat fluxes and the temperature field. The magnitude and direction of the measured turbulent heat flux also depend strongly on the direction of flow relative to surface heating. However, these spatial differences are smoothed by the shear layer at the canyon top. Buoyancy effects from the heated walls were not seen to have as large an impact on the measured flow field as has been shown in numerical experiments. At night canyon walls are shown to be the source of positive sensible heat fluxes. The measurements show that materials and their location, as well as geometry, play a role in regulating the heat exchange between the urban surface and atmosphere.     

Thermal bioclimate in idealized urban street canyons in Campinas, Brazil

Among several urban design parameters, the height-to-width ratio (H/W) and orientation are important parameters strongly affecting thermal conditions in cities. This paper quantifies changes in thermal comfort due to typical urban canyon configurations in Campinas, Brazil, and presents urban guidelines concerning H/W ratios and green spaces to adapt urban climate change. The study focuses on thermal comfort issues of humans in urban areas and performs evaluation in terms of physiologically equivalent temperature (PET), based on long-term data. Meteorological data of air temperature, relative humidity, wind speed and solar radiation over a 7-year period (2003–2010) were used. A 3D street canyon model was designed with RayMan Pro software to simulate the influence of urban configuration on urban thermal climate. The following configurations and setups were used. The model canyon was 500 m in length, with widths 9, 21, and 44 m. Its height varied in steps of 2.5 m, from 5 to 40 m. The canyon could be rotated in steps of 15°. The results show that urban design parameters such as width, height, and orientation modify thermal conditions within street canyons. A northeast–southwest orientation can reduce PET during daytime more than other scenarios. Forestry management and green areas are recommended to promote shade on pedestrian areas and on façades, and to improve bioclimate thermal stress, in particular for H/W ratio less than 0.5. The method and results can be applied by architects and urban planners interested in developing responsive guidelines for urban climate issues.     

A Vegetated Urban Canopy Model for Meteorological and Environmental Modelling

An urban canopy model is developed for use in mesoscale meteorological and environmental modelling. The urban geometry is composed of simple homogeneous buildings characterized by the canyon aspect ratio (h/w) as well as the canyon vegetation characterized by the leaf aspect ratio (σ _l ) and leaf area density profile. Five energy exchanging surfaces (roof, wall, road, leaf, soil) are considered in the model, and energy conservation relations are applied to each component. In addition, the temperature and specific humidity of canopy air are predicted without the assumption of thermal equilibrium. For radiative transfer within the canyon, multiple reflections for shortwave radiation and one reflection for longwave radiation are considered, while the shadowing and absorption of radiation due to the canyon vegetation are computed by using the transmissivity and the leaf area density profile function. The model is evaluated using field measurements in Vancouver, British Columbia and Marseille, France. Results show that the model quite well simulates the observations of surface temperatures, canopy air temperature and specific humidity, momentum flux, net radiation, and energy partitioning into turbulent fluxes and storage heat flux. Sensitivity tests show that the canyon vegetation has a large influence not only on surface temperatures but also on the partitioning of sensible and latent heat fluxes. In addition, the surface energy balance can be affected by soil moisture content and leaf area index as well as the fraction of vegetation. These results suggest that a proper parameterization of the canyon vegetation is prerequisite for urban modelling.     

Thermal comfort in an east–west oriented street canyon in Freiburg (Germany) under hot summer conditions

Summary Field-measurements were conducted in an urban street canyon with an east–west orientation, and a height-to-width ratio H/W = 1 during cloudless summer weather in 2003 in Freiburg, Germany. This experimental work adds to the knowledge available on the microclimate of an urban canyon and its impact on human comfort. Air temperature T _a , air humidity VP, wind speed v and direction dd were measured continuously. All short-wave and long-wave radiation fluxes from the 3D surroundings were also measured. The degree of comfort was defined in terms of physiologically equivalent temperature (PET). Furthermore, the data gathered within the canyon were compared to data collected by a permanent urban climate station with the aim of furthering the understanding of microclimatic changes due to street geometry. Changes in the meteorological variables T _a , v and dd in the canyon in comparison to an unobstructed roof level location were found to be in good agreement with previous studies, i.e., a small increase of T _a in the canyon adjacent to irradiated surfaces, and a good correlation of v and dd between canyon and roof levels. The daily dynamics of canyon facet irradiances and their impacts on the heat gained by a pedestrian were strongly dependent on street geometry and orientation. Thermal stress was mostly attributable to solar exposure. Under cloudless summer weather, a standing body was found to absorb, on average, 74% of heat in the form of long-wave irradiance and 26% as short-wave irradiance. Shading the pedestrian as well as the surrounding surfaces is, hence, the first strategy in mitigating heat stress in summer under hot conditions.     

A Laboratory Model for the Flow in Urban Street Canyons Induced by Bottom Heating

Water tank experiments are carried out to investigate the convection flow induced by bottom heating and the effects of the ambient wind on the flow in non-symmetrical urban street canyons based on the PIV (Particle Image Visualization) technique. Fluid experiments show that with calm ambient wind,the flows in the street canyon are completely driven by thermal force, and the convection can reach the upper atmosphere of the street canyon. Horizontal and vertical motions also appear above the roofs of the buildings. These are the conditions which favor the exchange of momentum and air mass between the street canyon and its environment. More than two vortices are induced by the convection, and the complex circulation pattern will vary with time in a wider street canyon. However， in a narrow street canyon， just one vortex appears. With a light ambient wind, the bottom heating and the associated convection result in just one main vortex. As the ambient wind speed increases, the vortex becomes more organized and its center shifts closer to the leeward building.     

A Physically-Based Scheme For The Urban Energy Budget In Atmospheric Models

An urban surface scheme for atmospheric mesoscale models ispresented. A generalization of local canyon geometry isdefined instead of the usual bare soil formulation currently usedto represent cities in atmospheric models. This allows refinement ofthe radiative budgets as well as momentum, turbulent heat and ground fluxes.The scheme is aimed to be as general as possible, in order to representany city in the world, for any time or weather condition(heat island cooling by night, urban wake, water evaporation after rainfalland snow effects).Two main parts of the scheme are validated against published data.Firstly, it is shown that the evolution of the model-predictedfluxes during a night with calm winds is satisfactory, considering both the longwave budget and the surface temperatures. Secondly, the original shortwave scheme is tested off-line and compared to the effective albedoof a canyon scale model. These two validations show that the radiative energy input to the urban surface model is realistic.Sensitivity tests of the model are performed for one-yearsimulation periods, for both oceanic and continental climates. The scheme has the ability to retrieve, without ad hoc assumptions, the diurnal hysteresis between the turbulent heat flux and ground heat flux. It reproduces the damping of the daytime turbulent heat flux by the heat storage flux observed in city centres. The latent heat flux is negligible on average,but can be large when short time scales are considered (especially afterrainfall). It also suggests that in densely built areas, domesticheating can overwhelm the net radiation, and supply a continuous turbulentheat flux towards the atmosphere. This becomes very important inwinter for continental climates. Finally, a comparison with a vegetation scheme shows that the suburban environment can be represented with a bare soil formulation for large temporal or spatial averages (typical of globalclimatic studies), but that a surface scheme dedicated to the urban surface is necessary when smaller scales are considered: town meteorological forecasts, mesoscale or local studies.     

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.     

Transfer processes in a simulated urban street canyon

The transfer processes within and above a simulated urban street canyon were investigated in a generic manner. Computational fluid dynamics (CFD) was used to aid understanding and to produce some simple operational parameterisations. In this study we addressed specifically the commonly met situation where buoyancy effects arising from elevated surface temperatures are not important, i.e. when mechanical forces outweigh buoyancy forces. In a geophysical context this requires that some suitably defined Richardson number is small. From an engineering perspective this is interpreted as the important case when heat transfer within and above urban street canyons is by forced convection. Surprisingly, this particular scenario (for which the heat transfer coefficient between buildings and the flow is largest), has been less well studied than the situation where buoyancy effects are important. The CFD technique was compared against wind-tunnel experiments to provide model evaluation. The height-to-width ratio of the canyon was varied through the range 0.5–5 and the flow was normal to the canyon axis. By setting the canyon’s facets to have the same or different temperatures or to have a partial temperature distribution, simulations were carried out to investigate: (a) the influence of geometry on the flow and mixing within the canyon and (b) the exchange processes within the canyon and across the canyon top interface. Results showed that the vortex-type circulation and turbulence developed within the canyon produced a temperature distribution that was, essentially, spatially uniform (apart from a relatively thin near-wall thermal boundary layer) This allowed the temperatures within the street canyon to be specified by just one value T _can , the canyon temperature. The variation of T _can with wind speed, surface temperatures and geometry was extensively studied. Finally, the exchange velocity u _E across the interface between the canyon and the flow above was calculated based on a heat flux balance within the canyon and between the canyon and the flow above. Results showed that u _E was approximately 1% of a characteristic wind velocity above the street canyon. The problem of radiative exchange is not addressed but it can, of course, be introduced analytically, or computationally, when necessary.