Lo que no se puede medir

Intentamos medir el área quemada de los veinte incendios más grandes del Perú. Lo conseguimos con cuatro. Ese fracaso, y no las hectáreas, es el resultado.

Gabriela Pozo · Reto C

Corte de datos: 21 de octubre de 2026 · Trabajo académico independiente, no reemplaza los reportes del COEN, INDECI, SERFOR ni SERNANP.

Cómo leer esta página

Foco de calor

Un píxel de unos 375 m donde un satélite VIIRS detectó fuego. Incluye quemas agrícolas y de pastos, no solo incendios forestales.

Incendio

Aquí, un grupo de focos vecinos: dos focos van al mismo incendio si están a 1 km o menos y a 48 h o menos entre sí.

Envolvente

El polígono que rodea los focos de un incendio. Mide la extensión del grupo, no el área quemada.

dNBR

La diferencia del índice NBR entre una imagen Sentinel-2 de antes y otra de después. Mide cambio en la vegetación, no fuego.

«Medible»

Que existe un par de imágenes Sentinel-2, antes y después, con poca nube sobre la zona del incendio: como máximo el 30 % tapado.

Colores de severidad

  • Baja · dNBR 0,10–0,27
  • Moderada-baja · 0,27–0,44
  • Moderada-alta · 0,44–0,66
  • Alta · 0,66 o más

Por debajo de 0,10 se cuenta como no quemado.

1El resultado principal: el 80 % no era medible

Sobre cada uno de los 20 incendios con más focos de calor entre el 14 y el 21 de octubre de 2026 se corrió exactamente el mismo procedimiento: buscar una imagen Sentinel-2 antes y otra después del fuego, enmascarar nubes, calcular el dNBR y cruzarlo con un mapa de cobertura. La regla para aceptar una imagen se fijó antes de mirar los datos: como máximo el 30 % de la zona tapada.

4 de 20

incendios medidos: una tasa de éxito del 20 %.

13

descartados por nubes sobre la zona; otros 2 sin ninguna pasada de Sentinel-2 y 1 por el prefiltro.

p = 0,62

sin diferencia significativa entre Andes y Amazonía en la tasa de éxito.

El sesgo no es geográfico sino de calendario. En la Amazonía se midió el 12,5 % de los incendios y en los Andes el 25,0 %, una diferencia que la prueba exacta de Fisher no distingue del azar (p = 0,62). En octubre empieza la temporada de lluvias en ambas regiones. La consecuencia es incómoda: la temporada de fuego coincide con la temporada en que no se ve el suelo, así que cualquier indicador nacional de área quemada en tiempo casi real mide solo el subconjunto de incendios que las nubes dejaron pasar.

IncendioRegiónFocosEnvolvente (ha)Quemado (ha)Resultado
Chichas (Arequipa)Andes1782 9571 689,5Medido · pastizal
Iparía (Ucayali)Amazonía108724255,4Medido · bosque
San Cristóbal (Moquegua)Andes59587609,9Medido · pastizal
Curahuasi (Apurímac)Andes45475306,0Medido · pastizal
Río Tambo (Junín)Amazonía2793 051—Nubes sobre la zona (55,4 % de la zona inválida)
Curahuasi (Apurímac)Andes2013 227—Sin pasada de Sentinel-2
Callería (Ucayali)Amazonía1161 530—Nubes sobre la zona (47,2 % de la zona inválida)
Huanoquite (Cusco)Andes801 220—Nubes sobre la zona (100,0 % de la zona inválida)
Nueva Requena (Ucayali)Amazonía65416—Nubes sobre la zona (77,7 % de la zona inválida)
Cayna (Huánuco)Andes63724—Nubes sobre la zona (100,0 % de la zona inválida)
Mazamari (Junín)Amazonía62793—Sin pasada de Sentinel-2
Palca (Junín)Andes61853—Nubes sobre la zona (100,0 % de la zona inválida)
Santa Cruz (Loreto)Amazonía50560—Nubes sobre la zona (100,0 % de la zona inválida)
Nueva Requena (Ucayali)Amazonía49315—Descartada por el prefiltro (18,2 % de la zona inválida)
Nueva Requena (Ucayali)Amazonía48308—Nubes sobre la zona (100,0 % de la zona inválida)
San Juan de Jarpa (Junín)Andes42497—Nubes sobre la zona (89,4 % de la zona inválida)
Caja (Huancavelica)Andes39181—Nubes sobre la zona (47,5 % de la zona inválida)
Ayahuanco (Ayacucho)Andes39400—Nubes sobre la zona (55,1 % de la zona inválida)
Ancahuasi (Cusco)Andes36259—Nubes sobre la zona (65,7 % de la zona inválida)
San Miguel (Ayacucho)Andes35264—Nubes sobre la zona (63,7 % de la zona inválida)

Nueva Requena (Ucayali, 49 focos): la zona estaba limpia en un 81,8 %, pero el pipeline descarta de entrada las escenas con más del 80 % de nubes en la escena completa. Es un caso medible perdido por una decisión propia. No se corrigió para no cambiar el método a mitad del análisis.

Se excluyó un grupo de focos en Paracas (Ica, 42 focos en 49 ha) con un criterio escrito y aplicado a los 40 primeros candidatos: no es un incendio si más de la mitad de las detecciones S-NPP de 2012-2025 en su zona son de tipo 2 en FIRMS (fuente estática terrestre, como una planta industrial). En Paracas lo son 3 770 de 3 838 (98,2 %), detectadas en 2 139 días distintos; en los demás candidatos, ninguna. Su lugar lo ocupa San Miguel (Ayacucho), el siguiente por número de focos.

Los cuatro que sí se pudieron medir

Mismo método, misma escala de colores: severidad (dNBR) sobre el infrarrojo cercano de la imagen posterior, con los focos VIIRS encima y el contorno de la zona analizada.

Mapa de severidad del incendio de Chichas (Arequipa): 1 689,5 hectáreas con dNBR de 0,10 o más sobre el infrarrojo cercano posterior, con 178 focos VIIRS y el contorno de la zona analizada.
Chichas · Arequipa · Andes
1 689,5 ha · pastizal dominante
Antes 15 oct (0,0 % de nubes en la zona) → después 20 oct (0,2 %)
Mapa de severidad del incendio de Iparía (Ucayali): 255,4 hectáreas con dNBR de 0,10 o más sobre el infrarrojo cercano posterior, con 108 focos VIIRS y el contorno de la zona analizada.
Iparía · Ucayali · Amazonía
255,4 ha · bosque dominante
Antes 23 sep (0,0 % de nubes en la zona) → después 20 oct (19,0 %)
Mapa de severidad del incendio de San Cristóbal (Moquegua): 609,9 hectáreas con dNBR de 0,10 o más sobre el infrarrojo cercano posterior, con 59 focos VIIRS y el contorno de la zona analizada.
San Cristóbal · Moquegua · Andes
609,9 ha · pastizal dominante
Antes 12 oct (0,0 % de nubes en la zona) → después 20 oct (0,0 %)
Mapa de severidad del incendio de Curahuasi (Apurímac): 306,0 hectáreas con dNBR de 0,10 o más sobre el infrarrojo cercano posterior, con 45 focos VIIRS y el contorno de la zona analizada.
Curahuasi · Apurímac · Andes
306,0 ha · pastizal dominante
Antes 13 oct (0,0 % de nubes en la zona) → después 18 oct (0,9 %)

Leyenda común: amarillo, severidad baja; naranja, moderada-baja; rojo, moderada-alta; granate, alta; celeste, sin dato por nube o agua; puntos cian, focos VIIRS. Cada mapa trae su barra de escala. Las zonas tienen tamaños distintos.

2El mapa nacional y cómo se comporta el fuego

Mapa del Perú con los 8 637 focos VIIRS del 14 al 21 de octubre de 2026 coloreados por potencia radiativa, concentrados en la selva central y los Andes del sur, y los 20 incendios con más focos marcados con círculos: cuatro verdes (medidos) y dieciséis grises (no medibles).
Focos VIIRS por potencia radiativa (FRP) y los 20 incendios con más focos: verde, los 4 medidos; gris, los 16 que no se pudieron medir.
79,9 %

de los 8 637 focos fue diurno.

89,5 %

de los 2 602 incendios duró menos de 12 h entre su primer y su último foco.

6,55 MW

de potencia radiativa (FRP) mediana por foco.

11

de 2 602 incendios superaron los 3 días.

El fuego peruano de esa semana fue mayoritariamente corto, de baja potencia y de día. VIIRS detecta mejor de noche, cuando el suelo está frío y el contraste térmico es mayor, así que el predominio diurno no se explica por el sensor.

1 345 de los 2 602 incendios tienen un solo foco; para ellos la duración es cero. La duración mide el tiempo entre la primera y la última detección, no cuánto ardió el fuego.

3El caso medido: Chichas, 1 690 ha donde nunca había ardido

Mapa del área quemada en Chichas, Arequipa: una mancha extensa, mayormente amarilla (severidad baja) con zonas naranjas (moderada-baja), con los 178 focos VIIRS como puntos celestes y la zona analizada rodeada por una línea negra.
Severidad (dNBR) sobre el infrarrojo cercano de la imagen Sentinel-2 del 20 de octubre de 2026. Puntos celestes: los 178 focos VIIRS. Línea negra: la zona analizada (focos + 1 km). Coordenadas UTM 18S.
Imagen Sentinel-2 en color real del 15 oct de 2026: laderas pardas y quebradas de la puna de Chichas, sin cicatriz. Imagen Sentinel-2 en color real del 20 oct de 2026 del mismo terreno: una gran mancha negra de quemado ocupa el centro; abajo a la derecha hay nubes blancas y sus sombras. 15 oct20 oct
Chichas antes y después, Sentinel-2 en color real (10 m), mismo recorte y mismo ajuste de color en las dos fechas. A la izquierda de la línea, el 15 oct; a la derecha, el 20 oct. Desliza para mover la línea. El recorte es el rectángulo que contiene la zona analizada; las manchas blancas y oscuras de la esquina inferior derecha del 20 oct son nubes y sus sombras, no quemado.
1 690 ha

quemadas (dNBR ≥ 0,10), el 31 % de las 5 399 ha analizadas.

4,3 %

del distrito de Chichas (38 976 ha según los límites GADM 4.1).

0

focos dentro del área quemada entre 2012 y 2025.

Recurrencia: catorce años sin un solo foco

dentro del área quemadaanillo de 1 km 0 10 20 30 40 50 60 70 2012: 0 focos dentro del área, 0 en el anillo de 1 km 0 2012 2013: 0 focos dentro del área, 1 en el anillo de 1 km 1 2013 2014: 0 focos dentro del área, 0 en el anillo de 1 km 0 2014 2015: 0 focos dentro del área, 0 en el anillo de 1 km 0 2015 2016: 0 focos dentro del área, 0 en el anillo de 1 km 0 2016 2017: 0 focos dentro del área, 0 en el anillo de 1 km 0 2017 2018: 0 focos dentro del área, 4 en el anillo de 1 km 4 2018 2019: 0 focos dentro del área, 0 en el anillo de 1 km 0 2019 2020: 0 focos dentro del área, 0 en el anillo de 1 km 0 2020 2021: 0 focos dentro del área, 0 en el anillo de 1 km 0 2021 2022: 0 focos dentro del área, 0 en el anillo de 1 km 0 2022 2023: 0 focos dentro del área, 0 en el anillo de 1 km 0 2023 2024: 0 focos dentro del área, 0 en el anillo de 1 km 0 2024 2025: 0 focos dentro del área, 0 en el anillo de 1 km 0 2025 2026*: 44 focos dentro del área, 17 en el anillo de 1 km 2026*: 44 focos dentro del área, 17 en el anillo de 1 km 61 2026*
Focos Suomi-NPP de confianza nominal o alta por año: oscuro, dentro del área que ardió en 2026; claro, en el anillo de 1 km alrededor. *2026: solo los ocho días del 14 al 21 de octubre.

Entre 2012 y 2025, con entre 346 y 366 días de registro al año, el área que ardió en 2026 no registró ni un solo foco. En el anillo de 1 km hubo uno en 2013 y cuatro en 2018. En los ocho días de registro de octubre de 2026 hubo 44 focos dentro y 61 contando el anillo.

Queda descartada una quema de manejo recurrente. Una zona de quema habitual aparece en el registro año tras año, y esta no aparece nunca. La salvedad es de escala: una quema muy pequeña o muy rápida puede ocurrir entre dos pasadas del satélite sin dejar foco.

Además, el 56,7 % de los focos de Chichas fue nocturno: la tercera firma más nocturna de los veinte, en un país donde el 79,9 % de los focos fue diurno.

4Reto C: qué tipo de suelo ardió

Cobertura (WorldCover 2021)Hectáreas% de lo quemado
Pastizal1 223,072,4 %
Matorral455,627,0 %
Suelo desnudo10,30,6 %
Agua0,40,0 %
Bosque0,10,0 %
Cultivo0,00,0 %
Total1 689,4100 %

Cero hectáreas de cultivo y una décima de bosque. El daño es pérdida de forraje y el afectado es la ganadería extensiva. Con recurrencia nula, además, no es forraje que se queme en un ciclo de manejo.

5Tres hipótesis que cayeron

a) «En pastizal la severidad sale baja y en bosque alta»: no se sostiene

Baja Moderada-baja Moderada-alta Alta 0 % 25 % 50 % 75 % 100 % Chichaspastizal Chichas · Baja: 85,7 % Chichas · Moderada-baja: 14,3 % Chichas · Moderada-alta: 0,1 % Curahuasipastizal Curahuasi · Baja: 56,0 % Curahuasi · Moderada-baja: 37,1 % Curahuasi · Moderada-alta: 6,6 % Curahuasi · Alta: 0,3 % Iparíabosque Iparía · Baja: 42,2 % Iparía · Moderada-baja: 40,0 % Iparía · Moderada-alta: 17,6 % Iparía · Alta: 0,1 % San Cristóbalpastizal San Cristóbal · Baja: 41,1 % San Cristóbal · Moderada-baja: 57,8 % San Cristóbal · Moderada-alta: 1,1 %
Distribución de severidad (escala de Key y Benson) del área quemada en los cuatro incendios medidos.

Los pastizales van del 14,4 % al 58,9 % de severidad moderada o superior, y el único bosque medible, Iparía, cae dentro de ese rango con 57,7 %. Chichas es el atípico. Con n = 4 no se puede concluir que la escala de Key y Benson esté mal calibrada para el pastizal andino; para eso hace falta severidad medida en campo. Y no encontramos ningún estudio de dNBR con validación de campo en puna o páramo: ese vacío es en sí mismo un hallazgo.

b) «Los focos sobrestiman de forma sistemática»: no

00 500500 1 0001 000 1 5001 500 2 0002 000 2 5002 500 3 0003 000 medida = envolventeEnvolvente de los focos (ha)Área medida (ha) Chichas: envolvente 2 957 ha, medida 1 689,5 haChichas (−42,9 %) Iparía: envolvente 724 ha, medida 255,4 haIparía (−64,7 %) San Cristóbal: envolvente 587 ha, medida 609,9 haSan Cristóbal (+3,9 %) Curahuasi: envolvente 475 ha, medida 306,0 haCurahuasi (−35,6 %)
Envolvente de los focos frente al área medida. Verde: pastizal dominante; ocre: bosque. Bajo la diagonal, la envolvente sobrestima.

−64,7 % en Iparía, −42,9 % en Chichas, −35,6 % en Curahuasi y +3,9 % en San Cristóbal. Tres sobrestiman y uno subestima, en un rango demasiado amplio para un factor de corrección. Lo que queda en pie es que contar focos no mide hectáreas.

c) «Las nubes castigan sobre todo a la Amazonía»: no

RegiónIncendiosMedidos% medidosNubesSin pasadaPrefiltro
Amazonía8112,5 %511
Andes12325,0 %810

Fisher, p = 0,62. Las nubes tumbaron el 62,5 % de los casos amazónicos y el 66,7 % de los andinos.

6Robustez

0 500 1 000 1 500 2 000 Umbral de dNBR(0,05–0,35) · 1 973,6 ha Umbral de dNBR: mínimo 33,9 ha Umbral de dNBR: máximo 2 007,5 ha 34 2 008 Par de imágenes(2.ª mejor) · 412,8 ha Par de imágenes: mínimo 1 276,7 ha 1 277 Máscara SCL(permisiva – estricta) · 216,1 ha Máscara SCL: mínimo 1 473,4 ha 1 473 Zona(500 m – 2 km) · 4,2 ha Zona: mínimo 1 688,9 ha Zona: máximo 1 693,1 ha 1 689 1 693
Rango de hectáreas quemadas en Chichas al cambiar una decisión por vez. La línea vertical es la base (1 689,5 ha).

El umbral recorre 1 973,6 ha; el par de imágenes, 412,8; la máscara, 216,1; la zona, 4,2. El umbral solo mueve 3,1 veces más que las otras tres decisiones juntas (633,1 ha).

0 250 500 750 1 000 1 250 1 500 1 750 2 000 2 250 0,05 0,08 0,10 0,15 0,20 0,27 0,35 Umbral de dNBR Umbral 0,05: 2 007,5 ha Umbral 0,08: 1 842,9 ha Umbral 0,10: 1 689,5 ha Umbral 0,15: 1 168,9 ha Umbral 0,20: 686,6 ha Umbral 0,27: 241,9 ha Umbral 0,35: 33,9 ha Umbral usado (0,10): 1 690 ha
Hectáreas quemadas en Chichas según el umbral de dNBR. La línea discontinua marca el umbral usado (0,10).

La curva colapsa rápido porque los píxeles quemados se apiñan cerca del umbral. Con 0,35, a Chichas le queda el 2 % del área medida con 0,10; a Iparía, el 39 %.

Más rigor aparente, menos señal. Endurecer la máscara de nubes reduce el área en 216,1 ha, pero no porque quite nubes: la capa SCL de Sentinel-2 clasifica 178,7 ha de la propia cicatriz como «zona oscura» y 37,4 ha como «sin clasificar», y la máscara estricta las descarta.

7Validación contra MCD64A1

MCD64A1, el producto mensual de área quemada de MODIS, aún no existe para octubre de 2026, así que se validó el método en incendios de 2024 con una regla fijada antes de ver resultados.

00 1 0001 000 2 0002 000 3 0003 000 4 0004 000 5 0005 000 6 0006 000 MCD64A1 (ha)Nuestro método (ha) 2024-23859 (bosque): nuestro 2 286,7 ha, MCD64A1 2 212,2 ha 2024-28172 (bosque): nuestro 3 049,5 ha, MCD64A1 2 547,8 ha 2024-25613 (bosque): nuestro 1 724,9 ha, MCD64A1 2 543,6 ha 2024-21152 (bosque): nuestro 2 676,0 ha, MCD64A1 2 047,0 ha 2024-7631 (bosque): nuestro 2 159,8 ha, MCD64A1 1 255,0 ha 2024-7629 (bosque): nuestro 2 509,6 ha, MCD64A1 1 349,4 ha 2024-9748 (pastizal): nuestro 1 754,8 ha, MCD64A1 1 498,7 ha 2024-27870 (bosque): nuestro 2 777,2 ha, MCD64A1 2 622,8 ha 2024-28664 (bosque): nuestro 2 737,4 ha, MCD64A1 2 036,7 ha 2024-23334 (bosque): nuestro 3 138,0 ha, MCD64A1 3 155,7 ha 2024-23382 (bosque): nuestro 2 257,6 ha, MCD64A1 2 199,1 ha 2024-8669 (bosque): nuestro 568,4 ha, MCD64A1 281,7 ha 2024-27840 (bosque): nuestro 2 308,7 ha, MCD64A1 2 192,4 ha 2024-7552 (bosque): nuestro 1 463,6 ha, MCD64A1 1 183,4 ha 2024-23352 (bosque): nuestro 941,3 ha, MCD64A1 752,0 ha 2024-15326 (bosque): nuestro 1 459,1 ha, MCD64A1 1 082,8 ha 2024-23796 (bosque): nuestro 1 900,8 ha, MCD64A1 2 198,0 ha 2024-28027 (bosque): nuestro 1 455,5 ha, MCD64A1 1 643,4 ha 2024-23394 (bosque): nuestro 5 816,2 ha, MCD64A1 4 854,5 ha 2024-21811 (bosque): nuestro 2 165,2 ha, MCD64A1 2 066,0 ha 2024-27844 (bosque): nuestro 5 679,1 ha, MCD64A1 5 590,0 ha Los 3 de la regla Pastizal dominante Otros (bosque)
Hectáreas de nuestro método frente a MCD64A1 en los 21 candidatos medibles. Sobre la diagonal, nuestro método da más.
MuestranRazón mediana (nuestro / MCD64A1)RangoChichas corregida (ha)
3 casos (regla fijada)31,0481,016 – 1,1981 612
21 candidatos medibles211,1710,678 – 2,0181 443
solo pastizal dominante11,1711,171 – 1,1711 443

La cifra principal sigue siendo 1 690 ha; las corregidas son alternativas, por tres razones: la validación es casi toda bosque amazónico y el caso es pastizal andino; MCD64A1 usa píxeles de 500 m y es otra estimación satelital, no un dato de campo; y varios de los 21 candidatos se solapan en el espacio, así que no son independientes.

8Concentración y MAUP

0 %0 % 25 %25 % 50 %50 % 75 %75 % 100 %100 % igualdad perfecta% de unidades, de menos a más focos% acumulado de focos distrito: Gini 0,906 distrito · Gini 0,906 provincia: Gini 0,758 provincia · Gini 0,758 grilla 10 x 10 km: Gini 0,954 grilla 10×10 km · Gini 0,954 grilla 10 x 10 km desplazada 5 km: Gini 0,953 grilla 10×10 km desplazada 5 km · Gini 0,953
Curvas de Lorenz de los 8 637 focos con cuatro particiones, contando también las unidades sin focos.
ParticiónUnidades con focosGini (solo con focos)Gini (todas)
distrito644 de 1 8150,7350,906
provincia180 de 1950,7370,758
grilla 10 × 10 km1 588 de 13 4640,6090,954
grilla 10 × 10 km desplazada 5 km1 595 de 13 4530,6070,953

El Gini va de 0,607 a 0,954 según la unidad y según si se cuentan las unidades sin fuego. Cualquier índice de concentración que no declare su unidad no significa nada. Lo robusto es otra cosa: la distribución es muy desigual en las cuatro particiones.

9Implicaciones

10Qué haría falta para ponerle precio

La teledetección acota qué y cuánto se quemó, pero no determina si fue una pérdida.

11Método

  1. Focos: VIIRS de NASA FIRMS (Suomi-NPP, NOAA-20 y NOAA-21, CSV de 7 días descargado el 21 de octubre de 2026), confianza nominal o alta, filtrados a Perú con GADM 4.1.
  2. Incendios: dos focos pertenecen al mismo incendio si están a ≤ 1 km y ≤ 48 h entre sí.
  3. Par de imágenes: Sentinel-2 L2A de Microsoft Planetary Computer; se elige la imagen con menos zona inválida, antes y después del fuego.
  4. dNBR y severidad: NBR = (B8A − B12)/(B8A + B12); dNBR = NBR antes − NBR después; clases de Key y Benson.
  5. Cobertura: ESA WorldCover 2021 cruzado con los píxeles quemados.
  6. Contraste: recurrencia con el archivo anual S-NPP 2012-2025, sensibilidad a cada decisión y validación contra MCD64A1.
ImagenFechaTileNubes en la zonaNubes en la escena
Antes2026-10-1518LYH0,0 %0,3 %
Después2026-10-2018LYH0,2 %11,4 %

Código, tablas y cómo reproducir: src/, outputs/ y el README del repositorio.

12Limitaciones y qué no se pudo hacer