UBC Faculty Research and Publications

Evaluation of the Town Energy Balance (TEB) scheme with direct measurements from dry districts in two… Masson, V.; Grimmond, C. S. B.; Oke, Timothy R. 2002

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata


52383-Oke_AMS_2002JAM1011.pdf [ 557.91kB ]
JSON: 52383-1.0041955.json
JSON-LD: 52383-1.0041955-ld.json
RDF/XML (Pretty): 52383-1.0041955-rdf.xml
RDF/JSON: 52383-1.0041955-rdf.json
Turtle: 52383-1.0041955-turtle.txt
N-Triples: 52383-1.0041955-rdf-ntriples.txt
Original Record: 52383-1.0041955-source.json
Full Text

Full Text

OCTOBER 2002  MASSON ET AL.  1011  Evaluation of the Town Energy Balance (TEB) Scheme with Direct Measurements from Dry Districts in Two Cities V. MASSON Centre National de Recherches Me´te´orologiques, Toulouse, France  C. S. B. GRIMMOND Atmospheric Science Program, Department of Geography, Indiana University, Bloomington, Indiana  T. R. OKE Department of Geography, University of British Columbia, Vancouver, British Columbia, Canada (Manuscript received 31 October 2001, in final form 9 May 2002) ABSTRACT The Town Energy Balance (TEB) model of Masson simulates turbulent fluxes for urban areas. It is forced with atmospheric data and radiation recorded above roof level and incorporates detailed representations of the urban surface (canyon geometry) to simulate energy balances for walls, roads, and roofs. Here the authors evaluate TEB using directly measured surface temperatures and local-scale energy balance and radiation fluxes for two ‘‘simple’’ urban sites: a downtown area within the historic core of Mexico City, Mexico (stone buildings five to six stories in height), and a light industrial site in Vancouver, British Columbia, Canada (flat-roofed, single-story warehouses). At both sites, vegetation cover is less than 5%, which permits direct evaluation of TEB in the absence of a coupled vegetation scheme. Following small modifications to TEB, notably to the aerodynamic resistance formulations, the model is shown to perform well overall. In Mexico City, with deep urban canyons and stone walls, almost two-thirds of the net radiation is partitioned into storage heat flux during the day, and this maintains large heat releases and an upward turbulent sensible heat flux at night. TEB simulates all of these features well. At both sites TEB correctly simulates the net radiation, surface temperatures, and the partitioning between the turbulent and storage heat fluxes. The composite wall temperature simulated by TEB is close to the average of the four measured wall temperatures. A sensitivity analysis of model parameters shows TEB is fairly robust; for the conditions considered here, TEB is most sensitive to roof characteristics and incoming solar radiation.  1. Introduction Although a general understanding of the interactions between the atmosphere and urbanized areas based on carefully conducted experimental studies is emerging (see, e.g., Rotach 1995; Grimmond and Oke 1999a), the complexity and diversity of cities around the world often means observational studies are limited, either to a particular site or a single physical process. Numerical models have been developed to simulate the urban climate, but to date most of these models have been formulated either at the microscale (building and/or urban canyon) (e.g., Mills 1997; Arnfield and Grimmond 1998; Arnfield et al. 1998), or at the mesoscale (whole city), using vegetation–atmosphere transfer models originally developed for other surface types, but with parameters modified for application to the urban environment (e.g., Corresponding author address: Dr. Vale´ry Masson, Centre National de Recherches Me´te´orologiques, 42 av Coriolis, 31057 Toulouse Cedex, France. E-mail: valery.masson@meteo.fr  ᭧ 2002 American Meteorological Society  Best 1998; Taha 1999). In an attempt to couple the two scales, Masson (2000) developed the Town Energy Balance (TEB) scheme, which can be used on its own for densely urbanized areas or with the Interactions between Soil, Biosphere, and Atmosphere (ISBA) model (Noilhan and Planton 1989) when vegetation is also present. TEB incorporates canyon geometry, with three typical surfaces—roof, wall, and road—in order to reproduce the effects produced by buildings. A specific energy balance is computed for each of these three surfaces. This approach is relatively simple, but it still allows most of the physical effects associated with the urban energy balance to be reproduced, including radiative trapping (longwave and shortwave), the momentum flux, the turbulent sensible and latent heat fluxes, heat storage uptake and release, and water and snow interception. The anthropogenic heat fluxes due to traffic and industry are prescribed, and they are implicitly computed for domestic heating by the use of a minimum internal building temperature. Masson (2000) validated the radiative portions of  1012  JOURNAL OF APPLIED METEOROLOGY  VOLUME 41  FIG. 1. Coupling of TEB within an atmospheric model with relation to measured and modeled variables. Dashed lines are the limits of the atmospheric grid boxes. Dotted lines are the middle of each grid box. For definition of the box see text. Here ⌬Q A is the advection flux below the measurement height.  TEB; however, only sensitivity experiments were conducted for the complete scheme. Here the objective is to present an independent evaluation of TEB using directly measured surface temperatures and surface energy balance fluxes for two ‘‘simple’’ urban sites. One is a portion of the historic core of Mexico City, Mexico, the other is a light industrial site in Vancouver, British Columbia, Canada. At both sites vegetation cover is less than 5% of the plan area (Grimmond and Oke 2002, their Fig. 5) but the sites are otherwise very different in terms of the structure and building construction. These tests allow direct evaluation of TEB and the appropriateness of the assigned parameter values and their sensitivity/robustness. Evaluation of the combined TEB–ISBA scheme, relevant to vegetated urban areas, is being pursued separately. TEB has been incorporated into MesoNH, the French community mesoscale model (Lafore et al. 1998), but here we consider the offline or stand-alone version. 2. The urban energy balance Fundamental to this paper is the concept of the urban surface energy balance, defined by (Oke 1988) Q* ϩ Q F ϭ Q H ϩ Q E ϩ ⌬Q S ϩ ⌬Q A [W mϪ2 ], (1) with Q* being the net all-wave radiation, Q F the an-  thropogenic heat flux, Q H the turbulent sensible heat flux, Q E the turbulent latent heat flux, ⌬Q S the storage heat flux, and ⌬Q A the net advective heat flux. In numerical models, such as TEB, each of the surface energy balance fluxes in (1) can be addressed independently. In this study, because the sites chosen have very little vegetation cover and there was no precipitation or irrigation, Q E is set to zero. The TEB–ISBA scheme, when operated in the offline mode, does not explicitly resolve the ⌬Q A term in (1). The model is forced with a temperature (T a ), humidity (q a ), and wind speed (U a ) measured in the inertial sublayer (Fig. 1). For TEB online, these are at the middle of the first atmospheric level, and for TEB offline they are at the measurement height. The rooftop level is the surface of the atmospheric model, so the TEB computed sensible heat flux (Q HTEB ) for the urban canopy is assigned at the base of the first atmospheric grid box (Figs. 1 and 2). The modeled output fluxes (Q H , Q E , outgoing longwave radiation L↑ and outgoing shortwave radiation K↑, when coupled to a 3D flow model, are assigned as input at the surface of the atmospheric model (the base of level 1 in Fig. 1), while forcing variables correspond to the midpoint of the level 1 grid. Note that the surface of the atmospheric model is located at the top of the mean topography and the building (roughness element) height for the individual grid square. In the 3D case,  OCTOBER 2002  MASSON ET AL.  the fluxes are then used to calculate the values of T a , q a , and U a at the next time step. These interact with the surrounding grids, so the ⌬Q A term can be resolved. The impact on the surface fluxes is returned by the forcing (T a , q a , and U a ) driving the next time step of the TEB–ISBA surface scheme. If ⌬Q A is positive (usually because of energy supplied by the surface), this energy comes from the Q H term convected from the surfaces, so Q H at a higher level (typically the middle of level 1) is smaller by an amount ⌬Q A . The converse holds for a negative advection flux (Q H increasing with height). In this manner the vertical profile of Q H is affected by this advection flux. When the energy balance is measured, however, (1) must be modified to read Q* ϭ Q H ϩ Q E ϩ ⌬Q S [W m Ϫ2 ].  (2)  The terms in (2) are evaluated here by measurements at the top of a ‘‘box,’’ following the concept of Oke (1987, 1988) illustrated in Fig. 1. The height of the box extends from above the roughness sublayer to a depth in the ground where there is no vertical flux over the period of concern. The height of the box is typically about two times the height of the buildings. The horizontal dimension of the box is sufficient to ensure that microscale inhomogeneity merges into a representative local-scale property. The height and length of the box must be sufficient to ensure that microscale variability of fluxes in the source area (or ‘‘footprint’’) is eliminated by turbulent mixing in the roughness sublayer, but not so great that the source area extends upwind into anomalous land cover. Thus, Q HOBS is located in the inertial layer, at a height sufficient to avoid microscale effects found in the roughness sublayer, but not so high that anomalous upwind surfaces are sampled. The ⌬Q A and Q F do not appear explicitly in (2). If the land cover in the source area of the sensors is approximately homogeneous at the local scale, ⌬Q A is minimized and can be neglected. This can also be achieved at less than completely homogeneous sites through filtering the data by wind direction to eliminate sectors with anomalous surface cover. Heat released by combustion of fuels (Q F ) is a net source of energy in cities; however, the instruments used to measure Q*, Q H , and Q E sense such anthropogenic contributions. Hence the total measured fluxes include this term, although the exact partitioning is not known. Adding an additional Q F term to (2) would therefore double-count this source. Here three components (Q*, Q H , and Q E ) are measured directly and ⌬Q S is determined as the residual of (2). In effect this means that storage changes in the volume of the box are expressed as an equivalent flux density through the top surface of the box. This has the effect that the net measurement errors incurred in estimating the other terms, including ⌬Q A , accumulate in ⌬Q S . It should be noted that closure of measured energy balances is rarely achieved, either over simple sites  1013  where micrometeorological theory is most likely to hold (e.g., Foken and Oncley 1995) or more complex sites such as forests (e.g., Lee 1998). This contributes to flux uncertainty, especially of ⌬Q S , when resolved as a residual. Given the differences between the modeled and measured energy balances [i.e., (1) and (2)], exact correspondence between measured and modeled fluxes is not to be expected. In particular, comparisons of the turbulent fluxes are affected by the assumption of zero evaporation by TEB (because ISBA is not run). This means TEB should overestimate Q H by the value of measured Q E . Here, Q H calculated by TEB is compared at the measurement height with the sum of measured Q H ϩ Q E . As discussed above, offline TEB Q H values at this height include ⌬Q A . Thus from a modeling perspective the grid cell is regarded as homogeneous, but from the measurement perspective the area is homogeneous only if the site has sufficiently extensive fetch to minimize local-scale ⌬Q A (Fig. 1). Therefore, in the evaluation of the model in the offline version, it is necessary to assume that the turbulent fluxes Q H and Q E from top of the urban canopy layer (where TEB fluxes are computed; Fig. 1) to the top of the box (where the measurements are made) remain the same. 3. Observations Here we present a brief overview of the observation methods used to gather data appropriate to the evaluation of TEB; more complete details are available in Grimmond and Oke (1999b), Oke et al. (1999), and Voogt and Grimmond (2000). a. The study sites Data were collected for Mexico City (referred to here as Me93) for a 7-day period during the dry season in December 1993 (Oke et al. 1999). The study site is located in the historic core of the city (see Fig. 1 of Oke et al. 1999) where the buildings are mainly institutional and commercial. In the vicinity of the measurement location, the mean building height is 18.4 Ϯ 6.6 m. The average complete-to-plan area ratio (or the three-dimensional aspect ratio; Voogt and Oke 1997), ␭ c , is 1.75. Average (3D) surface cover for the dominant source areas for the flux measurements is 25% impervious, 32% roofs, and 42% walls (Oke et al. 1999). Vegetation is negligible (ϳ1%). In the street canyons and alleys adjacent to the tower, the sky view factor varies from about 0.22 to 0.51. Based on an inventory of materials within the radius of 500 m around the Me93 tower site, walls are made of concrete (two thirds) of stone (one third), with roofs of concrete, tar, sheet metal, or tile. Building walls and roofs are usually grey or brown. The roads are paved with concrete or asphalt, or surfaced with tiles, cobblestones, or flagstones. Data were also collected for 15 days, in August 1992,  1014  JOURNAL OF APPLIED METEOROLOGY  at a light industrial site in Vancouver, British Columbia (hereinafter referred to as V192), during an extended period of drought. The area is characterized by flatroofed buildings one to three stories in height (average height 6.9 Ϯ 2.5 m, ␭ c , is 1.39), used for light industry and as warehouses. The buildings are arranged in city blocks, with east–west and north–south orientation (see Fig. 3 of Voogt and Oke 1997). The axis of the main blocks and alleys run east–west between the main streets. Less than 5% of the plan area is covered by vegetation. Most buildings are made of concrete, and roads and pathways are made of asphalt and concrete. b. Measured surface energy balance fluxes and forcing data At each site, instruments were mounted high enough above the surface to ensure that the measurement box (Fig. 1) is representative of the integrated local scale (horizontal length scales of ϳ10 2–10 4 m). In Mexico City and Vancouver, equipment was mounted at 28.4 and 28.5 m above ground, respectively. In both cases the turbulent fluxes were measured using eddy covariance techniques and radiometry was used to determine Q*. At the V192 site, ground heat flux plates were installed at a depth of 80 mm with temperature sensors above to account for flux divergence between the soil surface and the heat flux plates. In addition, standard meteorological data, needed as input for TEB, were collected at both sites. The primary meteorological variable that forces TEB, incoming solar radiation (K↓), was not measured at either site. However, K↓ data were collected by the Institute of Geophysics, Universidad Nacional Auto´noma de Me´xico at the university (ϳ13 km south of the Me93 site), and by the Canadian Atmospheric Environment Service at the Vancouver International Airport (ϳ8 km south of V192 site). Using previous studies on the spatial variation of K↓, and patterns of air pollution and flow in the two cities, we assessed the degree to which the measured solar forcing might differ from that at the validation site. We conclude this may amount to a few percent too low (typically less than 5%) in the afternoon at Me93, but it is likely to be negligible (less than 1%) at V192. Incoming longwave radiation (L↓) was not measured at either site but was calculated using observed air temperature and humidity following the method of Prata (1996). In Mexico City, the surface temperatures of a roof and two roads in canyons adjacent to the site were measured by infrared thermometers. In Vancouver, the surface temperatures of four walls (one facing each cardinal direction), the roof, and the inside of the east-facing wall of a building at the site were measured using similar instruments. These observations do not provide a representative sample of all facets found within the localscale turbulent source area (typically 15–40 ϫ 10 3 m 2 by day, larger at night) at either site. Rather, they rep-  VOLUME 41  resent fairly typical microscale surfaces to provide some guidance as to the thermal conditions of small areas (sensor field of view about 5–10 m 2 ). In addition, airborne infrared measurements were conducted at the Vancouver site during one day at three different times: morning, mid-, and late afternoon. These measurements give the surface temperatures of each type of facet (walls, roofs, sunlit and shaded roads) at the local scale, an area more comparable with that of the source area of the turbulent fluxes. Therefore, even though only one day of data is available, it is of great interest to the validation of TEB. Note that the airborne measured wall temperatures compare well with the wall measurements made from a vehicle (Voogt and Oke 1997). 4. Implementation of TEB Small modifications were made to the scheme originally presented in Masson (2000), notably to the aerodynamic resistance formulations (those indicated with thick lines in Fig. 2). • The road resistance here follows a roughness length formulation, instead of that of Rowley et al. (1930), although the latter is still used for walls. This modification is included because of analyses of the V192 site, where roads are important and where airborne measured road temperatures allow discrimination between the two formulations (see section 5b). The road roughness length is set to 50 mm. This is much larger than the commonly used roughness for asphalt (Ͻ1 mm). However, 50 mm is considered to be a reasonable value when obstacles of the order of 1 m are present, as is the case for roads with cars, pedestrians, signs, lamp posts, bollards, etc. Sample calculations using a morphometric formula (Grimmond and Oke 1999a) show the roughness length can even be larger than the values used here if numerous cars and trucks are present. • The temperatures of internal roof and wall layers are now linked to the internal temperature of the buildings (T ibld ) using a standard resistance (R i ) equal to 0.123 K m W Ϫ1 , based on building insulation values (McMullan 1992). In Masson (2000), T ibld was directly applied to the inner surfaces, whereas here R i mimics both the convection and the radiative interactions in the building. This modification has negligible effect, except for poorly insulated buildings, such as warehouses with metal roofs. In such cases the temperatures of the wall and roof inner layers now have greater diurnal amplitude, leading to larger amplitude for the outside surfaces. • A new temperature evolution equation is used to represent the annual cycle of the internal temperature inside the buildings (especially for warm temperatures in summer):  OCTOBER 2002  1015  MASSON ET AL.  FIG. 2. Schematic representation of the surfaces (roof, wall, road indicated by subscript R, w, and r, respectively), prognostic temperatures (T ), and aerodynamic resistances (R) used in TEB and the output fluxes. Resistances shown with thick lines have been altered from the original scheme (Masson 2000). Thus, Q HTEB ϭ f roof Q Hroof ϩ f road Q Htop ϩ Q Findustry and Q ETEB ϭ f roof Q Eroof ϩ f road Q Etop ϩ Q Eindustry .  T ϩi bld ϭ T Ϫi bld  ΂  ΃  ΂ ΃  ␶ Ϫ ⌬␶ ⌬t ϩ T* , ␶ ␶  where T ϩibld and are T Ϫibld the temperatures at the future and previous time step, respectively, ⌬t is the time step, ␶ is equal to 1 day, and T* is the average of the internal (ceiling, wall, and floor) surface temperatures. However, as in the initial version of the scheme, T ϩibld has a minimum value, in order to represent anthropogenic heating. The scheme is initialized with the meteorological data. No anthropogenic heat flux is prescribed. The parameters listed in Tables 1 and 2 are used to initialize the scheme in Me93 and V192. In both cases the morphometric parameters are weighted according to the wind direction frequency to be as consistent as possible with the observed fluxes. The roof, wall, and road fractions are calculated from Grimmond and Oke (1999b). The road is estimated to be composed of concrete (40%) and asphalt (60%) pavement (50 mm thick), over stone aggregate (0.2 m) and gravel and sandy clay soil. In general, the thermal parameters and emissivities are esti-  mated from data listed in the American Society of Heating, Refrigerating, and Air-Conditioning Engineers handbook (ASHRAE 1989). The albedo values are based on those listed for roofs, walls, and roads (Oke 1987). At V192 the building walls are made of uninsulated concrete. The warehouse, where the microscale surface temperatures were observed, had a roof with a steel layer overlaid by a thin gravel layer. However, as most of the other warehouse roofs within this vicinity are made of concrete, topped with a thin insulation layer and gravel, this configuration is used in TEB (Table 2). 5. Evaluation of TEB a. Mexico City 1) SURFACE  TEMPERATURES AT THE MICROSCALE  Since TEB computes only one road surface temperature, independent of direction, the TEB value is compared to the average of two road surface temperature measurements (one in an S–N, one in an E–W street). The general form of the diurnal evolution of roof and  1016  JOURNAL OF APPLIED METEOROLOGY  VOLUME 41  TABLE 1. Input parameters for the TEB scheme for the Vancouver light industrial (Vl92) and Mexico City historic core (Me93) sites. Parameters  Unit  Mexico City  Vancouver  Geometric parameters Building fractionb Building heightb Wall/plane area ratiob ␭C Canyon aspect ratiob : H/W Roughness lengthb  — m — — m  0.55 18.8 0.75 1.18 2.2  0.51 5.8 0.39 0.39 0.35  Radiative parametersc Roof albedo Wall albedo Road albedo Roof emissivity Wall emissivity Road emissivity  — — — — — —  0.20 0.25 0.08 0.90 0.85 0.95  0.12 0.50 0.08 0.92 0.90 0.95  a  Thermal parametersc (see Table 2 for values) Roof  Asphalt roll on concrete, good insulation Massive, thick, stone or concrete, insulation, 25% window Asphalt 60%, concrete 40% on gravel, sandy clay 0.15 0.05  Wall Road Roof roughness length Road roughness length  m m  Temperature initialization Inside building temperature Deep soil temperature  ЊC ЊC  Thin gravel over concrete, poor insulation Concrete without insulation Asphalt 60%, concrete 40% on gravel, sandy clay 0.15 0.05  20 22  23 20  Local-scale parameters. Source: Grimmond and Oke (1999a). c Microscale parameters of individual facets. a  b  road temperatures are correctly reproduced by TEB (Fig. 3, Table 3). However, TEB roof temperatures show larger diurnal amplitudes than those of the instrumented roof, with faster warming and cooling rates. This may be due to larger storage in the roof by TEB or to incorrect specification of the roof materials (see section  6 for a sensitivity test). A way to correct the daytime overestimation might be to increase the roughness length of the roof and hence the removal of heat via Q H . It must also be appreciated that the observations are only for one roof. The TEB road temperature has both the correct temporal course and diurnal amplitude,  TABLE 2. Thermal parameters for roofs, walls, and roads used in TEB for the Vancouver and Mexico City sites. Layer sequence: 1 is nearest to the surface. Here d is thickness of layer (m), C is heat capacity of layer (MJ mϪ3 KϪ1 ), and ␭ is thermal conductivity (W mϪ1 KϪ1 ). Mexico City 1 Roof layer Asphalt roll d C ␭  0.01 1.7 0.2  Vancouver  2  3  4  1  2  3  4  Concrete (stone)  Insulation  Gypsum  Gravel  Gravel  Insulation  Concrete  0.1 1.5 0.93  0.05 0.25 0.03  0.025 0.87 0.16  0.01 1.76 1.4  0.02 1.76 1.4  0.01 0.04 0.03  0.030 2.21 1.51  Concrete  Dense concrete  Wall layer Stone/window Stone/window Stone/window  Insulation/ window  Dense concrete Dense concrete  d C ␭  0.015 1.54 0.88  0.12 1.54 0.88  0.30 1.54 0.88  0.015 0.32 0.21  0.010 2.11 1.51  0.02 2.11 1.51  0.14 1.00 0.67  0.03 2.11 1.51  Road layer  Asphalt/ concrete  Asphalt/ concrete  Stone aggregate  Gravel and soil  Asphalt/ concrete  Asphalt/ concrete  Stone aggregate  Gravel and soil  d C ␭  0.01 1.74 0.82  0.04 1.74 0.82  0.20 2.00 2.1  1.00 1.40 0.4  0.01 1.74 0.82  0.04 1.74 0.82  0.20 2.00 2.1  1.00 1.40 0.4  OCTOBER 2002  1017  MASSON ET AL.  TABLE 3. Performance statistics for surface temperatures in Mexico City and a Vancouver light industrial site. Bias ϭ TEB Ϫ obs. Bias (K)  Rmse (K)  ϩ4 ϩ1.9  4.2 4.5  YD 223–236 (all) Wall temperature Roof temperature  ϩ2.3 ϩ2.5  3.0 7.4  YD 225–231 (period 1) Wall temperature Roof temperature  ϩ1.9 ϩ1.7  2.4 7.5  YD 232–236 (period 2) Wall temperature Roof temperature  ϩ2.7 ϩ3.6  3.7 7.2  Mexico City Average road temperature Roof temperature Vancouver  FIG. 3. Ensemble mean diurnal cycle for 6 days of measured and modeled surface temperatures for the Mexico City site: (top) roof surface temperature. and (bottom) road surface temperature.  but there is a bias of about ϩ4 K. Given the uncertainties in using a calculated L↓ (especially at night), the large difference between the observed temperatures in the two canyons (25 K during the day), and the microscale (rather than local scale) representativeness of the measurements, this is considered a good result. 2) ENERGY  FLUXES AT THE LOCAL SCALE  TEB is evaluated using the measured hourly energy fluxes: net radiation, sensible heat flux, and storage heat flux (Fig. 4 gives the ensemble mean results and Table 4 the summary statistics). During the night, the observed Q* mean is Ϫ103 W  m Ϫ2 , while TEB is only Ϫ82 W m Ϫ2 . The daytime Q* is well reproduced (Table 4). The observed daytime ⌬Q S /Q* ratio is very high (0.58); hence uptake of sensible heat into storage is more efficient than its convection to the boundary layer over this densely builtup district. This can probably be attributed to the thermal properties of the construction materials, the large 3D surface area exposed by this area of dense building and deep canyons, and the relatively light winds during the observation period (the mid- to late afternoon maximum speeds were Ͻ3.5 m s Ϫ1 ). TEB reproduces this result well; the simulated ⌬Q S /Q* ratio is 0.56 (Table 4). The observations show little hysteresis in the diurnal relation between the ⌬Q S and the Q*, and again TEB reproduces this well (Fig. 5). Similarly, TEB produces ⌬Q S fluxes of the correct magnitude, especially for the daytime maximum. The nocturnal measurements show the heat release from the urban fabric to be larger than the Q* drain (nighttime ⌬Q S /Q* ϭ 1.21). This imbalance maintains an average upward-directed Q H of 21 W m Ϫ2 . The model produces a nighttime ⌬Q S /Q* ratio of 1.17, but since the calculated Q* intensity is smaller than observed, the release of stored heat is also underestimated by 30 W m Ϫ2 . Therefore, the modeled Q H reaches only 12 W m Ϫ2 , on average. However, it is encouraging that TEB is able to create a positive Q H every night, including the cloudy night, as observed (not shown). The observed rapid increase of Q H (reaching between 50 and 100 W m Ϫ2 ) in the early morning is not simulated by the model. The observed increase in Q H might correspond to energy released from traffic (Q F ) during the morning commuter activity that is not included in these TEB runs. Ichinose et al. (1999), for example, report a traffic-induced heat flux in downtown Tokyo reaching 60 W m Ϫ2 . Thus in summary, TEB is able to reproduce correctly most of the behavior of the fluxes typical of the Me93 central city site, including the relatively low Q* in the  1018  JOURNAL OF APPLIED METEOROLOGY  VOLUME 41  TABLE 4. Performance statistics of TEB for heat fluxes (W m Ϫ2 ) at the Mexico City site. Average [observed (obs), modeled (TEB)], bias (TEB Ϫ obs), and rmse.  Overall period  Daytime  Nighttime  Heat flux  Q*  QH ϩ QE  ⌬Q S  Obs TEB Bias Rmse Obs TEB Bias Rmse Obs TEB Bias Rmse  45 55 10 32 257 252 Ϫ5 41 Ϫ103 Ϫ82 21 24  57 54 Ϫ3 25 108 113 5 32 21 12 Ϫ9 18  Ϫ12 1 13 39 149 139 Ϫ10 45 Ϫ125 Ϫ95 30 35  winter, the large daytime heat uptake by the urban fabric, and the positive Q H at night. b. Vancouver light industrial site 1) SURFACE  TEMPERATURES AT THE MICRO- AND LOCAL SCALES  The surface temperature observations include in situ observations of individual roof, wall, and road facets of a warehouse, that is, at the microscale; and airborne measurements at 1000, 1400, and 1700 LST on one day [yearday (YD) 228] of a large number of buildings and canyons, that is, at the local scale. The latter is most appropriate to use in the evaluation of TEB. For the rest of the observations our discussion focuses on two periods: period 1, six sunny days (YD 225–231) characterized by a sea breeze in the early afternoon; and period 2, five more cloudy days (YD 232–236) when the sea breeze is much weaker and occurs only in the late afternoon. The observed microscale roof temperatures show a very strong diurnal cycle, with variations of 50 K (Fig.  FIG. 4. Ensemble mean diurnal cycle of measured and modeled surface energy balance fluxes for 6 days at the Mexico City site: (top) net all-wave radiation, (middle) total turbulent flux, and (bottom) storage sensible heat flux. FIG. 5. Ensemble mean (6 days) measured and modeled (TEB) hysteresis loop relating the storage heat flux to the net radiation at the Mexico City site.  OCTOBER 2002  MASSON ET AL.  1019  FIG. 6. Ensemble mean diurnal cycle of observed (microscale) and modeled (TEB) surface temperature for the Vancouver site for (a) YD 225–231 and (b) YD 232–236: (top) roof and (bottom) individual N-, E-, S-, and W-facing walls and average, and simulated wall temperature.  6). TEB simulates a smaller amplitude than is observed (the reverse of the Me93 case). TEB does reproduce the general form of the nocturnal cooling pattern, but observed roof temperatures are low at night and the model overestimates them by 8 K, on average. The roof where  the surface temperature observations were conducted is unrepresentative of the local scale (see section 4), thus the poor comparisons are not unexpected. TEB is in somewhat better agreement with the airborne roof data of the morning and evening flights but underestimates  1020  JOURNAL OF APPLIED METEOROLOGY  VOLUME 41  FIG. 7. Time series of observed (micro- and local scale) and modeled (TEB) surface temperatures for YD 228: (left) roof, (middle) averaged microscale and local-scale wall temperature, and (right) local-scale road temperature. Simulated road surface temperature: Masson (2000) scheme [road resistance follows Rowley et al. (1930)] and new roughness length formulation.  the early afternoon value, which is near the daily maximum, by approximately 6 K (Fig. 7). The TEB wall temperature is evaluated against the average of the observed wall temperatures (north, east, south, and west facing). For sunny days in period 1, TEB wall temperatures are similar to this averaged wall temperature (Table 3). The bias and root-mean-square error (rmse) between TEB and the observations are small (1.9 and 2.4 K, respectively). Note that the difference is slightly larger when compared to the airborne measurements; the overestimation probably reaches 5 K in the daytime (Fig. 7), but the maximum wall temperature near 1600 LST is reproduced. However, a discrepancy appears for period 2. On the sunny days of period 1, the individual wall temperatures follow a standard daytime pattern and the mean wall temperature peaks at about 1600 LST (Fig. 6a). In period 2, on the other hand, there is significant cloud cover in the middle of the day. While three of the walls show little change from period 1, the temperature of the south-facing wall was considerably cooler (Fig. 6b). TEB does not capture the magnitude of this reduced warming (except for the cloudiest day, not shown) and model performance is correspondingly poorer. It seems the temporal history of warming of the individual walls over the day, and its disruption by cloud, is responsible. Such temporal specificity is not incorporated in the generic wall approach of TEB. If this is correct, this type of impact can be expected to be greatest in cases of urban geometry with low building height-to-street width ratios, as is the case at the V192 site. Overall, however, the good agreement between the wall temperature measurement and TEB output, especially when radiative forcing is symmetric through the day, supports one of the main simplifications of the  scheme: the use of only one wall temperature that is independent of orientation. For roads, the airborne-based measurements can be stratified into shaded and sunlit fractions. The sunlit fraction is much larger, because of the small height of the buildings (Voogt and Grimmond 2000). The road temperature, averaged according to sunlit and shaded areas, reaches 38ЊC in the afternoon of YD 228. TEB simulates only one road temperature, which is comparable to the local-scale observation. TEB simulates it well, with an overestimation of only 3–4 K at midday (Fig. 7). The simplification of TEB to consider only one road temperature is also supported by this comparison. To illustrate the improvement resulting from the new parameterizations (see section 4), TEB was run with the old formulation of the road aerodynamic resistance (Fig. 7). The main problem with the original TEB predictions is an increase of the temperature during the day, 6 K higher than with the new version. The new roughness length formulation creates a larger conductance, and thus a larger Q H , leading to better agreement between the modeled temperature and the observations. Since the road surface at the V192 site is large, and has an important role for the exchanges from the canyons (walls are low), this validation is particularly pertinent. 2) ENERGY  FLUXES AT THE LOCAL SCALE  TEB correctly reproduces Q* during the entire period (Fig. 8). As Table 5 shows this is true for the overall set (bias ϭ Ϫ9 W m Ϫ2 ) and for both day and night (biases of Ϫ17 and ϩ2 W m Ϫ2 , respectively). When comparing modeled and measured data it is important to consider the potential for measurement errors and uncertainties. Here, the observed fluxes and  OCTOBER 2002  MASSON ET AL.  1021  FIG. 8. Ensemble mean diurnal cycle of measured and modeled surface energy balance fluxes for the Vancouver site, for (a) YD 225–231 and (b) YD 232–236: (top) net radiation, (middle) total turbulent heat flux, and (bottom) storage sensible heat flux.  1022  JOURNAL OF APPLIED METEOROLOGY  VOLUME 41  TABLE 5. Performance statistics for TEB for heat fluxes at the Vancouver light industrial site (as in Table 4). Days 223–236  Overall period  Daytime  Nighttime  Obs TEB Bias Rmse Obs TEB Bias Rmse Obs TEB Bias Rmse  Days 225–231  Days 232–236  Q*  QH ϩ QE  ⌬Q S  Q*  QH ϩ QE  ⌬Q S  Q*  QH ϩ QE  ⌬Q S  150 141 Ϫ9 59 323 306 Ϫ17 76 Ϫ59 Ϫ57 2 24  95 133 38 76 168 234 66 103 7 11 4 12  55 8 Ϫ47 91 156 73 Ϫ83 121 Ϫ66 Ϫ68 Ϫ2 23  161 159 Ϫ2 57 318 313 Ϫ5 72 Ϫ62 Ϫ59 3 20  87 147 60 97 143 243 100 127 7 11 4 11  74 12 Ϫ62 105 175 70 Ϫ105 136 Ϫ69 Ϫ70 Ϫ1 19  136 126 Ϫ10 59 323 301 Ϫ21 77 Ϫ56 Ϫ55 1 27  103 119 16 50 196 225 29 69 8 11 3 13  32 6 Ϫ26 66 126 76 Ϫ50 89 Ϫ64 Ϫ66 Ϫ2 26  meteorological conditions at the V192 site were compared with concurrent observations made at a residential site (sunset, Vs92; Grimmond and Oke 1999b). In urban areas, Q* is expected to be fairly conservative spatially (Schmid et al. 1991; Oke 1997). Comparison of the data collected concurrently at the two Vancouver sites shows only small differences in Q*; the largest differences Ͻ50 W m Ϫ2 occur in the late morning when Q* at V192 is larger. This indicates reason for confidence in the V192 Q* data. During night, the measurements show relatively large ⌬Q S , with a peak release just after sunset. This allows a positive Q H to be maintained during the first hours of the night. Then ⌬Q S equilibrates close to the value of Q*, and Q H is small for the rest of the night. In terms of the observed storage heat flux, it is important to recall that it is computed here as the residual of the measured energy balance and thus errors in any of the other fluxes accumulate in ⌬Q S . However, during night, fluxes other than Q* and ⌬Q S are small (see the measured Q H and Q E ) or negligible (Q F or ⌬Q A ). The warehouses in the V192 area have minimal heating or air conditioning, few chimneys, and traffic is light in the day and negligible at night. The ⌬Q A is negligible, because nocturnal winds are weak except on one night. Therefore,  the observed nocturnal ⌬Q S is probably reliable. TEB simulates this nocturnal storage pattern well (bias ϭ Ϫ2 W m Ϫ2 , rmse Ͻ 23 W m Ϫ2 throughout the period; see Table 5). The net energy storage of the model during YD 226–231 is 0.4 MJ m Ϫ2 , corresponding to a mean storage heat flux of about ϩ0.8 W m Ϫ2 , or a 3 K heating of a 0.1-m layer of soil or concrete. This buildup of heat is consistent with results from heat flux plates installed at the site (not shown). Note that hourly values are not necessarily good, but the modeled integrated daily heat uptake is about right. During the day, observed Q H is large and Q E is small. In the drought conditions of the measurement period, TEB partitions Q* between ⌬Q S and Q H . For period 2 (Fig. 8b) and YD 224, 225 (not shown), there is reasonable correspondence between the predictions of TEB and the measured sensible heat fluxes (Table 5). The hysteresis behavior of ⌬Q S (Fig. 9) is slightly overestimated. For period 1, however, TEB significantly overpredicts Q H , while measurements show a much larger ⌬Q S than TEB (Fig. 8a). This discrepancy may be due to an overestimation of observed ⌬Q S [in (2)], through the neglect of ⌬Q A . During period 1, the wind direction is predominantly from the west so there may be a sea breeze (Oke and Hay 1994). Under these circumstances it is possible that the lack of ability to evaluate the advection term in both the measurements and the model results prevents a more appropriate comparison of the fluxes. As noted in section 2, the advection flux in the box ⌬Q A appears in the Q H term at the top of the box when modeled offline and in the ⌬Q S term in the observations. An improved estimate would come from running the model in the 3D mode. 6. Sensitivity of the scheme  FIG. 9. Ensemble mean measured and modeled (TEB) hysteresis loop relating the storage heat flux to the net radiation at the Vancouver light industrial site for YD 232–236.  The ability of the TEB scheme to reproduce most features of surface–atmosphere energy exchanges and the resultant surface temperatures of relatively dry urban districts has been demonstrated. However, many model input properties and some of the forcing data (e.g., radiation) are subject to uncertainties. Therefore, the sen-  OCTOBER 2002  1023  MASSON ET AL.  TABLE 6. Sensitivity analysis to varying TEB input parameters relative to the reference case for Mexico City. Max Q* (W mϪ2 )  Daytime ⌬Q S /Q*  Nocturnal Nocturnal Q* QH (W mϪ2 ) (W mϪ2 )  Reference  439  0.56  Ϫ82  Geometric parameters Building fraction ϩ 0.10 (→ H/W ϭ 1.42) Building fraction Ϫ 0.10 (→ H/W ϭ 1.03) Wall/plane area ratio ϩ 0.20 (→ H/W ϭ 1.40) Wall/plane area ratio Ϫ 0.20 (→ H/W ϭ 0.96)  ϩ1 ϩ1 ϩ6 Ϫ3  Ϫ0.01 ϩ0.01 ϩ0.01 Ϫ0.01  ϩ1 Ϫ1 Ϫ2 ϩ2  ϩ3 Ϫ3  Ϫ32  ϩ0.03  ϩ2  Ϫ1  Radiative parameters Roof albedo ϩ 0.10 Wall albedo ϩ 0.10 Road albedo ϩ 0.10 Roof emissivity Ϫ 0.05 Wall emissivity Ϫ 0.05 Road emissivity Ϫ 0.05 Thermal parameters Roof thickness ϫ 2 Wall thickness Ϭ by 2 Road thickness Ϭ by 2 Roughness lengths Town z 0 Ϭ by 2 Town z 0 ϫ 2 Road z 0 Ϭ by 5 Roof z 0 Ϭ by 5 Temperature initialization Internal building temperature ϩ 5 K Road and deep soil temperature ϩ 5 K Forcing Incoming infrared Ϫ 20 W mϪ2 Incoming solar Ϫ 10%  Ϫ6 Ϫ6 ϩ4  ϩ1  ϩ1  Ϫ8  Ϫ0.04  ϩ4  Ϫ2  ϩ2  ϩ0.01  Ϫ1  Ϫ1 ϩ1  Ϫ2 ϩ3 Ϫ1 Ϫ11  ϩ0.02 Ϫ0.02 ϩ0.01 ϩ0.05  Ϫ1 ϩ1 Ϫ2  Ϫ2 ϩ1  Ϫ1 Ϫ1  Ϫ0.01 Ϫ0.01  Ϫ1 Ϫ1  ϩ1 ϩ2  Ϫ12 Ϫ49  ϩ0.04 ϩ0.03  Ϫ10 ϩ2  Ϫ5 Ϫ3  Modification of resistance Without new road resistance  sitivity of the scheme to input parameters was tested via dedicated runs across a wide range of values for each parameter. Tables 6 and 7 list the modifications and the resulting sensitivity of the calculated fluxes and temperatures. Comparisons are made against a reference run. Daytime Q* is relatively insensitive to most parameters except for the albedo of horizontal surfaces, the roof roughness length, and K↓. Geometric parameters modify Q* to a lesser extent (canyon aspect ratios vary in the range Ϯ0.2 from the reference). The effect of roofs is always the most important, because there is no trapping effect for this surface. Because of the geometry of the site, the road albedo has a very small impact for Me93. Nighttime Q* also is not very sensitive, except to the calculated L↓, the roof thickness or its thermal structure (for V192), and the building fraction, if the  Twall bias (K)  Troad bias (K)  ϩ14  Ϫ1 Ϫ0.01  Troof bias (K)  Ϫ0.1  Ϫ1.3 (␴ ϭ 1) ϩ0.5  Ϫ0.4 (␴ ϭ 2.2)  Ϫ0.4 Ϫ0.2 Ϫ0.1  Ϫ0.1  ϩ0.1 Ϫ0.1 ϩ0.1 Ϫ0.1 (␴ ϭ 0.7)  ϩ0.1 Ϫ0.4 ϩ0.1  ϩ0.2 (␴ ϭ 1.2)  ϩ0.4 Ϫ0.6  ϩ0.5 Ϫ0.6 ϩ0.5  ϩ0.2  ϩ0.1 ϩ0.2  ϩ0.1 ϩ0.6  Ϫ1.9 Ϫ1.1 (␴ ϭ 0.8)  Ϫ0.7 Ϫ0.6  Ϫ1.1 Ϫ0.7  ϩ1.5 (␴ ϭ 1.2)  ϩ0.1  road and roof energy balances are very different (as is the case in V192). The ⌬Q S /Q* ratios for these sensitivity analyses, are in the range 0.26–0.31 and 0.52–0.61 for V192 and Me93, respectively. For most parameters, variations in input values result in differences of less than 2% in ⌬Q S / Q* between runs. Because only dry sites are involved this also means the scheme is stable with respect to the partitioning between the storage and turbulent sensible heat fluxes. For the V192 simulation, ⌬Q S /Q* is sensitive to roof roughness length (because it is a major control on the convective/conductive sensible heat sharing between the air and the roof ). A large impact is found for the choice of the road aerodynamic resistance. The road surface temperature measurements at V192 enabled us to discriminate between the two formulations (see section 5b). A too-small roughness length for roads  1024  JOURNAL OF APPLIED METEOROLOGY  VOLUME 41  TABLE 7. Sensitivity analysis to varying TEB input parameters relative to the reference case for Vancouver. Max Q* (W mϪ2 )  Daytime ⌬Q S /Q*  Nocturnal Nocturnal Q* QH (W mϪ2 ) (W mϪ2 )  Troof bias (K)  Twall bias (K)  Troad bias (K)  Ϫ0.1  Ϫ0.3 ϩ0.1 Ϫ0.4 ϩ0.4  Reference  520  0.28  Ϫ61  14  Geometric parameters Building fraction ϩ 0.10 (→ H/W ϭ 0.50) Building fraction Ϫ 0.10 (→ H/W ϭ 0.33) Wall/plane area ratio ϩ 0.20 (→ H/W ϭ 0.60) Wall/plane area ratio Ϫ 0.20 (→ H/W ϭ 0.19)  Ϫ5 ϩ3 ϩ1 Ϫ3  Ϫ0.02 ϩ0.02 ϩ0.02 Ϫ0.02  ϩ4 Ϫ4 Ϫ1 ϩ2  Ϫ2 ϩ1 ϩ2 Ϫ2  Ϫ34 Ϫ4 Ϫ21 ϩ4  ϩ0.01 ϩ1  Ϫ1 Ϫ1  ϩ1  ϩ1 ϩ3 ϩ0.8 (␴ ϭ 2.3) Ϫ2 ϩ4 Ϫ1 Ϫ0.7 (␴ ϭ 2.2)  Ϫ0.3 ϩ0.1  Ϫ0.1 ϩ0.2  ϩ0.6 Ϫ0.6  ϩ0.6 Ϫ0.6 ϩ1.9 (␴ ϭ 1.0)  ϩ0.1 ϩ0.1  ϩ0.1 ϩ0.3  Radiative parameters Roof albedo ϩ 0.10 Wall albedo ϩ 0.10 Road albedo ϩ 0.10 Roof emissivity Ϫ 0.05 Wall emissivity Ϫ 0.05 Road emissivity Ϫ 0.05 Thermal parameters Roof thickness ϫ 2 Wall thickness Ϭ by 2 Road thickness Ϭ by 2 Steel roof (see text) Roughness lengths Town z 0 Ϭ by 2 Town z 0 ϫ 2 Road z 0 Ϭ by 5 Roof z 0 Ϭ by 5 Temperature initialization Building internal temperature ϩ 5 K Surface and deep soil temperature ϩ 5 K Forcing Incoming infrared Ϫ 20 W mϪ2 Incoming solar Ϫ 10% Air temperature Ϫ 2 K Modification of resistance Without building room resistances Without new road resistance  ϩ2  ϩ0.1 Ϫ0.1 Ϫ1.0 (␴ ϭ 0.9) ϩ0.3  Ϫ0.6 Ϫ0.2 ϩ0.1  Ϫ0.8 ϩ0.2  ϩ2  ϩ0.02  ϩ2 Ϫ2  ϩ0.01 Ϫ0.02  Ϫ6 ϩ1 Ϫ2 ϩ5  Ϫ3 ϩ3 Ϫ8 Ϫ18  ϩ0.01 Ϫ0.01 ϩ0.02 ϩ0.03  Ϫ1 ϩ1 Ϫ2 Ϫ1  ϩ1 ϩ1  Ϫ1  Ϫ0.01  Ϫ1 Ϫ1  ϩ1  Ϫ15 Ϫ58 ϩ10  ϩ0.02 ϩ0.01 Ϫ0.02  Ϫ9 ϩ2 ϩ4  Ϫ4 Ϫ1.6 (␴ ϭ 0.9) Ϫ2 Ϫ0.8 (␴ ϭ 0.7) ϩ4 Ϫ0.9  Ϫ0.9 Ϫ0.9 Ϫ1.2  Ϫ1.0 Ϫ1.0 Ϫ1.3  ϩ2 Ϫ9  ϩ0.03 ϩ0.02  ϩ3 Ϫ2  Ϫ3 ϩ3  Ϫ1.0  Ϫ0.3 ϩ1.8 (␴ ϭ 1.4)  leads to an overestimate of the surface temperature during daytime (by 10 K, not shown), which in turn leads to too-large heat storage (Table 7). The impact of modifying the internal building resistance is significant, and it improves the roof temperature, especially at night (not shown). Some warehouses in the V192 study area have roofs made of gravel on steel, instead of gravel on steel and wood. Using such a roof in TEB (20 mm gravel, 0.50 mm steel, 10 mm wood) decreases the storage and gives a ⌬Q S /Q* ratio of 0.26. This dampens the thermal response of the roof surface (especially at night) compared to microscale measurements on this type of roof (not shown). This signals the need to take care in setting the value of the roof parameters for such cases. For the Me93 site, the roof specification again has a large impact, as does the radiative forcing. The road and wall temperatures simulated by TEB are not much affected by uncertainties in the radiative and thermal input parameters. In contrast, the roof temperature is sensitive, both on average (bias) and its diurnal evolution (standard deviation), particularly to the  ϩ1.9 (␴ ϭ 2.0) ϩ0.1  Ϫ0.5  albedo assigned. All temperatures are sensitive to incoming radiation and the roughness length. Summarizing the results of the sensitivity tests it seems that emissivities, temperature initializations, and wall characteristics have little impact. Similarly, the output is relatively insensitive to the overall roughness length for the city, even when large variations were considered. Road characteristics are not critical, except perhaps for albedo in areas with low buildings and wide roads. Errors in assigning geometric parameters have a moderate impact. Uncertainties in setting incoming radiative fluxes are significant. Thus we recommend that, if possible, measured values be used. Care must be taken not to use roughness values for horizontal surfaces that are too small and to take careful note of roof characteristics. However, it is our view that for the conditions tested, TEB is robust. To put this in perspective, across the range of sensitivity conditions tested here, the ⌬Q S /Q* ratio is always within 5% of the reference run.  OCTOBER 2002  1025  MASSON ET AL.  7. Concluding comments Evaluation of TEB using field observations from two relatively dry and sparsely vegetated urban areas, with very different urban structure, suggests that overall the model performs well. The evaluation process resulted in small modifications to the original TEB scheme; these were in the aerodynamic resistance formulations. The agreement between the measurements and TEB gives some support to one of the simplifications of the scheme: the use of only one wall temperature. The results of this evaluation are informative not only in assessing the ability of the model to calculate fluxes and thermal responses in these urban settings, but also in aiding interpretation of existing field data and to inform the better design of future observational programs. For example, the effect of the timing of clouds in sites with different canyon geometries, and the role of advection both in measurements and modeling need further investigation. Next we intend to extend the evaluation to include vegetated and moist urban districts using the combined TEB–ISBA scheme. This study is perhaps the first to attempt genuine validation of the output of an urban land surface model that is designed for use within a mesoscale atmospheric model. Several previous studies have provided comparisons with field data from standard climate stations, especially with regard to air temperature, but none has checked output against observed fluxes and temperatures together. Agreement against air temperature can be relatively easily achieved, but when agreement is between fluxes and surface temperatures, it shows that the energetic processes underlying the thermal estimates are also physically realistic. Significant difficulties remain with obtaining high quality measurements of surface fluxes. Energy balance closure is not achieved even in simple rural sites. In the urban case, uncertainties associated with anthropogenic releases of heat, water, and pollutants including radiatively active aerosols and gases from concentrated and sometimes organized source distributions, and microadvection that characterizes exchanges in among elements of the urban canopy layer contribute to uncertainty in fluxes, especially heat storage change, if it is resolved as a residual. Nevertheless, the estimates from several cities show similar features; this gives some confidence that anomalous effects do not dominate in the balance. Similarly, there are challenges for modeling. Not least is the need to strike a balance between the detail and complexity of the modeled phenomena on the one hand and the simplicity of input properties and computational efficiency on the other. Given certain conceptual and practical difficulties we also note the difficulty of getting model output and observations to apply to a common ‘‘surface’’ or atmospheric plane. The TEB model and the nature of the observational database used in the present study were not designed together. One of the issues that this research draws at-  tention to is the importance of scale and ensuring that the measurement and modeling communities are both cognizant of this. Better and more comprehensive observational data will emerge and both our understanding of the physics of the urban atmosphere and computational capacity will grow, but in our judgment the present comparison of observed and modeled climatic conditions represents a reasonable convergence of presentday capabilities of measurement and simulation. Acknowledgments. This work has been supported by grants to the authors by the Centre National de Recherches Me´te´orologiques, the Centre National de Recherches Scientifiques de France, the National Science Foundation, the Natural Sciences and Engineering Research Council of Canada, and the Canadian Foundation for Climate and Atmospheric Science. The authors are very grateful to the following colleagues who supplied data or provided assistance with interpretations: Stephanie Meyn of the University of British Columbia; Dr. Rachel Spronken-Smith of the University of Canterbury; Dr. Agustı´n Muhlia Vela´zquez, Coordinador del Observatorio de Radiacio´n Solar, Instituto de Geofisica, UNAM; and Dr. James Voogt of the University of Western Ontario. REFERENCES Arnfield, A. J., and C. S. B. Grimmond, 1998: An urban canyon energy balance model and its application to urban storage heat flux modelling. Energy Build., 27, 61–68. ——, J. M. Herbert, and G. T. Johnson, 1998: A numerical simulation investigation of urban canyon energy budget variations. Preprints, Second Symp. on the Urban Environment, Albuquerque, NM, Amer. Meteor. Soc., 2–5. ASHRAE, 1989: ASHRAE Handbook: 1989 Fundamentals. ASHRAE, 37 sections ϩ errata ϩ index, 797 pp. Best, M., 1998: Representing urban areas in numerical weather prediction models. Preprints, Second Symp. on the Urban Environment, Albuquerque, NM, Amer. Meteor. Soc., 148–151. Foken, T., and S. Oncley, 1995: Results of the workshop ‘‘Instrumental and Methodical Problems of Land Surface Flux Measurements.’’ Bull. Amer. Meteor. Soc., 76, 1191–1193. Grimmond, C. S. B., and T. R. Oke, 1999a: Aerodynamic properties of urban areas derived from analysis of surface form. J. Appl. Meteor., 38, 1262–1292. ——, and ——, 1999b: Heat storage in urban areas: Local-scale observations and evaluation of a simple model. J. Appl. Meteor., 38, 922–940. ——, and ——, 2002: Turbulent heat fluxes in urban areas: Observations and local-scale urban meteorological parameterization scheme (LUMPS). J. Appl. Meteor., 41, 792–810. Ichinose, T., K. Shimodozono, and K. Hanaki, 1999: Impact of anthropogenic heat on urban climate in Tokyo. Atmos. Environ., 33, 3897–3909. Lafore, J., and Coauthors, 1998: The Meso-NH atmospheric simulation system. Part I: Adiabatic formulation and control simulation. Ann. Geophys., 16, 90–109. Lee, X., 1998: On micrometeorological observations of surface-air exchange over tall vegetation. Agric. Forest Meteor., 91, 39– 49. Masson, V., 2000: A physically-based scheme for the urban energy balance in atmospheric models. Bound.-Layer Meteor., 94, 357– 397.  1026  JOURNAL OF APPLIED METEOROLOGY  McMullan, R., 1992: Environmental Science in Building. 3d ed. Macmillan, 332 pp. Mills, G., 1997: An urban canopy-layer climate model. Theor. Appl. Climatol., 57, 229–244. Noilhan, J., and S. Planton, 1989: A simple parameterization of land surface processes for meteorological models. Mon. Wea. Rev., 117, 536–549. Oke, T. R., 1987: Boundary Layer Climates. 2d ed. Routledge, 435 pp. ——, 1988: The urban energy balance. Prog. Phys. Geogr., 12, 471– 508. ——, 1997: Urban environments. The Surface Climates of Canada, W. G. Bailey, T. R. Oke, and W. R. Rouse, Eds., McGill–Queen’s University Press, 303–327. ——, and J. E. Hay, 1994: The Climate of Vancouver. B.C. Geographical Series, No. 50, 84 pp. ——, R. A. Spronken-Smith, E. Jauregui, and C. S. B. Grimmond, 1999: The energy balance of central Mexico City during the dry season. Atmos. Environ., 33, 3919–3930.  VOLUME 41  Prata, A. J., 1996: A new long-wave formula for estimating downward clear-sky radiation at the surface. Quart. J. Roy. Meteor. Soc., 122, 1127–1151. Rotach, M. W., 1995: Profiles of turbulence statistics in and above an urban street canyon. Atmos. Environ., 29, 1473–1486. Rowley, F. B., A. B. Algren, and J. L. Blackshaw, 1930: Surface conductances as affected by air velocity, temperature and character of surface. ASHRAE Trans., 36, 429–446. Schmid, H. P., H. A. Cleugh, C. S. B. Grimmond, and T. R. Oke, 1991: Spatial variability of energy fluxes in suburban terrain. Bound.-Layer Meteor., 54, 249–276. Taha, H., 1999: Modifying a mesoscale meteorological model to better incorporate urban heat storage: A bulk parameterization approach. J. Appl. Meteor., 38, 466–473. Voogt, J. A., and T. R. Oke, 1997: Complete urban surface temperatures. J. Appl. Meteor., 36, 1117–1132. ——, and C. S. B. Grimmond, 2000: Modeling surface sensible heat flux using surface radiative temperatures in a simple urban area. J. Appl. Meteor., 39, 1679–1699.  


Citation Scheme:


Citations by CSL (citeproc-js)

Usage Statistics



Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items