Geostatistics applied to swine facilities equipped with evaporative cooling system

The purpose of this study was to evaluate the spatial variation of temperature, relative humidity and black globe humidity temperature (BGHT) index in growing-finishing swine facilities, with and without evaporative air cooler. The experiment was conducted in a 60 x 9 m facility containing 32 pens of 3.20 x 4.00 m each, divided into two sectors. The treatments were without evaporative air cooler (control) and with evaporative air cooler. Temperature and relative humidity data were collected weekly in two periods (morning and afternoon) at 32 different sampling points, encompassing the center of each pen and the outer perimeter of the facility. Animal performance variables were: feed intake, feed conversion, weight out of the nursery, slaughter weight, weight gain and daily weight gain. According to the results, the analyzed variables showed spatial dependence, with lower temperature and higher relative humidity in the evaporative air cooler environment. However, evaporative air cooler system provided better environment thermal conditions inside the facilities, lowering BGHT values, which were lower in the adiabatic cooling sector. Although thermal conditions were better in the air cooler environment, swine performance was not affected. globo negro e umidade (ITGU) em instalações de suínos na fase de crescimento e terminação, com e sem resfriador de ar evaporativo. O experimento foi conduzido em uma edificação de 60 x 9 m, contendo 32 baias de 3,20 x 4,00 m cada uma, divididos em dois setores. Os tratamentos foram sem resfriador de ar evaporativo (controle) e com resfriador de ar evaporativo sendo os dados de temperatura e umidade relativa coletados semanalmente em dois períodos (manhã e tarde) em 32 pontos de amostragem diferentes abrangendo o centro de cada baia e o perímetro exterior das instalações. As variáveis de desempenho dos animais foram: consumo de ração, conversão alimentar, peso ao sair da creche, peso ao abate, ganho de peso e ganho de peso diário. Os resultados mostraram que as variáveis analisadas apresentaram dependência espacial com temperatura mais baixa e umidade do ar mais alta no setor com resfriador evaporativo. Desta forma, o resfriamento evaporativo melhorou as condições térmicas ambientais no interior da instalação, reduzindo os valores do ITGU, os quais foram menores no setor de resfriamento adiabático. Embora o conforto térmico do ambiente com resfriamento evaporativo tenha sido melhor, os suínos responderam a esta melhoria, em termos de desempenho.


Introduction
From the 70s, pig production began to be held in a confinement system aiming to improve the sanitary control, reduce energy loss from the animal and increase productivity (Sampaio et al., 2007). Under heat stress, the production efficiency, reproduction, feed conversion, health and welfare of pigs can be severely affected. Physiological and behavioral responses and thus swine production in tropical conditions can be favored by adopting evaporative cooling, as water evaporation to the atmosphere is usually an effective solution to alleviate heat stress (Haeussermann et al., 2007).
Air temperature and relative humidity (RH) inside facilities influence animal health and performance. By using evaporative cooling system, environment thermal conditions may improve (Ponciano et al., 2013;Liang et al., 2014). Sartor et al. (2003) reported a reduction in the black globe humidity index inside the facility from 83.5 to 82.4 during the critical period of the day (2 p.m.), besides reduction in feed conversion from 3.27 to 2.88, and increase in weight gain, from 0.97 to 1.02 kg a day.
Currently, traditional statistical analysis based on means has been replaced by geostatistical analyses, in which semivariograms consider the spatial dependence between observed samples. According to Vieira (2000), when the studied attribute varies from a place to another, geostatistics is a tool that enables the spatial vision of this behavior based on the spatial dependence. The geostatistics concept has been applied in different fields as a tool to study the variability of environment attributes. Also in breeding facilities, it is of great importance in interpreting and understanding results, based on the structure of its natural variability, mainly the interaction between physical attributes of the environment (Wang et al., 2002).
In this context, this research aimed to evaluate the spatial variability of air temperature and relative humidity (RH) in growingfinishing swine facilities with and without air cooler system, besides evaluating swine performance under these conditions.

Material and Methods
The experiment was developed in a farm destined for growing-finishing swine production, at Santo Antônio district, Concórdia, Santa Catarina state, Brazil, during summer period in the Southern hemisphere from November to February.
The facility area was 60 x 9 m, covering 32 pens measuring 3.20 x 4.00 m, central corridor of 1 m, ceiling height of 2.5 m, divided into two sectors separated by 50 micra black canvas curtains, containing one treatment each. The facility was opened with natural thermal conditioning and management of yellow polyethylene curtain. At ceiling height, liner was placed with the same material as the curtain. The roof consisted of 6-mm fiber cement and masonry gables. These treatments consisted of: environment without evaporative air cooler (control) and with evaporative air cooler. Sixteen pens housed control animals, submitted to conventional system, an open facility in natural thermal conditions without any mechanical ventilation and cooling system. The remaining sixteen pens were equipped with evaporative air cooler system. In this system, four Joape 777 coolers were used, two at each row of pen, distributed longitudinally at the center of the pen at 1.5 m height ( Figure  1). These coolers had the following characteristics: 18,000 m³ air flow, 74 dB noise level, 430 W consumption and 1.0 x 1.07 x 0.37 m dimensions; and were activated by automatic thermostat (± 0.1 °C accuracy), adjusted to turn on and off when temperature reached 27.0 °C and 26.9 °C, respectively.
Thirty-two sampling points were examined in the experiment, according to Figure 1; half of sampling was performed around the facility at the highest point of each pen, and the other half measurements were taken inside the facility at the center of each pen for both systems. At each sampling point, temperature and RH were measured using a multifunction device (TESTO  ), weekly at 9 a.m. and 2 p.m. throughout the 16 experimental weeks, totalizing 1,024 samples each treatment.
A total of 352 animals, half male and half female, from industrial Agroceres PIC crossbreeding, were utilized. Each pen housed 11 animals aged 58 days and with 23 kg of mean weight. Animal performance was evaluated through the following variables: feed intake, feed conversion, weight out of the nursery, slaughter weight, weight gain, and daily weight gain in growing-finishing periods (112 days). These variables were evaluated through analyses of variance using the GLM procedure of SAS, taking into account the effects of sex, block (initial weight and pen location) within the sex, treatment and interaction between treatment and sex.
The black globe humidity temperature (BGHT) index is the most precise indicative of thermal comfort in animal production according to Buffington et al. (1981), who proposed BGHT calculation through data collected at each point, following Eq. 1: and from 52.55 to 58.77% in the afternoon. BGHT index in the morning varied between 72.91 and 74.97, and in the afternoon, between 76.39 and 78.64. Throughout the experimental time frame, the outside conditions were: 34.6 °C the maximum temperature average and 38% the minimum RH average.
In order to allow scrutiny of variables, histograms were built to visualize the variable behavior in relation to data concentration trend (symmetric or asymmetric). This trend especially for the non-spatial analysis of data can drive different analytical procedures. In this way, data distributions regarding temperature, RH and BGHT index in the morning and afternoon have not obeyed a normal distribution (Figure 2).
Semivariograms are illustrated in Figure 3, while the parameters used in the analysis are shown in Table 2. There was a strong spatial dependence of data in both studied periods, with the best fit spherical model represented by the high value of the coefficient of determination (R 2 ) and low value of RSS.
The spherical model is one of the most commonly used models in geostatistical analysis according to results presented by Campos et al. (2009)  According to temperature and RH values in the morning and afternoon at each point, the spatial analysis of data was carried out using GS + software. Data were submitted to geostatistical analysis, where the model was chosen based on the smallest residual sum of squares (RSS), identifying parameters as the nugget effect (Co), sill (Co + C 1 ) and range (Ao). In order to verify the existence of spatial dependence, the spatial dependence index (SDI), which is the ratio representing the percentage of data variability explained by spatial dependence, was adapted from Kitamura et al. (2007) according to Eq. 3. It may be classified as weak (SDI ≤ 25%), moderate (25% < SDI < 75%) and strong (SDI ≥ 75%).

Parameters
Temperature morning where: SDI -evaluator index of the spatial dependence ratio; and, C 1 -spatial variance; and Co + C 1 -total variance (sill) for the semivariogram.

Results and Discussion
In total 1,024 recordings were obtained during the experimental period (Table 1). The temperature varied from 24.0 to 25.9 °C in the morning and from 27.3 to 29.3 °C in the afternoon. The RH ranged from 62.04 to 68.82% in the morning    Table 2. Parameters of the theoretical semivariogram from temperature and humidity in the morning and afternoon periods via ordinary kriging of temperature through the facility in both treatments was observed ( Figure 4A). As already expected, higher values of temperature in the afternoon were accompanied by lower RH, unlike observed in the morning (Figures 4A and B). In an environment with evaporative cooler, air temperature was lower and more homogeneous in both day periods than in the environment without cooler. In the afternoon the difference between systems was only 2 °C.
Lower temperature and higher RH values in the environment with coolers were observed close to the place where coolers were set. These were effective in maintaining the environment temperature pre-set in the thermostat. However, this was not homogeneous throughout the facility. Air RH was greater in the environment with cooler either in the morning or afternoon ( Figure 4B).
As the experiment was conducted in the Southern hemisphere, the northern side of the facility presented higher values of temperature and lower values of RH in comparison to the southern side at both periods. This condition was due to solar output along the day on the facility oriented in the east-west direction.
According to Haeussermann et al. (2007), evaporative coolers are efficient to relieve heat stress in swine. However, the effect of increasing RH by evaporative cooling in closed facilities is still controversial. These same authors concluded that the evaporative fraction reached 100% during steady state, and 63% under warm and dry environment conditions (28 °C and 53% RH). Under moderate conditions (21 °C and 69%) the evaporative fraction reached 89% and under cold conditions (13 °C and 83%), the evaporative fraction reached 65%.
As shown in Figure 4B, the maximum RH value was 58% in the environment with cooler in the afternoon, which is the most stressful period for animals, evidencing that the air still had evaporative capacity. Sartor et al. (2003) reported that evaporative cooling system improved thermal conditions inside the facilities for finishing swine, by reducing BGHT index, as well as providing better feed conversion and greater weight gain of animals.
BGHT index values characterized as thermoneutral environment by Orlando et al. (2006) were close to 69.1 and 68.7 for swine kept under thermal comfort. As observed in Figure 5, these values were not achieved in the environment with and without evaporative cooler, showing that, despite the temperature decreases and RH increases using coolers, they were not effective in promoting animal thermal comfort.
No significant effect (p > 0.05) was observed for evaporative cooling system regarding feed intake, feed conversion, weight gain out of the nursery, slaughter weight, weight gain and daily weight gain (Table 3).
Evaporative cooling system improved physiological, behavioral and productive responses from growing swine (Huynh et al., 2006), environmental thermal condition (Sartor et al., 2003) and provided an increase in weaning piglet weight (Justino et al., 2015). Otherwise, Collin et al. (2001) reported that weight gain and voluntary feed intake of swine were reduced by up to 37 and 30% respectively at 33 °C when compared to 23 ºC.
Nevertheless, little is known about the benefits of evaporative cooling system at finishing phase, as heat stress relief, as well as management procedures required by this system, once the efficiency of the evaporative cooling is directly influenced by local microclimatic conditions, and by animal physiological and behavioral responses.

Conclusions
1. Higher RH values and lower temperatures were found with evaporative cooler, providing better environment conditions for pigs compared with the environment without cooler.
2. The environment with evaporative cooler showed that temperature and RH conditions are prone to increase the equipment efficiency.
3. The spatialization of temperature, RH and BGHT index provided more detailed information on the operation of the evaporative cooler, although this equipment does not provide thermal comfort for swine.
4. Even modifying thermal conditions using evaporative air cooler, it was not enough to enhance pig performance; also, daily weight gain did not show significant difference (p > 0.05) between sex.
5. By psychometric air conditions presented within the facility, it is urged to adopt more air coolers or replace them by other evaporative cooling system, so that animals respond to environment change. Table 3. Means and standard errors of performance variables according to sex and treatment (cooling system)