New Algorithm For Down Scaling Temperature Values Chow2007BSERT
New Algorithm For Down Scaling Temperature Values Chow2007BSERT
New Algorithm For Down Scaling Temperature Values Chow2007BSERT
net/publication/245383326
CITATIONS READS
35 2,946
2 authors:
Some of the authors of this publication are also working on these related projects:
The development and assessment of the residential buildings energy code in the Kingdom of Saudi Arabia View project
All content following this page was uploaded by D.H.C. Chow on 03 March 2017.
New algorithm for generating hourly temperature values using daily maximum, minimum and
average values from climate models
D.H.C. Chow and Geoff J. Levermore
BUILDING SERV ENG RES TECHNOL 2007 28: 237
DOI: 10.1177/0143624407078642
Published by:
http://www.sagepublications.com
On behalf of:
Additional services and information for Building Services Engineering Research and Technology can be found at:
Subscriptions: http://bse.sagepub.com/subscriptions
Reprints: http://www.sagepub.com/journalsReprints.nav
Permissions: http://www.sagepub.com/journalsPermissions.nav
Citations: http://bse.sagepub.com/content/28/3/237.refs.html
The use of building simulation programmes for predicting building performance is increasing
all the time especially with the advent of cheap, fast computers. Hourly weather data,
in particular outdoor dry bulb temperature (DBT) and solar radiation values, are required for
simulation programmes for building performance. When hourly values are not available
there are algorithms for generating hourly temperature values from daily values. These use
the daily maximum temperature TMAX, and daily minimum temperature, TMIN. However,
climate prediction models, such as HadCM3 and HadRM3 also provide the daily
average dry-bulb temperature TAVE as well as the daily maximum and minimum.
The average temperature is important for selecting weather years and also because the
average temperature is often different from the average of the maximum and minimum,
assumed in the simpler algorithms. Buildings being designed now will need to perform
under future weather conditions with climate change, so the downscaling of daily values
from climate prediction models to hourly values is required. This paper describes a new,
more accurate algorithm for generating hourly temperature values in the UK that uses all
three temperature parameters from climate change models, and demonstrates the
improvement of the quality of the generated values against traditional algorithms that use
just the daily maxima and minima.
Practical application: The proposed algorithm for generating hourly DBT values from
daily maximum, minimum and average values is intended primarily for deriving hourly data
for running building simulation programmes, as some weather stations and future climate
prediction models only provide daily values of weather parameters.
Climate change is affecting all aspects of human life, and as well as being affected by
climate change, buildings can also affect the degree of climate change, since well-designed
buildings will require less energy to run, thus minimising the amount of carbon dioxide
emitted. As the climate is predicted to change significantly in the next 100 years, if buildings
are designed to last, it is important for designers to know how buildings will respond and
perform then. Building simulation programmes are useful for this, but they require hourly
weather data which are not provided by most climate prediction models. By having a
quality-assured algorithm for down-scaling the raw daily values to hourly values, data from
climate prediction models, such as HadCM3 and HadRM3 can be used for building
simulations for any location in the world.
The hourly DBTs derived from using the algorithm suggested in this paper may also be used
in conjunction with other hourly weather data, such as wet-bulb temperature, solar
irradiance, cloud cover and wind speed, derived from the same dataset.
3.00
2.00
1.00
Max pos. error
Error (°C)
−2.00
−3.00
1 2 3 4 5 6 7 8 9 10 11 12
Month
Figure 1 Monthly errors of predicting average daily temperature using the average of TMAX and TMIN for Heathrow 1976–1995
20.00
Temperature (°C)
16.00
14.00
12.00
Known Tmin value
10.00
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Hour
(See Table 1) and uses two sinusoidal curves where: For t5tmin:
to fit the data.
Hourly temperatures, T(t) are calculated cosððtmin tÞ=24 þ tmin tmax Þ þ 1
f1 ¼ : ð2Þ
using Equation (1). 2
For tmin5t5tmax:
TðtÞ ¼ f1 TMIN þ f2 TMAX ð1Þ
cosððt tmin Þ=tmax tmin Þ þ 1
f1 ¼ : ð3Þ
where f1 and f2 are factors given in a table in 2
the guide, and TMIN and TMAX are the daily For tmin5t:
minimum and maximum temperatures respec-
tively. These are related to sinusoidal inter- cosðð24 þ tmin tÞ=24 þ tmin tmax Þ þ 1
f1 ¼
polations, and can be expressed more 2
mathematically in Equations (2)–(5). ð4Þ
The shaded area in Figure 4, Shade 2 is given The temperature between and tmax is
in Equation (10). generated as shown in Equation (11):
Z
TMAX TMIN TMAX þ TMIN
Shade 2 ¼ TðtÞ ¼ TMIN
2 2
tmin
ðt tmax Þ
cos þ 1 TMIN
ðt 2 þ tmin Þ 2ðtmax Þ
cos þ1
2ð tmin Þ TMAX TMIN
¼
2
TMAX TMIN
¼ ðt tmax Þ
2 cos þ 1 þ TMIN : ð11Þ
2ðtmax Þ
2 The shaded area in Figure 5, Shade 3, is given
ð tmin Þ 1 : ð10Þ in Equation (12).
Z tmax
The situation from to tmax is shown in TMAX TMIN
Figure 5. Shade 3 ¼
2
ðt tmax Þ
cos þ1
2ðtmax Þ
TMAX TMIN
T MAX ¼ ðtmax Þ
2
T MAX + T MIN
2
2 1þ : ð12Þ
T MIN
Figure 6 shows the difference between the
half-sine wave and the quarter-sine wave
methods. Using the quarter-sine wave and
tmin l tmax
varying the position of allow the correct
Figure 4 Quarter-sine wave between tmin and TAVE to be achieved with the given TMAX and
TMIN.
T MAX
T MAX
Quarter-sines
T MAX + T MIN
T MAX + T MIN
2
Half-sine
2
T MIN
T MIN
tmin l tmax
tmin l tmax
Figure 6 Difference between half-sine wave and quarter-sine
Figure 5 Quarter-sine wave between and tmax waves
observed data. Absolute errors at each whole and Q-Sin algorithms for each individual year
degree for the cumulative distribution of the between 1976 and 1995. The performance of
generated data from real observed data are the new proposed Q-Sin algorithm is signifi-
summed and compared for each month. cantly better than the other algorithms for all
20 years, across all three analyses. It should
4.1.3 Comparing averaged hourly data also be noted that the algorithms where the
The third analysis of the algorithms days are linked together rather than treated
uses temperature at hour 1, hour 2, separately also performed better, having
hour 3, . . . hour 24 for each month (January, smaller errors.
February, March, . . . December). These are
compared with the real monthly averaged
values, and the sum of the errors from each 4.3 Effect of fitting average temperature
hour compared. Figures 10–12 show the performances of
the algorithms for each month, by averaging
4.2 Effect of linking days the errors for all 20 years for each month. The
Figures 7–9 show the average errors for the analysis of comparing cumulative plots
three tests for linked days and individual days (Figure 11) still shows the new proposed
1.80
CIBSE individual
1.70 days
Yearly average error (°C)
1.60
14R-1 individual
1.50 days
1.40 CIBSE linking days
1.30
1.20 14R-1 linking days
1.10
1.00 Q-Sin
1976
1977
1978
1979
1980
1981
1982
1983
1984
1985
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
Year
Figure 7 Yearly RMSE analysis of hourly temperature algorithms for Heathrow 1976–1995
35%
Yearly cumulatve percentage error
CIBSE individual
30% days
5%
Q-Sin
0%
1976
1977
1978
1979
1980
1981
1982
1983
1984
1985
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
Year
Figure 8 Yearly cumulative error analysis of hourly temperature algorithms for Heathrow 1976–1995
Figure 9 Yearly average hourly error analysis of hourly temperature algorithms for Heathrow 1976–1995
1.00
0.95 CIBSE individual
days
Monthly average error (°C)
0.90
0.85 14R-1 individual
days
0.80
0.75 CIBSE linking days
0.70
0.65 14R-1 linking days
0.60
0.55 Q-Sin
0.50
1 2 3 4 5 6 7 8 9 10 11 12
Month
Figure 10 Monthly RMSE analysis of hourly temperature algorithms for Heathrow 1976–1995
20.00%
Monthly cumulative percentage error
Figure 11 Monthly cumulative error analysis of hourly temperature algorithms for Heathrow 1976–1995
14R-1 individual
0.40 days
0.10
Q-Sin
0.00
1 2 3 4 5 6 7 8 9 10 11 12
Month
Figure 12 Monthly average hourly error analysis of hourly temperature algorithms for Heathrow 1976–1995
2.10
CIBSE individual days
1.90
Monthly average error (°C)
1.50
CIBSE linking days
1.30
14R-1 linking days
1.10
0.70
Q-Sin (CIBSE TmaxR-1)
0.50
1 2 3 4 5 6 7 8 9 10 11 12
Month
Figure 13 Revised monthly RMSE analysis of hourly temperature algorithms for Heathrow 1976–1995
algorithm having significantly fewer errors. changed to the CIBSE values for the Q-Sin
However, Figures 10 and 12 suggest that for algorithm, then results may improve.
two summer months (June and July), the For May, June, July and August, the times
CIBSE linked and individual days algorithm for TMAX were changed from hour 14 to 15,
performed better. This could be due to the 15, 16 and 15 respectively. The results are
difference in the designated times for the daily shown in Figures 13–15. This change appears
maximum temperature, TMAX. to give an over-riding success to the new
Figures 10 and 12 show that the 14R-1 proposed Q-Sin algorithm for downscaling
linking days algorithm generally perform daily temperature values to hourly values.
better than the CIBSE linking days algorithm, Figures 16 shows the maximum, minimum
with the exception of May, June, July and and average error range of hourly tempera-
August. The reason could be due to the ture for each month using the Q-Sin algo-
designation of TMAX times. If these were rithm. On average, the generated
35.00%
CIBSE individual days
Monthly cumulative percentage error
30.00%
14R-1 individual days
25.00%
10.00%
Q-Sin (14R-1)
5.00%
Q-Sin (CIBSE TmaxR-1)
0.00%
1 2 3 4 5 6 7 8 9 10 11 12
Month
Figure 14 Revised monthly cumulative error analysis of hourly temperature algorithms for Heathrow 1976–1995
1.20
CIBSE individual days
Monthly average hourly error (°C)
1.00
14R-1 individual days
0.80
CIBSE linking days
0.60
14R-1 linking days
0.40
Q-Sin (14R-1)
0.20
Q-Sin (CIBSE TmaxR-1)
0.00
1 2 3 4 5 6 7 8 9 10 11 12
Month
Figure 15 Revised monthly average hourly error analysis of hourly temperature algorithms for Heathrow 1976–1995
1.5
1.0
Difference in (°C)
0.5 Max
Min
0.0 Average
Average abs err +
−0.5 Average abs err −
−1.0
−1.5
−2.0
1 2 3 4 5 6 7 8 9 10 11 12
Month
Figure 16 Maximum, minimum and average error of generated hourly temperatures using Q-Sin algorithm for Heathrow 1976–1995
those described by the CIBSE Guide in require methods of downscaling from daily
Table 1. values to hourly values. Current and pro-
These revised times are not specific to the posed algorithms for these are to be covered
Q-Sin algorithm, and can also be applied to in a future paper as this paper concentrates on
improve tradition algorithms where only daily the more complex algorithm of downscaling
maximum and minimum temperature values daily DBT values to hourly values.
are available.
The use of this new algorithm is useful
particularly for generating TRYs and DSYs
for future time-slices, as most climate predic- References
tion models do not provide hourly values for
temperature, but the daily maximum, mini- 1 CIBSE Guide J Weather Solar and Illuminance
mum and average. Since the selection of Data, CIBSE, London, March 2002.
2 Hulme M et al. Climate change scenarios for
TRYs and DSYs use daily average tempera-
the United Kingdom: The UKCIP02 Scientific
tures, it would be useful to keep this variable Report, Tyndall Centre for Climate Change
when generating the resulting hourly values. Research, School of Environmental Sciences,
Since the improvement in the results is University of East Anglia, Norwich, 2002.
significant, it would provide a more accurate 3 CIBSE Guide A2 Weather & Solar Data, 1982.
simulation of the effect of climate change on 4 ASHRAE Handbook of Fundamentals,
buildings. Chapter 26, page 6, ASHRAE, Atlanta, 1993.
Building simulation programmes also 5 Sharmer K, Greif J. The European Solar
require hourly values of other weather para- Radiation Atlas Vol 2: Database and
meters such as wet-bulb temperature, solar Exploitation Software 2000; Paris, Les Presses
irradiance and wind speed. These will also de l’Ecole des Mines.