Atributos físicos de Cambissolo Háplico em vinhedos submetidos a intensidades de tráfego

Physical attributes of a Haplic Cambisol in vineyards subjected to different traffic intensities

DIVISÃO 3 - USO E MANEJO DO SOLO

COMISSÃO 3.3 - MANEJO E CONSERVAÇÃO DO SOLO E DA ÁGUA

Atributos físicos de Cambissolo Háplico em vinhedos submetidos a intensidades de tráfego(¹)

Physical attributes of a Haplic Cambisol in vineyards subjected to different traffic intensities

Cassio Marques de Valois^I; Jucinei José Comin^II; Milton da Veiga^III; Gustavo Brunetto^IV; Jamil Abdalla Fayad^V; Fernando Cesar Bauer^VI; Arcângelo Loss^VI; Marcel Pires de Moraes^VII; Cleiton Junior Ribeiro Lazzari^VII

^IMestrando do Programa de Pós-Graduação em Agroecossistemas - PGA, Centro de Ciências Agrárias - CCA, Universidade Federal de Santa Catarina - UFSC. Campus Universitário Reitor João David Ferreira Lima Trindade, Bairro Itacorubi. CEP 88040-900 Florianópolis (SC). Bolsista da CAPES. E-mail: [e-mail]

^IIProfessor Associado VI, Departamento de Engenharia Rural (ENR), CCA, UFSC. Rodovia Admar Gonzaga, 1346, Bairro Itacorubi. CEP 888034-000 Florianópolis (SC). E-mail: [e-mail]

^IIIPesquisador, Estação Experimental da Empresa de Pesquisa Agropecuária e Extensão Rural de Santa Catarina (EPAGRI). BR 282, km 342. Caixa Postal 116. CEP 89620-000 Campos Novos (SC). E-mail: [e-mail]

^IVProfessor Adjunto II, Departamento de Solos, Centro de Ciências Rurais, Universidade Federal de Santa Maria. Av. Roraima, 1000, Bairro Camobi. CEP 97105-900 Santa Maria (RS). E-mail: [e-mail]

^VExtensionista, EPAGRI Regional da Grande Florianópolis. Caixa Postal 502. CEP 88.034-901 Florianópolis (SC). E-mail: [e-mail]

^VIProfessor Adjunto, ENR, CCA, UFSC. E-mail: [e-mail], [e-mail]

^VIIAcadêmico do Curso de Agronomia, ENR, CCA, UFSC. E-mail: [e-mail], [e-mail]

RESUMO

Em vinhedos, o tráfego de máquinas agrícolas causa compactação do solo que pode ser mais acentuada nas entrelinhas pela maior intensidade de tráfego. Este trabalho objetivou avaliar a variabilidade dos atributos físicos do solo em razão da posição de amostragem em relação à linha de plantio, em vinhedos comerciais da variedade Bordô em pé franco, submetidos a diferentes manejos da cultura. Em fevereiro de 2012, três vinhedos foram selecionados no município de Major Gercino (SC), sendo um manejado sem a utilização de máquinas agrícolas (vinhedo 1: com seis anos de implantação) e dois que utilizam máquinas nas entrelinhas (vinhedos 2 e 3: com oito e 17 anos, respectivamente). Na linha de plantio (Lp), na linha do rodado (Lr) e entre a linha de plantio e a linha do rodado (El), foram coletadas amostras indeformadas de solo nas camadas de 0,00-0,05, 0,05-0,10, 0,10-0,15 e 0,15-0,20 m de profundidade. Avaliaram-se a densidade do solo (Ds), a resistência à penetração (Rp), as classes de poros (macroporos - MaP, mesoporos - MeP, microporos - MiP e porosidade total - PT), os índices de agregação e de estabilidade de agregados (Diâmetro médio geométrico dos agregados secos ao ar e estáveis em água - DMGsa e DMGea, índice de estabilidade de agregados - IEA), a umidade gravimétrica (Ug), a densidade das partículas e a granulometria do solo. O estado de compactação aumentou no sentido da Lp para a Lr. A ocorrência de tráfego, em geral, influenciou negativamente os atributos avaliados. Entre as posições de amostragem, os maiores valores de Ds e Rp ocorreram na Lr, em que a Rp atingiu valores críticos ao desenvolvimento de raízes; nos vinhedos mecanizados, foram encontrados elevados valores de Rp e baixos de MaP nas camadas de 0,00-0,05 e 0,05-0,10 m. Observaram-se agregados com maior diâmetro e menores IEA, o que indica deterioração da qualidade física do solo. No vinhedo 1, houve menor variabilidade dos atributos físicos, enquanto no vinhedo 2 verificou-se influência negativa do tráfego de máquinas, principalmente para os atributos relacionados à agregação e elevada MiP. No vinhedo 3, observaram-se os maiores índices de Rp e de Ds, bem como o menor valor de MaP. Portanto, em relação aos vinhedos estudados, verificou-se aumento da degradação física do solo proporcionado pelo tráfego de máquinas (vinhedos 2 e 3). Dentre as posições de coletas e os vinhedos, evidenciou-se, por meio da análise de componentes principais, a separação da posição Lr do vinhedo 3, em razão dos atributos Ds e Rp, permitindo inferir que o manejo adotado na posição Lr está ocasionando maior degradação da qualidade física do solo.

Termos de indexação: variabilidade de atributos físicos do solo, Vitis labrusca, análise de componentes principais.

SUMMARY

Agricultural machine traffic causes soil compaction in vineyards, which may be accentuated between the rows due to greater intensity of traffic. The aim of this study was to evaluate the variability of soil physical properties as a function of the sampling position in relation to the plant row in commercial vineyards of the Bordô variety planted directly in the soil and subjected to different management practices. In February 2012, three vineyards were selected in the municipality of Major Gercino, southern Brazil - one without the use of agricultural machinery (vineyard 1 with six years of age) and two that use machines between the rows (vineyards 2 and 3, with eight and 17 years of age, respectively). Undisturbed soil samples were collected in the crop row (Cr), the wheel row (Wr), and between the crop row and the wheel row (Br) at the 0.00-0.05, 0.05-0.10, 0.10-0.15, and 0.15-0.20 m depths. We evaluated bulk density (Bd), resistance to penetration (Rp), pore classes (macropores - MaP, mesopores - MeP, micropores - MiP, and total porosity - TP), aggregation and aggregate stability indices (mean geometric diameter of air-dried aggregates - MGDad, and water stable aggregates - MGDws, and aggregate stability index - ASI), gravimetric moisture (Gm), and particle density and particle size distribution in the soil. Compaction increased in the direction from the Cr to the Wr. The occurrence of traffic, in general, negatively affected the properties evaluated. Among the sampling positions, the greatest values of Bd and Rp occurred in the Wr, where the Rp reached critical values for root development. In the mechanized vineyards, high values of Rp and low values of MaP were found in the 0.00-0.05 and 0.05-0.10 m layers. Aggregates of greater diameter and the lowest ASI were also found, which indicates deterioration of soil physical quality. In vineyard 1, there was less variability of the physical properties; while in vineyard 2 the negative effect of machine traffic was observed, especially for the properties related to aggregation and high MiP. In vineyard 3, greater indices of Rp and Bd were observed, as well as the lowest value of MaP. Therefore, in relation to the vineyards studied, an increase in soil physical degradation was caused by machine traffic (vineyards 2 and 3). Among the soil collection positions and the vineyards, separation of the Wr position of vineyard 3 was seen by means of principal component analysis, stemming from the Bd and Rp properties, from which it may be inferred that the management adopted in the Wr position of vineyard 3 is leading to greater degradation of soil physical quality.

Index terms: variability of soil physical properties, Vitis labrusca, principal component analysis.

INTRODUÇÃO

O Estado de Santa Catarina possui a quinta maior área cultivada com videira do Brasil, destacando-se a região do Vale do Rio Tijucas como uma das principais regiões vitivinícolas. Cultivam-se videiras americanas (Vitis labrusca), e as uvas são destinadas à elaboração de sucos e vinhos de mesa. Os vinhedos são instalados, predominantemente, em Cambissolos localizados em relevo plano a suave ondulado.

Anualmente, por causa da elevada e constante umidade relativa do ar e altas temperaturas, especialmente durante o período vegetativo e reprodutivo, as videiras são submetidas a contínuas aplicações de fungicidas para o controle preventivo de doenças fúngicas foliares. As aplicações são realizadas usando-se máquinas e equipamentos agrícolas, como tratores e pulverizadores, que trafegam nas entrelinhas dos vinhedos, inclusive quando o solo possui teor de água acima da capacidade de campo. Isso causa a aproximação das partículas do solo pela pressão dos rodados e, consequentemente, pode provocar a compactação do solo nas entrelinhas (van Dijck & van Asch, 2002; Lima et al., 2004; Rosa et al., 2009; Braida et al., 2011).

A compactação do solo em entrelinhas de vinhedos (van Dijck & van Asch, 2002; Ferrero et al., 2005), citrus (Lima et al., 2004) e pereira (Vocanson et al., 2006) pode modificar os atributos físicos do solo, como a densidade do solo, a resistência à penetração, a porosidade, a estabilidade de agregados, a estrutura e o grau de floculação das argilas (Silva et al., 2006; Echenique et al., 2007; Portugal et al., 2008; Rosa et al., 2009), comparativamente aos atributos físicos na linha de plantio dos pomares onde não é realizado o tráfego. Entretanto, a modificação dos atributos depende da intensidade do tráfego a que o solo é submetido (Lima et al., 2004; Elaoud & Chehaibi, 2011). Essas alterações reduzem o suprimento de oxigênio e promovem restrições físicas ao crescimento das raízes das plantas, diminuindo o volume de solo explorado e, por consequência, podem reduzir a quantidade de água e nutrientes absorvidos, refletindo negativamente na taxa fotossintética, no crescimento e no rendimento das culturas (De Maria et al., 1999; Beutler & Centurion, 2004; Beutler et al., 2006; Fidalski et al., 2010; Costa et al., 2012; Luciano, 2012).

Valores de densidade (Ds) em solos arenosos considerados normais estão compreendidos entre 1,20 e 1,90 kg dm^-3; o valor médio para solos compactados que redunda em restrição ao crescimento radicular (Ds crítica) encontra-se na faixa de 1,65 kg dm^-3 (Reinert & Reichert, 2006). Já o valor proposto por Reichert et al. (2003) de Ds crítica para solos com textura média (20 a 55 % de argila) foi de 1,55 kg dm^-3. Como a porosidade total é inversamente proporcional à densidade do solo, e a macroporosidade é a maior responsável pela aeração do solo e drenagem interna do perfil, é necessário que essa esteja numa faixa mínima de 10 % para suprir a demanda de oxigênio no solo (Xu et al., 1992).

Para diversas culturas, os valores de resistência à penetração (Rp) acima dos 2,0 MPa são considerados limitantes para o crescimento das raízes (Taylor et al., 1966; Nesmith, 1987). Para a cultura da videira, em solos com horizontes muito coesos até 0,35 m de profundidade, onde os valores de Rp ultrapassam os 3,0 MPa, Echenique et al. (2007) verificaram redução no rendimento, na área foliar e na massa de material podado da videira. Em geral, o aumento da Rp tem sido verificado em solos cultivados em razão do tráfego de máquinas agrícolas, onde há tendência de aumento da Rp com o aumento da intensidade do tráfego (Elaoud & Chehaibi, 2011). Além disso, a Rp também é utilizada para comparar diferentes manejos e locais em sistemas de produção, evidenciando a degradação física de solos submetidos à aplicação de pressão na superfície pelo tráfego de máquinas, como o caso de compactação de pomares (Echenique et al., 2007; Portugal et al., 2008; 2010), reflorestamento (Seixas & Souza, 2007) e culturas anuais (Fuentes Llanillo et al., 2006; Andrade et al., 2010), sendo também associada a um menor crescimento radicular (Costa et al., 2012).

No Brasil, a maioria dos estudos que abordam a compactação e problemas estruturais de solo trata de culturas anuais e solos que são revolvidos anualmente, enquanto poucos trabalhos relatam situações de solos não revolvidos sob culturas perenes (Reichert et al., 2007), como é o caso dos vinhedos. Esses autores consideram que, para uma avaliação confiável e que possibilite a inferência sobre a qualidade física do solo, é necessário determinar vários atributos físicos e interpretá-los em conjunto para evitar uma avaliação errônea do estado de compactação de determinado solo.

Considerando a constante utilização de tráfego de máquinas agrícolas nas entrelinhas dos vinhedos instalados em Cambissolos, que é uma das classes de maior ocorrência em SC (Embrapa, 1998), este trabalho objetivou avaliar a variabilidade dos atributos físicos do solo em razão da posição de amostragem em relação à linha de plantio, em três vinhedos comerciais da variedade Bordô em pé franco, submetidos a diferentes manejos.

MATERIAL E MÉTODOS

Este trabalho foi realizado em 2012 em três vinhedos comerciais, no município de Major Gercino, situado no Alto Vale do Rio Tijucas, SC, nas coordenadas geográficas 27^º 26' 42" sul e 49^º 06' 27" oeste (Vinhedo 1 - 1,5 ha), 27^º 25' 16" sul e 49^º 04' 31" oeste (Vinhedo 2 - 6 ha), 27^º 25' 18" sul e 49^º 05' 53" oeste (Vinhedo 3 - 3,4 ha), com altitudes de 530 a 545 m e distância entre os vinhedos de 1,5 e 2,5 km. O clima do local é temperado úmido, com verões quentes e sem estação seca definida, do tipo Cfa, de acordo com a classificação de Köppen, com temperatura média anual de 20 ^ºC e precipitação pluvial entre 1.500 e 1.700 mm bem distribuída ao longo do ano (Pandolfo et al., 2002). Os solos dos vinhedos foram classificados como Cambissolo Háplico Tb eutrófico típico, textura média, pouco cascalhento a cascalhento, A moderado, relevo plano, substrato sedimentos recentes de Granito da Suíte Valsungana (Embrapa, 2006). Os atributos químicos determinados em amostras de solo coletadas em cada vinhedo na camada de 0,00-0,20 m são apresentados no quadro 1.

Os vinhedos foram compostos pela variedade Bordô (Vitis labrusca) em pé franco, implantados nos anos de 2006, 2004 e 1995 (Vinhedos 1, 2 e 3, respectivamente). Anteriormente à implantação desses, as áreas eram cultivadas com fumo (Nicotiana tabacum) em sistema de manejo com revolvimento intensivo do solo (aração e gradagem) e uso constante de agrotóxicos e adubações com fertilizantes de síntese química, o que explica a característica de epi-eutrofismo desses solos (Quadro 1), que originalmente eram classificados como álicos (Embrapa, 1998).

O sistema de condução utilizado nos três vinhedos é do tipo latada, com espaçamento de 3 m entre fileiras e 2 m entre plantas. Em todos os vinhedos, foram semeadas a lanço, sem incorporação, aveia-preta (Avena strigosa) e ervilhaca (Vicia sativa), sendo roçadas quatro vezes por ano; e a matéria seca foi depositada na superfície do solo.

No vinhedo 1, considerada área-controle, não foi utilizado maquinário agrícola e todas as práticas de manejo foram realizadas manualmente. O tráfego de máquinas ocorreu externamente ao vinhedo. No vinhedo 2, o tráfego foi realizado 20 vezes por ano, com um trator marca Tramontini^® 4 × 2 modelo T3025-4 (950 kg de massa total, rodado dianteiro 6'' × 12'' e traseiro 9.5'' × 24'', com pressão de inflação de 25 e 40 lb, respectivamente), com taxa de compressão estimada de 0,387 kg cm^-2 (Mialhe, 1980). Já no vinhedo 3, o tráfego foi realizado 30 vezes por ano com um trator de rabiças marca Yanmar^® modelo Cultivador Motorizado TC14 Super (433 kg de massa total, rodado 6'' × 12'' com pressão de inflação de 25 lb), com taxa de compressão estimada de 0,928 kg cm^-2. Utilizou-se um pulverizador de turbina com capacidade para 300 L, nos vinhedos 2 e 3.

Para determinação dos atributos físicos do solo, foram coletadas amostras com estrutura preservada utilizando-se anéis volumétricos de aço inox com 0,05 m de altura e 0,07 m de diâmetro interno. As coletas foram realizadas em quatro camadas (0,00-0,05, 0,05-0,10, 0,10-0,15 e 0,15-0,20 m de profundidade) e em três locais distintos nos vinhedos: na linha de plantio (Lp); na linha do rodado (Lr); e entre a linha do rodado e a linha de plantio (El), com três repetições por local de coleta. As amostras foram coletadas em três subparcelas homogêneas dos vinhedos com a seleção do local de coleta feita aleatoriamente com o uso de um marcador lançado em cada uma delas. Os pontos de coleta se distanciaram aproximadamente 40 m entre si.

Depois de realizada a limpeza e preparação das amostras, utilizaram-se os métodos de rotina para a determinação dos atributos físicos do solo descritos em Veiga (2011). Para calcular o volume de água retida nas classes de poros e a densidade do solo (Ds), foi determinada a massa das amostras após sua saturação por capilaridade e, sequencialmente, após aplicadas tensões de 6 kPa (box de areia) e 600 kPa (Extrator de Richards), bem como após efetuada a secagem em estufa a 105 ^ºC, nesse caso, de uma fração da amostra até massa constante. A Ds (kg dm^-3) corresponde à relação entre a massa de solo seco em estufa e o volume interno do respectivo anel volumétrico. O volume total de poros (PT), macro (MaP, poros com Ø >50 µm), meso (MeP, poros com Ø entre 50-0,5 μm) e microporos (MiP, poros com Ø < 0,5 μm) correspondem, respectivamente, à relação entre o volume de água retido no solo após sua saturação, entre a saturação e a tensão de 6 kPa, entre as tensões de 6 e 600 kPa e com umidade equilibrada na tensão de 600 kPa, todos calculados em relação ao volume interno do respectivo anel volumétrico.

A resistência do solo à penetração foi determinada em laboratório no centro das mesmas amostras utilizadas para discriminar a porosidade e Ds, após estabilização da umidade na tensão de 600 kPa. Para essa determinação, foi utilizado um penetrômetro de bancada marca Marconi, modelo MA933, equipado com pino de penetração de 80 mm de altura e ponta cônica angular de 30º, ajustado para velocidade de penetração de 1,0 mm s^-1 e registros em intervalos de 1,0 mm. O valor para cada amostra correspondeu à média dos 30 valores obtidos entre 11 e 40 mm de profundidade. Feita a leitura no penetrômetro, foi pesado aproximadamente 20 g do solo para determinar o fator de correção da umidade gravimétrica (Ug) na tensão de 600 kPa para solo seco.

Após definida a Rp, iniciou-se a análise de distribuição de tamanhos e estabilidades dos agregados, em que foram analisados os agregados secos ao ar e estáveis em água. O solo utilizado foi retirado dos anéis e cuidadosamente destorroado, separando-se os agregados nos planos de fraqueza, até que toda a amostra passasse pela peneira com abertura de malha de 8 mm. As amostras foram secas ao ar livre e à sombra por uma semana.

Para se determinar a distribuição de tamanho dos agregados secos ao ar, foram pesados de 25 a 30 g da amostra e, em seguida, colocadas na sequência de peneiras com 20 cm de diâmetro e 5 cm de altura, com aberturas de malha de 4,00; 2,00; 1,00; 0,50 mm; e fundo, sequencialmente. Com a amostra nas peneiras, foram realizados 12 movimentos de ida e vinda em um sentido, girado o conjunto em 90º e repetidos os movimentos, com a posterior pesagem do material retido em cada peneira e no fundo, calculando-se o diâmetro médio geométrico de agregados secos ao ar (DMGsa).

Essa amostra utilizada foi recomposta, com a exclusão do material retido no fundo para determinar a estabilidade dos agregados em água. Para essa definição foi utilizada a mesma sequência de peneiras acopladas a um oscilador mecânico, com amplitude de oscilação de aproximadamente 4 cm e frequência de oscilação de 30 ciclos min^-1 dentro de recipiente com água, similar ao aparelho de Yoder. As amostras foram acondicionadas no fundo da primeira peneira e ficaram em descanso por 10 min (umedecendo). Decorrido esse tempo, foi iniciada, por mais 10 min, a agitação das peneiras para separação do material e, em seguida, o solo retido nas diferentes peneiras foi seco em estufa a 105 °C, calculando-se o diâmetro médio geométrico de agregados estáveis em água (DMGea). De posse desses valores, foi realizado o cálculo do índice de estabilidade de agregados (IEA) pela fórmula IEA = DMGea/DMGsa.

A densidade das partículas (Dp) foi determinada utilizando-se o método do balão volumétrico e álcool, sendo os cálculos em base de massa (Gubiani et al., 2006). A granulometria foi discriminada pelo método da pipeta, adaptado por Veiga (2011) para análise de rotina.

Os dados provenientes das análises físicas do solo foram avaliados quanto à homocedastia, pelo teste de Cochran (Snedecor & Cochran, 1989) e distribuição normal dos resíduos, pelo teste de Lilliefors (Campos, 1979). Posteriormente, foram submetidos à análise de variância com aplicação do teste F e dos valores médios, quando significativos, comparados entre si pelo teste de Tukey a 5 %. Não foi necessária nenhuma transformação de dados. Cada vinhedo foi considerado uma área, sendo os tratamentos considerados os locais de coleta (Lp, El, Lr), com três repetições por profundidade. Também foi realizada uma análise de componentes principais (ACP), por meio do programa XLSTAT (Addinsoft^®, 2010), para cada camada avaliada e entre os locais de coleta, cujas variáveis selecionadas foram as que apresentaram correlações de Pearson maiores que 0,6 entre si. A ACP é utilizada para reduzir as dimensões dos dados e, consequentemente, facilitar a análise por meio do gráfico do círculo de correlações (Herlihy & McCarthy, 2006).

RESULTADOS E DISCUSSÃO

Atributos físicos do vinhedo 1

No vinhedo 1, onde não há tráfego de máquinas, quando o solo foi influenciado, decorreu do pisoteio (humano) na região compreendida pelas posições El e Lr. Para MaP e PT, verificaram-se diferenças, na profundidade de 0,05-0,10 m, entre a Lp e a Lr, com maiores valores para Lp, intermediários para a El e menores para Lr (Quadro 2). Com maior valor de MaP na Lp, a PT também apresentou o mesmo padrão, pois a MaP é a primeira que sofre interferência no processo de compactação (Reichert et al., 2007). Isso também pode ter influenciado a Ug, que foi superior na Lp em relação à Lr. Aparentemente, a reduzida pressão causada pelo pisoteio foi suficiente para alterar a MaP. Em estudo conduzido em um Cambissolo Háplico em Bento Gonçalves, RS, também se observou redução da MaP na entrelinha em relação à linha de plantio de um vinhedo, mas essa diferenciação ocorreu em razão do tráfego de máquinas agrícolas (Rosa et al., 2009). Tanto a MeP como a MiP não foram influenciadas pelo manejo do vinhedo 1 (Quadro 2).

A Ds foi superior na El em relação à Lp na camada de 0,10-0,15 m, ambas não diferindo da Lr, podendo indicar o efeito do manejo para as podas, colheitas e pulverizações com aplicação de forças no solo mais próximas às plantas de videira. No entanto, esses valores encontrados são menores do que o limite para considerar o solo compactado (Reinert & Reichert, 2006).

Adicionalmente, os valores de MaP são adequados e suficientes para permitir a drenagem do solo em todas as posições, sempre próximos dos 0,200 m³ m^-3 (Xu et al., 1992; van Lier, 2001). Os atributos relacionados à agregação e estabilidade dos agregados (DMGsa, DMGea e IEA) também foram semelhantes, indicando que, em relação à estrutura do solo, o manejo não causou efeitos nos diferentes pontos de amostragem, submetidos ou não à ação antrópica por meio do pisoteio.

No vinhedo 1, não foram verificados efeitos das posições de amostragem sobre a Rp, mas todos os valores determinados foram 2,5 MPa. Nas profundidades de 0,05-0,10 e 0,10-0,15 m, verificaram-se valores superiores a 3,0 MPa para a El e Lp, respectivamente (Quadro 2), indicando existência de condições restritivas ao desenvolvimento radicular para diversas culturas, inclusive para a videira (Taylor et al., 1966; Echenique et al., 2007). O histórico de uso do solo, anteriormente à implantação do vinhedo, associado ao intenso revolvimento do solo, refletiu diretamente no atributo físico Rp.

Dessa forma, em razão da ausência de tráfego de máquinas e por ser um vinhedo novo, para a maioria dos atributos físicos avaliados, não foram detectadas diferenças entre posições de amostragem, indicando pequeno efeito do pisoteio humano sobre os atributos físicos do solo.

Atributos físicos do vinhedo 2

No vinhedo 2, que sofre tráfego de máquinas agrícolas, observaram-se diferenças entre as posições de coleta para a Rp na camada de 0,00-0,05 m (mesmo sem diferenças para Ug), onde verificou-se o maior valor na Lr (4,2 MPa), seguido da El (2,9 MPa) e com o menor valor na Lp (1,7 MPa) (Quadro 2). Esse atributo evidenciou uma variação do estado físico do solo causado pela pressão dos rodados no local de tráfego, onde a Rp atingiu valor considerado restritivo ao crescimento radicular (Echenique et al., 2007). Além disso, na Lr, observou-se aumento de 11 % para a Ds e redução de 38 % para a MaP, em comparação à Lp (Quadro 2). Na posição El, a Rp apresentou valor limitante, mas não restritivo, e na Lp não verificou-se resistência considerada crítica ao desenvolvimento radicular para nenhuma cultura (Taylor et al., 1966; Nesmith, 1987).

Na camada de 0,10-0,15 m, a Ds apresentou diferença entre a Lp e a El, com valores intermediários na Lr. O esperado seria essa diferença ocorrer entre Lp e Lr, mas o maior valor foi verificado na El, que pode em certas ocasiões também ter sofrido tráfego, bem como influência desse, em razão da proximidade com a Lr, que se apresenta mais compactada em relação à Lp.

Para a MeP, somente verificaram-se maiores valores na Lp em relação às outras posições na camada de 0,15-0,20 m. Para as demais camadas, não houve diferenças significativas (Quadro 2). Isso pode ser por causa da maior destruição da MeP ocasionada pela pressão das máquinas na Lr e El, transformando-a em poros menores, o que é corroborado pela correlação alta e negativa entre a MeP dessa camada com a Rp da de 0,00-0,05 m (r = -0,75, p<0,05) e a Ds da de 0,10-0,15 m (r = -0,73, p<0,05). Também foram verificadas correlações negativas entre a MeP com a Ds para a camada de 0,00-0,05 m (r = -0,91, p<0,05), entre a MeP da de 0,10-0,15 m e a Rp da de 0,0-0,05 (r = -0,69, p<0,05); e correlações positivas entre a Ds da camada de 0,10-0,15 m e a Rp das camadas de 0,00-0,05 m (r = 0,75, p<0,05) e 0,10-0,15 m (r = 0,67, p<0,05). Esse padrão demonstra o efeito negativo do tráfego de máquinas, acarretando em aumento dos valores de Ds e Rp, com consequências negativas para a MeP.

A MaP no vinhedo 2 apresentou valores inferiores (proporcionalmente) àqueles encontrados no vinhedo 1, o que pode estar atrelado à diferença granulométrica entre os dois (Quadro 1), sendo verificado, na camada de 0,00-0,20 m, valores de correlação de r = 0,60; r = -0,45; e r = -0,52, respectivamente, para areia, silte e argila (p<0,05) e a MaP. A redução da MaP também é por causa da compactação do solo causada pelo tráfego no vinhedo 2, o que é corroborado pelos valores de correlações (altos e negativos) entre a MaP (camada de 0,00-0,05 m) e a Ds das camadas de 0,10-0,15 m (r = -0,86, p<0,05) e 0,15-0,20 m (r = -0,72, p<0,05), e a Rp da camada de 0,05-0,10 m (r = -0,72, p<0,05). Para a camada de 0,00-0,05 m, é possível aferir redução de cerca de 30 % na MaP da Lr, em relação à Lp (Quadro 2), ao mesmo tempo em que essa chega em valor próximo àquele tido como base de volume de poros que se apresentem livres de água para permitir o fluxo de gases no solo (van Lier, 2001; Klein et al., 2008).

A baixa permeabilidade à água e ao ar em solos com predominância da MiP, como no caso do vinhedo 2, que apresenta em todos locais e profundidades mais de 60 % da PT composta pela MiP e chegando a 75 % na camada 0,00-0,05 m na Lr , determina elevada capacidade de retenção de água e prejudica o sistema radicular, em razão da escassez de oxigênio às raízes (Dalbó et al., 2011). E, ainda, quando presente em regiões com elevada temperatura e com elevados índices de pluviosidade, aumenta a propensão à ocorrência de doenças fúngicas, o que pode prejudicar ainda mais a produtividade da cultura (Protas, 2003).

Em relação ao atributo DMGsa, é possível verificar valores proporcionalmente maiores, se comparado com os encontrados no vinhedo 1. No entanto, esses maiores agregados secos ao ar, que apresentaram r = 0,63 (p<0,05) com os teores de argila (camada de 0,00-0,20 m), não se mantiveram estáveis em água, sendo verificado correlação negativa (r = -0,39, p<0,05) entre DMGea e os teores de argila. Isso pode estar relacionado ao efeito do tráfego sobre a estrutura do solo, onde os agregados naturais são destruídos e ocorre aumento daqueles formados por efeito de compactação, os quais são menos estáveis, fato evidenciado pelo menor IEA (r = -0,60, p<0,05, com os teores de argila na camada de 0,00-0,20 m), que não chegou a atingir 0,50 no vinhedo 2 (Quadro 2) (Carpenedo & Mielniczuk, 1990; Albuquerque et al., 2005; Salton et al., 2008), mesmo com maior teor de argila e, portanto, mais agente cimentante que o vinhedo 1.

Portanto, no vinhedo 2, tem-se a formação de agregados maiores em razão da destruição de agregados naturais pela compactação causada pelo tráfego, que forma agregados de maior tamanho (Bronick & Lal, 2005). Esse comportamento pode ser explicado pela aproximação das partículas causada pela pressão das máquinas agrícolas, resultando em agregados maiores (maiores índices de DMGsa), porém sem o processo de estabilização desses, já que apenas a entrada de água é suficiente para desestabilizá-los (DMGea consideravelmente menor no vinhedo 2, proporcionalmente aos valores verificados para o vinhedo 1). Isso porque não se formaram pela ação dos agentes cimentantes permanentes e temporários, que estão ligados aos processos biológicos de formação e estabilização desses agregados (Tisdall & Oades, 1982; Carpenedo & Mielniczuk, 1990; Salton et al., 2008).

A correlação alta e negativa entre Ds com DMGea, sendo r = -0,69 (p<0,05) e r = -0,76 (p<0,05), respectivamente, para as camadas de 0,05-0,10 e 0,10-0,15 m, corrobora o efeito negativo do tráfego de máquinas, acarretando aumento dos valores de Ds e diminuição dos de DMGea. Também foram observados valores de correlação negativos entre IEA com a Ds (r = -0,60, p<0,05) e a Rp (r = -0,36, p<0,05), ambos para a camada de 0,00-0,20 m.

Atributos físicos do vinhedo 3

Para o vinhedo 3, verificou-se efeito significativo das posições de amostragem sobre os atributos Ds e Rp nas camadas até 0,10 m, com valores variando entre 0,70 e 1,49 kg dm^-3 para a Ds e entre 1,51 e 4,50 MPa para a Rp (Quadro 2). Os maiores valores de Ds foram encontrados na posição de coleta Lr, mesmo local onde também foram observados os maiores valores para Rp, demonstrando alta correlação entre esses dois atributos (r = 0,75, r = 0,67, p<0,05, para as camadas superficiais do solo), consequência da maior intensidade do tráfego e compressão do solo pelas máquinas associada ao maior tempo de implantação desse vinhedo. Também verificou-se correlação negativa da PT com a Ds na camada superficial do solo (r = -0,79, p<0,05) e não foram verificadas diferenças estatísticas entre as posições de coleta para PT. Mesmo sem apresentar diferenças na PT (0,00-0,05 m), observou-se redução de 30 % na Lr em relação à Lp, o que corrobora com os maiores valores de Ds e Rp na Lr.

Na camada superficial, a Ds foi maior na Lr e menor na El, enquanto na camada de 0,05-0,10 m a El não diferiu das demais, onde essa posição é intermediária à Lp (que sofreu menor alteração dos atributos físicos) e à Lr (mais alterada). Esse é o padrão observado nos vinhedos que sofrem tráfego de máquinas agrícolas. A Rp, que possui alta e positiva correlação com a Ds, como mencionado anteriormente, evidenciou a variação ocorrida no solo na posição Lr, apresentando valor superior a 100 % em relação às posições El e Lp, que não diferiram entre si. Ao estudar um Cambissolo Húmico com 443 g kg^-1 de argila sob diferentes usos em Lages, SC, com e sem preparo do solo, Andrade et al. (2010) observaram valores de Ds entre 1,13 e 1,38 kg dm^-3, com os maiores valores a uma profundidade próxima dos 0,15 m. A Rp nesse solo acompanhou o aumento da Ds, com valores entre 0,13 e 1,54 MPa, confirmando a boa relação entre esses atributos para inferir sobre as condições para o crescimento radicular. Também em um estudo em um Cambissolo Háplico em Bento Gonçalves, RS, observou-se na entrelinha em relação à linha de plantio de um vinhedo, alteração em atributos físicos, como aumento da Ds na entrelinha, decorrente do tráfego de máquinas agrícolas (Rosa et al., 2009).

De acordo com o valor de Ds crítica proposta por Reichert et al. (2003) de 1,55 kg dm^-3 para solos com textura média (20 a 55 % de argila), os valores verificados na Lr do vinhedo 3, bem como na Lr do vinhedo 2, ambos para a camada de 0,05-0,10 m, causam problemas para o crescimento radicular, o que também foi verificado por Suzuki (2005), em um Argissolo Vermelho distrófico latossólico (278 g kg^-1 de argila, 331 g kg^-1 de silte e 391 g kg^-1 de areia), em Alto das Canas, RS, causando redução no crescimento radicular da soja.

Juntamente a esses valores de Ds, é importante salientar que os valores de Rp, observados nos solos do vinhedo 2 e 3 (Quadro 2), foram superiores a 3,0 MPa na posição Lr, sendo considerados limitantes ou mesmo restritivos ao crescimento radicular da videira e das plantas de cobertura (Taylor et al., 1966; Nesmith, 1987; Merotto Jr. & Mundstock, 1999; Echenique et al., 2007). Essa limitação ainda pode ser maior ou menor dependendo da umidade do solo, como exposto por Genro Junior et al. (2004), os quais verificaram diferentes valores de Rp ao realizarem amostragens em um Latossolo em Ijuí, RS, que variaram inversamente com a umidade volumétrica do solo. Neste trabalho, foram encontradas correlações negativas entre a Rp da camada de 0,00-0,05 m e a Ug da de 0,05-0,10 m (r = -0,69, p<0,05) e entre a Rp da de 0,05-0,10 m e a Ug da de 0,15-0,20 m (r = -0,69, p<0,05).

O vinhedo 3, na camada de 0,05-0,10 m, mesmo na posição Lp, apresentou elevados valores de Rp, o que pode interferir negativamente no crescimento radicular e, consequentemente, na produção das videiras (Echerique et al., 2007). Nesse vinhedo, apesar de o solo ter sofrido maior intensidade de tráfego e maior compressão do solo, foram observados valores de DMGsa, DMGea e IEA intermediários, em relação aos vinhedos 1 e 2, o que pode estar relacionado ao maior teor de matéria orgânica do solo (MOS) nessa área (r = 0,78, p<0,05, com DMGea), comparativamente ao solo dos demais vinhedos (Quadro 1), o que pode aumentar o tamanho dos agregados estáveis em água. Esses maiores teores de MOS são decorrentes do uso de espécies de plantas de cobertura do solo (aveia e ervilhaca), que, em razão do maior tempo de implantação do vinhedo 3, pode ter favorecido o aumento dos teores de MOS. Os maiores teores de MOS contribuem na agregação, na consistência e na retenção de umidade dos solos, além de aumentar o número de ligações entre as partículas, proporcionando melhores condições à formação e estabilidade de macroagregados (Tisdall & Oades, 1982; Salton et al., 2008; Santos, 2008; Braida et al., 2011). Os maiores teores de MOS associados ao sistema radicular das espécies de plantas de cobertura do solo podem ter auxiliado na redução da ação negativa do tráfego sobre os agregados, em relação ao vinhedo 2, fato evidenciado pelo IEA, que apresentou valores mais similares ao vinhedo 1.

Análise dos componentes principais (ACP)

A distribuição das variáveis analisadas na ACP é apresentada na figura 1, com variância acumulada de 67,99 % (eixos F1 e F2), sendo, desse total, maior efeito expresso pelo eixo F1 (43,37 %). Verificou-se que o vinhedo 1 separou-se dos demais, formando um agrupamento das posições de coleta, principalmente por causa dos atributos MaP, IEA e DMGea, que ficaram todos em um mesmo quadrante, sendo a maior ênfase desses atributos dada à posição de Lp, em comparação às posições Lr e El, que apresentaram maior proximidade do atributo Ds (0,00-0,05 m). Esse padrão indicou que nas posições de Lr e El, por serem aquelas que sofrem maior pisoteio humano para os tratos culturais do vinhedo, apresentaram maiores valores de Ds, corroborando com os dados apresentados no quadro 2. Para os atributos relacionadas à posição Lp (MaP, IEA e DMGea), verificou-se que esses remetem às boas e desejadas características físicas do solo, por evidenciarem condições favoráveis à drenagem e representarem a melhor agregação do solo do vinhedo 1, em relação aos vinhedos 2 e 3.

Na posição Lp do vinhedo 1, a maior influência foi atribuída aos atributos IEA e MaP, portanto demonstrando nesse local, e mesmo no vinhedo, características de um solo com boa qualidade física e, por meio dessa, diferenciando o manejo realizado nesse vinhedo em relação aos demais vinhedos.

Próximo ao grupamento das posições do vinhedo 1, tem-se a posição Lp do vinhedo 3. Essa aproximação é decorrente do atributo DMGea ao longo das camadas analisadas, o que pode estar relacionado ao maior teor de MOS, comparativamente ao solo dos demais vinhedos (Quadro 1), mantendo o tamanho dos agregados estáveis em água, com valores similares entre o vinhedo 1 e o 3, conforme verificado no quadro 2. No entanto, no vinhedo 3 a posição Lr se destacou das demais posições do próprio vinhedo, assim como dos vinhedos 1 e 2, em razão dos atributos Rp (camadas de 0,00-0,10 m) e Ds (camadas de 0,10-0,20 m). Esse vinhedo é o que tem a maior intensidade de tráfego e há mais tempo sujeito a ele, com maiores valores de Ds e Rp, considerados críticos ao desenvolvimento da videira (Echenique et al., 2007).

A mesma posição (Lr) no vinhedo 2 também se aproximou de atributos que remetem à degradação física do solo, principalmente a Rp (0,15-0,20 m), a MeP (0,10-0,15 m), a MiP (0,10-0,20 m) e o DMGsa (0,00-0,20 m), e se distanciou de atributos como o DMGea e o IEA, o que indicou que os agregados maiores (maior DMGsa) foram formados pela compactação causada pelo tráfego às custas da destruição dos agregados formados por processos naturais (Bronick & Lal, 2005). Porém, esses agregados maiores do vinhedo 2 não são resistentes às forças de desagregação pelo rápido umedecimento e peneiramento em água, fato evidenciado pela maior distância do valor de DMGea e do IEA da posição Lr, e, consequentemente, maior aproximação da posição Lp do vinhedo 1. Esse padrão corroborou com os índices DMGea e IEA dos vinhedos 1 e 2 (Quadro 2). O tráfego também confere uma elevação do valor da MiP e redução do valor da MaP em razão da pressão exercida pelos rodados na superfície, pois essa é a primeira classe de poros influenciada no processo de compactação (Reichert et al., 2007).

As posições Lr e El do vinhedo 2 e El do vinhedo 3 também ficaram mais próximas do atributo Rp (0,10-0,15 e 0,15-0,20 m) do que a Lr do vinhedo 3, bem como estão mais próximas dos atributos PT, MeP, MiP e DMGsa, que por sua vez, excluindo a MeP, inferem em uma qualidade física mais depreciada nesses dois vinhedos, em comparação ao vinhedo 1.

A posição Lp do vinhedo 2 se diferenciou das demais, principalmente em relação à PT, MeP e Ug. No entanto, bem como o caso da posição Lp do vinhedo 3, isso está mais condicionado a uma característica inerente ao solo do vinhedo do que causada pelo manejo, pois o solo desse vinhedo apresentou valor mais elevado de silte e argila (Quadro 1), portanto causando diferenciação granulométrica em relação aos vinhedos 1 e 3.

As máquinas utilizadas nos vinhedos 2 e 3 aparentam não ter causado expressivas diferenças entre si quanto à deterioração do solo, mas há possibilidade de maior restrição no vinhedo 3, por causa do maior tempo de uso e da maior intensidade de tráfego e carga aplicada, como também observado por Elaoud & Chehaibi (2011), que mensuraram maior compactação do solo em razão da maior intensidade de tráfego. Nos vinhedos 2 e 3, principalmente em relação à posição Lr, e mesmo na El, verificou-se que intensivo uso de máquinas agrícolas causou degradação física do solo.

Dentre as posições de coletas e os vinhedos estudados, evidenciou-se, por meio da ACP, a separação da posição Lr do vinhedo 3 decorrente dos atributos Ds e Rp. Isso se deve ao maior tempo de uso, maior intensidade de tráfego e compressão aplicada ao solo, como também observado por Elaoud & Chehaibi (2011), que mensuraram maior compactação do solo em razão da maior intensidade de tráfego. Assim, pode-se inferir que o manejo adotado no vinhedo 3 na posição Lr está ocasionando maior degradação da qualidade física do solo, quando comparado aos demais vinhedos.

CONCLUSÕES

1. O maior estado de compactação do solo ocorreu no sentido da linha de plantio para a do rodado, com maiores alterações nos atributos relacionados à porosidade do solo nos dois vinhedos submetidos ao tráfego de máquinas, à agregação e estabilidade de agregados no vinhedo 2 e à densidade do solo e da resistência à penetração no vinhedo 3.

2. No vinhedo 1, por meio da ACP, evidenciaram-se as melhores condições físicas do solo; entretanto, nos vinhedos 2 e 3, encontrou-se a degradação física. A posição na linha do rodado do vinhedo 3 diferenciou-se das demais posições e dos vinhedos, caracterizando maior degradação da qualidade física do solo.

AGRADECIMENTOS

À UFSC e EPAGRI, pela disponibilização de infraestrutura, materiais e apoio na pesquisa. Ao MDA/SAF/CNPq, edital Nº 58/2010, pelo auxílio financeiro. À CAPES, pela bolsa de mestrado concedida ao primeiro autor.

LITERATURA CITADA

Recebido para publicação em 25 de junho de 2013

Aprovado em 19 de maio de 2014

• ADDINSOFT. XLSTAT statistical analysis software. Version 2013.2. Disponível em: <www.xlstat.com>. Acesso em: 10 de Jul. 2010.
• ALBUQUERQUE, J.A.; ARGENTON, J.; BAYER, C.; WILDNER, L.P. & KUNTZE, M.A.G. Relação de atributos do solo com a agregação de um Latossolo Vermelho sob sistemas de preparo e plantas de verão para cobertura do solo. R. Bras. Ci. Solo, 29:415-424, 2005.
• ANDRADE, A.P.; MAFRA, A.L.; BALDO, G.R.; PICCOLLA, C.D.; BERTOL, I. & ALBUQUERQUE, J.A. Physical properties of a Humic Cambisol under tillage and cropping systems after 12 years. R. Bras. Ci. Solo, 34:219-226, 2010.
• BEUTLER, A.N.; CENTURION, J.F.; CENTURION, M.A.P.C. & SILVA, A.P. Efeito da compactação na produtividade de cultivares de soja em Latossolo Vermelho. R. Bras. Ci. Solo, 30:787-794, 2006.
• BEUTLER, A.N. & CENTURION, J.F. Soil compaction and fertilization in soybean productivity. Sci. Agríc., 61:626-631, 2004.
• BRAIDA, J.A.; BAYER, C.; ALBUQUERQUE, A. & REICHERT, J.M. Matéria orgânica e seu efeito na física do solo. In: KLAUBERG FILHO, O., ed. Tópicos em ciência do solo. Viçosa, MG, Sociedade Brasileira de Ciência do Solo, 2011. v.7. p.221-278.
• BRONICK, C.J. & LAL, R. Soil structure and management: A review. Geoderma, 124:3-22, 2005.
• CAMPOS, H. Estatística experimental não paramétrica. 3.ed. Piracicaba, Escola Superior de Agricultura Luiz de Queiroz, 1979. 373p.
• CARPENEDO, V. & MIELNICZUK, J. Estado de agregação e qualidade de agregados de Latossolos Roxos, submetidos a diferentes sistemas de manejo. R. Bras. Ci. Solo, 14:99-105, 1990.
• COSTA, M.A.T.; TORMENA, C.A.; LUGÃO, S.M.B.; FIDALSKI, J.; NASCIMENTO, W.G. & MEDEIROS, F.M. Resistência do solo à penetração e produção de raízes e de forragem em diferentes níveis de intensificação do pastejo. R. Bras. Ci. Solo, 36:993-1004, 2012.
• DALBÓ, M.A.; VEIGA, M. & GARDIN, J.P.P. Desenvolvimento do sistema radicular da videira em função de porta-enxertos e de atributos físicos e químicos do solo. Agropec. Catarinense, 24:59-63, 2011.
• DE MARIA, I.C.; CASTRO, O.M. & DIAS, H.S. Atributos físicos do solo e crescimento radicular de soja em Latossolo Roxo sob diferentes métodos de preparo do solo. R. Bras. Ci. Solo, 23:703-709, 1999.
• ECHENIQUE, M.C.; APCARIAN, A.; REEB, P. & ARUANI, M.C. Growth-yield relationship of grapevine cultivars on soils with hardened layers, Alto Valle of the Río Negro, southern wine-growing region of Argentina. Agric. Técn., 67:262-270, 2007.
• ELAOUD, A. & CHEHAIBI, S. Soil compaction due to tractor traffic. J. Failure Anal. Prev., 11:539-545, 2011.
• EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Centro Nacional de Pesquisa do Solo. Levantamento de reconhecimento dos solos do Estado de Santa Catarina. Rio de Janeiro, 1998. CD-ROM. (Boletim de Pesquisa, 6)
• EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Centro Nacional de Pesquisa do Solo. Sistema brasileiro de classificação de solos. 2.ed. Brasília, 2006. 306p.
• FERRERO, A.; USOWICZ, B. & LIPIEC, J. Effects of tractor traffic on spatial variability of soil strength and water content in grass covered and cultivated sloping vineyard. Soil Till. Res., 84:127-138, 2005.
• FIDALSKI, J.; TORMENA, C.A. & SILVA, A.P. Least limiting water range and physical quality of soil under groundcover management systems in citrus. Sci. Agríc., 67:448-453, 2010.
• FUENTES LLANILLO, R.; RICHART, A.; TAVARES FILHO, J.; GUIMARÃES, M. F. & FERREIRA, R.R.M. Evolução de propriedades físicas do solo em função dos sistemas de manejo em culturas anuais. Semina: Ci. Agron., 27:205-220, 2006.
• GENRO JUNIOR, S.A.; REINERT, D.J. & REICHERT, J.M. Variabilidade temporal da resistência à penetração de um Latossolo Argiloso sob semeadura direta com rotação de culturas. R. Bras. Ci. Solo, 28:477-484, 2004.
• GUBIANI, P.I.; REINERT, D.J. & REICHERT, J.M. Método alternativo para a determinação da densidade de partículas do solo - exatidão, precisão e tempo de processamento. Ci. Rural, 36:664-668, 2006.
• HERLIHY, M. & McCARTHY, J. Association of soil test phosphorus with phosphorus fractions and adsorption characteristics. Nutr. Cycl. Agroecosyst., 75:79-90, 2006.
• KLEIN, V.A.; VIEIRA, M.L.; DURIGON, F.F.; MASSING, J.P. & FÁVERO, F. Porosidade de aeração de um Latossolo Vermelho e rendimento de trigo em plantio direto escarificado. Ci. Rural, 38:365-371, 2008.
• LIMA, C.L.R.; SILVA, A.P.; IMHOFF, S.; LIMA, H.V. & LEÃO, T.P. Heterogeneidade da compactação de um Latossolo Vermelho-Amarelo sob pomar de laranja. R. Bras. Ci. Solo, 28:409-414, 2004.
• LUCIANO, R.V. Variabilidade espacial e temporal de atributos do solo e relação com a qualidade da uva para vinificação no Planalto Catarinense. Lages, Universidade do Estado de Santa Catarina, 2012. 161p. (Tese de Doutorado)
• NESMITH, D.S. Soil compaction in double cropped wheat and soybean on Ultissol. Soil Sci. Soc. Am. J., 51:183-186, 1987.
• MEROTTO JUNIOR, A. & MUNDSTOCK, C.M. Wheat root growth as affected by soil strength. R. Bras. Ci. Solo, 23:197-202, 1999.
• MIALHE, L.G. Máquinas motoras na agricultura. São Paulo, Universidade de São Paulo, 1980. v.2. 367p.
• PANDOLFO, C.; BRAGA, H.J.; SILVA JÚNIOR, V.P.; MASSIGNAM, A.M.; PEREIRA, E.S. & THOMÉ, V.M.R. Atlas climático digital do Estado de Santa Catarina. Florianópolis, EPAGRI, 2002. CD-ROM
• PORTUGAL, A.F.; COSTA, O.D.A.; COSTA, L.M. & SANTOS, B.C.M. Atributos químicos e físicos de um Cambissolo Háplico Tb distrófico sob diferentes usos na Zona da Mata Mineira. R. Bras. Ci. Solo, 32:249-258, 2008.
• PORTUGAL, A.F.; COSTA, O.D.V. & COSTA, L.M. Propriedades físicas e químicas do solo em áreas com sistemas produtivos e mata na região da Zona da Mata mineira. R. Bras. Ci. Solo, 34:575-585, 2010.
• PROTAS, J.F.S. Uvas viníferas para processamento em regiões de clima temperado. Bento Gonçalves, Embrapa-CNPUV, 2003. (Sistema de Produção, 4)
• REICHERT, J.M.; REINERT, D.J. & BRAIDA, J.A. Qualidade dos solos e sustentabilidade de sistemas agrícolas. R. Ci. Amb., 27:29-48, 2003.
• REICHERT, J.M.; SUZUKI, L.E.A.S. & REINERT, D.J. Compactação do solo em sistemas agropecuários e florestais: Identificação, efeitos, limites críticos e mitigação. In: CERETTA, C.A.; SILVA, L.S. & REICHERT, J.M., eds. Tópicos em ciência solo. Viçosa, MG, Sociedade Brasileira de Ciência do Solo, 2007. v.5, p.49-134.
• REINERT, D.J. & REICHERT, J.M. Propriedades físicas do solo. Santa Maria, Universidade Federal Santa Maria, 2006. 18p.
• ROSA, J.D.; MAFRA, A.L.; NOHATTO, M.; FERREIRA, E.Z.; OLIVEIRA, O.L.P.; MIQUELUTTI, D.J.; CASSOL, P.C. & MEDEIROS, J.C. Atributos químicos do solo e produtividade de videiras alterados pelo manejo de coberturas verdes na Serra Gaúcha. R. Bras. Ci. Solo, 33:179-187, 2009.
• SALTON, J.C.; MIELNICZUK, J.; BAYER, C.; BOENI, M.; CONCEIÇÃO, P.C.; FABRÍCIO, A.C.; MACEDO, M.C.M. & BROCH, D.L. Agregação e estabilidade de agregados do solo em sistemas agropecuários em Mato Grosso do Sul. R. Bras. Ci. Solo, 32:11-21, 2008.
• SANTOS, G.A.S. Fundamentos da matéria orgânica do solo: Ecossistemas tropicais & subtropicais. In: SANTOS, G.A.; SILVA, L.S.; CANELLAS, L.P. & CAMARGO, F.A.O., eds. Dinâmica e função da matéria orgânica. 2.ed. Porto Alegre, Metrópole, 2008. p.7-16.
• SEIXAS, F. & SOUZA, C.R. Avaliação e efeito da compactação do solo, devido à frequência de tráfego, na produção de madeira de eucalipto. R. Árvore, 31:1047-1052, 2007.
• SILVA, M.A.S.; MAFRA, A.L.; ALBUQUERQUE, J.A.; ROSA, J.D.; BAYER, C. & MIELNICZUK, J. Propriedades físicas e teor de carbono orgânico de um Argissolo Vermelho sob distintos sistemas de uso e manejo. R. Bras. Ci. Solo, 30:327-337, 2006.
• SNEDECOR, W. & COCHRAN, W.G. Statistical methods. 8.ed. Ames, Iowa State University Press, 1989. 502p.
• SUZUKI, L.E.A.S. Compactação do solo e sua influência nas propriedades físicas. Santa Maria, Universidade Federal de Santa Maria, 2005. 149p. (Dissertação de Mestrado)
• TAYLOR, H.M.; ROBERTSON, G.M. & PARKER, J.J. Soil strength root penetration relations for medium to coarse textured soil materials. Soil Sci., 102:18-22, 1966.
• TEDESCO, M.J.; GIANELLO, C.; BISSANI, C.; BOHNEN, H. & VOLKWEISS, S.J. Análise de solo, plantas e outros materiais. 2.ed. Porto Alegre, Universidade Federal do Rio Grande do Sul, 1995. 174p. (Boletim Técnico, 5)
• TISDALL, J.M. & OADES, J.M. Organic matter and water-stable aggregates in soils. J. Soil Sci., 33:141-164, 1982.
• van DIJCK, S.J.E. & van ASCH, Th.W.J. Compaction of loamy soils due to tractor traffic in vineyards and orchards and its effect on infiltration in southern France. Soil Till. Res.; 63:141-153, 2002.
• van LIER, J.Q. Oxigenação do sistema radicular: Uma abordagem física. R. Bras.Ci. Solo, 25:233-238, 2001.
• VEIGA, M. Metodologias para coleta de amostras e análises físicas do solo. Florianópolis, EPAGRI, 2011. 52p. (Boletim Técnico, 156)
• VOCANSON A.; ROGER-ESTRADE J.; BOIZARD H. & JEUFFROY M.H. Effects of soil structure on root development of pea (Pisum sativum L.) according to the sowing date and the cultivar. Plant Soil, 281:121-135, 2006.
• XU, X.; NIEBER, J.L. & GUPTA, S.C. Compaction effects on the gas diffusion coefficients in soils. Soil Sci. Soc. Am. J., 56:1743-1750, 1992.

(1

) Parte da Dissertação de Mestrado do primeiro autor.

Datas de Publicação

• Publicação nesta coleção
  02 Out 2014
• Data do Fascículo
  Ago 2014

Histórico

• Recebido
  25 Jun 2013
• Aceito
  19 Maio 2014