A reproducible algorithm and R implementation for classifying ecological flood pulse phases from flood-area time series

Resumo Técnico: Um Algoritmo Reproduzível e Implementação em R para Classificação de Fases de Pulso de Inundação Ecológica

Definição do Problema
Estudos ecológicos de ecossistemas de rios e planícies de inundação dependem do Conceito de Pulso de Inundação, que postula que a expansão e contração hidrológicas periódicas regulam a disponibilidade de habitat, a troca de nutrientes e os processos biológicos. No entanto, existe uma lacuna metodológica na pesquisa empírica: os dados hidrológicos são tipicamente trajetórias contínuas, enquanto as observações biológicas são frequentemente agregadas em categorias discretas e arbitrárias (por exemplo, "estação úmida", "nível alto" ou meses específicos). Essas categorias simples falham em capturar a direção funcional de um pulso de inundação. Uma inundação de uma magnitude específica pode ter implicações ecológicas distintas dependendo se o sistema está em expansão, estável ou em recessão. Há uma carência de um método reproduzível e transparente para converter séries temporais hidrológicas contínuas em "Fases de Pulso de Inundação Ecológica" (FPPI) específicas que considerem tanto a magnitude quanto a trajetória.

Metodologia
O artigo introduz o Índice de Discriminação de Fase de Pulso de Inundação (FPDI), um fluxo de trabalho projetado para traduzir séries temporais contínuas de magnitude de pulso de inundação em cinco categorias distintas de FPPI:

1. Água Baixa (Low Water)
2. Expansão de Inundação (Flood Expansion)
3. Acúmulo de Inundação (Flood Build-up)
4. Recessão de Inundação (Flood Recession)
5. Contração de Inundação (Flood Contraction)

O algoritmo opera através das seguintes etapas:

• Padronização de Entrada: A entrada é uma série temporal cronologicamente ordenada de magnitude de pulso de inundação (ex: área de inundação, nível da água). O sistema realiza uma auditoria de valores ausentes, duplicatas e intervalos de datas.
• Cálculo da Direção Hidrográfica: A mudança entre observações consecutivas ($\Delta_t = X_t - X_{t-1}$) determina a direção. Um limiar definido pelo usuário (padrão de 0,50 km²) classifica as observações como Ascendente ($\Delta > 0,50$), Descendente ($\Delta < -0,50$) ou Estável. A primeira observação é classificada como Inicial.
• Normalização Anual: Para considerar a variabilidade interanual na magnitude absoluta da inundação, os valores de inundação são normalizados dentro de cada ano civil: $FloodNorm_t = (X_t - min(X_y)) / (max(X_y) - min(X_y))$.
• Atribuição de Fase: Os rótulos de FPPI são atribuídos combinando a direção hidrográfica com a magnitude da inundação normalizada usando limiares definidos pelo usuário (limiar padrão de água baixa = 0,35; limiar de inundação alta = 0,60).
  • Água Baixa: Estável e abaixo do limiar de água baixa.
  • Expansão de Inundação: Ascendente e abaixo do limiar de inundação alta.
  • Acúmulo de Inundação: Ascendente e no nível ou acima do limiar de inundação alta.
  • Recessão de Inundação: Descendente e no nível ou acima do limiar de inundação alta.
  • Contração de Inundação: Descendente e abaixo do limiar de inundação alta.
• Consolidação Mensal: Para compatibilidade com amostragens biológicas, a série ao nível de observação é consolidada em um único rótulo por mês civil, retendo a última observação dentro daquele mês.

Implementação e Validação
O fluxo de trabalho é implementado como um pacote R leve e puro chamado fpdi. Ele não requer código compilado e inclui funções para auditoria de dados, classificação e validação hidrológica.

• Estudo de Caso: O método foi demonstrado usando uma série temporal de área de inundação (2017–2020) para o Lago Yardi, Etiópia. As áreas de inundação foram derivadas de imagens Sentinel-2 usando o Índice de Água de Diferença Normalizada Modificado (MNDWI).
• Validação Hidrológica: A ANOVA de uma via confirmou que as fases atribuídas separaram significativamente os dados hidrológicos ($F = 19,71, p = 3,37 \times 10^{-9}$), com as fases explicando aproximadamente 65% da variação da área de inundação (eta-quadrado = 0,65).
• Análise de Sensibilidade: O sistema foi testado através de uma grade de limiares de direção, água baixa e inundação alta. Os resultados indicaram que a estrutura de fase foi robusta; a variação dos limiares dentro dos intervalos testados resultou em mudanças mínimas nos rótulos (0,0% a 4,65% para rótulos mensais), e a distinção hidrológica das fases permaneceu substancial.

Principais Contribuições

1. Ponte Operacional: O FPDI fornece um método reproduzível e transparente para ligar as dinâmicas contínuas de pulso de inundação com unidades analíticas ecológicas discretas, sem propor uma nova teoria ecológica.
2. Sensibilidade Direcional: Ao contrário das categorias estáticas, o FPPI distingue entre estados de magnitude semelhante, mas trajetórias opostas (ex: Acúmulo de Inundação vs. Recessão de Inundação), o que é crítico para entender processos ecológicos como desova ou migração.
3. Reprodutibilidade: O pacote R garante que a lógica de classificação seja auditável, permitindo que pesquisadores inspecionem limiares, regras de consolidação e transformações de dados de entrada.
4. Modularidade: O pacote separa a lógica central de classificação do processamento de sensoriamento remoto a montante, tornando-o acessível a ecólogos que já possuem séries temporais hidrológicas, permitindo também a futura integração com ferramentas de dados de satélite.

Significância e Alegações
O autor afirma explicitamente que o FPDI não é uma alternativa ao Conceito de Pulso de Inundação, mas sim um "método operacional aninhado" para aplicar a teoria a sistemas específicos, particularmente pequenos lagos conectados a planícies de inundação e sistemas de oxbow (braços mortos).

• O método é projetado para tornar as dinâmicas locais de inundação "auditáveis, reproduzíveis e analiticamente utilizáveis" na pesquisa ecológica e de pesca.
• Ele aborda a limitação das categorias baseadas no calendário ou simples categorias de nível alto/baixo ao definir fases baseadas no estado funcional real do sistema.
• O autor reconhece que, embora o FPDI identifique a "oportunidade de fase" hidrológica, ele não garante respostas biológicas específicas, que dependem de histórias de vida de espécies específicas e outros fatores ambientais.
• O artigo conclui que a demonstração no Lago Yardi produziu com sucesso fases interpretáveis e hidrologicamente distintas, além de uma implementação de software robusta, validando a utilidade do fluxo de trabalho para transferir a teoria do pulso de inundação para a análise prática de dados.