Acessibilidade / Reportar erro

ESTIMAÇÃO DE VOLUME DE MADEIRA DE EUCALIPTO POR COKRIGAGEM, KRIGAGEM E REGRESSÃO

PREDICTIONG OF EUCALYPTUS WOOD BY COKRIGING, KRIGING AND REGRESSION

Resumos

No semiárido pernambucano na região do Polo Gesseiro do Araripe, onde é produzido 97% do gesso consumido no Brasil, um experimento florestal com 1875 árvores de eucaliptos foi cortado e todas as árvores foram cubadas rigorosamente pelo método de Smalian. A localização de cada árvore foi demarcada num plano cartesiano, e uma amostra de 200 árvores foi retirada pelo processo inteiramente aleatório. Nas 200 unidades amostrais foram utilizados três métodos de estimação para a variável volume de madeira, análise de regressão, krigagem e cokrigagem. Para a cokrigagem a variável secundária foi o DAP (Diâmetro a Altura do Peito) e para a regressão foi utilizado o modelo de Spurr ou da variável combinada que utiliza duas variáveis explicativas a altura total da árvore (H) e o DAP. As variáveis volume e DAP apresentaram dependência espacial. Para comparar os métodos foram utilizados o coeficiente de determinação (R2), a distribuição residual dos erros (valores observados x estimados). Os melhores resultados foram conseguidos com a regressão R2 = 0,82, volume total = 166,25 m3, a cokrigagem forneceu R2 = 0,72 e volume total = 164,14 m3 e a krigagem teve R2 = 0,32 e o volume total = 163,21 m3, o volume real do plantio era de 166,14 m3.

Inventário Florestal; Volume de Madeira; Geoestatística


In the Gypsum Pole of Araripe, semiarid zone of Pernambuco, where is produces 97% of the plaster consumed in Brazil, a forest experiment with 1875 eucalyptus was cut off and all the trees were rigorously cubed by the Smalian method. The location of each tree was marked on a Cartesian plane, and a sample of 200 trees was removed by entirely random process. In the 200 sample trees, three estimation methods for variable volume timber, regression analysis, kriging and cokriging were used. To cokriging method, the secondary variable was the DBH (Diameter at Breast Height), and for the regression model of Spurr or the combined variable, it uses two explanatory variables: total height of the tree (H) and the DBH. The variables volume and DBH showed spatial dependency. To compare de methods it was used the coefficient of determination (R2) and the residual distribution of the errors (real x estimated data). The best results were achieved with the Spurr equation R2= 0.82 and total volume estimated 166.25 m3. The cokriging provided and R2 = 0.72 with total volume estimated of 164.14 m3 and kriging had R2 = 0.32 and the total volume estimated of 163.21 m3. The real volume of the experiment was 166.14 m3.

Forest Inventory; Wood Volume; Geostatistics


INTRODUÇÃO

A estimação do volume de madeira de um plantio de eucalipto é um processo que demanda tempo, trabalho e dinheiro e que depende entre outros fatores do número de amostras retiradas. A geoestatística surge como uma possibilidade de se realizar boas estimações com um número menor de coletas de dados em campo do que as técnicas que usam a estatística clássica e tendo a vantagem de localizar espacialmente o volume de madeira dentro de um plantio, (Carvalho et al. 2002CARVALHO, J. R. P. de; SILVEIRA, P. M. da; VIEIRA, S. R. Geoestatística na determinação da variabilidade espacial de características químicas do solo sob diferentes preparos. Pesq. Agropec. Bras, Brasília, v. 37, n. 8, p. 1151 - 1159, Ago. 2002.).

Na região do Araripe - PE encontra-se o pólo gesseiro do Araripe, um pólo industrial que concentra uma grande quantidade de indústrias gesseiras, que consomem muita energia. Essa energia em sua grande maioria é proveniente da derrubada da Caatinga (Araújo, 2004ARAÚJO, S. M. S. O pólo gesseiro do Araripe: Unidades geo-ambientais e impactos da mineração. 2004. 276 p. Tese (Doutorado em Ciências, Área de Administração e Política de Recursos Minerais), Universidade Estadual de Campinas, Campinas - SP2004.. ). Uma alternativa possível para essa demanda energética é o plantio de eucalipto para fornecimento de lenha. A determinação do volume produzido pela floresta plantada, segundo Leite e Andrade (2002LEITE, H. G.; ANDRADE, V. C. L. Um método para condução de inventários florestais sem o uso de equações volumétricas., R. ÁrvoreViçosa-MG, v.26, n.3, p.321-328, 2002.), é necessária para se planejar o corte e a comercialização do produto final.

A necessidade de estimação do volume produzido pela floresta plantada segundo Leite e Andrade, (2002LEITE, H. G.; ANDRADE, V. C. L. Um método para condução de inventários florestais sem o uso de equações volumétricas., R. ÁrvoreViçosa-MG, v.26, n.3, p.321-328, 2002.) é informação necessária para se planejar o corte e a comercialização do produto final.

Silva (2008SILVA, S. M. F. S. da. Comparação entre equações volumétricas regionais e equações baseadas em volumes da primeira tora em clones de eucalyptus na chapada do Araripe- PE. 2008, 60 p. Dissertação (Mestrado em Ciências Florestais). Universidade Federal Rural de Pernambuco. Recife - PE. Brasil.) comparou seis modelos de regressão para volume de eucalipto na região do Araripe no mesmo plantio em que a atual pesquisa foi desenvolvida, ele considerou como o melhor deles o modelo de Silva et al. (1993SILVA, J.A.A.; MEUNIER, I.M.J.; BORDERS, B.E.; FARIAS, G.G.A.; ASSUNÇÃO, E. P. Equação volumétrica para Eucalyptus camaldulensis, na região de Barbalha, Ceará, usando o volume da primeira tora como variável independente. Revista Árvore, v.17, p.30-37, 1993b.) que utiliza quatro variáveis explicativas. Guimarães et al. (1996GUIMARÃES, D. P. e LEITE, H. G. Influência do número de árvores na determinação de equação volumétrica para Eucalyptus grandis. Rev. Scientia Forestalis. n. 50, p. 37-42, dez, 1996.) mostraram a influência do número de amostras em equações volumétricas para o eucalipto. Leite e Andrade (2003LEITE, H. G. e ANDRADE, V. C. L. de. Importância das variáveis altura dominante e altura total em equações hipsométricas e volumétricas. R. Árvore, Viçosa-MG, v.27, n.3, p.301-310, 2003.) afirmam que a altura dominante deve ser considerada em estimativas volumétricas de plantios de eucaliptos e que as equações devem ser regionalizadas.

Já a geoestatística usa geralmente, duas técnicas básicas, descritas por Landin, (2002)LANDIM, P. M. B.; STURARO, J. R.; MONTEIRO, R. C. Exemplos de Aplicação de cokrigagem, UNESP / Campus de Rio Claro. Departamento de Geologia Aplicada - IGCE. Texto Didático - 09. 2002., a krigagem, em que a própria variável resposta depende apenas dela mesma e da distância entre suas medições. E a cokrigagem, que além da variável resposta (variável primaria, de difícil medição) utiliza outra variável explicativa (variável secundária, de fácil medição e fortemente correlacionada com a primária) além de utilizar a distância entre suas medições.

No uso da geoestatística é condição impositiva que exista dependência espacial, vários trabalhos mostram que existe dependência espacial para alguns atributos florestais tais como DAP, altura, volume, área basal etc. Rufino et al. (2006RUFINO, T.M.C; THIERSCH, C.R.; FERREIRA, S.O.; KANEGAE JÚNIOR, H. e FAIS, D. Uso da Geoestatística no Estudo da Relação entre Variáveis Dentrométricas de Povoamentos de Eucalyptus sp. e Atributos do Solo. Ambiência Guarapuava, PR v.2 p. 83-93 Edição Especial 1, 2006.) verificaram dependência espacial para altura dominante, área basal, DAP e número de fustes. Mello et al. (2005MELLO, J. M. de.; BATISTA, J.L.F.; RIBEIRO JÚNIOR, P.J. e OLIVEIRA, M.S. de. Ajuste e seleção de modelos espaciais de semivariograma visando à estimativa volumétrica de Eucalyptus grandis., Rev. Scientia Forestalisn. 69. p. 25-37. dez/2005b.b) encontrou dependência espacial para o volume de madeira de eucalipto.

Para a estimativa de volume de madeira de eucalipto encontram-se as condições ótimas para o uso da cokrigagem, pois é de conhecimento geral em ciências florestais a existência de forte correlação entre o DAP (variável de fácil medição) e o volume de madeira (variável de difícil medição). Lundgren (2011LUNDGREN, W. J. C. Uso da geoestatística na estimativa volumétrica de clones de eucalyptus sp no Pólo Gesseiro do Araripe. 2011. 165 p. Tese (Doutorado em Ciências Florestais) - Universidade Federal Rural de Pernambuco, Recife, Disponível em: http://www.bibliotecaflorestal.ufv.br/bitstream/handle/123456789/4729/Tese_Wellington-Jorge-Cavalcanti-Lundgren.pdf?sequence=1. Acesso em 28/09/2014.
http://www.bibliotecaflorestal.ufv.br/bi...
) encontrou forte correlação entre essas duas variáveis, além da já verificada dependência espacial dos dois atributos.

Essa pesquisa mostra que a cokrigagem é uma técnica promissora para o manejo florestal, e que as técnicas geoestatísticas demandam um número menor de medições em campo do que a regressão conseguindo resultados com erros menores que 5%. Nesse trabalho o erro cometido por cada umas das técnicas é conhecido, e não estimado, pois o volume de cada árvore foi medido e também o volume total de madeira existente na floresta.

O objetivo principal é verificar se o DAP fornece estimação precisa de volume de madeira em floresta plantada de eucalipto, além de comparar as técnicas de estimação krigagem, cokrigagem e regressão.

MATERIAL E MÉTODO

A pesquisa foi desenvolvida no município de Araripina - PE que está situado no extremo oeste do estado, no entorno da Chapada do Araripe e faz parte do Pólo Gesseiro do Araripe - PE, possuindo coordenadas geográficas de 07°27'37'' S e 40°24'36'' W e altitude de 831 metros (BARROS, 2009BARROS, B. C. de. Volumetria, densidade, carbono e calorimetria de espécies nativas e exóticas no pólo gesseiro do Araripe - PE. 2009. 67 p. Dissertação (Mestrado em Ciência Florestais). Universidade Federal Rural de Pernambuco, Recife2009..).·.

A precipitação média anual é de aproximadamente 650 mm, concentrando-se entre os meses de dezembro a março em que ocorrem mais de 70% do total anual ficando os demais meses com o restante o que provoca deficiências hídricas que se acumulam a partir do mês de maio até o mês de janeiro. O clima é do tipo BShw, semiárido, quente, com chuvas de verão-outono, pela classificação de Koppen. A temperatura média anual é de 24°C (ARAÚJO, 2004ARAÚJO, S. M. S. O pólo gesseiro do Araripe: Unidades geo-ambientais e impactos da mineração. 2004. 276 p. Tese (Doutorado em Ciências, Área de Administração e Política de Recursos Minerais), Universidade Estadual de Campinas, Campinas - SP2004.. ). Segundo EMBRAPA (2007)EMBRAPA, Síntese do Levantamento de Reconhecimento de Baixa e Média Intensidade dos Solos da Folha Crato. Comunicado técnico 48. dez/2007. Disponível em: http://www.cnps.embrapa.br/solosbr/pdfs/comtec48_2007folha_crato.pdf. Acesso em: 11/03/2011.
http://www.cnps.embrapa.br/solosbr/pdfs/...
, o solo é do tipo Latossolo Vermelho-Amarelo.

No plantio o espaçamento utilizado foi de 2 x 3 m em 4 blocos, cada bloco foi plantado inicialmente com 735 árvores. Os blocos estavam separados entre si por uma distância de 5 metros, cada bloco continha inicialmente 15 (quinze) parcelas constituídas de um único clone de 7 X 7 árvores, totalizando uma área por parcela de 216 metros quadrados. Os clones de eucaliptos foram fornecidos pela empresa Suzano, os clones utilizados foram: E. tereticornis (Clone: C49), Híbrido de E. urophylla x E. tereticornis (Clones: C51, C158, C80, C315, C101, C78 e C156), E. urophylla (Clones: C39, C41 e C33), E. brassiana (Clones: C27, C31, C25 e C11).

O sistema de coordenadas foi definido tendo o eixo X na direção do espaçamento de 3 metros indo de 0 a 174 metros e o eixo Y na direção do espaçamento de 2 metros indo de 0 a 138 metros. Esquema apresentado na Figura 1.

Figura 1:
Posição das árvores do plantio no plano cartesiano com 1875 árvores (Esquerda) e posição da amostra com 200 árvores (Direita). |Figure 1: Position tree planting in the Cartesian plane trees in 1875 (Left,) and position of the sample with 200 trees (Right).

Cada árvore no plantio foi cubada rigorosamente pelo método de Smalian: V = ((A1 + A2)/2) x L, em que V = volume da tora, A1 = Área da base maior da tora, A2 = Área da base menor da tora e L = comprimento da tora. Cada árvore foi derrubada com motosserra a uma altura de 20 cm do solo e medidos os diâmetros de 30 em 30 cm até a altura de 2,10 m, a partir dessa altura os diâmetros foram medidos de metro em metro com suta digital. Um total de 1875 árvores foram cubadas rigorosamente. Conhecendo-se o volume total do povoamento e o volume individual de cada árvore. Em seguida uma amostra sistemática desencontrada de tamanho igual a 200 foi retirada. A amostra cobriu toda a área de forma consistente, permitindo a informação do atributo o que facilitará a construção dos semivariogramas. Para os dois atributos foram testadas as normalidades pelo método de Kolmogorov-Smirnov usando o software GS+ (GAMMA, 2011GAMMA DESIGN SOFTWARE. GS+: Geostatistics for the Environmental Sciences. Gamma Design Software, Plainwell, Michigan USA, 2011. Disponível em: http://www.gammadesign.com/. Acesso em: 11/03/2011.
http://www.gammadesign.com/...
).

Foi realizada a análise de tendência para os dois atributos pelo método apresentado por Camargo et al. (2001CAMARGO, E. C. G.; DRUCK, S.; CÂMARA, G. Análise de Superfícies por Geoestatística Linear. In: Análise Espacial de Dados Geográficos. 2001. Disponível em: http://www.dpi.inpe.br/gilberto/livro/analise/cap3-superficies.pdf. Acesso em: 28/09/2014.
http://www.dpi.inpe.br/gilberto/livro/an...
). Foi colocado no plano cartesiano o eixo X versus os atributos e também o eixo Y versus os atributos, Se o coeficiente de determinação (R2) for baixo podemos considerar a inexistência de tendência. A correlação entre os dois atributos foi calculada pelo método de Pearson. Uma correlação forte é útil tanto para a regressão como para a cokrigagem.

A existência de dependência espacial para os atributos foi verificada pelo semivariograma em que Y*(h) = semivariância amostral, h = distância entre os dois valores medidos, N(h) = número de pares de pontos medidos, Z(xi) = valor medido no ponto xi , Z(xi + h) = valor do atributo medido a uma distância h do ponto Z(xi), xi = local no plano cartesiano (X,Y).

Construído o semivariograma amostral, existem alguns modelos de semivariogramas teóricos que podem ser ajustados a ele, os mais utilizados são:

em que: Y(h) = semivariância teórica, C0 = efeito pepita, h = distância, C = Patamar e A = alcance. O modelo que obtiver maior coeficiente de determinação (R2) será o escolhido. (LUNDGREN, 2011LUNDGREN, W. J. C. Uso da geoestatística na estimativa volumétrica de clones de eucalyptus sp no Pólo Gesseiro do Araripe. 2011. 165 p. Tese (Doutorado em Ciências Florestais) - Universidade Federal Rural de Pernambuco, Recife, Disponível em: http://www.bibliotecaflorestal.ufv.br/bitstream/handle/123456789/4729/Tese_Wellington-Jorge-Cavalcanti-Lundgren.pdf?sequence=1. Acesso em 28/09/2014.
http://www.bibliotecaflorestal.ufv.br/bi...
).

Quando a semivariância se comporta de maneira diferente dependendo da direção em que é observada tem-se a anisotropia, fato que obriga a utilização de diferentes semivariogramas para cada direção. Para a análise da existência de anisotropia, quatro diferentes semivariogramas foram construídos nas direções 0°, 45°, 90° e 135° (LANDIN, 1998LANDIM, P. M. B. Análise Estatística de Dados Geológicos, São Paulo, ed. Fundação Editora da UNESP, 1998, p. 63.). Escolhido o modelo teórico e retirada a tendência se ela existir pode-se estimar o valor do atributo em qualquer local dentro do polígono em estudo usando-se a krigagem descrita a seguir: com a condição de que em que: Z*x0 = valor do atributo estimado no ponto x0 , N = número de pontos, Zxi = valor do atributo no ponto xi, λi = peso da interpolação.

A semivariância cruzada verifica a existência de dependência espacial conjunta para os dois atributos e a sua existência permite a efetivação da cokrigagem (LANDIN et al. 2002LANDIM, P. M. B.; STURARO, J. R.; MONTEIRO, R. C. Exemplos de Aplicação de cokrigagem, UNESP / Campus de Rio Claro. Departamento de Geologia Aplicada - IGCE. Texto Didático - 09. 2002.). A semivariância cruzada é obtida pela equação em que Y*12(h) = semivariância cruzada amostral, Z1(xi) = valor da variável primária medido no ponto xi, Z1(xi + h) = valor da variável primaria medido a uma distância h do ponto Z1(xi), xi = local no plano cartesiano (X,Y), Z2(xi) = valor da variável primária medido no ponto xi, Z2(xi + h) = valor da variável primaria medido a uma distância h do ponto Z2(xi). Escolhido o modelo teórico que melhor se ajusta a semivariância cruzada amostral, pode-se utilizar a cokrigagem para estimar o valor do atributo em qualquer ponto dentro do polígono em estudo, a cokrigagem é descrita por Landin (1998) como: com a imposição de que e em que: Z1* (x0) = valor estimado do atributo Z1 no ponto x0, Z1(x1i) = valor atributo Z1 no ponto xi, Z2j(x1j) = valor do atributo Z2 no ponto xj, N = número de pontos, λ1i e λ1j = pesos das interpolações para os atributos Z1 e Z2, x1i e x1j = pontos onde os atributos Z1 e Z2 foram medidos.

Para a construção dos semivariogramas e mapas de krigagem, cokrigagem e estatísticas clássicas foi utilizado o software GS+ da Gamma Desing, (Robertson, 1998ROBERTSON, G.P. GS+: Geostatistics for the environmental sciences. Versão 5.03 Beta, Plainwell, Gamma Design Software, 1998. 152p. ).

Landin (1998)LANDIM, P. M. B. Análise Estatística de Dados Geológicos, São Paulo, ed. Fundação Editora da UNESP, 1998, p. 63. apresenta um avaliador de dependência espacial denominado de ADE = (patamar - efeito pepita)/patamar, em que, Patamar = o valor da semivariância quando a distância deixa de influenciá-la, efeito pepita = microestruturas não captadas pela distância amostrada, por erros de amostragem, e ou erros de análise (VALENTE, 1989VALENTE, J. M. G. P. Geomatemática - Lições de Geoestatística. Ouro Preto. Ed. da Fundação Gorciex. 25 p. v. 8. 1989.). Quando o ADE < 25%, a dependência espacial é considerada fraca; ADE entre 25 e 75% é moderada, e ADE > 75% possui forte dependência espacial.

O modelo de regressão utilizado foi o de Spurr. Esse modelo utiliza duas variáveis explicativas, o DAP e a altura, sendo descrito pela equação V = a + b(DAP2 x H), em que, V = volume, a e b = coeficientes da equação, DAP = diâmetro a altura do peito e H = altura da árvore. Os parâmetros do modelo foram encontrados pelo softwareSTATISTICA (Statsoft 1995). A escolha desse modelo se deve ao fato de que é comumente utilizado para estimação de volume em Ciências Florestais. Nesse trabalho usamos as alturas reais de cada uma das árvores, o que induz a uma melhora considerável para as estimativas do volume pela equação de regressão.

Para a comparação entre as técnicas foi usado o erro real cometido pela estimação para cada árvore usando o gráfico de valores observados versus valores estimados e o valor do coeficiente de correlação (R), e sendo conhecido o volume total real, também foi comparada a melhor estimativa total de volume pela diferença entre valor real total e valor total estimado. Para comparação entre as técnicas krigagem e cokrigagem também foram construídos mapas de erros usando o software GS+, (GAMMA, 2011GAMMA DESIGN SOFTWARE. GS+: Geostatistics for the Environmental Sciences. Gamma Design Software, Plainwell, Michigan USA, 2011. Disponível em: http://www.gammadesign.com/. Acesso em: 11/03/2011.
http://www.gammadesign.com/...
).

RESULTADOS E DISCUSSÃO

Nenhuma das 3 variáveis apresentaram normalidade pelo teste de Kolmogorov-Sminorv. Mello (2004MELLO, J. M. de. Geoestatística aplicada ao inventário florestal. 2004. 110 p. Tese (Doutorado em Recurso Florestais). Escola Superior de Agricultura Luiz de Queiroz. Universidade de São Paulo. Piracicaba. Estado de São Paulo - Brasil. ) afirma que "não é necessário que o atributo tenha distribuição Normal para que as técnicas geoestatísticas sejam utilizadas, porém essa característica é desejável, pois assim ela ganha propriedades estatísticas ótimas". Os valores das estatísticas descritivas das variáveis analisadas tanto para os dados populacionais quanto para os dados amostrais são apresentados na Tabela 1.

Tabela 1:
Estatísticas descritivas para a população e amostra das variáveis volume (V), DAP e altura (H). CV = coeficiente de variação, N = número de dados, D. padrão = Desvio padrão. |Table 1: Descriptive statistics for population and sample of variables volume (V), DBH and height (H). CV = coefficient of variation, N = number of data, D. padrão = standart deviation.

Não houve uma grande diferença entre os valores das estatísticas amostrais e populacionais, indicando que a amostra representa bem a população. O CV do volume foi regular enquanto os CV do DAP e da altura foram baixos. A correlação entre a variável primaria (V) e as variáveis explicativas foram 0,88 para (V - DAP) e 0,71 para (V - H) ambas consideradas altas.

A equação volumétrica para o modelo de regressão foi: V = 0,016021 + 0,00003(DAP2.H). O gráfico dos valores observados versus os estimados pela regressão encontram-se na figura 04. O R2 para a regressão foi de 0,83, o que é considerado alto. Pela regressão a soma dos volumes estimados foi de 166,25 m3 enquanto o volume total verdadeiro é de 166,14 m3 o que representa um erro de 0,11 m3, ou seja, um erro de 0,07 %.

Para a geoestatística os dados dos dois atributos não apresentaram tendência nem anisotropia, podendo ser construído o semivariograma omnidirecional descrito na equação 1. A existência de dependência espacial foi confirmada pelos semivariogramas, tanto para o volume de madeira quanto para o DAP, fato que coincide com os trabalhos de Mello (2004MELLO, J. M. de. Geoestatística aplicada ao inventário florestal. 2004. 110 p. Tese (Doutorado em Recurso Florestais). Escola Superior de Agricultura Luiz de Queiroz. Universidade de São Paulo. Piracicaba. Estado de São Paulo - Brasil. ) e Lima et al. (2010LIMA, C. G. da R.; CARVALHO, M. de P. e; NARIMATSU, K.C.P.; SILVA, M.G.da e QUEIROZ, H.A. de. Atributos físico-químicos de um latossolo do cerrado brasileiro e sua relação com características dendrométricas do eucalipto. Revista Brasileira de Ciências do Solo, vol. 34. p. 163-173, 2010.), o semivariograma cruzado também apresentou dependência espacial. A Figura 2 mostra esses semivariogramas.

Figura 2:
Semivariograma omnidirecional do volume (A), semivariograma omnidirecional do DAP (B) e semivariograma cruzado (V - DAP) (C). |Figure 2:Omnidirecional Semivariogram for volume (A), Omnidirecional semivariogram for DAP (B) and cross semivariogram (V - DAP) (C).

O modelo que apresentou maior coeficiente de determinação (R2) foi o exponencial. As estimativas dos parâmetros para cada um dos semivariogramas experimentais estão na Tabela 2. O ADE pode ser considerado alto para o semivariograma e DAP. O ADE pode ser considerado alto para o semivariograma cruzado e moderado para os semivariogramas omnidirecionais do volume e DAP.

Tabela 2:
Estimativas dos parâmetros dos semivariogramas omnidirecional e cruzado em que Vol = volume, Exp. = Exponencial. | Table 2: Parameters estimates of the omnidirectional and cross semivariograms Vol = volume, Exp = Exponential.

O alcance teórico do volume no semivariograma foi de 230 m. Este resultado corrobora os resultados de Mello (2005a)MELLO J. M. de, BATISTA J. L. F., OLIVEIRA, M. S. de, JÚNIOR, P. J. R. Estudo da dependência espacial de características dendrométricas para Eucalyptus grandis. Cerne, Lavras, v. 11, n. 2, p. 113-126, abr./jun. 2005a. que encontrou alcance para plantio de eucaliptos de 250 m e Mello et al. (2009)MELLO, J. M. de; DINIZ, F.S.; OLIVEIRA, A.D. de; MELLO, C.R. de; SCOLFORO, J.R.S. e ACERBI JÚNIOR, F.W. Continuidade espacial para características dendrométricas (número de fustes e volume) em plantios de Eucalyptus grandis., R. ÁrvoreViçosa-MG, v.33, n.1, p.185-194, 2009. encontrou alcance de 280 m também para plantio de eucalipto. O alcance do DAP foi menor (108 m), Rufino et al.(2006RUFINO, T.M.C; THIERSCH, C.R.; FERREIRA, S.O.; KANEGAE JÚNIOR, H. e FAIS, D. Uso da Geoestatística no Estudo da Relação entre Variáveis Dentrométricas de Povoamentos de Eucalyptus sp. e Atributos do Solo. Ambiência Guarapuava, PR v.2 p. 83-93 Edição Especial 1, 2006.) afirma ter encontrado dependência espacial para o DAP, porém não fornece o alcance do atributo.

A Tabela 3 mostra que a estimativa total da regressão foi de 166,25 m3, a krigagem estimou o volume total em 163,21 m3 e a cokrigagem estimou o volume total em 164,18 m3. Assim a melhor estimativa de volume total foi obtida pela técnica regressão já que o volume real era de 166,14 m3.

Tabela 3:
Estimativas e erros para as três técnicas. | Table 3:Estimates and errors for the three techniques.

Os mapas da krigagem e cokrigagem para o volume foram construídos e estão na Figura 3. Os dois mapas de volume apresentaram estrutura espaciais semelhantes. Visualmente percebe-se uma maior suavidade nas linhas de contorno para o mapa da cokrigagem.

Figura 3:
Mapas da krigagem (A) e desvios padrão da krigagem (B), cokrigagem (C) e desvios padrão da cokrigagem (D). |Figure 3: Maps kriging (A) and kriging standard deviations (above and right), cokriging (C) and cokriging standard deviations (D).

Os mapas de desvios mostram que a cokrigagem conseguiu resultados mais confiáveis que a krigagem. A Figura 3 apresenta esses mapas

A Figura 4 mostra os gráficos dos valores observados versus os estimados pela regressão, krigagem e cokrigagem, o R2 da cokrigagem foi maior do que o da krigagem e o da regressão foi maior do que o da cokrigagem, a informação adicional do DAP melhorou a estimativa dos volumes para a construção dos mapas, a regressão forneceu o melhor resultado.

Figura 4:
Gráficos dos valores observados x estimados para regressão (A), krigagem (B) e cokrigagem (C). |Figure 4: Plots of the estimated x observed values estimated for regression (A), kriging (B) and cokriging (C).

A regressão foi a técnica que forneceu menores erros, porém a regressão tem a exigência de que para a estimação do volume de uma árvore ainda em campo, seja necessário que se vá ao campo novamente e se realize as medições ou estimativas das variáveis explicativas. Na geoestatística esse trabalho não é mais necessário, para se estimar o volume de uma árvore ainda no campo, é apenas necessário se conhecer sua localização no plano cartesiano.

CONCLUSÕES

1 - A Regressão foi a melhor técnica para estimativa de volume.

2 - A krigagem forneceu um baixo valor de R, porém a estimativa total teve um erro percentual menor que 5% valor comumente considerado baixo para inventários florestais, apesar do R baixo, o mapa da krigagem mostrou-se útil para previsões localizadas geograficamente.

3 - A cokrigagem foi melhor que a krigagem na estimação do volume e pode ser usada em inventários florestais.

4 - O DAP mostrou -se uma variável explicativa útil tanto para estimativas realizadas pela regressão como pela cokrigagem.

  • ARAÚJO, S. M. S. O pólo gesseiro do Araripe: Unidades geo-ambientais e impactos da mineração. 2004. 276 p. Tese (Doutorado em Ciências, Área de Administração e Política de Recursos Minerais), Universidade Estadual de Campinas, Campinas - SP2004..
  • BARROS, B. C. de. Volumetria, densidade, carbono e calorimetria de espécies nativas e exóticas no pólo gesseiro do Araripe - PE. 2009. 67 p. Dissertação (Mestrado em Ciência Florestais). Universidade Federal Rural de Pernambuco, Recife2009..
  • CAMARGO, E. C. G.; DRUCK, S.; CÂMARA, G. Análise de Superfícies por Geoestatística Linear. In: Análise Espacial de Dados Geográficos. 2001. Disponível em: http://www.dpi.inpe.br/gilberto/livro/analise/cap3-superficies.pdf. Acesso em: 28/09/2014.
    » http://www.dpi.inpe.br/gilberto/livro/analise/cap3-superficies.pdf
  • CARVALHO, J. R. P. de; SILVEIRA, P. M. da; VIEIRA, S. R. Geoestatística na determinação da variabilidade espacial de características químicas do solo sob diferentes preparos. Pesq. Agropec. Bras, Brasília, v. 37, n. 8, p. 1151 - 1159, Ago. 2002.
  • EMBRAPA, Síntese do Levantamento de Reconhecimento de Baixa e Média Intensidade dos Solos da Folha Crato. Comunicado técnico 48. dez/2007. Disponível em: http://www.cnps.embrapa.br/solosbr/pdfs/comtec48_2007folha_crato.pdf. Acesso em: 11/03/2011.
    » http://www.cnps.embrapa.br/solosbr/pdfs/comtec48_2007folha_crato.pdf
  • GAMMA DESIGN SOFTWARE. GS+: Geostatistics for the Environmental Sciences. Gamma Design Software, Plainwell, Michigan USA, 2011. Disponível em: http://www.gammadesign.com/. Acesso em: 11/03/2011.
    » http://www.gammadesign.com/
  • GUIMARÃES, D. P. e LEITE, H. G. Influência do número de árvores na determinação de equação volumétrica para Eucalyptus grandis. Rev. Scientia Forestalis. n. 50, p. 37-42, dez, 1996.
  • LANDIM, P. M. B. Análise Estatística de Dados Geológicos, São Paulo, ed. Fundação Editora da UNESP, 1998, p. 63.
  • LANDIM, P. M. B.; STURARO, J. R.; MONTEIRO, R. C. Exemplos de Aplicação de cokrigagem, UNESP / Campus de Rio Claro. Departamento de Geologia Aplicada - IGCE. Texto Didático - 09. 2002.
  • LEITE, H. G. e ANDRADE, V. C. L. de. Importância das variáveis altura dominante e altura total em equações hipsométricas e volumétricas. R. Árvore, Viçosa-MG, v.27, n.3, p.301-310, 2003.
  • LEITE, H. G.; ANDRADE, V. C. L. Um método para condução de inventários florestais sem o uso de equações volumétricas., R. ÁrvoreViçosa-MG, v.26, n.3, p.321-328, 2002.
  • LIMA, C. G. da R.; CARVALHO, M. de P. e; NARIMATSU, K.C.P.; SILVA, M.G.da e QUEIROZ, H.A. de. Atributos físico-químicos de um latossolo do cerrado brasileiro e sua relação com características dendrométricas do eucalipto. Revista Brasileira de Ciências do Solo, vol. 34. p. 163-173, 2010.
  • LUNDGREN, W. J. C. Uso da geoestatística na estimativa volumétrica de clones de eucalyptus sp no Pólo Gesseiro do Araripe. 2011. 165 p. Tese (Doutorado em Ciências Florestais) - Universidade Federal Rural de Pernambuco, Recife, Disponível em: http://www.bibliotecaflorestal.ufv.br/bitstream/handle/123456789/4729/Tese_Wellington-Jorge-Cavalcanti-Lundgren.pdf?sequence=1. Acesso em 28/09/2014.
    » http://www.bibliotecaflorestal.ufv.br/bitstream/handle/123456789/4729/Tese_Wellington-Jorge-Cavalcanti-Lundgren.pdf?sequence=1
  • MELLO, J. M. de. Geoestatística aplicada ao inventário florestal. 2004. 110 p. Tese (Doutorado em Recurso Florestais). Escola Superior de Agricultura Luiz de Queiroz. Universidade de São Paulo. Piracicaba. Estado de São Paulo - Brasil.
  • MELLO J. M. de, BATISTA J. L. F., OLIVEIRA, M. S. de, JÚNIOR, P. J. R. Estudo da dependência espacial de características dendrométricas para Eucalyptus grandis. Cerne, Lavras, v. 11, n. 2, p. 113-126, abr./jun. 2005a.
  • MELLO, J. M. de.; BATISTA, J.L.F.; RIBEIRO JÚNIOR, P.J. e OLIVEIRA, M.S. de. Ajuste e seleção de modelos espaciais de semivariograma visando à estimativa volumétrica de Eucalyptus grandis., Rev. Scientia Forestalisn. 69. p. 25-37. dez/2005b.
  • MELLO, J. M. de; DINIZ, F.S.; OLIVEIRA, A.D. de; MELLO, C.R. de; SCOLFORO, J.R.S. e ACERBI JÚNIOR, F.W. Continuidade espacial para características dendrométricas (número de fustes e volume) em plantios de Eucalyptus grandis., R. ÁrvoreViçosa-MG, v.33, n.1, p.185-194, 2009.
  • ROBERTSON, G.P. GS+: Geostatistics for the environmental sciences. Versão 5.03 Beta, Plainwell, Gamma Design Software, 1998. 152p.
  • RUFINO, T.M.C; THIERSCH, C.R.; FERREIRA, S.O.; KANEGAE JÚNIOR, H. e FAIS, D. Uso da Geoestatística no Estudo da Relação entre Variáveis Dentrométricas de Povoamentos de Eucalyptus sp. e Atributos do Solo. Ambiência Guarapuava, PR v.2 p. 83-93 Edição Especial 1, 2006.
  • SILVA, J.A.A.; MEUNIER, I.M.J.; BORDERS, B.E.; FARIAS, G.G.A.; ASSUNÇÃO, E. P. Equação volumétrica para Eucalyptus camaldulensis, na região de Barbalha, Ceará, usando o volume da primeira tora como variável independente. Revista Árvore, v.17, p.30-37, 1993b.
  • SILVA, S. M. F. S. da. Comparação entre equações volumétricas regionais e equações baseadas em volumes da primeira tora em clones de eucalyptus na chapada do Araripe- PE. 2008, 60 p. Dissertação (Mestrado em Ciências Florestais). Universidade Federal Rural de Pernambuco. Recife - PE. Brasil.
  • STATSOFT. STATISTICA for Windows [Computer Program Manual]. Inc. (1995). Tulsa, OK: StatSoft, Inc. 2325 East 13th Street, Tulsa, OK, 74104.
  • VALENTE, J. M. G. P. Geomatemática - Lições de Geoestatística. Ouro Preto. Ed. da Fundação Gorciex. 25 p. v. 8. 1989.

Datas de Publicação

  • Publicação nesta coleção
    Apr-Jun 2015

Histórico

  • Recebido
    07 Maio 2012
  • Aceito
    21 Jan 2015
UFLA - Universidade Federal de Lavras Universidade Federal de Lavras - Departamento de Ciências Florestais - Cx. P. 3037, 37200-000 Lavras - MG Brasil, Tel.: (55 35) 3829-1706, Fax: (55 35) 3829-1411 - Lavras - MG - Brazil
E-mail: cerne@dcf.ufla.br