Resumen
| Humedal | Última fecha válida | Agua (ha) | NDCI | NDTI | Alertas |
|---|
| Doñana | 2026-10-12 | 7.815 | -0.163 | -0.449 | ninguna |
| Mar Menor | 2026-09-21 | 13.210 | -0.144 | -0.430 | ninguna |
| Tablas de Daimiel | 2026-09-24 | 189 | 0.156 | -0.011 | ninguna |
| L'Albufera de València | 2026-10-08 | 9.948 | 0.050 | -0.675 | ninguna |
| Laguna de Fuente de Piedra | 2026-08-18 | 330 | 0.052 | 0.034 | ninguna |
| Laguna de Gallocanta | 2026-10-16 | 434 | 0.502 | -0.431 | eutrofización |
Doñana (Andalucía · Natura 2000 ES0000024)
Marismas del Guadalquivir. Inundación estacional; desecación por sobreexplotación del acuífero.
Última fecha válida2026-10-12
Agua7.815 ha17 % del área inundable
NDCI medio-0.163
NDTI medio-0.449
Nubes en el sitio2 %
Observaciones válidas354de 609 fechas
Sin alertas: valores dentro del rango de referencia.
Área inundable 46.709 ha de las 128.265 ha del sitio Natura 2000, de las cuales 7.688 ha con agua casi siempre; medida acumulando 35 fechas de meses húmedos. Calado y lluvia del panel inferior: ICTS-Doñana (Estación Biológica de Doñana, CSIC) — CC BY 4.0.
Mar Menor (Murcia · Natura 2000 ES6200030)
Laguna hipersalina. Crisis de eutrofización recurrentes (2016, 2019, 2021).
Última fecha válida2026-09-21
Agua13.210 ha98 % del área inundable
NDCI medio-0.144
NDTI medio-0.430
Nubes en el sitio0 %
Observaciones válidas373de 562 fechas
Sin alertas: valores dentro del rango de referencia.
Área inundable 13.459 ha de las 13.583 ha del sitio Natura 2000, de las cuales 13.062 ha con agua casi siempre; medida acumulando 33 fechas de meses húmedos.
Tablas de Daimiel (Castilla-La Mancha · Natura 2000 ES0000013)
Humedal fluvial dependiente del acuífero 23. Se seca casi por completo en sequías.
Última fecha válida2026-09-24
Agua189 ha40 % del área inundable
NDCI medio0.156
NDTI medio-0.011
Nubes en el sitio0 %
Observaciones válidas266de 486 fechas
Sin alertas: valores dentro del rango de referencia.
Área inundable 475 ha de las 3.100 ha del sitio Natura 2000, de las cuales 30 ha con agua casi siempre; medida acumulando 24 fechas de meses húmedos.
L'Albufera de València (C. Valenciana · Natura 2000 ES0000023)
Laguna costera eutrófica rodeada de arrozal.
Última fecha válida2026-10-08
Agua9.948 ha49 % del área inundable
NDCI medio0.050
NDTI medio-0.675
Nubes en el sitio0 %
Observaciones válidas359de 502 fechas
Sin alertas: valores dentro del rango de referencia.
Área inundable 20.192 ha de las 27.538 ha del sitio Natura 2000, de las cuales 8.927 ha con agua casi siempre; medida acumulando 38 fechas de meses húmedos.
Laguna de Fuente de Piedra (Andalucía · Natura 2000 ES0000033)
Laguna salina endorreica. Colonia de flamenco; se seca en verano.
Última fecha válida2026-08-18
Agua330 ha26 % del área inundable
NDCI medio0.052
NDTI medio0.034
Nubes en el sitio6 %
Observaciones válidas427de 978 fechas
Sin alertas: valores dentro del rango de referencia.
Área inundable 1.284 ha de las 8.668 ha del sitio Natura 2000, de las cuales 649 ha con agua casi siempre; medida acumulando 34 fechas de meses húmedos.
Laguna de Gallocanta (Aragón · Natura 2000 ES2430043)
Laguna salina endorreica. Invernada de grulla; ciclos plurianuales de inundación.
Última fecha válida2026-10-16
Agua434 ha36 % del área inundable
NDCI medio0.502
NDTI medio-0.431
Nubes en el sitio3 %
Observaciones válidas513de 904 fechas
EUTROFIZACIÓN · alta
Pico de NDCI 0.513 en los últimos 20 días, por encima de los picos habituales en estas fechas (0.431 es el percentil 90 de 56 ventanas de años anteriores), con el 92 % de la lámina por encima de 0.2. Supera además el umbral de floración algal de la literatura.
Área inundable 1.206 ha de las 2.814 ha del sitio Natura 2000, de las cuales 688 ha con agua casi siempre; medida acumulando 34 fechas de meses húmedos.
Mapa
Cómo se calcula
- Imagen. Cada pocos días se descargan las imágenes Sentinel-2 L2A del catálogo público Earth Search, recortadas al polígono Natura 2000 de cada humedal y pasadas a metros (EPSG:25830) con píxel de 20 m; Doñana a 40 m.
- Fechas válidas. Se tira la fecha si queda nublado más del 20 % del humedal, si las escenas cubren menos del 95 %, si hay neblina o si el espectro del agua delata una corrección atmosférica fallida.
- Agua. Un píxel es agua si su MNDWI pasa del corte calculado ese día con el método de Otsu (o de 0 si el histograma no lo permite). Es agua libre si además su NDVI es menor que 0,15, y vegetación inundada si no. Contando píxeles salen las hectáreas.
- Turbidez y clorofila. NDTI y NDCI se promedian solo sobre el agua libre.
- Valor de ahora. La mediana de las fechas válidas de los últimos 20 días, si hay al menos 2; si no, la última.
- Comparación histórica. Se compara con las fechas de años anteriores a ±30 días del mismo día del año, desde julio de 2017, con un mínimo de 5 observaciones de referencia.
- Alerta. Desecación si la lámina queda por debajo del percentil 10 de esa época; descenso brusco si cae por debajo del 70 % de la mediana de los 45 días previos y además más que el 90 % de las caídas de esa época; eutrofización si el pico de NDCI supera el percentil 90 de los picos de la época más 0,02; turbidez si el NDTI supera su percentil 90.
Parámetros
| Parámetro | Valor | Uso |
|---|
| Píxel de trabajo | 20 m | Doñana a 40 m. |
| Nubes en el humedal | ≤ 20 % | Si hay más, la fecha se tira. |
| Cobertura de las escenas | ≥ 95 % | Si hay menos, la fecha se tira. |
| Azul de nube | > 0,22 | Reflectancia azul a partir de la cual un píxel es nube. |
| Azul de neblina | > 0,12 | Azul mediano del humedal; la fecha se tira. |
| Infrarrojo / verde sobre el agua | > 1,1 | Corrección atmosférica fallida; la fecha se tira. |
| Otsu: píxeles válidos | ≥ 500 | Para poder calcular el corte del día. |
| Otsu: separabilidad | ≥ 0,6 | Si no, no hay dos grupos que separar. |
| Otsu: banda del corte | -0,35 a 0,35 | Fuera de ella se usa el corte de reserva. |
| Corte de reserva de MNDWI | 0 | Cuando falla alguna de las tres condiciones. |
| NDVI del agua libre | < 0,15 | Por encima, vegetación inundada. |
| Valor actual | 20 días | Mediana de las fechas válidas. |
| Ventana estacional | ±30 días | En años anteriores. |
| Referencia mínima | 5 obs. | Sin ellas no hay alerta relativa. |
| Desecación | percentil 10 | De la lámina en esa época. |
| Descenso brusco | < 70 % | De la mediana de los 45 días previos, y peor que el 90 % de las caídas de la época. |
| Eutrofización | percentil 90 + 0,02 | Del pico de NDCI; alta si además pasa de 0,2. |
| Turbidez | percentil 90 | Del NDTI en esa época. |
Validación
El área inundada que mide el satélite en Doñana (agua libre más vegetación inundada) se ha comparado con el calado medido en el suelo de la marisma por la ICTS-Doñana: la correlación es de 0,94 (Pearson; 0,78 de Spearman) sobre las 168 fechas de 2022 a 2026 en que hay las dos medidas, y el área sube con cada tramo de calado. Las hectáreas absolutas también se sostienen: la lámina permanente del Mar Menor sale entre 12.950 y 13.340 ha de mediana anual en nueve años, frente a las 13.500 ha que se citan habitualmente para la laguna.
En detalle: Qué se mide, y con qué · El umbral de agua se calcula en cada fecha · Qué fechas se tiran, y por qué · Con qué se compara la superficie medida · Cómo se decide que algo es una alerta · El contraste con medidas de campo · Lo que este informe no puede decir · Glosario de siglas
Qué se mide, y con qué
Cada pocos días un satélite pasa por encima del humedal y mide cuánta luz devuelve el suelo en
varios colores. El agua tiene una firma muy reconocible: refleja algo en el verde y absorbe casi
todo el infrarrojo. Comparando esas dos bandas píxel a píxel se sabe qué parte del humedal está
inundada, y contando píxeles eso se convierte en hectáreas.
De cada fecha válida salen tres medidas:
- Agua libre: píxeles con MNDWI por encima del umbral del día y NDVI <
0.15. La lámina abierta, sin plantas emergiendo.
- Vegetación inundada: el mismo umbral de agua pero con NDVI por encima. Carrizo,
masiega, arrozal: agua que no se ve como espejo pero cuenta hidrológicamente. Es lo que en
Doñana correlaciona con el calado medido en el suelo.
- Turbidez (NDTI) y clorofila (NDCI), promediadas solo sobre el agua libre: sobre
vegetación no significarían nada.
Todo se recorta con el polígono oficial Natura 2000 del humedal, reproyectado a
EPSG:25830 para trabajar en metros, y se muestrea a 20 m de píxel;
Doñana a 40 m, porque necesita 4,5 escenas por fecha y a 20 m tardaba cuatro veces más sin cambiar el resultado.
El umbral de agua se calcula en cada fecha
La receta de manual es «hay agua donde MNDWI > 0». Aquí no se usa, porque medía
sistemáticamente un 15 % menos de agua que la propia clasificación de la ESA. El motivo es
físico: en una orilla somera y turbia el fondo se ve a través del agua, y el índice se queda por
debajo de cero sin dejar de ser agua. Un corte fijo, además, no puede valer igual en una laguna
salina de fondo claro y en una marisma cargada de materia orgánica.
Así que el corte se calcula en cada fecha y para cada humedal con el método de
Otsu sobre el histograma de ese día: se busca el valor que mejor separa las dos poblaciones
—agua y no agua— que la imagen contiene. Con tres cautelas, porque un método automático también se
equivoca:
- hacen falta al menos 500 píxeles válidos para que el histograma
signifique algo;
- la separabilidad (cuánta de la varianza total explica la partición) debe llegar a
0.6: por debajo de eso no hay dos grupos, hay uno, y partirlo
sería inventarse una orilla;
- el corte tiene que caer entre -0.35 y 0.35; fuera de esa
banda no puede ser una orilla.
Si alguna de las tres falla, esa fecha usa el umbral de reserva 0.0 en lugar de
un valor absurdo.
Qué fechas se tiran, y por qué
De las 4.041 fechas descargadas solo 2.292 —un 57 %— llegan al informe. Ese filtro es la mitad del trabajo: una nube fina o una corrección atmosférica fallida no dan un error, dan un
número plausible y equivocado, que es mucho peor.
- Nubes: las clases de nube de la SCL, más todo píxel con azul por encima de
0.22. Si queda nublado más del 20 % del humedal, fuera.
- Cobertura: las escenas del día tienen que cubrir al menos el
95 % del humedal. Medir media laguna y compararla con la laguna entera
sería un descenso brusco de mentira.
- Neblina: si el azul mediano del sitio pasa de 0.12, fuera. Las calimas
finas no las marca la SCL y llegaban a triplicar el agua detectada.
- Desacuerdo con la ESA: si la SCL ve agua de sobra y el índice detecta menos del
25 % de esa cantidad, algo va mal en esa escena.
- Forma del espectro sobre el agua: sobre agua real el infrarrojo cercano queda por
debajo del verde. Si la razón infrarrojo/verde pasa de 1.1, o el
azul se desploma por debajo de 0.002, lo que falló fue la corrección
atmosférica, no la laguna.
Ese último control resolvió un falso culpable. El satélite más nuevo, Sentinel-2C, parecía tener
un sesgo propio; en realidad el problema eran correcciones atmosféricas fallidas, que afectan a los
tres satélites y solo son mucho más frecuentes en él. Detectar la causa permitió seguir usándolo en
vez de vetarlo.
Con qué se compara la superficie medida
Decir «el humedal está al 16 % de su superficie» invita a una pregunta: ¿el 16 % de qué? El
polígono Natura 2000 es un límite administrativo. El de Doñana son 128.265 ha que incluyen
terreno que no se inunda jamás, así que medir la lámina contra ese total no dice
nada hidrológico.
Por eso el denominador de este informe es el área inundable: la parte que ha llegado a
tener agua alguna vez en nueve años, medida acumulando las fechas válidas de los meses húmedos.
En Doñana son 46.709 ha, 2,7 veces menos que el polígono. El método se valida a sí mismo: al Mar Menor
le sale un 99 % de área inundable y un 96 % con agua permanente, que es exactamente lo que debe
salir en una laguna costera.
Cómo se decide que algo es una alerta
Ninguna alerta compara con un valor absoluto ni con otro humedal. Todas comparan cada humedal
consigo mismo en la misma época del año: se toma la ventana de ±30 días
alrededor de la fecha en los años anteriores, y hacen falta al menos 5
observaciones para que exista referencia. El «valor de ahora» es la mediana de las observaciones
válidas de los últimos 20 días, no la última, que podría ser un mal día.
- Desecación (grave): la lámina actual por debajo del percentil 10 de lo que ese humedal
da en estas fechas. No se aplica a las lagunas de agua permanente.
- Descenso brusco: caída por debajo del 70 % de la mediana de los
45 días anteriores y además peor que el
90 % de las caídas registradas en esa misma época en años
anteriores. Sin esa segunda condición la regla avisaba de que las lagunas se secan en verano,
que es lo que llevan haciendo desde siempre: en Gallocanta saltaba en 100 de sus 494 fechas.
- Eutrofización: el pico de NDCI de la ventana actual por encima del percentil
90 % de los picos de esa época en años anteriores, más un margen de
0.02. Pico y no mediana, porque una floración de algas dura días y una mediana de
tres semanas la borra.
- Turbidez: NDTI actual por encima del percentil 90 de la época.
El umbral de floración algal que da la literatura (NDCI > 0.2, unos 40 mg/m³
de clorofila-a) no sirve para disparar en estas lagunas: se calibró en aguas continentales
profundas, y en lagunas someras, salinas y de fondo claro el índice está inflado de forma crónica
—Gallocanta lo supera en las 494 fechas que tiene—. Aquí solo se usa para graduar la
gravedad: una alerta que además pasa ese umbral es alta; si no, media.
El contraste con medidas de campo
El satélite dice cuánta superficie hay inundada, no por qué. Donde existe una red de medida en el
suelo se añade al pie de la gráfica, y en Doñana la hay: ICTS-Doñana (Estación Biológica de Doñana, CSIC) — CC BY 4.0, con cinco estaciones que
publican calado dentro de la marisma y siete que miden lluvia.
Sirve de dos maneras. Como explicación, porque la lluvia acumulada de cada mes va detrás
del calado y se ve qué crecidas vienen de qué precipitación. Y como validación independiente:
el área inundada que mide el satélite y el calado que miden los sensores en el suelo correlacionan
0,94 sobre las 168 fechas que comparten, y la relación es monótona por tramos de calado. Son dos
instrumentos que no tienen nada que ver midiendo lo mismo, y coinciden.
Lo que este informe no puede decir
- NDTI y NDCI son proxies sin calibrar: no hay unidades de turbidez ni miligramos de
clorofila detrás, y sus umbrales de literatura no son transferibles a lagunas someras y
salinas. Valen para comparar un humedal consigo mismo, no para dar cifras absolutas.
- El agua abierta y la vegetación inundada se intercambian: En Tablas de Daimiel la lámina detectada pasó de 80 a 200 y a 26 ha en diez días de agosto de 2026, con imágenes en color natural casi idénticas. Esas
hectáreas no desaparecen, cambian de clase: el agua bajo vegetación, o con algas en
superficie, queda en el filo del umbral y pasa a contarse como vegetación inundada, que es la
línea discontinua de la gráfica. Cuando una baja y la otra sube a la vez, no ha pasado nada
en el humedal. Por eso las alertas usan la mediana de las últimas observaciones y no la
última, y por eso las hectáreas de una fecha suelta no son una medida fiable de nada.
- Nueve años son pocos: la referencia empieza en julio de 2017, cuando empieza el
archivo de Sentinel-2 L2A. Un año excepcional puede parecer normal si solo hay ocho con los
que compararlo.
- Cuidado con los años incompletos: 2017 arranca en julio y el año en curso acaba hoy,
así que sus medias anuales no cubren el mismo periodo que las demás. Comparándolas sin más
apareció una falsa tendencia de eutrofización en el Mar Menor que era, simplemente, el otoño
que aún no ha ocurrido. Las alertas no caen en esa trampa porque siempre comparan con la
misma época del año.
- No distingue la causa: una desecación puede venir de que no ha llovido, de que el
acuífero ha bajado o de que alguien ha abierto una compuerta. El satélite no lo sabe; para eso
está el contraste con las medidas de campo, que hoy solo existe en Doñana.
- Una alerta es un aviso, no un diagnóstico: dice que algo se sale de lo que ese humedal
suele hacer en estas fechas, y que merece que alguien lo mire.
Glosario de siglas
- Sentinel-2
- Pareja de satélites de la Agencia Espacial Europea, más un tercero desde 2024, que fotografían Europa entera cada 5 días con píxeles de 10 a 20 m. Son gratuitos y públicos, y de ahí sale todo lo que hay en este informe.
- L2A Level-2A
- El nivel de proceso de la imagen: reflectancia de superficie, ya descontada la atmósfera. El nivel anterior, L1C, mide lo que ve el satélite; el L2A intenta medir lo que hay en el suelo. Esa corrección a veces falla, y eso da trabajo (ver más abajo).
- banda
- Cada uno de los colores que mide el satélite por separado, incluidos los que el ojo no ve. Aquí se usan seis: azul (B02), verde (B03), rojo (B04), borde del rojo (B05), infrarrojo cercano (B08) e infrarrojo de onda corta (B11).
- NIR / SWIR Near Infrared / Short-Wave Infrared
- Infrarrojo cercano y de onda corta. Importan porque el agua los absorbe casi por completo: en una imagen infrarroja una laguna es una mancha negra, y ahí está el truco para medirla.
- índice normalizado
- Todos los índices que siguen se calculan igual, (A − B)/(A + B) con dos bandas, y por eso valen siempre entre −1 y +1. Dividir por la suma cancela buena parte del efecto de la luz: el mismo sitio da un valor parecido en verano y en invierno, o con el sol más bajo.
- MNDWI Modified Normalized Difference Water Index
- (verde − SWIR)/(verde + SWIR). Es el detector de agua: sube donde hay lámina y baja en suelo y vegetación. Es el índice del que sale la superficie inundada.
- NDWI Normalized Difference Water Index
- (verde − NIR)/(verde + NIR). El hermano mayor del anterior, aquí solo como contraste: si los dos discrepan mucho, la fecha es sospechosa.
- NDVI Normalized Difference Vegetation Index
- (NIR − rojo)/(NIR + rojo). Mide cuánta vegetación viva hay. Sirve para separar el agua abierta del carrizo, la masiega o el arrozal, que también están inundados pero se ven verdes desde arriba.
- NDTI Normalized Difference Turbidity Index
- (rojo − verde)/(rojo + verde). Proxy de turbidez: el agua con sedimento en suspensión refleja más en el rojo. No da unidades; indica si esta fecha está más turbia de lo normal.
- NDCI Normalized Difference Chlorophyll Index
- (borde del rojo − rojo)/(borde del rojo + rojo). Proxy de clorofila, es decir de cuántas algas hay en el agua. Es el índice con el que se vigila la eutrofización.
- SCL Scene Classification Layer
- Un mapa que la propia ESA publica junto a la imagen y que etiqueta cada píxel: nube, sombra, agua, vegetación, suelo… Se usa para tirar las fechas con nubes y, además, como segunda opinión independiente sobre dónde hay agua.
- Otsu
- Método de 1979 (Nobuyuki Otsu) para elegir el corte que mejor parte un histograma en dos grupos. Aquí decide, en cada fecha y para cada humedal, a partir de qué valor de MNDWI un píxel cuenta como agua.
- STAC / COG SpatioTemporal Asset Catalog / Cloud-Optimized GeoTIFF
- El catálogo por el que se buscan las imágenes y el formato en el que están guardadas. El COG permite descargar solo el recorte del humedal en vez de la escena entera, que son cientos de megas.
- Natura 2000
- La red europea de espacios protegidos. Sus polígonos oficiales son el contorno con el que se recorta cada humedal; cada sitio tiene un código como ES0000024.
- Ramsar
- Convenio internacional de 1971 sobre humedales de importancia internacional. España tiene 76 sitios inscritos; los seis de este informe están entre ellos.
- ETRS89 / UTM 30N / EPSG:25830
- El sistema de coordenadas en el que se trabaja: metros, no grados. Es lo que permite hablar de hectáreas con sentido, porque cada píxel mide lo mismo en toda la imagen.
- ha hectárea
- 10.000 m², una parcela de 100 × 100 m. A 20 m de píxel cada píxel son 0,04 ha; a 40 m, 0,16 ha.
- percentil
- El percentil 10 de una serie es el valor que solo el 10 % de las observaciones baja. Casi todas las alertas de este informe se definen así: no con un umbral absoluto, sino preguntando si el valor de hoy está entre los peores de lo que ese humedal suele dar en estas fechas.
- ICTS Infraestructura Científica y Técnica Singular
- La categoría oficial de las grandes instalaciones científicas españolas. La ICTS-Doñana, de la Estación Biológica de Doñana (CSIC), es la que publica el calado medido dentro de la marisma que aparece en el panel inferior de las gráficas de Doñana.
- CSIC Consejo Superior de Investigaciones Científicas
- El mayor organismo público de investigación de España, del que depende la Estación Biológica de Doñana.
- CC BY 4.0 Creative Commons Attribution 4.0
- Licencia que permite reutilizar los datos, incluso comercialmente, con la única condición de citar la fuente.