International Journal of Science, Environment ISSN 2278-3687 (O) and Technology, Vol. 6, No 4, 2017, 2216 – 2223 2277-663X (P)

SPATIAL AND TEMPORAL CLASSIFICATION OF BASED ON TEMPERATURE HUMIDITY INDEX (THI) Prasad A1*, Kannan A2, Gleeja V.L3, Anil K. S4, Joseph Mathew5 and Girish Varma G6 1Assistant Professor, Department of Livestock Production Management, College of Veterinary and Animal Sciences, Mannuthy, Thrissur, 2Associate Professor and Head, LRS, Thiruvizhamkunnu, Kerala, 3Assistant Professor, Department Statistics, College of Veterinary and Animal Sciences, Mannuthy, Thrissur, 4Professor and Head, ULF Mannuthy, Thrissur, 5Professor, Department of Livestock Production Management, College of Veterinary and Animal Sciences, Mannuthy, Thrissur, 6Professor, Veterinary Physiology, College of Veterinary and Animal Sciences, Mannuthy E-mail: [email protected] (*Corresponding Author)

Abstract: The study was conducted with the objectives to delineate THI zones in Kerala with respect to thermal comfort of dairy cattle. Daily meteorological data for the past ten years (from 2002 to 2011) was collected from thirteen different weather stations located in different agro ecological zones of Kerala. The ten year average of daily THI, temperature and humidity of each weather station and of the entire state of Kerala were calculated. Moving averages were used to find periods. The 60 day moving average of THI at all the thirteen weather stations for the past ten years were calculated. Altitude was found to have significant effect on THI (p<0.05). The present study suggests following linear equation to predict THI from altitude. THI = -0.0151 X Altitude in meters + 91.923 (R² = 0.9339) There is similarity in pattern of THI in all zones and the highest THI period in Kerala is from 19-Feb to May 27(98 days). Keywords: Thermal stress, THI, Climatic classification, Kerala climate.

INTRODUCTION

There is no place in the world where climate is perfectly conducive at all times for dairy cattle. Climatic constraints exist at somewhere or at some time in every geographical entity. In a hot humid tropical environment thermal stress is the major climatic constraint. Kerala is tropical humid state in having high percentage of crossbreds among cattle population. For cattle, increased air temperature and humidity measured as Temperature Humidity Index (THI) above critical thresholds are related to low dry matter intake (DMI) and to reduced efficiency of milk production (West, 2003). The Temperature Humidity Index (THI) combines the values of temperature and humidity to measure intensity of heat exposure through a formula to evaluate a unique thermo physiological effect (Jendritzky et al., 2002). There were climatic classification attempts with dairy cattle interpretations using temperature Received June 5, 2017 * Published Aug 2, 2017 * www.ijset.net

2217 Prasad A, Kannan A, Gleeja V.L, Anil K. S, Joseph Mathew and Girish Varma G and humidity alone and using thermal indices. THI based zonation was done in Mediterranean area by Gantner et al 2011. Hence the study was conducted with the objectives to delineate THI zones in Kerala with respect to thermal comfort of dairy cattle. MATERIALS AND METHODS Daily meteorological data for the past ten years (from 2002 to 2011) was collected from thirteen different weather stations located in different agro ecological zones of Kerala, India. Which were maintained by Indian Meteorological Department (IMD) and Kerala Agricultural University (KAU). The following method was used for statistical analysis to find out the spatial and temporal classification of Kerala based on THI. From the daily maximum and minimum temperature, the daily average temperature was calculated and from the morning and evening relative humidity values average relative humidity was calculated. Ten years daily average of both temperature and humidity was used to find out the average THI. THI calculated using maximum temperature was designated as Maximum THI. From among the different THI formulae, the THI was calculated using the formula THI = db ◦F- {(0.55-0.55RH) (db ◦F-58)} (LPHSI 1991). Each weather station was assumed as representing the meso climatological jurisdiction for evolving a zonation based on THI (Rao 2008). The ten year average of daily THI, temperature and humidity of each weather station and of the entire state of Kerala were calculated. The average THI of each station was compared using paired T test in SPSS software for spatial classification of THI zones in Kerala and (GIS) software was used to prepare THI based map of Kerala. The comparative relation between average THI altitude was also analysed. Temporal Classification Moving averages were used to find periods. The 60 day moving average of THI at all the thirteen weather stations for the past ten years were calculated. The highest 60 day moving average of THI was taken as the most probable highest THI period in a year. RESULTS AND DISCUSSION Classification of Kerala based on THI The weather stations across the state in the decreasing order of their THI based on ten years daily data and corresponding average THI are given in Table 1 and Fig 1. The table also contains the latitude, longitude, altitude, period of highest THI of that particular station and average of maximum THI and average THI of respective stations. Paired t tests conducted on

Spatial and Temporal Classification of Kerala Based on …. 2218 average THI delineated the weather stations into six statistically significantly different THI groups (p<0.05). Average THI delineated the weather stations into six statistically significantly different THI groups (p<0.05). Altitude was found to have significant effect on THI (p<0.05). The figure 1 shows the variation in THI with relation to altitude in Kerala. The present study suggests following linear equation to predict THI from altitude. THI = -0.0151 X Altitude in meters + 91.923 (R² = 0.9339). Table 1. THI Based Spatial and Temporal Details of Weather Stations of Kerala (2002-2011) Highest 60 day THI Latitude Maximum Longitude Altitude Period Average Stations in THI in degrees (m) Calendar THI degrees day March May Panniyur 12.07 75.40 14 18 17 101.08 95.32

Punalur 9.01 76.93 56 2 1 96.25 92.61

Cochin 10.15 76.40 9 16 15 96.03 92.43

Thiruvananthapuram 8.48 76.92 5 21 20 95.80 91.42

Vellayani 8.43 76.99 21 20 19 95.25 91.26

Pattambi 10.80 76.18 62 18 17 95.69 90.95

Kozhikode 11.13 75.95 88 3 2 95.25 90.91

Vellanikkara 10.54 76.28 22.5 17 16 95.52 90.80

Kottayam 9.59 76.52 3 15 14 94.67 90.75

Peelikode 12.19 75.14 10 27 26 94.94 90.18

Alleppy 9.49 76.33 2 3 2 93.58 90.09

Wayanad 11.61 76.21 921 18 17 84.17 78.79

Pampadumpara 9.79 77.15 1040 19 18 81.82 75.48

Same data was used to prepare the THI map of Kerala using GIS (Fig1).

2219 Prasad A, Kannan A, Gleeja V.L, Anil K. S, Joseph Mathew and Girish Varma G

Fig 1. THI map of Kerala

Spatial and Temporal Classification of Kerala Based on …. 2220

Table 2. Grouping of Weather Stations based on THI

Statistical grouping Stations under each group Average THI of weather stations

LI Pampadumpara 75.52

L2 Wayanad 78.73

Alappy H1 89.87 Peelikkode

Kozhikkode H2 Kottayam 90.66 Vellanikkra

Pattambi H3 Trivandrum 91.07 Vellayani

Cochin H4 93.40 Panniyur

Figure 2 shows the THI of each Julian day for the different weather stations. There is similarity in pattern of THI in all stations of Kerala. Figure 3 shows the average of averages of all thirteen weather stations depicted as the Average THI of the state of Kerala. Highest THI period in Kerala is from 19-Feb to May 27(98 days). From the 60 day moving average of THI of Kerala (Fig 4), the highest 60 day period of THI was found to be between March 17 and May 16. Daily averages of THI of all the thirteen weather stations for ten years plotted graphically illustrated the most probable average THI of each Julian day in the respective area (Fig 5). The figure also show that though there was difference in THI values among different stations, showed a similar pattern of the THI graph in a year.

2221 Prasad A, Kannan A, Gleeja V.L, Anil K. S, Joseph Mathew and Girish Varma G

Fig. 2. Daily average THI of different locations in Kerala

Fig. 3. Average THI of Kerala in a year

Many earlier works had classified the climate of Kerala into seasons. (Somanathan, 1980; Nishant, 2009; Sowmya, 2012). These were based on monthly average temperatures. English calendar months are not apposite to the monsoon based climate of Kerala. moreover, in thermal stress study in which even the minute to minute changes in thermal factors affecting the physiology of the animal; monthly or weekly averages are far from accurate and scientifically implausible over simplifications.

Spatial and Temporal Classification of Kerala Based on …. 2222

Fig. 4. Sixty day moving average THI of Kerala

96

94 92 90 88

86 84 Moving Average THI 82 1 31 61 91 121 151 181 211 241 271 301 331 361 Days in a calender year

The present study relied upon daily averages to accurately find out which are the high stress and which are period of low stress. So rather than seasonality classification this work has identified high thermal stress periods which are relevant in stress management of dairy cattle. Conclusion This classification can be used for adaptability and stress evaluation studies of animals especially dairy cattle. So, the spatial classification and stress status between the zones could be summarised as, 1. Extreme high stress zone (H4) - Pockets of severe stress 2. High stress zone (H1, H2, H3) – Major areas of high THI zone in the coastal and midland of the state of Kerala 3. Moderately comfortable zones (L2) – Comfortable THI zones with differing comfort levels at different altitudes closely related to high range 4. Comfortable zones (L1) –Comfortable THI zones of the high range References [1] Gantner, V., Mijic P., Kutorovac, K., Solic, D and Gantner, R 2011. Daily production of dairy cattle. Mljekarstvo 61(1) : 56-63. [2] Jendritzky, G., Maarouf, A., Fiala, D. and Staiger, H. 2002. An update on the development of a Universal Thermal Climatic Index. In Proc. 15 th AMS Conf. on Biomet and Aeriobiology / 16 th Intl. Cong. Biomet. 129-133, Boston, Mass: AMS. [3] LPHSI. 1990. Livestock and Poultry Indices, Agriculture Engineering Guide, Clemenson University, Clemenson SC, 29634, USA. 2223 Prasad A, Kannan A, Gleeja V.L, Anil K. S, Joseph Mathew and Girish Varma G

[4] Nishant, P. 2009. Management protocol for averting drop of milk production in dairy cattle during summer. M.V.Sc. thesis, Kerala Agricultural University, Thrissur, 65p. [5] Rao, G.S.L.H.V.P. 2003. Agricultural Meteorology. Kerala Agricultural University, Thrissur, Kerala, India, p. 326. [6] Somanathan, V.L. 1980. Bioclimatological studies on dry matter intake and water consumption of growing livestock. M.V.Sc. thesis, Kerala Agricultural University. [7] Sowmya, T.C. 2012. Evaluation of adaptability and production performance of high yielding crossbred cattle in Kerala. M.V.Sc. thesis, Kerala Veterinary Animal Sciences University pp. 70. [8] West, J.W. 2003. Effects of heat stress on production in dairy cattle. J. Dairy Sci. 86: 2131-2144.