ORIGINAL PAPER Models to predict both sensible and latent heat transfer in the respiratory tract of Morada Nova sheep under semiarid tropical environment Vinícius Carvalho Fonseca1 & Edilson Paes Saraiva1 & Alex Sandro Campos Maia2 & Carolina Cardoso Nagib Nascimento2 & Josinaldo Araújo da Silva1 & Walter Esfraim Pereira1 & Edgard Cavalcanti Pimenta Filho1 & Maria Elivânia Vieira Almeida1 Received: 11 January 2016 /Revised: 6 September 2016 /Accepted: 22 September 2016 /Published online: 10 October 2016 # ISB 2016 Abstract The aim of this study was to build a prediction model both sensible and latent heat transfer by respiratory tract for Morada Nova sheep under field conditions in a semi- arid tropical environment, using easily measured physiologi- cal and environmental parameters. Twelve dry Morada Nova ewes with an average of 3 ± 1.2 years old and average body weight of 32.76 ± 3.72 kg were used in a Latin square design 12 × 12 (12 days of records and 12 schedules). Tidal volume, respiratory rate, expired air temperature, and partial vapor pres- sure of the expired air were obtained from the respiratory facial mask and using a physiological measurement system. Ewes were evaluated from 0700 to 1900 h in each day under shade. A simple nonlinear model to estimate tidal volume as a function of respiratory rate was developed. Equation to estimate the expired air temperature was built, and the ambient air temper- ature was the best predictor together with relative humidity and ambient vapor pressure. In naturalized Morada Nova sheep, respiratory convection seems to be amechanism of heat transfer of minor importance even under mild air temperature. Evaporation from the respiratory system increased together with ambient air temperature. At ambient air temperature, up to 35 °C respiratory evaporation accounted 90 % of the total heat lost by respiratory system, on average. Models presented here allow to estimate the heat flow from the respiratory tract for Morada Nova sheep bred in tropical region, using easily measured physiological and environmental traits as respiratory rate, ambient air temperature, and relative humidity. Keywords Genetic adaptation . Naturalized breed . Predictionmodel . Tidal volume Introduction Mathematical models have been developed to evaluate heat exchange at the respiratory surfaces in cattle and sheep (Stevens 1981; da Silva et al., 2002; Berman, 2005; Maia et al., 2005). These models were conceived using easily mea- sured environmental and physiological parameters to estimate heat dissipation through the respiratory tract via convection and evaporation under field conditions.Models of evaporation in the respiratory system of animals generally assume that expired air is saturated at a temperature lower than deep body temperature, evincing that environmental air temperature, hu- midity, and body temperature are important factors for the determination of evaporation in the respiratory tract (Stevens 1981; da Silva et al., 2002; Maia et al., 2005). Equations to estimate heat loss via the respiratory tract should take factors like species, breed, and environment which they were selected. Maia et al. (2005) showed that using functions estimated from animals bred in temperate cli- mates were inappropriate for animals bred under tropical con- ditions. They studied Holstein cows acclimatized to a tropical environment and observed that measures like tidal volume were underestimated when they used equations developed by Stevens (1981) for Holstein cows bred in a temperate zone. Electronic supplementary material The online version of this article (doi:10.1007/s00484-016-1255-3) contains supplementary material, which is available to authorized users. * Vinícius Carvalho Fonseca vinicius_fonseca86@hotmail.com 1 Universidade Federal da Paraiba, Paraiba, Areia, Brazil 2 Universidade Estadual Paulista, São Paulo, Jaboticabal, Brazil Int J Biometeorol (2017) 61:777–784 DOI 10.1007/s00484-016-1255-3 http://dx.doi.org/10.1007/s00484-016-1255-3 http://crossmark.crossref.org/dialog/?doi=10.1007/s00484-016-1255-3&domain=pdf In the case of sheep, da Silva et al. (2002) developed functions with Corriedale sheep, but these are evolutionarily distinct from breeds naturalized to the semiarid tropical environment. Morada Nova sheep are known to be well adapted to the semiarid environment of Brazil (Egito et al., 2002; Fonsêca et al., 2014). By contrast, Corriedale sheep were selected for temperate climates and are characterized by a wool coat. Studies have shown that evaporation in the respiratory tract is the main avenue of heat dissipation at high ambient temper- atures in wool sheep (Alexander and Williams, 1962; Hofmeyr et al., 1969; Hofman and Riegle, 1977). Starling et al. (2002) showed that at high ambient temperature, respi- ratory heat transfer accounted for approximately 70 % of total heat loss in Corriedale bred in a tropical zone. On the other hand, in a study of naturalized goats from a semiarid tropical zone, Maia et al. (2015) documented that evaporative heat loss by the respiratory tract was only 7 % of total heat loss under conditions of high temperature. In sum, specific physiological responses in animals are related to the environment that they are adapted to. Therefore, mathematical models of heat transfer should take this into account. No attempt has been made to evaluate the characteristics of the respiratory system in heat exchange in naturalized sheep bred in a semiarid tropical environment. Thus, the aim of this study was to build predictive models of both sensible and latent heat transfer in the respiratory tract of Morada Nova sheep under field conditions in a semiarid trop- ical environment, using easily measured physiological and environmental traits. Materials and methods Twelve dry Morada Nova, Ovis aries, ewes averaging 3 ± 1.2 years in age and 32.76 ± 3.72 kg in body weight were studied in a 12 × 12 Latin square design (12 days of records and 12 schedules). Observations were made from April to May 2015. Ewes were evaluated under ambient conditions in São João do Cariri, PB, Brazil (07° 23′ 27″ S, 36° 31′ 58″ W, 458 m altitude), from 0700 to 1900 h each day under shade. Physiological parameters were recorded in the first an- imal between 0700 and 0800 h, the second between 0800 and 0900 h, and so on. Thus, the twelfth animal was evaluated between 1800 and 1900 h. By the end of the trial, every animal had been evaluated on each of the 12 days in all the schedules. Ewes were fed twice a day (0530 and 1910 h) with forage (elephant grass hay) and concentrate (corn, soybean, wheat, and mineral supplementation) in a ratio of 74:26. Water was provided ad libitum.Daily feed consumption was estimated at approximately 1.5 % of average body weight. Sensible and latent heat flows from the respiratory tract were determined according to the methods described by Maia et al. (2005) and da Silva and Maia (2013). The sensible heat flow (CR; Wm−2) from the respiratory system is given by CR ¼ A−1VTpcP F TE−Tarð Þ ð1Þ and the latent heat flow (ER; W m−2) by ER ¼ A−1λVTF ψE−ψatmð Þ ð2Þ where Tar and TE are the temperatures of the ambient and expired air (°K), respectively; p is the air density (kg m−3), cP is the specific heat of the air (J kg−1 °C−1), λ is the latent heat of water vaporization (J g −1), VT is the tidal volume (L breaths−1), F is the respiratory rate (breaths s−1), and A is the body surface area of the animal, which is estimated as (Bennett, 1973): A ¼ 0:171BW −0:05025 ð3Þ where, BW is body weight of the animal (Kg). The absolute humidity of the ambient air (ψatm; g m−3) and of the expired air (ψE; g m−3) are given by ψatm ¼ Tar−12166:87e Tarð Þ ð4Þ ψE ¼ TE −12166:87e TEð Þ ð5Þ which e(Tar) and e(TE) are the partial vapor pressures (kPa) of the ambient and expired air, respectively. Thermal properties of the air (p, λ, and cP) were obtained from a simple function of Tar (Da Silva, 2008). The respiratory parameters tidal volume (L breaths−1), re- spiratory rate, expired air temperature (°C), and partial vapor pressure of the expired air (kPa) were obtained using a respi- ratory facial mask and a physiological measurement system (Maia et al., 2016). Ambient air temperature (°C; model HOBO, onset), relative humidity (%; model HOBO, onset), partial vapor pressure of the air (kPa), rectal temperature (°C), and wind speed (m s−1; model APM-360, USA) were also collected. In addition, black globe temperature (°C) was mea- sured with a thermocouple (Model type T, Salcas, Brazil. Accuracy 0.2 °C) inserted into the centre of a hollow 0.15- m-diameter copper sphere, matt black painted, which was placed 50 cm above the ground near the animals. The mean radiant temperature (K) was obtained according to da Silva et al. (2010). Prior to the study, the ewes were habituated and conditioned to the use of the respiratory mask using operant conditioning principles employing feed as a positive rein- forcement. This habituation period extended from December 2014 to March 2015. 778 Int J Biometeorol (2017) 61:777–784 During each breath, the inlet flow (inspired air) and the outlet flow (expired air) were carried through two valves of the measurement system (Fig. 1). The expired air leaving the facial mask was directed through a tracheal tube (MLA1015 Breathing Tube, ADInstruments, Australia) to a flow head (MLT1.000, ADInstruments, Australia). The flow head was directed to a gas mixing-chamber (MLA246, ADInstruments, Aus t r a l i a ) connec ted to a sp i romete r (ML141 , ADInstruments, Australia) providing the ventilation (VR; L s−1) and respiratory rates (min−1). From these measure- ments, we calculated the tidal volume (VT; L breaths−1). The gas-mixing chamber was connected to a Field Metabolic System (FMS-1201-05, Sable System, USA) through a plastic tube (Bevaline Tubing, Sable System, USA). Inside the tube, a sample of expired air (150 mL/min) was continuously pumped into a gas analyzer (H2O) tomeasure the partial vapor pressure of the expired air (kPa). These measurements were performed every time a facial mask was placed on an animal’s muzzle. Atmospheric vapor pressure was continuously recorded by an external H2O vapor analyzer (PV, kPa; RH-300, Sable System, USA) connected to a pump (SS4, Sable System, USA), which continuously aspirated atmospheric air (150 ml min−1) close to the mask’s inlet valve. The expired air temperature was recorded with a nasal temper- ature probe (MLT 415/AL, ADInstruments, Australia, ac- curacy ±0.3 °C) placed inside the facial mask. The rectal temperature (°C) was continuously recorded with a rectal thermocouple (MLT 1404, ADInstruments, Australia, ac- curacy ±0.3 °C). The FMS, spirometer, H2O vapor ana- lyzer, and thermocouples were connected to a data acqui- si t ion system (PowerLab 16/32, ADInstruments, Australia) and linked to a computer that continuously and simultaneously recorded all measurements at a rate of one observation per second. Data were analyzed using general linear models and least squares methods. Nonlinear regression was used to estimate VT as a function of the respiratory rate. Multiple regression was used to develop a model for F, TE, CR, and ER, as a function of several variables, given the initial model: yijkl ¼ αþ ai þ t j þ dk þ b1 Tarð Þ þ b2 RHð Þ þ b3e Tarð Þ þ b4 TBð Þ þ εijkl where yijkl is the l th observation of TE, F, CR, and ER in the ith animal, jth time of the day and kth day of records, α is the intercept, ai is the effect of the ith animal (i = 1,...,12), tj is the effect of jth time of the day (k = 0700,...,1800 h), and dk is the effect of the kth day of records (k = 1,...,12). The values of b1, b2, b3, and b4 are partial linear regression coefficients and ɛijkl is the random error associated with the lth observation on the ith animal, jth time of the day, and kth day of record. Starting with Tar, independent variables were successively added to the model until the largest coefficient of determination (R2) was reached. At each step, a set of simultaneous equations was built and then solved by stepwise selection. The statistical package used for all analyses was SAS version 9.3 (SAS Inst. Inc., Cary, NC). Results and discussion From the data obtained in our study, we estimated the rela- tionship between tidal volume and respiratory rate for wool- less Morada Nova sheep bred in semiarid environment (R2 = 0.95; Fig. 2) as VT L breath−1 � � ¼ 10:0351F−1:0376; ð6Þ VT ranged from 0.345 to 0.644 L/breath under ambient conditions of air temperature (25.07 to 35.56 °C), radiant mean temperature (25.22 to 46.39 °C), relative humidity (35.84 % to 78.47 %), and vapor pressure (1.42 to 2.82 kPa; Table 1). Previous findings of sheep described similar Fig. 1 Physiological measurement system using facial mask developed for sheep Int J Biometeorol (2017) 61:777–784 779 variation for tidal volume on resting conditions (Mitchell 1972; Hofman and Riegle 1977; Schmitd-Nielsen 1991; da Silva et al. 2002). For example, Schmidt-Nielsen (1991) mea- sured a tidal volume of 0.403 L/breath for a 60-kg sheep. Da Silva et al. (2002), working with Corriedale ewes under am- bient temperatures ranging from 17.5 to 40 °C, documented tidal volumes from 0.200 to 0.700 L/breath. Similarly, shorn and unshorn Dorset ewes presented tidal volumes from 0.150 to 0.400 L/breath (Hofman and Riegle, 1977). According to this model, an increase in F (breaths min−1) decreased the VT, consistent with other studies of sheep (Hofman and Riegle, 1977; da Silva et al. 2002; Silva, 2013) and cattle (Stevens, 1981; Maia et al., 2005). Variation in the respiratory rate as a function of ambient air temperature was described by the equation (R2 = 0.61): F ¼ 49:86−2:472Tar þ 0:047Tar2; ð7Þ F ranged from 16.26 to 26.83 (Fig. 3). When other vari- ables such as ambient vapor pressure, relative humidity, and rectal temperature were included in the model, R2 did not increase. In the present study, the range of ambient and mean radiant temperature was not sufficient to cause panting in the ewes. Working with wool-less Santa Inês sheep bred in a tropical environment, Neiva et al. (2004) documented close range for F under both shady and sunny conditions. However, large variation in F under a similar range of ambient temper- ature was observed in wool sheep under panting state (Hofman and Reagle, 1977; Johnson, 1991; da Silva and Starling, 2003; Silva, 2013). With Merino ewes, Johnson (1991) described variation in the respiration rate from 65 to 120 breaths per minute. Similarly, Silva (2013) observed re- spiratory rates of 155 to 180 breaths per minute in male Corriedale maintained in a climatic chamber (average ambient temperature of 30 °C). When we plotted the data of the present study using the nonlinear function of da Silva et al. (2002; VT (m3/ breath) = 0.0496F−0.463), the curve deviated from the da- ta, showing a large overestimate for the tidal volume (Fig. 2). Maia et al. (2005) worked with Holstein cows acclimatized to a hot environment and observed that when they used a function estimated by Stevens (1981) for Holstein cows bred in a temperate climate, the estimates 0,300 0,500 0,700 0,900 1,100 1,300 1,500 1,700 1,900 15 17 19 21 23 25 27 ht a er b L , e m ul o vl a di T -1 Respiration rate, min-1 Da Silva et al. (2002) equation present equation Fig. 2 Tidal volume as a function of the respiration rate Table 1 Average, minimum, and maximum values of air temperature (Tar), relative humidity (RH), mean radiant temperature (TRM), wind speed (WS), and ambient vapor pressure (e[Tar]) observed in São João do Cariri from April to May, 2015 Variable n Average Minimum Maximum Tar, °C 748 31.84 ± 0.27 25.07 37.02 RH, % 743 49.04 ± 1.24 25.22 42.39 TRM, °C 702 32.40 ± 0.30 35.84 78.47 WS, m s−1 452 0.32 ± 0.04 0 1.21 e(Tar), kPa 740 1.94 ± 0.05 1.42 2.82 16 18 20 22 24 26 28 24 26 28 30 32 34 36 R e sp ir a ti o n r a te , m in -1 Ambient air temperature, °C Fig. 3 Respiration rate as a function of ambient air temperature 780 Int J Biometeorol (2017) 61:777–784 deviated from the data, underestimating tidal volume. A similar pattern was found by Silva and Maia (2011) using a function provided by Gatenby (1986) to estimate cuta- neous evaporation in cattle as a function of surface tem- perature. These authors worked with Holstein cows in a tropical environment while Gatenby (1986) developed the equation for crossbred beef cattle. Correlation coefficients of TE with respect to ambient air temperature (Tar), ambient vapor pressure e(Tar), rectal tem- perature (TR), and relative humidity (RH) are shown in Table 2. TE is more a function of the ambient air temperature, relative humidity, and ambient vapor pressure than rectal tem- perature (Fig. 4a–d). Based on these relationships, we chose two models that estimate TE with good accuracy. The first (R2 = 0.9778) included Tar and RH and the second (R2 = 0.9764) included Tar and e(Tar). The inclusion of e(Tar) in the first function and RH in the second did not increase the R2. When we included TR in both equations, the estimate did not improve. Thus, the models are the following: TE Å 1 � � ¼ 14:56þ 0:64Tar−0:08RH þ 0:0007RH 2; ð8Þ TE Å 2 � � ¼ 12:31þ 0:65Tar−0:16e Tarð Þ: ð9Þ However, in both models, the inclusion of only ambient air temperature gave an R2 of 0.96. Thus, the equation presented in Fig. 4a accurately predicted TE as a function of Tar. Da Silva et al. (2002) described that the best model (R2 = 0.96) of TE for Corriedale sheep included Tar, RH, and deep body temperature (TR). Stevens (1981) reported that the best model of TE included ambient air temperature and relative humidity. Meanwhile, Maia et al. (2005) reported the best model of TE was a function of Tar (R2 = 0.91). Deep body temperature (TR) ranged from 38 to 39.18. The minimum and maximum values were observed at 0700 and 1800 h, respectively. Da Silva and Minomo (1995) reported a similar pattern for Corriedale sheep bred in a tropical environ- ment. According to our investigation, the amount of body ther- mal energy did not coincide with the highest values of expired air temperature or respiratory rate, which may be an indication that the great part of heat was dissipated in other ways (Fig. 5). Table 2 Correlation coefficients among ambient air temperature (Tar), partial air vapor pressure (Pvatm), rectal temperature (TR), relative humidity (RH), and expired air temperature (TE) Variable Tar e(Tar) TR RH TE P value Tar 1 −0.8270 0.3485 −0.9442 0.9588 0.001 e(Tar) – 1 −0.3854 0.9491 −0.7799 0.001 TR – – 1 −0.3820 0.4083 0.001 RH – – – 1 −0.8917 0.001 TE – – – – 1 0.001 24 25 26 27 28 29 30 31 32 33 34 35 36 25 26 27 28 29 30 31 32 33 34 35 36 37 E x p ir ed a ir t em p er at u re , °C Air temperature, °C T E = 11.53 + 0.666Tar R 2 = 0.976 30 35 40 45 50 55 60 65 70 75 80 85 Relative humidity, % T E = 48.08 - 0.445R H + 0.002R H 2 R 2 = 0.952 1,0 1,2 1,4 1,6 1,8 2,0 2,2 2,4 2,6 2,8 3,0 25 26 27 28 29 30 31 32 33 34 35 36 37 E x p ir ed a ir t em p er at u re , °C E x p ir ed a ir t em p er at u re , °C E x p ir ed a ir t em p er at u re , °C Ambient vapor pressure, kPa T E = 43.14 - 5.3487e(Tar) R 2 = 0.828 37,8 38,0 38,2 38,4 38,6 38,8 39,0 39,2 39,4 Rectal temperature, °C T E = -28691 + 1484.15T R - 19.17T R 2 R 2 = 0.692 25 26 27 28 29 30 31 32 33 34 35 36 37 25 26 27 28 29 30 31 32 33 34 35 36 37 Fig. 4 Expired air temperature as a function of ambient air temperature (a), relative humidity (b), ambient air vapor pressure (c), and rectal temperature (d) Int J Biometeorol (2017) 61:777–784 781 As previously reported in sheep and cattle (Da Silva et al., 2002; Maia et al., 2005), the sensible (W m−2) heat loss de- creased with increasing ambient temperature, while evapora- tive heat loss increased together with Tar (Fig. 6): CR ¼ 33:61−2:995Tar þ 0:092Tar2−0:001Tar3; R2 ¼ 0:76 � � ð10Þ ER ¼ −295:8þ 28:49Tar−0:909Tar2 þ 0:009Tar3 R2 ¼ 0:90 � � : ð11Þ Sensible heat loss by the respiratory tract in Morada Nova sheep seems to be of little significance for thermal balance, even under mild air temperatures (25 °C). In a similar situa- tion, da Silva et al. (2002) reported values of respiratory con- vection in Corriedale sheep that ranged 7 to 9Wm−2 when the thermal gradient (Δt) between ambient and expired air tem- perature reached 7 °C. In the present study, at 25 °C, the thermal gradient reached 2 or 3 °C, on average. Our study shows that latent heat loss is the main mechanism of respira- tory heat loss. At ambient temperatures above 35 °C, respira- tory cooling via evaporation accounted for 90 % of total heat lost. As expected, the amount of latent heat lost in Morada Nova sheep is small when compared with the findings of da Silva et al. (2002). There is large differences between Morada Nova (weight = 33 kg on average) and Corriedale (weight = 75 kg on average) regarding thermal energy content. In other words, heavier animals have a greater capacity for heat storage than lighter animals. It is evident that there are large differences between woolly and wool-less naturalized sheep bred in a semiarid tropical environment in terms of heat transfer by the respiratory tract. Maia et al. (2005) and Silva and Maia (2011) had already emphasized that the use of the same models for animals bred under tropical and temperate conditions was inappropriate. As a recent example, in work with Morada Nova ewes in a semiarid tropical region, Oliveira and Costa (2013) used models developed by da Silva et al. (2002) to calculate the heat loss by the respiratory system. They estimated that evap- orative heat transfer ranged from 12 to 38 W m−2, while re- spiratory convection was insignificant, ranging from −0.002 to 0.004 W m−2. However, in the present study, we observed values of 0.65–8.94 W m−2 for evaporative heat loss and 0.74–1.59 W m−2 for respiratory convection. Therefore, these results confirm that the use of the same functions without considering the environment where animals were selected could cause large errors and the results would not represent the animal’s thermal balance. Differences between the tropical semiarid zone and tem- perate regions revolve around the challenge of maintaining animal homeothermy. The main contrast between tropical and temperate latitudes is the radiant environment (Da Silva et al. 2012). In tropical regions, the mean radiant temperature (TMR) is usually close to or higher than the air and body temperatures, so animals gain heat by thermal radiation. In contrast, heat loss by radiation almost invariably occurs in animals in temperate latitudes. In other words, genetic adap- tation in the tropics gave rise to morphological and physiolog- ical characteristics that facilitate heat dissipation and protec- tion from solar radiation (pigmented skin). By contrast, ani- mals adapted to temperate regions have characteristics that help retain heat (great thermal insulation by the presence of 38 38,5 39 39,5 10 13 16 19 22 25 7 8 9 10 11 12 13 14 15 16 17 18 ○ R e c ta l te m p e ra tu re , ° C ■ R e sp ir a to ri o n r a te , m in -1 Time of the day 27 29 31 33 35 37 ∆ E x p ir e d a ir t e m p e ra tu re , °C P < 0.001 Fig. 5 Respiration rate, rectal temperature, and expired air temperature (mean ± SEM) during daylight 0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10 25 27 29 31 33 35 x S e n sib le h e a t, W m -² ° L a te n t h ea t , W m -² Ambient a ir temperature,°C Fig. 6 Heat loss from the convection and evaporation in the respiratory tract of Morada Nova sheep as a function of ambient air temperature 782 Int J Biometeorol (2017) 61:777–784 wool). These statements are notable when we observe the animal body surface area and coat traits. Silanikove (2000) pointed out that the larger body surface areas of smaller races make them more efficient heat dissipaters in warm climates while a smaller body surface area may help to conserve heat in cold climates. Coat surface characteristics can drive differences between individuals and determine the best avenue for heat dissipation under high temperature. The presence of long wool difficult the sensible (main surface convection) and latent heat loss from the skin surface. Several studies have shown that some long-wool sheep breeds selected for temperate zones, the re- spiratory system is the main avenue for heat transfer at high temperatures, but results are still conflicting (Cheviot and Scottish Blackface, Brockway et al., 1965; Merino, Hofmeyer et al. 1969; Dorset, Hofman and Riegle, 1977; Corriedale, Starling et al., 2002; Silva and Starling, 2003). Starling et al. (2002), working with Corriedale sheep in a tropical environment, reported that respiratory heat loss accounted for 70 % of total heat loss, on average. Similar results were reported by Brockway et al. (1965) and Hofmeyer et al. (1969) in a study with Merino sheep. However, Silva and Starling (2003) described that the heat loss by the respiratory system represented 37 % of total evap- oration in Corriedale sheep. Despite these differences, it is clear that the respiratory tract plays a key role in thermal equi- librium of wool sheep selected for cold environments. On the other hand, cutaneous evaporation seems to be the principal way to eliminate heat loss in animals raised on high ambient temperature as in tropical region. In study with Anglo Nubian goats acclimatized on tropical condition, Maia et al. (2016) described that under ambient temperature close to 30 °C, thermal equilibrium was reached with 15 W m−2, 35 W m−2 and 40 W m−2 of respiratory evaporation, sensible heat loss (radiation and free convection), and cutaneous evap- oration, respectively. Ligeiro et al. (2006) working with cross- bred goats acclimatized to a tropical environment and pure- bred breeds adapted to temperate regions (Saanen and Alpine) observed similar results and greater cutaneous evaporation in crossbred animals. In addition, Maia et al. (2015) reported that under high ambient temperature (36 °C) cutaneous evapora- tion represented over the 93 % of the total heat loss for Bundefined breed goats^ protected from direct and diffuse solar radiation in a semiarid tropical region. According to our findings, the respiratory system is a sig- nificant avenue for latent heat loss under high temperatures, but no prior attempts had been made to identify what the main way for latent heat dissipation in wool-less sheep bred in trop- ical environment. Considering that the heat produced by me- tabolism in Morada Nova sheep can reach 45 W m−2 (Fonsêca, unpublished data) and heat loss (convection and evaporation) by respiratory tract of 9 W m−2 on average, the remaining 36 W m−2 should be lost in the animal surface by sensible way (radiation and convection) and cutaneous evap- oration for the thermal equilibrium to be reached. In semiarid tropical environment, even in shade, during times of greatest thermal load, low thermal gradient (animal coat surface—am- bient) and high radiant mean temperature can impair the elim- ination of heat by sensible avenues, and probably, sweat evap- oration should assume this role. However, more studies are needed to investigate better this hypothesis because some part of the heat can be storage instead dissipated to the environ- ment. The fact is that by the characteristics of the hair coat of Morada Nova sheep, the latent heat can be more easily dissi- pated at the skin than if it were covered by long wool. Conclusions Models presented here allow estimation of heat flow from the respiratory tract in Morada Nova sheep bred in a tropical re- gion and use easily measured physiological and environmen- tal parameters, such as the respiratory rate, ambient air tem- perature, and relative humidity. Acknowledgments We gratefully acknowledge the support of the Animal Science Department and Animal Biometeorology Laboratory of the Universidade Estadual de São Paulo, Brazil. This study was supported by Research Foundation of São Paulo State (FAPESP; Process number: 2011/17388-6). We thank all members of the Behaviour and Biometeorology Group of the Universidade Federal da Paraíba (BIOET) and the members of the Animal Biometeorology Innovation Group of the Universidade Estadual de São Paulo (INOBIO). References Alexander G, Williams D (1962) Temperature regulation in the newborn lamb. VI. Heat exchanges in lambs in a hot environment 13:122– 143 Bennett JW (1973) Regional body surface area of sheep. J Agric Sci 81: 429–432 Berman A (2005) Estimates of heat stress relief needs for Holstein dairy cows. J Anim Sci 83:1377–1384 Brockway JM, McDonald JD, Pullar JD (1965) Evaporative heat loss mechanisms in sheep. J Physiol Lond 179:554–568 Da Silva RG (2008) Biofísica ambiental: Os animais e seu ambiente. São Paulo, Brazil Da Silva RG, Maia ASC (2011) Evaporative cooling and cutaneous sur- face temperature of Holstein cows in tropical conditions. Braz. J Anim Sci 40:1143–1147 Da Silva RG, Maia AS (2013) Principles of animal biometeorology. New York Da Silva RG, Minomo FR (1995) Cicardian and seasonal variation of the body temperature of sheep in a tropical environment. Int J Biometeorol 39:69–73 Da Silva RG, Starling JMC (2003) Cutaneous and respiratory evaporation in sheep under high ambient temperatures. Braz. J Anim Sci 32: 1956–1961 Da Silva RG, La Scala N Jr, Lima Filho AE, Catharin MC (2002) Respiratory heat loss in the sheep: a comprehensive model. Int J Biometeorol 46:136–140 Int J Biometeorol (2017) 61:777–784 783 Da Silva RG, GuilherminoMM,Morais DAEF (2010) Thermal radiation absorbed by dairy cows in pasture. Int J Biometeorol 54:5–11 Da Silva RG, Maia AS, Costa LLM, Queiroz JPAF (2012) Latent heat loss of dairy cows in an equatorial semi-arid environment. Int J Biometeorol 56:927–932 Egito AA,Mariante AS, Albuquerque MSM (2002) Brazilian program of animal genetic resources conservation. Arch Zootec 51:39–52 Fonsêca VFC, Saraiva EP, Pimenta Filho EC, Furtado DA, Mariz TAM, Silva ALS, Almeida MEV, Pinheiro AC (2014) Influence of the climatic environment and duration of labor on the mother- offspring interaction in Morada Nova Sheep. J Anim Sci 92:4123– 4129 Gatenby RM (1986) Exponential relation between sweat rate and skin temperature in hot climates. J Agric Sci 106:175–183 Hofman WF, Riegle GD (1977) Thermo respiratory responses of shorn and unshorn sheep to mild heat stress. Resp Physiol 30:327–338 Hofmeyer HS, Guidry AJ, Waltz FA (1969) Effects of temperature and wool length on surface and respiratory evaporative losses of sheep. J Appl Physiol 26:517–523 Johnson KG (1991) Body temperatures and respiratory rates of free- ranging Merino sheep in and out of shade during summer. Aust J Agric Res 42:1347–1357 Ligeiro EC, Campos Maia AS, daSilva RG, Loureiro CMB (2006) Heat loss by cutaneous evaporation associated with hair coat morpholog- ical characteristics of dairy goats bred in tropical environment. Braz J Anim Sci 35:544–549 Maia AS, da Silva RG, Loureiro CMB (2005) Respiratory heat loss of Holstein cows in a tropical environment. Int J Biometeorol 49:332– 336 Maia ASC, da Silva RG, Nascimento ST, Nascimento CCN, Pedroza HP, Domingos HGT (2015) Thermoregulatory responses of goats in hot environments. Int J Biometeorol 59:1025–1033 Maia ASC, Nascimento ST, Nascimento CCN, Gebremedhin KG (2016) Thermal equilibrium of goats. J Therm Biol 58:43–49 Mitchell B (1972) Relationship between body weight and tidal volume during general anaesthesia in sheep, cattle, pigs and horses. Vet Anaest Analg 3:56–60 Neiva JNM, Teixeira M, Turco SHN, Oliveira SMP, Moura AAAN (2004) Effect of thermal stress on the physiological and performance parameters of sheep maintained in confinement in the northeast of Brazil. Braz J Anim Sci 33:668–678. Oliveira SEO, Costa CCM (2013) Respiratory heat loss in Morada Nova sheep in Brazilian semi-arid regions. J Anim Behav Biometeorolol 1:17–20 Schmidt-Nielsen K (1991) Animal physiology: adaptation and environ- ment, 4th edn. Cambridge University Press, Cambridge Silanikove N (2000) The physiological basis of adaptation in goats to harsh environments. Sma Rum Res 35:181–193 Silva RB (2013) Effect of environment in the physiological variables and methane emission associated to production and heat loss in sheep. Dissertation, Universidade Estadual Paulista Starling JMC, daSilva RG, Cerón-Munhoz M, Barbosa GSSC, Paranhos da Costa MJR (2002) Analyses of some physiological variables for evaluation of adaptation level on sheep submitted to heat stress. Braz J Anim Sci 31:2070–2077 Stevens DG (1981) A model of respiratory vapor loss in Holstein dairy cattle. Trans ASAE 24:151–158 784 Int J Biometeorol (2017) 61:777–784 Models... Abstract Introduction Materials and methods Results and discussion Conclusions References