Metodología

De dónde parte

Que la cercanía entre plantaciones y viviendas aumenta el riesgo de incendio está documentado en Chile. McWethy et al. (2018) mostraron que entre 2001 y 2017 el fuego ocurrió de preferencia en plantaciones exóticas. Jaque Castillo et al. (2021) identificaron la cercanía entre la expansión urbana y las plantaciones como el factor principal del riesgo en el Gran Concepción. Hora et al. (2025) encontraron más plantación a menos de 30 m de las viviendas destruidas en el incendio de 2023 en el Biobío. Y al presentar MapBiomas Fuego, sus autores informaron que cerca del 30 % de los incendios afectó áreas próximas a ciudades (Universidad de Chile, 2026).

También existen mapas nacionales de interfaz urbano-forestal. CONAF publica un análisis nacional de riesgo y zonas de interfaz (2021, actualizado en 2024), con resolución de 10 m y un visor público de riesgo en las zonas de interfaz. Según su descripción, no separa la plantación forestal de la vegetación nativa ni compara años. Hay además mapas de interfaz para todo el país en Schug et al. (2023) y Gelabert et al. (2026). Sarricolea et al. (2020) estimaron que cerca de 3 millones de personas viven en la interfaz urbano-forestal entre Valparaíso y La Araucanía, con resultados por región, y Hernández et al. (2024) midieron con MapBiomas el crecimiento de la infraestructura urbana del país entre 2000 y 2022. El Ministerio de Educación (visor de comunidades educativas) y SENAPRED (Visor Chile Preparado) ya muestran colegios, jardines y centros de salud sobre la amenaza de incendio de CONAF, así que Límite Vivo no repite esa capa.

En nuestra revisión no encontramos una herramienta pública que, para cada comuna del país y con las series de MapBiomas, cuantifique la plantación forestal a 100, 200 y 500 m de zonas pobladas, la siga desde el año 2000, separe si fue la ciudad la que llegó a la plantación o la plantación la que llegó a la ciudad, cuente a las personas que viven al lado con el Censo 2024 y lo contraste con el área quemada dentro de la misma comuna y la misma franja. Eso es lo que agrega Límite Vivo, como complemento de los trabajos anteriores.

Hay además una necesidad práctica. La Ley 21.364 (artículo 28) exige a cada comuna un Plan Comunal para la Reducción del Riesgo de Desastres, y según Corporación Ciudades, con datos de SENAPRED al 6 de julio de 2026, 127 de 345 comunas no lo tenían (nota de prensa). El proyecto de ley de incendios forestales y rurales (boletín 16.335-14), que propone zonas de interfaz y de amortiguación, seguía en trámite en octubre de 2026.

Datos

• Cobertura y uso del suelo de MapBiomas Chile, Colección 1, años 2000 y 2022, y 2000 a 2005 para saber qué plantación ya existía al inicio. Es la base de las cifras principales. Se leyó del asset público de la Colección 1 en Google Earth Engine, en su grilla nativa de 30 m; da exactamente las mismas cifras que los GeoTIFF anuales del almacenamiento público de MapBiomas, que se usaron primero y que el 30 de octubre de 2026 dejaron de estar disponibles para descarga directa.
• Cobertura y uso del suelo de MapBiomas Chile, Colección 2, años 2000, 2022 y 2024, leída en Google Earth Engine y descargada en la misma grilla que la Colección 1. Se usa como análisis de sensibilidad (ver control de calidad).
• Área quemada anual de MapBiomas Fuego Chile, Colección 1, años 2023, 2024 y 2025, y frecuencia de quema 2013-2025. Cubre de la Región de Valparaíso a la de La Araucanía, así que todo lo relacionado con incendios solo se calcula en esas regiones.
• Límites comunales de MapBiomas Chile (nivel político 4, 343 comunas), los mismos que usa MapBiomas en sus estadísticas. Están publicados como ráster en Earth Engine; se vectorizaron y se nombraron cruzándolos con la capa oficial de comunas del MOP (servicio DPA), con la que cada comuna coincide en más del 92 % de su superficie.
• Cartografía del Censo de Población y Vivienda 2024 (INE): manzanas urbanas y de aldeas, con personas, viviendas ocupadas, grupos de edad y personas con dificultad para caminar o subir escaleras.
• Estadísticas oficiales por comuna de las Colecciones 1 y 2, para el control de calidad y el contexto de cada ficha.

Método

1. Zonas pobladas. Se parte de la clase infraestructura urbana (código 24), que según el equipo de MapBiomas Chile incluye también las carreteras principales (Hernández et al., 2024). Para quedarse con ciudades, pueblos y villorrios se eliminan las formas de menos de 3 píxeles de ancho (unos 75 a 90 m) con una apertura morfológica, y luego las manchas de menos de 2 ha.
2. Distancia. Para cada píxel se calcula la distancia en metros hasta el borde de la zona poblada más cercana, con el tamaño real del píxel en esa latitud (unos 30 m de alto y entre 17 y 29 m de ancho). Se consideran zonas pobladas hasta cerca de 1 km fuera del límite comunal, porque una ciudad vecina también deja expuesta la plantación de la comuna.
3. Franjas. Se suman las hectáreas de plantación forestal (código 9) dentro de la comuna a menos de 100, 200 y 500 m de una zona poblada. El área de cada píxel depende de su latitud.
4. Comparación en el tiempo. El cálculo se repite para cada año disponible en cada colección.
5. Sectores. Con la cobertura de 2022 (Colección 1) se marca el borde de cada zona poblada que tiene plantación a menos de 100 m. Los tramos de borde expuesto que quedan a menos de unos 90 m entre sí se unen, y luego se cortan en celdas de 0,01° (cerca de 1 km) para que cada sector sea un trozo de borde que se pueda recorrer. Cada píxel de plantación a menos de 100 m se asigna al punto de zona poblada más cercano, y así se suman las hectáreas de cada sector. El largo del borde es aproximado, porque se cuenta en píxeles de 30 m. Un sector puede incluir plantación ubicada en una comuna vecina, porque lo que importa es el borde donde vive la gente.
6. Nombres de los sectores. Cada sector toma el nombre del lugar más cercano registrado en OpenStreetMap (barrio, población, villa o localidad; © colaboradores de OpenStreetMap, licencia ODbL). Si ese lugar está a más de 500 m se indica la distancia, y si está a más de 2 km el sector se identifica solo por sus coordenadas. Es una referencia para ubicarse, no un límite oficial.
7. Búsqueda por dirección. Las zonas pobladas y un margen de 200 m se dividen en celdas de 0,001° (unos 100 m). Para cada celda se guarda la distancia hasta el borde de la plantación más cercana (con tope de 3 km), las hectáreas de plantación en un radio de 500 m y si hubo área quemada a menos de 500 m entre 2023 y 2025. La dirección se ubica con sugerencias del buscador Photon, que usa datos de OpenStreetMap (lo escrito se envía a ese servicio), o con la ubicación del navegador, que no se envía a ningún servicio externo. Como OpenStreetMap no tiene el número de muchas casas en Chile, la persona puede ajustar el punto en un mapa opcional, que carga sus imágenes desde OpenStreetMap solo si se abre.
8. Quién se acercó a quién. La plantación que en 2022 está a menos de 100 o 500 m de una zona poblada se separa, píxel a píxel, en tres partes que suman exacto: la que ya estaba junto a las casas en 2000, la que ya existía pero más lejos y quedó cerca porque la zona poblada creció hasta ella, y la que se plantó después. Una plantación se considera existente si aparece como plantación en algún año entre 2000 y 2005; sin esa ventana, un rodal cosechado justo en 2000 se contaría como plantación nueva al replantarse. Se calcula con las dos colecciones.
9. Personas. Cada manzana del Censo 2024 se lleva a la grilla de MapBiomas (cada píxel queda con la manzana que contiene su centro; las manzanas más chicas que un píxel se miden en un punto interior) y se mide su distancia mínima al borde de la plantación de 2022. Si alguna parte de la manzana está a menos de 100 o 500 m, cuentan todas sus personas, viviendas y grupos de edad. Como en manzanas grandes eso incluye a quienes viven algo más lejos, se calcula también un criterio estricto, que cuenta la manzana solo si su punto interior está dentro de la distancia, y se informan las dos cifras como rango.
10. Historial de quema de los sectores. Para la franja de 100 m junto a cada sector, cualquiera sea su cobertura, se suma la superficie que se quemó al menos una vez entre 2013 y 2025 y la que se quemó dos o más veces, con la frecuencia de quema de MapBiomas Fuego.

Control de calidad

• La plantación y la infraestructura totales de cada comuna, calculadas desde los mapas de la Colección 1, se compararon con la estadística oficial de MapBiomas. En 960 comparaciones la diferencia mediana fue de 0,04 % y la máxima de 0,31 %.
• Con la Colección 2 la plantación también calza (diferencia mediana de 0,07 %), pero la infraestructura difiere cerca de 3 % de forma sistemática. La imagen disponible en Earth Engine no es exactamente la versión con que se publicaron las estadísticas, y por eso la Colección 2 se usa como análisis de sensibilidad y no como base.
• Las franjas se recalcularon en seis comunas con un método independiente (polígonos, buffer en metros en proyección UTM y área elipsoidal). En la franja de 500 m los dos métodos difieren menos de 1 %, en la de 200 m menos de 3 % y en la de 100 m hasta 9 %, porque a esa distancia el tamaño del píxel pesa más. En las seis comunas el método principal dio la cifra más baja, así que las hectáreas de la franja de 100 m tienden a quedarse cortas, no a exagerar.
• Hay pruebas automáticas para el filtro de zonas pobladas (un camino de 1 o 2 píxeles se elimina, una ciudad con un camino conserva la ciudad) y para el cálculo de Mantel-Haenszel. Las comunas sin plantación, como Macul, dan cero.
• La suma de la plantación de los 3.118 sectores (9.680 ha) difiere en 0,68 % de la cifra nacional a menos de 100 m (9.746 ha), porque se descartan tramos muy cortos y porque parte de la plantación se asigna a bordes de comunas vecinas. Ninguna comuna con 5 ha o más a menos de 100 m queda sin sectores, y la ubicación de los sectores se revisó sobre mapas de control en Tomé, Viña del Mar y Cabrero.
• La distancia guardada en 12 celdas al azar de la búsqueda por dirección se recalculó directamente sobre el ráster, con una diferencia máxima de 2 m.
• En cada comuna, las tres partes de "quién se acercó a quién" suman exactamente la franja del cálculo principal. Hay pruebas automáticas que comprueban la suma, que cada movimiento se reconoce y que un rodal replantado no se cuenta como plantación nueva.
• Todas las manzanas del censo quedaron medidas. Solo faltan Juan Fernández e Isla de Pascua, que no están entre las 343 comunas de MapBiomas. En los sectores, la superficie quemada entre 2013 y 2025 nunca es menor que la plantación quemada entre 2023 y 2025 de la misma franja, como debe ser.
• Revisando mapas de control de Viña del Mar, Santa Juana y Tomé se detectó que la infraestructura incluía carreteras, lo que llevó al filtro descrito arriba. La verificación independiente detectó que la primera versión medía la distancia de centro a centro de píxel, lo que subestimaba la franja de 100 m en cerca de 15 %. Ambos errores se corrigieron antes de publicar.

Contraste con los incendios de 2023 a 2025

Se usó la cobertura de 2022, anterior a los incendios, y el área quemada de 2023, 2024 y 2025. En la franja de 500 m de cada comuna se comparó la proporción quemada de la plantación con la del bosque nativo, el matorral, el pastizal y los humedales de la misma franja, y también con la de la agricultura.

La comparación se hace dentro de cada comuna para no mezclar lugares distintos, porque cada comuna tiene su clima, su relieve y sus propios incendios, y lo que interesa es qué se quemó más estando en la misma franja y bajo el mismo fuego. Los resultados por comuna se combinaron con la razón de Mantel-Haenszel, y el intervalo de confianza de 95 % se obtuvo con 2.000 remuestreos de comunas. Entran las comunas con al menos 20 ha de plantación en la franja y al menos 10 ha quemadas en ella.

Con las dos colecciones la plantación cercana a zonas pobladas se quemó en mayor proporción que la vegetación natural de la misma franja. Con la Colección 2 el intervalo queda claramente sobre 1. Con la Colección 1 el límite inferior llega a 0,98, así que con esa colección el resultado no alcanza a ser concluyente por sí solo.

El resultado no depende de los umbrales de inclusión. Con la Colección 1, al cambiar el mínimo de plantación y de área quemada la razón queda entre 1,74 y 1,82 veces, y al sacar una comuna a la vez queda entre 1,54 y 2,25 veces. Hay comunas donde pasó lo contrario, como Viña del Mar en 2024, y se muestran tal cual en sus fichas.

Limitaciones

• Mide exposición, no probabilidad de incendio. No considera dónde se inician los fuegos, el viento, la pendiente, la sequía ni el manejo de cada plantación.
• Depende de la clasificación de MapBiomas. Las dos colecciones difieren, sobre todo en la infraestructura, y eso cambia cuánto creció la franja entre 2000 y 2022 o 2024. Por eso se informan las dos.
• Con píxeles de 30 m, la franja de 100 m tiene una incertidumbre de método de hasta 10 %. No se consideran casas aisladas ni zonas pobladas de menos de 2 ha.
• El contraste con incendios cubre tres temporadas y solo las regiones de Valparaíso a La Araucanía. Con más años podría cambiar.
• Es una asociación, no una relación causal. La plantación podría quemarse más también por dónde está dentro de la franja, por ejemplo en laderas.
• "Quién se acercó a quién" depende de que la clasificación de zonas pobladas y plantación sea comparable entre años. Las dos colecciones dan proporciones distintas, por eso se informa el rango y no una sola cifra.
• Las personas son las de manzanas urbanas y de aldeas (89 % de la población del país): quedan fuera las viviendas rurales dispersas, que el INE publica en polígonos de referencia estadística. El censo es de 2024 y la plantación de 2022, y una manzana cuenta completa aunque solo una parte esté dentro de la distancia.
• La búsqueda de trabajos previos no fue exhaustiva. Si conoces un trabajo que haga lo mismo, lo citaremos.

Código y datos

Todo el procesamiento está en Python (rasterio, numpy, scipy y geopandas) y el código es público. Los resultados por comuna están disponibles en un archivo CSV.

Descargar datos de todas las comunas(CSV, 58 kB)

Fuentes de datos

• MapBiomas – Colección 1 de la Serie Anual de Mapas de Cobertura y Uso del Suelo de Chile, consultada el 30/10/2026 a través del enlace: https://chile.mapbiomas.org
• MapBiomas – Colección 2 de la Serie Anual de Mapas de Cobertura y Uso del Suelo de Chile, consultada el 30/10/2026 a través del enlace: https://chile.mapbiomas.org
• MapBiomas – Colección 1 de MapBiomas Fuego Chile (área quemada anual y frecuencia de quema), consultada el 30/10/2026 a través del enlace: https://chile.mapbiomas.org
• Instituto Nacional de Estadísticas (INE) – Cartografía del Censo de Población y Vivienda 2024, manzanas, consultada el 30/10/2026 en https://censo2024.ine.gob.cl/resultados/