Spatiotemporal Heterogeneity of Daily Rainfall Persistence across the Maluku Archipelago: A First-Order Markov Chain Analysis of ERA5 Reanalysis Data (2000–2025)

Resumo Técnico: Heterogeneidade Espaciotemporal da Persistência de Precipitação Diária no Arquipélago de Maluku

Problema
Arquipélagos tropicais, particularmente o Arquipélago de Maluku dentro do Continente Marítimo Indonésio, exibem uma extrema heterogeneidade hidroclimática impulsionada pela interação entre a complexa topografia insular, mares quentes circundantes e a circulação atmosférica de múltiplas escalas. Embora a climatologia média de precipitação e os índices de precipitação extrema tenham sido extensivamente estudados nesta região, a persistência temporal da precipitação diária — a tendência de condições úmidas ou secas continuarem de um dia para o outro — permanece mal caracterizada. Esta lacuna é crítica para a gestão de riscos hidrológicos em sistemas insulares pequenos e com limitação de recursos hídricos, onde a duração de períodos secos ou úmidos influencia diretamente a suscetibilidade à seca e o potencial de inundações repentinas. Além disso, a escassez de observações terrestres contínuas e de longo prazo na região tem dificultado a caracterização sistemática dessas dinâmicas.

Metodologia
Este estudo aborda a lacuna de dados acoplando dados de reanálise ERA5 de alta resolução (2000–2025) a uma estrutura de cadeia de Markov de primeira ordem de dois estados. A análise foca em oito zonas climáticas representativas que abrangem diversos cenários fisiográficos: interiores montanhosos (ex: Ilha de Seram), planícies costeiras, ilhas de oceano aberto e a fronteira sul seca (Ilhas de Tanimbar).

As principais etapas metodológicas incluíram:

1. Processamento de Dados: Totais de precipitação diária foram extraídos de células de grade ERA5 próximas a locais de referência. Um limiar fixo de 1 mm dia⁻¹ foi utilizado para classificar os dias em estados discretos "secos" (estado 0) e "úmidos" (estado 1).
2. Modelagem de Markov: Uma cadeia de Markov de primeira ordem foi aplicada para calcular matrizes de probabilidade de transição mensal ($P_{DD}$, $P_{DW}$, $P_{WD}$, $P_{WW}$). Estas matrizes quantificam a probabilidade de transição entre estados secos e úmidos.
3. Métricas Derivadas: A partir das probabilidades de transição, o estudo derivou:
   • Probabilidades estacionárias ($\pi_0$, $\pi_1$) representando a probabilidade de longo prazo de estados secos ou úmidos.
   • Durações esperadas de episódios ($E(L_d)$, $E(L_w)$) calculadas como o inverso das probabilidades de transição ($1/P_{DW}$ e $1/P_{WD}$).
   • Durações de episódios extremos definidas pelo percentil 90 ($P_{90}$) dos comprimentos de episódios contínuos.
   • Limiares de precipitação extrema ($P_{90}$, $P_{95}$) baseados na distribuição da precipitação em dias úmidos.
4. Validação: A validade do modelo foi avaliada comparando o ciclo anual reconstruído de precipitação (via simulação estocástica) contra o ciclo anual observado do ERA5, confirmando que a estrutura captura adequadamente a evolução sazonal.

Principais Resultados
A análise revela uma pronunciada heterogeneidade espacial na persistência da precipitação através do arquipélago, governada pela interação entre a orografia e a circulação monçônica:

• Gradientes Espaciais:
  • Interiores Montanhosos (Zona 3, Seram): Registraram a maior precipitação diária média (11,85 mm dia⁻¹) e o maior limiar de precipitação extrema ($P_{95} = 72,8$ mm). Esta zona exibiu a mais forte persistência úmido-úmido, consistente com o reforço orográfico da chuva convectiva.
  • Ilhas do Sul (Zona 8, Saumlaki): Registraram a menor precipitação diária média (4,45 mm dia⁻¹) e uma mediana próxima a zero. Esta zona mostrou a maior persistência seco-seco, impulsionada pela intrusão sazonal de ar continental seco do norte da Austrália durante a Monção de Sudeste.
• Dinâmicas Sazonais:
  • Episódios Úmidos: As durações esperadas de episódios úmidos atingiram o pico em junho (média de 8,87 dias; $P_{90} = 21,3$ dias), coincidindo com o transporte ativo de umidade monçônica.
  • Episódios Secos: As durações esperadas de episódios secos atingiram o pico em julho (média de 6,72 dias; $P_{90} = 14,8$ dias), correspondendo ao pico da estação seca.
  • Variabilidade: A variabilidade mensal na precipitação foi menor em janeiro e abril, mas aumentou marcadamente durante o meio e a parte final do ano, atingindo um máximo primário em julho e um máximo secundário em dezembro. Estes meses de elevada variabilidade indicam maior flutuação interanual e redução da previsibilidade durante o pico da estação seca e o início da estação úmida, consistentes com a maior instabilidade na transição do monção regional.

Significância e Alegações
O artigo alega que a persistência da precipitação no Arquipélago de Maluku não é um fenômeno aleatório ou específico de uma ilha, mas está organizada ao longo de um gradiente espacial coerente e fisicamente explicável. O estudo demonstra que esta persistência é regida conjuntamente pela circulação monçônica e pela complexidade fisiográfica em escala de ilha, em vez de apenas um desses fatores agindo isoladamente.

A principal significância do trabalho reside na sua contribuição metodológica: estabelece uma abordagem transferível e reproduzível para quantificar a memória atmosférica em regiões de ilhas tropicais com escassez de dados. Ao utilizar uma estrutura de Markov de primeira ordem parcimoniosa aplicada a dados de reanálise disponíveis gratuitamente, o estudo fornece uma alternativa viável às descrições climatológicas convencionais onde as redes terrestres são insuficientes.

Implicações Práticas
Os achados oferecem relevância direta para:

• Monitoramento de Secas: Identificação das ilhas do sul como zonas de alto risco para episódios secos prolongados e escassez de água.
• Alerta Precoce de Inundações Repentinas: Destaque para as bacias montanhosas como áreas propensas a episódios úmidos estendidos e perigos de precipitação cumulativa.
• Adaptação Climática: Argumentação de que uma estratégia de adaptação única para todo o arquipélago é inadequada devido à natureza assimétrica dos perigos (persistência de episódios úmidos nas montanhas vs. persistência de episódios secos no sul).

Os autores mantêm uma postura modesta, reconhecendo que a resolução de 0,25° do ERA5 pode não resolver circulações locais de pequena escala (ex: brisas marítimas) e que a suposição de Markov de primeira ordem simplifica a memória atmosférica. No entanto, eles afirmam que a estrutura recupera com sucesso sinais de persistência espacialmente resolvidos e fisicamente coerentes que, de outra forma, seriam inacessíveis nesta região.