Relation between potential sheep heat-stress and meteorological

Similar documents
URBAN HEAT ISLAND IN SEOUL

Mesonet Cattle Comfort Advisor Description Edited May 26, 2016

Climate Indices and Indicators for Agriculture

Urban heat island in the metropolitan area of São Paulo and the influence of warm and dry air masses during summer

NOAA National Centers for Environmental Information State Summaries 149-FL. Observed and Projected Temperature Change

Seasonal & Diurnal Temp Variations. Earth-Sun Distance. Eccentricity 2/2/2010. ATS351 Lecture 3

Dr. Haritini Tsangari Associate Professor of Statistics University of Nicosia, Cyprus

What is Climate? Understanding and predicting climatic changes are the basic goals of climatology.

Content. Climate Parameters used for the measurement and recording of weather phenomena The Psychrometric Chart

IGAD Climate Prediction and and Applications Centre Monthly Bulletin, August May 2015

ESCI 344 Tropical Meteorology Lesson 7 Temperature, Clouds, and Rain

Climate and the Atmosphere

Chapter 3. Materials and Methods

Chapter 3 Section 3 World Climate Regions In-Depth Resources: Unit 1

Analytical Report. Drought in the Horn of Africa February Executive summary. Geographical context. Likelihood of drought impact (LDI)

National Wildland Significant Fire Potential Outlook

Wind: Global Systems Chapter 10

Temporal and Spatial Distribution of Tourism Climate Comfort in Isfahan Province

Average Weather In March For Fukuoka, Japan

LAB J - WORLD CLIMATE ZONES

DROUGHT IN MAINLAND PORTUGAL

September 2016 No. ICPAC/02/293 Bulletin Issue October 2016 Issue Number: ICPAC/02/294 IGAD Climate Prediction and Applications Centre Monthly Bulleti

INVESTIGATION FOR A POSSIBLE INFLUENCE OF IOANNINA AND METSOVO LAKES (EPIRUS, NW GREECE), ON PRECIPITATION, DURING THE WARM PERIOD OF THE YEAR

3) What is the difference between latitude and longitude and what is their affect on local and world weather and climate?

Why the Earth has seasons. Why the Earth has seasons 1/20/11

Mozambique. General Climate. UNDP Climate Change Country Profiles. C. McSweeney 1, M. New 1,2 and G. Lizcano 1

EXPLORING THE URBAN HEAT ISLAND INTENSITY OF DUTCH CITIES

Energy Systems, Structures and Processes Essential Standard: Analyze patterns of global climate change over time Learning Objective: Differentiate

Manfred A. Lange Energy, Environment and Water Research Center The Cyprus Institute. M. A. Lange 11/26/2008 1

IGAD Climate Prediction and Applications Centre Monthly Bulletin, August 2014

ATMOSPHERIC CIRCULATION AND WIND

SEASONAL AND DAILY TEMPERATURES

WINTER NIGHTTIME TEMPERATURE INVERSIONS AND THEIR RELATIONSHIP WITH THE SYNOPTIC-SCALE ATMOSPHERIC CIRCULATION

Application and verification of ECMWF products 2013

2012 Meteorology Summary

LECTURE 28. The Planetary Boundary Layer

1990 Intergovernmental Panel on Climate Change Impacts Assessment

Seasonal prediction of extreme events

Climatic Classification of an Industrial Area of Eastern Mediterranean (Thriassio Plain: Greece)

West Africa: The 2015 Season

Impacts of the climate change on the precipitation regime on the island of Cyprus

Outdoor Thermal Comfort and Local Climate Change: Exploring Connections

Summary and Conclusions

Model Calculation of Odour Sensation in the Vicinity of a Livestock Building: a Meteorological Analyses when Odour Occur

PROPOSAL OF SEVEN-DAY DESIGN WEATHER DATA FOR HVAC PEAK LOAD CALCULATION

CLIMATE-SMART AGRICULTURE: THE POTENTIAL ROLE OF MODELLING AS PART OF

RELATION BETWEEN AIR POLLUTION EPISODES AND DISCOMFORT INDEX IN THE GREATER ATHENS AREA, GREECE

FINAL EXAM PRACTICE #3: Meteorology, Climate, and Ecology

Ganbat.B, Agro meteorology Section

Urban-rural humidity and temperature differences in the Beijing area

Fog Hazards in Punjab

Analysis of meteorological measurements made over three rainy seasons in Sinazongwe District, Zambia.

Factors That Affect Climate

Meteorology. Chapter 15 Worksheet 1

Impact on Agriculture

Application and verification of ECMWF products 2014

The Atmosphere and Atmospheric Energy Chapter 3 and 4

Monthly Overview. Rainfall

Climate. Annual Temperature (Last 30 Years) January Temperature. July Temperature. Average Precipitation (Last 30 Years)

Extremes Events in Climate Change Projections Jana Sillmann

CHAPTER 3. The sun and the seasons. Locating the position of the sun

Topic 1 The Atmosphere and Atmospheric Variables

Prentice Hall EARTH SCIENCE

IGAD CLIMATE PREDICTION AND APPLICATIONS CENTRE (ICPAC) UPDATE OF THE ICPAC CLIMATE WATCH REF: ICPAC/CW/NO. 24, AUGUST 2011

Explain the impact of location, climate, natural resources, and population distribution on Europe. a. Compare how the location, climate, and natural

Diurnal Variation of Rainfall Associated with Tropical Depression in South China and its Relationship to Land-Sea Contrast and Topography

Using the Mesonet Cattle Comfort Advisor

Climate concepts & climate zones

Climate Change Impact on Air Temperature, Daily Temperature Range, Growing Degree Days, and Spring and Fall Frost Dates In Nebraska

Average Weather For Coeur d'alene, Idaho, USA

Sudan Seasonal Monitor 1

La circulation atmosphérique dans le sud du Portugal et l'effet du réservoir d'alqueva. Cas d étude ALEX 2014

Climatic Zoning for Energy Efficient Buildings in the Palestinian Territories (the West Bank and Gaza)

Areas of the World with high Insolation

CH. 3: Climate and Vegetation

Monthly overview. Rainfall

WMO technical conference on climate services building on clips legacy

3 Temperate and Polar Zones

UPDATE OF REGIONAL WEATHER AND SMOKE HAZE FOR MAY 2015

Zambia. General Climate. Recent Climate Trends. UNDP Climate Change Country Profiles. Temperature. C. McSweeney 1, M. New 1,2 and G.

Using Temperature and Dew Point to Aid Forecasting Springtime Radiational Frost and/or Freezing Temperatures in the NWS La Crosse Service Area

International Journal of Science, Environment and Technology, Vol. 6, No 4, 2017,

2015: A YEAR IN REVIEW F.S. ANSLOW

Torrential events on the Spanish Mediterranean coast (Valencian Region). Spatial precipitation patterns and their relation to synoptic circulation

Chapter 9 Atmosphere Study Guide

SUDAN METEOROLOGCIAL AUTHORITY

Lesson- 6 WEATHER AND CLIMATE

Application and verification of ECMWF products: 2010

CLIMATE CLIMATE STUDYING CLIMATE. PART II: People and their 3/20/2018. GEOG 101 Part II People and their Physical Environment R E M I N D E R S

Suriname. General Climate. Recent Climate Trends. UNDP Climate Change Country Profiles. Temperature. C. McSweeney 1, M. New 1,2 and G.

CLIMATE DEFINED CLIMATE INFLUENCES STUDYING CLIMATE R E M I N D E R S. PART II: People and their 10/19/2018

Making Rain on Arid Regions The GESHEM Rain System

Studies on adaptation capacity of Carpathian ecosystems/landscape to climate change

Environmental Science Chapter 13 Atmosphere and Climate Change Review

NOTES AND CORRESPONDENCE. Seasonal Variation of the Diurnal Cycle of Rainfall in Southern Contiguous China

Weather and Climate Summary and Forecast Summer 2017

C Y P A D A P T. M. Petrakis C. Giannakopoulos G. Lemesios.

Analysis of Rainfall and Other Weather Parameters under Climatic Variability of Parbhani ( )

Flood hazard mapping in Urban Council limit, Vavuniya District, Sri Lanka- A GIS approach

The Atmosphere. Importance of our. 4 Layers of the Atmosphere. Introduction to atmosphere, weather, and climate. What makes up the atmosphere?

Transcription:

Ref: C628 Relation between potential sheep heat-stress and meteorological conditions Dimitris K. Papanastasiou and Thomas Bartzanas, Laboratory of Agricultural Engineering and Environment, Institute for Research and Technology of Thessaly, Centre for Research and Technology Hellas, Dimitriados St 95 & Pavlou Mela St, 3rd floor, GR-38333 Volos, Greece Constantinos Kittas, Laboratory of Agricultural Constructions and Environmental Control, Department of Agricultural Crop Production and Rural Environment, University of Thessaly, Fytokou St., GR-38446 N. Ionia, Magnisia, Greece Abstract It is well known that heat-stress affects the body growth, the biological functions and the productive and reproductive characteristics of sheep. Heat-stress levels are significantly determined by the prevailing meteorological conditions. The aim of this paper is to investigate qualitatively and quantitatively the impact of some meteorological parameters (i.e. incoming solar radiation, wind speed, wind direction, rainfall) on sheep s thermal comfort conditions during summer. Data recorded during the summers of the period 27 212 by an automated meteorological station that operates in a rural area near the east coast of central Greece is analyzed. Potential sheep s heat-stress is assessed using the Temperature Humidity Index (THI). Daytime (8: 2: local time) and night time (21: 7: local time) hours were examined separately. During daytime sheep were exposed to potential heat-stress (THI 22.2) in 84% of the hours. In more than half of the hours (53%), THI exceeded the extreme severe heat-stress threshold (i.e. 25.6). During night time sheep were exposed to potential heat-stress in 32% of the hours. The hourly mean of THI values remained higher than the extreme severe heat-stress threshold approximately during the hours 11: 2: (local time). As expected, solar radiation is related to adverse heat-stress conditions. Even low irradiance values were associated with high hourly THI values. 8% of the hourly THI values exceeded the severe heat-stress threshold (i.e. 23.3) when the irradiance was over 421 W m -2, while when irradiance was over 841 W m -2 the corresponding value was 94%. Wind speed was related to unfavorable thermal comfort conditions both during daytime and night time hours. This fact could be mainly attributed to mesoscale circulations that develop in the greater area. Additionally, THI levels were much lower during rainy days than during non rainy days both during daytime and night time hours. Keywords: sheep, heat-stress, temperature humidity index, meteorological conditions 1 Introduction One of the major factors that affect animal health and welfare is climate conditions, as unfavourable conditions can trigger heat-stress to ruminants. Heat-stress affects crucially the body growth, the biological functions and the productive and reproductive characteristics of sheep, therefore, it is considered as a serious threat to their short- or long-term thermal welfare. Silanikove (2) stated that in the long-term growth, milk production and reproduction of ruminants are impaired under heat-stress due to changes in biological functions. Sevi et Proceedings International Conference of Agricultural Engineering, Zurich, 6-1.7.214 www.eurageng.eu 1/8

al., 21) stated that high temperatures induce adverse effects on the thermal and energy balance, the mineral metabolism, the immune function, the udder health, and the milk production of lactating ewes during summer under the Mediterranean climate. Also, Sevi et al. (22) reported that high temperatures may induce a worsening of nutritional properties associated with the fatty acid profile of ewe milk. Finocchiaro et al. (25) found that milk production yields of Mediterranean dairy sheep are affected by heat- stress. Panagakis and Chronopoulou (21) found that under abnormally hot summer conditions, dairy ewes had respiration rates above normal altering their shade seeking behavior. Sevi and Caroprese (212) stated that a reduction of thermal stress in dairy sheep results in cumulative beneficial effects such as lengthening of lactation, maintenance of good processing features of milk and a reduction in veterinary costs. The influence of meteorological factors on animals heat-stress has been studied by many authors. A widely accepted method to study animals heat stress is the application of the THI, which is estimated by exploiting temperature and relative humidity values (Marai et al., 27). It is well documented that a combination of high ambient temperatures and high relative humidity is unsuitable for sheep (Thwaites, 1985). A poor thermal environment can affect the incidence and severity of certain endemic diseases, as well as the animals thermal comfort, growth rate and milk yield (Papanastasiou et al., 213; Seedorf et al., 1998). Some authors have also incorporated measurements of solar radiation and wind speed to assess heat-stress conditions (Eigenberg et al., 25; Mader et al., 26). The economic importance of the sheep sector must be pointed out (Theocharopoulos et al., 27). Sheep farming is the largest livestock sector in Greece, accounting for 3% of the total value of livestock output. Sheep milk and meat are among the major agricultural commodities with a share of ~9% of the total value of agricultural production. Greece is among the main producers of sheep milk and meat in EU as 67 Ktonnes of sheep milk and 8 Ktonnes of sheep meat are produced annually, accounting for the 1% of total EU production. The aim of this study is to assess sheep s thermal comfort conditions during summer and to investigate qualitatively and quantitatively the impact of some meteorological parameters (i.e. incoming solar radiation, wind speed, wind direction, rainfall) on them. 2 Materials and methods 2.1 Data Hourly averaged values of temperature, relative humidity, incoming solar radiation, wind speed, wind direction and rainfall observed during the summers of the period 27 212 were exploited in this study. Data was recorded by an automated meteorological station that operates in Velestino, a rural area near the east coast of central Greece (figure 1), by the Laboratory of Agricultural Engineering and Environment, Institute for Research and Technology of Thessaly, Centre for Research and Technology Hellas. 2.2 Estimation of sheep s potential heat-stress Marai et al. (27) suggested that an appropriate climatic index to estimate the severity of sheep heat-stress is the THI given in equation 1, where T is the dry-bulb temperature ( C) and RH is the relative humidity (%). The same authors defined four heat-stress categories which are presented in Table 1. The hourly averaged values of temperature and relative humidity were exploited to calculate the hourly values of THI that were analyzed in this study. THI T.31.31 RH T 14.4 (1) Proceedings International Conference of Agricultural Engineering, Zurich, 6-1.7.214 www.eurageng.eu 2/8

Figure 1: Location, where measurements were recorded (Velestino, 22 o 45 E, 39 o 24 N). Left map: map of Greece; right map: focused on the greater Velestino area. Table 1: Definition of heat-stress categories according to THI values Heat-stress category THI < 22.2 absence of heat-stress 22.2 THI < 23.3 moderate heat-stress 23.3 THI < 25.6 severe heat-stress THI 25.6 extreme severe heat-stress 3 Results and Discussion Daytime (8: 2: local time) and night time (21: 7: local time) hours were examined separately in this study. 7 Frequency of occurence (%) 6 5 4 3 2 1 THI < 22.2 22.2 <= THI < 23.3 23.3 <= THI < 25.6 THI >= 25.6 Figure 2: Frequency of occurrence (%) of es during daytime (gray bars) and night time (black bars). The distribution of the hourly THI values during daytime and night time is presented in figure 2. Figure 2 shows that during daytime sheep were exposed to potential heat-stress Proceedings International Conference of Agricultural Engineering, Zurich, 6-1.7.214 www.eurageng.eu 3/8

(THI 22.2) in 84% of the hours. In more than half of the hours (53%), THI exceeded the extreme severe heat-stress threshold (i.e. 25.6). During night time sheep were exposed to potential heat-stress in 32% of the hours. Figure 3 shows that sheep were exposed to potential extreme severe heat-stress approximately during the hours 11: 2: (local time). On the contrary, the hourly average of THI values remained lower than the heat-stress threshold (i.e. 22.2) approximately during the hours 23: 9:. This fact shows that the nocturnal exceddances of the heat-stress threshold identified in figure 2 were spread all night long, being not often observed during specific hours. 29 27 THI value 25 23 21 19 17 2 4 6 8 1 12 14 16 18 2 22 24 Local Time Figure 3: Average diurnal variation of THI values. The horizontal lines correspond to THI thresholds (i.e. 22.2, 23.3 and 25.6) Table 2: Frequency of occurrence (%) of es per solar radiation classes Incoming solar radiation class (W m -2 ) - 21 211-42 421-63 631-84 841-15 THI < 22.2 4 24 11 4 1 22.2 <= THI < 23.3 8 12 9 5 5 23.3 <= THI < 25.6 21 24 31 22 21 THI >= 25.6 31 4 49 69 73 Table 2 verifies that sheep s heat-stress levels were strongly determined by the incoming solar radiation. Even when irradiance was extremely low (i.e. 21 W m -2 ), 6% and 31% of the hourly THI values were higher than the heat-stress threshold and the extreme severe heat-stress threshold, respectively. These percentage increased as irradiance increased, reaching 99% and 73%, respectively, when irradiance was very high (i.e. 841 15 W m -2 ). Figure 4 shows the distribution of the hourly THI values during daytime and night time in relation to the prevailing wind speed. This figure shows that all wind speed classes could be associated with unfavourable es both during daytime and night time hours. However, this relation was more pronounced during daytime. This fact shows that wind speed was not a factor that strongly affected heat-stress levels. However, this conclusion should be treated with caution and should be examined jointly with the variation of the wind direction. Figure 5 shows the prevailing wind directions per wind speed class during daytime. Figure 5a shows that when wind speed was very low (i.e. 2.1 m s -1 ) the wind blew from north directions. When wind speed increased (figures 5b and 5c), the wind also blew from southeast directions (i.e. 2.1 6.3 m s -1 ) and from west directions (i.e. 4.2 6.3 m s -1 ), while when speed was very high (i.e. 6.3 8.3 m s -1 ) (figure 5d) wind blew mainly from west directions. Taking into account the topography of the greater area (figure 1, right map), the following conclusions could be drawn: Proceedings International Conference of Agricultural Engineering, Zurich, 6-1.7.214 www.eurageng.eu 4/8

Frequency of occurence (%) 8 7 6 5 4 3 2 1 a Frequency of occurence (%) 8 7 6 5 4 3 2 1 b THI < 22.2 22.2 <= THI < 23.3 23.3 <= THI < 25.6 THI >= 25.6 THI < 22.2 22.2 <= THI < 23.3 23.3 <= THI < 25.6 THI >= 25.6 Figure 4: Frequency of occurrence (%) of es per WS classes during daytime (a) and night time (b). Forward slashed bars: < WS 2.1 m s -1 ; black bars: 2.1 m s -1 < WS 4.2 m s -1 ; gray bars: 4.2 m s -1 < WS 6.3 m s -1 ; backward slashed bars: 6.3 m s -1 < WS 8.4 m s -1. a. North directions could be associated to the transport of relatively cold air masses from the neighbouring mountains to the study area. Consequently, when wind speed was very low, non heat-stress conditions prevailed in the 27% of the examined hours (figure 4a). b. Southeast directions could be associated to the sea breeze circulation that develops frequently in the area in the summer (Papanastasiou & Melas, 29). When sea breeze develops, cool sea air is advected over the land, triggering a reduction in temperature and an increase in relative humidity. The temperature reduction triggers a decrease in apparent temperature, while the increase in relative humidity triggers an increase in apparent temperature. The fact that heat-stress conditions deteriorated when wind speed also blew from southeast directions (figure 4a) shows that the effect of relative humidity on THI was more pronounced than the effect of temperature. A similar conclusion has been drawn by Papanastasiou et al. (21) when studying human comfort levels in the area. c. West directions could be associated with the advection of warm air from the Thessaly plain, where noon temperatures remain at high levels during summer. This fact could explain the deterioration of heat-stress conditions when wind speed was very high (figure 4a). 337.5 22.5 337.5 22.5 315 45 315 45 292.5 a 67.5 292.5 b 67.5 27 5 1 15 2 25 3 9 27 5 1 15 2 25 3 9 247.5 112.5 247.5 112.5 225 135 225 135 22.5 18 157.5 22.5 18 157.5 Proceedings International Conference of Agricultural Engineering, Zurich, 6-1.7.214 www.eurageng.eu 5/8

337.5 22.5 337.5 22.5 315 45 315 45 292.5 c 67.5 292.5 d 67.5 27 5 1 15 2 25 3 9 27 5 1 15 2 25 3 9 247.5 112.5 247.5 112.5 225 135 225 135 22.5 18 157.5 Figure 5: Frequency of occurrence (radius axis, %) of wind direction (angle axis, o ) during daytime. a: < WS 2.1 m s -1 ; b: 2.1 m s -1 < WS 4.2 m s -1 ; c: 4.2 m s -1 < WS 6.3 m s -1 ; d: 6.3 m s -1 < WS 8.4 m s -1. 22.5 18 157.5 Table 3: Frequency of occurrence (%) of es during rainy and non rainy days, during daytime and night time Daytime Night time Rainy days Non rainy days Rainy days Non rainy days THI < 22.2 65 15 88 68 22.2 <= THI < 23.3 6 7 7 1 23.3 <= THI < 25.6 23 24 4 15 THI >= 25.6 6 54 1 7 Table 3 reveals that heat-stress conditions were worse during non rainy days than during rainy days, both during daytime and night time hours. During daytime and during non rainy days, the hourly THI value surpassed the severe and the extreme severe heat-stress threshold in 78% and 54% of the hours, respectively. The corresponding percentages for rainy days were 29% and 6%, respectively. During the daytime hours of the rainy days, cloudiness is increased, causing a reduction in the amount of solar radiation that reaches the ground. As it is above mentioned, the incoming solar radiation affected decisively the heatstress levels. Therefore, THI levels were diminished. 4 Conclusions This paper studied the impact of incoming solar radiation, wind speed, wind direction and rainfall on sheep s thermal comfort conditions that prevail during summer in a rural area near the east coast of central Greece. Sheep s potential heat-stress was assessed by means of THI. Daytime and night time hours were examined separately. The analysis showed that sheep were exposed to potential moderate, severe and extreme severe heat-stress in 8%, 24% and 53% of the daytime hours, while heat-stress was absent in 15% of the daytime hours. The corresponding percentages for the night time hours were 1%, 15%, 7% and 68%, respectively. Sheep were usually exposed to potential extreme severe heat-stress during the vast majority of daytime hours (approximately during the hours 11: 2:, local time), while they usually were not exposed to heat-stress during the night time hours. The analysis revealed that the major meteorological factor that affects THI levels was the incoming solar radiation, followed by the wind direction. Even low irradiance values were related to unfavorable heat-stress conditions. Additionally, low irradiance values were asso- Proceedings International Conference of Agricultural Engineering, Zurich, 6-1.7.214 www.eurageng.eu 6/8

ciated with increased cloudiness, a fact that contributed to the reduction of the hourly THI values during rainy days. Wind direction affected THI levels as either warm air masses or cold air masses rich in moisture could be advected to the study area. It was also found that the effect of relative humidity on THI was more pronounced than the effect of temperature, when cold air rich in moisture was advected to the area. The analysis showed that wind speed did not have a significant impact on THI levels 5 Acknowledgement This work was supported by the Joint Call of ERANET ICT-Agri C-2 projects Smart Integrated Livestock Farming: integrating user-centric & ICT-based decision support platforms SILF 6 References Eigenberg, R. A., Brown-Brandl, T. M., Nienaber, J. A., & Hahn, G. L. (25). Dynamic Response Indicators of Heat Stress in Shaded and Non-shaded Feedlot Cattle, Part 2: Predictive Relationships. Biosystems Engineering, 91 (1), 111-118. Finocchiaro, R., Van Kaam, J. B. C. H. M., Portolano, B., & Misztal, I. (25). Effect of heat stress on production of Mediterranean dairy sheep. Journal of Dairy Science, 88 (5), 1855-864. Mader, T. L., Davis, M. S., & Brown-Brandl, T. (26). Environmental factors influencing heat stress in feedlot cattle. Journal of Animal Science, 84, 712-719. Marai, I. F. M., El-Darawany, A. A., Fadiel, A., & Abdel-Hafez, M. A. M. (27). Physiological traits as affected by heat stress in sheep-a review. Small Ruminant Research, 71 (1-3), 1-12. Panagakis, P., & Chronopoulou, E. (21). Preliminary evaluation of the apparent short term heat-stress of dairy ewes reared under hot summer conditions. Applied Engineering in Agriculture, 26(6), 135-142. Papanastasiou, D. K., & Melas, D. (29) Climatology and impact on air quality of sea breeze in an urban coastal environment. International Journal of Climatology, 29, 35-315. Papanastasiou, D. K., Melas, D., Bartzanas, T., & Kittas, C., (21). Temperature, comfort and pollution levels during heat waves and the role of sea breeze. International Journal of Biometeorology, 54, 37-317. Papanastasiou, D. K., Zhang, G., Bartzanas, T., Panagakis, P., Norton, T., & Kittas, C., (213). Development of models to assess potential sheep heat-stress during heat waves. Joint European Conference on Precision Livestock Farming (EC-PLF 213), 1 12 September, Leuven, Belgium. In: Berckmans, D, & Vandermeulen, J. (eds). Precision Livestock Farming 13. 391-396. Seedorf, J., Hartung, J., Schröder, M., Linkert, K. H., Pedersen, S., Takai, H., Johnsen, J. O., Metz, J. H. M., Groot Koerkamp, P. W. G., Uenk, G. H., Phillips, V. R., Holden, M. R., Sneath, R. W., Short, J. L., White, R. P., & Wathes, C. M. (1998). Temperature and moisture conditions in livestock buildings in Northern Europe. Journal of Agricultural Engineering Research, 7 (1), 49-57. Sevi, A., Annicchiarico, G., Albenzio, M., Taibi, L., Muscio, A., & Dell Aquila, S. (21). Effects of solar radiation and feeding time on behavior, immune response and production of lactating ewes under high ambient temperature. Journal of Dairy Science, 84(3), 629-64. Proceedings International Conference of Agricultural Engineering, Zurich, 6-1.7.214 www.eurageng.eu 7/8

Sevi, A., & Caroprese, M. (212). Impact of heat stress on milk production, immunity and udder health in sheep: A critical review. Small Ruminant Research, 17(1), 1-7. Sevi, A., Rotunno, T., Di Caterina, R., & Muscio, A. (22). Fatty acid composition of ewe milk as affected by solar radiation and high ambient temperature. Journal of Dairy Research, 69(2), 181-194. Silanikove, N. (2). Effects of heat stress on the welfare of extensively managed domestic ruminants. Livestock Production Science, 67 (1-2), 1-18. Theocharopoulos, A., Melfou, K., & Papanagiotou, E., (27). A microeconomic approach for agricultural development: a dea application to Greek sheep farms. NEW MEDIT, 4, 48-53. Thwaites, C. J. (1985). Chapt. 2: Physiological responses and productivity in sheep. Vol. II: Ungulates. In: Yousef, M. K. (ed.). Stress Physiology in Livestock. 25-38, Boca Raton, FL: CRC Press. Proceedings International Conference of Agricultural Engineering, Zurich, 6-1.7.214 www.eurageng.eu 8/8