Microsoft Word - 42-Agra_18096 390 Original Article Biosci. J., Uberlandia, v. 30, supplement 1, p. 390-400, June/14 PRODUTIVIDADE DA CANA-DE-AÇÚCAR: VARIABILIDADE LINEAR E ESPACIAL ENTRE COMPONENTES TECNOLÓGICOS E DA PRODUÇÃO PRODUCTIVITY OF CANE SUGAR: LINEAR AND SPATIAL VARIABILITY BETWEEN COMPONENTS AND TECHNOLOGICAL PRODUCTION Flávio Carlos DALCHIAVON1; Morel de Passos e CARVALHO2; Rafael MONTANARI2; Marcelo ANDREOTTI2; Alan Rodrigo PANOSSO2 1. Professor Efetivo do Instituto Federal de Mato Grosso, Campus Campo Novo do Parecis, Departamento de Agronomia, Campo Novo do Parecis, MT, Brasil. flavio.dalchiavon@cnp.ifmt.edu.br; 2. Professore, Doutor, Faculdade de Engenharia de Ilha Solteira - UNESP, Ilha Solteira, SP. RESUMO: Na atualidade, a cultura da cana-de-açúcar desempenha frente à realidade brasileira, sobretudo no aspecto relacionado às fontes energéticas alternativas, importantíssimo papel. No ano de 2009, no município de Suzanápolis (SP), no Cerrado Brasileiro, foi instalado um experimento com a cultura da cana-de-açúcar em um Argissolo Vermelho eutrófico, com o objetivo de selecionar, por meio de coeficientes de correlação de Pearson, da modelagem de regressões lineares simples e de múltiplas e da correlação espacial, os melhores componentes tecnológicos e os de produção, para explicar a variabilidade da produtividade da cana-de-açúcar. Foi instalada a malha geoestatística, para a coleta de dados, com 120 pontos amostrais, numa área de 14,53 ha. Para as regressões lineares simples, a população de plantas é o componente da produção que apresenta a melhor correlação quadrática com a produtividade da cana-de-açúcar, dada por: PRO = -0,553**xPOP2+16,14*xPOP-15,77. Entretanto, para as regressões lineares múltiplas, a equação PRO = - 21,11+4,92xPOP**+0,76xPUR** é a que melhor se apresenta para estimar a referida produtividade. Espacialmente, a melhor correlação com a produtividade da cana-de-açúcar é determinada, também, pelo componente da produção da população de plantas. PALAVRAS-CHAVE: Energia renovável. Geoestatistica. Sacchharum spp. Regressões. INTRODUÇÃO A importância socioeconômica da cultura da cana-de-açúcar (Sacchharum spp.) para o Brasil é inquestionável, sendo a principal matéria-prima para a produção do etanol utilizado pelos veículos automotores, além da sua destinação para a produção do açúcar, um dos principais produtos de exportação nacional. De acordo com estimativa da Conab (2012), na atual safra, o setor sucroalcooleiro nacional processará uma produção de 602,179 milhões de toneladas de colmos em 8,567 milhões de hectares, tendo o estado de São Paulo como o maior produtor, participando com 53,6% desta (323,131 milhões de toneladas), numa área de 4,426 milhões de hectares, e produtividade média em torno 73,0 t ha-1. O conhecimento da produtividade da cana- de-açúcar pode ser obtidopelos seus componentes tecnológicos e da produção, utilizando técnicas estatísticas que permitam estimá-la. Dewey e Lu (1959) foram os pesquisadores pioneiros a utilizar esta metodologia em plantas de trigo, destacando que os métodos de correlação linear e regressão múltipla foram os primeiros métodos a serem utilizados nas análises dos componentes da produção. Gioia (2011) e Lima (2012), mediante observações em seus trabalhos, afirmaram quea utilização dos componentes tecnológicos (açúcares totais recuperáveis, percentual de sólidos solúveis totais, percentual de sacarose, pureza e percentual de fibra) e da produção (volume de colmos por hectare e população de plantas por metro quadrado), enquanto variáveis independentes, tanto linear quanto espacialmente, para estimar a produtividade da cana-de-açúcar (variável dependente), são ferramentas extremamente importantes, das quais os agricultores dispõem, uma vez que alterações nesses componentes são responsáveis diretos pelo ajuste da produtividade. Entretanto, sabe-se que existe variabilidade em tais componentes entre indivíduos de uma população em função de vários fatores, que vão desde a implantação da cultura (influenciando o número de plantas por área), até a disponibilidade de assimilados (modificada pelo arranjo de plantas). Assim, para avaliar o comportamento espacial desses componentes e da produtividade, pode-se utilizar de técnicas geoestatísticas, que, por meio de semivariogramas, permitem definir o alcance e o grau da dependência espacial, fornecendo parâmetros para a elaboração de mapas de variabilidade pela krigagem e/ou cokrigagem (MONTANARI et al., 2010; ZANÃO JÚNIOR et al., 2010). O semivariograma é um gráfico que caracteriza a estrutura da variância do atributo estudado como variável regionalizada em função da Received: Accepted: 391 Produtividade da cana-de-açúcar… DALCHIAVON, F. C. et al Biosci. J., Uberlandia, v. 30, supplement 1, p. 390-400, June/14 distância no terreno, tendo o alcance da dependência espacial como principal parâmetro fornecido pela geoestatística, representando a distância dentro da qual os valores de um determinado atributo são iguais entre si (DALCHIAVON et al., 2011a). Por outro lado, a cokrigagem pode ser mais precisa do que a krigagem de um atributo simples quando o semivariograma cruzado apresentar dependência entre dois atributos quaisquer (MONTANARI et al., 2012). Além do que, exclusivamente em relação à variável primária, que normalmente é sub- amostrada, essa técnica possibilita efetuar a estimativa dos seus valores não amostrados. Neste contexto, é reconhecido o uso e contribuição destas técnicas na identificação e mapeamento da variabilidade intrínseca aos componentes de produção dos cultivos agrícolas. Segundo Pires (2002), estas técnicas podem permitir avaliar a variabilidade espacial para características como produtividade, bem como a ocorrência de padrões (manchas) de áreas de lavoura associadas, por exemplo, com os componentes da planta. Este tipo de estudo permite verificar, mesmo havendo diferenças entre plantas em determinados componentes relacionados à produtividade, se há possibilidade de se estabelecer áreas com comportamento similar, ou seja, se existe dependência espacial de outros componentes, que não a produtividade de colmos, comumente avaliada em nível de lavoura. Contudo, a variabilidade linear e espacial dos componentes, tecnológicos e da produção, e sua correlação com a produtividade da cana-de-açúcar têm sido pouco pesquisadas. Diante do exposto, o presente trabalho teve por objetivo, através de coeficientes de correlação de Pearson, da modelagem de regressões lineares simples e múltiplas e da correlação espacial, selecionar os melhores componentes, tecnológicos e da produção, para explicar a variabilidade da produtividade da cana-de-açúcar. MATERIAL E MÉTODOS O experimento foi desenvolvido no ano de 2009, na Fazenda Caiçara, pertencente à Usina Vale do Paraná S/A Álcool e Açúcar, no município de Suzanápolis (SP) no noroeste do estado de São Paulo, cujas coordenadas são 20°28'10'' S e 50°49'20'' W.O solo no qual a malha experimental foi instalada foi classificado como Argissolo Vermelho eutrófico típico textura arenosa/média. Na Figura 1 estão apresentadas as precipitações e temperaturas médias durante o período experimental, cujos valores foram: precipitação acumulada de 2881 mm, temperatura máxima, média e mínima de 31,4; 25,1 e 19,9ºC, respectivamente. Antes da instalação do experimento, foram coletadas 20 amostras simples de solo, nas profundidades de 0-0,20 m e 0,20-0,40 m, das quais obteve uma amostra composta representativa para cada profundidade, para determinação da sua fertilidade, cujos resultados estão contidos na Tabela 1. Figura 1. Precipitação pluvial e médias térmicas obtidas durante o período experimental (set/2007 a mai/2009). 392 Produtividade da cana-de-açúcar… DALCHIAVON, F. C. et al Biosci. J., Uberlandia, v. 30, supplement 1, p. 390-400, June/14 Tabela 1. Caracterização inicial da fertilidade do solo nas profundidades de 0 -0,20 m e 0, 20 - 0,40 m. Atributo químico (a) pH Complexo sortivo Prof. MO CaCl2 P K + Ca2+ Mg2+ H++Al3+ Al3+ SB CTC V% m% m g dm-3 mg dm-3 ----------------------------------mmolc dm -3-------------------- ------------- -------%----- --- 0-0,20 16 5,1 4,3 1,4 13,9 7,3 16,6 0,7 22,6 39,2 57,6 3,0 0,20-0,40 14 5,2 4,1 0,9 12,5 7,6 16,4 0,7 21,0 37,4 56,1 3,2 (a) MO = matéria orgânica, V% = índice de saturação por bases, m% = índice de saturação por alumínio. Para a implantação da cana-de-açúcar (variedade SP 79 1011), em área sob pastagem há 10 anos, foi dessecada a braquiária com Roundup WG (glifosato), utilizando 2,5 kg ha-1 p.c.. No preparo do solo foi realizada uma gradagem pesada e duas intermediárias. Aplicou-se a lanço 2,0 t ha-1 de calcário, incorporando-o a 0,35 m de profundidade com um arado de aivecas. Para finalizá-lo, utilizou-se uma grade niveladora. O espaçamento de plantio (20/jun/2006) foi de 1,50 m entre sulcos, utilizando-se 500 kg ha-1 da fórmula 06-30-24 (NPK). Após o primeiro corte da cultura (set/2007 - colheita mecanizada) foi aplicada 1 t ha-1 de gesso agrícola e 500 kg ha-1 da fórmula 18-00-27 na adubação da soqueira. Na coleta de dados, a cultura encontrava-se no segundo corte e bisada (21 meses em campo), assim realizou-se a colheita manual após queima da cana-de-açúcar (mai/2009). Na alocação da malha experimental foram definidas as direções x e y do sistema de coordenadas cartesianas, sendo constituída de 120 pontos amostrais, distribuídos numa área de 14,53 ha (418,46 m x 349,00 m). O espaçamento utilizado entre os pontos na grande malha foi de 42 m x 43 m, enquanto que na malha de refinamento, construída com a finalidade de detectar alcances da dependência espacial para espaçamentos menores do que 40 m, composta por 3 pontos adjacentes ao ponto central, no espaçamento de 9 m (Figura 2). Figura 2. Croqui da malha experimental de amostragem. Os atributos pesquisados foram: a) produtividade de colmos por hectare (PRO), b) volume de colmos por hectare (VOL), c) população de plantas por metro quadrado (POP), d) açúcares totais recuperáveis (ATR), e) percentual de sólidos solúveis totais (BRI), f) percentual de sacarose (POL), g) pureza (PUR) e h) percentual de fibra (FIB). A PRO (t ha-1) foi obtida pela pesagem dos colmos contidos numa área útil de 9 m2, sendo duas linhas de 3 m (3m x 3m) por ponto amostral. Para o VOL (m3 ha-1) foram coletados cinco colmos, sendo medidos o comprimento e os diâmetros (da base, do 393 Produtividade da cana-de-açúcar… DALCHIAVON, F. C. et al Biosci. J., Uberlandia, v. 30, supplement 1, p. 390-400, June/14 meio e do ápice), obtendo o diâmetro médio utilizado no cálculo final do VOL, a partir da POP (pl. m-2), obtida por contagem, do respectivo ponto amostral, com área útil de 9 m2. Os atributos tecnológicos ATR (kg t-1), BRI, POL, PUR e FIB (%), foram determinados no laboratório da Usina Vale do Paraná S/A Álcool e Açúcar, conforme metodologia do Consecana (2006). Para cada atributo foi efetuada a análise descritiva com o auxilio da estatística clássica, utilizando o Software SAS. Seguidamente, montou- se a matriz de correlação, objetivando efetuar as correlações lineares simples para as combinações entre os atributos estudados, e modelaram-se as regressões lineares simples com a planilha de cálculos Excel, conforme Dalchiavon e Carvalho (2012). Por outro lado, utilizando-se o pacote computacional SAS, foram efetuadas as regressões lineares múltiplas da PRO (variável dependente) em função dos demais atributos pesquisados (variáveis independentes), no intuito de selecionar aqueles que, nos devidos casos, proporcionariam as melhores equações para estimar a PRO, por intermédio do stepwise, a 10% de probabilidade. O desempenho das equações foi analisado pelo coeficiente de determinação (R2) e pela contribuição de cada componente no modelo, estimada pelo aumento observado na soma de quadrados devido à regressão, quando da adição de cada variável no modelo, conforme Maia et al. (2001). Para cada atributo foi analisada a dependência espacial pelo cálculo do semivariograma simples, com base nos pressupostos de estacionaridade da hipótese intrínseca pelo uso do pacote Gamma Design Software 7.0 (Gs+, 2004). Os ajustes dos semivariogramas simples, em função de seus modelos, foram efetuados conforme Dalchiavon et al. (2012), retirando-se a tendência dos dados, por meio da técnica da regressão múltipla polinomial, para os atributos que não apresentaram dependência espacial. A decisão final do modelo que representou o ajuste foi atestada pela validação cruzada, assim como para a definição do tamanho da vizinhança que proporcionou a melhor malha de krigagem e/ou cokrigagem. Para cada atributo foram relacionados o efeito pepita (Co), o patamar (Co+C) e o alcance (Ao). A análise do ADE foi efetuada conforme a seguinte expressão (Gs+, 2004): ADE = [C/(C+Co)].100......................................(1) onde: ADE é o avaliador da dependência espacial; C, a variância estrutural; e C+Co, o patamar, adotando-se a interpretação proposta por Dalchiavon et al. (2012). RESULTADOS E DISCUSSÃO Os atributos pesquisados que apresentaram distribuição de frequência do tipo normal foram a PRO, VOL, POP, BRI e FIB (Tabela 2), sendo significativos pelo teste de normalidade de Shapiro e Wilk a 5%, uma vez que suas probabilidades oscilaram entre 0,066 (FIB) e 0,403 (VOL), indicando que as medidas de tendência central não refletem valores atípicos na distribuição. A PRO obtida (89,2 t ha-1) está de acordo com a média brasileira (Unica, 2009), entretanto, esperava-se maior PRO, haja vista que se trata de cana com 21 meses de ciclo. Tabela 2. Análise descritiva da produtividade e de componentes tecnológicos da cana-de-açúcar em um Argissolo Vermelho eutrófico da Fazenda Caiçara (Suzanápolis-SP). Atributo (a) Medidas estatísticas descritivas Média Mediana Valor Desvio Padrão Coeficiente Probabilidade do teste (b) Mínimo Máximo Variação (%) Curtose Assimetria Pr110,4 t ha-1 (PRO), 10,0 a > 112,0 pl. m-2 (£POP) e 14,4 a >16,1% (BRI). Por outro lado, em pequenas áreas dos mapas foram observados os menores valores, que ficaram entre 55,8-83,1 t ha-1 (PRO), 8,0-10,0 pl. m-2 (£POP) e 12,7-14,4% (BRI). O mapa da PRO apresentou manchas com produtividades extremamente variadas, fato que também fora notado por Molin et al. (2007), Dalchiavon et al. (2011b) e Dalchiavon e Carvalho (2012), ao pesquisarem a produtividadedo milho safrinha, do feijoeiro invernal e da soja, respectivamente. Na Tabela 6 podem ser averiguados ainda os parâmetros dos semivariogramas cruzados, bem como na Figura 4 seus respectivos semivariogramas cruzados e mapas de cokrigagem. Portanto, a melhor cokrigagem foi da PRO=f(#POP), pois explicou 78,1% da variabilidade da PRO. Assim, do ponto de vista espacial da área pesquisada, nos sítios em que a variedade de cana-de-açúcar SP 79 1011 apresentou £POP variando entre 10,0 e >12,0 pl. m-2 (Figura 3d), a PRO esteve compreendida entre 77,0 e >120,0 t ha-1 (Figura 4b). Em contrapartida, naqueles onde a £POP esteve entre 8,0 e 10,0 pl. m- 2, a PRO ficou compreendida entre 33,0 e 77,0 t ha-1. Para a presente cokrigagem, o modelo ajustado foi o exponencial direto, com alcance de 357,3 m e ADE médio (50,0%). Da mesma forma, as Figuras 4c,d evidenciam a cokrigagem da PRO=f(BRI), cujo modelo semivariográfico foi o gaussiano direto, com 73,0 m de alcance da dependência espacial e ADE alta (63,6%), explicando 23,4% da variabilidade espacial da PRO (Tabela 6), ficando plenamente de acordo com a cokrigagem 397 Produtividade da cana-de-açúcar… DALCHIAVON, F. C. et al Biosci. J., Uberlandia, v. 30, supplement 1, p. 390-400, June/14 apresentada por Lima (2012). Portanto, nos sítios em que o BRI ficou entre 14,4 e > 16,1% (Figura 3f),a PRO esteve compreendida entre 78,0 e > 120,0 t ha-1 (Figura 4d), ao passo que naqueles onde o BRI esteve entre 12,7 e 14,4%, a PRO ficou compreendida entre 35,0 e 78,0 t ha-1. Figura 3. Semivariogramas simples e mapas de krigagem de atributos da cana-de-açúcar em um Argissolo Vermelho eutrófico da Fazenda Caiçara (Suzanápolis-SP) (B) (A) (C) (D) (E) (F) 398 Produtividade da cana-de-açúcar… DALCHIAVON, F. C. et al Biosci. J., Uberlandia, v. 30, supplement 1, p. 390-400, June/14 Figura 4. Semivariogramas cruzados e mapas de cokrigagem de atributos da cana-de-açúcar em um Argissolo Vermelho eutrófico da Fazenda Caiçara (Suzanápolis-SP) CONCLUSÕES Para as regressões simples, a população de plantas é o componente da produção que apresenta a melhor correlação quadrática com a produtividade da cana-de-açúcar, dada por: PRO = - 0,553**.POP2+16,14*.POP-15,77. Entretanto, para as regressões lineares múltiplas, a equação PRO = - 21,11+4,92.POP**+0,76.PUR** é a que melhor se apresenta para estimar a referida produtividade. Do ponto de vista espacial, a melhor correlação com a produtividade da cana-de-açúcar é determinada, também, pelo componente da produção população de plantas. ABSTRACT: Nowadays, the culture of the sugarcane plays an important role regarding the Brazilian reality, especially in the aspect related to the alternative energy sources. In 2009, the municipality of Suzanápolis (SP), in the Brazilian Cerrado, an experiment was conducted with the culture of the sugarcane in a Red eutrophic, with the aim of selecting, using Pearson correlation coefficients, modeling, simple, linear and multiple regressions and spatial correlation, and also the best technological and productive components, to explain the variability of the productivity of the sugarcane. The geostatistical grid was installed in order to collect the data, with 120 sampling points, in an area of 14.53 ha. For the simple linear regressions, the plants population is the component of production that presents the best quadratic correlation with the productivity of the sugarcane, given by: PRO = -0.553**xPOP2+16.14*xPOP-15.77. However, for multiple linear regressions, the equation PRO = -21.11+4.92xPOP**+0.76xPUR** is the one that best presents in order to estimate that productivity. Spatially, the best correlation with yield of the sugarcane is also determined by the component of the production population of plants. KEYWORDS: Renewable energy. Geostatistical. Sacchharum spp. Regressions. (B) (A) (C) (D) 399 Produtividade da cana-de-açúcar… DALCHIAVON, F. C. et al Biosci. J., Uberlandia, v. 30, supplement 1, p. 390-400, June/14 REFERÊNCIAS CHAVES, L. H. G.; FARIAS, C. H. A. Variabilidade espacial do estoque de carbono nos Tabuleiros Costeiros da Paraíba: Solo cultivado com cana-de-açúcar. Revista Brasileira de Ciências Agrárias, Recife, v. 3, n. 1, p.20-25, 2008. COMPANHIA NACIONAL DE ABASTECIMENTO - CONAB. Acompanhamento da safra brasileira: Cana-de-açúcar. Primeiro Levantamento, abril/2012, Brasília: CONAB, 2012. 19p. CONSELHO DOS PRODUTORES DE CANA DE AÇÚCAR, AÇÚCAR E ÁLCOOL DO ESTADO DE SÃO PAULO - CONSECANA. Manual de instruções. 5.ed. Piracicaba, 2006. 112p. CORÁ, J. E.; ARAUJO, A. V.; PEREIRA, G. T.; BERALDO, J. M. G. Variabilidade espacial de atributos do solo para adoção do sistema de agricultura de precisão na cultura de cana-de-açúcar. Revista Brasileira de Ciência do Solo, Viçosa, v. 28, n. 6, p.1013-1021, 2004. DALCHIAVON, F. C.; CARVALHO, M. P.: FREDDI, O. S.; ANDREOTTI, M.; MONTANARI, R. Variabilidade espacial da produtividade do feijoeiro correlacionada com atributos químicos de um Latossolo Vermelho Distroférrico sob sistema de semeadura direta. Bragantia, Campinas, v. 70, n. 4, p. 908-916, 2011b. DALCHIAVON, F. C.; CARVALHO, M. P.; ANDREOTTI, M.; MONTANARI, R. Variabilidade espacial de atributos da fertilidade de um Latossolo Vermelho Distroférrico sob Sistema Plantio Direto. Revista Ciência Agronômica, Fortaleza, v. 43, n. 3, p. 453-461, 2012. DALCHIAVON, F. C.; CARVALHO, M. P. Correlação linear e espacial dos componentes de produção e produtividade da soja. Semina: Ciências Agrárias, Londrina, v. 33, n. 2, p. 541-552, 2012. DALCHIAVON, F. C.; CARVALHO, M. P.; NOGUEIRA, D.C.; ROMANO, D.; ABRANTES, F. L.; ASSIS, J. T.; OLIVEIRA, M. S. Produtividade da soja e resistência mecânica à penetração do solo sob sistema plantio direto no cerrado brasileiro. PesquisaAgropecuária Tropical, Goiania, v. 41, n. 1, p. 8-19, 2011a. DEWEY, D. R.; LU, K. H. A correlation path coefficient analysis of components of crested wheatgrass seed production. AgronomyJournal, Madison, v. 51, n. 2, p. 515-518, 1959. GIOIA, M. T. Produtividade e qualidade da cana soca correlacionadas com atributos físico-químicos de um Argissolo Vermelho distrófico do noroeste Paulista. 2011. 127f. Dissertação (Mestrado em Agronomia) - Universidade Estadual Paulista, Ilha Solteira, 2011. GS+: Geostatistics for environmental sciences. 7. ed. Michigan, Plainwell: Gamma Desing Software, 2004. 159p. LANDELL, M. G. A.; SILVA, M. A. Manual do experimentador: melhoramento da cana-de-açúcar. In: Metodologia de Experimentação: ensaios de competição em cana-de-açúcar. Pindorama: Instituto Agronômico, 2004. Não paginado. (Apostila de Treinamento Interno) LAVANHOLI, M. G. D. P. Qualidade da cana-de-açúcar como matéria prima para produção de açúcar e álcool. In: DINARDO-MIRANDA, L. L.; VASCONCELOS, A. C. M.; LANDELL, M. G. A. (Ed.). Cana-de-açúcar. Campinas: Instituto Agronômico e Fundação: IAC, 2008. Cap.32, p. 697-722. LIMA, C. G. R. Análise linear e espacial entre alguns atributos produtivos e tecnológicos da cana-de- açúcar com o pH de um Argissolo Vermelho de Suzanápolis (SP). 2012. 101f. Tese (Doutorado em Agronomia) - Universidade Estadual Paulista, Ilha Solteira, 2012. MAIA, C. E.; MORAIS, E. R. C.; OLIVEIRA, M. Classificação da composição iônica da água de irrigação usando regressão linear múltipla. Revista Brasileira de Engenharia Agrícola e Ambiental, Campina Grande, v. 5, n. 1, p. 55-59, 2001. 400 Produtividade da cana-de-açúcar… DALCHIAVON, F. C. et al Biosci. J., Uberlandia, v. 30, supplement 1, p. 390-400, June/14 MOLIN, J. P. VIEIRA JUNIOR, P. A.; DOURADO NETO, D.; FAULIN, G. C.; MASCARIN, L. Variação espacial na produtividade de milho safrinha devido aos macronutrientes e à população de plantas. Revista Brasileira de Milho e Sorgo, Sete Lagoas, v. 6, n. 3, p. 309-324, 2007. MONTANARI, R.; CARVALHO, M. P.; ANDREOTTI, M.; DALCHIAVON, F. C.; LOVERA, L. H.; HONORATO, M. A. O. Aspectos da produtividade do feijão correlacionados com atributos físicos do solo sob elevado nível tecnológico de manejo. Revista Brasileira de Ciência do Solo, Viçosa, v. 34, n. 6, p. 1811-1822, 2010. MONTANARI, R.; ZAMBIANCO, E. C.; CORRÊA, A. R.; PELLIN, D. M. P.; CARVALHO, M. P.; DALCHIAVON, F. C. Atributos físicos de um Latossolo Vermelho correlacionados linear e espacialmente com a consorciação de guandu com milheto. Revista Ceres, Viçosa, v. 59, n. 1, p. 407-410, 2012. PIRES, J. L. F. Estimativa do potencial produtivo da soja e variabilidade espacial de área de produção. 2002. Tese (Doutorado em Fitotecnia) - Faculdade de Agronomia. Universidade Federal do Rio Grande do Sul, Porto Alegre, 2002. RAPASSI, R. M. A.; TARSITANO, M. A. A.; BOLONHEZI, A. C. Avaliação técnica e econômica de sistemas de produção da cana-de-açúcar (Saccharum spp) na região oeste do estado de São Paulo. Informações Econômicas, São Paulo, v. 39, n. 1, p. 11-21, 2009. RIPOLI, T. C. C.; RIPOLI, M. L. C. Biomassa de cana-de-açúcar: colheita, energia e ambiente. Piracicaba: Barros e Marques Ed. Eletrônica, 2004. 302p. SCAPARI, M. S.; BEUCLAIR, E. G. F. Anatomia e botânica. In: DINARDO-MIRANDA, L. L.; VASCONCELOS, A. C. M.; LANDELL, M. G. A. (Ed.). Cana-de-açúcar. Campinas: Instituto Agronômico e Fundação- IAC, 2008. Cap. 2, p. 47 - 56. SILVA, M. A. Interação genótipo x ambiente e estabilidade fenotípica de cana-de-açúcar em ciclo de cana de ano. Bragantia, Campinas, v. 67, n. 1, p. 109-117, 2008. SOUZA, Z. M.; CERRI, D. G. P.; COLET, M. J.; RODRIGUES, L. H. A.; MAGALHÃES, P. S. G; MANDONI, R. J. A. Análise dos atributos do solo e da produtividade da cultura de cana-de-açúcar com o uso da geoestatística e árvore de decisão. Ciência Rural, Santa Maria, v. 40, n. 4, p. 840-847, 2010. SOUZA, Z. M.; CERRI, P. D. G.; MAGALHÃES, P. G.; CAMPOS, M. C. C. Correlação dos atributos físicos e químicos do solo com a produtividade de cana-de-açúcar. Revista de Biologia e Ciências da Terra, Campina Grande, v. 8, n. 2, p. 183-190, 2008. UNIÃO DA INDÚSTRIA DE CANA-DE-AÇÚCAR-UNICA. Produção de cana-de-açúcar no Brasil. Disponível em: . Acesso em 20 Set. 2009. WATANABE, R. T.; FIORETTO, R. A.; HERMANN, R. Propriedades químicas do solo e produtividade da cana-de-açúcar em função da adição da palhada de colheita, calcário e vinhaça em superfície (sem mobilização). Semina: Ciências Agrárias, Londrina, v. 25, n. 1, p. 93-100, 2004. ZANÃO JÚNIOR, L. A.; LANA, R. M. Q.; CARVALHO-ZANÃO, M. P.; GUIMARÃES, E. C. Variabilidade espacial de atributos químicos em diferentes profundidades em um Latossolo em sistema de plantio direto. Revista Ceres, Viçosa, v. 57, n. 3, p. 429-438, 2010.