RELAÇÃO ENTRE ATRIBUTOS FÍSICOS E QUÍMICOS DOS SOLOS E A PRODUTIVIDADE DE CAPIM DOURADO NA REGIÃO DO JALAPÃO, TO

RELATIONSHIP BETWEEN SOIL PHYSICAL AND CHEMICAL PROPERTIES AND YIELD OF GOLDEN GRASS IN THE JALAPÃO REGION, TO, BRAZIL

Resumos

Na região do Jalapão, no Estado do Tocantins, existem grandes áreas contínuas de Neossolos Quartzarênicos, com desenvolvimento do capim dourado (Syngonanthus nitens), espécie nativa utilizada na confecção de artesanatos nessa região, possuindo destacada importância na economia local. O objetivo deste estudo foi verificar a relação entre os atributos físicos e químicos dos solos e a produtividade de capim dourado nos Campos Limpos Úmidos associados às Veredas no Parque Estadual do Jalapão. Foram coletadas amostras de solo e escapos de capim dourado e realizadas medições de declividade do terreno ao longo de caminhamento ou transeções perpendiculares às linhas de drenagem das Veredas selecionadas para o estudo. Foram identificados três tipos principais de solos ao longo das transeções: Neossolo Quartzarênico órtico típico, Neossolo Quartzarênico hidromórfico típico e Organossolo Háplico fíbrico típico. A produtividade de capim dourado foi analisada por meio de modelos de regressão múltipla, envolvendo 26 variáveis dos solos. Em razão da significância, o modelo final ficou constituído por seis dessas variáveis: pH, Mg²⁺, Ca²⁺, areia fina, argila e umidade do solo. O modelo de regressão quadrática múltipla envolvendo essas seis variáveis contemplou R² de 0,85 e R² ajustado de 0,81 da relação avaliada, onde a umidade dos solos foi a variável mais importante.

Syngonanthus nitens; vereda; regressão múltipla; solos arenosos

In the Jalapão region, State of Tocantins, Brazil, there are large contiguous areas of Neossolos Quartzarênicos (Tropsamments). Golden grass (Syngonanthus nitens) is a native species widely used as material for handicrafts in the region, and it is therefore of great importance for the local economy. The aim of this study was to analyze the relationship between soil physical and chemical properties and golden grass yield in Campos Limpos Úmidos within Jalapão State Park (Parque Estadual do Jalapão). We collected soil samples and golden grass scapes and made slope measurements along the transects perpendicular to the drainage lines. We identified three soil units: Neossolo Quartzarênico órtico típico, Neossolo Quartzarênico hidromórfico típico, and Organossolo Háplico fíbrico típico (Orthpsamment, Hydropsamment, and Hydrofibrist, respectively). Multiple regression models involving 26 variables further portrayed golden grass yield. Six of these variables subsequently proved to be of greatest significance and constituted the final model: pH, Mg²⁺, Ca²⁺, fine sand, clay, and soil moisture. The multiple quadratic regression model involving these six variables considered an R² of 0.85 (adjusted R²of 0.81) of the relationship evaluated. Soil moisture was the most important variable.

Syngonanthus nitens; multivariate regression; sandy soils

INTRODUÇÃO

Os Neossolos Quartzarênicos (RQ) ocupam 30 milhões de hectares no bioma Cerrado (Caetano et al., 2013); no entanto, observam-se poucos estudos envolvendo a caracterização desses solos sob ambientes nativos. No Estado do Tocantins, os RQ ocupam aproximadamente 18,52 % do seu território (Seplan, 2012).

Na região do Jalapão, Estado do Tocantins, ocorrem grandes áreas contínuas com desenvolvimento de RQ, que apresentam intensa rede de drenagem e por causa das características físicas e químicas dessa classe de solo, além do relevo plano, ocorre o desenvolvimento de Campos Limpos Úmidos (Santana et al., 2014). A cobertura vegetal natural na região encontra-se preservada (Santiago e Pereira Júnior, 2010), ou seja, verifica-se baixa conversão da vegetação natural em uso antrópico. Em razão dessa preservação, da riqueza em biodiversidade, das limitações do uso agrícola, das fragilidades dessa classe de solo, entre outras variáveis, foram criadas diversas unidades de conservação, o que torna a região um local propício para pesquisas com espécies de ocorrências naturais do Cerrado.

Nos locais onde os RQ predominam, na região do Jalapão, destaca-se o extrativismo do capim dourado, Syngonanthus nitens (Bong.) Ruhland, da família Eriocaulaceae. O capim dourado é uma sempre-viva, ou seja, mantém a forma e a coloração dos escapos e as inflorescências inalteradas após serem coletados, como se tais estruturas estivessem vivas (Giulietti et al., 1996); essas características fazem dessa espécie a matéria-prima para a produção de artesanato. Além da ocorrência do capim dourado nos Campos Limpos Úmidos da região do Jalapão, TO (Schmidt et al., 2007; Schmidt e Ticktin, 2012;Santana et al., 2014), a ocorrência dessa espécie ainda é relatada na Chapada dos Veadeiros, GO, no Distrito Federal (Munhoz e Felfili, 2007) e no Triângulo Mineiro, MG (Ramos et al., 2006), e está sempre associada à fitofisionomia de Veredas. De acordo com Carvalho (1991), as Veredas constituem subsistemas úmidos que participam do controle do fluxo do lençol freático, desempenhando papel fundamental no equilíbrio hidrológico dos cursos d’água no ambiente do Cerrado. Schmidt e Ticktin (2012), Pacifico et al. (2011), Schmidt et al. (2011) caracterizaram o capim dourado como um importante produto de exploração não florestal, no ambiente de Veredas na região do Jalapão.

Estudos relacionados à produtividade vegetal e aos atributos físicos ou químicos do solo são comuns em espécies cultivadas; entretanto, tais estudos em espécies de ocorrência natural são pouco comuns, em razão principalmente ao grau de dispersão das espécies que ocorrem em mosaicos vegetacionais, como acontece no bioma Cerrado, sem que haja padrão definido de ocorrência da espécie. Os estudos de produtividade em espécies de ocorrência natural no Cerrado têm-se concentrado nas espécies frutíferas, tanto para consumo in natura como para industrialização; são os casos do baru (Dipteryx alata) (Sano e Simon, 2008) e do araticum (Annona spp.) (Mesquita et al., 2007), entre outros.

Diante do exposto, este trabalho teve como objetivo verificar a relação entre os atributos físicos e químicos dos solos e a produtividade de capim dourado nos Campos Limpos Úmidos associados às Veredas no Parque Estadual do Jalapão, TO.

MATERIAL E MÉTODOS

A área de estudo situa-se no Parque Estadual do Jalapão (PEJ), no município de Mateiros, TO, entre as latitudes 10º 08’ 73” S e 10º 35’ 09” S e as longitudes 47º 04’ 01” W e 47º 35’ 06” W. O PEJ foi criado pela Lei Estadual nº. 1.203, de 12 de janeiro de 2001, possui 158.885 ha e ocupa 16 % da área municipal de Mateiros. No interior do PEJ, ocorrem várias nascentes que alimentam a bacia do rio Sono e o sistema hidrográfico Tocantins-Araguaia. O relevo dominante é plano a suave ondulado, com altitudes variando de 400 a 500 m. Destacam-se as escarpas erosivas representadas pelas serras do Espírito Santo, Jalapinha e Porco, que atingem altitudes entre 770 e 798 m. Quanto aos solos da região, predominam os Neossolos Quartzarênicos órticos, originados de rochas sedimentares areníticas do Grupo Urucuia (Seplan, 2012).

O PEJ abriga diferentes formações vegetacionais savânicas, com predominância de espécies campestres (Sano et al., 2008; Santiago e Pereira Júnior, 2010). O clima na região, segundo o critério de Thornthwaite, é classificado como C1dA’a’, que corresponde a clima subúmido seco com moderada deficiência hídrica no inverno e precipitação pluvial média anual entre 1.500 e 1.600 mm. O período chuvoso compreende os meses de novembro a abril, enquanto o período de seca inicia-se em maio, terminando em outubro. A temperatura média anual é de 27 ºC (Seplan, 2012).

Para o estudo, foram selecionadas 10 Veredas no PEJ, consideradas representativas do Parque em relação às suas dimensões e proporção das suas áreas de inundação, proporcionando o desenvolvimento de classes de solos similares e ocorrência de vegetação nativa preservada com ocorrência do capim dourado. Essas informações foram extraídas por meio de análises prévias em fotografias aéreas, na escala de 1:25.000, obtidas em novembro de 2002 pelo Projeto AASA 92/02 (Naturatins, 2002) e de imagens do satélite RapidEye de 16 de setembro e 14 de outubro de 2010, utilizadas por Santana et al. (2014).

Como o mês de setembro representa a época em que o capim dourado atinge a plena maturação (Naturatins, 2007), as campanhas de campo para a coleta de amostras para este estudo foram realizadas de 27 agosto a 1º de setembro de 2011, período sem precipitação pluvial. Durante as atividades de campo, foram coletadas amostras de solos, de capim dourado, além de medições de declividades dos terrenos ocupados pelos Campos Limpos Úmidos associados às Veredas avaliadas.

Para a amostragem, nas Veredas selecionadas para o estudo, foram estabelecidos três caminhamentos ou transeções perpendiculares às linhas de drenagem, com 100 m de espaçamento entre elas. As transeções partiam-se das áreas mais afastadas da rede de drenagem, onde se verificava o início da ocorrência de Campos Limpos Úmidos, até as proximidades da área inundada, próxima à rede de drenagem, marcada pela presença da palmácea buriti (Mauritia flexuosa).

As transeções foram divididas em três segmentos de mesma extensão em relação ao maior comprimento de cada uma delas, denominadas de S1 (segmento mais afastado da rede de drenagem, com ocorrência de solos mais secos), S2 (segmento intermediário, com solos com umidade intermediária) e S3 (segmento próximo à rede de drenagem, com solos mais úmidos), como demonstrado na figura 1.

Figura 1
Esquema amostral de solos, capim dourado e de medições de declividade do terreno nas Veredas selecionadas para o estudo. T: transeções; S: segmentos; MG: mata de galeria; e CLU: campo limpo úmido.

Em cada segmento, foram coletadas amostras indeformadas e deformadas de solos na profundidade de 0-10 cm, ambas com três repetições; amostras aleatórias da população de capim dourado em uma área de 0,5 × 0,5 m, com três repetições. Foram feitas leituras de diferença de nível para estimar a declividade do terreno, utilizando-se de uma mangueira de plástico transparente (Guimarães et al., 2002), anotando-se a diferença de cota a cada posição no terreno no sentido de S1 para S3. Assim, ao final da amostragem, para as 10 Veredas selecionadas totalizaram-se 90 amostras de solos (deformadas e indeformadas), 90 amostras de capim dourado e 90 medições de declividade (três transeções × três segmentos).

Com as amostras indeformadas dos solos, determinaram-se a umidade volumétrica e densidade dos solos (Ds). As amostras deformadas foram secas ao ar, destorroadas e passadas em peneiras de malha de 2 mm, para obter a fração terra fina seca ao ar (TFSA). Com essas amostras, foram realizadas análises granulométricas (areia fina, areia grossa, areia total, argila e silte) e químicas. As análises químicas constaram de determinações de pH em água; acidez potencial (H+Al) (Acetato de cálcio 0,5 mol L^-1, pH 7,0); P e K (extrator Mehlich-1); Ca²⁺ e Mg²⁺ (KCl 1 mol L^-1). A matéria orgânica (MO) foi extraída por via úmida pelo método Walkey-Black e o carbono orgânico (CO) foi obtido pelo equivalente de 58 % da MO. Foram calculados os valores de somas de bases (S); capacidade de troca catiônica (CTC) (T); CTC efetiva (t); saturação por bases (V) e saturação por alumínio (m), por meio das amostras deformadas. As análises físicas e químicas dos solos seguiram métodos descritos emEmbrapa (1997). Todas as análises laboratoriais de solos foram realizadas nos Laboratórios de Física e Química de Solos da Embrapa Cerrados.

Os escapos coletados de capim dourado foram acondicionados em sacos de papel e posteriormente foram encaminhados ao laboratório, onde foram secos em estufa de circulação forçada de ar a uma temperatura de 55 °C por 48 h ou até peso constante. Em seguida, foram pesados em balança analítica de precisão (quatro casas decimais) no Laboratório de Física de Solos da Embrapa Cerrados.

Dessa forma, foram avaliadas 23 variáveis envolvendo atributos físicos e químicos dos solos e declividade do terreno. De posse desses dados, foi criada uma matriz de dados com as variáveis analisadas mais três variáveis categóricas (segmentos 1, 2 e 3), para investigar e modelar a relação entre produtividade de escapos de capim dourado (variável dependente) e as variáveis estudadas consideradas independentes, utilizando a análise de regressão linear múltipla.

Para realizar as análises estatísticas, foi utilizado o Statistical Package for Social Sciences (SPSS), versão 19.0 (SPSS, 2008), em que as variáveis independentes utilizadas para a construção do modelo foram selecionadas pelo procedimentobackward disponível nesse programa. Neste estudo, foi ajustado um modelo de regressão polinomial múltipla de segunda ordem, também chamado de regressão quadrática múltipla.

Para definir o modelo, visto que os resultados das variáveis independentes analisadas evidenciaram pequena variação, para cada segmento das transeções foi realizada a média aritmética dos dados. Com isso, reduziu-se o número de amostras de solos, do capim dourado e da declividade de 90 para 30. Utilizou-se a transformação logarítmica da variável dependente Y e estabeleceu-se nova variável dependente ln(Y_i), que corresponde ao “logaritmo natural de Y+1”, a fim de reduzir a assimetria da distribuição e aproximar sua curvatura da distribuição normal. Essa transformação foi utilizada, ainda, para reduzir o problema da heterocedasticidade, ou seja, a distribuição desigual das variâncias dos erros aleatórios da função de regressão, procedimento semelhante ao realizado por Lúcio et al. (2010, 2011), com o diferencial que os autores usaram transformação λy = ln(x + 0,5).

Foi realizada ainda uma análise exploratória para verificar a existência de possíveis valores atípicos (outliers) por meio dos resíduos padronizados maiores que +2 ou menores que -2. A existência de tais dados pode distorcer os resultados do estudo; logo, a sua presença deve ser verificada, identificada e desconsiderada na análise. A análise de variância foi utilizada para avaliar o modelo de regressão linear múltipla pelo método dos mínimos quadrados. O teste t avaliou a significância estatística dos coeficientes de regressão estimados para as variáveis independentes em nível de significância de 5 %. A análise de multicolinearidade foi avaliada no modelo de regressão para verificar a correlação entre duas ou mais variáveis independentes ou entre uma delas e as demais incluídas no modelo, pelo fator de inflação da variância (Mayers, 2000).

Com esses procedimentos, foram avaliados os quatro pressupostos para a regressão (linearidade, independência dos erros, normalidade dos erros e igualdade de variâncias) para verificar a existência de comportamento linear entre as variáveis, a existência de homocedasticidade nos dados, a distribuição normal dos resíduos, a multicolinearidade das variáveis e a ocorrência de outliers.

Para avaliar o grau de desvio do modelo, em relação à linha de ajuste perfeito (linha 1:1), foi utilizado o recurso da regressão linear simples, aplicado aos pares dos valores estimado e observado. A qualidade preditiva do modelo foi avaliada quanto à precisão, à exatidão e ao desempenho. A precisão foi avaliada pelo coeficiente de correlação de Pearson (r) entre a produtividade observada e a estimada, indicando o grau de dispersão dos dados obtidos em relação à média, ou seja, o erro aleatório. A exatidão foi quantificada pelo índice de concordância (d) de Willmott (Willmott, 1982). Essa exatidão está relacionada ao afastamento dos valores estimados em relação aos observados. Seus valores variam de zero (nenhuma concordância) a 1 (concordância perfeita). O desempenho do modelo foi quantificado pelo índice de consistência c = r × d (Camargo e Sentelhas, 1997), que, de acordo com os critérios de interpretação proposto pelos autores, o índice c varia de zero a 1, correspondendo a desempenho péssimo a ótimo, respectivamente.

RESULTADOS E DISCUSSÃO

Segmento S1 das transeções

Os solos desse segmento apresentaram textura arenosa 988,68 g kg^-1 de areia total e 11,31 g kg^-1 de soma dos teores de argila e silte, representando solos com grãos simples, sem estruturação (Quadro 1). Santana et al. (2010) encontraram valores de 896,70 g kg^-1 de areia total e a soma dos teores de argila e silte com 103,30 g kg^-1, em RQ no município de Lagoa do Tocantins, região do Jalapão. Nesse segmento, os solos foram classificados como Neossolos Quartzarênicos órticos típicos, segundo o Sistema Brasileiro de Classificação de Solos (Embrapa, 2013).

Quadro 1
Média das análises físicas e químicas dos solos, declividade do terreno e produtividade de capim dourado dos segmentos S1, S2 e S3 das Veredas avaliadas no PEJ

Atributo	Segmento 1	Segmento 2	Segmento 3
Argila (g kg-1)(1)	4,66	29,66	95,33
Silte (g kg-1)(1)	6,66	24,33	109,66
Areia grossa (g kg-1)(1)	439,33	473,33	273,66
Areia fina (g kg-1)(1)	549,33	472,66	521,33
Areia total (g kg-1)(1)	988,66	945,99	794,99
Densidade do solo - Ds (kg dm-3)(2)	1,49	1,09	0,33
Umidade (cm3 cm-3)(2)	0,03	0,40	0,73
Declividade (%)	5,27	5,36	5,42
pH(H2O)(3)	5,52	5,47	5,22
K+ (cmolcdm-3)(4)	0,01	0,04	0,11
Mg2+ (cmolcdm-3)(5)	0,32	0,35	0,37
Ca2+ (cmolcdm-3)(5)	0,11	0,12	0,16
Al3+ (cmolcdm-3)(5)	0,47	1,02	1,77
H+Al (cmolc dm-3)(6)	2,75	6,69	12,98
P (mg kg-1)(4)	1,36	1,71	4,31
Soma de bases - S (cmolc dm-3)	0,44	0,52	0,64
CTC efetiva - t (cmolc dm-3)	0,91	1,52	2,42
CTC a pH 7 - T (cmolc dm-3)	3,19	7,21	13,63
Saturação por Al - m (%)	51,53	63,10	72,85
Saturação por bases - V (%)	14,72	9,05	4,99
Carbono orgânico - CO (g kg-1)(7)	6,91	46,08	195,06
Matéria orgânica - MO (g kg-1)	4,01	26,72	113,13
Produtividade (g m-2)	1,74	3,36	2,48

• ⁽¹⁾ Método da pipeta (Embrapa, 1997); ⁽²⁾Embrapa (1997); ⁽³⁾ pH em água, relação solo:solução 1:2,5; ⁽⁴⁾ Extrator Mehlich-1; ⁽⁵⁾Extrator KCl 1 mol L^-1; ⁽⁶⁾ Extraído com solução de acetato de cálcio 0,5 mol L^-1 a pH 7,0;⁽⁷⁾ Método Walkey-Black modificado (Embrapa, 1997). CTC: capacidade de troca catiônica.

Nesse segmento, foi encontrado o menor teor de água no solo (0,03 cm³cm^-3). A Ds foi a mais alta em relação aos demais segmentos, com valor médio de 1,49 kg dm^-3. Já o teor de carbono orgânico (CO) foi o menor, com 4,01 g kg^-1. O baixo teor de CO ocorreu em razão da constituição mineral arenosa dos solos, com baixa CTC, que confere baixa fixação da matéria orgânica, além da elevada temperatura média anual da região (27^oC), que também contribui para elevadas taxas de oxidação biológica do CO do solo (Bayer et al., 2000). Santana et al. (2010)também encontraram baixos teores de CO em RQ na região do Jalapão, onde os teores de MO variaram de 6,57 a 7,69 g kg^-1. Carbono orgânico, S, Al³⁺, acidez potencial (H+Al) e, consequentemente, valores de T e CTC efetiva apresentaram menores valores em relação aos demais segmentos. Os solos desse segmento apresentaram baixo pH, classificando-se como medianamente ácido.

Em relação à declividade, não foram verificadas grandes amplitudes entre S1 e S3. A média variou de 5,27 a 5,42 %, classificando-se o relevo como suave ondulado.

A saturação por bases (V) apresentou valores muito inferiores a 50 %, e a saturação por Al (m), maior que 50 %, caracterizando os solos como distróficos e álicos. O teor de P também foi baixo, semelhante aos encontrados por Silva et al. (2010) em solos dessa região do Jalapão. Nesse segmento, a produtividade do capim dourado foi a mais baixa, com 1,74 g m^-2. Essa menor produtividade pode ser explicada pela baixa umidade do solo (Quadro 1). SegundoOliveira et al. (2010), os baixos teores de umidade no solo não propiciam condições favoráveis para a absorção dos nutrientes, que se dá por fluxo de massa ou difusão.

Segmento S2 das transeções

Os solos apresentaram textura arenosa, com o teor de areia total em torno de 946,00 g kg^-1 e com 53,99 g kg^-1 de soma dos teores de argila e silte, sendo classificados como Neossolos Quartzarênicos hidromórficos típicos (Embrapa, 2013). Nesse segmento, foram encontrados os menores e os maiores teores de areia fina e areia grossa, respectivamente, dentre os três segmentos avaliados. A Ds evidenciou valor mediano entre as densidades dos solos nos segmentos S1 e S3, com valor de 1,09 kg dm^-3. O relevo também se enquadrou na classe de suave ondulado (5,36 % de declividade).

O valor do pH foi baixo, menor que em S1; entretanto, continuou sendo classificado como medianamente ácido. O teor de Al³⁺ teve acréscimo, traduzindo-se em maior acidez potencial (H+Al) em relação aos valores encontrados nos solos em S1. A soma de bases continuou baixa, o que implicou em teores de V muito inferiores a 50 % e teor de m superior a 50 %, ou seja, os solos são distróficos e álicos.

Em relação ao teor de P, houve pequeno acréscimo em relação aos solos em S1. Já o teor de CO, teve aumento considerável em relação aos solos em S1 (Quadro 1). O acúmulo de material orgânico foi, provavelmente, ocasionado pelo ambiente anaeróbico nesse segmento, que, ao final da época da seca, atingiu valores médios de umidade de 0,40 cm³cm^-3. Esse alto teor de umidade, no final da época seca, pode indicar que nesse ambiente ocorra saturação do solo na época das chuvas, pois ocupa o terço médio da transeção, uma vez que o ambiente hidromórfico é caracterizado por saturação permanente ou periodicamente (Embrapa, 2013).

A produtividade do capim dourado foi a mais elevada nesse segmento em relação aos demais, atingindo 3,36 g m^-2. A sua alta produtividade pode ser explicada pela umidade encontrada, o que propiciou as condições mais favoráveis ao desenvolvimento e à produção de escapos da espécie em relação aos segmentos S1 e S3. Segundo Schmidt et al. (2008), a emissão dos escapos ocorre nos meses de abril e maio; sua floração inicia-se a partir de julho. Nos meses de abril e maio, inicia-se o período em que cessam as chuvas na região do Jalapão e começa a época da seca.

Segmento S3 das transeções

A textura do solo nesse segmento se enquadrou em areia franca, com teor de areia total em torno de 795,00 g kg^-1 e com teores de argila e silte de 95,33 e 109,66 g kg^-1, respectivamente. Os solos foram classificados como Organossolos Háplicos fíbricos típicos (Embrapa, 2013). O teor de areia grossa foi o menor em relação a S1 e S2 (Quadro 1). Em S3, foi encontrada a maior umidade no solo (0,73 cm³ cm^-3); nessa região, o solo se encontra saturado, em razão da sua posição mais baixa no relevo em relação aos segmentos anteriores, aliada à proximidade com o curso d’água. A Ds foi a menor (0,33 kg dm^-3) entre os segmentos avaliados, uma vez que essa diminui com o incremento de MO, em razão da baixa densidade da MO em relação à fração mineral do solo (Araújo et al., 2007; Valladares, 2008). Esse valor de Ds está dentro da amplitude encontradao por Valladares et al. (2008), em estudo sobre Organossolos e solos afins. A declividade atingiu 5,42 %, continuando na mesma classe de relevo dos demais segmentos (suave ondulado).

Nesse segmento, o solo apresentou o mais baixo valor de pH(H₂O), mas permaneceu classificado como medianamente ácido. O Al³⁺ teve acréscimo considerável em relação aos solos em S1 e S2, refletindo aumento na acidez potencial (H+Al), provavelmente por causa do aumento de CO. A diminuição dos valores de pH(H₂O) e V se deve ao maior teor de CO no solo, pois no processo de mineralização e formação de substâncias húmicas há a liberação de nitratos e íons H⁺ na solução do solo (Meurer, 2006). A soma de bases (S) continuou baixa, o que implicou em baixos teores de V (<50 %) e elevado teor de m (>50 %), classificando-se como solos distróficos e álicos. Ocorreu aumento nos teores de T em relação aos solos em S1 e S2, possivelmente em razão do aumento no teor de CO. Nesse ambiente, foi detectado acúmulo de material orgânico originado de ambiente anaeróbico com umidade elevada, com teor de CO maior que 80,0 g kg^-1, configurando-se como solos com predomínio de material orgânico. O teor de CO encontrado nesse segmento (113,13 g kg^-1) está próximo ao intervalo encontrado por Ebeling et al. (2011), que obtiveram valores entre 116,66 e 745,83 g kg^-1 em Organossolos de diversas regiões do Brasil. Os autores atribuíram a variabilidade em decorrência da cobertura vegetal, do clima e da umidade ao longo do tempo. Em geral, os Organossolos apresentam características como elevada acidez e baixa saturação por bases (Ebeling et al., 2008; Valladares, 2008), como verificado neste trabalho. Foi constatado aumento no teor de P em relação aos teores encontrados em S1 e S2, com valores superiores aos encontrados porSilva et al. (2010) em Veredas da região do Jalapão, possivelmente associado ao incremento de CO. A produtividade do capim dourado foi intermediária (2,48 g m^-2), em relação aos segmentos S1 e S2, o que deve estar associado à saturação de água do solo.

Comportamento das variáveis ao longo dos Segmentos S1, S2 e S3

As médias dos resultados das variáveis avaliadas nas 10 Veredas do PEJ evidenciaram que, para os segmentos S1, S2 e S3, os valores obtidos foram crescentes para K⁺, Mg²⁺, Ca²⁺, P, CO, S, T, Al³⁺, H+Al, t, m, argila, silte, umidade e declividade no sentido de S1 para S3. Já os valores de pH, V, areia total e Ds apresentaram decréscimo.

Já a produtividade do capim dourado apresentou teores maiores em S2.

Relação entre atributos dos solos e a produtividade do capim dourado

Os resultados da análise da regressão quadrática múltipla entre a produtividade de escapos de capim dourado e as variáveis independentes indicaram que o modelo adotado foi estatisticamente significativo, melhorando com a retirada de uma amostra considerada outlier. A variabilidade explicada pelo modelo (soma dos quadrados do modelo - SS_M) foi de 11,815, e a variabilidade não explicada (soma dos quadrados dos resíduos - SS_R) foi de 2,121, correspondendo a uma variação total (SS_T) de 13,935. A comparação entre a média dos quadrados da regressão (MQ_M = 1,969) e a média dos quadrados dos resíduos (MQ_R = 0,096) confirmou que as variações na produtividade não são consequência do acaso (resíduos) e podem ser explicadas pelas diferenças existentes entre os valores das variáveis independentes incluídas no modelo.

A análise de variância apresentou valor de F = 20,429 (significativo a p<0,001), aceitando-se a hipótese de que existe relação linear entre a variável dependente em pelo menos uma das variáveis independentes, ou seja, existe relação de causa-efeito entre a variável dependente e as independentes. O modelo apresentou valor de R² de 0,85 e R² ajustado de 0,81.

Mediante processo de redução do número de variáveis independentes a serem consideradas no modelo, utilizando o método backward, o teste t a 5 % de significância selecionou as variáveis pH (x₁)², Mg²⁺ (x₂)², Ca²⁺(x₃)², areia fina - AF (x₄)², argila (x₅) e umidade (x₆). As estimativas dos coeficientes das variáveis independentes e do intercepto para o modelo de regressão quadrática múltipla encontram-se no quadro 2.

Quadro 2
Coeficientes estimados e padronizados relativos ao teste t, no intervalo de confiança de 95 % e fator inflação da variância (VIF) para os coeficientes da regressão quadrática múltipla para a produtividade de capim dourado nos Campos Limpos Úmidos associados às Veredas avaliadas do Parque Estadual do Jalapão, Tocantins

Variável	Coeficiente			Teste t	Sig.	Intervalo de confiança de 95 % para β		VIF
	Não padronizado		Padronizado			Limite
	β	Erro-padrão	β			Inferior	Superior
X0 = intercepto	6,38	1,42		4,48	0,00	3,43	9,34
(x1)2 = pH × pH	-0,12	0,04	-0,31	-2,65	0,01	-0,21	-0,03	2,01
(x2)2 = Mg × Mg	-9,46	1,73	-0,70	-5,46	0,00	-13,06	-5,87	2,40
(x3)2 = Ca × Ca	-13,47	4,35	-0,42	-3,10	0,01	-22,49	-4,45	2,68
(x4)2 = AF × AF	-2,71×10-6	0,00	-0,61	-5,40	0,00	0,00	0,00	1,86
(x5) = Argila	-0,01	0,00	-0,645	-3,67	0,00	-0,02	-0,00	4,52
(x6) = Umidade	1,94	0,47	0,81	4,13	0,00	0,96	2,91	5,54

• Variável dependente: ln(produtividade); AF: areia fina; e Sig.: significância.

Dos coeficientes estimados, apenas o intercepto X₀ (valor da variável dependente quando todas as variáveis independentes são iguais a zero) e a umidade apresentaram interação positiva com a produtividade de escapos. O coeficiente argila (x₅) evidenciou interação negativa, e as demais variáveis apresentaram interações quadráticas; os valores dos coeficientes são negativos.

Essas interações explicam a adaptação do capim dourado ao ambiente natural de ocorrência do bioma Cerrado no Jalapão, uma espécie nativa em Campo Limpos Úmidos associados às Veredas. Essas interações corroboram o estudo de Haridasan (2008), o qual constatou que espécies nativas que ocorrem nesse bioma são resistentes e tolerantes às condições edáficas consideradas desfavoráveis às plantas cultivadas, mas sua distribuição e frequência nas comunidades nativas, crescimento e produtividade são determinados pela disponibilidade de nutrientes, regime hídrico do solo e outros fatores edáficos.

O ambiente de ocorrência do capim dourado apresentou solos medianamente ácidos e distróficos, que não influenciaram o seu desenvolvimento, concordando com Haridasan (2008). Entretanto, o capim dourado é uma espécie exigente em teor de água no solo.

De acordo com Haridasan (2008), algumas espécies nativas do bioma Cerrado crescem em solos ácidos e são tolerantes ou resistentes ao Al, porque sua capacidade de absorção de nutrientes essenciais, crescimento e reprodução não são prejudicados pelas altas concentrações de Al no solo. Esse fato é observado pela elevação da saturação por Al do S1 em relação ao S2 e mesmo assim a produtividade teve acréscimo próximo a 100 %. O mesmo é comprovado pelo fato de m não fazer parte da equação que explica a relação entre atributos físicos e químicos do solo e a produtividade de capim dourado.

Todos os coeficientes de regressão foram significativos (p<0,05), aceitando-se a hipótese de existência de uma relação linear significativa entre as variáveis independentes consideradas no modelo na contribuição da variável dependente (produtividade) (Quadro 2).

Os coeficientes de β padronizados são todos mensurados em termos de unidades de desvios-padrão e são diretamente comparáveis; a maior contribuição das variáveis independentes no modelo foi dada pela umidade, seguida por pH, Ca, areia fina, argila e Mg. O fator inflação da variância (VIF) encontrado no modelo para as variáveis independentes está de acordo com Mayers (2000), o qual sugeriu que a partir de 10 é um valor preocupante.

Assim, o modelo de regressão quadrática múltipla ajustado para a produtividade de capim dourado encontra-se apresentada na equação 1:

Sendo ŷ: produtividade estimada; e: expoente; (x₀) = 6,383 (intercepto); (x₁)² = -0,116 (pH × pH); (x₂)² = -9,464 (Mg × Mg); (x₃)²= -13,470 (Ca × Ca); (x₄)² = -2,712×10-6 (AF × AF); (x₅) = -0,009 (Argila); e (x₆) = 1,935 (Umidade).

A dispersão dos 29 valores estimados e observados da produtividade do capim dourado no PEJ é evidenciada na figura 2, cuja média calculada por segmento encontra-se no quadro 1.

Figura 2
Diagrama de dispersão das amostras da produtividade do capim dourado observada (eixo X) e a estimada (eixo Y) estabelecida pelo modelo proposto para o Parque Estadual do Jalapão. (R2: coeficiente de determinação; r: coeficiente de correlação de Pearson; d: coeficiente de Willmott; e c: índice de desempenho).

O valor de r = 0,84 indicou que a estimativa da produtividade pelo modelo proposto foi boa. A exatidão, representada pelo índice d = 0,91, demonstrou que o modelo foi capaz de estimar a produtividade de capim dourado com relativa exatidão. O desempenho do modelo, representado pelo valor de c = 0,77, é interpretado como muito boa, segundo o critério proposto por Camargo e Sentelhas (1997).

CONCLUSÕES

A produtividade de capim dourado do PEJ possui interação com os seguintes atributos do solo: pH(H₂O), Mg²⁺, Ca²⁺, areia fina, argila e umidade.

A umidade do solo foi a variável mais importante no modelo de regressão quadrática múltipla estabelecida para a produtividade do capim dourado no PEJ no estádio de maturação dessa espécie, ou seja, na época propícia para a coleta de escapos.

REFERÊNCIAS

• Araújo R, Goebert WJ, Lacerda MPC. Qualidade de um solo sob diferentes usos e sob cerrado nativo. R Bras Ci Solo. 2007;31:1099-108.
• Bayer C, Mielniczuk J, Martin-Neto L. Efeitos de sistemas de preparo e de cultura na dinâmica da matéria orgânica e na mitigação das emissões CO₂ R Bras Ci Solo. 2000;24:559-607.
• Caetano JO; Benites VM, Silva GP, Silva IR, Assis RL, C Filho A. Dinâmica da matéria orgânica de um Neossolo Quartzarênico de cerrado convertido para cultivo em sucessão de soja e milheto. R Bras Ci Solo. 2013;37:1245-55.
• Camargo AP, Sentelhas PC. Avaliação de desempenho de diferentes métodos de estimativa da evapotranspiração potencial no Estado de São Paulo, Brasil. R Bras Agrometeorol. 1997;5:89-97.
• Carvalho PGS. As veredas e sua importância no domínio dos cerrados. Inf Agropec. 1991;15:54-6.
• Empresa Brasileira de Pesquisa Agropecuária - Embrapa. Sistema brasileiro de classificação de solos. 3ª.ed. Brasília: Embrapa Solos; 2013.
• Empresa Brasileira de Pesquisa Agropecuária - Embrapa. Centro Nacional de Pesquisa de Solo. Manual de métodos de análise de solos. Rio de Janeiro: 1997.
• Giulietti AM, Wanderley MGL, Longhi-Wagner HM, Pirani JR, Parra LR. Estudos em “sempre vivas”: Taxonomia com ênfase nas espécies de Minas Gerais, Brasil. Acta Bot Bras. 1996;10:329-83.
• Guimarães AJ, Araújo GM, Correira GF. Estrutura fitossociológica em área natural e antropizada de uma vereda em Uberlândia, MG. Acta Bot Bras. 2002;16:317-29.
• Haridasan M. Nutritional adaptations of native plants of the cerrado biome and acid soils. Braz J Plant Physiol. 2008;20:183-95.
• Instituto Natureza do Tocantins. Naturatins. Levantamento aerofotogramétrico, na escala 1:25.000, do Parque Estadual do Jalapão. Projeto AASA 92/2002. Curitiba: Aeroimagem; 2002.
• Instituto Natureza do Tocantins. Naturatins. Portaria 362/2007 que adota medidas de ordenamento à coleta e ao manejo do capim dourado. Palmas: 2007.
• Lúcio ADC, Couto MRM, Lopes SJ, Storck L. Transformação Box-Cox em experimentos com pimentão em ambiente protegido. Hortic Bras. 2011;29:38-42.
• Lúcio ADC, Couto MRM, Trevisan JN, Martins GAK, Lopes SJ. Excesso de zeros nas variáveis observadas: Estudo de caso em experimento com brócolis. Bragantia. 2010;69:1035-46.
• Mayers RH. Classical and modern regression with applications. 2^nded. Boston: Duxbury Classic; 2000.
• Mesquita MAM, Naves RV, Souza EB, Bernardes TG, Silva LBE. Caracterização de ambientes com alta ocorrência natural de araticum (Annona crassiflora Mart.) no estado de Goiás. R Bras Frutic. 2007;29:15-9.
• Meurer EJ. Fundamentos de química do solo. 3ª.ed. Porto Alegre: Evangraf; 2006.
• Munhoz CBR, Felfile JM. Florística do estrato herbáceo-subarbustivo de um Campo Limpo Úmido em Brasília, Brasil. Biota Neotrop. 2007;7:205-15.
• Oliveira EMM, Ruiz HA, Alvarez V VH, Ferreira PA, Costa FO, Almeida ICC. Nutrient supply by mass flow and diffusion to maize plants in response to soli aggregate size and water potential. R Bras Ci Solo. 2010;34:317-27.
• Pacifico M, Napolitano A, Masullo M, Hilario F, Vilegas W, Piacente S, Santos LC. Metabolite fingerprint of “capim dourado” (Syngonanthus nitens), a basis of Brazilian handcrafts. Ind. Crops Produc. 2011;33:488-96.
• Ramos MVV, Curi N, Mota PEF, Vitorino ACT, Ferreira MN, Silva MLN. Veredas do Triângulo Mineiro: solos, água e uso. Ci Agrotec. 2006;30:283-93.
• Sano EE, Rosa R, Brito JLS, Ferreira LG. Mapeamento semidetalhado do uso da terra do bioma Cerrado (notas científicas). Pesq Agropec Bras. 2008;43:153-6.
• Sano SM, Simon MF. Produtividade do baru (Dipteryx alata Vong.) em ambientes modificados, durante 10 anos. In: Anais do 9º Simpósio Nacional do Cerrado, 2008; Brasília. Planaltina: Embrapa Cerrados; 2008.
• Santana HMP, Lacerda MPC, Barros MA, Barbosa IO. Unidades pedoambientais da região de Santa Tereza, estado do TO. Pesq Agropec Trop. 2010;40:8-19.
• Santana HMP, Sano EE, Lacerda MPC. Levantamento dos campos limpos úmidos no Parque Estadual do Jalapão, TO, por meio de fotografias aéreas analógicas: Abordagem metodológica. R Bras Cartogr. 2014;66:71-81.
• Santiago AR, Pereira Júnior AC. Mapeamento da cobertura da terra dos Parques Estaduais do Jalapão (PEC), Cantão (PEC) e município de Itaguatins (Tocantins). Ambiência. 2010;6:109-24.
• Schmidt IB, Figueiredo IB, Borghetti F, Scariot A. Produção e germinação de sementes de “capim dourado” Syngonanthus nitens(Bong.) Ruland (Euriocaulaceae): Implicações para o manejo. Acta Bot Bras. 2008;22:37-42.
• Schmidt IB, Figueiredo IB, Scariot A. Ethnobotany and efects of harvesting on the population ecology of Syngonanthus nitens(Bong.) Ruland (Euriocaulaceae), a NTFP from Jalapão region, central Brazil. Econ Bot. 2007;61:73-85.
• Schmidt IB, Mandle L, Ticktin T, Gaoue O. What do matrix population models reveal about the sustainability of non-timber forest product harvest? J Appl Ecol. 2011;48:815-26.
• Schmidt IB, Ticktin T. When lessons from population models and local ecological knowledge coincide - Effects of flower stalk harvesting in the Brazilian savanna. Biol Conserv. 2012;152:187-95.
• Secretaria de Planejamento e da Gestão Pública - Seplan. Atlas do Tocantins: Subsídios ao planejamento da gestão territorial. Palmas: 2012.
• Silva RB, Santos AC, Silva JEC, Oliveira OBT, Aráujo AS. Diagnóstico do solo em áreas de ocorrência de capim dourado na região do Jalapão, estado do Tocantins. Amazônia: Ci Desenv. 2010;6:91-101.
• Statistical Package for Social Sciences. SSPS - Release 19.0 for Windows. Chicago: 2008.
• Valladares GS, Gomes EG, Mello JCCBS, Pereira MG, Anjos LHC, Ebeling AG, Benites VM. Análise dos componentes principais e métodos multicritérios de estudo de Organossolos e solos afins. R Bras Ci Solo. 2008;32:285-96.
• Valladares SG. Caracterização de solos brasileiros com elevados teores de material orgânico. Magistra. 2008;20:95-104.
• Willmott CT. Some comments on the evaluation of model performance. Bull Am Meteorol Soc. 1982;63:1309-13.

Datas de Publicação

• Publicação nesta coleção
  Jul-Aug 2015

Histórico

• Recebido
  25 Jun 2014
• Aceito
  20 Mar 2015