Caracterização estrutural e florística da área de estudo

Um levantamento florístico e estrutural foi efetuado sobre indivíduos com mais de 10 cm de diâmetro à altura do peito (DAP). O trabalho de campo foi realizado entre os dias 11 e 14 de setembro de 2013, consistindo no levantamento dos indivíduos contidos em 7 parcelas amostrais de 20 m por 100 m, distribuídas ao longo da estrada que corta o PET (Figura 1). Foram coletadas as seguintes informações: espécie, DAP, altura total e posição na parcela. Esses dados foram utilizados como auxílio para a interpretação das informações espectrais e dos índices de vegetação.

Aquisição e pré-processamento dos produtos de sensoriamento remoto

Foram adquiridas duas cenas do sensor RapidEye: 28 de junho e 17 de outubro de 2012. Estas datas correspondem aos períodos de baixo e alto índice de área foliar (IAF), como também de alto e baixo valor de AZS, respectivamente. Todas as imagens foram obtidas próximo a visada do nadir e adquiridas com correção geométrica e ortorretificação (produto RapidEye nível 3A). Os valores de número digital dos pixels foram transformados em valores de reflectância de superfície por meio do modelo Fast Line-of-sight Atmospheric Analysis of Hypercubes (FLAASH) (^{EXELIS, 2014}), que permite a correção atmosférica baseada em modelos físicos. Para a correção atmosférica foi escolhido um modelo atmosférico subtropical com modelo de aerossol rural e visibilidade padrão (50 km). O modelo digital de elevação Advanced Spaceborne Thermal Emission and Reflection Radiometer - Global Digital Elevation Map version 2 (ASTER GDEM 2) foi utilizado para a modelagem topográfica. Não foi aplicada a correção do efeito da topografia sobre os dados de reflectância de superfície e índices de vegetação.

Modelagem topográfica e índices de vegetação

A partir do modelo ASTER GDEM 2 foi elaborado o modelo de relevo sombreado considerando as características locais do relevo e os parâmetros da geometria de iluminação das cenas RapidEye. Para o mês de junho, utilizou-se o valor do ângulo azimutal solar de 16,37° e o do ângulo de elevação solar de 37,93°. Para o mês de outubro, o valor do ângulo azimutal solar utilizado foi de 22,31° e o ângulo de elevação solar foi de 71,25°. O relevo sombreado foi classificado empiricamente em sombreado, intermediário e iluminado para os seguintes intervalos de valores de cosseno do ângulo de incidência i: 0 - 0,78; 0,78 - 0,87 e 0,87 - 1,0, respectivamente. Dada a alta correlação entre o relevo sombreado e o fator de iluminação para a avaliação dos dados de reflectância e índices de vegetação, o primeiro foi selecionado. A altimetria relativa foi extraída diretamente do modelo ASTER GDEM 2 por meio do fatiamento da amplitude de elevação do PET em 9 classes com intervalos de 40 metros.

Os índices de vegetação EVI e NDVI foram obtidos por meio da reflectância das bandas do sensor RapidEye. O EVI foi obtido com o uso da Equação 1, apresentada por ^{Huete et al. (2002)}:

Onde: G corresponde ao fator de ganho, o qual se atribui o valor de 2,5; C₁ e C₂ são coeficientes de correção do espalhamento atmosférico causado pelos aerossóis e assumem os valores de 6,0 e 7,5, respectivamente; L é o fator de ajustamento para correção do efeito de solos, com valor de 1,0.

O NDVI foi obtido com o uso da Equação 2, apresentada por ^{Rouse et al. (1973)}:

Onde: ρ_NIR representa a banda correspondente à região do infravermelho próximo do espectro eletromagnético; p_red representa a banda correspondente à região do vermelho do espectro eletromagnético.

Análise dos dados espectrais em função da modelagem topográfica

A análise da influencia do efeito topográfico sobre a reflectância e índices de vegetação foi feita a partir da leitura de pixels contendo valores de reflectância e de índices de vegetação, dentro de cada classe da modelagem topográfica. Devido à grande variação de tamanho entre as classes da modelagem topográfica, foi realizada uma amostragem aleatória de 1000 pixels para cada classe. Os dados resultantes foram organizados em planilhas e analisados estatisticamente e graficamente.

Afim de permitir a comparação entre as bandas, foi realizada a normalização dos valores amostrados em função de cada classe das variáveis topográficas. Os valores amostrados em função da orientação de vertentes foram normalizados a partir dos valores obtidos para a classe de orientação sul (a priori mais baixos). Para a altimetria foi escolhida a classe alti-métrica mais baixa, com altitudes entre 120 e 160 metros. Já os valores referentes ao relevo sombreado foram amostrados em função da classe de condição intermediária de iluminação. A análise dos dados foi feita por meio de gráficos.

Para avaliar de forma direta as variações decorrentes da influencia das condições de iluminação, foram elaborados dois gráficos tridimensionais a partir de dois transectos orientados nas direções Norte-Sul e Leste-Oeste. Os gráficos permitiram avaliar de forma mais eficaz o comportamento da reflectância e dos índices de vegetação sobre diferentes condições de iluminação determinadas pelo relevo.

Análise dos dados do inventário florestal

Foram identificadas 31 famílias botânicas, 65 gêneros e 74 espécies nas 14 parcelas estudadas. Esses resultados concordam com os reportados na literatura para a mesma área de estudo (^{Ruschel et al., 2005}; ^{Ruschel et al., 2007}; ^{Breunig et al., 2015}). A maioria das espécies pertence ao estrato de arvoretas, o que compreende 77,13% das espécies, enquanto 21,26% das espécies pertencem ao estrato superior e apenas 1,61% ao estrato dominante. Os estratos de arvoretas e superior foram verificados em todas as sete parcelas, enquanto que espécies com maior altura características do estrato dominante foram verificadas apenas nas parcelas 1, 3 e 7 (Figura 1). Cerca de 20% dos indivíduos identificados no levantamento estão localizados na parcela 7 (maior diversidade). Em relação ao aspecto topográfico, esta parcela diferencia-se das demais devido à grande declividade do terreno, além de estar localizada mais próxima ao Rio Uruguai.

Modelagem topográfica

Os resultados obtidos a partir da modelagem topográfica permitiram identificar as principais características do relevo local. A área do PET apresenta uma amplitude altimétrica de 360 m, variando de 120 m nas áreas próximas ao rio Uruguai, até 480m nas regiões mais elevadas (Figura 2a). Essa amplitude mostra que, em algumas áreas, o PET apresenta um relevo dissecado com encostas íngremes, associadas aos vales. A orientação de vertentes (Figura 2b) mostrou que as faces com vertentes orientadas na direção Norte ocupam 19,55% da área de estudo, com as vertentes a Leste apresentando menor área ocupada (7,68%). As demais direções apresentação uma cobertura média de 12,50% (±3,78%). A direção de orientação é parâmetro fundamental para o cálculo do relevo sombreado e afeta a quantidade de radiação disponível para a vegetação (entre outros fatores). Dessa forma, pode alterar os valores da reflectância em imagens de satélite. Considerando a trajetória do sol em um dia no Hemisfério Meridional, as vertentes Leste, Nordeste, Norte, Noroeste e Oeste, que recebem maior incidência de radiação solar, correspondem a 64,29% da área total do PET. As vertentes direcionadas a Sudoeste, Sul e Sudeste somam 35,72% da área do PET, sendo aquelas com menor incidência de radiação.

Figura 2 Resultados da modelagem topográfica do modelo ASTER GDEM v2 para (a) altimetria e (b) orientação de vertentes. 

Figure 2 Results of the topographic modelling with ASTER GDEM v2 for the (a) altimetry and (b) aspect. 

Devido a latitude intermediária, o PET está sujeito a uma grande amplitude do AZS ao longo do ano, que pode ultrapassar 30° de julho para dezembro. Tamanha variação leva a grandes modificações da imagem de relevo sombreado (Figura 3). Verifica-se uma área sombreada muito maior no inverno (Figura 3a) quando comparado com as cenas adquiridas no final de outubro, com o Sol mais próximo do zênite (Figura 3b). As mudanças ocorridas entre as duas datas estão diretamente ligadas a diminuição do AZS de 52° (junho) para de 19° (outubro), alcançando uma amplitude solar de 33°. Essa variação leva a um maior efeito do sombreamento sobre a resposta espectral da vegetação nas cenas RapidEye adquiridas em junho.

Figura 3 Relevo sombreado para (a) 28 de junho e (b) 17 de outubro, com ângulo zenital solar (AZS) de 52° e 19°, respectivamente. 

Figure 3 Shaded relief for (a) June 28 and (b) October 17, 2012 with solar zenith angle (AZS) of 52º and 19°, respectively. 

Análise espectral em função da modelagem topográfica

As variações da reflectância de superfície estão associadas com as condições de sombreamento ao longo do ano e ao efeito da perda de folhas das espécies dominantes (^{SEMA/RS, 2005}; ^{Breunig et al., 2015}). Este efeito combinado sobre os produtos de sensoriamento remoto foi discutido por ^{Breunig et al. (2015)}, com base na avaliação de imagens do sensor Multi-angle Imaging Spectroradiometer (MISR) e com modelos de transferência radiativa. Os resultados reportados são concordantes com os vistos na inspeção dos espectros de reflectância média do RapidEye (Figura 4), adquiridos para o período de maior AZS (Figura 4a) e menor AZS (Figura 4b). Nas duas datas foi observada a influencia das condições topográficas locais, na qual as áreas sujeitas a maior iluminação apresentam maior reflectância. As áreas sombreadas apresentam reflectância mais baixa, conforme esperado. Contudo, quando observada a variação no inverno, esses efeitos são muito mais pronunciados quando comparados com os vistos em meados de outubro de 2012.

Figura 4 Espectros de relectância de loresta estacional decidual para (a) 27 de junho e (b) 17 de outubro de 2012. Cada espectro refere-se a média de 1000 pixels aleatoriamente obtidos na área de cada classe de iluminação. As barras de desvio padrão são indicadas. 

Figure 4 Reflectance spectra of seasonal deciduous forest for (a) June 27 and (b) October 17, 2012. Each spectrum refers to the average of 1000 pixels randomly selected over each illumination class. The standard deviation bars are indicated. 

A diferença de valores de reflectância entre junho e outubro deve considerar o caráter decidual da floresta combinado ao efeito da topografia local. Em junho, o dossel apresenta perda de folhas e um sombreamento bastante elevado, levando a grandes variações nos espectros (Figura 4a). Com o início do brotamento de folhas na primavera, ocorre uma homogeneização do dossel das espécies dominantes de forma a reduzir a anisotropia do dossel, associado a um aumento de IAF (^{Epiphanio et al., 1995}; ^{Walter-Shea et al., 1997}). Em paralelo, ocorre uma redução de aproximadamente 33° no AZS de junho para outubro, levando a uma redução do efeito de iluminação e tornando a cena mais homogênea, mesmo em áreas sujeitas a sombreamento em outubro (Figura 4b).

A distribuição da reflectância das bandas do vermelho e infravermelho próximo (IVP) na área do PET pode ser comparada (Figura 5). Em todas as bandas analisadas, a cena da data de 28 de junho (Figuras 5a e 5c, respectivamente) parece ser mais influenciada pelas condições topográficas locais do que a cena de 17 de outubro (Figuras 5b e 5d). As bandas mais influenciadas por este efeito foram a do vermelho e do IVP. Além das características do relevo, é importante considerar o efeito causado pelas variações causadas pela caducifólia de algumas espécies dominantes do PET. O trabalho de campo mostrou que, em junho, a floresta apresenta um menor índice de área foliar que em outubro. Ademais, as características florísticas e fitossociológicas introduziram variações nos padrões de reflectância para as distintas parcelas.

Figura 5 Distribuição espacial da relectância do vermelho (a e b) e IVP (c e d) da floresta estacional decidual, em função da geometria de iluminação solar para as datas de 28 de junho e 17 de outubro de 2012. 

Figure 5 Spatial distribution of red (a and b) and near infrared (NIR) (c and d) reflectance of the seasonal deciduous forest, as a function of the illumination geometry for June 28 and October 17, 2012. 

No que tange a orientação das vertentes, aquelas orientadas na direção Norte/Noroeste (faces iluminadas) apresentaram menores variações de reflectância entre junho e outubro, quando comparadas as vertentes Sul-Sudeste (faces sombreadas). Nas vertentes sul/sudoeste, o sombreamento predomina. As variações na reflectância entre datas são explicadas pelas mudanças fenológicas e florístico-estruturais da floresta estacional decidual do PET e pelo efeito da topografia local.

A avaliação dos gráficos dos valores médios amostrados de reflectância em função das variáveis topográficas (Figura 6) mostrou a existência de grandes tendências admitidas pela reflectância à medida que os parâmetros da modelagem topográfica são alterados. Os gráficos das Figuras 6a e 6b permitiram verificar que as vertentes dispostas nas direções contrárias à da incidência direta de radiação solar apresentaram valores de reflectância consideravelmente menores em relação as vertentes dispostas nas direções com incidência direta de radiação solar.

Figura 6 Valores de reflectância normalizada das bandas RapidEye do azul, vermelho e infravermelho próximo, mostrados em função da orientação de vertentes, altimetria e relevo sombreado em (a), (c) e (e) para 28 de junho e em (b) (d) e (f) para 17 de outubro de 2012, respectivamente. 

Figure 6 Normalized reflectance values of the blue, red and near infrared RapidEye bands, displayed as a function of the aspect, altimetry and shaded relief in (a), (c) and (e) for June 28, and in (b), (d) and (f ) for October 17, 2012, respectively. 

Considerando as bandas utilizadas na elaboração dos índices de vegetação EVI e NDVI, verificou-se que os valores de reflectância apresentaram maior amplitude de variação entre vertentes dispostas em condições distintas de iluminação/sombreamento na data de junho, sendo verificadas amplitudes de reflectância de 0,57%, 0,84% e 14,65%, respectivamente, para as bandas do azul, vermelho e infravermelho próximo. Na data de outubro, para as mesmas bandas, foram verificadas variações de reflectância de 0,17%, 0,27% e 5,99%. Desta forma, pode ser verificado uma maior isotropia do dossel em outubro. Os valores médios de reflectância amostrados em função da altimetria do PET (Figura 6c e 6d) permitiram verificar que as bandas do azul e vermelho apresentaram um comportamento de redução da reflectância em relação as áreas de maior altimetria do PET. No entanto, a banda do IVP apresentou um padrão contrário, com valores mais elevados de reflectância nas áreas mais altas. Essas duas tendências podem estar associadas a maior absorção por pigmentos em função do maior adensamento do dossel nas áreas mais elevadas. As tendências foram identificadas nas duas datas avaliadas, porém, em outubro com menor variação em função da altitude. É importante destacar que não existe uma relação direta entre altitude e reflectância, mas sim uma relação entre as características da vegetação para distintas altitudes (p.ex., distribuição de espécies), que induzem variações na reflectância.

O cálculo da imagem do cosseno da iluminação (relevo sombreado) mostrou que nas duas datas avaliadas (Figuras 6e e 6f), os maiores valores de reflectância encontram-se nas áreas iluminadas por meio de incidência direta de radiação. Por sua vez, os valores mais baixos foram verificados nas áreas com condições de iluminação intermediárias ou iluminadas de forma indireta. Os baixos valores de reflectância restringiram-se as áreas com grandes condições de sombreamento, conforme esperado. No entanto, devido as diferenças de AZS entre as duas datas influenciar diretamente na redução das condições de sombreamento, a amplitude de variação dos valores de reflectância sofreu uma redução entre as datas, sendo verificadas amplitudes de 0,39%, 0,52% e 9,22%, respectivamente, para as bandas do azul, vermelho e infravermelho próximo. Na data de outubro, as variações foram de 0,1% para a banda do azul, 0,19% para a banda do vermelho e 4,36% para a banda do IVP. Na data de outubro, devido ao maior brotamento de folhas (aumento da clorofila), foi verificada uma maior absorção por parte das bandas do visível, condicionando uma menor amplitude de variação, em virtude da baixa reflectância. A banda do IVP, no entanto, apresentou elevados valores de reflectância, tendo uma menor influência de baixos valores de reflectância de áreas sombreadas, uma vez que em outubro as áreas sombreadas restringiram-se apenas aos fundos dos vales.

A análise das três bandas avaliadas mostrou que as bandas do vermelho e do IVP foram as que apresentaram maior sensibilidade tanto em relação as condições de estado fenológico da vegetação, quanto as variações de iluminação. Neste sentido, a elaboração dos transectos a partir de vertentes dispostas nos sentidos Norte-Sul e Leste-Oeste mostraram que a orientação das vertentes tende a condicionar as variações da reflectância em função das condições de iluminação associadas a cada direção das vertentes (^{Matsushita et al., 2007}).

Em relação a banda do vermelho, os valores de reflectância para a data de junho (Figura 7a) apresentaram uma grande amplitude de variação associada a direção das vertentes. Na data de junho, os valores de reflectância tendem a ser mais elevados em relação aos valores analisados na data de outubro. Na primeira data ocorrem condições geométricas que propiciam o aumento do sombreamento e a caducifólia das espécies que constituem o estrato arbóreo dominante, contribuindo para elevar a reflectância na banda do vermelho. Em outubro, em função da absorção da clorofila na banda do vermelho, os valores de reflectância foram baixos em sua maioria (Figura 7c). No entanto, da mesma forma como verificado na data de junho, a vertente orientada ao norte apresentou os valores mais elevados de reflectância, sendo que a vertente com direção sul concentrou os valores mais baixos. Em geral, foi verificado um forte condicionamento da reflectância em função da direção das vertentes.

Figura 7 Variação da reflectância das bandas do vermelho e infravermelho próximo do RapidEye em função das características do relevo para um transecto Norte-Sul na data de junho (a, b) e outubro (c, d) de 2012. 

Figure 7 Reflectance variation of the red and near infrared RapidEye bands as a function of relief characteristics over a North-South transect for June (a, b) and October (c, d), 2012. 

Na data de junho, o transecto contendo a banda do IVP (Figura 7b) mostrou grande variação nos valores de reflectância. Esta variação mostrou-se relacionada diretamente às condições de iluminação condicionadas pelo relevo, sendo que a vertente com direção norte, que recebe iluminação direta, apresentou os valores mais elevados. Por sua vez, a vertente orientada na direção sul, na qual a incidência de radiação solar é baixa, faz com que ocorra sombreamento ou menor resposta espectral da vegetação. Desta forma, a análise da reflectância das bandas do vermelho e do IVP aponta que para junho ocorre um condicionamento de reflectância pelas condições de iluminação, sendo este um fator importante na resposta da reflectância da vegetação, além de seus próprios atributos biofísicos. Em outubro, a reflectância da banda do IVP (Figura 7d), embora apresente-se elevada e com menor amplitude de variação em relação a data anterior, possui valores com a mesma tendência em relação a orientação das vertentes e ao condicionamento da incidência de iluminação causada por estas.

Embora as condições de iluminação tendam a mostrar-se como fatores dominantes sobre as características da reflectância, os valores encontrados nos transectos permitiram identificar a ocorrência de uma redução na amplitude de variação das medidas de reflectância para as duas bandas. No IVP, esta redução pode estar relacionada ao estado fenológico da vegetação, visto que em outubro, em função do período de brotamento das folhas, em especial das espécies caducifólias, tende a ocorrer uma maior isotropia em relação ao cenário de junho. Em relação aos transectos elaborados no sentido Leste-Oeste (Figura 8), foi encontrado um condicionamento menor da iluminação sobre os valores de reflectância, quando comparado com a variação encontrada nas vertentes Norte-Sul.

Figura 8 Variação da reflectância das bandas do vermelho e infravermelho próximo do RapidEye em função das características do relevo em um transecto Leste-Oeste em junho (a, b) e outubro (c, d) de 2012. 

Figure 8 Reflectance variation of the red and near infrared RapidEye bands as a function of relief characteristics over a East-West transect for June (a, b) and October (c, d), 2012. 

Analisando os valores da reflectância da banda do IVP foi identificada uma variação maior entre as medidas, sendo que na data de junho (Figura 8b) foi possível identificar de forma mais clara as diferenças de medidas de reflectância entre as vertentes. As vertentes orientadas na direção leste mostraram reflectância mais elevada, enquanto que a vertente disposta no sentido oeste apresentou valores mais baixos de reflectância, conforme esperado. Na data de outubro (Figura 8d), assim como ocorreu no transecto Norte-Sul, foi identificada uma redução nas amplitudes dos valores pelo estado fenológico da vegetação. Da mesma forma, como verificado na data de junho, as medidas de reflectância em outubro apresentaram-se dependentes das condições de iluminação condicionadas pelo efeito topográfico, porém, em menor magnitude. Analisando a reflectância do vermelho (junho-Figura 8a e outubro- Figura 8c), fica clara sua menor dependência do sombreamento quando comparado com o IVP.

Avaliação dos índices de vegetação

Visto que a floresta subtropical decídua predomina no PET, verifica-se que as caraterísticas locais de topografia controlam em parte as variações dos índices de vegetação (Figura 9). Da mesma forma como verificado por ^{Matsushita et al. (2007)}, o EVI foi mais afetado pelas condições locais de iluminação. Apesar de ser mais sensível às variações estruturais da vegetação, o EVI não normaliza os efeitos externos (topográficos) como faz o NDVI, que é um índice normalizado. De fato, a alta correlação do EVI com a banda do IVP (^{Galvão et al., 2011}) mostrou que o índice não normaliza os dados, mas sim, representa um índice da diferença entre a reflectância do vermelho e do IVP (^{Peng et al., 2011}; ^{Breunig et al., 2015}).

Figura 9 Variações espaciais do NDVI e EVI para as datas de 28 de junho (a, c) e 17 de outubro (b, d) de 2012. 

Figure 9 Spatial variation of NDVI and EVI for June 28 (a, c) and October 17 (b, d), 2012. 

Embora o estudo feito por ^{Matsushita et al. (2007)} comprove que o EVI apresenta grande dependência em relação ao efeito topográfico (Figura 9c e 9d), para a data de junho foi verificada uma dependência do índice NDVI em relação as áreas sombreadas do PET (Figura 9a). No entanto, esta variação foi menor em relação a apresentada pelo EVI mesmo em áreas com elevado sombreamento. Em relação a data de outubro (Figura 9b), foi identificada uma redução da amplitude dos valores da reflectância, indicando que, à medida em que se avança para a primavera e verão, ocorre uma forte redução da dependência do efeito topográfico devido a redução do AZS e brotamento de folhas novas (redução da anisotropia do dossel). O EVI, que apresentou elevada variação de medidas, conforme as diferentes características de iluminação presentes na data de junho (Figura 9c), ainda se mostrou fortemente controlado pelas condições locais de iluminação na data de outubro (Figura 9d).

Cabe destacar que o efeito do brotamento de folhas não foi avaliado para quantificar a interferência de cada variável sobre o comportamento do NDVI e EVI, bem como o total de produtividade primária liquida. A distribuição espacial do índice EVI para a área do PET apresentou as mesmas características descritas por ^{Matsushita et al. (2007)} sobre a sua relação com a influência das condições de iluminação da área no momento da aquisição das imagens. As áreas dos vales e encostas iluminadas são nitidamente visíveis em ambas as datas analisadas.

Tanto o NDVI, quanto o EVI, amostrados em função das variáveis topográficas de orientação de vertentes, altimetria e relevo sombreado, apresentaram variações em relação às condições de sombreamento geradas em função das caraterísticas do relevo do PET (Figura 10). Embora ocorra uma variação natural destes índices entre as duas datas avaliadas em função da mudança nas condições fenológicas da vegetação do PET, as diferenças entre as áreas sombreadas e iluminadas são nítidas e são consequência das variações locais de iluminação, deciduidade-brotamento das espécies dominantes e variações locais da distribuição de espécies. Em geral, as diferenças dos índices de vegetação diminuem de junho para outubro devido a um efeito combinado de redução do AZS e o aumento do índice de área foliar (^{Breunig et al., 2015}).

Figura 10 Variações normalizadas do NDVI (a, c, e) e EVI (b, d, f) em função da orientação de vertentes, altimetria e relevo sombreado para 28 de junho e 17 de outubro de 2012. 

Figure 10 Normalized variations of NDVI (a, c, e) and EVI (b, d, f ) as a function of the aspect, altimetry and shaded relief for June 28 and October 17, 2012. 

O NDVI (Figura 10a) e EVI (Figura 10b) apresentaram menor dependência da orientação de vertentes em outubro com relação a junho devido a redução do sombreamento e ao brotamento de novas folhas. O EVI novamente foi mais afetado, principalmente no inverno. Apesar de não haver relações de causa-efeito entre altimetria e índices de vegetação, foi identificado um padrão estável para o NDVI em junho (Figura 10c) e outubro (Figura 10d). Para o EVI, uma tendência de aumento do índice com a altitude foi encontrada para ambas as imagens RapidEye. Esse padrão pode indicar uma maior sensibilidade do EVI ao brotamento de folhas e a variações na estrutura do dossel (^{Huete et al., 2002}).

A avaliação do NDVI em função do relevo sombreado (Figura 10e) mostrou uma maior amplitude de variação entre cada classe de iluminação na data de junho. Neste sentido, as variações das condições de iluminação tendem a causar uma maior segmentação entre os valores de NDVI associados a cada condição de iluminação. No entanto, na data de outubro verificou-se que as diferenças entre os valores médios amostrados ocorreram de forma menos acentuada, o que mostra que a maior incidência de radiação direta provoca um aumento da reflectância, bem como a redução das áreas sombreadas, reduzindo o efeito topográfico sobre o NDVI. Por meio da análise dos valores médios amostrados de EVI em função do relevo sombreado (Figura 10f) foi verificado que este índice apresentou uma amplitude de variação maior que o NDVI. Os valores de junho apresentaram grande amplitude de variação associados as classes de iluminação, sendo o valor médio das áreas sombreadas o que apresentou maior variação. Na data de outubro, a amplitude dos valores das classes diminuiu consideravelmente, embora ainda a amplitude entre os valores médios das classes seja maior em comparação com os valores obtidos a partir do índice NDVI. O efeito do sombreamento apresentou um controle mais efetivo sobre os valores do que as áreas com iluminação direta. De forma sintética, os perfis identificados pelos símbolos vazados da Figura 10 tendem a variar pouco ao redor de 1, enquanto que os dados de junho (símbolos sólidos) apresentam grandes variações.

De acordo com os transectos elaborados sobre vertentes orientadas na direção Norte-Sul, os índices de vegetação NDVI e EVI foram influenciados pelas condições de iluminação (Figura 11). Em tese, os índices de vegetação deveriam ressaltar a vegetação de acordo com as suas características fenológicas e eliminar os efeitos externos. No entanto, quando analisados sobre condições diferentes de iluminação sobre o relevo, os índices NDVI e EVI apresentaram variações de acordo com a diferença de incidência da radiação solar em relação a orientação das vertentes. A Figura 11a mostra que o NDVI de junho apresentou uma variação de acordo com a orientação das vertentes no transecto. A vertente com orientação na direção norte apresentou os valores mais elevados de NDVI, enquanto que valores mais baixos foram encontrados na vertente oposta (sul), caracterizada por apresentar condições de sombreamento causadas pela relação entre o efeito topográfico e a geometria de iluminação. Na data de outubro (Figura 11c), no entanto, a variação verificada foi menor do que na data anterior, ocorrendo valores muito próximos em ambas as vertentes.

Figura 11 Variação dos índices de vegetação em função das características do relevo, em um transecto norte-sul para o (a) NDVI e (b) EVI de junho e (c) NDVI e (d) EVI de outubro. 

Figure 11 Variation of vegetation indices as a function of relief characteristics in a North-South transect for the (a) NDVI and (b) EVI of June and (c) NDVI and (d) EVI of October. 

O EVI analisado no transecto Norte-Sul (Figuras 11b e 11d) mostrou ser mais sensível as condições de sombreamento, em decorrência da grande amplitude de variação verificada entre as vertentes. Na vertente orientada na direção norte, os valores do EVI mantiveram-se elevados, enquanto que na vertente oposta, os valores encontrados foram mais baixos. Embora se tratando de uma área coberta por floresta, na qual os índices de vegetação deveriam apresentar valores elevados, os efeitos das condições de iluminação afetam a resposta do EVImesmo na data de outubro (Figura 11d), em que há uma redução dos efeitos topográficos causada pela diminuição do AZS.

Em relação ao transecto Leste-Oeste (Figura não apresentada), o NDVI apresentou-se mais estável em relação ao condicionamento da iluminação causada pelo efeito topográfico. Na data de junho foram identificados baixos valores do índice associados a vertente oeste, enquanto que na data de outubro os valores apresentaram-se de forma quase que homogenia entre as vertentes, não havendo variações consideráveis. Novamente, o índice EVI apresentou grande dependência das condições de iluminação dispostas pelas vertentes, sendo que os maiores valores se concentraram em ambas as datas nas vertentes com orientação leste, as quais apresentam maior incidência de radiação solar em relação as vertentes orientadas na direção oeste. De maneira geral, entre as duas datas analisadas, as mudanças ocorridas nos valores amostrados de reflectância e nos índices de vegetação estão diretamente relacionadas a variação sazonal das condições de iluminação do PET em função do decréscimo do AZS verificado entre as duas datas. De acordo em ^{Ayoade (2011)}, a variação do ângulo de incidência solar associada à latitude do lugar é considerado o principal fator que condiciona as variações de iluminação das áreas. Cabe destacar que as caracteristicas sazonais da vegetação podem influenciar a resposta espectral, resultando em um efeito combinado da geometria e iluminação e fenologia (^{Breunig et al., 2015}).

Com o levantamento florestal das sete parcelas amostrais foi possivel perceber que embora o PET seja totalmente coberto por vegetação, o dossel não é homogêneo. Isso se deve a sua estratificação em três estratos arbóreos principais, mais um arbustivo e outro herbáceo, o que configura uma característica de grande sombreamento interno do dossel vegetativo implicando em alterações diretas na resposta espectral da vegetação.

Devido a estas condições, podem ser verificadas variações na resposta das bandas, em especial na banda do vermelho, causadas pelas estruturas não foliares, pelo efeito de substrato (background) e pelo sombreamento sobre os estratos inferiores do dossel (causado pelas árvores de maior porte). Estas variações influenciam na resposta dos índices de vegetação, nos quais o aumento da resposta na banda do vermelho tende a reduzir os valores dos índices.

Embora exista a ocorrência destas características de sombreamento interno do dossel, o sombreamento causado pelo relevo é o que mais influenciou na resposta tanto da reflectância das bandas, quanto na dos valores dos índices de vegetação, formando áreas de maior e menor reflectância. O sombreamento interno afeta de forma mais direta as variações entre pixel vizinhos, especialmente considerando sensores de resolução espacial alta como o RapidEye. Novos estudos, com um maior controle de campo, podem auxiliar no entendimento do efeito do sombreamento interno em florestas caducifólias.

Em geral, o resultado mostra que os índices de vegetação avaliados podem ser influenciados pelas variações locais de topografia, alterando as estimativas de parâmetros biofísicos e bioquímicos da vegetação (^{Galvão et al., 2016}).