SciELO - Brazil - Geoestatística para a avaliação do controle ambiental do sistema de ventilação em instalações comerciais para frangos de corte Geoestatística para a avaliação do controle ambiental do sistema de ventilação em instalações comerciais para frangos de corte
Artigos Científicos • Eng. Agríc. 34 (6) • Dez 2014 • https://doi.org/10.1590/S0100-69162014000600004 copiar
Geoestatística para a avaliação do controle ambiental do sistema de ventilação em instalações comerciais para frangos de corte
Geostatistic to evaluete the environmental control in different ventilation systems in broiler houses
Zigomar M. Souza
avicultura; conforto térmico; qualidade do ar; sistema de ventilação
air quality; poultry; thermal welfare; ventilation system
Thayla M. R. de C. CuriI; Rimena do A. VercellinoII; Juliana M. MassariIII; Zigomar M. SouzaIV; Daniella J. de MouraV
IEngª Agrícola, Dra., Faculdade de Engenharia Agrícola, UNICAMP/Campinas – SP, Fone: (19) 98162-3050, thaylamrcarvalho@hotmail.com
IIZootecnista, Mestre, Faculdade de Engenharia Agrícola, UNICAMP/Campinas – SP, rimena.vercellino@gmail.com
IIIMédica Veterinária, Mestre, Faculdade de Engenharia Agrícola, UNICAMP/Campinas – SP, jujumassari@hotmail.com
IVEngº Agrícola, Prof. Dr., Faculdade de Engenharia Agrícola, UNICAMP/Campinas – SP, zigomarms@feagri.unicamp.br
VEngª Agrônoma, Profa. Dra., Faculdade de Engenharia Agrícola, UNICAMP/Campinas – SP, daniella.moura@feagri.unicamp.br
Palavras-chave: avicultura, conforto térmico, qualidade do ar, sistema de ventilação
The ventilation and cooling systems are important to environmental control inside broilers house enabling the bird to express the maximum potential of production. The aim of this study was to evaluate, using geostatistical analysis, the influence of the ventilation system and the different typologies of building in relation to thermal welfare and air quality in broilers house with birds in the final phase. The broiler facilities evaluated were: Blue House – ventilation systems with exhaust fans and nebulizers, Dark House – ventilation systems with exhaust fans and pad cooling and Solid Wall – ventilation systems with exhaust fans and pad cooling. The attributes evaluated were dry bulb temperature, relative humidity, air velocity, concentration of carbon dioxide and ammonia, during the summer time at 2PM on 52 equidistant points inside the building. The data analysis was by statistical and geostatistical. Geostatistical analysis was an efficient tool to evaluate distribution of environmental conditions in the broilers houses and was able to identify the influence of presence or absence of evaporative panel inside conditions.
Keywords: air quality, poultry, thermal welfare, ventilation system
A avicultura de corte tem buscado o progresso por meio de pesquisas nas áreas de genética, instalações, nutrição, manejo, sanidade e conforto ambiental. Todos esses estudos têm como propósito compreender melhor os fatores que influenciam o desenvolvimento e o desempenho de frangos de corte (NAZARENO et al., 2009). Dentre os diversos fatores que influenciam na produção de frangos de corte, os fatores ambientais, como a temperatura, umidade relativa, iluminância, radiação, entre outros, assumem relevante importância no processo de criação dos animais por comprometer a função vital mais importante das aves, a homeotermia (AMARAL et al., 2011).
É sabido que os animais atingem sua produtividade ótima quando são mantidos em ambiente termoneutro, ou seja, quando a energia do alimento não é desviada para compensar desvios térmicos em relação ao intervalo de termoneutralidade para eliminar ou manter seu calor. Neste sentido, o ambiente gerado pelo sistema de ventilação nas instalações de frango de corte desempenha papel fundamental no controle de vários fatores ambientais, diluindo o ar interno por meio do ar fresco externo, uma vez que apresenta função sanitária para a retirada do excesso da umidade, da poeira, dos odores e gases, e para prover oxigênio às aves, e função térmica para a retirada do calor produzido pelas aves, equipamentos, ganho de calor pela radiação solar e iluminação (BUCKLIN et al., 2009; MOSTAFA et al., 2012; CORKERY et al., 2013). Neste sentido, aves na fase final de produção estão em conforto térmico quando a temperatura se apresenta entre 15 e 25 ºC, e a umidade relativa do ar, entre 50 a 70% (NICHOLSON et al., 2004; MELUZZI & SIRRI, 2009). A velocidade do ar é um parâmetro que mais influencia na perda de calor sensível, com valor ideal de 2 m s-1 para frangos de corte até a fase final de crescimento (DOZIER et al., 2006). Em relação à qualidade do ar, algumas instituições, como o GLOBALGAP (2012), estabeleceram critérios para o ambiente aéreo de aviários, sendo o limite máximo de 20 ppm para a concentração de amônia e 5.000 ppm para a concentração de dióxido de carbono.
Ambientes considerados quentes e úmidos, mesmo sem haver mortalidade apresentaram redução em 67% no ganho de peso diário. Por esta razão, é fundamental explorar a eficiência do sistema de ventilação a fim de melhorar todo o projeto de construção e tipologia, melhorando o ambiente interno (BUSTAMANTE et al., 2013). Atualmente, a necessidade do controle efetivo do ambiente interno (térmico e aéreo) nas instalações para criação de frangos de corte é a principal razão do uso de diferentes métodos de simulação para profissionais da área, tais como a geoestatística (BARBOSA FILHO et al., 2009; BORGES et al., 2010; YANAGI JÚNIOR et al., 2011; CARVALHO et al., 2011; CARVALHO et al., 2012; PEREIRA et al., 2012; SILVA et al., 2012; CHOWDHURY et al., 2013; SILVA et al., 2013). A geoestatística é uma ferramenta que permite a interpretação dos resultados baseada na estrutura de sua variabilidade natural, por meio do cálculo da dependência espacial (VIEIRA, 2000). Portanto, o conhecimento da variabilidade espacial dos atributos ambientais e a construção de mapas de krigagem podem auxiliar no manejo correto da ventilação. Assim, o presente trabalho teve como objetivo avaliar e estudar os sistemas de ventilação em relação à qualidade do ambiente térmico e aéreo no interior das instalações comerciais de frangos de corte, por meio da geoestatística.
O presente trabalho foi conduzido em três aviários comerciais, localizados na região de Amparo-SP, com clima Cwa com duas estações bem definidas (conforme classificação de Köeppen).
O registro dos dados foi realizado semanalmente, quando as aves da linhagem Cobb Vantress 500 apresentavam 28; 35 e 42 dias de idade, às 14h, no período de verão do ano de 2012.
O aviário do tipo Blue House (BH), aviário com sistema de ventilação artificial do tipo túnel de pressão negativa com isolamento lateral por cortinas azuis na face interna e prata na face externa, para o controle da luminosidade, estava localizado em Amparo-SP (latitude 22º45'43" sul, longitude 46º47'43" oeste, altitude de 706 m), com orientação cartográfica noroeste-sudeste, dimensões de 17 x 90 x 2,45 m (largura x comprimento x pé-direito) e densidade de 15 aves m-2. O sistema de ventilação era artificial de pressão negativa, contendo dez exaustores, com diâmetro de 1,38 m, hélice com três pás e motor de 735,5 W e vazão nominal de 10,56 m3 s-1 (à pressão de zero Pascal), sem uso de painel evaporativo e vedação com cortinas laterais na cor azul na face interna e prata na face externa.
O aviário do tipo Dark House (DH), aviário com sistema de ventilação artificial do tipo túnel de pressão negativa com isolamento lateral por cortinas pretas na face interna e prata na face externa, para o controle da luminosidade, estava localizado em Amparo-SP (latitude 22º44'01" sul, longitude 46º44'02" oeste, altitude de 706 m), com orientação cartográfica noroeste-sudeste, dimensões de 20 x 150 x 2,90 m (largura x comprimento x pé-direito) e densidade de 12 aves m-2. O sistema de ventilação era artificial de pressão negativa, contendo dezesseis exaustores, com diâmetro de 1,30 m, hélice com três pás e motor de 1.103,25 W e vazão nominal de 10,56 m3 s-1 (à pressão de zero Pascal), com uso de painel evaporativo feito de tijolo cerâmico e vedação com cortinas laterais na cor preta na face interna e prata na face externa.
O aviário do tipo Solid Wall (SW), aviário com sistema de ventilação artificial do tipo túnel de pressão negativa com isolamento lateral por paredes de alvenaria para controle da luminosidade, estava localizado em Amparo-SP (latitude 22º45'37" sul, longitude 46º46'42" oeste, altitude de 706 m), com orientação cartográfica norte-sul, dimensões de 20 x 120 x 3 m (largura x comprimento x pé-direito) e densidade de 13 aves m-2. O sistema de ventilação era artificial de pressão negativa, contendo dezesseis exaustores, com diâmetro de 1,38 m, hélice com três pás e motor de 735,5 W e vazão nominal de 11,42 m3 s-1 (à pressão de zero Pascal), com uso de painel evaporativo feito de celulose e vedação de alvenaria com parede de blocos de concreto.
O registro dos dados de temperatura (Tbs, ºC) e velocidade do ar (Var, m s-1) foi realizado por meio de um anemômetro de fio quente (VelociCalc®, TSITM), com escala de leitura para temperatura de -18 a 93 ºC e resolução 0,1 ºC, e velocidade do ar com amplitude de 0 a 30 m s-1, resolução de velocidade de 0,01 m s-1 e exatidão de ± 0,015 m s-1. Para os dados de umidade relativa (UR, %), foram coletados por meio do equipamento THDL 400 (Instrument®), com amplitude de 25 a 95% e precisão de ±5%. Os dados de concentração de amônia (NH3, ppm) e concentração de dióxido de carbono (CO2, ppm) foram registrados por meio de um medidor instantâneo de concentração de gases (GasAlertMIcro 5, BW Technologies®).
Nos dados de temperatura de bulbo seco e umidade relativa, foram realizadas correções numéricas adaptadas da metodologia proposta por BARBOSA FILHO et al. (2009) que consiste em corrigir a defasagem de tempo entre a coleta do dado no primeiro e no último ponto, em decorrência da coleta manual dos dados.
A variabilidade dos atributos foi verificada baseada na estatística descritiva dos dados registrados (Tbs, UR, Var, CO2 e NH3), obtendo-se valores da média, mediana, coeficiente de variação e o teste de normalidade de Kolmogorov-Smirnov. A análise estatística foi realizada, utilizando o programa computacional estatístico SAS®-.
A dependência espacial foi analisada por meio de ajustes de semivariogramas (VIEIRA, 2000), com base na pressuposição de estacionariedade da hipótese intrínseca, a qual é estimada pela [eq. (1)].
N (h) - número de pares experimentais;
Z(x1) e Z (x1 + h) - observações separadas por uma distância h (m).
O semivariograma é representado no gráfico pela semivariância versus h. Do ajuste de um modelo matemático aos valores calculados da semivariância, são estimados os coeficientes do modelo teórico para o semivariograma (o efeito pepita - C0, m; patamar, C0+C1, m; e o alcance, a, m). Para analisar o grau da dependência espacial dos atributos em estudo, utilizou-se da classificação de CAMBARDELLA et al. (1994), em que são considerados como grau da dependência espacial forte os semivariogramas que têm um efeito pepita < 25% do patamar, moderada quando está entre 25 e 75% e fraca > 75%.
Na determinação da existência ou não da dependência espacial, utilizou-se do exame de semivariogramas, através do programa GS+. A elaboração dos mapas de distribuição espacial das variáveis foi realizada no programa SURFER v.9. Com o objetivo de facilitar a comparação do grau de dependência espacial das variáveis em estudo, o efeito pepita foi expresso como percentagem do patamar (TRANGMAR et al., 1985).
De acordo com o teste de Kolmogorov-Smirnov, as variáveis de Tbs e UR para o aviário SW, Var para os aviários BH e DH, e concentração de NH3 para os aviários BH e SW, aos 28 dias de idade das aves, UR para o aviário DH, Var para os aviários BH e SW, e a concentração de NH3 para o aviário SW, aos 35 dias de idade, e Tbs e Var para o DH, UR para o DH e SW, Var para o BH e DH, e para a concentração de NH3 para os três aviários avaliados, aos 42 dias de idade das aves, apresentaram normalidade (Tabela 1). Com exceção da concentração de CO2, as variáveis avaliadas apresentaram valores próximos de média e mediana, o que permite assumir que as distribuições são suficientemente simétricas. Segundo LITTLE & HILLS (1978), quando o valor da média e da mediana estão próximos, os dados aproximam-se da distribuição normal. Isso pode ser um indicativo de que as medidas de tendência central não são dominadas por valores atípicos na distribuição (CAMBARDELLA et al., 1994).
Em relação à variabilidade dos dados ambientais classificada em baixa (coeficiente de variação abaixo de 12%), média (coeficiente de variação entre 12 e 24%) e alta (coeficiente de variação acima de 24%), de acordo com a literatura (WARRICK & NIELSEN, 1980), foi possível observar que, aos 28 dias de idade das aves, a Tbs para os três aviários avaliados e para a UR no SW, aos 35 e aos 42 dias de idade, as variáveis de Tbs e UR para os três aviários avaliados, ocorreu baixa variabilidade dos dados. A média variabilidade dos dados foi observada para a UR nos aviário BH e DH aos 28 dias de idade e para a concentração de NH3 no aviário DH aos 42 dias de idade das aves. A Var, concentração de CO2 e NH3 para os aviários BH, DH e SW apresentaram alta variabilidade aos 28 e 42 dias de idade das aves, da mesma forma que a Var e concentração de CO2 para os três aviários e concentração de NH3 para o BH e SW aos 35 dias de idade das aves. A avaliação da homogeneidade dos parâmetros climáticos é essencial para assegurar a qualidade do manejo do ambiente térmico e aéreo, a fim de minimizar a ocorrência de regiões quentes e com pouca renovação de ar (pontos mortos), melhorando a qualidade e a produção do lote (BOURNET & BOULARD, 2010).
Aos 28 dias de idade das aves, os semivariogramas da Tbs e da concentração de NH3 no aviário BH e da UR no BH e SW foram ajustados no modelo gaussiano, o modelo esférico ajustou- -se para a Tbs nos aviários DH e SW, para a Var nos três aviários avaliados, concentração de CO2 e NH3 para SW e DH, respectivamente. A UR para DH, a concentração de CO2 para os aviários BH e DH e a concentração de NH3 para o SW tiveram seus semivariogramas ajustados no modelo exponencial. A Tbs e a UR para os aviários BH e DH, e a concentração NH3 para o BH foram ajustadas no modelo gaussiano aos 35 dias de idade das aves, assim como os semivariogramas para a Tbs e a UR em SW,Var para BH e SW e concentração de NH3 para DH foram ajustados no modelo esférico, e no modelo exponencial para a Var em DH, concentração de CO2 e NH3 para os três aviários e para SW, respectivamente. Os semivariogramas para aves com 42 dias de idades foram ajustados em gaussiano para a Tbs em SW, UR para BH, DH e SW, concentração de CO2 no SW e concentração de NH3 no SW, esférico para Tbs no DH e Var para os três aviários, e exponencial para a Tbs em Bh e para a concentração de CO2 e NH3 para o aviário DH.
A análise da relação C0/(C0+C1) mostrou que as variáveis apresentaram grau de dependência espacial forte e/ou moderada (Tabela 2). Quanto maior a dependência espacial, menor será a contribuição do efeito pepita na variabilidade dos dados, consequentemente melhores serão as estimativas na krigagem. Portanto, as distribuições dos atributos microclimáticos e aéreos no espaço não são aleatórias, uma vez que todos apresentaram valores moderados ou fortes para o grau de dependência espacial (CARVALHO et al., 2012).
O alcance representa a distância limite de amostragem em que os pontos coletados estão correlacionados, contribuindo para o planejamento amostral e a definição dos locais onde as variáveis climáticas serão amostradas (McBRATNEY & WEBSTER, 1983). Os dados apresentaram valores de alcance variando de 12 a 90 m.
O valor médio da temperatura esteve além do recomendado para os aviários estudados, ou seja, acima do intervalo entre 19,1 e 20,2 ºC, capaz de melhorar as características de desempenho zootécnico (ARAUJO et al., 2011). O aviário DH apresentou os menores valores de Tbs, indicando que o sistema de resfriamento foi mais eficaz comparado com os outros aviários (FIGURA 2). O aviário SW, com 28 dias de idade, apresentou áreas com temperaturas maiores (26 ºC) e áreas com temperaturas menores (24 ºC), remetendo a necessidade de maior controle do sistema de ventilação para manter o ambiente mais homogêneo.
O aviário BH apresentou os maiores valores de Tbs com maior potencial de perda de calor por meios evaporativos, pelos valores de UR se apresentarem próximos de 50%, sendo que a literatura recomenda entre 50 e 70% (COBB, 2013). Foi possível observar ilhas do centro dos aviários com valores de Tbs mais amenas em função do sistema de resfriamento por painel evaporativo no DH e SW (SILVA et al., 2013), e os maiores valores estavam próximos às entradas de ar através do painel evaporativo (FIGURA 3). Sendo que os aviários DH e SW apresentaram valores de UR acima do recomendado em função do sistema de ventilação tipo túnel, que tende a apresentar maiores níveis de UR em relação ao recomendado (CARVALHO et al., 2012). A Tbs alta, associada à UR do ar alta, dificulta as perdas de calor do animal (ROCHA et al., 2010).
Os valores registrados da Var nos aviários BH e DH apresentaram-se no intervalo considerado ideal, ou seja, próximo de 1,5 m s-1 (MEDEIROS et al., 2005), e os menores valores foram observados no aviário SW, indicando necessidade de intensificação do sistema de ventilação. Os maiores valores da velocidade do ar foram observados próximo às entradas de ar, aos 28 dias de idade das aves no aviário BH e aos 42 dias de idade das aves nos aviário DH e SW (FIGURA 4).
A concentração de CO2 esteve abaixo do limite máximo considerado ideal para todos os aviários, abaixo de 3.000 ppm (GLOBALGAP, 2012), indicando que o sistema de ventilação utilizado pelos aviários estudados foi efetivo para a renovação do ar. Ainda é possível observar maiores valores de CO2 e NH3 próximo aos exaustores, em função do sistema de ventilação tipo túnel de pressão negativa, que faz o encaminhamento do ar no sentido longitudinal do aviário – FIGURAS 5 e 6 (CARVALHO et al., 2012; SILVA et al., 2013). Os níveis de dióxido de carbono podem ser utilizados para a definição das taxas de ventilação necessárias em aviários de frango de corte (CALVET, 2011).
A concentração de NH3 esteve abaixo do limite máximo considerado ideal para todos os aviários, abaixo de 20 ppm (GLOBALGAP, 2012), embora os aviários SW, aos 42 dias de idade das aves, estivesse acima do limite máximo de 10 ppm preconizado pelo manual de frango de corte da linhagem (COBB, 2013). Os mapas de distribuição da concentração de NH3 apresentaram níveis mais altos próximos aos exaustores, de forma análoga aos mapas de concentração de CO2 (FIGURA 6).
Por meio da análise geoestatística, observou-se que o sistema de ventilação utilizado nos diferentes aviários foi capaz de manter a qualidade do ar em relação à concentração de amônia e dióxido de carbono, embora não suficiente para a manutenção da ambiência em condições adequadas à criação das aves. A temperatura e a umidade relativa apresentaram-se além do adequado nos aviários estudados, com exceção do aviário BH, que apresentou umidade aquém. A velocidade do ar esteve em níveis adequados para os aviários BH e DH, e aquém para o SW, remetendo que a presença ou a ausência do painel evaporativo influencia na condição-padrão de funcionamento do sistema de ventilação. O sistema de ventilação com sistema de resfriamento por painel evaporativo auxilia no resfriamento do ar no centro do aviário, enquanto nas paredes laterais ocorrem temperaturas acima do recomendado pela literatura.
Recebido pelo Conselho Editorial em: 02-7-2014
Aprovado pelo Conselho Editorial em: 24-9-2014
AMARAL, A.G.; YANAGI JUNIOR, T.; LIMA, R.R.; TEIXEIRA, V.H.; SCHIASSI, L. Effect of the production environment on sexed broilers reared in a commercial house. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, Belo Horizonte, v.63, n.3, p.649-658, jun. 2011.
ARAUJO, R.B.; FARIA, D.E.; FARIA FILHO, D.E.; LIMA, C.G.; TREVISAN, R.B.; SOUZA, K.M.R.; SAKAMOTO, M.I.; SOUZA, V.N. Modelos de superfície de resposta para predição do desempenho de frangos e elaboração de análise econômica. Revista Brasileira de Saúde e Produção Animal, Salvador, v.12, n.3, p.770-783, 2011.
BARBOSA FILHO, J.A.D.; VIEIRAI, F.M.C.; SILVA, I.J.O.; GARCIA, D.B.; SILVA, M.A.N.; FONSECA, B.H.F. Transporte de frangos: caracterização do microclima na carga durante o inverno. Revista Brasileira de Zootecnia, Viçosa, MG, v.38, n.12, dez. 2009.
BORGES, G.; MIRANDA, K.O.S.; RODRIGUES, V.C.; RISI, N. Uso da geoestatística para avaliar a captação automática dos níveis de pressão sonora em instalações de creche para suínos. Engenharia Agrícola, Jaboticabal, v.30, n.3, p.377-385, 2010.
BOURNET, P.E.; BOULARD, T. Effect of ventilator configuration on the distributed climate of greenhouses: A review of experimental and CFD studies. Computers and Electronics in Agriculture, Amsterdam, n.74, p.195217, 2010.
BUCKLIN, R.A.; JACOB, J.P.; MATHER, F.B.; LEARY, J.D.; NÄÄS, I.A. Tunnel ventilation of broiler houses Gainesville: Institute of food and Agricultural Sciences, 2009. (FACTSHEET PS 46).
BUSTAMANTE, E.; GARCÍA-DIEGO, F.J.; CALVET, S.; ESTELLÉS, F.; BELTRÁN, P.; HOSPITALER, A.; TORRES, A.G. Exploring ventilation efficiency in poultry buildings: The validation of computational fluid dynamics (CFD) in a cross-mechanically ventilated broiler farm. Energies, Singapura, v.6, p.2605-2623, abril 2013.
CALVET, S.; ESTELLÉS, F.; CAMBRA-LÓPEZ, M.; TORRES, A.G.; VAN DEN WEGHE, H.F. The influence of broiler activity, growth rate, and litter on carbon dioxide balances for the determination of ventilation flow rates in broiler production. Poultry Science, Oxford, v.90, p.24492458, 2011.
CAMBARDELLA, C.A.; MOORMAN, T.B.; NOVAK, J.M. et al. Field scale variability of soil properties in Central Iowa soils. Soil Science Society of America Journal, Madison, v.58, n.5, p.1501-1511, 1994.
CARVALHO, T.M.R.; MOURA, D.J.; SOUZA, Z.M.; SOUZA, G.S.; BUENO, L.G.F. Litter and air quality in different broiler housing conditions. Pesquisa Agropecuária Brasileira, Brasília, v.46, n.4, p. 351-361, 2011.
CARVALHO, T.M.R., MOURA, D.J.; SOUZA, Z.M.; SOUZA, G.S.; BUENO, L.G.F.; LIMA, K.A.O. Use of geostatistics on broiler production for evaluation of different minimum ventilation systems during brooding phase. Revista Brasileira Zootecnia, Viçosa, MG, v.41, n.1, p.194-202, jan.2012.
CHOWDHURY, S.; THEMUDO, G.E.; SANDBERG, M.; ERSBØLL, A.K. Spatio-temporal patterns of Campylobacter colonization in Danish broilers. Epidemiology and Infection, Cambridge, v.141, n.5, p.997-1008, 2013.
COBB. Suplemento de crescimento e nutrição para frangos de corte Massachusetts, 2013. 14 p. Disponível em: <http://www.cobb-vantress.com/docs/default-source/cobb-700-guides/cobb700-broiler-performance-and-nutrition-supplement---portuguese.pdf?sfvrsn=0>. Acesso em: 15 jul. 2014.
CORKERY, G.; WARD, S.; KENNY, C.; HEMMINGWAY, P. Monitoring Environmental Parameters in Poultry Production Facilities. Computer Aided Process Engineering, Graz, v. 1, p.1-10, 2013. Disponível em <http://hdl.handle.net/10197/4257>. Acesso em: 24 nov. de 2013.
DOZIER, W.A.; PURSWELL, J.L.; BRANTON, S.L. Growth responses of male broilers subjected to high air velocity for either twelve or twenty-four hours from thirty-seven to fifty-one days of age. The journal of applied poultry research, Oxford, v.15, p.362-366, 2006.
GLOBALG.A.P. Integrated farm assurance / poultry. Controls Points and compliance criteria. Cologne, 2012. 54 p. Disponível em <www.globalgap.org>. Acesso em: 24 fev. de 2012.
LITTLE, T.M.; HILLS, F.J. Agricultural experimentation: design and analysis. New York: John Wiley & Sons, 1978. 368p.
McBRATNEY, A.B.; WEBSTER, R. How many observations are needed for regional estimation of soil properties? Soil Science, Philadelphia, v.135, n.3, p.177-183, 1983.
MEDEIROS, C.M.; BAÊTA, F.C.; OLIVEIRA, R.F.M.; TINÔCO, I.F.F.; ALBINO, L.F.T.; CECON, P.R. Índice térmico ambiental de produtividade para frangos de corte. Revista Brasileira de Engenharia Agrícola e Ambiental, Campina Grande, v.9, n.4, p. 660-665, 2005.
MELUZZI, A.; SIRRI, F. Welfare of broiler chickens. Italy Journal of Animal Science, Pavia, v.8, n.1, p.161-173, mar. 2009.
MOSTAFA, E. et al. Computational fluid dynamics simulationof air temperature distribution inside broiler building fitted with duct ventilation system. Biosystems Engineering, Queensland, v.112, p.293-303, 2012.
NAZARENO, A.C.; PANDORFI, H.; ALMEIDA, G.L.P.; GIONGO, P.R.; PEDROSA, E.M.R.; GUISELINI, C. Avaliação do conforto térmico e desempenho de frangos de corte sob regime de criação diferenciado. Revista Brasileira de Engenharia Agrícola e Ambiental, Campina Grande, v.13, n.6, p.802-808, 2009.
NICHOLSON, F. A.; CHAMBERS, B. J.; WALKER, A. W. Ammonia Emissions from Broiler Litter and Laying Hen Manure Management Systems. Biosystems Engineering, Queensland, v.89, n.2, p.175185, 2004.
PEREIRA, P.A.; YANAGI JUNIOR, T.; SILVA, J.P.; LIMA, R.R.; CAMPOS, A.T.; ABREU, L.H.P. Technical evaluation of artificial ligthing systems for broiler houses. Engenharia Agrícola, Jaboticabal, v.32, n.6, p.1011-1024, nov./dez. 2012.
ROCHA, H.P.; FURTADO, D.A.; NASCIMENTO, J.W.B.; SILVA, J.H.V. Índices bioclimáticos e produtivos em diferentes galpões avícolas no semiárido paraibano. Revista Brasileira de Engenharia Agrícola e Ambiental, Campina Grande, v.14, n.12, p.1330-1336, 2010.
SILVA, E.G.; SANTOS, A,C.; FERREIRA, C.L.S.; SOUSA, J.P.L.; ROCHA, J.M.L.; SILVEIRA JÚNIOR, O. Variabilidade espacial das características ambientais e peso de frangos de corte em galpão de ventilação negativa. Revista Brasileira de Saúde e Produção Animal, Salvador, v.14, n.1, 2013.
SAS INSTITUTE. Statistical analysis system for Windows: computer program manual. Cary, 1995. 705p.
SILVA, I. M.; PANDORFI, H.; ALMEIDA, G. L. P.; GUISELINI, C.; CALDAS, A. M.; JACOB, A. L. Análise espacial das condições térmicas do ambiente pré-ordenha de bovinos leiteiros sob regimes de climatização. Revista Brasileira de Engenharia Agrícola e Ambiental, Campina Grande, v.16, n.8, 2012.
TRANGMAR, B. B.; YOST, R. S.; UEHARA, G. Application of geostatistics to spatial studies of soil properties. Advances in Agronomy, Newark, v.38, n.1, p.45-93, 1985.
VIEIRA, S.R. Geoestatística em estudos de variabilidade espacial do solo. In: NOVAIS, R.F.; ALVARES, V.H.; SCHAEFER, C.E.G.R. Tópicos em ciência do solo Viçosa: Sociedade Brasileira de Ciência do Solo, 2000. v.1, p.1-53.
WARRICK, A.W.; NIELSEN, D.R. Spatial variability of soil physical properties in the field. In: HILLEL, D. Applications of soil physics New York: Academic,1980. p.319-344, 1980.
YANAGI JÚNIOR, T.; AMARA, A.G.; TEIXEIRA, V.H.; LIMA, R.R. Spatial characterization of thermoacoustic and luminance environment in commercial broiler chicken house. Engenharia Agrícola, Jaboticabal, v.31, n.1, p.1-12, jan./feb. 2011.