Beta — piloto fechado · dados podem conter inconsistências
🌊 Como tá o Mar?

Metodologia científica

Descrição técnica do sistema de previsão de turbidez costeira.

Conteúdo

  1. Fundamentação física
  2. Dados utilizados
  3. Algoritmo de turbidez Sentinel-2
  4. Modelo preditivo
  5. Índice de limpeza
  6. Limitações conhecidas
  7. Referências

1. Fundamentação física

A turbidez das praias de Marataízes é controlada por dois mecanismos principais:

O gradiente espacial observado — maior turbidez na Praia da Colônia (4 km da foz) do que na Restinga (7,5 km) — reflete a atenuação da pluma sedimentar com a distância à desembocadura.

2. Dados utilizados

Fonte Variáveis Período
Sentinel-2 (ESA) Turbidez observada (FNU), cobertura de nuvens 2015–presente, ~5 dias
Open-Meteo Chuva 24h/72h/7d, temperatura ar, pressão 2015–presente, diário
ERA5 (ECMWF) Altura/período/direção de ondas, SST, vento 2015–2024, horário
ANA (Hidroweb) Chuva na bacia hidrográfica (Cachoeiro de Itapemirim) 2015–2024, diário

Em produção (pós-2024), ERA5 é substituído por previsão GFS/Open-Meteo Marine API para dados de ondas em tempo quase-real.

3. Algoritmo de turbidez Sentinel-2

A turbidez é estimada a partir da reflectância de superfície na banda vermelha (B4, 665 nm) e, para turbidez muito alta, na banda NIR (B5, 705 nm), seguindo o algoritmo de comutação de Dogliotti et al. (2015):

se ρ_B4 < 0.05:

T = A_low × ρ_B4 / (1 − ρ_B4 / C_low) # Banda vermelha

senão se ρ_B4 > 0.07:

T = A_high × ρ_B5 / (1 − ρ_B5 / C_high) # Banda NIR

senão:

T = interpolação linear entre as duas estimativas

Onde A_low = 228,1, C_low = 0,1641, A_high = 3078,9, C_high = 0,2112 (calibração original para águas costeiras europeias, aplicada como referência).

Pixels com cobertura de nuvens, sombra de nuvens ou adjacentes à costa (< 3 pixels de distância) são mascarados. A turbidez final por praia é a mediana dos pixels válidos na zona de batimento de onda (buffer de 300 m centrado no ponto de referência).

4. Modelo preditivo

O modelo é um ensemble de dois estimadores LightGBM treinados conjuntamente:

Regressor

Prevê turbidez contínua em FNU. Otimizado por RMSE.

Classificador

Prevê classe {limpa, media, barrenta}. Usado para probabilidade.

Features principais (22 no total): chuva acumulada em múltiplas janelas (24h, 72h, 7d), chuva na bacia hidrográfica com defasagem, velocidade e direção do vento, componente meridional do vento, altura/período/direção de ondas, SST, temperatura do ar, maré (amplitude e fase), índice de sizígia, e distância à foz.

Validação: split temporal estrito. Treino 2017–2024 (n = 1.371 amostras), validação 2025–2026 (n = 197 amostras). Métricas: R² = 0,59, MAE = 8,3 FNU, AUC = 0,86 para classe limpa.

Horizonte de previsão: D+0 a D+6. A incerteza cresce com o horizonte porque depende de previsão meteorológica, não de observação.

5. Índice de limpeza (0–100)

A turbidez em FNU é mapeada para um índice intuitivo de 0 (barrento) a 100 (cristalino) via função sigmóide inversa calibrada para Marataízes:

Índice = 100 / (1 + exp(0,10 × (FNU − 24,5)))

Parâmetros calibrados (k = 0,10; midpoint = 24,5 FNU) por minimização do erro quadrático em cinco âncoras derivadas da distribuição empírica de turbidez observada em Marataízes (2015–2024) e da sensibilidade humana à transparência da água:

FNU Índice Significado
5~88Limpa
10~81Fronteira limpa/boa
18~66Fronteira boa/turva
30~37Fronteira turva/barrenta
60~3Barrenta forte

6. Limitações conhecidas

7. Referências