© Borgis - Balneologia Polska 3-4/2005, s. 106-110
Grzegorz Durło
Wyznaczanie dziennej wartości wielkości ochładzającej powietrza
Calculation of the mean diurnal bioclimatic cooling power
Department of Forest Protection and Forest Climatology
Agricultural University in Cracow
Streszczenie
Praca przedstawia metodę obliczania średniego dziennego ochładzania bioklimatycznego na podstawie standardowych pomiarów meteorologicznych. Opracowanie wykonano w oparciu dane, zgromadzone w okresie 1995-2004 na stacjach meteorologicznych i posterunkach mikroklimatycznych, położonych w południowej części kraju. Wartość dziennego ochładzania bioklimatycznego wyznaczono na podstawie średniej ważonej z trzech terminów obserwacyjnych, w której wagą były wyniki pomiarów z II terminu. Przeciętna różnica pomiędzy wartościami, średniej dziennej (13 godzin, od 7:00 do 19:00) i średniej ważonej, ochładzania bioklimatycznego suchego, dla serii rocznych, wynosiła 0,4%. Między wartością średniej z II terminu i średnią dzienną 1,8%, między II terminem a średnią ważoną 2,2%. Te same różnice dla serii miesięcznych wynosiły odpowiednio: 1,4%, 7,6% oraz 6,2%. Dla serii pięciodniowych jeszcze zwiększyły się i wynosiły odpowiednio 5,0 %, 18,5%, i 19,3%. Wartość dziennego ochładzania bioklimatycznego wyznaczonego w oparciu o średnią ważoną z trzech terminów obserwacyjnych jest najbardziej zbliżona do średniej rzeczywistej wartości wskaźnika, wyznaczonej w oparciu o dane cogodzinne.
Polecane
książki z księgarni medycznej BORGIS:
Introduction
The magnitude of cooling down of air is one of the physical indexes of heat exchange between the human body and the environment. It is very useful for the assessment of a thermal sensibility of man being in motion and dressed according to season of the year and actual atmospheric conditions. The magnitude of cooling down of air is an important element for the typological elaborations, and one of the main criteria for the estimation of weather conditions for the needs of bioclimatology (Kozłowska-Szczęsna et all 1997).
In most of bioclimatic papers the authors use the magnitude of cooling down of air calculated for the observation time-limit II, motivating that it represents that part of the day in which man works most actively, takes climatotherapeutic treatments, or rests. It has turned out, however, that at this time of the day there is dinner being served in many sanatoriums, convalescent hospitals, and rest-houses in health resorts of Poland. Moreover the period between 12:00 and 15:00 is considered to be the most unfavorable due to a high level of solar radiation and a high hazard of an excessive sunburn. This is especially important in open areas such as beaches, meadows, recreational forest clearings, and areas surrounding water reservoirs. It should also be noticed that human activity is not limited to only one form of spending time out of doors. This is the case when work takes place in the morning and forenoon, while rest and recreation in the afternoon and evening. As a rule, with the exception of some climatotherapeutic treatments especially arranged for a definite group of patients with a defective thermoregulation system, man stays out of doors from sunrise to sunset. Therefore, would the magnitude of cooling down of air determined on the basis of the observation time-limit II only, be adequate to characterize the conditions of thermal sensibility of a human being during the whole day? It seems that it would not be adequate. Besides, if we have values of meteorological elements from three observation time-limits at our disposal why not to use them?
The purpose of this study was to prove that to calculate the mean diurnal magnitude of bioclimatic cooling down it is necessary to use all values of magnitude of air cooling down, or meteorological elements used to its calculation measured in three observation time-limits.
It was assumed in the research hypothesis that using the formula for calculation of the weighted arithmetic mean, with the magnitude of bioclimatic cooling down from the measuring time-limit II as the weight factor, is the most correct method of calculation of the diurnal magnitude of air cooling down H in stations measuring meteorological elements in three observation time-limits.
Methods
The series of meteorological data gathered at meteorological stations and microclimatic posts situated in southern Poland (Beskid Śląski Mts, Beskid Sądecki Mts, Beskid Niski Mts, and Kraków-Częstochowa Upland) were used in this study (Tab. 1). The meteorological measurements were carried out by the HUGER automatic meteorological station with data recording every 60 minutes. At the meteorological station in Wisła Wyrchczadeczka the measurements were taken every 30 minutes obtaining 2880 results in November 2004 (Table 12). Together, there were gathered over 160 000 data on hourly values of air temperature and wind velocity 200 cm above the ground. Moreover, 34 series of microclimatic measurements, carried out during 1995–2004 in the warm part of the year in the area mentioned above, were used. Using this method a total of 12 000 data on air temperature, wind velocity, and magnitude of air cooling down (dry bioclimatic cooling down) measured with the Hill´s cathetometer, were gathered. The manual microclimatic measurements were carried out using the psyhrometer TB91AI, hand anemometer, and Hill´s cathetometer. The microclimatic measurements were conducted from May to August, while the automatic measurements during the entire year.
H1 = (36,5 - t) x (0,20 + 0,40 x √v)
when v Ł 1.0 m/s
H1 = (36,5 - t) x (0,13 + 0,47 x √v)
when v > 1.0 m/s
where: t – air temperature (oC), v – wind speed (mĚs-1).
Table 1.Localization of meteorological station1 and microclimatic post2.
| Name of station/post | Geographical coordinates | Altitude m a.s.l | Landform |
| Jaworzyna Krynicka1 | 49o26´ N
20o58´ E | 1114,0 | summit |
| Krynica Dolna2 | 49o26´ N 20o58´ E | 480,0 | lower parts of gentle valley slopes |
| Krynica Czarny Potok2 | 49o26´ N
20o´ E | 560,0 | valley bottom |
| Kopciowa1 | 49o27´ N 20o58´ E | 720,0 | upper parts of gentle valley slopes |
| Wojkowa2 | 49o25´ N 20o58´ E | 820,0 | hilltop |
| Złatna2 | 49o31´ N 19o14´ E | 700,0 | hilltop |
| Istebna Wyrchczadeczka1 | 49o39´ N
18o53´ E | 680,0 | mountain ridge |
| Zelków2 | 50o09´ N 19o48´ E | 260,0 | valley bottom |
| Kobylany2 | 50o08´ N 19o46´ E | 270,0 | valley mouth |
| Będkowice2 | 50o10´ N 19o45´ E | 400,0 | summit |
The magnitude of air cooling down was calculated by Hill´s formulae.
In the remaining cases the direct measurements of the biometeorological index were used. For each twenty-four hour period the following were computed:
a. mean magnitude of air cooling down (arithmetic average) on the basis of hourly measurements taken from 07:00 to 19:00 (13 hours);
b. magnitude of air cooling down on the basis of measurements from the observation time-limit II (single value);
c. mean magnitude of air cooling down (weighted average) on the basis of measurements taken at the I, II and III time-limits:
Hśr = (HI + 2HII + HIII)/4
It was assumed that the mean calculated from all hourly values from 07:00 to 19:00 (13 values) or 30 min values (25 values) of the same period of time represents in the best way the conditions of bioclimatic cooling down during the day. It reflects the whole day variation of this index, especially of its most variable element, the wind. In connection with this the following indexes were compared with each other:
a. value of bioclimatic cooling down from the observation time-limit II with value of the index calculated on the basis of measurements taken at 13 (25) hours, from 07:00 to 19:00;
b. value of bioclimatic cooling down calculated by the formula [3] as proposed by the author with value of the index calculated on the basis of measurements taken at 13 (25) hours. Values of index H, shown in tables, are expressed in W?m-2.
The indexes were calculated using the program BioKlima 2.3 (Błażejczyk 1999). The significance of differences between means of individual measurement series was estimated using t-test for dependent variables. In two cases the measurement series did not have normal distribution. In those cases the Wilcoxon´s test was used as a non-parametric equivalent of t-test for dependent variables. Calculations were made using STATISTICA 6.1 (StatSoft, Inc. 2004)
Results
Generally the mean value of a dry bioclimatic cooling down was lower than the diurnal mean, as well as the mean calculated with formula [3] (Tables 2-11). Such a relationship occurred regardless of the length of the measurement series, season of the year or weather conditions. The effect of weather conditions was only expressed by the fact that differences were smaller in the case of great cloudiness during the whole day and lack of advection. They were most evident in the case of minimum values (Tables 2-11). A different situation was observed in the case of standard deviation. In all cases it reached the highest value for the cooling down index from the time-limit II. This confirmed the hypothesis that magnitude of air cooling down in the time-limit II was characterized by a great variation, and this made the search for analogy to diurnal distribution of this index difficult. The relationship between maximum and minimum values confirmed this characteristic, since in all annual and monthly series the amplitudes of the bioclimatic cooling down for the time-limit II were distinctly higher than for the means (Tables 2, 11 and 12).
Table 2. The cooling power index on Jaworzyna Krynicka Meteorological station.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 2001-12-01 do 2001-12-31 | II time limits | 365,0 | 1044,3 | 156,6 | 2818,2 | 495,6 |
| day | 365,0 | 1051,6 | 312,9 | 2746,7 | 437,5 |
| formula [3] | 365,0 | 1055,0 | 260,4 | 2859,0 | 430,3 |
| 2002-01-01 do 2002-01-31 | II time limits | 365,0 | 1001,4 | 122,3 | 2462,4 | 498,7 |
| day | 365,0 | 1029,1 | 163,2 | 2585,7 | 456,8 |
| formula [3] | 365,0 | 1035,2 | 165,8 | 2501,4 | 466,3 |
| 2003-02-01 do 2003-02-28 | II time limits | 365,0 | 1025,6 | 197,3 | 2722,2 | 491,4 |
| day | 365,0 | 1045,1 | 235,7 | 2610,1 | 478,4 |
| formula [3] | 365,0 | 1049,4 | 221,3 | 2580,0 | 472,9 |
where: n – number of days; Hśr, Hmin, Hmax – average, minimum, maximum value of cooling power index (WĚm-2); ó – standard deviation (WĚm-2).
Table 3. The cooling power index on Krynica Dolna microclimatic post.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 2000-06-23 do 2000-06-27 | II time limits | 120,0 | 253,5 | 202,2 | 850,2 | 377,2 |
| day | 120,0 | 635,5 | 334,0 | 903,1 | 285,6 |
| formula [3] | 120,0 | 510,5 | 321,6 | 936,8 | 263,5 |
Table 4. The cooling power index on Ujsoły Złatna microclimatic post.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 1997-06-03 do 1997-06-09 | II time limits | 120 | 645,4 | 460,5 | 752,3 | 168,6 |
| day | 120 | 499,9 | 576,5 | 696,5 | 103,8 |
| formula [3] | 120 | 446,5 | 565,6 | 684,3 | 91,2 |
Table 5. The cooling power index on Wojkowa microclimatic post.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 1995-07-07 do 1995-07-11 | II time limits | 120 | 718,4 | 426,4 | 894,2 | 209,8 |
| day | 120 | 611,6 | 444,6 | 721,3 | 147,5 |
| formula [3] | 120 | 602,5 | 425,1 | 759,3 | 168,6 |
Table 6. The cooling power index on Kynica Czarny Potok microclimatic post.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 2000-06-03 do 2000-06-09 | II time limits | 120 | 411,1 | 220,7 | 777,0 | 317,2 |
| day | 120 | 453,0 | 262,3 | 683,2 | 213,5 |
| formula [3] | 120 | 447,7 | 270,1 | 654,3 | 193,5 |
Table 7. The cooling power index on Będkowice microclimatic post.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 2003-06-03 do 2003-06-09 | II time limits | 120 | 490,0 | 211,5 | 658,8 | 215,5 |
| day | 120 | 513,2 | 281,5 | 731,3 | 210,1 |
| formula [3] | 120 | 526,1 | 312,0 | 780,7 | 205,3 |
Table 8. The cooling power index on Kobylany microclimatic post.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 2004-06-03 do 2004-06-09 | II time limits | 120 | 422,6 | 178,2 | 940,3 | 390,4 |
| day | 120 | 458,6 | 284,0 | 707,5 | 239,3 |
| formula [3] | 120 | 477,3 | 307,0 | 726,7 | 232,1 |
Table 9. The cooling power index on Zelków microclimatic post.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 2003-07-03 do 2003-07-09 | II time limits | 120 | 564,4 | 396,6 | 722,3 | 128,0 |
| day | 120 | 575,5 | 415,5 | 686,9 | 115,0 |
| formula [3] | 120 | 581,1 | 402,9 | 656,3 | 101,0 |
Table 10. The cooling power index on Istebna Wyrchczadeczka meteorological station.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 1997-08-11 do 1997-08-15 | II time limits | 48 | 310,0 | 221,1 | 500,5 | 98,4 |
| day | 48 | 350,4 | 235,4 | 435,0 | 49,9 |
| formula [3] | 48 | 368,1 | 264,6 | 415,3 | 59,6 |
Table 11. The cooling power index on Krynica Kopciowa meteorological station.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 2002-11-01 do 2002-11-30 | II time limits | 744 | 1089,5 | 191,1 | 2452,2 | 655,5 |
| day | 744 | 1008,1 | 378,7 | 2074,8 | 460,2 |
| formula [3] | 744 | 1019,5 | 265,6 | 2134,1 | 502,0 |
| 2002-12-01 do 2002-12-31 | II time limits | 744 | 965,5 | 281,3 | 2152,7 | 565,5 |
| day | 744 | 1013,0 | 477,7 | 2389,2 | 392,2 |
| formula [3] | 744 | 997,6 | 381,1 | 2304,9 | 396,8 |
| 2003-01-01 do 2003-01-31 | II time limits | 744 | 1220,7 | 303,1 | 2451,6 | 728,4 |
| day | 744 | 1101,6 | 338,1 | 2135,8 | 567,8 |
| formula [3] | 744 | 1121,1 | 336,2 | 2106,8 | 555,3 |
Table 12. The cooling power index on Istebna Wyrchczadeczka meteorological station.
| Date | Formula | n | Hśr | Hmin | Hmax | s |
| 2004-11-01 do 2004-11-30 | II time limits | 720 | 940,7 | 236,1 | 2170,7 | 433,1 |
| day | 720 | 917,9 | 263,3 | 1836,6 | 327,0 |
| formula [3] | 720 | 903,3 | 309,2 | 1816,5 | 313,5 |
Also in series of measurements taken every 30 minutes the amplitude of the magnitude of air cooling down turned out to be considerably higher than in the remaining ones. The situation in shorter series was similar, although there were cases when the diurnal mean had the highest amplitude (table 7).
The average difference between values of the diurnal mean and weighted mean of the cooling down index for annual series was 0.4%. The difference between the value from the time-limit II and the diurnal mean was 1.8%, and between the value from the time-limit II and the weighted mean was 2.2%. These differences for the monthly series were 1.4%, 7.6% and 6.2% respectively. In series of several days these differences had increased. Between the means they were 5.0% on the average, and between the means from the time-limit II and the diurnal means 18.5%, and between the values of the means from the time-limit II and the weighted means 19.3%. In single cases the differences between the indexes from the time-limit II and the means were as much as 84 %. A detailed analysis of the results of this study did not confirm the relationship between values of the index calculated by different methods and the location of the weather station.
Conclusions
The variation of the index of bioclimatic cooling down during the day is often very high, and this was confirmed by the results of hourly measurements of this element in all series of measurements. To identify the diurnal, or twenty-four hour value of a dry bioclimatic cooling down with the value measured or calculated for the observation time-limit II only is an oversimplification. In single cases the differences reach even 300 %. During the year (annual mean) the magnitude of bioclimatic cooling down from the observation time-limit II has values distinctly lower than the means computed on the basis of hourly values, and it should not be used in climatological elaborations, with the exception of those stations which have no data from morning and evening observation time-limits. The results of this study confirmed the hypothesis that the formula proposed by the author to calculate the weighted arithmetic mean with the magnitude of bioclimatic cooling down from the observation time-limit II as the weight factor is the most proper method of determination of the diurnal magnitude of air cooling power index H, at stations taking meteorological measurements in three observation time-limits.
Polecane
książki z księgarni medycznej BORGIS:
Piśmiennictwo
1. Błażejczyk K, Błażejczyk M., 1999, Program komputerowy BioKlima dla Windows. Wersja 2.2.
2. Kozłowska-Szczęsna T., Błażejczyk K., Krawczyk B., 1990, Bioklimatologia człowieka. Metody i ich zastosowanie w badaniach bioklimatu Polski. Monografie, 1. IGiPZ PAN. Warszawa, 1997.
3. StatSoft, Inc., 2004, STATISTICA (data analysis software system), version 6.1 www.statsoft.com.

Adres do korespondencji:
dr inż. Grzegorz Durło
Akademia Rolnicza im. H. Kołłątaja w Krakowie
Katedra Ochrony Lasu i Klimatologii Leśnej,
Al. 29 Listopada 46, 31-425 Kraków, Polska
telefon: (012) 411-91-44 wew. 481
e-mail: rldurlo@cyf-kr.edu.pl
Balneologia Polska 3-4/2005Pozostałe artykuły z numeru 3-4/2005: