Acessibilidade / Reportar erro

Cor do solo na identificação de áreas com diferentes potenciais produtivos e qualidade de café

Soil color for the identification of areas with different yield potential and coffee quality

Resumo

O objetivo deste trabalho foi verificar a eficácia da cor do solo na identificação de áreas com diferentes potenciais produtivos e qualidade de café, em Argissolo Vermelho-Amarelo desenvolvido de arenito, no Planalto Ocidental Paulista, no Estado de São Paulo. Coletaram-se 173 amostras em 39 ha, para a determinação de matiz, valor, croma e índice de avermelhamento do solo, com base nos dados da espectroscopia de reflectância difusa (ERD). Avaliaram-se, também, a produtividade e a qualidade do fruto de café, e calculou-se o índice de diagnose foliar (DRIS). A maior correlação espacial observada pelo variograma cruzado ocorreu entre matiz e DRIS, em pontos distantes em até 497 m. A menor correlação espacial foi observada em croma e produção, em pontos distantes em até 207 m. As áreas com o maior potencial de produtividade (20 a 33 sacas por hectare) apresentaram matiz entre 6,99 e 7,06, croma entre 5,00 e 5,08 e valor entre 4,40 e 4,45. As áreas com o menor potencial de produtividade (3 a 7 sacas por hectare) apresentaram matiz entre 7,27 e 7,68, croma entre 5,20 e 5,28 e valor entre 4,51 e 4,67. Nestes locais, o teor de argila foi de 16 g kg-1, os teores de P variaram entre 65 e 75 mg dm-3, e a soma de bases estava entre 56 e 58 mmolc dm-3. O uso da cor determinada por ERD é eficaz para identificar áreas com diferentes potenciais produtivos e qualidade de café, com precisão de 61 a 97%.

Termos para indexação:
Coffea arabica; cor do solo; geoestatística; pedometria

Abstract

The objective of this work was to verify soil color effectiveness for the identification of areas with different yield potentials and coffee quality, in an Ultisol developed from sandstone, in Planalto Ocidental Paulista, in the state of São Paulo, Brazil. One hundred seventy-three samples were collected from 39 ha, for the determination of hue, value, chroma, and the soil redness index, based on data of diffuse reflectance spectroscopy (DRS). Productivity and quality of coffee fruit were also evaluated, and the leaf diagnosis index (DRIS) was calculated. The highest spatial correlation, observed in the cross variogram, occurred between hue and DRIS in points up to 497 m apart. The lowest spatial correlation was observed for chroma and production, in points spaced by up to 207 m. The areas with the highest yield potential (20 to 33 bags per hectare) showed hue between 6.99 and 7.06, croma between 5.00 and 5.08, and value between 4.40 and 4.45. The areas with the lowest productivity potential (3 to 7 bags per hectare) showed hue between 7.27 and 7.68, chroma between 5.20 and 5.28, and value between 4.51 and 4.67. In these places, clay content was of 16 g kg-1, P levels varied between 65 and 75 mg dm-3, and the sum of bases was between 56 and 58 mmolc dm-3. The use of the color determined by DRS is effective to identify areas with different productive potentials and coffee quality, with 61 to 97% accuracy.

Index terms:
Coffea arabica; soil color; geostatistics; pedometrics

Introdução

O café é classificado com base em parâmetros qualitativos, como defeitos no grão e qualidade da bebida. A variabilidade natural do solo e de seus atributos, em nível detalhado (<1:5.000), é uma das causas da variação quantitativa e qualitativa da resposta do café, expressa em diferentes produtividades e qualidade de grãos. Entre as dificuldades de elaboração de modelos de previsão para a cultura do café está a caracterização da variabilidade natural dos solos e seus atributos (Sanchez et al., 2013SANCHEZ, M.G.B.; MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; PEREIRA, G.T. Delineation of specific management areas for coffee cultivation based on the soil-relief relationship and numerical classification. Precision Agriculture, v.14, p.201-214, 2013. DOI: 10.1007/s11119-012-9288-z.
https://doi.org/10.1007/s11119-012-9288-...
).

Na mesma área, podem ser encontrados cultivos com baixa produtividade e qualidade de grão e que, no momento da colheita, são agrupados aos de alta produção e qualidade (Sanchez et al., 2013SANCHEZ, M.G.B.; MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; PEREIRA, G.T. Delineation of specific management areas for coffee cultivation based on the soil-relief relationship and numerical classification. Precision Agriculture, v.14, p.201-214, 2013. DOI: 10.1007/s11119-012-9288-z.
https://doi.org/10.1007/s11119-012-9288-...
). Essa mistura pode provocar uma diminuição do rendimento final de até 2,5 vezes, considerando-se o ganho das sacas beneficiadas e o custo com adubação (Pollo, 2013POLLO, G.Z. Suscetibilidade magnética, atributos do solo e da planta na discriminação de áreas de manejo específico na cultura do café. 2013. 58p. Dissertação (Mestrado) - Universidade Estadual Paulista, Jaboticabal.). Uma alternativa para auxiliar o planejamento estratégico e a tomada de decisão do cafeicultor pode ser a definição de áreas de manejo específico, que apresentem máxima homogeneidade dos atributos do solo.

Ferramentas computacionais - como o modelo digital de elevação (MDE), a geoestatística, a análise multivariada e a integração destas ferramentas - têm sido cada vez mais utilizadas em ciências agrárias para auxiliar a identificação de áreas homogêneas (Siqueira et al., 2010SIQUEIRA, D.S.; MARQUES JR., J.; PEREIRA, G.T. The use of landforms to predict the variability of soil and orange attributes. Geoderma, v.155, p.55-66, 2010. DOI: 10.1016/j.geoderma.2009.11.024.
https://doi.org/10.1016/j.geoderma.2009....
; Sanchez et al., 2013SANCHEZ, M.G.B.; MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; PEREIRA, G.T. Delineation of specific management areas for coffee cultivation based on the soil-relief relationship and numerical classification. Precision Agriculture, v.14, p.201-214, 2013. DOI: 10.1007/s11119-012-9288-z.
https://doi.org/10.1007/s11119-012-9288-...
).

A evolução computacional e a formação de pessoas para utilizar essas ferramentas no planejamento agrícola ocorre em paralelo a outras análises de campo e laboratório, que têm sido testadas para a identificação de áreas de manejo específico, entre as quais a quantificação da cor do solo por espectroscopia de refletância difusa (ERD) (Marques Jr. et al., 2014MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; TEIXEIRA, D.D.B.; BARRÓN, V.; TORRENT, J. Magnetic susceptibility and diffuse reflectance spectroscopy to characterize the spatial variability of soil properties in a Brazilian Haplustalf. Geoderma, v.219/220, p.63-71, 2014. DOI: 10.1016/j.geoderma.2013.12.007.
https://doi.org/10.1016/j.geoderma.2013....
). Com à evolução dos sensores e de técnicas mais refinadas de quantificação, a cor ganha novas aplicações e perspectivas para a caracterização da variabilidade dos atributos do solo de caráter agronômico, como o teor de argila, o fósforo adsorvido e a saturação por bases, por exemplo (Marques Jr. et al., 2014MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; TEIXEIRA, D.D.B.; BARRÓN, V.; TORRENT, J. Magnetic susceptibility and diffuse reflectance spectroscopy to characterize the spatial variability of soil properties in a Brazilian Haplustalf. Geoderma, v.219/220, p.63-71, 2014. DOI: 10.1016/j.geoderma.2013.12.007.
https://doi.org/10.1016/j.geoderma.2013....
; Camargo et al., 2015CAMARGO, L.A.; MARQUES JÚNIOR ., J.; BARRÓN, V.; ALLEONI, L.R.F.; BARBOSA, R.S.; PEREIRA, G.T. Mapping of clay, iron oxide and adsorbed phosphate in Oxisols using diffuse reflectance spectroscopy. Geoderma, v.251/252, p.124-132, 2015. DOI: 10.1016/j.geoderma.2015.03.027.
https://doi.org/10.1016/j.geoderma.2015....
). Os parâmetros da cor do solo (matiz, valor, croma), quando determinados pela leitura da reflectância do solo na faixa do visível (350 a 750 nm), podem ser utilizados na estimativa do teor de óxidos de ferro, como hematita e goethita (Viscarra Rossel et al., 2010VISCARRA ROSSEL, R.A., BUI, E.N.; CARITAT, P. de; MCKENZIE, N.J. Mapping iron oxides and the color of Australian soil using visible-near-infrared reflectance spectra. Journal of Geophysical Research, v.115, p.1-13, 2010. DOI: 10.1029/2009JF001645.
https://doi.org/10.1029/2009JF001645...
; Bahia et al., 2015BAHIA, A.S.R. de S.; MARQUES JR., J.; SIQUEIRA, D.S. Procedures using diffuse reflectance spectroscopy for estimating hematite and goethite in Oxisols of São Paulo, Brazil. Geoderma Regional, v.5, p.150-156, 2015. DOI: 10.1016/j.geodrs.2015.04.006.), carbono total (Xie et al., 2011XIE, H.T.; YANG, X.M.; DRURY, C.F.; YANG, J.Y.; ZHANG, X.D. Predicting soil organic carbon and total nitrogen using mid- and near-infrared spectra for Brookston clay loam soil in Southwestern Ontario, Canada. Canadian Journal of Soil Science, v.91, p.53-63, 2011. DOI: 10.4141/cjss10029.
https://doi.org/10.4141/cjss10029...
) e outras importantes propriedades do solo que covariam com ela, tanto do ponto de vista pedológico quanto do agronômico (Barrón et al., 2000BARRÓN, V.; MELLO, J. W. V.; TORRENT, J. Caracterização de óxidos de ferro em solos por espectroscopia de reflectância difusa. In: NOVAIS, R.F. de; ALVAREZ V., V.H.; SCHAEFER, C.E.G.R., (Ed.). Tópicos em ciência do solo. Viçosa: Sociedade Brasileira de Ciência do Solo, 2000. v.1, p.139-162.; Barrón & Torrent, 2013BARRÓN, V.; TORRENT, J. Iron, manganese and aluminium oxides and oxyhydroxides. In: NIETO, F.; LIVI, K.J.T.; OBERTI, R. (Ed.). Minerals at the nanoscale. London: European Mineralogical Union: Mineralogical Society of Great Britain & Ireland, 2013. p.297-336. (European Mineralogical Union Notes in Mineralogy, v.14). DOI: 10.1180/EMU-notes.14.9.
https://doi.org/10.1180/EMU-notes.14.9...
). A cor pode ser determinada de diferentes maneiras (Torrent & Barrón, 2008TORRENT, J.; BARRÓN, V. Diffuse reflectance spectroscopy. In: ULERY, A.L.; DREES, L.R. (Ed.). Methods of soil analysis: part 5: mineralogical methods. Madison: Soil Science Society of America, 2008. p.367-387. (SSSA Book Series, 5). DOI: 10.2136/sssabookser5.5.c13.
https://doi.org/10.2136/sssabookser5.5.c...
). Em ciências agrárias, é comum utilizar-se a comparação visual de padrões de cor conforme a caderneta de Munsell (Soil Color Charts) que, no entanto, apresenta alguma dificuldade de comparação dos resultados, em razão das diferenças de percepção visual entre os avaliadores (Barrón et al., 2000BARRÓN, V.; MELLO, J. W. V.; TORRENT, J. Caracterização de óxidos de ferro em solos por espectroscopia de reflectância difusa. In: NOVAIS, R.F. de; ALVAREZ V., V.H.; SCHAEFER, C.E.G.R., (Ed.). Tópicos em ciência do solo. Viçosa: Sociedade Brasileira de Ciência do Solo, 2000. v.1, p.139-162.; Campos & Demattê, 2004CAMPOS, R.C.; DEMATTÊ, J.A.M. Cor do solo: uma abordagem da forma convencional de obtenção em oposição à automatização do método para fins de classificação de solos. Revista Brasileira de Ciência do Solo, v.28, p.853-863, 2004. DOI: 10.1590/S0100-06832004000500008.
https://doi.org/10.1590/S0100-0683200400...
). Apesar de o sistema Munsell ter sido proposto inicialmente em ciências agrárias, para fins pedológicos na classificação de solos, ele pode ser utilizado para a definição de áreas de manejo específico para as culturas da cana-de-açúcar (Marques Jr. et al., 2014MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; TEIXEIRA, D.D.B.; BARRÓN, V.; TORRENT, J. Magnetic susceptibility and diffuse reflectance spectroscopy to characterize the spatial variability of soil properties in a Brazilian Haplustalf. Geoderma, v.219/220, p.63-71, 2014. DOI: 10.1016/j.geoderma.2013.12.007.
https://doi.org/10.1016/j.geoderma.2013....
) e da soja (Resende et al., 2014RESENDE, J.M. do A.; MARQUES JÚNIOR, J.; MARTINS FILHO, M.V.; DANTAS, J.S.; SIQUEIRA, D.S.; TEIXEIRA, D.D.B. Variabilidade espacial de atributos de solos coesos do leste Maranhense. Revista Brasileira de Ciência do Solo, v.38, p.1077-1090, 2014.), estudos de gênese de solo (Dantas et al., 2014aDANTAS, J.S.; MARQUES JÚNIOR, J.; MARTINS FILHO, M.V.; RESENDE, J.M. do A.; CAMARGO, L.A.; BARBOSA, R.S. Gênese de solos coesos do leste maranhense: relação solo-paisagem. Revista Brasileira de Ciência do Solo, v.38, p.1039-1050, 2014a. DOI: 10.1590/S0100-06832014000400001.
https://doi.org/10.1590/S0100-0683201400...
), quantificação do potencial de erosão (Dantas et al., 2014aDANTAS, J.S.; MARQUES JÚNIOR, J.; MARTINS FILHO, M.V.; RESENDE, J.M. do A.; CAMARGO, L.A.; BARBOSA, R.S. Gênese de solos coesos do leste maranhense: relação solo-paisagem. Revista Brasileira de Ciência do Solo, v.38, p.1039-1050, 2014a. DOI: 10.1590/S0100-06832014000400001.
https://doi.org/10.1590/S0100-0683201400...
), e identificação de áreas com diferentes potenciais de adsorção de fósforo (Camargo et al., 2015CAMARGO, L.A.; MARQUES JÚNIOR ., J.; BARRÓN, V.; ALLEONI, L.R.F.; BARBOSA, R.S.; PEREIRA, G.T. Mapping of clay, iron oxide and adsorbed phosphate in Oxisols using diffuse reflectance spectroscopy. Geoderma, v.251/252, p.124-132, 2015. DOI: 10.1016/j.geoderma.2015.03.027.
https://doi.org/10.1016/j.geoderma.2015....
; Peluco et al., 2015PELUCO, R.G.; MARQUES JÚNIOR, J.; SIQUEIRA, D.S.; PEREIRA, G.T; BARBOSA, R.S.; TEIXEIRA, D. de B. Mapeamento do fósforo adsorvido por meio da cor e da suscetibilidade magnética do solo. Pesquisa Agropecuária Brasileira, v.50, p.259-266, 2015. DOI: 10.1590/S0100-204X2015000300010.
https://doi.org/10.1590/S0100-204X201500...
).

Nesse contexto, o emprego da ERD para a determinação dos parâmetros da cor tem-se mostrado útil para a caracterização quantitativa e mais precisa dos solos e mostrado que é possível explorar melhor esta propriedade a partir de um instrumento não muito complexo, além de promover o desenvolvimento científico e aplicações neste campo de estudo em ciência do solo (Barrón et al., 2000BARRÓN, V.; MELLO, J. W. V.; TORRENT, J. Caracterização de óxidos de ferro em solos por espectroscopia de reflectância difusa. In: NOVAIS, R.F. de; ALVAREZ V., V.H.; SCHAEFER, C.E.G.R., (Ed.). Tópicos em ciência do solo. Viçosa: Sociedade Brasileira de Ciência do Solo, 2000. v.1, p.139-162.; Torrent & Barrón, 2008TORRENT, J.; BARRÓN, V. Diffuse reflectance spectroscopy. In: ULERY, A.L.; DREES, L.R. (Ed.). Methods of soil analysis: part 5: mineralogical methods. Madison: Soil Science Society of America, 2008. p.367-387. (SSSA Book Series, 5). DOI: 10.2136/sssabookser5.5.c13.
https://doi.org/10.2136/sssabookser5.5.c...
; Barrón & Torrent, 2013BARRÓN, V.; TORRENT, J. Iron, manganese and aluminium oxides and oxyhydroxides. In: NIETO, F.; LIVI, K.J.T.; OBERTI, R. (Ed.). Minerals at the nanoscale. London: European Mineralogical Union: Mineralogical Society of Great Britain & Ireland, 2013. p.297-336. (European Mineralogical Union Notes in Mineralogy, v.14). DOI: 10.1180/EMU-notes.14.9.
https://doi.org/10.1180/EMU-notes.14.9...
).

O objetivo deste trabalho foi verificar a cor de solo, quantificada por espectroscopia de reflectância difusa (ERD), para a identificação de áreas com diferentes potenciais produtivos e qualidade de café, em um Argissolo Vermelho-Amarelo desenvolvido de arenito, no Planalto Ocidental Paulista.

Material e Métodos

O trabalho foi realizado em uma área de 39 ha, no Município de Matão, no sudoeste do Estado de São Paulo (Figura 1) - a 21º36'S e 48º29'W, com altitude entre 590 e 615 m. Segundo a classificação de Thornthwaite (1948)THORNTHWAITE, C.W. An approach toward a rational classification of climate. Geographical Review, v.38, p.55-94, 1948. DOI: 10.2307/210739.
https://doi.org/10.2307/210739...
, o clima local é definido como B1rB'4a', tipo mesotérmico úmido, com pequena deficiência hídrica, com evapotranspiração de verão abaixo de 70% da evapotranspiração anual.

Figura 1
Localização da área e malha amostral (A); corte para representação do perfil de elevação e espessura do horizonte A+E (B); modelo digital de elevação que representa a variabilidade espacial da espessura do horizonte A+E (C).

O material geológico está relacionado aos arenitos do Grupo Bauru. Geomorfologicamente, a área está inserida no Planalto Ocidental Paulista, que ocupa aproximadamente metade da extensão territorial do Estado de São Paulo - cerca de 13 milhões de hectares (IPT, 1981IPT. INSTITUTO DE PESQUISAS TECNOLÓGICAS DO ESTADO DE SÃO PAULO. Mapa geomorfológico do Estado de São Paulo: mapas. São Paulo, 1981. 94p. (IPT Publicação, 1183; IPT Monografia, 5).).

O solo da área de estudo foi classificado como Argissolo Vermelho-Amarelo distrófico (Santos et al., 2006SANTOS, H.G. dos; JACOMINE, P.K.T.; ANJOS, L.H.C. dos; OLIVEIRA, V.A. de; OLIVEIRA, J.B. de; COELHO, M.R.; LUMBRERAS, J.F.; CUNHA, T.J.F. (Ed.). Sistema brasileiro de classificação de solos. 2.ed. Brasília: Embrapa, 2006. 306p.), com teor médio de Fe2O3 (obtido por ataque sulfúrico) de 30 g kg-1. À profundidade de 0,0-0,2 m, o teor médio de argila é de 160 g kg-1, a soma de bases é de 52 mmolc dm-3, e a suscetibilidade magnética é de 1,50 10-6 m3 kg-1. Realizou-se, também, a avaliação da espessura dos horizontes A+E, por meio de três tradagens profundas sob a projeção da copa da planta central da unidade amostral (Figura 1 A) e, posteriormente, calculou-se a média das determinações. Estes valores foram utilizados na modelagem do variograma e elaboração do mapa de espessura do horizonte A+E, para complementar a caracterização da área e auxiliar a compreensão das relações de causa e efeito entre solo e paisagem. A mudança textural em decorrência do maior teor de argila, na camada mais profunda, foi utilizada como parâmetro indicativo da profundidade dos horizontes A+E, conforme realizado por Sanchez et al. (2009SANCHEZ, R.B.; MARQUES JÚNIOR, J.; SOUZA, Z.M. de; PEREIRA, G.T.; MARTINS FILHO, M.V. Variabilidade espacial de atributos do solo e de fatores de erosão em diferentes pedoformas. Bragantia, v.68, p.1095-1103, 2009. DOI: 10.1590/S0006-87052009000400030.
https://doi.org/10.1590/S0006-8705200900...
, 2013)SANCHEZ, M.G.B.; MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; PEREIRA, G.T. Delineation of specific management areas for coffee cultivation based on the soil-relief relationship and numerical classification. Precision Agriculture, v.14, p.201-214, 2013. DOI: 10.1007/s11119-012-9288-z.
https://doi.org/10.1007/s11119-012-9288-...
.

A área é cultivada há 25 anos com Coffea arabica 'Catuaí Amarelo IAC 62'. À época do estudo, as plantas estavam com 7 anos de idade, dispostas em espaçamento 3,5 x 0,5 m, cultivadas desde o início em sistema de fertirrigação por gotejamento. Coletaram-se amostras de solo à profundidade de 0,0-0,20 m, a intervalos regulares de 45 m, sob a projeção da copa da planta central, num total de 173 pontos de amostragem. Para a avaliação da produtividade e da qualidade do fruto, definiram-se unidades amostrais de nove plantas, tendo-se avaliado o total de 1.557 plantas.

As amostras de solo coletadas na malha e nas trincheiras foram submetidas à análise granulométrica, em solução de NaOH 0,1 N como dispersante químico, com agitação em aparato de baixa rotação (Day, 1965DAY, P. R. Particle fraction and particle-size analysis. In: BLACK, C.A. (Ed.). Methods of soil analysis. Madison: American Society of Agronomy, 1965. v.1, p.545-566. (ASA. Agronomy, 9).). A argila foi determinada pelo método da pipeta. Cálcio, magnésio e potássio trocáveis e o fósforo disponível foram extraídos pelo método da resina trocadora de íons (Raij et al., 1987RAIJ, B. van; QUAGGIO, J.A.; CANTARELLA, J.; FERREIRA, M.E.; LOPES, A.S.; BATAGLIA, O.C. Análise química do solo para fins de fertilidade. Campinas: Fundação Cargill, 1987. 170p.), e o carbono orgânico foi extraído conforme Claessen (1997)CLAESSEN, M.E.C. (Org.). Manual de métodos de análise de solo. 2.ed. rev. e atual. Rio de Janeiro,: EMBRAPA-CNPS, 1997. 212p..

Para a obtenção dos espectros de reflectância difusa utilizou-se o sensor Lambda 950 UV/VIS/NIR. Usou-se aproximadamente 0,5 g de terra fina seca ao ar (TFSA), moída em almofariz de ágata até a obtenção de coloração constante. Os valores de reflectância foram determinados a cada 0,5 nm na faixa do visível, no intervalo de 380 a 780 nm (Torrent & Barrón, 2008TORRENT, J.; BARRÓN, V. Diffuse reflectance spectroscopy. In: ULERY, A.L.; DREES, L.R. (Ed.). Methods of soil analysis: part 5: mineralogical methods. Madison: Soil Science Society of America, 2008. p.367-387. (SSSA Book Series, 5). DOI: 10.2136/sssabookser5.5.c13.
https://doi.org/10.2136/sssabookser5.5.c...
).

Após a obtenção dos espectros de reflectância das amostras de solo, determinaram-se os valores triestímulo XYZ, definidos pela Comisión Internacional de L'Eclairage - CIE (Wyszecki & Stiles, 1982WYSZECKI, G.; STILES, W.S. Color science: concepts and methods, quantitative data and formulae. 2.ed. New York: John Wiley & Sons, 1982. 968p.) e, a partir destas coordenadas, foram deduzidos os valores de Munsell referentes a matiz (H), croma (C) e valor (V), pelo programa Munsell Conversion, versão 6.4. Os valores de H, C e V, conforme método de Munsell, foram utilizados no cálculo do índice de avermelhamento (IAV), pela equação 1 (Torrent & Barrón, 2008TORRENT, J.; BARRÓN, V. Diffuse reflectance spectroscopy. In: ULERY, A.L.; DREES, L.R. (Ed.). Methods of soil analysis: part 5: mineralogical methods. Madison: Soil Science Society of America, 2008. p.367-387. (SSSA Book Series, 5). DOI: 10.2136/sssabookser5.5.c13.
https://doi.org/10.2136/sssabookser5.5.c...
), a seguir: IAV = (10 - H)×C/V, em que: V é o valor Munsell; C é o valor numérico do croma; e H é o valor numérico que corresponde ao matiz.

A produtividade foi estimada com base na produção das plantas da unidade amostral. A classificação do tamanho dos grãos foi realizada por um jogo de peneiras que os separa pelo tamanho e pela forma. Consideraram-se para a análise os grãos retidos na peneira com abertura 17 mm.

Os resultados das análises foliares do cafeeiro foram submetidos à metodologia do sistema integrado de diagnose e recomendação (DRIS). Esta técnica compara o balanço nutricional dos nutrientes na folha ao padrão nutricional ótimo, baseado em altos índices de produtividade (Beaufils, 1973BEAUFILS, E.R. Diagnosis and recommendation integrated system (DRIS): a general scheme of experimentation and calibration based on principles developed from research in plant nutrition. Pietermaritzburg: University of Natal, 1973. 132p. (Soil Science Bulletin, 1).). Para o cálculo pelo DRIS, utilizaram-se informações sobre os teores de N, P, K, Ca, Mg, S, B, Cu, Fe, Mn e Zn, avaliados nas folhas de 1.557 plantas da unidade amostral.

Calcularam-se a média, o desvio-padrão (DP) e o coeficiente de variação (CV) para as variáveis do solo e da planta. As médias das variáveis do solo e da planta, em relação à espessura do horizonte A+E, foram comparadas pelo teste de Tukey, a 5% de probabilidade.

A dependência espacial foi avaliada pelo variograma simples e por mapas de krigagem (Isaaks & Srivastava, 1989ISAAKS, E.H.; SRIVASTAVA, R.M. An introduction to applied geoestatistics. New York: Oxford University Press, 1989. 560p.). A correlação espacial entre os parâmetros da cor e as variáveis da planta de café foi avaliada pelo variograma cruzado, que é ajustado somente para as variáveis que apresentam correlação entre si (Deutsch & Journel, 1998DEUTSCH, C.V.; JOURNEL, A.G. GSLIB: Geostatistical Software Library and User's Guide. 2.ed. Oxford: Oxford University Press, 1998. 366p.). A técnica do variograma cruzado foi utilizada para avaliar a semelhança dos modelos matemáticos (variogramas) que descrevem a variabilidade espacial das variáveis do solo e da cultura. A dependência espacial entre as variáveis foi estimado pela equação 2 (Deutsch & Journel, 1998DEUTSCH, C.V.; JOURNEL, A.G. GSLIB: Geostatistical Software Library and User's Guide. 2.ed. Oxford: Oxford University Press, 1998. 366p.), a seguir:

γ ̂ zy ( h ) = 1 2 N ( h ) i 1 N ( h ) [ z ( x i ) z ( x i + h ) ] [ y ( x i ) y ( x i + h ) ] ,

em que: γ̂zy (h) é a semivariância experimental cruzada, para uma distância de separação h; z(xi) é o valor da variável principal (aquela a ser estimada), no ponto i; y(xi) é o valor da variável secundária no ponto i; e N(h) é o número de pares de pontos separados pela distância h.

O variograma simples é um caso particular do variograma cruzado, em que a semivariância é calculada somente para uma propriedade e, assim, é considerado como uma medida da autocorrelação espacial da variável.

Os dados referentes aos 173 pontos de amostragem de solo foram hierarquizados em grupos, com base no percentil da produtividade e do parâmetro qualitativo percentagem de grãos retidos na peneira de 17 mm (Tabela 1).

Tabela 1
Valores médios dos atributos da cor resultantes dos grupos de resposta da planta.

A média dos resultados referentes ao valor, croma, matiz e IAV, para cada um dos grupos de resposta da planta, foi utilizada para construir gráficos de dispersão e representar a relação das variáveis dependentes (resposta da planta), em função das variáveis independentes (atributos da cor). Utilizou-se o método dos mínimos quadrados, para ajustar o modelo de regressão ao conjunto de pontos das variáveis dependentes e independentes.

Construiu-se o triângulo de cor para identificar locais com diferentes potenciais de produtividade, com base no conceito do diagrama de cromaticidade elaborado a partir dos valores triestímulos (XYZ). A construção do triângulo foi feita em duas etapas, conforme a seguir. Na primeira etapa, os intervalos dos modelos de regressão para valor (4,40 a 4,67), croma (5,00 a 5,28) e matiz (6,99 a 7,68) foram subdivididos em quatro grupos, com base no intervalo do modelo de regressão da produtividade - 33 a 20 sacas ha-1, 20 a 13 sacas ha-1, 13 a 7 sacas ha-1 e 7 a 3 sacas ha-1 (Figura 2).

Figura 2
Triângulos de cor que relacionam a amplitude encontrada nos modelos de regressão da produtividade com o valor, o croma e o matiz.

Na segunda etapa, foi feita a integração dos triângulos para valor, croma e matiz. Cada intersecção das componentes da cor corresponde aos locais com diferentes características físicas e químicas que, por sua vez, condicionam diferentes potenciais de produtividade de café na área estudada. Esse tipo de representação gráfica elaborada a partir da ERD pode mostrar a interação existente entre os três parâmetros da cor, no sistema Munsell, e uma quarta variável de interesse agronômico. Assim, esta é uma alternativa rápida e viável, para a elaboração de triângulos de cor para diferentes regiões, com o objetivo de auxiliar a identificação de áreas de manejo no campo.

Resultados e Discussão

Os locais com menor espessura do horizonte A+E (entre 39 e 45 cm) apresentaram os maiores valores de DRIS (8,28), IAV (3,62) e produtividade (26,13 sc ha-1), e os menores valores para matiz (6,94), valor (4,33), croma (5,10), percentagem de grãos retidos na peneira de 17 mm (31,33%), argila (15,66%), carbono orgânico (28,44 g dm-3), fósforo (49,04 mg dm-3) e soma de bases (46,52 mmolc dm-3). O maior valor do IAV (3,62), nos locais com menor espessura dos horizontes A+E, está relacionado ao maior teor de argila nessa camada, o que indica que estes solos foram relativamente mais erodidos e, assim, estão mais próximos do horizonte B textural, que é mais argiloso.

Com exceção do croma, a média de todas as outras variáveis relacionadas à cor do solo e à resposta da planta diferiram em relação às áreas e à espessura dos horizontes A+E. Observou-se maior desbalanço nutricional na parte superior da área de estudo, com maior IAV e menor espessura do horizonte A+E; nesse local, o valor de DRIS está mais distante de zero (8,28) do que na parte inferior da área de estudo (6,88), o que indica que parte dos nutrientes foi exportada para o desenvolvimento do grão. A maior produtividade na parte superior da área (26,13 sc ha-1), em comparação à parte inferior (13,76 sc ha-1), está de acordo com essa explicação.

A maioria das raízes do café, com maior potencial de absorção (menos de 3 mm de diâmetro), está localizada sob o dossel das plantas em até 0,20 m de profundidade (Motta et al., 2006MOTTA, A.C.V.; NICK, J.A.; YORINORI, G.T.; SERRAT, B.M. Distribuição horizontal e vertical da fertilidade do solo e das raízes de cafeeiro (Coffea arabica L.) cultivar Catuaí. Acta Scientiarum. Agronomy, v.28, p.455-463, 2006. DOI: 10.4025/actasciagron.v28i4.758.
https://doi.org/10.4025/actasciagron.v28...
). Mediante os resultados expostos, os locais na área com solo mais avermelhado (maior IAV) indicam menor espessura do horizonte A+E, ambiente onde há melhor absorção de água e nutrientes, o que aumenta o potencial produtivo. Nos locais em que a espessura do horizonte A+E é maior ocorre maior percentagem de grãos retidos na peneira de 17 mm, o que pode ser atribuído aos resultados da soma de bases (57,58 mmolc dm-3), carbono orgânico (40,22 g dm-3) e fósforo (69,35 mg dm-3), que confirmam a maior fertilidade do solo nesses locais.

Os locais menos avermelhados (maior espessura do horizonte A+E) apresentaram os maiores valores de CV dos atributos estudados, com exceção do valor, percentagem de grãos retidos na peneira de 17 mm, fósforo e carbono orgânico. Sanchez et al. (2013)SANCHEZ, M.G.B.; MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; PEREIRA, G.T. Delineation of specific management areas for coffee cultivation based on the soil-relief relationship and numerical classification. Precision Agriculture, v.14, p.201-214, 2013. DOI: 10.1007/s11119-012-9288-z.
https://doi.org/10.1007/s11119-012-9288-...
também relatam melhores ambientes de produção para a cultura do café, em locais com menor espessura do horizonte A+E.

O modelo esférico do variograma foi o que melhor se ajustou aos atributos matiz, croma, IAV, produtividade, DRIS e argila, enquanto o modelo exponencial foi o que melhor se ajustou aos atributos valor, percentagem de grãos retidos na peneira de 17 mm, fósforo, carbono orgânico e soma de bases (Figura 3). O modelo esférico indica mudanças bruscas do padrão de variabilidade entre pontos mais distantes (Cambardella et al., 1994CAMBARDELLA, C.A.; MOORMAN, T.B.; NOVAK, J.M.; PARKIN, T.B.; KARLEN, D.L.; TURCO, R.F.; KONOPKA, A.E. Field-scale variability of soil properties in Central Iowa soils. Soil Science Society of America Journal, v.58, p.1501-1511, 1994. DOI: 10.2136/sssaj1994.03615995005800050033x.
https://doi.org/10.2136/sssaj1994.036159...
); o exponencial indica padrões de variabilidade mais irregulares em distâncias curtas (Isaaks & Srivastava, 1989ISAAKS, E.H.; SRIVASTAVA, R.M. An introduction to applied geoestatistics. New York: Oxford University Press, 1989. 560p.).

Figura 3
Variogramas e mapas de krigagem para matiz, valor, croma e índice de avermelhamento (IAV), à profundidade de 0,0-0,20 m, e percentagem de grãos retidos na peneira de 17 mm, produtividade e DRIS. Esf., modelo esférico; Exp., modelo exponencial; Co, efeito pepita; Co + C1, patamar; a, alcance; C0/(C0+C1), grau de dependência espacial; R2, coeficiente de determinação.

Os valores de alcance variaram de 150 a 479 m, para os componentes da cor do solo, e de 150 a 180 m para a resposta da planta. Autores que realizaram estudos quanto à variabilidade espacial, em Argissolo originado de arenito do Planalto Ocidental Paulista, também encontraram valores de alcance dentro da mesma faixa de variação - 150 a 479 m -, em área de citros (284 a 350 m) (Siqueira et al., 2010SIQUEIRA, D.S.; MARQUES JR., J.; PEREIRA, G.T. The use of landforms to predict the variability of soil and orange attributes. Geoderma, v.155, p.55-66, 2010. DOI: 10.1016/j.geoderma.2009.11.024.
https://doi.org/10.1016/j.geoderma.2009....
) e cana-de-açúcar (160 a 410 m) (Sanchez et al., 2009SANCHEZ, R.B.; MARQUES JÚNIOR, J.; SOUZA, Z.M. de; PEREIRA, G.T.; MARTINS FILHO, M.V. Variabilidade espacial de atributos do solo e de fatores de erosão em diferentes pedoformas. Bragantia, v.68, p.1095-1103, 2009. DOI: 10.1590/S0006-87052009000400030.
https://doi.org/10.1590/S0006-8705200900...
). Essa similaridade do padrão de distribuição espacial pode estar associada à variabilidade inicial causada pelo material de origem (Schemiko et al., 2014SCHEMIKO, D.C.B.; VESELY, F.F.; FERNANDES, L.A.; SOWEK, G.A. Distinção dos elementos deposicionais fluviais, eólicos e lacustres do Membro Morro Pelado, Permiano Superior da Bacia do Paraná. Geologia USP. Série Científica, v.14, p.29-46, 2014. DOI: 10.5327/Z1519-874X201400030003.
https://doi.org/10.5327/Z1519-874X201400...
; Fernandes & Ribeiro, 2015FERNANDES, L.A.; RIBEIRO, C.M.M. Evolution and palaeoenvironment of the Bauru Basin (Upper Cretaceous, Brazil). Journal of South American Earth Sciences, v.61, p.71-90, 2015. DOI: 10.1016/j.jsames.2014.11.007.
https://doi.org/10.1016/j.jsames.2014.11...
) comum a estes trabalhos.

Observando-se os mapas de matiz, valor e croma, nota-se que os menores valores estão localizados na área com menor espessura do horizonte A+E (Figura 1 B e C). Estes resultados são complementares aos da Tabela 2 e indicam os limites entre os ambientes com diferentes características e as classes de variabilidade. Produtividade e percentagem de grãos de café retidos na peneira de 17 mm apresentaram padrão inverso de variabilidade. Os locais com menor espessura do horizonte (39 a 45 cm) tiveram maior produtividade (36 a 60 sacas ha-1) e de frutos menores (menos de 16% dos frutos foram retidos na peneira de 17 mm).

Tabela 2
Valores médios, desvio-padrão e coeficiente de variação (CV) dos componentes da cor do solo, à profundidade de 0,0-0,2 m, e características da cultura do café (produtividade e parâmetros qualitativos, n = 173) relacionados à espessura dos horizontes A+E(1) (1) Médias seguidas de letras iguais, nas linhas, não diferem entre si, pelo teste de Tukey, a 5% de probabilidade. .

Observaram-se as seguintes relações de dependência positiva: do matiz em relação à percentagem de grãos retidos na peneira de 17 mm (449 m); do valor, à percentagem de grãos retidos na peneira de 17 mm (410 m); do croma, à produtividade (207 m); do IAV, à produtividade (259 m); e do IAV, à argila (473 m). Os variogramas cruzados dos pares - matiz-produtividade, matiz-DRIS, IAV-percentagem de grãos retidos na peneira de 17 mm, IAV-carbono orgânico, IAV-fósforo, IAV-soma de bases - apresentaram dependência espacial negativa (Figura 4). Estes resultados indicam que os componentes da cor obtidos pela ERD têm potencial para seu uso como covariáveis na cokrigagem e para auxiliar o detalhamento de mapas de produtividade e do tamanho do grão na cultura do café.

Figura 4
Variogramas cruzados dos atributos de cor, à profundidade de 0,0-0,20 m, produtividade, qualidade do grão de café e atributos físicos e químicos do solo. Esf., modelo esférico; Exp., modelo exponencial; Gauss, modelo gaussiano; Co, efeito pepita; Co + C1: patamar; a, alcance; C0/(C0+C1), grau de dependência espacial; R2, coeficiente de determinação.

Modelos quadráticos foram encontrados para os atributos analisados. A produtividade pode ser estimada pelos componentes da cor do solo à profundidade de 0,0-0,2 m, com precisão de 66 a 95% (Figura 5). Quanto às características do grão, os atributos da cor estimam o tamanho do grão com precisão de 88%. Para o DRIS, a precisão é de 82%. Em solos originados de sedimentos do Grupo Barreiras, os componentes da cor do solo foram utilizados em relações diretas para estimar a fertilidade, com precisão de 85% (Resende, 2013RESENDE, J.M. do A. Caracterização pedométrica de atributos de Argissolos coesos do leste maranhense. 2013. 83p. Tese (Doutorado) - Universidade Estadual Paulista, Jaboticabal.); a erodibilidade foi estimada com precisão de 70 a 73%; e a acurácia, com 53% (Dantas et al., 2014bDANTAS, J.S.; MARTINS FILHO, M.V.; MARQUES JÚNIOR, J.; RESENDE, J.M. do A.; TEIXEIRA, D.D.B.; BARBOSA, R.S.; SIQUEIRA, D.S. Coeficiente de erodibilidade em sulcos e entressulcos de Argissolos coesos estimado pela cor do solo. Pesquisa Agropecuária Brasileira, v.49, p.700-707, 2014b. DOI: 10.1590/S0100-204X2014000900006.
https://doi.org/10.1590/S0100-204X201400...
). Em solos com teor de Fe2O3 de 20 a 200 g kg-1, o IAV apresentou baixa precisão (4%) na estimativa de áreas com diferentes potenciais de adsorção de fósforo (Peluco et al., 2015PELUCO, R.G.; MARQUES JÚNIOR, J.; SIQUEIRA, D.S.; PEREIRA, G.T; BARBOSA, R.S.; TEIXEIRA, D. de B. Mapeamento do fósforo adsorvido por meio da cor e da suscetibilidade magnética do solo. Pesquisa Agropecuária Brasileira, v.50, p.259-266, 2015. DOI: 10.1590/S0100-204X2015000300010.
https://doi.org/10.1590/S0100-204X201500...
).

Figura 5
Modelos de regressão dos atributos da cor do solo, à profundidade de 0,0-0,20 m, com a produtividade, qualidade do grão de café e atributos físicos e químicos do solo.

A representação gráfica com as três componentes da cor mostra que as áreas com maior potencial de produtividade (20 a 22 sacas ha-1) estão relacionadas ao valor entre 4,45 e 4,40, matiz entre 7,09 e 6,99 e croma entre 5,0 e 5,08. Estas variações dos atributos da cor expressam indiretamente a variação do teor de argila entre 7 e 35%, teores de fósforo entre 8 e 187 mg dm-3 e soma de bases entre 14 e 216 mmolc dm-3. Os locais com menor produtividade (3 a 7 sacas ha-1) estão relacionados com o valor entre 4,51 e 4,67, matiz entre 7,27 e 7,68 e croma entre 5,20 e 5,28. Nestes locais, o teor de argila é de 16%, os teores de fósforo variam entre 65 e 75 mg dm-3, e a soma de bases, entre 56 e 58 mmolc dm-3 (Figura 6). A faixa espectral na região do visível - 370 a 750 nm - já é utilizada no detalhamento de áreas com diferentes potenciais produtivos para a cultura da cana-de-açúcar (Marques Jr. et al., 2014MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; TEIXEIRA, D.D.B.; BARRÓN, V.; TORRENT, J. Magnetic susceptibility and diffuse reflectance spectroscopy to characterize the spatial variability of soil properties in a Brazilian Haplustalf. Geoderma, v.219/220, p.63-71, 2014. DOI: 10.1016/j.geoderma.2013.12.007.
https://doi.org/10.1016/j.geoderma.2013....
).

Figura 6
Produtividade do cafeeiro "Coffea arabica 'Catuaí Amarelo IAC 62' em diferentes ambientes de Argissolo Vermelho-Amarelo com baixo teor de ferro total caracterizados pelos componentes da cor obtida por ERD.

Conclusões

  1. A cor do solo, obtida pela espectroscopia de reflectância difusa (ERD), é eficiente tanto na identificação de áreas com diferentes potenciais de produtividade quanto de qualidade do grão de café.

  2. O método indireto de caracterização da cor pela ERD pode ser uma alternativa para viabilizar futuros trabalhos, que busquem identificar áreas quanto à produtividade e qualidade de grão de café em grandes áreas.

Agradecimentos

À Fundação de Amparo à Pesquisa do Estado do São Paulo (Fapesp, processo no. 11/07549-2), pelo apoio financeiro; e ao Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), pela concessão de bolsa de produtividade.

Referências

  • BAHIA, A.S.R. de S.; MARQUES JR., J.; SIQUEIRA, D.S. Procedures using diffuse reflectance spectroscopy for estimating hematite and goethite in Oxisols of São Paulo, Brazil. Geoderma Regional, v.5, p.150-156, 2015. DOI: 10.1016/j.geodrs.2015.04.006.
  • BARRÓN, V.; MELLO, J. W. V.; TORRENT, J. Caracterização de óxidos de ferro em solos por espectroscopia de reflectância difusa. In: NOVAIS, R.F. de; ALVAREZ V., V.H.; SCHAEFER, C.E.G.R., (Ed.). Tópicos em ciência do solo Viçosa: Sociedade Brasileira de Ciência do Solo, 2000. v.1, p.139-162.
  • BARRÓN, V.; TORRENT, J. Iron, manganese and aluminium oxides and oxyhydroxides. In: NIETO, F.; LIVI, K.J.T.; OBERTI, R. (Ed.). Minerals at the nanoscale London: European Mineralogical Union: Mineralogical Society of Great Britain & Ireland, 2013. p.297-336. (European Mineralogical Union Notes in Mineralogy, v.14). DOI: 10.1180/EMU-notes.14.9.
    » https://doi.org/10.1180/EMU-notes.14.9
  • BEAUFILS, E.R. Diagnosis and recommendation integrated system (DRIS): a general scheme of experimentation and calibration based on principles developed from research in plant nutrition. Pietermaritzburg: University of Natal, 1973. 132p. (Soil Science Bulletin, 1).
  • CAMARGO, L.A.; MARQUES JÚNIOR ., J.; BARRÓN, V.; ALLEONI, L.R.F.; BARBOSA, R.S.; PEREIRA, G.T. Mapping of clay, iron oxide and adsorbed phosphate in Oxisols using diffuse reflectance spectroscopy. Geoderma, v.251/252, p.124-132, 2015. DOI: 10.1016/j.geoderma.2015.03.027.
    » https://doi.org/10.1016/j.geoderma.2015.03.027
  • CAMBARDELLA, C.A.; MOORMAN, T.B.; NOVAK, J.M.; PARKIN, T.B.; KARLEN, D.L.; TURCO, R.F.; KONOPKA, A.E. Field-scale variability of soil properties in Central Iowa soils. Soil Science Society of America Journal, v.58, p.1501-1511, 1994. DOI: 10.2136/sssaj1994.03615995005800050033x.
    » https://doi.org/10.2136/sssaj1994.03615995005800050033x
  • CAMPOS, R.C.; DEMATTÊ, J.A.M. Cor do solo: uma abordagem da forma convencional de obtenção em oposição à automatização do método para fins de classificação de solos. Revista Brasileira de Ciência do Solo, v.28, p.853-863, 2004. DOI: 10.1590/S0100-06832004000500008.
    » https://doi.org/10.1590/S0100-06832004000500008
  • CLAESSEN, M.E.C. (Org.). Manual de métodos de análise de solo 2.ed. rev. e atual. Rio de Janeiro,: EMBRAPA-CNPS, 1997. 212p.
  • DANTAS, J.S.; MARQUES JÚNIOR, J.; MARTINS FILHO, M.V.; RESENDE, J.M. do A.; CAMARGO, L.A.; BARBOSA, R.S. Gênese de solos coesos do leste maranhense: relação solo-paisagem. Revista Brasileira de Ciência do Solo, v.38, p.1039-1050, 2014a. DOI: 10.1590/S0100-06832014000400001.
    » https://doi.org/10.1590/S0100-06832014000400001
  • DANTAS, J.S.; MARTINS FILHO, M.V.; MARQUES JÚNIOR, J.; RESENDE, J.M. do A.; TEIXEIRA, D.D.B.; BARBOSA, R.S.; SIQUEIRA, D.S. Coeficiente de erodibilidade em sulcos e entressulcos de Argissolos coesos estimado pela cor do solo. Pesquisa Agropecuária Brasileira, v.49, p.700-707, 2014b. DOI: 10.1590/S0100-204X2014000900006.
    » https://doi.org/10.1590/S0100-204X2014000900006
  • DAY, P. R. Particle fraction and particle-size analysis. In: BLACK, C.A. (Ed.). Methods of soil analysis Madison: American Society of Agronomy, 1965. v.1, p.545-566. (ASA. Agronomy, 9).
  • DEUTSCH, C.V.; JOURNEL, A.G. GSLIB: Geostatistical Software Library and User's Guide. 2.ed. Oxford: Oxford University Press, 1998. 366p.
  • FERNANDES, L.A.; RIBEIRO, C.M.M. Evolution and palaeoenvironment of the Bauru Basin (Upper Cretaceous, Brazil). Journal of South American Earth Sciences, v.61, p.71-90, 2015. DOI: 10.1016/j.jsames.2014.11.007.
    » https://doi.org/10.1016/j.jsames.2014.11.007
  • IPT. INSTITUTO DE PESQUISAS TECNOLÓGICAS DO ESTADO DE SÃO PAULO. Mapa geomorfológico do Estado de São Paulo: mapas. São Paulo, 1981. 94p. (IPT Publicação, 1183; IPT Monografia, 5).
  • ISAAKS, E.H.; SRIVASTAVA, R.M. An introduction to applied geoestatistics New York: Oxford University Press, 1989. 560p.
  • MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; TEIXEIRA, D.D.B.; BARRÓN, V.; TORRENT, J. Magnetic susceptibility and diffuse reflectance spectroscopy to characterize the spatial variability of soil properties in a Brazilian Haplustalf. Geoderma, v.219/220, p.63-71, 2014. DOI: 10.1016/j.geoderma.2013.12.007.
    » https://doi.org/10.1016/j.geoderma.2013.12.007
  • MOTTA, A.C.V.; NICK, J.A.; YORINORI, G.T.; SERRAT, B.M. Distribuição horizontal e vertical da fertilidade do solo e das raízes de cafeeiro (Coffea arabica L.) cultivar Catuaí. Acta Scientiarum. Agronomy, v.28, p.455-463, 2006. DOI: 10.4025/actasciagron.v28i4.758.
    » https://doi.org/10.4025/actasciagron.v28i4.758
  • PELUCO, R.G.; MARQUES JÚNIOR, J.; SIQUEIRA, D.S.; PEREIRA, G.T; BARBOSA, R.S.; TEIXEIRA, D. de B. Mapeamento do fósforo adsorvido por meio da cor e da suscetibilidade magnética do solo. Pesquisa Agropecuária Brasileira, v.50, p.259-266, 2015. DOI: 10.1590/S0100-204X2015000300010.
    » https://doi.org/10.1590/S0100-204X2015000300010
  • POLLO, G.Z. Suscetibilidade magnética, atributos do solo e da planta na discriminação de áreas de manejo específico na cultura do café 2013. 58p. Dissertação (Mestrado) - Universidade Estadual Paulista, Jaboticabal.
  • RAIJ, B. van; QUAGGIO, J.A.; CANTARELLA, J.; FERREIRA, M.E.; LOPES, A.S.; BATAGLIA, O.C. Análise química do solo para fins de fertilidade Campinas: Fundação Cargill, 1987. 170p.
  • RESENDE, J.M. do A. Caracterização pedométrica de atributos de Argissolos coesos do leste maranhense 2013. 83p. Tese (Doutorado) - Universidade Estadual Paulista, Jaboticabal.
  • RESENDE, J.M. do A.; MARQUES JÚNIOR, J.; MARTINS FILHO, M.V.; DANTAS, J.S.; SIQUEIRA, D.S.; TEIXEIRA, D.D.B. Variabilidade espacial de atributos de solos coesos do leste Maranhense. Revista Brasileira de Ciência do Solo, v.38, p.1077-1090, 2014.
  • SANCHEZ, M.G.B.; MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; PEREIRA, G.T. Delineation of specific management areas for coffee cultivation based on the soil-relief relationship and numerical classification. Precision Agriculture, v.14, p.201-214, 2013. DOI: 10.1007/s11119-012-9288-z.
    » https://doi.org/10.1007/s11119-012-9288-z
  • SANCHEZ, R.B.; MARQUES JÚNIOR, J.; SOUZA, Z.M. de; PEREIRA, G.T.; MARTINS FILHO, M.V. Variabilidade espacial de atributos do solo e de fatores de erosão em diferentes pedoformas. Bragantia, v.68, p.1095-1103, 2009. DOI: 10.1590/S0006-87052009000400030.
    » https://doi.org/10.1590/S0006-87052009000400030
  • SANTOS, H.G. dos; JACOMINE, P.K.T.; ANJOS, L.H.C. dos; OLIVEIRA, V.A. de; OLIVEIRA, J.B. de; COELHO, M.R.; LUMBRERAS, J.F.; CUNHA, T.J.F. (Ed.). Sistema brasileiro de classificação de solos 2.ed. Brasília: Embrapa, 2006. 306p.
  • SCHEMIKO, D.C.B.; VESELY, F.F.; FERNANDES, L.A.; SOWEK, G.A. Distinção dos elementos deposicionais fluviais, eólicos e lacustres do Membro Morro Pelado, Permiano Superior da Bacia do Paraná. Geologia USP. Série Científica, v.14, p.29-46, 2014. DOI: 10.5327/Z1519-874X201400030003.
    » https://doi.org/10.5327/Z1519-874X201400030003
  • SIQUEIRA, D.S.; MARQUES JR., J.; PEREIRA, G.T. The use of landforms to predict the variability of soil and orange attributes. Geoderma, v.155, p.55-66, 2010. DOI: 10.1016/j.geoderma.2009.11.024.
    » https://doi.org/10.1016/j.geoderma.2009.11.024
  • THORNTHWAITE, C.W. An approach toward a rational classification of climate. Geographical Review, v.38, p.55-94, 1948. DOI: 10.2307/210739.
    » https://doi.org/10.2307/210739
  • TORRENT, J.; BARRÓN, V. Diffuse reflectance spectroscopy. In: ULERY, A.L.; DREES, L.R. (Ed.). Methods of soil analysis: part 5: mineralogical methods. Madison: Soil Science Society of America, 2008. p.367-387. (SSSA Book Series, 5). DOI: 10.2136/sssabookser5.5.c13.
    » https://doi.org/10.2136/sssabookser5.5.c13
  • VISCARRA ROSSEL, R.A., BUI, E.N.; CARITAT, P. de; MCKENZIE, N.J. Mapping iron oxides and the color of Australian soil using visible-near-infrared reflectance spectra. Journal of Geophysical Research, v.115, p.1-13, 2010. DOI: 10.1029/2009JF001645.
    » https://doi.org/10.1029/2009JF001645
  • WYSZECKI, G.; STILES, W.S. Color science: concepts and methods, quantitative data and formulae. 2.ed. New York: John Wiley & Sons, 1982. 968p.
  • XIE, H.T.; YANG, X.M.; DRURY, C.F.; YANG, J.Y.; ZHANG, X.D. Predicting soil organic carbon and total nitrogen using mid- and near-infrared spectra for Brookston clay loam soil in Southwestern Ontario, Canada. Canadian Journal of Soil Science, v.91, p.53-63, 2011. DOI: 10.4141/cjss10029.
    » https://doi.org/10.4141/cjss10029

Datas de Publicação

  • Publicação nesta coleção
    Set 2016

Histórico

  • Recebido
    02 Jun 2015
  • Aceito
    24 Nov 2016
Embrapa Secretaria de Pesquisa e Desenvolvimento; Pesquisa Agropecuária Brasileira Caixa Postal 040315, 70770-901 Brasília DF Brazil, Tel. +55 61 3448-1813, Fax +55 61 3340-5483 - Brasília - DF - Brazil
E-mail: pab@embrapa.br