Scalar fluxes from urban street canyons. Part I: Laboratory simulation

Flow over urban surfaces depends on surface morphology and interaction with the boundary layer above. However, the effect of the flow on scalar fluxes is hard to quantify. The naphthalene sublimation technique was used to quantify scalar vertical fluxes out of a street canyon under neutral conditions. For an array of eight canyons with aspect ratio H/W=0.75 (here, H is building height and W is the street width), increased flux was observed in the first two or three canyons for moderate and low roughness upstream. This is consistent with predictions of the length scale for initial adjustment of flow to an urban canopy. The flux was constant after the initial adjustment region and thus dependent only on local geometry. For a street canyon in the equilibrium part of the array, each facet of the street canyon was coated with naphthalene to simulate scalar release from street, walls and roof, to evaluate the effect of street canyon geometry on fluxes for H/W=0.25, 0.6, 1 and 2. Fluxes from the roof and downstream wall were considerably larger than fluxes from the street and upstream wall, and only the flux from the downstream wall exhibited a simple decrease with H/W. For each H/W there was a monotonic decrease between downstream wall, street and upstream wall transfer. This suggests that flow decelerates around the recirculation region in the lee of the upstream building, i.e. a recirculating jet rather than a symmetrical vortex. The addition of a second source within the street canyon resulted in reduced fluxes from each facet for H/W>0.25, due to increased concentration of naphthalene in the canyon air.     

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.     

Pollutant Concentrations in Street Canyons of Different Aspect Ratio with Avenues of Trees for Various Wind Directions

This study summarizes the effects of avenues of trees in urban street canyons on traffic pollutant dispersion. We describe various wind-tunnel experiments with different tree-avenue models in combination with variations in street-canyon aspect ratio W/H (with W the street-canyon width and H the building height) and approaching wind direction. Compared to tree-free street canyons, in general, higher pollutant concentrations are found. Avenues of trees do not suppress canyon vortices, although the air ventilation in canyons is hindered significantly. For a perpendicular wind direction, increases in wall-average and wall-maximum concentrations at the leeward canyon wall and decreases in wall-average concentrations at the windward wall are found. For oblique and perpendicular wind directions, increases at both canyon walls are obtained. The strongest effects of avenues of trees on traffic pollutant dispersion are observed for oblique wind directions for which also the largest concentrations at the canyon walls are found. Thus, the prevailing assumption that attributes the most harmful dispersion conditions to a perpendicular wind direction does not hold for street canyons with avenues of trees. Furthermore, following dimensional analysis, an estimate of the normalized wall-maximum traffic pollutant concentration in street canyons with avenues of trees is derived.     

On the vertical profile of surface radiative fluxes in southwest Germany

Summary Vertical profile of surface radiative fluxes in an area of heterogeneous terrain in south-west Germany is presented. Main data sets utilized for the study were recorded during the REgio KLIma Projekt (REKLIP). Supporting observational data were provided by the German weather service and German geophysical consultant service. Elevation of the study sites ranged from 212 m a.s.l. to 1489 m a.s.l. From May to September, monthly mean albedo was generally low at the study sites, ranging from 19% to 24%. For the other months, monthly mean albedo lie between 22% and 25% at the lowland site but extended between 27% and 71% at the highly elevated mountain site. Following the altitudinal increase in surface albedo, net radiative flux and radiation efficiency declined with elevation at an annual mean of 1.15 Wm⁻²/100 m and 0.008/100 m respectively. Absorbed shortwave radiation and effective terrestrial radiation showed mean decline of 1.54 Wm⁻²/100 m and 0.34 Wm⁻²/100 m, respectively, with the mean sky-to-earth radiation deficit amounting to about 52 Wm⁻² for the lowland site and 73 Wm⁻² for the highest elevated site. Some empirical models which express shortwave and longwave radiative fluxes in terms of meteorological variables have been validated for the lowland and mountain sites. Monthly mean daily total estimates of solar radiation obtained from ?ngst?m-Prescott relation were quite consistent with observed values. Parameterisation of downward atmospheric radiation under all sky condition was achieved by extending Brutsaert clear sky atmospheric model. Relationship between outgoing longwave radiation and screen temperature at the study sites was best described by an exponential function unlike the linear relationship proposed by Monteith and Unsworth. Net radiative flux for the lowland and mountain sites has been expressed in terms of absorbed shortwave radiation, cloud amount and screen temperature. Received March 5, 2001 Revised October 29, 2001     

Mass Transfer Velocity and Momentum Vertical Exchange in Simulated Deep Street Canyons

A box model to simulate mass transfer inside deep street canyons and with atmospheric flow above is introduced and discussed. Two ideal deep street canyons with aspect ratios of 3 and 5 (the aspect ratio being the ratio between building height and street width H/W) are considered. This range of aspect ratios, found in many densely populated historical centres in Mediterranean cities as well as in other cities around the world, potentially creates high air pollutant concentration levels. Our model is based on a combination of analytical solutions and computation fluid dynamics (CFD) simulations using carbon monoxide (CO) as a tracer pollutant. The analytical part of the model is based on mass transfer velocity concepts while CFD simulations are used both for a preliminary validation of the physical hypothesis underlying the model (steady-state simulations) and to evaluate the concentration pattern with time (transient or wash-out simulations). Wash-out simulation curves were fitted by model curves, and mass transfer velocities were evaluated through a best-fitting procedure. Upon introducing into the model the contribution of traffic-produced turbulence, the modelled CO concentration levels became comparable with those obtained in real-world monitoring campaigns. The mass transfer rate between the canyon and the above atmosphere was then expressed in terms of an overall mass transfer velocity, which directly allows the evaluation of the mass transfer rate between the bottom volume of the canyon (pedestrian level) with the above atmosphere. Overall mass transfer velocities are reported as a function of the operating conditions studied (H/W = 3–5 and wind speeds = 2–8 ms⁻¹). Finally, a simple expression is reported for determining pollutant concentrations at the pedestrian level based on the overall mass transfer velocity defined.     

Temporal variations in heat fluxes over a central European city centre

Summary Energy fluxes have been measured over an area near the centre of the city of Łódź, Poland, since November 2000. The site was selected because the building style (surface cover and morphology) is typical of European cities, yet distinct from the majority of cities where energy balance observations have been studied thus far. The multi-year dataset permits consideration of temporal changes in energy balance partitioning over a wide range of seasonal and synoptic conditions and of the role of heat storage and anthropogenic fluxes in the energy balance. Partitioning of net radiation into the turbulent fluxes is consistent in the two years, with the largest differences occurring due to differing precipitation. The monthly ensemble diurnal cycles of the turbulent fluxes over the two years are similar. The largest differences occur during the July–September period, and are attributable to greater net radiation and lower rainfall in 2002. The latent heat flux accounts for approximately 40% of the turbulent heat transfer on an annual basis. The average daily daytime Bowen ratio and its variability are slightly reduced during the summer (growing) season. Anthropogenic heat is a significant input to the urban energy balance in the winter. The fluxes observed in this study are consistent with results from other urban sites.     

Flow in a Street Canyon for any External Wind Direction

An analytical model has been developed for the flow along a street canyon (of height H and width W), generated by an external wind blowing at any angle relative to the axis of the street. Initially, we consider the special case of a wind blowing parallel to the street. The interior of the street is decomposed into three regions, and the flow within each region is assumed to depend only on the external wind and the distance to the closest solid boundary. This decomposition leads to two different flow regimes: one for narrow streets (H/W > 1/2) and one for wide streets (H/W < 1/2). The theoretical model agrees well with results obtained from numerical simulations using a Reynolds-Averaged Navier–Stokes model. We then generalize the model to the case of arbitrary wind direction. Numerical solutions show that the streamlines of the mean flow in the street have a spiral form, and for most angles of incidence, the mass flux along the street scales on the component of the external wind resolved parallel to the street. We use this result to generalize the model derived for wind blowing parallel to the street, and the results from this model agree well with the numerical simulations. The model that has been developed can be evaluated rapidly using only very modest computing power, so it is suitable for use as an operational tool.     

A Double-Canyon Radiation Scheme for Multi-Layer Urban Canopy Models

We develop a double-canyon radiation scheme (DCEP) for urban canopy models embedded in mesoscale numerical models based on the Building Effect Parametrization (BEP). The new scheme calculates the incoming and outgoing longwave and shortwave radiation for roof, wall and ground surfaces for an urban street canyon characterized by its street and building width, canyon length, and the building height distribution. The scheme introduces the radiative interaction of two neighbouring urban canyons allowing the full inclusion of roofs into the radiation exchange both inside the canyon and with the sky. In contrast to BEP, we also treat direct and diffuse shortwave radiation from the sky independently, thus allowing calculation of the effective parameters representing the urban diffuse and direct shortwave radiation budget inside the mesoscale model. Furthermore, we close the energy balance of incoming longwave and diffuse shortwave radiation from the sky, so that the new scheme is physically more consistent than the BEP scheme. Sensitivity tests show that these modifications are important for urban regions with a large variety of building heights. The evaluation against data from the Basel Urban Boundary Layer Experiment indicates a good performance of the DCEP when coupled with the regional weather and climate model COSMO-CLM.     

Climate feedbacks in a general circulation model incorporating prognostic clouds

 This study performs a comprehensive feedback analysis on the Bureau of Meteorology Research Centre General Circulation Model, quantifying all important feedbacks operating under an increase in atmospheric CO₂. The individual feedbacks are analysed in detail, using an offline radiation perturbation method, looking at long- and shortwave components, latitudinal distributions, cloud impacts, non-linearities under 2xCO₂ and 4xCO₂ warmings and at interannual variability. The water vapour feedback is divided into terms due to moisture height and amount changes. The net cloud feedback is separated into terms due to cloud amount, height, water content, water phase, physical thickness and convective cloud fraction. Globally the most important feedbacks were found to be (from strongest positive to strongest negative) those due to water vapour, clouds, surface albedo, lapse rate and surface temperature. For the longwave (LW) response the most important term of the cloud ‘optical property’ feedbacks is due to the water content. In the shortwave (SW), both water content and water phase changes are important. Cloud amount and height terms are also important for both LW and SW. Feedbacks due to physical cloud thickness and convective cloud fraction are found to be relatively small. All cloud component feedbacks (other than height) produce conflicting LW/SW feedbacks in the model. Furthermore, the optical property and cloud fraction feedbacks are also of opposite sign. The result is that the net cloud feedback is the (relatively small) product of conflicting physical processes. Non-linearities in the feedbacks are found to be relatively small for all but the surface albedo response and some cloud component contributions. The cloud impact on non-cloud feedbacks is also discussed: greatest impact is on the surface albedo, but impact on water vapour feedback is also significant. The analysis method here proves to be a␣powerful tool for detailing the contributions from different model processes (and particularly those of the clouds) to the final climate model sensitivity. Received: 15 June 2000 / Accepted: 10 January 2001     

A Multi-layer Radiation Model for Urban Neighbourhoods with Trees

A neighbourhood-scale multi-layer urban canopy model of shortwave and longwave radiation exchange that explicitly includes the radiative effects of tall vegetation (trees) is presented. Tree foliage is permitted both between and above buildings, and mutual shading, emission and reflection between buildings and trees are included. The basic geometry is a two-dimensional canyon with leaf area density profiles and probabilistic variation of building height. Furthermore, the model accounts for three-dimensional path lengths through the foliage. Ray tracing determines the receipt of direct shortwave irradiance by building and foliage elements. View factors for longwave and shortwave diffuse radiation exchange are computed once at the start of the simulation using a Monte Carlo ray tracing approach; for subsequent model timesteps, matrix inversion rapidly solves infinite reflections and interception of emitted longwave between all elements. The model is designed to simulate any combination of shortwave and longwave radiation frequency bands, and to be portable to any neighbourhood-scale urban canopy geometry based on the urban canyon. Additionally, the model is sufficiently flexible to represent forest and forest-clearing scenarios. Model sensitivity tests demonstrate the model is robust and computationally feasible, and highlight the importance of vertical resolution to the performance of urban canopy radiation models. Full model evaluation is limited by the paucity of within-canyon radiation measurements in urban neighbourhoods with trees. Where appropriate model components are tested against analytic relations and results from an independent urban radiation transfer model. Furthermore, system response tests demonstrate the ability of the model to realistically distribute shortwave radiation among urban elements as a function of built form, solar angle and tree foliage height, density and clumping. Separate modelling of photosynthetically-active and near-infrared shortwave bands is shown to be important in some cases. Increased canyon height-to-width ratio and/or tree cover diminishes the net longwave radiation loss of individual canyon elements (e.g., floor, walls), but, notably, has little effect on the net longwave loss of the whole urban canopy. When combined with parametrizations for the impacts of trees on airflow and hydrological processes in the urban surface layer, the new radiation model extends the applicability of urban canopy models and permits more robust assessment of trees as tools to manage urban climate, air quality, human comfort and building energy loads.     

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.     

The sea surface temperature deviation and the heat flow at the sea-air interface

The deviation of the sea surface temperature from the water temperature below is calculated as a function of the heat flow through the air-sea interface, using wind tunnel determinations of the effective thermal diffusivity in a boundary layer. The influence ofQ, shortwave radiation, andH, latent and sensible heat transfer plus effective back radiation, and U, wind speed, can be described by:T ₀ –T _w =C ₁ ·H/U +C ₂ ·Q/U. The calculated coefficients vary slightly with reference depth, Tables II and III. They are in good agreement with independent observations.On leave at Department of Oceanography, Oregon State University, Corvallis, Oregon in 1969–70.     

Urban–rural contrasts of meteorological parameters in Łódź

Summary Data from two automatic stations in Łódź (one urban and one rural) for the period 1997–2002 are analyzed to reveal urban–rural contrasts of such parameters as air temperature, relative humidity, water vapour pressure and wind speed. Under favourable weather conditions the highest temperature differences between the urban and rural station exceeds 8 °C. Relative humidity is lower in the town, sometimes by more than 40%. Water vapour pressure differences can be either positive (up to 5 hPa) or negative (up to −4 hPa). Wind speed at the urban station is on average lower by about 34% in night and 39% during daytime. Regression analysis shows that for rural winds lower than 1.13 m s⁻¹ urban winds can be stronger than rural speeds. Attention has also been paid to singularities in the course of the analyzed parameters over 24 hour periods. It is shown that the typical course of the urban heat island intensity under favourable conditions is similar in all season. Four stages of this course have been distinguished. Wind speed differences also seem to change in a typical way. Case studies show that humidity contrasts, unlike temperature, can evolve in different ways under fine weather conditions. Types of relative humidity evolution are proposed.     

Rapid determination of canyon geometry parameters for use in surface radiation budgets

Summary  The parameterisation of canyon geometry as a sky-view factor is common practice when studying surface radiation budgets. However, the accurate computation of sky-view factors can prove to be slow and has often resulted in studies involving the sky-view factor to be limited to just a few study sites. A fully automated digital approach to calculating sky-view factors direct from digital imagery is presented along with a simple methodology to determine the impact of surface geometry on the nature of incoming shortwave radiation. The validity of these new automated techniques is tested by a comparative study with alternative methodologies, producing nearly identical results. It is hoped that the rapidity of the techniques should allow for the impacts of canyon geometry on surface radiation budgets to be studied on a scale not previously possible. Received July 6, 2000 Revised December 11, 2000     

Investigation about the dependence of spectral diffuse-to-direct-beam irradiance ratio on atmospheric turbidity and solar zenith angle

Summary The modifications of the solar spectral diffuse and direct-beam irradiances as well as the diffuse-to-direct-beam ratio, E_dλ/E_bλ, as a function of the aerosol optical depth, AOD, and solar zenith angle, SZA, is investigated. The E_dλ/E_bλ ratios decrease rapidly with wavelength and exponential curves in the form E_dλ/E_bλ = aλ^−b can be fitted with a great accuracy. These curves are strongly modified by the solar spectrum distribution, which is affected by the aerosol loading, aerosol optical properties and SZA. The spectral dependence of the above E_dλ/E_bλ ratios in logarithmic coordinates does not yield a straight line, while a significant departure from the linearity is revealed. The reasons for this departure are investigated in detail and it is established that the aerosol physical properties such as single scattering albedo and size distribution along with the effect of SZA are responsible. These parameters strongly affect the scattering processes in the atmosphere and as a consequence the diffuse spectral distribution. The E_dλ/E_bλ ratio, which is an indicator of the atmospheric transmittance (King, 1979), exhibits a strong wavelength and aerosol-loading dependence. The observed differences between turbid and clear atmospheres constitute a manifestation of contrasting air properties and influence solar irradiance spectra. The present work aims at investigating the effect of atmospheric turbidity and SZA on the E_dλ/E_bλ ratio. For this reason, two distinct cases are examined: one having different atmospheric turbidity conditions but same SZA and a second having different SZAs and same atmospheric turbidity levels.     

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.     

Use of atmospheric radiation measurement program data from Barrow, Alaska, for evaluation and development of snow-albedo parameterizations

Summary Snow albedo is determined from the ratio of out-going to incoming solar radiation using three years of broadband shortwave radiometer data obtained from the Barrow, Alaska, Atmospheric Radiation Measurement (ARM) site. These data are used for the evaluation of various types of snow-albedo parameterizations applied in numerical weather prediction or climate models. These snow-albedo parameterizations are based on environmental conditions (e.g., air or snow temperature), snow related characteristics (e.g., snow depth, snow age), or combinations of both. The ARM data proved to be well suited for snow-albedo evaluation purposes for a low-precipitation tundra environment. The evaluation confirms that snow-age dependent parameterizations of snow albedo work well during snowmelt, while parameterizations considering meteorological conditions often perform better during snow accumulation. Current difficulties in parameterizing snow albedo occur for long episodes of snow-event free conditions and episodes with a high frequency of snow events or strong snowfall. In a further step, the first two years of the ARM albedo dataset is used to develop a snow-albedo parameterization, and the third year’s data serves for its evaluation. This parameterization considers snow depth, wind speed, and air temperature which are found to be significant parameters for snow-albedo modeling under various conditions. Comparison of all evaluated snow-albedo parameterizations with this new parameterization shows improved snow-albedo prediction. Correspondence: Nicole M?lders, Geophysical Institute and College of Natural Science and Mathematics, University of Alaska Fairbanks, 903 Koyukuk Drive, P.O. Box 757320, Fairbanks, AK 99775-7320, USA     

All-Sky Surface Radiation and Clear-Sky Surface Energy Budgets: Summer to Freeze-Up in the Western Maritime Arctic

The 2009 ArcticNet expedition was a field campaign in the Amundsen Gulf–eastern Beaufort Sea region from mid-July to the beginning of November aboard the CCGS Amundsen that provided an opportunity to describe the all-sky surface radiation and the clear-sky surface energy budgets from summer to freeze-up in the data sparse western maritime Arctic. Because the fractional area of open water was generally larger than the fractional area of ice floes, the net radiation at the water surface controlled the radiation budget. Because the water albedo is much less than the albedo of the ice floes, the extent and duration of open water in summer is an important albedo feedback mechanism. From summer to freeze-up, the net all-sky shortwave radiation declined steadily as the solar angle lowered, while coincidently the net all-sky longwave radiation became increasingly negative. The all-sky net surface radiation switched from positive in summer to negative during the freeze-up period. From summer to freeze-up, both upward and downward turbulent heat fluxes occurred. In summer, a positive surface energy budget residual contributed to the melting of ice floes and/or to the warming of the Arctic Ocean's mixed layer. During the freeze-up period, with temperatures below approximately ?5°C, the residuals were mainly negative suggesting that heat loss from the ocean's mixed layer and heat released by the phase change of water were significant components of the energy budget's residual.