2011 International Conference on Alternative Energy in Developing Countries and Emerging Economies
- 324 -
for sea surface temperature [7]. The model was run for
seven days plus one day for spin-up, repeated over an
entire year to create an annual series of hourly output.
Land cover data sets from the Global Land Cover 2000
Database [8] and the Land Development Department of
Thailand (LDD) were combined and converted to 24
USGS categories used in default MM5 inputs. The
amount of urban land surface area in the final domain
was under-represented in the default data, with just 5% of
grid cells classified as urban compared to 45% in the
modified data. The single urban category was also split
into three subcategories to account for varying building
morphology across the city. These subcategories (Fig. 2)
are identical in all aspects except aerodynamic surface
roughness, with 0.5 m for low-density and height urban
areas, 1.5 m for medium- to high-density and height
urban areas, and 3.0 m for high-rise areas. These values
were chosen from the undertaking of visual surveys,
examination of satellite imagery, and comparison with
aerodynamic surface roughness from other studies [9,
10].
Fig. 2. Urban subcategories designed for
this study, and location of the PCD Tower.
Default values for fraction of vegetation cover over the
land surface in MM5 are of low resolution and thus vary
little over the final domain. These were therefore
replaced with data originating from MODIS satellite-
derived normalized difference vegetation index data,
converted to vegetation fraction and averaged to give
monthly values, to more accurately represent vegetation
fraction.
A MODIS satellite-derived albedo product was also
incorporated, but due to a high proportion of missing data
in rainy months, tabulated rather than gridded values
were used. Seasonal averages were calculated for each
land cover category and inserted into the land surface
parameter table in MM5. Modified values over
categories containing cropland, pasture, grassland and
shrubland are lower and have less seasonal variation than
default MM5 values. This is to be expected as Bangkok
is a tropical area with a significant percentage of land
area covered with dense vegetation that varies relatively
little over the seasons. In forested and urban areas, which
both have similar vegetation cover worldwide, there was
little difference in albedo.
To test the response of the model to land surface
inputs, the model was run over two weeks in each of
January and July, both with default land surface inputs
and with the modified inputs described above. Hourly
near-surface temperature and wind speed at 100 m agl
were extracted from model output to the location of the
Bangkok Pollution Control Department (PCD)
meteorological tower (location shown in Fig. 2). These
were compared with hourly data from observations
acquired from the PCD tower, and mean bias
MB
calculated for each using
N
OP
MB
N
i
i
i
¦
)
(
,
(1)
where
P
i
and
O
i
are the predicted (modeled) and observed
values respectively at hour
i
, and
N
is the total number of
hourly prediction-observation pairs. Comparison of
default and modified inputs in Fig. 3 and Fig. 4 shows
improved temperature results with modified inputs,
particularly in the dry season (January). Although mean
bias for wind speed has little change, we used the
modified inputs for the annual simulation as they better
represent the land surface of Bangkok.
Fig. 3. Mean bias of near-surface temperature at the PCD
tower, for default and modified land surface model inputs.
Fig. 4. Mean bias of wind speed at 100 m agl at the PCD
tower, for default and modified land surface inputs.