Avaliação de risco à produção de frango de corte do Estado de São Paulo em função da temperatura ambiente

Risk evaluation of broiler production in the state of São Paulo as function of ambient temperature

ARTIGOS CIENTÍFICOS

CONSTRUÇÕES RURAIS E AMBIÊNCIA

Douglas D. Salgado^I; Irenilza de A. Nääs^II

^IEstatístico, Doutor em Engenharia Agrícola, [e-mail]

^IIEngª Civil, Professora Colaboradora, FEAGRI-UNICAMP, [e-mail]

RESUMO

As perdas produtivas na avicultura de corte, provenientes de climas com temperaturas diárias elevadas, são potencialmente de grande magnitude, pois abrangem perdas diretas e indiretas. Esta pesquisa teve a finalidade de apontar os municípios localizados onde há mais risco de temperaturas extremas diárias. Os dados históricos meteorológicos de temperaturas diárias, máximas e mínimas, fornecidos pelo CEPAGRI-UNICAMP, foram adaptados para análise estatística, descritiva e exploratória. As temperaturas diárias foram categorizadas em forma binária, estimando-se probabilidades condicionais (riscos) por meio das frequências relativas, obtidas em tabelas de contingência. Utilizando o software Minitab 15^®, foi realizada análise descritiva, análise de risco e, por fim, avaliação de associação. Este trabalho apontou que os municípios da região oeste do Estado de São Paulo foram os mais suscetíveis a apresentarem perda na produção avícola devido às temperaturas ambientais, sendo recomendado maior cuidado com o excesso de calor nos alojamentos das granjas. Também foi constatado que os valores médios e medianos das temperaturas mínimas são bons preditores do risco, devido à alta associação entre o risco e essas variáveis.

PAlavras-chave: bem-estar animal, temperaturas extremas, risco.

ABSTRACT

Broiler production productive losses due to weather with daily high temperature are potentially of great direct and indirect productive losses. This research aimed to develop a climatologic profile of counties of the State of São Paulo, under the conditions that involved risk of losses in broiler production. A set of historical data given by the CEPAGRI-UNICAMP were adapted to the statistical analysis, and furthermore the descriptive and exploratory analysis were performed with the temperatures. Daily temperature data were categorized in binary form to estimate the conditional probabilities (risk) through the relative frequencies, obtained from tables of contingency. Using the software Minitab 15^® both descriptive and risk analysis were performed, and at last a connection evaluation. These researches showed that the counties from the West region of the state of São Paulo are the most susceptible to present loss in poultry production due to ambient temperatures. For these counties it is recommended high concern from the producer with the excess of heat inside the houses in the farms. It was also found that mean and median values of minimum temperature are good predictors of risk due to the high association between risk and these variables.

Key words: : animal welfare, extreme temperatures, risk.

INTRODUÇÃO

O desenvolvimento tecnológico brasileiro, voltado à produção zootécnica, vem aperfeiçoando-se ao longo do tempo, equiparando-se aos grandes produtores internacionais, aumentando a produção para exportação, sob os elevados níveis de qualidade exigidos no mercado externo (SILVA, 2000; BRUNOZILLI & SOUZA, 2004; ROPPA, 2004). Para tanto, justificou-se a otimização do sistema de produção, automação e utilização de recursos de zootecnia de precisão (NÄÄS & CURTO, 2001). Contudo, faz-se também necessária a redução de perdas produtivas, principalmente nas fases finais de criação, sobretudo quando provocada pelo estresse térmico, reduzindo os índices zootécnicos (HAHN, 1982; YAHAV et al., 1995; COOPER & WASHBURN, 1998; TANKSON et al., 2001). Grande parte do equilíbrio energético na granja de produção de frangos de corte está na manutenção do ambiente interno dentro da faixa de conforto, o que ocorre em função da temperatura ambiente externa ao aviário, além da idade das aves (MACARI & GONZALES, 1990; DEATON et al., 1997; MOURA, 2001). Sendo a avicultura comercial uma atividade intensiva, as condições adversas de clima podem em parte ser eliminadas ou atenuadas, quando os galpões são localizados em regiões mais favoráveis (ARADAS et al., 2005; MIRAGLIOTTA et al., 2006; VALE et al., 2008).

A temperatura ideal para criação de frangos de corte está próxima aos 20 ºC (MACARI & GONZALES, 1990; COOPER & WASHBURN, 1998; TANKSON et al., 2001), a temperatura média da superfície da pele é ao redor de 33 ºC, e a temperatura interna, de 41 ºC; assim, como existe diferença de temperatura entre a superfície da pele e o ambiente, ocorrerá fluxo de energia (calor) do local mais quente (pele - 33 ºC) para o mais frio (ambiente - 25 a 28 ºC), promovendo perda de calor sensível ao ambiente. Estudos dos efeitos de variadas condições térmicas sobre as aves, por meios experimentais controlados, encontraram faixas de alto risco de perda produtiva atingindo até a mortalidade, decorrente, direta e indiretamente, da exposição ao estresse calórico (TANKSON et al., 2001; NIENABER & HAHN, 2004; PEREIRA, 2005; SALGADO, 2006; VALE et al., 2008).

As análises, descritiva e exploratória, são importantes meios de pré-avaliar os dados de clima, pois fornecem subsídios para o conhecimento mais concreto a respeito das variáveis registradas, o que permite o levantamento de diversas hipóteses importantes para a investigação científica. De maneira geral essas hipóteses podem ser avaliadas por meio de análises confirmatórias, como testes de hipóteses, modelos, etc. (SALGADO, 2006). Estudos de zoneamento de risco têm sido utilizados na agricultura e, na área de climatologia, a aplicação de modelos de predição tem sido dificultada pelas múltiplas interações de fatores ambientais naturais com abruptas variações não lineares, aumentando-se a incerteza sobre a previsibilidade de eventos extremos (JONES, 2003; SCHNEIDER, 2004).

O Brasil possui um rebanho de aproximadamente 822 milhões de aves produtoras de carnes e ovos (IBGE, 2006), cerca de nove milhões de toneladas produzidas por ano, o que faz do País o maior exportador mundial desse alimento. A região Sudeste é a segunda principal região produtora no País, com 28% de aves alojadas, sobretudo no Estado de São Paulo, com 18% (UBA, 2008). Nesse cenário, esta pesquisa teve a finalidade de elaborar uma descrição das temperaturas máximas e mínimas em 38 municípios do Estado de São Paulo, avaliar o risco de temperaturas diárias extremas (MACHADO, 1950; SARTORI, 2000 e ROSSATO et al., 2003), apontando-se as localidades que tendem a apresentar condições que são mais propícias às perdas na produção avícola, por excesso de calor, além de associar as medidas estatísticas de descrição ao risco de temperatura extrema.

MATERIAL E MÉTODOS

Os dados iniciais foram cedidos pelo Centro de Pesquisas Meteorológicas e Climáticas aplicadas à Agricultura da Universidade Estadual de Campinas (CEPAGRI, 2008). Os dados consistiam em temperaturas diárias (temperatura de bulbo seco) máximas e mínimas, registradas por um período de até 17 anos, em 40 estações meteorológicas distribuídas pelo Estado de São Paulo.

Preparação dos dados históricos e classificação das variáveis

Os dados foram dispostos da seguinte forma na planilha: 1) Nas diferentes colunas, estavam presentes as variáveis de controle - local (V1), dia (V2), mês (V3) e ano (V4), e as variáveis principais - temperatura máxima diária (V5) e temperatura mínima diária (V6), e 2). As linhas representam os registros diários consecutivos. Desse modo, a matriz dos dados foi composta por 86 mil registros diários (linhas) e cinco variáveis (colunas). Com base em MACHADO (1950), SARTORI (2000) e ROSSATO et al. (2003), foi elaborada uma variável de classificação binária, denominada "Temperatura diária extrema" (V7). A classificação foi estabelecida a partir das temperaturas máximas e mínimas diárias, simultaneamente, da seguinte maneira, com organização binária: a) Quando a temperatura diária máxima estava acima de 32° C e a mínima acima de 22° C, então a temperatura diária era considerada extrema e recebia o valor "um" (1), e 2) Caso contrário, a temperatura diária era não alarmante, recebendo o valor zero (0).

Foram detectadas as seguintes deficiências nos dados registrados: 1) Falta de sincronismo - os dados não foram registrados para cada cidade a partir do mesmo ano; 2) Desbalanceamento - dados coletados por diferentes períodos, resultando em muitos anos completos de coletas para alguns municípios, enquanto, para outros, não havia mais de um ano completo, e 3) Ausência sistemática de coletas de dados - para alguns anos, não havia registro nos meses finais, ou só existiam nos primeiros meses. Problemas como falta de sincronismo não puderam ser solucionados, pois se fossem aparados o início e o fim do conjunto de dados de períodos diferentes, somente dados de dez municípios seriam aproveitados. Isso corresponderia a uma redução no banco de dados, sobretudo quanto à variabilidade espacial das variáveis no Estado de São Paulo. Logo, sob o ponto de vista dos objetivos deste trabalho, optou-se por manter o maior número de diferentes pontos de registro, mesmo não havendo total sincronismo entre os pontos. Analogamente, pelo mesmo motivo, não se pôde fazer o balanceamento dos registros por município, pois a variação de diferentes períodos de registro era grande. Contudo, foi necessário extrair dados de alguns municípios em que, em determinados anos, havia uma ausência de registros consecutivos. Essa conduta foi necessária para evitar o "vício" nas estimativas. Assim, o banco de dados composto passou a ter dois municípios a menos e, por conseguinte, houve uma ligeira redução de linhas. No total, foram computados 38 municípios com dados diários, contabilizando o total de 77.687 linhas.

Análise descritiva das temperaturas

As análises inicialmente foram efetuadas com base nas variáveis contínuas, temperaturas mínimas e máximas diárias, para cada município. Assim puderam-se realizar as análises descritivas por meio das medidas de posição: Média, Mediana, Valor máximo, Valor mínimo; e de dispersão: Desvio-Padrão e Erro-Padrão.

Mineração de dados

Foram realizadas as análises sobre a variável já categorizada, denominada "Temperatura diária extrema", avaliando-se o risco de ocorrência desse fenômeno, para cada município. Por meio da variável V7 (Temperatura Diária Extrema), foi obtida uma tabela de contingência de modo a se mensurarem as probabilidades condicionais, que, por sua vez, foram estimadas pela frequência relativa. As frequências relativas, consequentemente, estimam o risco, pois esse é definido como probabilidade condicional (FLEISS, 2003). Posteriormente, foi avaliada a existência de associação entre as medidas-resumo das temperaturas (mínimas e máximas) e os valores do risco de temperaturas extremas. Para os cálculos de estatísticas descritivas e exploratórias, foi utilizado o software MINITAB 15^® (MINITAB, 2006).

RESULTADOS E DISCUSSÃO

Análise das temperaturas mínimas e máximas diárias

Os resultados das análises descritivas sobre as temperaturas (mínimas, médias e máximas) documentaram resumidamente, as características dos dados, permitindo uma descrição e comparação das temperaturas observadas nos municípios avaliados (Tabelas 1 e 2). A suposição de que as regiões de temperaturas mais elevadas tenderiam a ter maiores índices de risco, tem base em SOARES & DIAS (1986), que descrevem, em seus resultados, o efeito denominado de persistência térmica temporal, isto é, a probabilidade de ocorrência de calor num determinado dia é maior quando houve calor no dia anterior. Assim, a persistência temporal de temperaturas elevadas, em regiões mais quentes, favorece o aumento da frequência de registros de temperaturas extremas diárias, refletindo no aumento do risco.

Os municípios de Mirante do Paranapanema e São José do Rio Preto foram aqueles com temperaturas mínimas em patamares mais elevados, visto que as médias das temperaturas mínimas já alcançaram valores próximos de uma temperatura mínima de alerta (MACHADO, 1950). Por outro lado, o município de Piedade foi aquele que apresentou a menor média das temperaturas mínimas entre todas as cidades estudadas, sendo menor a possibilidade de ocorrência de um dia de temperatura extrema.

Considerada a referência de temperatura máxima de 32 °C (MACHADO, 1950; SARTORI, 2000 e ROSSATO et al., 2003), foi possível verificar que os municípios de Barretos e Tupi Paulista são passíveis de alcançar temperaturas máximas diárias críticas, que afetam negativamente o desempenho produtivo das aves (DEATON et al., 1997; ARADAS et al., 2005; MIRAGLIOTTA et al., 2006 e VALE et al., 2008). Contudo, na região onde se localiza o município de Piedade, houve menor valor médio de temperatura máxima, sugerindo que, nessa região, não haverá risco de ocorrência de temperatura extrema.

Análise do risco de temperatura extrema diária

Os riscos foram obtidos pela estimativa das probabilidades condicionais e foram calculados a partir de uma tabela de contingência, que foi configurada especialmente para a avaliação e comparação desses valores. A obtenção da tabela de contingência foi proveniente de um cruzamento entre duas variáveis categóricas que possuíam diferentes classes (FLEISS, 2003). A Tabela 3 apresenta os riscos de ocorrência de temperaturas extremas, dado um determinado município. Verifica-se que Andradina possui um risco de 11,22%, ou seja, valor que corresponde à probabilidade de haver, em qualquer dia, uma temperatura extrema nesse município. Foram encontrados os seguintes municípios com maior risco de incidência de temperatura extrema diária, que provocaria problemas de estresse térmico na produção de frangos de corte, conforme indicam YAHAV et al. (1995) e COOPER & WASHBURN (1998): Tupã (29,04%), Araçatuba (23,06%), São José do Rio Preto (21,90%), Osvaldo Cruz (19,84%) e Penápolis (19,28%). Também foram apontados os municípios de risco zero (0,00%): Atibaia, Monte Alegre do Sul e Piedade.

Comparando-se as Tabelas 1; 2 e 3, constatou-se que a maioria dos municípios que possuem temperaturas mais elevadas, pertence à região oeste do Estado de São Paulo. Verificou-se também que a classificação obtida pelos riscos de temperaturas extremas está mais vinculada com o perfil da temperatura mínima, do que o da temperatura máxima.

Os municípios de Araçatuba e de São José do Rio Preto foram os que apresentaram temperaturas ambientais extremas mais altas e com elevados riscos de temperaturas extremas diárias, o que os torna vulneráveis ao impacto negativo do efeito térmico sobre a criação de aves, de acordo com o preconizado por DEATON et al.(1997), MIRAGLIOTTA et al. (2006) e VALE et al. (2008). Nesse contexto, é dada maior atenção à região de São José do Rio Preto, pois é um dos municípios de intensa produção avícola de corte (APA, 2008). Assim, os produtores devem ter mais investimentos e alternativas para mitigar as elevadas temperaturas ambientais, às quais seus alojamentos estão submetidos.

As evidências da relação entre determinadas características das temperaturas e os riscos de temperaturas extremas foram encontradas por meio da análise de Componentes Principais (Figura 1). Essa análise permitiu a visualização conjunta da relação entre as medidas estatísticas das temperaturas e o risco de temperaturas extremas. As medidas referentes às temperaturas mínimas estão acompanhadas de "(T-mín)", enquanto aquelas referentes às temperaturas máximas estão acompanhadas de "(T-máx)". Associações dessa natureza foram encontradas por PEREIRA (2005) e SALGADO (2006), quando utilizaram metodologia semelhante.

CONCLUSÕES

A partir de dados meteorológicos históricos, foi possível apontar as regiões onde se localizam os municípios do Estado de São Paulo com maiores temperaturas (mínimas e máximas), bem como calcular o risco de ocorrência de temperaturas extremas. Fazendo-se um comparativo dos resultados, foram destacados os municípios de Araçatuba, Osvaldo Cruz, Mirante do Paranapanema, São José do Rio Preto e Tupã. Dentre as regiões de maior risco à ocorrência de temperaturas diárias extremas, destaca-se a de São José do Rio Preto (aproximadamente 22% de risco), que é de intensa produção avícola. Verificou-se que a média e, principalmente, as medianas das temperaturas mínimas têm forte relação com os valores encontrados de risco de temperaturas extremas, isso pode resultar em diminuição em custos operacionais e computacionais quando se deseja estimar esse risco.

Recebido pelo Conselho Editorial em: 20-11-2008

Aprovado pelo Conselho Editorial em: 4-11-2009

• APA. ASSOCIAÇÃO PAULISTA DOS AVICULTORES. Municípios de intensa produção avícola no Estado de São Paulo Disponível em: <http://www.apa.com.br>. Acesso em: 15 fev. 2008.
• ARADAS, M.E.C.; NÄÄS, I.A.; SALGADO, D.D. Comparing the thermal environment in broiler houses using two bird densities under tropical conditions. CIGR E-Journal, College Station, v.7, p.1-12, 2005.
• BRUNOZILLI, J.; SOUZA, D.P. Produtividade regional. Avicultura Industrial, Itu, v.1119, n.2, p.22-24, 2004.
• CEPAGRI. CENTRO DE PESQUISAS METEOROLÓGICAS E CLIMÁTICAS APLICADAS À AGRICULTURA. Clima das cidades do Estado de São Paulo Disponível em: <http://www.cpa. unicamp.br/produção/trabalhos-em-anais-de-eventos/2003 /tema1_2.pdf/download.pdf>. Acesso em: 24 set. 2008.
• COOPER, M.A.; WASHBURN, K.W. The relationships of body temperature to weight gain, feed consumption, and feed utilization in broilers under heat stress. Poultry Science, Stanford, v.77, n.2, p.237-242, 1998.
• DEATON, J.W.; REECE, F.N.; MCNAUGHTON, N. The effect of temperature during the growing period on broiler performance. Poultry Science, Ithaca, v.57, n.1, p.1.070-1.074, 1997.
• FLEISS, J.L.; LEVIN, B.; PAIK, M.C. Statistical methods for rates and proportions New York: John Wiley & Sons, 2003. 328 p. (Series in Probability and Statistics)
• HAHN, G.L. Compensatory performance in livestock: influences on environmental criteria. In: INTERNATIONAL LIVESTOCK ENVIRONMENT SYMPOSIUM, 2., 1982. Proceedings.. St. Joseph: ASABE, 1982. p.285-294.
• IBGE. INSTITUTO BRASILEIRO DE GEOGRAFIA E ESTATÍSTICA. Produção agropecuária 2006 Disponível em: <http://www.ibge.gov.br/home/presidencia/noticias/21052004 biomashtml. shtm>. Acesso em: 10 jan. 2008.
• JONES, R. Managing climate risks. In: OECD WORKSHOP ON THE BENEFITS OF CLIMATE CHANGE: IMPROVING INFORMATION FOR POLICY MAKERS, 1., 2003, Paris. Proceedings.. Paris: OECD, 2003. p.37.
• MACARI, M.; GONZALES, E. Fisiopatogenia da síndrome da morte súbita em frangos de corte. In: CONFERÊNCIA APINCO DE CIÊNCIA E TECNOLOGIA AVÍCOLAS, 1990, Campinas. Anais.. Campinas: APINCO, 1990. p.65-73.
• MACHADO, F.P. Contribuição ao estudo do clima do Rio Grande do Sul Rio de Janeiro: IBGE, 1950. p.91.
• MINITAB INC. Minitab^® 15.1.0.0 State College, 2006.
• MIRAGLIOTTA, M.Y.; NÄÄS, I. de A.; MANZIONE, R.L.; NASCIMENTO, F.F. Spatial analysis of stress conditions inside broiler house under tunnel ventilation. Scientia Agricola, Piracicaba, v.63, n.5, p.426-432, 2006.
• MOURA, D.J. Ambiência na produção de aves de corte In: SILVA, I.J.O. Ambiência na produção de aves em clima tropical. Jaboticabal: Associação Brasileira de Engenharia Agrícola, 2001. v.2, p.185-198. (Série Engenharia Agrícola)
• NÄÄS, I.A.; CURTO, F.P. Avicultura de precisão. In: DA SILVA, I. J. O. Ambiência na produção de aves em clima tropical. Jaboticabal: Associação Brasileira de Engenharia Agrícola, 2001. v.2, p.1-30. (Série Engenharia Agrícola)
• NIENABER, J.A.; HAHN, G.L. Engineering and management practices to ameliorate livestock heat stress. In: INTERNATIONAL SYMPOSIUM OF THE CIGR. NEW TRENDS IN FARM BUILDINGS, 1., 2004, Évora. Proceedings.. Évora: CIGR, 2004. 1 CD-ROM.
• PEREIRA, D.F. Metodologia para estimativa de bem-estar de matrizes de frango de corte, utilizando monitoramento digital e construção de modelos de simulação 2005. 155 f. Tese (Doutorado em Construções Rurais e Ambiência) - Faculdade de Engenharia Agrícola, Universidade Estadual de Campinas, Campinas, 2005.
• ROPPA, L. Evolução do plantel total de suínos Disponível em: <http://www. porkworld.com.br> Acesso em: 30 mar. 2004.
• ROSSATO, P.S.; SARTORI, M.G.B.; MISSIO, L.R. As ondas de calor na região central do RS entre os meses de maio a outubro. In: SIMPÓSIO BRASILEIRO DE GEOGRAFIA FÍSICA APLICADA, 10., 2003, Rio de Janeiro. Anais... Rio de Janeiro: DGFA, 2003. 1 CD-ROM.
• SALGADO, D.D. Modelo estatístico para predição de bem-estar de reprodutoras de frango de corte baseado em dados de ambiente e análise do comportamento 2006. 126 f. Dissertação (Mestrado em Construções Rurais e Ambiência) - Faculdade de Engenharia Agrícola, Universidade Estadual de Campinas, Campinas, 2006.
• SARTORI, M.G.B. Clima e percepção 2000. 138 f. Tese (Doutorado) - Faculdade de Filosofia, Letras e Ciências Humanas, Universidade de São Paulo, São Paulo, 2000.
• SCHNEIDER, S.H. Abrupt non-linear climate change, irreversibility and surprise. Global Environmental Change, East Anglia, v.14, n.1, p.245-258, 2004.
• SILVA, R.G. Introdução à bioclimatologia animal São Paulo: Nobel, 2000. 285 p.
• SOARES, J.R.; DIAS, M.A.F.S. Probabilidade de ocorrência de alguns eventos meteorológicos extremos na cidade de São Paulo. Revista Brasileira de Meteorologia, Rio de Janeiro, v.1, n.2, p.67-75, 1986.
• TANKSON, J.D.; VIZZIER-THAXTON, Y.; THAXTON, J.P.; MAY, J.D.; CAMERON, J.A. Stress and nutritional quality of broilers. Poultry Science, Stanford, v.80, n.9, p.1.384-1.389, 2001.
• UBA. UNIÃO BRASILEIRA DE AVICULTURA. Relatório 2007-2008 Disponível em: <http://www.uba.org.br/index.htm>. Acesso em: 20 set. 2008.
• VALE, M.M.; MOURA, D.J.; NÄÄS, I.A.; OLIVEIRA, S.R.M.; RODRIGUES, L.H.A. Data mining to estimate broiler mortality when exposed to heat wave. Scientia Agricola, Piracicaba, v.65, n.3, p.223-229, 2008.
• YAHAV, S.; GOLDFELD, S.; PLAVNIK, I.; HURWITZ, S. Physiological response of chickens and turkeys to relative humidity during exposure to high ambient temperature. Journal of Thermal Biology, Durham, v.20, n.3, p.245-253, 1995.

Avaliação de risco à produção de frango de corte do Estado de São Paulo em função da temperatura ambiente

Risk evaluation of broiler production in the state of São Paulo as function of ambient temperature

Datas de Publicação

• Publicação nesta coleção
  20 Jul 2010
• Data do Fascículo
  Jun 2010

Histórico

• Recebido
  20 Nov 2008
• Aceito
  04 Nov 2009