IMPLEMENTATION OF AN URBAN CANOPY PARAMETERIZATION FOR FINE-SCALE SIMULATIONS
Tanya L. Otte1'* and Avraham Lacser2
1 Atmospheric Sciences Modeling Division, NOAA ARL, Research Triangle Park, North Carolina
'Israel Institute for Biological Research, Ness Ziona, Israel
"On assignment to the National Exposure Research Laboratory, U.S. EPA
1.	INTRODUCTION
The Pennsylvania State University/National Center
(or Atmospheric Research Mesoscale Model (MM5)
(Grell et al. 1994) has been modified to include an urban
canopy parameterization (UCP) tor fine-scale urban
simulations (-1-km horizontal grid spacing). The UCP
accounts for drag exerted by urban structures, thf
enhancement of turbulent kinetic energy (TKE)
especially near the tops of the buildings, and the
modification of the energy budget within the urban
canopy (i.e., from the surface to the tops of buildings).
This UCP is applied to grid cells in MM5 that have a
non-zero fraction of urban land use. This refinement o1
MM5 is targeted to enable the Community Multiscale Air
Quality (CMAQ) Modeling System (Byun and Ching
1999) to capture the details of pollutant spatial
distributions in urban areas, Preliminary results will be
presented below.
2.	URBAN CANOPY PARAMETERIZATION
2.1 Momentum
The horizontal components of the momentum
equations are modified to account for the area average
effect of the sub-grid urban elements following Brown
(2000). The modifications are implemented in MM5 via
the TKE-based Gayno-Seaman planetary bounoary
layer (PBL) parameterization scheme (e.g., Shatrari
et al. 2000). The momentum equations accounting tor
the urban elements are:
dU	/ 2
= Fa - 0.5furbCd A z) U (U + y2)
9 V	! 2 2\0'5
— = Fy~Q.5furbCtfj4(z)V(lT +V )
0j "" =f'TKF +0.5TurbCdMz)|^S + + W2)
where the Fs are the peneral forcing terms in each
equation; U, V, and W are the wind components; and
TKE is the turbulent kinetic energy. In this lormulation.
it is assumed that the builoings affeel the flow but oo noi
take up any volume within the grid cell. Ca is a drag
coefficient. The urban fraction of the grid cell it
described by /urt. A(z) is the canopy area density, oi
the surface area of the obstacle (e.g., building;
perpendicular to the wind, per unit volume of the urban
canopy, expressed in rrf1. There are several
approaches to Describe A.i) (e.g., Uno et al. 198t
Brown 2000) where the function reaches its maximum at
the ground level, but vanishes at the top of thf
obstacles (so the drag term vanishes also at that level).
The integral of fi(x) from the ground level to the lops of
the tallest buildings (H) is Kf, which corresponds to the
ratio of the frontal area to the total surface area of the
buildings. In general /(z) is a function of the location
within the domain as it depends on building morphology.
A.I) can be estimated from and H assuming some
functional form for here we use a linear function.
To solve the modified momentum equations (the
added new term), we follow the analytical solution
suggested by Byun and Arya (1886). The TKE equation
is solved explicitly. To take proper account of the
influence of A(z), the vertical resolution in MM5 is
increased in the domain such that several prognostic
layers are below H.
2.2	Energy Budget
To account for the impact of urban settings on the
energy budget, the anthropogenic heat flux is included
in the heat equation and not in the surface energy
budget (e.g., Chin et al. 2000). The anthropogenic heat
flux is set as a function of urban land use subcategory,
and it has a temporal weighting function (e.g., Taha
1999), The heat equation also includes the heat
contribution of the city canyons following Yamada
(1982); the contribution due to rooftops has not yet been
implemented. The surface energy balance includes the
shadowing/trapping effect of the net radiation reaching
the ground in the city canyons modified by extinction of
the radiation through the urban canopy using a simple
exponential function (e.g., Brown 2000).
2.3	Urban Morphology
Parameters that are required for the UCP (e.g., H
and Xf) can be extracted from digital imagery (e.g.,
Rafti et al. 2001) which is commercially available from
several vendors tor various cities and with different
degrees of accuracy and precision. We did not hove
access to a true urban morphology database for our
area of interest, so we modified the MM5 land use
database for our domain to have seven subcategories of
urban areas adapted from R. E let sen (personal
communication 2001). These categories crudely
represent urban zones such as high-rise, industrial, and
urban residential. These urban zones are used to
create pseudo-morphology lor the region of interest, and
the physical characteristics of these areas art used in
the UCP. Each of the urban subcategories has a
different value for H. kf, canyon fraction, and maximum
anthropogenic heal fiux.
"Corresponding author address; Tanya Otte, U.S. EPA.
MD E243-03, RTF, NC 27711; otte.tanya@epa.gov

-------
3. PRELIMINARY RESULTS
The unmodified MM5 was run in a one-way nested
mode for several days in July 1995 during which there
was a high-ozone episode in "the Northeastern U.S. that
coincided with a photochemical field study. The MMt
domains included a five-domain configuration (108, 36,
12, 4, and 1.33 km grid spacing). The first four domains'
were run with 30 vertical layers (about 12 layers in the
PBL, and lowest level at 19 m) and physics options
appropriate for each resolution. To include the influence
of smaller obstacles, the UCP was used on the 1.33-km
domain (Figure 1} with 40 layers that included ten new
layers in the lowest 100 m (lowest level at 2 m).
Several simulations have been made with the
1.33-km domain to determine the impact of the UCF
and the morphology on the simulation. Four 1.33-km
simulations are discussed here. The first experiment
nocanSO, is a 30-layer experiment without the UCf!
The second experiment, nocan40, is a 40-iavei
experiment without the UCF. Both noean30 and
nocan40 are shown so that the impact of the UCP or,
the simulation can be evaluated independent of the
increase in vertical resolution. Expt. morph40 includes
the UCP with urban morphology. Expt. homog40
includes the UCF and a homogeneous representation ot
the city defined by the ' average" of the morphology.
Figure 2 is a comparison of the surface (~2-m;
temperature for Philadelphia International Airport (PHL)
trom the four 1.33-nm experiments compared to
observations lor 14 July 1995. Figure 2 shows that the
two experiments with the UCP improved the simulated
nighttime temperatures by ~2°C and improved the
maximum temperatures by -0.5'C. In the mid-
afternoon, the addition of the ten new layers improver
the verification, as the 40-layer experiments werf
superior to nocanSO Overall, the diurnal temperature
patterns of the simulations with the UCP compare morf
favorably with observations than the two experiment,'
without the UCF
As expected tor a non-urban location, Millville. New
Jersey (MIV), the use of the UCP had little effect on tht
diurnal pattern of the temperature (not shown). The lou-
simulations are very similar lor nighttime temperature.'
In addition, the same mid-aflernoon impact ot the ten
new layers occurs si MIV as well as PHL. Thi:
suggests that there are no detrimental effects of thf
UCP on areas beyond the urban core.
Verification statistics were calculated against thf
National Weather Service standard observation site:
shown in Figure 1. Root-mean-square error (RMSE)
statistics reflect the nighttime improvements in iht
temperature in the urban areas in the UCP simulations:
Over the 24-h simulation period, the two experiment:
with the UCF outperformed the two experiments
without. Expt. homoo40 has an RMSE of 1.9'C, and tht
RMSE in morpMO is 2.0'C. Expt. nocar>30 has ar
RMSE of 2.2*"C, while nocan40 has an RMSE of 2.4'C
The index of agreement (e.g., Willmott, 1982) for thf
same period reflects the same ordinal ranking of thj
experiments and a laroe separation between the UCF
and no-UCP simulations. These statistics tend to show
that the use of the UCP generally improves thf
temperature simulations throughout the domain
compared to the "standard" 30-layer configuration, and
the improvements are not solely due to the increase in
vertical resolution from 30 to 40 layers.
The hourly RMSE comparison for moisture (not
shown) indicates thai the 40-iayer experiments are
significantly better than nocan30. However, morph40
verified best of the four experiments (2.11 g/kg), while
tiomog40 verified third (2.21 g/kg).
The wind speed timeseries for various sites in the
domain do not show a clear "winner", potentially due to
the more erratic hourly changes in wind speed.
However, there is clearly an effect on the wind speed
timeseries both in the urban areas (e.g., PNE, not
shown) and downstream (e.g., WRI, not shown), ' The
RMSE over the 24-h simulation (not shown) indicates
that morph40 is the best of the experiments (1.32 m s"').
Expt. homog4Q was the worst of the four experiments
11.56 m s } for wind speed. The bias (mode! minus
observation) for the same time period (Figure 3) clearly
shows that morph40 is the best of the experiments
f-0,31 m s"1), followed by nocanSO (+0.75 ms"1),
nocari40 (-0,91 m s"'), and homog40 (-0,96 m s'1). This
is consistent with the RMSE scores discussed above.
The UCP had very little effect on wind direction in the
timeseries and the statistical scores over the same
simulation period, Thus the UCP with morphology has a
positive impact on the simulations of winds at 1.33-km.
4. CLOSING THOUGHTS
Initial tests with the UCP for simulations on 14 July
1995 show that the UCP at 1.33-km tends to produce
the desired effects on the wind, temperature, and TKE
fields. The temperature simulations with the UCP tend
to be superior to 30-layer and 40-layer simulations
without the UCP. In addition, the UCP simulations with
pseudo-morphology tend to be superior to the UCP
simulations with an "average' homogeneous
morphology (i.e., no urban subcategories), particularly
for wind speed. Further evaluation ot the UCP in MM5
on additional days is underway
Future plans include expanding the UCP to add the
root energy budget and to use more detailed land use
databases and morphology databases so that the urban
a.'Las can be more accurately characterized. Also, the
UCP will be coupled with a land-surlace model and
urban soil model.
Disclaimer. The information in this manuscript has been
piepared under funding by tht- United States
Environmental Protection Agency. It has been
subjected to Agency review and approved for
publication. Mention of trade names or commercial
P'oducts does not constitute endorsement or
recommendation for use.
4. REFERENCES
Biown, M, J., 2000: Urban paramelerizations for
mesoscaie meteorological models Mesoscate
Atmospheric Dispersion. Ed.. 2. Boybeyi.
byun, D. W., and S. P. S. An,a. 1986: A study of mixed
layer momentum evolution. Atmos. Environ 20
715-728.

-------
Byun, D. W., and J. K, S. Ching, 1999; Science
algorithms of the EPA Mccels-3 Community
Multiscale Air Quality (CMAQ) Modeling System,
EPA-600/R-99/030, U.S. EPA.
Chin, H. S., M. L Leach, and M J. Brown, 2000: A
sensitivity study of the urban effect on a regional
scale model: an idealized case. Preprints, Third
Sym. on the Urban Environment, Davis, CA, Arner.
Meteor. Soc.
Grell, G„ J. Dudhia, and D R. Stauffer, 1994: A
description of the fifth-generation Penn State/NCAR
Mesoscale Model (MM5). NCAR Tech. Note,
NCAR/TN398+STR, 138 pp.
Ratti, C., S. Di Sabatino, R. E. Britter, M. J. Brown, F.
Caton and S. Burian, 2001: Analysis of 3-D urban
databases with respect to pollution dispersion for &
number of European and American cities..
Proceedings, Third Int. Conf. on Urban Air Quality,
Loutraki, Greece.
Shafran, P. €.. N. L Seaman, and G. A. Gayno, 2000:
Evaluation of numerical predictions of boundary
layer structure during the Lake Michigan Ozone
Study. J. Appl. Meteor., 38, 412-426
Taha, H , 19S9: Modifying a mesoscale meteorological
mode! to better incorporate urban heat storage: A
bulk parameterization approach. J. Appl. Meteor.,
38, 466-473.
Uno, I., H. Ueda, and S. Wakamalsu, 1989: Numerical
modeling of the nocturnal urban boundary layer,
Bound-Layer Meteor., 49, 77-98.
Willmott, C. J., 1S82: Some comments on the evaluation
of model performance. Bull Amer. Meteor. Soc.,
63, 1309-1313.
Yamada, T., 1982: A numerical model study of turbulent
airflow in and above a forest canopy. J. Met. Soc.
Japan, 60 (1), 439—454.
TTN
NEL
PNI
WRi
PHL
ACY
MIV,
Figure 1. The stancard National Weather Service-
observation sites used tor verification on the 1.33-km
domain.
2
0
•2
•4
0
24
time (h)
Figure 2. Temperature timeseries (model minus
observation) of four 1.33-km experiments compared to
observations for Philadelphia international Airport (PHL)
tor 14 July 1995. Time is in UTC.
0.5

-0.5

¦1.5
time (h)
hpure 3. Wind speed bias for tour experiments verified
acains! ihe observation sites shown in Figure 1.
Verification is for 14 July 1995 1 ime is in UTC

-------
TECHNICAL REPORT DATA
1. REPORT NO.
2.
3 .
4. TITLE AND SUBTITLE
Implementation of an urban canopy parameterization
for fine-scale simulations
5.REPORT DATE
6.PERFORMING ORGANIZATION CODE
7. AUTHORCS)
T.L. Otte, and A. Lacser
8.PERFORMING ORGANIZATION REPORT NO.
9. PERFORMING ORGANIZATION NAME AND ADDRESS
Same as Block 12
10.PROGRAM ELEMENT NO.
11, CONTRACT/GRANT NO.
12. SPONSORING AGENCY NAME AND ADDRESS
U.S. Environmental Protection Agency
Office of Research and Development
National Exposure Research Laboratory
Research Triangle Park, NC 27711
13.TYPE OF REPORT AND PERIOD COVERED
14. SPONSORING AGENCY CODE
EPA/600/9
15. SUPPLEMENTARY NOTES
16. ABSTRACT
The Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model
(MMS) {Grell et al. 1994) has been modified to include an urban canopy parameterization (UCP)
for fine-scale urban simulations {~1-km horizontal grid spacing). The UCP accounts for drag
exerted by urban structures, the enhancement of turbulent kinetic energy (TKE) especially
near the tops of the buildings, and the modification of the energy budget within the urban
canopy (i.e., from the surface to the tops of buildings). This UCP is applied to grid cells
in MM5 that have an non-zero fraction of urban land use. This refinement of MM5 is targeted
to enable the Community Multiscale Air Quality (CMAQ) Modeling System (Byun and Ching 1999)
to capture the details of pollutant spatial distributions in urban areas. Preliminary results
will be presented below.
17. KEY WORDS AND DOCUMENT ANALYSIS
a. DESCRIPTORS
b.IDENTIFIERS/ OPEN ENDED
TERMS
c.COSATI



18. DISTRIBUTION STATEMENT
RELEASE TO PUBLIC
19. SECURITY CLASS (This
Report)
UNCLASSIFIED
21,NO, CH PAGES
20. SECURITY CLASS (This
Page)
UNCLASSIFIED
22, PRICE

-------