Resumos
O presente trabalho teve por objetivo analisar as respostas do capim 'Tifton 85' (Cynodon spp.) à adubação nitrogenada, em termos de área foliar verde e de massa seca da parte aérea, utilizando índices calculados a partir de informações contidas em imagens multiespectrais, em comparação com indicadores obtidos por métodos diretos. Os tratamentos consistiram de quatro doses de N-ureia (0; 100; 150 e 200 kg ha-1), dispostos em um desenho experimental inteiramente casualizado, com quatro repetições. Aos 26 dias, após um corte de uniformização, foram determinados os seguintes atributos do dossel: o Índice de Área Foliar (IAF), a massa seca da parte aérea (MSPA), o índice de cobertura verde (ICV), utilizando uma câmara digital na faixa do visível e o Índice de Vegetação por Diferença Normalizada (NDVI), utilizando uma câmara multiespectral. Os quatro atributos do dossel vegetativo responderam positivamente à adubação nitrogenada. Enquanto o ICV mostrou relações curvilineares com o IAF e a MSPA (R² IAF=0,89; R² MSPA=0,87), o NDVI mostrou relações lineares (r² IAF=0,62; r2 MSPA=0,64), sem evidenciar tendência à saturação. Foi concluído que é possível avaliar a adubação nitrogenada em pastagens de Cynodon spp., utilizando análise de imagens, e que tanto o NDVI quanto o ICV podem ser utilizados como indicadores de produtividade do capim Tifton 85
Cynodon spp.; índice de cobertura verde; índice de área foliar; NDVI; refletância do dossel
The purpose of this study was to verify the responses of 'Tifton 85' bermudagrass (Cynodon spp.) to nitrogen fertilization in terms of green leaf area and aerial dry mass, using indices calculated from spectral information contained in multispectral images, in comparison with indicators obtained by direct methods. The bermudagrass was managed under four N treatments (0; 100; 150 and 200 kg ha-1) using a completely randomized experimental design, with four repetitions. Twenty six days after a standardization cut, four canopy traits were determined: leaf area index (LAI), aerial dry mass (ADM), green coverage index (GCI), using a RGB camera and the Normalized Difference Vegetation Index (NDVI), using a multispectral camera. All the traits responded positively to nitrogen fertilization. However, while the GCI showed curvilinear relationships with LAI (R²=0.89) and ADM (R²=0.87), the NDVI showed linear relationships with both (r² LAI=0.62; r² ADM=0.64) without evidences of saturated responses. It was concluded that it is possible to evaluate nitrogen fertilizer on pastures of Cynodon spp., using image analysis. NDVI and ICV can be used as productivity indicators of 'Tifton 85' bermudagrass
canopy reflectance; Cynodon spp.; Green Coverage Index; leaf area index; NDVI
INTRODUÇÃO:
Nos últimos anos, têm ganhado destaque os estudos sobre parametrização de atributos
morfofisiológicos dos dosséis vegetais, feitos com métodos baseados nas propriedades
ópticas das plantas (OLLINGER, 2011OLLINGER, S.V. Sources of variability in canopy reflectance and the
convergent properties of plants. New Phytologist, v.189, p.375-394, 2011. Disponível
em:
<http://onlinelibrary.wiley.com/doi/10.1111/j.1469-8137.2010.03536.x/abstract>.
Acesso em: 28 jan. 2014. doi: 10.1111/j.1469-8137.2010.03536.x.
http://onlinelibrary.wiley.com/doi/10.11...
). Tais
estudos envolvem o uso de técnicas radiométricas, que utilizam bandas espectrais
específicas, ou sistemas baseados no processamento de imagens digitais obtidas por
câmeras fotográficas (SENA JÚNIOR et al., 2008SENA JÚNIOR, D.G. et al. Discriminação entre estágios nutricionais na
cultura de trigo com técnicas de visão artificial e medidor portátil de clorofila.
Engenharia Agrícola, v.28, p.187-195, 2008. Disponível em:
<http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-69162008000100019>.
Acesso em: 5 dez. 2013. doi: 10.1590/69162008000100019.
http://www.scielo.br/scielo.php?script=s...
).
Diversos índices gerados pela combinação de valores de refletância entre bandas
espectrais específicas são utilizados como indicadores da quantidade de vegetação e
outros atributos dos dosséis, sendo o Índice de Vegetação por Diferença Normalizada
(NDVI, do inglês Normalized Difference Vegetation Index), o mais conhecido (OLLINGER, 2011OLLINGER, S.V. Sources of variability in canopy reflectance and the
convergent properties of plants. New Phytologist, v.189, p.375-394, 2011. Disponível
em:
<http://onlinelibrary.wiley.com/doi/10.1111/j.1469-8137.2010.03536.x/abstract>.
Acesso em: 28 jan. 2014. doi: 10.1111/j.1469-8137.2010.03536.x.
http://onlinelibrary.wiley.com/doi/10.11...
). Em relação às câmaras digitais,
embora com menor resolução espectral em relação aos radiômetros (SENA JÚNIOR et al., 2008SENA JÚNIOR, D.G. et al. Discriminação entre estágios nutricionais na
cultura de trigo com técnicas de visão artificial e medidor portátil de clorofila.
Engenharia Agrícola, v.28, p.187-195, 2008. Disponível em:
<http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-69162008000100019>.
Acesso em: 5 dez. 2013. doi: 10.1590/69162008000100019.
http://www.scielo.br/scielo.php?script=s...
), o seu uso tem aumentado nos últimos anos,
dada sua acessibilidade e facilidades operativas, com um mínimo de treinamento do
usuário. O processamento analítico das imagens captadas com as câmeras permite a
construção de diversos índices, relacionados com a cobertura vegetal, o acúmulo de
biomassa, e os teores de nitrogênio em diversas espécies cultivadas (SENA JÚNIOR et al. 2008SENA JÚNIOR, D.G. et al. Discriminação entre estágios nutricionais na
cultura de trigo com técnicas de visão artificial e medidor portátil de clorofila.
Engenharia Agrícola, v.28, p.187-195, 2008. Disponível em:
<http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-69162008000100019>.
Acesso em: 5 dez. 2013. doi: 10.1590/69162008000100019.
http://www.scielo.br/scielo.php?script=s...
; LI et al., 2010LI , Y. et al. Estimating the nitrogen status of crops using a digital
camera. Field Crops Research, v.118, p.221-227, 2010. Disponível em:
<http://www.sciencedirect.com/science/article/pii/S0378429010001358>. Acesso
em: 3 jan. 2014. doi: 10.1007/s13593-011-0041-1.
http://www.sciencedirect.com/science/art...
). Todavia, no caso das espécies forrageiras, a
validação dessas técnicas, dependerá das características estruturais e botânicas da
vegetação em foco, o que implica no conhecimento da variabilidade de atributos como grau
de cobertura, índice de área foliar (IAF) e acúmulo de biomassa dos dosséis (SCHUT & KETELAARS, 2003SCHUT, A.G.T.; KETELAARS, J.J.M.H. Monitoring grass swards using imaging
spectroscopy. Grass and Forage Science, v.58, p.276-286, 2003. Disponível em:
<http://www.researchgate.net/publication/229924945_Monitoring_grass_swards_sing_imaging_spectroscopy>.
Acesso em: 29 out. 2013. doi: 10.1046/j.1365-2494.2003.00379.x.
http://www.researchgate.net/publication/...
; ROTZ et al., 2008ROTZ, J. D. et al. Classification of digital photography for measuring
productive ground cover. Rangeland Ecology & Management, v.61, p.245-248, 2008.
Disponível em:
<http://www.researchgate.net/publication/225076678_Classification_of_digital_photography_for_measuring_productive_ground_cover>.
Acesso em: 24 set. 2014. doi: 10.2111/07-011.1.
http://www.researchgate.net/publication/...
).
O conhecimento da quantidade e qualidade da forragem disponível é essencial para a
gestão da produção animal em pastagens (ZHAO et al.,
2007ZHAO, D. et al. Assessment of forage biomass and quality parameters of
bermudagrass using proximal sensing of pasture canopy reflectance. Grassland Science,
v.53, p.39-49, 2007. Disponível em:
<http://onlinelibrary.wiley.com/doi/10.1111/j.1744-697X.2007.00072.x/abstract>.
Acesso em: 21 jan. 2014. doi: 10.1111/j.1744-697X.2007.00072.x.
http://onlinelibrary.wiley.com/doi/10.11...
). O acesso a essas informações tem sido feito tradicionalmente por corte
e secagem da biomassa, seguida de análises laboratoriais (MANNETJE, 2000MANNETJE, L.'t. Measuring biomass of grassland vegetation. In: MANNETJE,
L.'t e JONES, R.M. (Eds.). Field and laboratory methods for grassland and animal
production research. Wallingford: CABI Publishing, 2000. p.151-177. Disponível em:
<http://www.scielo.br/scielo.php?script=sci_nlinks&ref=000101&pid=S0100-204X200300080001100028&lng=es>.
Acesso em: 2 set. 2013. doi: 10.1079/9780851993515.0000.
http://www.scielo.br/scielo.php?script=s...
). Tais procedimentos, embora precisos, possuem limitações,
relativas a custos e tempo de execução, o que destaca a importância das metodologias
indiretas, que permitem efetuar tais avaliações em tempo real, a exemplo da utilização
de informações espectrais obtidas de câmeras digitais e multiespectrais.
O gênero Cynodon é constituído por espécies e cultivares de ampla adaptação
às condições próprias de ambientes tropicais e subtropicais (VILELA et al., 2005VILELA, D. et al. (Eds) Cynodon, forrageiras que estão revolucionando a
pecuária brasileira. Juiz de Fora: Embrapa Gado de Leite, 2005. 250p. Disponível em:
<http://www.bdpa.cnptia.embrapa.br/busca?b=ad&id=594835&biblioteca=CNPAB&busca=autoria:%22%22&qFacets=autoria:%22%22&sort=&paginacao=t&paginaAtual=120>.
Acesso em: 22 nov. 2013. doi: 10.1590/S0103-84782012005000149.
http://www.bdpa.cnptia.embrapa.br/busca?...
). Entre os diversos híbridos oriundos de
cruzamentos intra ou interespecíficos, destaca-se a cultivar 'Tifton 85' (Cynodon
spp.), por combinar elevado potencial de produção de forragem com alta
digestibilidade (QUARESMA et al., 2011QUARESMA, J.P. S. et al. Produção e composição bromatológica do
capim-tifton 85 (Cynodon spp.) submetido a doses de nitrogênio. Acta Scientiarum:
Animal Sciences, v.33, p.45-150, 2011. Disponível em:
<http://www.alice.cnptia.embrapa.br/handle/doc/901620>. Acesso em: 16 out.
2013. doi: 10.4025/actascianimsci.v33i2.9261.
http://www.alice.cnptia.embrapa.br/handl...
). O
nitrogênio (N) é essencial à expressão da produtividade do capim 'Tifton 85', tanto por
efeitos diretos do nutriente sobre o IAF e a consequente interceptação da radiação solar
(OLIVEIRA et al., 2010LI , Y. et al. Estimating the nitrogen status of crops using a digital
camera. Field Crops Research, v.118, p.221-227, 2010. Disponível em:
<http://www.sciencedirect.com/science/article/pii/S0378429010001358>. Acesso
em: 3 jan. 2014. doi: 10.1007/s13593-011-0041-1.
http://www.sciencedirect.com/science/art...
), quanto por efeitos
indiretos, ligados à absorção e assimilação fotossintética do carbono (ALDERMAN et al., 2011ALDERMAN, P.D. et al. Carbohydrate and nitrogen reserves relative to
regrowth dynamics of 'Tifton 85' bermudagrass as affected by nitrogen fertilization.
Crop Science, v.51, p.1727-1738, 2011. Disponível em:
<https://www.crops.org/publications/cs/abstracts/51/4/1727>. Acesso em: 12 dez.
2013. doi: 10.2135/cropsci2010.09.0516.
https://www.crops.org/publications/cs/ab...
; SILVA et al., 2011SILVA, R.V.M.M. et al. Uso de clorofilometro na avaliação da adubação
nitrogenada e potássica no capim Tifton 85. Revista Brasileira de Saúde e Produção
Animal, v.12, p.828-839, 2011. Disponível em:
<http://www.scielo.br/pdf/cr/2012nahead/a35512cr4895.pdf>. Acesso em: 21 jan.
2014. doi: 10.1590/S0103-84782012005000109.
http://www.scielo.br/pdf/cr/2012nahead/a...
).
Dadas as possibilidades oferecidas pelas técnicas de análise de imagens provenientes de câmeras digitais, a hipótese deste trabalho foi a de que é possível avaliar os efeitos da adubação nitrogenada no capim 'Tifton 85' através dessas ferramentas. Nesse sentido, os objetivos foram: a) verificar o comportamento de dois índices calculados a partir de imagens digitais, a saber: o Índice de Cobertura Verde (ICV), obtido com câmera fotográfica comum e o NDVI, obtido com câmera multiespectral, em relação à adubação nitrogenada; b) determinar se tais índices podem servir como indicadores de produtividade do capim 'Tifton 85'.
MATERIAL E MÉTODOS:
O experimento foi conduzido na Estação Experimental de Parasitologia "Wilhelm Otto
Neitz" do Instituto de Veterinária da UFRRJ, Seropédica, RJ (22º45'48''S; 43º41'56''W,
33 m s.n.m), em área de ocorrência de Planossolo Háplico, textura superficial arenosa
(EMBRAPA, 2006EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA (EMBRAPA). Sistema
Brasileiro de Classificação de Solos. 2.ed. Rio de Janeiro: EMBRAPA/CNPS, 2006. 306p.
Disponível em:
<http://www.agrolink.com.br/downloads/sistema-brasileiro-de-classificacao-dos-solos2006.pdf>.
Acesso em: 27 de nov. 2013.
http://www.agrolink.com.br/downloads/sis...
). O plantio do capim 'Tifton 85'
foi realizado em 2006, conforme os procedimentos detalhados por Oliveira et al. (2010OLIVEIRA, A.P.P. et al. Respostas do capim-Tifton 85 à aplicação de
nitrogênio: cobertura do solo, índice de área foliar e interceptação da radiação
solar. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.62, p.429-438, 2010.
Disponível em:
<http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0102>. Acesso em:
9 out. 2013. doi: 10.1590/09352010000200024.
http://www.scielo.br/scielo.php?script=s...
). As unidades experimentais são constituídas
por parcelas de 16 m2, sendo a parte central (2 m x 2 m) considerada a área
útil para finalidade de amostragem. No mês de dezembro de 2011 foi realizado um corte de
uniformização, a uma altura de 0,05 m do solo, seguido de uma aplicação basal de cloreto
de potássio equivalente a 160 kg K2O ha-1, em todas as unidades
experimentais.
O experimento foi desenvolvido entre os meses de janeiro e março de 2012. O delineamento foi inteiramente casualizado, com quatro repetições. Os tratamentos consistiram em quatro doses de N (0; 50; 100 e 200 kg ha-1) na forma de sulfato de amônio, aplicadas no início do mês de fevereiro, após um corte de uniformização. Em 29 de fevereiro de 2012, aos 26 dias após o corte (26 DAC), foram obtidas as imagens digitais necessárias à determinação dos índices de vegetação estudados, em conjunto com as amostragens de biomassa e IAF. Para a determinação do índice de Cobertura Verde (ICV), foram adquiridas imagens com uma câmera digital (CD) marca Sony, modelo DSC-W50, entre as 7h30min e 9h horas. Um suporte de PVC, em formato de "L", foi feito para o acoplamento ortogonal da câmera, a uma altura padrão de 1,60 m. As imagens obtidas (duas por parcela) foram transferidas para um computador e armazenadas no formato TIFF. Subsequentemente foram coletadas amostras da biomassa em pé, dentro de um quadrante com 0,04 m2, por corte rente ao solo, as quais foram mantidas em caixa isotérmica, à sombra. A coleta de imagens para a estimativa da refletância do dossel foi realizada entre as 11 e 13 horas, sob condições de céu claro. Utilizou-se uma câmera multiespectral (CME) marca FluxData, modelo FD-1665-MS 5 (Channel Camera, Flux Data, Inc., NY), configurada nas seguintes bandas: azul (400-510nm); verde (500-580nm); vermelho (585-640nm), "red-edge" (720-730nm) e infravermelho próximo (770- 830nm). A CME foi posicionada a uma altura de 3,0 m e direcionada a nadir, utilizando uma estrutura metálica, sendo conectada a um computador, para permitir a visualização das imagens no momento da sua obtenção. Em cada parcela, foram tomadas três imagens, com resolução espacial de 0,33 mm/pixel. No campo de visão da imagem, era posicionado um painel padrão Spectralon (Labsphere, Inc., North Sutton, NH, USA), contendo quatro tons de cinza: branco difuso, cinza claro, cinza médio e cinza escuro, com refletâncias de 99%, 75%, 40% e 10%, respectivamente, para permitir a conversão dos valores de pixel em valores de refletância. O nível de refletância do dossel, dentro da faixa 400-700 nm, foi também avaliado com dois sensores de quantum (LI-190) conectados a radiômetros portáteis (marca LI-COR, modelo LI-250A) (LI-COR, Inc., Lincoln, NE), posicionados a uma altura de 0,5 m acima da superfície do relvado e apontados para cima e para abaixo, consecutivamente.
As imagens obtidas com CD foram processadas e analisadas com auxílio do programa MATLAB
v. 9 (MATHWORKS, 2000MATHWORKS. Image Processing Toolbox for use with MATLAB Release 12.
Image Processing Toolbox User's Guide Revised for Version 2.2.2, Massachusetts, 2000.
511p.). Na etapa de
pré-processamento, foram escolhidos pontos representativos de pixels de vegetação verde
ou de solo nu ou com vegetação seca, para criar um arquivo com o vetor classificador,
baseado na função discriminante de Fisher (JOHNSON &
WICHERN, 1999JOHNSON, R.A.; WICHERN, D.W. Applied multivariate statistical analysis.
4 ed. Upper Saddle River, New Jersey: Prentice-Hall, 1999. 815p. Disponível em:
<http://esminfo.prenhall.com/math/johnsonwichern/data.html>. Acesso em: 2 set.
2013.
http://esminfo.prenhall.com/math/johnson...
). Este procedimento permitiu a binarização da imagem em branco e
preto, em que os pontos pretos representaram a vegetação verde. As imagens foram, então,
recortadas de '2816C x 2112L' para '2040C x 2040L', correspondente a um campo de visão
de 1m² e resolução espacial de 0,49 mm/pixel, e processadas para a obtenção do ICV, que
corresponde ao porcentual de pixels verdes em relação ao número total de pixels da
imagem. As imagens de pixels obtidas com CME foram convertidas em imagens de
refletância, baseadas em cinco equações de calibração, uma para cada banda espectral,
resultantes de regressão linear simples, entre os valores de pixels das imagens
referentes a cada tom de cinza do painel Spectralon e os respectivos valores de
refletância fornecidos pelo fabricante, de acordo com a equação (1):
(1)
em que: Pijk = valor do pixel na linha i e coluna j da
banda k, adimensional; Ŷijk = refletância na linha i e coluna j da banda
k, em %; a e b = coeficientes de regressão. Procedimento similar foi
descrito por Abrahão et al. (2013ABRAHÃO, S.A. et al. Determination of nitrogen and chlorophyll levels in
bean-plant leaves by using spectral vegetation bands and indices. Revista Ciência
Agronômica, v.44, p.464-473, 2013. Disponível em:
<http://www.scielo.br/scielo.php?pid=S1806-66902013000300007&script=sci_arttext>.
Acesso em: 25 fev. 2014. doi: 10.1590/S1806-66902013000300007.
http://www.scielo.br/scielo.php?pid=S180...
). O
processamento das imagens foi realizado com auxílio do programa MATLAB v. 9 (MATHWORKS, 2000MATHWORKS. Image Processing Toolbox for use with MATLAB Release 12.
Image Processing Toolbox User's Guide Revised for Version 2.2.2, Massachusetts, 2000.
511p.). No presente estudo, foi utilizado o
NDVI, definido como (2):
(2)
em que: IPV = valor de refletância na banda do infravermelho próximo; V = valor de refletância na banda do vermelho.
A biomassa colhida no campo foi imediatamente fracionada no laboratório, separando-se o material verde do material morto ou senescente. O material verde foi separado em lâminas foliares e colmo+bainhas. Sub-amostras, contendo de 10 a 20 % da massa fresca de lâminas foliares, foram separadas e digitalizadas diretamente em scanner, com uma resolução de 200 dpi. As imagens obtidas foram arquivadas no formato JPEG e analisadas com auxílio do software SIARCS(r), para a determinação da sua área superficial. Subsequentemente, o material foi posto a secar em estufa com circulação forçada de ar (65ºC), por 48 h. Após pesagem, os dados foram utilizados para computar a massa seca de lâminas foliares (MSLF, kg m -2), a área foliar específica (AFE, m² folha kg MS folha -1), o IAF (calculado como AFE x MSLF) e a massa seca da parte aérea (MSPA, kg m-2) do capim.
Os dados meteorológicos relativos ao período experimental foram obtidos junto ao Posto
da Estação Experimental de Seropédica, da PESAGRO-RJ. Durante o mês de fevereiro de
2012, a precipitação pluviométrica total foi de 42,6 mm, ou seja, um terço da normal
climática regional (MATTOS et al., 1998MATTOS, C.C.L.V. et al. Boletim Agrometeorológico UFRRJ 1996. Floresta e
Ambiente, v.5, p.208-215, 1998. Disponível em:
<http://www.geocities.ws/floramrural/0208.pdf>. Acesso em: 5 set. 2013. doi:
10.4322/floram.2011.058.
http://www.geocities.ws/floramrural/0208...
). Houve
elevação das temperaturas médias do ar, com mínima de 21,1º e máxima de 31,3º. Esse
quadro ambiental prejudicou a velocidade de rebrota do capim, sendo atenuado por regas,
durante os dias de maior demanda evapotranspirativa.
Os dados foram submetidos à análise estatística para verificar a significância do fator dose de N (teste F, P≤0,05). Em atenção à natureza quantitativa da variável independente, usaram-se técnicas de regressão para a parametrização das respostas à adubação nitrogenada. Da mesma forma, em consonância com os objetivos deste estudo, as interrelações entre as variáveis experimentais, medidas ou estimadas, foram analisadas por regressão, com auxílio do software GraphPad Prism v. 5.0. Quando foi necessário optar entre modelos de regressão, utilizou-se o Critério de Informação de Akaike, já incluído no software, para a escolha do modelo mais provável.
RESULTADOS E DISCUSSÃO:
Aos 26 DAC, os dosséis já tinham acumulado uma quantidade apreciável de cobertura verde,
de forma que os níveis de radiação fotossinteticamente ativa interceptada (400-700nm)
oscilaram entre 35,8±3,2 (testemunha) e 71,1±2,8 % (200 kg N ha-1). Esses
valores foram menores que os observados na mesma área, em anos anteriores (OLIVEIRA et al., 2010LI , Y. et al. Estimating the nitrogen status of crops using a digital
camera. Field Crops Research, v.118, p.221-227, 2010. Disponível em:
<http://www.sciencedirect.com/science/article/pii/S0378429010001358>. Acesso
em: 3 jan. 2014. doi: 10.1007/s13593-011-0041-1.
http://www.sciencedirect.com/science/art...
), possivelmente em razão da
baixa precipitação pluviométrica no ciclo de rebrota estudado. As variáveis estudadas
(IAF, ICV, MSPA e NDVI) responderam positivamente à adubação nitrogenada. Tanto o IAF
(Figura 1A) quanto o ICV (Figura 1B), guardaram uma relação linear com as doses de N aplicadas.
O alto grau de associação entre essas variáveis deveu-se, em parte, à própria
metodologia de determinação, a qual excluiu a área foliar senescente ou morta,
manualmente, no caso do IAF, ou através da aplicação de algoritmos específicos, no caso
da classificação das imagens.
A) Índice de área foliar (IAF), B) Índice de Cobertura Verde (ICV), C) Índice de Vegetação por Diferença Normalizada (NDVI) e D)Matéria seca da parte aérea (MSPA) de dosséis do capim 'Tifton 85' em função de doses de nitrogênio aplicadas. As linhas pontilhadas representam o intervalo de confiança da regressão, ao nível de 95%. Cada ponto é uma média de quatro repetições. Seropédica, RJ, 29 de fevereiro de 2012.
Os dados coletados com os sensores quânticos indicaram que a refletância média do
dossel, dentro da faixa fotossinteticamente ativa (400-700nm), decresceu
significativamente (P<0,05), de 8,4 para 5,3 %, entre as doses de zero e 200kg N
há-1, respectivamente. A maior absorção de RFA pela clorofila a
está associada aos comprimentos de onda na faixa do vermelho (ABRAHÃO et al. , 2013ABRAHÃO, S.A. et al. Determination of nitrogen and chlorophyll levels in
bean-plant leaves by using spectral vegetation bands and indices. Revista Ciência
Agronômica, v.44, p.464-473, 2013. Disponível em:
<http://www.scielo.br/scielo.php?pid=S1806-66902013000300007&script=sci_arttext>.
Acesso em: 25 fev. 2014. doi: 10.1590/S1806-66902013000300007.
http://www.scielo.br/scielo.php?pid=S180...
), justificando as reduções da
refletância observadas dentro da faixa visível do espectro (XIONG et al., 2007XIONG, X. et al. Bermudagrass seasonal responses to nitrogen
fertilization and irrigation detected using optical sensing., Crop Science v.47,
p.1603-1610, 2007. Disponível em:
<http://agris.fao.org/agris-search/search.do?recordID=US201300831667>. Acesso
em: 17 nov. 2013. doi: 10.2135/cropsci2006.06.0400.
http://agris.fao.org/agris-search/search...
). Por outro lado, a refletância na banda do
infravermelho aumenta de forma contínua com os aumentos em área foliar ou MSPA (OLLINGER, 2011OLLINGER, S.V. Sources of variability in canopy reflectance and the
convergent properties of plants. New Phytologist, v.189, p.375-394, 2011. Disponível
em:
<http://onlinelibrary.wiley.com/doi/10.1111/j.1469-8137.2010.03536.x/abstract>.
Acesso em: 28 jan. 2014. doi: 10.1111/j.1469-8137.2010.03536.x.
http://onlinelibrary.wiley.com/doi/10.11...
). A consequência desse comportamento
espectral da vegetação é o aumento do NDVI com a dose de N aplicada (Figura 1C). Todavia, o modelo linear obtido entre NDVI
e dose de N teve um coeficiente de determinação menor (r2=0,63) que os
correspondentes ao IAF e ao ICV (Figura 1A e B),
evidenciando um aumento da variabilidade dentro de tratamentos. Parte dessa
variabilidade pode ser atribuída a flutuações nos teores de N disponível, entre
repetições, as quais afetam componentes importantes da geometria do dossel no capim
'Tifton 85', como altura e distribuição angular de folhas e estolões (ROSSIELLO & ANTUNES, 2012ROSSIELLO, R.O.P.; ANTUNES, M.A.H. Solar radiation utilization by
tropical forage grasses: light interception and use efficiency. In: BABATUNDE, E.B.
(Ed.) Solar radiation. Rijeka: In Tech, 2012. p.221-244. Disponível em:
<http://www.intechopen.com/books/solar-radiation/solar-radiation-utilization-by-tropical-forage-grasses-light-interception-and-use-efficienc>.
Acesso em: 7 fev. 2013. doi: 10.5772/34321.
http://www.intechopen.com/books/solar-ra...
). Em relação à
resposta da produção da MSPA à adubação nitrogenada, ela evidenciou um padrão de
resposta quadrático (Figura 1D), similar ao
observado em experimentos anteriores (SILVA et al.,
2011SILVA, R.V.M.M. et al. Uso de clorofilometro na avaliação da adubação
nitrogenada e potássica no capim Tifton 85. Revista Brasileira de Saúde e Produção
Animal, v.12, p.828-839, 2011. Disponível em:
<http://www.scielo.br/pdf/cr/2012nahead/a35512cr4895.pdf>. Acesso em: 21 jan.
2014. doi: 10.1590/S0103-84782012005000109.
http://www.scielo.br/pdf/cr/2012nahead/a...
).
Foi verificado que o ICV relacionou-se de forma curvilinear tanto com o IAF (Figura 2A) quanto com a MSPA (Figura 2B). No caso da relação com o IAF, o modelo de melhor ajuste
foi o sigmoidal, que indicou, para valores de IAF > 3,0, uma tendência crescente à
estabilização do ICV (Figura 2A). Isso significa
que incrementos adicionais de IAF implicarão em incrementos progressivamente menores do
ICV. A implicação é que um atributo como o ICV deva guardar uma relação assintótica com
o IAF, tal como acontece com a fração da RFA incidente interceptada pelo dossel (OLIVEIRA et al., 2010LI , Y. et al. Estimating the nitrogen status of crops using a digital
camera. Field Crops Research, v.118, p.221-227, 2010. Disponível em:
<http://www.sciencedirect.com/science/article/pii/S0378429010001358>. Acesso
em: 3 jan. 2014. doi: 10.1007/s13593-011-0041-1.
http://www.sciencedirect.com/science/art...
). Já com relação à variação de
MSPA induzida pela adubação nitrogenada, o ICV relacionou-se de forma exponencial (Figura 2B) sem evidenciar tendência à saturação, em
função de valores de biomassa não excessivamente elevados. Na figura 2, também são
apresentadas as relações entre o NDVI e o IAF (Figura
2C) ou a MSPA (Figura 2D), ambas de
natureza linear, mas com coeficientes de determinação menores
(0,62<r2<0,65) do que os indicados na figura 2A e 2B. Esse aumento na
variabilidade experimental pode ser consequência de que as estimativas de refletância
foram baseadas na estrutura intacta do dossel, enquanto as de área foliar ou biomassa
tomaram por base valores obtidos de forma destrutiva. Entretanto, subsiste o fato de
que, com a utilização da informação proveniente da banda do infravermelho próximo, as
relações entre os atributos do dossel mostraram um comportamento linear, diferentemente
das obtidas com a CD, na banda visível do espectro solar.
Relações entre o Índice de Cobertura Verde (ICV) e o Índice de Vegetação por Diferença Normalizada (NDVI), com valores de Índice de Área Foliar (IAF) (A e C) ou com a Massa seca da Parte Aérea (MSPA) (B e D). Cada ponto representa uma dose de adubação nitrogenada e é a média de quatro repetições. Barras verticais ou horizontais representam os erros padrão das estimativas para cada variável. Sy.x indica o desvio padrão da regressão. Seropédica, RJ, 29 de fevereiro de 2012.
Não foram observadas as tendências de saturação relatadas para outras espécies vegetais
com o uso do NDVI (OLLINGER, 2011OLLINGER, S.V. Sources of variability in canopy reflectance and the
convergent properties of plants. New Phytologist, v.189, p.375-394, 2011. Disponível
em:
<http://onlinelibrary.wiley.com/doi/10.1111/j.1469-8137.2010.03536.x/abstract>.
Acesso em: 28 jan. 2014. doi: 10.1111/j.1469-8137.2010.03536.x.
http://onlinelibrary.wiley.com/doi/10.11...
). Então, uma
vantagem associada à determinação "in situ" do NDVI seria a atenuação das tendências de
resposta não lineares observadas com o uso de câmeras digitais comuns ou com outros
métodos manuais de quantificação da cobertura verde, como a frequência de ocorrência
(OLIVEIRA et al., 2010LI , Y. et al. Estimating the nitrogen status of crops using a digital
camera. Field Crops Research, v.118, p.221-227, 2010. Disponível em:
<http://www.sciencedirect.com/science/article/pii/S0378429010001358>. Acesso
em: 3 jan. 2014. doi: 10.1007/s13593-011-0041-1.
http://www.sciencedirect.com/science/art...
). Isto pode ser uma
propriedade interessante, pensando na possibilidade de uso do NDVI como um estimador da
produção de matéria seca do capim, frente às diferentes doses de adubação nitrogenada.
SCHUT & KETELAARS (2003SCHUT, A.G.T.; KETELAARS, J.J.M.H. Monitoring grass swards using imaging
spectroscopy. Grass and Forage Science, v.58, p.276-286, 2003. Disponível em:
<http://www.researchgate.net/publication/229924945_Monitoring_grass_swards_sing_imaging_spectroscopy>.
Acesso em: 29 out. 2013. doi: 10.1046/j.1365-2494.2003.00379.x.
http://www.researchgate.net/publication/...
) observaram, em
mini-dosséis de Lolium perenne, que um índice obtido a partir de imagens
classificadas da cobertura vegetal foi um estimador da interceptação luminosa melhor do
que o próprio IAF. Esse resultado sugere que os índices de vegetação poderiam ser
utilizados como indicadores do IAF, dependendo da finalidade com o qual este último é
determinado. Em manejo de pastagens, o conceito de IAF "crítico" corresponde à condição
em que o dossel intercepta 95% da luz incidente (BROUGHAM, 1958BROUGHAM, R.W. Interception of light by the foliage of pure and mixed
stands of pasture plants. Australian Journal of Agricultural Research, v.9, p.39-52,
1958. Disponível em:
<http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-204X2001000100023>.
Acesso em: 21 jan. 2014. doi: 10.1071/AR9580039.
http://www.scielo.br/scielo.php?script=s...
), estando associado à maximização da taxa de acúmulo líquido
de forragem, e, portanto, ao momento mais adequado para a interrupção da rebrota, por
corte ou pastejo. Seria interessante, portanto, comparar o desempenho do ICV e do NDVI
com os dos analisadores ópticos do dossel, em valores de IAF mais próximos do crítico,
que, para o 'Tifton 85', poderia se situar entre 4,0 (OLIVEIRA et al., 2010LI , Y. et al. Estimating the nitrogen status of crops using a digital
camera. Field Crops Research, v.118, p.221-227, 2010. Disponível em:
<http://www.sciencedirect.com/science/article/pii/S0378429010001358>. Acesso
em: 3 jan. 2014. doi: 10.1007/s13593-011-0041-1.
http://www.sciencedirect.com/science/art...
) e 4,9 (OLIVEIRA et al.,
2000OLIVEIRA, M.A. et al. Análise de crescimento do capim bermuda "Tifton
85" (Cynodon spp.). Revista Brasileira de Zootecnia, v.29, p.1930-1938, 2000.
Disponível em:
<http://www.revistasbz.org.br/scripts/revista/sbz1/Artigos/2825.pdf>. Acesso
em: 3 nov. 2013.
http://www.revistasbz.org.br/scripts/rev...
), dependendo da duração do período de rebrota. XIONG et al. (2007XIONG, X. et al. Bermudagrass seasonal responses to nitrogen
fertilization and irrigation detected using optical sensing., Crop Science v.47,
p.1603-1610, 2007. Disponível em:
<http://agris.fao.org/agris-search/search.do?recordID=US201300831667>. Acesso
em: 17 nov. 2013. doi: 10.2135/cropsci2006.06.0400.
http://agris.fao.org/agris-search/search...
) comparou a efetividade de quatro índices ópticos
e concluiu que o NDVI foi o melhor indicador dos efeitos da irrigação, sazonalidade e
adubação nitrogenada sobre a qualidade de gramados de Cynodon dactylon L.
Esse levantamento das potencialidades do NDVI em relação às gramas bermuda suporta os
resultados indicados nas figuras 1C, 2C e 2D, se bem que, por serem produtos de uma
única data de amostragem, devam ser interpretados com cautela.
CONCLUSÃO:
É possível avaliar a adubação nitrogenada em pastagens de Cynodon spp., utilizando análise de imagens dos dosséis. Tanto o NDVI quanto o ICV podem ser utilizados como indicadores de produtividade do capim 'Tifton 85'.
- ABRAHÃO, S.A. et al. Determination of nitrogen and chlorophyll levels in bean-plant leaves by using spectral vegetation bands and indices. Revista Ciência Agronômica, v.44, p.464-473, 2013. Disponível em: <http://www.scielo.br/scielo.php?pid=S1806-66902013000300007&script=sci_arttext>. Acesso em: 25 fev. 2014. doi: 10.1590/S1806-66902013000300007.
» https://doi.org/10.1590/S1806-66902013000300007» http://www.scielo.br/scielo.php?pid=S1806-66902013000300007&script=sci_arttext - ALDERMAN, P.D. et al. Carbohydrate and nitrogen reserves relative to regrowth dynamics of 'Tifton 85' bermudagrass as affected by nitrogen fertilization. Crop Science, v.51, p.1727-1738, 2011. Disponível em: <https://www.crops.org/publications/cs/abstracts/51/4/1727>. Acesso em: 12 dez. 2013. doi: 10.2135/cropsci2010.09.0516.
» https://doi.org/10.2135/cropsci2010.09.0516» https://www.crops.org/publications/cs/abstracts/51/4/1727 - BROUGHAM, R.W. Interception of light by the foliage of pure and mixed stands of pasture plants. Australian Journal of Agricultural Research, v.9, p.39-52, 1958. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-204X2001000100023>. Acesso em: 21 jan. 2014. doi: 10.1071/AR9580039.
» https://doi.org/10.1071/AR9580039» http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-204X2001000100023 - EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA (EMBRAPA). Sistema Brasileiro de Classificação de Solos. 2.ed. Rio de Janeiro: EMBRAPA/CNPS, 2006. 306p. Disponível em: <http://www.agrolink.com.br/downloads/sistema-brasileiro-de-classificacao-dos-solos2006.pdf>. Acesso em: 27 de nov. 2013.
» http://www.agrolink.com.br/downloads/sistema-brasileiro-de-classificacao-dos-solos2006.pdf - JOHNSON, R.A.; WICHERN, D.W. Applied multivariate statistical analysis. 4 ed. Upper Saddle River, New Jersey: Prentice-Hall, 1999. 815p. Disponível em: <http://esminfo.prenhall.com/math/johnsonwichern/data.html>. Acesso em: 2 set. 2013.
» http://esminfo.prenhall.com/math/johnsonwichern/data.html - LI , Y. et al. Estimating the nitrogen status of crops using a digital camera. Field Crops Research, v.118, p.221-227, 2010. Disponível em: <http://www.sciencedirect.com/science/article/pii/S0378429010001358>. Acesso em: 3 jan. 2014. doi: 10.1007/s13593-011-0041-1.
» https://doi.org/10.1007/s13593-011-0041-1» http://www.sciencedirect.com/science/article/pii/S0378429010001358 - MANNETJE, L.'t. Measuring biomass of grassland vegetation. In: MANNETJE, L.'t e JONES, R.M. (Eds.). Field and laboratory methods for grassland and animal production research. Wallingford: CABI Publishing, 2000. p.151-177. Disponível em: <http://www.scielo.br/scielo.php?script=sci_nlinks&ref=000101&pid=S0100-204X200300080001100028&lng=es>. Acesso em: 2 set. 2013. doi: 10.1079/9780851993515.0000.
» https://doi.org/10.1079/9780851993515.0000» http://www.scielo.br/scielo.php?script=sci_nlinks&ref=000101&pid=S0100-204X200300080001100028&lng=es - MATTOS, C.C.L.V. et al. Boletim Agrometeorológico UFRRJ 1996. Floresta e Ambiente, v.5, p.208-215, 1998. Disponível em: <http://www.geocities.ws/floramrural/0208.pdf>. Acesso em: 5 set. 2013. doi: 10.4322/floram.2011.058.
» https://doi.org/10.4322/floram.2011.058» http://www.geocities.ws/floramrural/0208.pdf - MATHWORKS. Image Processing Toolbox for use with MATLAB Release 12. Image Processing Toolbox User's Guide Revised for Version 2.2.2, Massachusetts, 2000. 511p.
- OLIVEIRA, M.A. et al. Análise de crescimento do capim bermuda "Tifton 85" (Cynodon spp.). Revista Brasileira de Zootecnia, v.29, p.1930-1938, 2000. Disponível em: <http://www.revistasbz.org.br/scripts/revista/sbz1/Artigos/2825.pdf>. Acesso em: 3 nov. 2013.
» http://www.revistasbz.org.br/scripts/revista/sbz1/Artigos/2825.pdf - OLIVEIRA, A.P.P. et al. Respostas do capim-Tifton 85 à aplicação de nitrogênio: cobertura do solo, índice de área foliar e interceptação da radiação solar. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.62, p.429-438, 2010. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0102>. Acesso em: 9 out. 2013. doi: 10.1590/09352010000200024.
» https://doi.org/10.1590/09352010000200024» http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0102 - OLLINGER, S.V. Sources of variability in canopy reflectance and the convergent properties of plants. New Phytologist, v.189, p.375-394, 2011. Disponível em: <http://onlinelibrary.wiley.com/doi/10.1111/j.1469-8137.2010.03536.x/abstract>. Acesso em: 28 jan. 2014. doi: 10.1111/j.1469-8137.2010.03536.x.
» https://doi.org/10.1111/j.1469-8137.2010.03536.x» http://onlinelibrary.wiley.com/doi/10.1111/j.1469-8137.2010.03536.x/abstract - QUARESMA, J.P. S. et al. Produção e composição bromatológica do capim-tifton 85 (Cynodon spp.) submetido a doses de nitrogênio. Acta Scientiarum: Animal Sciences, v.33, p.45-150, 2011. Disponível em: <http://www.alice.cnptia.embrapa.br/handle/doc/901620>. Acesso em: 16 out. 2013. doi: 10.4025/actascianimsci.v33i2.9261.
» https://doi.org/10.4025/actascianimsci.v33i2.9261» http://www.alice.cnptia.embrapa.br/handle/doc/901620 - ROSSIELLO, R.O.P.; ANTUNES, M.A.H. Solar radiation utilization by tropical forage grasses: light interception and use efficiency. In: BABATUNDE, E.B. (Ed.) Solar radiation. Rijeka: In Tech, 2012. p.221-244. Disponível em: <http://www.intechopen.com/books/solar-radiation/solar-radiation-utilization-by-tropical-forage-grasses-light-interception-and-use-efficienc>. Acesso em: 7 fev. 2013. doi: 10.5772/34321.
» https://doi.org/10.5772/34321» http://www.intechopen.com/books/solar-radiation/solar-radiation-utilization-by-tropical-forage-grasses-light-interception-and-use-efficienc - ROTZ, J. D. et al. Classification of digital photography for measuring productive ground cover. Rangeland Ecology & Management, v.61, p.245-248, 2008. Disponível em: <http://www.researchgate.net/publication/225076678_Classification_of_digital_photography_for_measuring_productive_ground_cover>. Acesso em: 24 set. 2014. doi: 10.2111/07-011.1.
» https://doi.org/10.2111/07-011.1» http://www.researchgate.net/publication/225076678_Classification_of_digital_photography_for_measuring_productive_ground_cover - SCHUT, A.G.T.; KETELAARS, J.J.M.H. Monitoring grass swards using imaging spectroscopy. Grass and Forage Science, v.58, p.276-286, 2003. Disponível em: <http://www.researchgate.net/publication/229924945_Monitoring_grass_swards_sing_imaging_spectroscopy>. Acesso em: 29 out. 2013. doi: 10.1046/j.1365-2494.2003.00379.x.
» https://doi.org/10.1046/j.1365-2494.2003.00379.x» http://www.researchgate.net/publication/229924945_Monitoring_grass_swards_sing_imaging_spectroscopy - SENA JÚNIOR, D.G. et al. Discriminação entre estágios nutricionais na cultura de trigo com técnicas de visão artificial e medidor portátil de clorofila. Engenharia Agrícola, v.28, p.187-195, 2008. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-69162008000100019>. Acesso em: 5 dez. 2013. doi: 10.1590/69162008000100019.
» https://doi.org/10.1590/69162008000100019» http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-69162008000100019 - SILVA, R.V.M.M. et al. Uso de clorofilometro na avaliação da adubação nitrogenada e potássica no capim Tifton 85. Revista Brasileira de Saúde e Produção Animal, v.12, p.828-839, 2011. Disponível em: <http://www.scielo.br/pdf/cr/2012nahead/a35512cr4895.pdf>. Acesso em: 21 jan. 2014. doi: 10.1590/S0103-84782012005000109.
» https://doi.org/10.1590/S0103-84782012005000109» http://www.scielo.br/pdf/cr/2012nahead/a35512cr4895.pdf - VILELA, D. et al. (Eds) Cynodon, forrageiras que estão revolucionando a pecuária brasileira. Juiz de Fora: Embrapa Gado de Leite, 2005. 250p. Disponível em: <http://www.bdpa.cnptia.embrapa.br/busca?b=ad&id=594835&biblioteca=CNPAB&busca=autoria:%22%22&qFacets=autoria:%22%22&sort=&paginacao=t&paginaAtual=120>. Acesso em: 22 nov. 2013. doi: 10.1590/S0103-84782012005000149.
» https://doi.org/10.1590/S0103-84782012005000149» http://www.bdpa.cnptia.embrapa.br/busca?b=ad&id=594835&biblioteca=CNPAB&busca=autoria:%22%22&qFacets=autoria:%22%22&sort=&paginacao=t&paginaAtual=120 - XIONG, X. et al. Bermudagrass seasonal responses to nitrogen fertilization and irrigation detected using optical sensing., Crop Science v.47, p.1603-1610, 2007. Disponível em: <http://agris.fao.org/agris-search/search.do?recordID=US201300831667>. Acesso em: 17 nov. 2013. doi: 10.2135/cropsci2006.06.0400.
» https://doi.org/10.2135/cropsci2006.06.0400» http://agris.fao.org/agris-search/search.do?recordID=US201300831667 - ZHAO, D. et al. Assessment of forage biomass and quality parameters of bermudagrass using proximal sensing of pasture canopy reflectance. Grassland Science, v.53, p.39-49, 2007. Disponível em: <http://onlinelibrary.wiley.com/doi/10.1111/j.1744-697X.2007.00072.x/abstract>. Acesso em: 21 jan. 2014. doi: 10.1111/j.1744-697X.2007.00072.x.
» https://doi.org/10.1111/j.1744-697X.2007.00072.x» http://onlinelibrary.wiley.com/doi/10.1111/j.1744-697X.2007.00072.x/abstract
Datas de Publicação
-
Publicação nesta coleção
Abr 2015
Histórico
-
Recebido
25 Out 2013 -
Aceito
24 Jul 2014