Acessibilidade / Reportar erro

Uso de imagens de radar para o cálculo da produção primária de plantas aquáticas nas várzeas da Amazônia

Use of radar imagery for estimating net primary productivity of aquatic vegetation in the Amazon floodplain

Resumos

A área da planície de inundação da Amazônia é estimada em 300 000km² e sua produtividade primária em 1,17 x 10(14) g C yr-1. Deste total de área e produtividade, estimativas sugerem que 43% e 62%, respectivamente, são atribuídos às plantas aquáticas. Estas estimativas variam de acordo com o pulso de inundação. Por exemplo, durante o período de seca as plantas terrestres (herbáceas) geralmente ocupam áreas que apresentam plantas aquáticas na cheia. A área e a produtividade destes ecossistemas são informações essenciais para a compreensão da dinâmica biogeoquímica da Amazônia. Imagens de satélites (radar) combinadas com amostragem de campo foram utilizadas para estimar a biomassa e mapear a área de cobertura de plantas aquáticas emergentes para calcular a produção primária de plantas aquáticas na várzea do baixo Amazonas. A combinação de bandas C e L forneceu a melhor correlação (r=0,82) e um ponto de saturação de biomassa intermediário (620 gm-2) para estimar biomassa aérea. O método de segmentação e classificação por região foi utilizado para classificar combinações de bandas C e L para cada período de nível de água, e forneceu uma precisão de mapeamento maior que 95% para determinação espacial de áreas cobertas por plantas aquáticas. Combinando a distribuição espacial de plantas aquáticas, o modelo para estimativa de biomassa aérea e a porcentagem de biomassa submersa, estimou-se espacialmente uma produção primária líquida anual de 1.9x10(12) g C yr-1 (±28%) para as plantas aquáticas em uma área de 394km².

Imagens de satélite radar; produção primária; plantas aquáticas; várzea


Estimates suggest that 43% and 63% of the 300,000km² of Amazon floodplain area and of the 1.17 x 10(14) g C yr-1 of primary productivity, respectively, are attributed to the aquatic vegetation. These numbers change according to the flood pulse. For instance, in the dry season terrestrial herbaceous plants generally occupy areas that have lost their aquatic vegetation. The areal extent and productivity of these ecosystems are essential to construe even an initial understanding of the biogeochemistry of the Amazon. Field measurements were combined with synthetic aperture radar images to evaluate the use of RADARSAT and JERS-1 for estimating biomass changes and mapping of aquatic vegetation, and subsequently estimating the net primary productivity of aquatic vegetation in the lower Amazon. The combination of C and L bands provides the best correlation (r =0.82) and an intermediate saturation point (620 gm-2) for estimating above water biomass of aquatic vegetation. A combination of RADARSAT and JERS-1 images from each water period was classified using a region growing algorithm, and yielded accuracy higher than 95% for the seasonal vegetated areas of the floodplain. The combination of the seasonal mapped area of aquatic vegetation with the statistical SAR-algorithm for estimating above water biomass and the percentage of below water biomass yielded a total annual NPP of 1.9x10(12) g C yr-1 (±28%) for aquatic vegetation.

Radar satellite imagery; primary production; aquatic vegetation; várzea


PROJETO LBA

Uso de imagens de radar para o cálculo da produção primária de plantas aquáticas nas várzeas da Amazônia

Use of radar imagery for estimating net primary productivity of aquatic vegetation in the Amazon floodplain

Maycira Costa

University of Victoria, Canada. Department of Geography. P.O. Box 3050, Victoria, BC, Canada. V8W 3P5. email:maycira@office.geog.uvic.ca Instituto Nacional de Pesquisas Espaciais - INPE - Divisão de Sensoriamento Remoto

RESUMO

A área da planície de inundação da Amazônia é estimada em 300 000km2 e sua produtividade primária em 1,17 x 1014 g C yr-1. Deste total de área e produtividade, estimativas sugerem que 43% e 62%, respectivamente, são atribuídos às plantas aquáticas. Estas estimativas variam de acordo com o pulso de inundação. Por exemplo, durante o período de seca as plantas terrestres (herbáceas) geralmente ocupam áreas que apresentam plantas aquáticas na cheia. A área e a produtividade destes ecossistemas são informações essenciais para a compreensão da dinâmica biogeoquímica da Amazônia. Imagens de satélites (radar) combinadas com amostragem de campo foram utilizadas para estimar a biomassa e mapear a área de cobertura de plantas aquáticas emergentes para calcular a produção primária de plantas aquáticas na várzea do baixo Amazonas. A combinação de bandas C e L forneceu a melhor correlação (r=0,82) e um ponto de saturação de biomassa intermediário (620 gm-2) para estimar biomassa aérea. O método de segmentação e classificação por região foi utilizado para classificar combinações de bandas C e L para cada período de nível de água, e forneceu uma precisão de mapeamento maior que 95% para determinação espacial de áreas cobertas por plantas aquáticas. Combinando a distribuição espacial de plantas aquáticas, o modelo para estimativa de biomassa aérea e a porcentagem de biomassa submersa, estimou-se espacialmente uma produção primária líquida anual de 1.9x1012 g C yr-1 (±28%) para as plantas aquáticas em uma área de 394km2.

Palavras-chave: Imagens de satélite radar, produção primária, plantas aquáticas, várzea.

ABSTRACT

Estimates suggest that 43% and 63% of the 300,000km2 of Amazon floodplain area and of the 1.17 x 1014 g C yr-1 of primary productivity, respectively, are attributed to the aquatic vegetation. These numbers change according to the flood pulse. For instance, in the dry season terrestrial herbaceous plants generally occupy areas that have lost their aquatic vegetation. The areal extent and productivity of these ecosystems are essential to construe even an initial understanding of the biogeochemistry of the Amazon. Field measurements were combined with synthetic aperture radar images to evaluate the use of RADARSAT and JERS-1 for estimating biomass changes and mapping of aquatic vegetation, and subsequently estimating the net primary productivity of aquatic vegetation in the lower Amazon. The combination of C and L bands provides the best correlation (r =0.82) and an intermediate saturation point (620 gm-2) for estimating above water biomass of aquatic vegetation. A combination of RADARSAT and JERS-1 images from each water period was classified using a region growing algorithm, and yielded accuracy higher than 95% for the seasonal vegetated areas of the floodplain. The combination of the seasonal mapped area of aquatic vegetation with the statistical SAR-algorithm for estimating above water biomass and the percentage of below water biomass yielded a total annual NPP of 1.9x1012 g C yr-1 (±28%) for aquatic vegetation.

Key words: Radar satellite imagery, primary production, aquatic vegetation, várzea.

INTRODUÇÃO

A planície de inundação do Rio Amazonas faz parte da vasta área de regiões alagadas tropicais e subtropicais, as quais compõem aproximadamente 56% das áreas alagadas do globo (Mitsch & Gosselink, 2000). A área da planície de inundação da Amazônia é estimada em 300 000km2 (Junk, 1993) e sua produtividade primária em 1,17 x 1014 g C yr-1 (Melack & Forsberg, 2000). Deste total de área e produtividade, estimativas sugerem que 43% e 62%, respectivamente, são atribuídos as plantas aquáticas(Melack & Forsberg, 2000), as quais se proliferam predominantemente na regiao de várzea (Piedade et al. 1991). Esta alta produtividade é decorrente da alta eficiência destas plantas em converterem a energia solar interceptada em biomassa, característica essencial para a sua sobrevivência em um ambiente de contínua variação de nível da água devido ao pulso de inundação (Piedade et al. 1991).

Devido à heterogeneidade espacial da planície de inundação, estimativas mais refinadas de área e produtividade só são possíveis com uma intensa amostragem de campo ou com o uso de imagens de satélite associado à amostragem de campo. A primeira opção se torna inviável devido às dimensões espaciais da Amazônia. A segunda opção é apresentada nesta pesquisa. Outros autores (Potter et al. 1993; Sellers et al. 1995) utilizaram com sucesso dados de satélites ópticos em modelos de estimativa de produção primária em estudos regionais e globais. Contudo, imagens de satélites ópticos sem cobertura de nuvens não estão disponíveis na frequência necessária para a Amazônia. Por exemplo, Novo et al. (1997) utilizaram imagens do satélite Landsat adquiridas durante 10 anos para construir um mosaico sem nuvens da várzea do Amazonas. Para a região geográfica que estamos estudando, em um levantamento de quatro anos de aquisição de imagens Landsat (um total de 38 imagens), somente 3% das imagens não apresentam nuvens; as outras apresentam aproximadamente 100% de cobertura de nuvens. Portanto, o uso de imagens de satélite radar é uma das únicas soluções para estudos temporais na Amazônia.

Este trabalho relata um método eficiente que utiliza imagens de satélites radares para mapeamento sazonal e estimativa de produção primária líquida (PPL) de plantas aquáticas na Amazônia.

MÉTODOS

Modelo para estimativa de PPL usando imagens de radar

A PPL foi estimada baseada na seguinte equação:

onde Bt e Bt+1 correspondem à biomassa seca total no começo e final do intervalo temporal de medida e Lt+1 é a porcentagem de biomassa seca perdida durante o mesmo intervalo. Todos os parâmetros são calculados para cada célula de resolução espacial (pixel) da imagem de radar correspondente à área coberta por uma determinada planta aquática. Esta mesma equação foi originalmente utilizada por Junk & Piedade (1993), porém, os autores basearam-se unicamente em amostragem de campo.

Os parâmetros necessários para o uso da equação 1 são:

(i)A área espacial ocupada pela plantas aquáticas

(ii)Medidas mensais de biomassa e perda de biomassa

A explicação detalhada da estimativa destes parâmetros utilizando imagens de radar é apresentada nos próximos itens.

Área de estudo

A área de estudo, Lago Grande de Monte Alegre (LGMA), (Figura 1), está localizada nas coordenadas 2º00'S/54º00'O e 2º30'S/54º30'O, próximo à cidade de Santarém. A vegetação predominante na planície de inundação é constituída por grandes bancos de plantas aquáticas herbáceas, arbustos pioneiros, e floresta. Na região de "terra firme" (regiões marginais não alagáveis), a vegetação é composta predominantemente por pequenos campos agrícolas, cerrados e floresta secundária (RADAMBRASIL, 1976). Hymenachne amplexicaulis (rabo de rato), Echinochloa polystachya (canarana) e Paspalum repens (premenbeca) são as espécies dominantes de plantas aquáticas na planície de inundação. Este estudo é focada na espécie H. amplexicaulis por ser a mais comum na região (corresponde a 65% das amostras coletadas). Esta espécie forma grandes bancos mono-específicos na zona litoral do lago, onde a profundidade média é aproximadamente 3,5±0,5m em maio (água alta) e 0m em outubro (água baixa). Esta espécie apresenta uma fase terrestre e está distribuída largamente em toda a várzea Amazônica, principalmente nos lagos e áreas susceptíveis à sedimentação (Junk & Piedade, 1993).


A figura 2 apresenta o ciclo monomodal do Rio Amazonas, na altura de Óbidos, estação mais próxima da área de estudo, onde os níveis mais altos de água ocorrem em maio/junho e os níveis mais baixos em outubro/novembro. As curvas do ciclo hidrológico também mostram variações similares do nível da água para a mesma época do ano para os três anos de aquisição de dados de campo e satélite. Devido a essa variação similar nos níveis de água foi possível assumir que os dados adquiridos (satélite e campo) representam um único ciclo hidrológico contínuo, isto é, estamos assumindo que os dados foram adquiridos durante um ano do ciclo hidrológico. Esta consideração é importante porque a variação temporal da altura da água é de importância fundamental no ciclo de crescimento das plantas aquáticas na Amazônia (Junk & Piedade, 1997).


Dados de campo

Um total de cinco campanhas de campo foram realizadas coincidentes com o período de água alta (maio 1996), descida de água (agosto 1996), água baixa (novembro 1996), subida de água (abril 1997), e novamente, água alta (junho 1999). Pontos amostrais foram estabelecidos em regiões de grandes bancos homogêneos de plantas aquáticas utilizando o método de estratificação aleatória (Figura 1 indica as regiões amostrais onde um total de 96 pontos amostrais foram determinados). Nos pontos amostrais foram coletadas amostras para determinação de biomassa aérea e submersa usando um quadrado de 0,25m2 de área. Outros parâmetros como altura das plantas, espécie dominante, estágio fenológico e profundidade local foram determinados. As amostras de biomassa foram apropriadamente secas e pesadas para determinação de biomassa seca (Dowing & Anderson, 1985). O valor de biomassa seca submersa total foi estimado considerando a biomassa submersa de três plantas individuais, o comprimento de talos submersos individuais, o número de plantas por quadrante e a profundidade local. A correlação entre comprimento das plantas e profundidade local é 0,95, indicando que a medida de profundidade pode ser utilizada para estimar o comprimento total das hastes submersas (Costa, 2000).

Imagens de satélite

Imagens dos satélites Radarsat e JERS foram adquiridas coincidentes com as campanhas de campo. As características destas imagens são apresentadas de forma sucinta na tabela 1. Condições de nenhuma chuva ou pouca chuva foram observadas durante três dias anteriores e durante o período de aquisição das imagens, de forma que as condições meteorológicas não influenciaram a qualidade radiometrica das mesmas.

Dois objetivos determinaram as etapas necessárias para processamento das imagens: (i) estimativa de biomassa a partir de valores de retroespalhamento do radar e (ii) determinação da distribuição espacial de plantas aquáticas a partir das imagens de radar. As etapas de processamento são apresentadas a seguir.

Calibração radiométrica

O processo de calibração radiométrica implica no uso de constantes determinadas para as imagens do Radarsat e JERS. As equações a seguir foram utilizadas nesta etapa:

Onde sº é o valor de retroespalhamento calibrado radiometricamente em unidade linear a partir do valor de nível digital (ND), Ao e Agj são os valores de offset e ganho, èj é o ângulo de incidência levando em consideração o elipsoide da Terra. k é a constante temporal aplicada (-68.2). A equação 1 foi aplicada nas imagens do Radarsat (Shepherd, 1997) e a equação 2 nas imagens JERS (Shimada, 1996).

Correção geométrica

Utilizou-se o método de orto-retificação para correção geométrica e o algoritmo de convolução cúbica na reamostragem. A precisão média final dos pontos teste é de 21 m. Após as etapas de correção, dois processamentos distintos foram realizadas nas imagens (i) determinação dos valores de sº em unidade logarítmica (dB) a partir dos valores em unidade linear e (ii) processamento para segmentação e classificação.

Determinação dos valores de sº

Valores médios de sº foram determinados para as estacões de coleta de campo. O número médio de pixels em cada amostra é de 208 (±89SD), o que significa que as médias amostrais apresentam um nível de significância de 95% de acordo com a estatística para radar e considerando o efeito do "speckel" (Laur et al. 1996).

A partir dos valores de sº, o "índice de biomassa" para cada amostra foi relacionado estatisticamente com os valores de biomassa aérea. Subsequentemente, o modelo estatístico foi avaliado de acordo com e erro quadrático médio (EQM) (Neter et al. 1990).

Segmentação e Classificação - determinação espacial das áreas de plantas aquáticas

Antes da classificação as imagens foram (i) filtradas com o filtro Gamma (três interações recursivas com uma janela 5x5) para minimizar o efeito do "speckel" (Costa et al. 2002).

Subsequentemente, as imagens foram (ii) convertidas de 32 bit para 8 bit de acordo com a média e três vezes o desvio padrão da média (99% da distribuição de níveis digitais originais) para reduzir o tamanho dos arquivos de imagens. A seguir, (iii) máscaras previamente determinadas correspondentes à superfície de água e "terra firme" (região acima da planície de inundação) foram subtraídas das imagens. Estas máscaras foram determinadas através da classificação das imagens adquiridas no período de água baixa e alta, respectivamente. Subsequentemente, (iv) pares de imagens sazonais (JERS e Radarsat) foram submetidos a um processo de segmentação automática de acordo o algoritmo de crescimento de regiões (Bins et al. 1996). O processo de segmentação requer a definição do limiar de área (mínimo número de pixels em uma região) e do limiar de similaridade de nível digital (diferença média do valor de nível digital entre duas regiões) para separação de duas regiões. O limiar de área foi determinado em 50 pixels devido à (i) necessidade de garantir que os dados estatísticos de uma região (definida pelo número de pixels) estão de acordo com os critérios de uma população com distribuição normal para imagens de radar (Laur et al. 1996) e (ii) evitar a segmentação exagerada das imagens. O limiar de similaridade de nível digital foi estabelecido de acordo com o método da diferença significativa de duas médias (Snedecor & Cochran, 1980), e resultou em um valor de 30. Finalmente, (v) os pares de imagens segmentadas foram submetidos a uma classificação supervisionada de regiões baseada no método da distância de Bahattacharrya (Richards, 1986). A tabela 2 apresenta os números de regiões de treinamento e teste utilizadas para a classificação supervisionada e a matriz de erro da classificação, respectivamente. Detalhes do processo de seleção de regiões podem ser encontrados em Costa (2004).

RESULTADOS

Biomassa de plantas aquáticas

A tabela 3 apresenta os valores medidos de biomassa seca aérea, submersa e total para H. amplexicaulis, a espécie mais comum na área de estudo (65% das amostras). A biomassa submersa corresponde à biomassa do talo submerso e a biomassa da raiz (5% da biomassa total de acordo com Junk & Piedade (1997).

De uma forma geral, a biomassa das plantas aumentou de forma constante de novembro até agosto, seguindo a subida e descida dos níveis de água. De fato, as diferenças sazonais de biomassa são principalmente devido ao aumento da biomassa submersa. A biomassa submersa contribui para a biomassa total o correspondente a 52%, 88%, 86%, 81% e 83%, respectivamente, para novembro, abril, maio, junho e agosto. Esta alta porcentagem é uma estratégia de sobrevivência desta espécie para acompanhar a subida da água. Estratégia similar é observada para outras espécies de gramíneas fixas ao substrato da Amazônia (Piedade et al. 1991).

Imagens de radar

Determinação dos valores de retroespalhamento para estimativas de biomassa

O modelo estatístico gerado, log B a = 4,002 + 0,175 H sº índice, resultou em um coeficiente de correlação igual a 0,82. Os resultados do modelo estatístico são apresentados na tabela 4.

O teste do modelo estatístico foi realizado com 25% das amostras do modelo e avaliado de acordo com o erro quadrático médio (EQM) (Neter et al. 1990). Os valores de EQM (modelo) e EQM* (teste do modelo) são similares, o que demostra que o modelo estatístico é estável e pode ser utilizado para as estimativas de biomassa.

A curva do modelo estatístico é apresentada na figura 3. A curva sugere a existência de um ponto de saturação para os valores de sº. Baseado na solução da segunda derivada da curva do polinômio de terceira ordem, definiu-se o ponto de saturação da biomassa aérea em 620 gm-2. A definição deste ponto é importante para entender a capacidade de estimativa do modelo estatístico, porque acima deste valor as estimativas de biomassa aérea através do modelo não são confiáveis.


Pontos de saturação similares (400gm-2) foram determinados para estimativa de biomassa aérea de arroz em campos alagados utilizando imagens de radar do ERS-1 (Le Toan et al. 1997). Contudo, estimativas utilizando imagens do Radarsat adquirida com ângulo de incidência menor (23º) resultou em um ponto de saturação inferior (200gm-2), para o mesmo tipo de cultura de arroz (Ribbes & Le Toan, 1999). Estas diferenças podem estar relacionadas com a estrutura de plantio de arroz quando comparado com o processo natural de crescimento de plantas aquáticas na Amazônia. O plantio de arroz estudado por Ribbes & Le Toan (1999) apresenta uma distribuição de 20 agregados de plantas por m-2, distribuídos em fileiras. Na área de estudo, H. amplexicaulis não forma agregados, a densidade é em média 111 talos m-2 (Costa, et al. 2002), e não existe uma orientação horizontal preferencial de crescimento. Também, a complementaridade das bandas C e L (este estudo) deve influenciar na determinação de um ponto de saturação mais alto para biomassa aérea.

Segmentação e Classificação - determinação espacial das áreas coberta por plantas aquáticas

De uma forma geral, o método de classificação resultou em uma precisão média maior que 95% para os pares de imagens (Radarsat e JERS) sazonais (Tabela 5). As áreas de planta aquática foram corretamente classificadas em mais de 95% dos casos para qualquer estação do ano. A maior confusão na classificação (11%) foi observada entre as classes de floresta inundada e floresta não inundada em novembro (inundada sazonalmente). Possivelmente o alto conteúdo de umidade do solo, devido ao começo da subida de água, resultou em valores de retroespalhamento altos, os quais favorecem a confusão com a classe de floresta inundada. Os mapas temáticos sazonais são apresentados na figura 4.


Cálculo da produção primária líquida

A PPL foi calculada para cada célula de resolução da imagem (pixel) considerando (i) a área classificada de plantas aquáticas (item 3.2.3), (ii) o algoritmo estatístico para determinação de biomassa seca aérea (item 3.2.1), (iii) a porcentagem de biomassa seca submersa na biomassa seca total (item 3.1) e (iv) a produção de biomassa líquida mensal. A estratégia para determinação de cada um destes parâmetros é apresentada a seguir.

(i)Separação de áreas sazonais correspondentes a plantas aquáticas

A partir dos mapas temáticos resultantes da classificação, máscaras sazonais correspondentes a plantas aquáticas foram determinadas e sobrepostas nos pares sazonais de imagens de radar.

(ii)Geração de imagens de biomassa seca aérea (Ba)

Imagens de biomassa seca aérea foram calculadas aplicando-se o modelo estatístico nas imagens sazonais do Radarsat e JERS correspondentes a área coberta por plantas aquáticas (item anterior).

(iii)Geração de imagens de biomassa seca total (Bt)

A biomassa seca total foi calculada para cada pixel da imagem de biomassa seca aérea (item anterior) considerando a porcentagem de biomassa seca submersa na biomassa seca total. Imagens de biomassa total foram geradas para novembro, abril, maio, junho e agosto. O cálculo da biomassa total para cada mês (BT) é o somatório da biomassa de todos os pixels da imagem de biomassa total do respectivo mês. A tabela 6 apresenta os valores de BT calculados para novembro, abril, maio, junho e agosto. Os valores de BT calculados são estatisticamente semelhantes (p<0,05) aos BT calculados por meio de medidas de campo; contudo, observa-se uma diferença de 28% entre estes valores. Note que esta diferença é influenciada principalmente pela maior discrepância entre as BT de novembro (54% de diferença), quando os valores de BT são os mais baixos e, portanto, menos significantes para a estimativa de produção primária líquida total. Desconsiderando-se a discrepância observada para novembro, a diferença total entre a BT medida em campo e a BT estimada com imagens de satélite diminui para 21%.

Os resultados mostram que a BT aumenta continuamente a partir de novembro, quando a nova geração de plantas aquáticas começa a se desenvolver resultando em uma biomassa seca total de 0,1x1012g. As folhas e talos se formam e crescem rapidamente mantendo as estruturas fotossintéticas acima do nível da água. Em maio, no máximo nível da água, a biomassa seca total é 1,45x1012g em uma área de 395km2. A partir de maio o nível da água começa a descer arrastando algumas plantas para o Rio Amazonas. Consequentemente em agosto, a área ocupada pelas plantas e a biomassa seca total decrescem para 281 km2 e 0,95x1012g, respectivamente. Valores similares de biomassa seca total mensal são observados para outras espécies em outras regiões da Amazônia (Piedade et al. 1992).

(iv)Estimativa de biomassa mensal

Como reportado anteriormente, nossos dados foram adquiridos em cinco meses distintos do ciclo de crescimento das plantas, o qual começa em novembro/dezembro, com a subida das águas, maio/junho, com o pico de máximo nível de água, até setembro, quando as plantas estão na fase de senescência. Contudo, como já havia salientado Piedade et al. (1991), estimativas mensais são importantes para estudos de PPL de plantas aquáticas na Amazônia devido ao crescimento contínuo das plantas em um ciclo hidrológico. Portanto, a biomassa total (BT) dos meses nos quais não se obteve dados (campo e satélite) foi estimada estatisticamente a partir da biomassa total calculada com dados de campo e satélite. Utilizou-se a curva de mínimos quadrados determinada com os BT medidos (R2=0,99; n=5; p<0,01) para determinação da biomassa total estimada correspondente aos meses de dezembro, janeiro, fevereiro, março e julho. Este procedimento foi considerado satisfatório porque um contínuo aumento da biomassa de plantas aquáticas é observado para várias regiões da Amazônia (Piedade et al. 1992; Junk & Piedade, 1997). A mesma relação estatística foi utilizada nas imagens de satélite (imagens Bt item anterior) para estimar imagens de Bt para dezembro, janeiro, fevereiro, março e julho. A tabela 6 apresenta os valores de biomassa total (BT) estimados.

(v)Estimativa de PPL

A produção primária líquida é o incremento da biomassa total e biomassa perdida entre dois meses consecutivos para cada pixel da imagem. Caso não haja incremento da biomassa total (diminuição ou constância) entre dois meses consecutivos, a PPL é determinada apenas em função da biomassa perdida entre dois meses consecutivos (Junk & Piedade, 1993). Matematicamente, a PPL é calculada por meio das seguintes equações:

Onde PPLmi é a produtividade mensal, Bti e B(t+1)i são as biomassas totais de dois meses consecutivos, t e (t+1), L(t+1) é a biomassa perdida durante dois meses consecutivos e corresponde a 25% de B(t+1)i (Junk & Piedade, 1993) e "Bti é o incremento na biomassa total entre dois meses consecutivos.

Calculada a PPLmi, a produtividade primária líquida total (PPLi), para cada pixels da imagem, é a soma das PPLmi de cada pixel. A tabela 7 apresenta os valores de PPL média e total por ano estimada a partir de dados de satélite. A partir da imagem de PPLi e da porcentagem de carbono em H. amplexicaules (46% da biomassa total, (Costa, 2000)), calculou-se a imagem de carbono total assimilado em um ano do ciclo hidrológico (Figure 5).


DISCUSSÃO

A tabela 8 mostra valores de máxima biomassa para plantas aquáticas enraizadas da Amazônia reportados por vários autores. De uma forma geral os valores variam de aproximadamente 5000 a 8000 gm-2. Os valores mais altos são para E. polystachya (8000gm-2), uma espécie que suporta condições de alto stress ambiental (Piedade et al. 1991). Valores relativamente mais baixos são reportados para Paspalum fasciculatum e H. amplexicaulis, espécies com um menor porte quando comparadas com E. polystachya. O método baseado em dados de radar estimou uma biomassa máxima de 4500gm-2 para H. amplexicaulis, a qual esta dentro dos valores medidos para outras espécies da Amazônia utilizando intensa amostragem de campo. Esta mesma comparação é valida para a estimativa de PPL pelo método de radar e estimativas baseadas em coletas de campo, como reportado no item anterior.

Os valores de PPL média estimada (6350gm-2yr-1) ou uma PPL anual de 4,1x1012 gyr-1, para uma área máxima de plantas aquáticas de 394km2 (Tabela 7) são semelhantes às estimativas apresentadas para outras espécies de plantas aquáticas na Amazônia. Por exemplo, 9930 gm-2 yr-1 para E. polystachya (Piedade et al. 1991) e 2700 gm-2 em quatro meses para Oryza perennis (Junk & Piedade, 1993). As diferenças em PPL para cada espécie estão possivelmente relacionada com as estratégias necessárias para sobreviverem em um determinado habitat e as características da espécie. Por exemplo, E. polystachya coloniza regiões de alto stress ambiental devido à correnteza e grandes variações de nível de água (8m). Para sobreviver nestas condições esta espécie apresenta rápido crescimento como resultado da uma alta eficiência em converter a radiação solar interceptada em biomassa. Esta e uma característica típica de plantas C4 (Piedade et al., 1991). Por outro lado, áreas mais protegidas e com menor variação de nível de agua (4m), como no caso da região de estudo, são colonizadas por H. amplexicaulis, uma espécie que apesar do menor porte e do aparente modo menos eficiente de fotossintetizar por ser uma planta C3 (Costa, 2000), apresenta valores de PPL altos. Esta alta PPL ocorre parcialmente devido à alta densidade destas plantas (111 plantas m-2, Costa et al., 2002) comparada com a densidade de E. polystachya (30 plantas m-2, Costa et al. 2002) e uma possível adaptação para condições de máxima luz solar resultando em um aumento nas taxas fotossintéticas, o que caracteriza as plantas C4. Juncus roemerianus, uma espécie C3 característica da América do Norte, apresenta este mesmo tipo de adaptação (Mitsch & Gosselink, 2000). Portanto, na planície de inundação da Amazônia, tanto as plantas aquáticas C3 (este estudo) como as C4 (Piedade et al. 1991), de formas distintas, apresentam estratégias para alcançar alta produtividade.

CONCLUSÕES

Imagens dos satélites Radarsat e JERS foram utilizadas eficientemente para o estimar a produtividade primária de plantas aquáticas na planície de inundação da Amazônia. A produção primária é estimada a partir da determinação espacial da área ocupada por plantas aquáticas e da relação estatística entre o "índice de biomassa" () e valores de biomassa amostrados em campo. As áreas de plantas aquáticas foram mapeadas com precisão maior que 95% pelo método de classificação supervisionada por regiões. Este método (segmentação/classificação) foi pela primeira vez utilizado para mapeamento da planície de inundação da Amazônia. A relação estatística para determinação de biomassa aérea apresentou um coeficiente de correlação de 0,82 (p<0,05) e um ponto de saturação de biomassa aérea em 620gm-2. Acima do valor de saturação de biomassa as estimativas do modelo estatístico não são confiáveis. Contudo, a máxima biomassa média aérea observada (agosto, 1996) está abaixo do valor de saturação calculado, o que torna o modelo estatístico confiável para estimativas de biomassa aérea para a região de estudo.

De uma forma geral, a estimativa de PPL pelo método de radar apresenta uma boa concordância com os valores medidos em campo e valores reportados por outros autores, e permitiu estimar uma produção de carbono anual de 2920 g C m-2 yr-1 ou um total de 1,9x1012 g C yr-1 para a área de estudo coberta por plantas aquáticas (máximo em maio - 394km2). Espacialmente (Figure 5), os valores mais baixos de carbono (<900gCm-2) são observados em regiões onde as plantas se desenvolveram somente no começo do ciclo de crescimento. Valores mais altos de carbono (>5000 gCm-2) são observados em regiões próximas ao Rio Amazonas onde a influência de água rica em nutrientes é maior e, portanto uma alta produtividade é observada.

O método baseado em imagens de radar é uma alternativa confiável para estudos de produtividade de plantas aquáticas emergentes em escala local e regional, porque fornece informação com a resolução espacial (pixel) das imagens de satélite comparado com as extrapolações espaciais necessárias quando se utiliza método de amostragem de campo tradicional ou modelo baseado em processos fisiológicos.

AGRADECIMENTOS

A Canadian International Space Agency e Fundação de Amparo a Pesquisa do Estado de São Paulo pelos recursos financeiros para conduzir os trabalhos de campo. Contribuições financeiras do projeto LBA coordenado pelo Dr. John Melack são também reconhecidas. A Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) forneceu a bolsa de estudo para a autora. A Japonesa Space Agency e a Canadian Space Agency forneceram as imagens dos satélites JERS-1 e Radarsat, respectivamente. Gostaria de agradecer aos vários colegas e amigos que ajudaram durante o trabalho de campo, Dra. Novo, Dr. Mantovani e Gilson. Agradeço também a revisão detalhada deste manuscrito por parte de um dos revisores.

BIBLIOGRAFIA CITADA

Recebido em 04/07/2003

Aceito em 23/05/2005

  • Bins, L.S.; Fonseca, L.M.G.; Eethal, G.J. ;Mitsuo, F. 1996. Satellite imagery segmentation: a region growing approach. Proceedings. CD In: Simpósio Brasileiro de Sensoriamento Remoto. Salvador, Bahia, Brasil.
  • Costa, M.P.F. 2000. Net Primary Productivity of Aquatic Vegetation of the Amazon Floodplain: a Multi-SAR Satellite Approach Ph.D. Thesis, University of Victoria, Victoria, Canada. 233p.
  • Costa, M.P.F. 2004. Use of SAR satellite for mapping zonation of vegetation communities in the Amazon floodplain. International Journal of Remote Sensing 25(10):1817-1835.
  • Costa, M.P.F., Niemann, O., Novo, E., Ahern, F. 2002. Biophysical Properties and mapping of aquatic vegetation during the hydrological cycle of the Amazon floodplain using JERS-1 and RADARSAT, 2001. International Journal of Remote Sensing 23(7):1401-1260.
  • Downing, J.A.; Anderson, M.R. 1985. Estimating the standing biomass of aquatic macrophytes. Canadian Journal of Fisheries and Aquatic Science, 42:1860-1869.
  • Junk, W.J. 1993. Wetlands of Tropical South America. In: Whigham D.F., Hejny, S., Dykyjova, D. (Eds). Wetlands of the world IC. Kluwer, Dordrecht, pp: 679-739.
  • Junk, W.J.; Piedade, M.T. 1993. Biomass and primary productivity of herbaceous plant communities in the Amazon floodplain. Hydrobiology, 263: 155-162.
  • Junk, W.J.; Piedade, M.T. 1997. Plant life in the floodplain with special reference to herbaceous plants. In: W.J. Junk (Eds.). The Central Amazon Floodplain: Ecology of a Pulsing System. Springer. Berlin. pp.147-181.
  • Laur, H., Bally, P., Meadows, P., Sanchez, J., Schaettler, B.,; Lopinto, E. 1996. Derivation of the backscattering coefficient sş in ERS.SAR.PRI products. ESA/ESRIN, Issue 2, Rev. 2. 55p.
  • Le Toan, T., Ribbes, F., Wang, L., Floury, N., Ding, K.; Kong, J.A.,; Fujita, M. 1997. Rice crop mapping and monitoring using ERS-1 data based on experiment and modeling results. IEEE Transaction of Geoscience and Remote Sensing, 35(1): 41-55.
  • Melack, J.M.; Forsberg, B.R. 2000. Biogeochemistry of Amazon floodplain lakes and associated wetlands. In: M.E. McClain, R.L. Victoria, J.E. Richey (Eds.) Biogeochemistry of the Amazon Basin Oxford University Press. Oxford. pp.210-260.
  • Mitsch, J.W.; Gosselink, J.G. 2000. Wetlands John Wiley & Sons, Inc, New York, USA. 920p.
  • Neter, J.; Wasserman, W.; Kutner, M.H. 1990. Applied Linear Statistical Models IRWIN. USA. 1181p.
  • Novo, E.M.L.M.; Leite, F.A.; Avila, J.; Ballester, V.; Melack, J.M. 1997. Assessment of Amazon floodplain habitats using TM/Landsat data. Ciência e Cultura 49(4):280-284.
  • Piedade, M.T.F.; Junk, W.J.; Long, S.P. 1991. The productivity of the C4 grass Echinochloa polystachya on the Amazon floodplain. Ecology 72(4), 1456-1463.
  • Piedade, M.T.F.; Junk, W.J.; Mello, J.A.N. 1992. A floodplain grassland on the central Amazon. In: B.J.Jones and M.J. Roberts (Eds). Primary Productivity of Grass Ecosystems of the Tropics and Sub-tropics Chapman & Hall/UNEP. New York. pp.127-158.
  • Potter, C.S.; Randerson, J.T.; Field, C.B.; Matsaon, P.A.; Vitousek, P.M.; Mooney, H.A.; Klooster, S.A. 1993. Terrestrial ecosystem production: a process model based on global satellite and surface data. Global Biogeochemical cycles 7(4):811-841.
  • RADAMBRASIL. 1976. Folha SA.21-Santarém. Levantamento de Recursos Naturais. Vol.10. Rio de Janeiro, Brasil. 200p.
  • Ribbes, F.; Le Toan, T. 1999. Rice field mapping and monitoring with RADARSAT data. International Journal of Remote Sensing, 20(4), 745-765.
  • Richards, J.A. 1986. Remote Sensing Digital Analysis Springer-Verlag. Berlin. 281p.
  • Sellers, P.J.; Meeson, B.W.; Hall, F.G.; Asrar, G; Murphy, R.E.; Schiffer, R.A; Bretherton, F.P.; Dickinson, R.E.; Ellingson, R.G.; Field, C.B.; Huemmrich, K.F.; Justice, C.O.; Melack, J.M.; Roulet, N.T.; Schimel, D.S.; & Try, P.D. 1995. Remote Sensing of the land surface for studies of global change: models-algorithms-experiments. Remote Sensing of Environment 51(1):3-26.
  • Shepherd, N. 1997. Output scaling look up tables in RADARSAT CDPF products. In: Proceedings of a Workshop on RADARSAT Data Quality Quebec, Canada. 20p.
  • Shimada, M. 1996. Radiometric and geometric calibration of JERS-1 SAR. Advances in Space Research, 17: 79-88.
  • Snedecor, G.W.; Cochran, W.G. 1980. Statistical Methods Iowa State University Press. Iowa, USA. 507p.

Datas de Publicação

  • Publicação nesta coleção
    17 Out 2005
  • Data do Fascículo
    Jun 2005

Histórico

  • Aceito
    23 Maio 2005
  • Recebido
    04 Jul 2003
Instituto Nacional de Pesquisas da Amazônia Av. André Araujo, 2936 Aleixo, 69060-001 Manaus AM Brasil, Tel.: +55 92 3643-3030, Fax: +55 92 643-3223 - Manaus - AM - Brazil
E-mail: acta@inpa.gov.br