Afección a Red Natura 2000 en Galicia por Incendios Forestales (2001-2023)

Trabajo sobre SIG
Título	Afección a Red Natura 2000 en Galicia por Incendios Forestales (2001-2023)
Autores	Rebeca Iglesias Buendía Damián Lorente García Victor M. Landázuri Salazar
Asignatura	Sistemas de Información Geográfica Aplicados a la Ingeniería Civil
Curso	Curso 24/25
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

Este trabajo tiene como objetivo el análisis de las afecciones por incendios forestales a los espacios protegidos de la Red Natura 2000 en Galicia, concretamente a las zonas ZEPA, ZEC/LIC, y ZEPA/ZEC/LIC, de tal manera que dicho análisis se ha realizado por medio de un software de Sistemas de Información Geográfica (SIG), concretamente QGIS, herramienta que ha permitido integrar, manipular y presentar datos espaciales y temporales.

Para ello, por un lado se han obtenido datos espaciales de incendios forestales a nivel europeo entre los años 2001 y 2023, y por otro lado datos espaciales de las zonas ZEPA, ZEC/LIC, y ZEPA/ZEC/LIC de la Red Natura 2000 en Galicia, y a partir de dichos datos espaciales se han obtenido valores numéricos que cuantifican la afección a dichas zonas protegidas de la Red Natura 2000 en Galicia, los cuales han sido representados convenientemente en múltiples cartografías distintas.

Como punto final, se ha incluido una conclusión en la cual se juzga, en función de los datos obtenidos, la justificación de la necesidad de actuación en las áreas estudiadas, debido a sus grados de afección, y en función de indicadores empleados para ello.

1 Introducción

Galicia es una de las comunidades autónomas de España con mayor riqueza natural. Su biodiversidad y variedad de ecosistemas han motivado la declaración de espacios protegidos, integrados en la Red Natura 2000 bajo las figuras de LIC (Lugar de Importancia Comunitaria), ZEC (Zona Especial de Conservación) y ZEPA (Zona de Especial Protección para las Aves). Actualmente, está compuesta por 59 ZECs y 16 ZEPAs, que supone un 38,9% de la superficie total de Galicia.

Como podremos observar más adelante, en el proyecto figuran 73 zonas protegidas Esto se debe a que hay dos zonas (Islas Cíes (ES0000001) y Serra de Enciña da Lastra (ES1130009)) que son a la vez ZEC/LIC y ZEPA.

Esta riqueza ambiental se ve constantemente amenazada por incendios forestales cuya frecuencia e intensidad son un grave riesgo ambiental, económico y social, impactando directamente en la vegetación, suelo, fauna…

A pesar de los esfuerzos institucionales para prevenir y controlar dichos incendios, el problema persiste y se ve agravado en algunos casos debido al abandono rural, acumulación de biomasas…
En este contexto, surge la necesidad de analizar en detalle el grado de afección que los incendios forestales han tenido sobre las zonas protegidas de Galicia en los últimos años.

El proyecto consiste en abordar esta problemática utilizando Sistemas de Información Geográfica (SIG), integrando las zonas protegidas y los incendios desde el año 2001 hasta 2023. El objetivo es cuantificar la superficie afectada de RED NATURA 200 debido a incendios forestales en Galicia, obteniendo una base objetiva de información para una mejor planificación, prevención, conservación, gestión y restauración del territorio.

2 Metodología

Lo primero se seleccionó dos capas ráster: una correspondiente a los incendios forestales (por años) y otra a las áreas protegidas (ZEPA, ZEC/LIC, ZEC LIC/ZEPA).

Recortamos ambas capas a la región de Galicia para optimizar el procesamiento: Barra de herramientas > Ráster > Extracción > Cortar ráster por extensión. (Ver Anejo I)

Las capas recortadas se poligonizaron usando la siguiente herramienta: Barra de herramientas > Ráster > Conversión > Poligonizar.

A continuación, recortamos los polígonos que se extendían fuera de Galicia (Portugal, Asturias y Castilla y León). Seleccionamos los polígonos y seleccionamos: Barra de herramientas > Editar > Editar geometría > Dividir objetos espaciales y eliminamos los elementos sobrantes tras el recorte.

A continuación, calculamos la superficie de áreas protegidas afectadas por incendios mediante la intersección de ambas capas. Durante este proceso, surgió un problema: algunas geometrías no estaban correctamente cerradas, impidiendo realizar la intersección. Para solucionarlo, validamos las geometrías mediante (Ver Anejo I) : Barra de herramientas > Vectorial > Herramientas de geometría > Comprobar validez. El análisis indicó que había dos geometrías inválidas. Para corregirlo, buscamos la herramienta "Corregir geometrías” Y volvimos a comprobar la validez. (Ver Anejo I)

Una vez corregidas las geometrías, volvimos a ejecutar la operación de intersección utilizando: Barra de herramientas > Vectorial > Herramientas de geoproceso > Intersección. Esto nos generó una nueva capa que representaba las zonas protegidas afectadas por incendios.

Para calcular el área de las zonas protegidas no afectadas, accedimos a la tabla de atributos de la nueva capa y utilizamos la calculadora de campo para obtener el área de cada polígono. (Ver Anejo I)

Los datos de la capa de zonas protegidas afectadas fueron extraídos y trabajados en Excel, se obtuvieron las áreas de las zonas afectadas. (Ver Anejo I)

Una vez obtenido el porcentaje de superficie afectada, cargamos en el proyecto el archivo Excel que habíamos elaborado previamente, en formato .XLSX, con el fin de vincularlo a la capa de zonas protegidas no afectadas.

Para ello, en la caja de capas, seleccionamos la capa correspondiente y accedimos a Propiedades. Dentro de la ventana de propiedades, nos dirigimos al apartado Uniones y allí configuramos la unión con el archivo Excel, estableciendo los campos comunes entre ambas tablas para vincular correctamente los datos. (Ver Anejo I)

Una vez unido el Excel representamos la capa de zonas afectadas en función de la gravedad de afección.

3 Resultados

En base al procedimiento de trabajo, se han obtenido los siguientes datos descriptivo-espaciales:

1. Cartografía temática de espacios de la Red Natura 2000 clasificados como zonas ZEPA, ZEC/LIC, y ZEPA/ZEC/LIC. (Ver Anejo III)
   
   
2. Cartografía temática de la distribución de incendios forestales en Galicia entre los años 2001 y 2023. (Ver Anejo III).
   
   
3. Cartografía específica del porcentaje de áreas afectadas por incendios forestales entre los años 2001 y 2023 de las zonas ZEPA, ZEC/LIC, y ZEPA/ZEC/LIC de la Red Natura 2000 en Galicia. (Ver Anejo III)
   
   
4. Tabla de identificación, caracterización y estado de las zonas ZEPA, ZEC/LIC, y ZEPA/ZEC/LIC de la Red Natura 2000 en Galicia (información complementaria al resto de cartografía ver Anejo III). Los hábitats y especies protegidos incluidos en esta tabla han sido identificados a través del portal oficial de la Unión Europea sobre la Red Natura 2000 biodiversity.europa.eu).
   
   
5. Cartografía específica de la distribución de áreas afectadas por incendios forestales entre los años 2001 y 2023 en las zonas ZEPA, ZEC/LIC, y ZEPA/ZEC/LIC de la Red Natura 2000 en Galicia.
   
   
6. Cartografía especifica de gravedad de afeccion por incendios en las zonas ZEPA, ZEC/LIC, y ZEPA/ZEC/LIC de la Red Natura 2000 en Galicia.

Distribucion de incendios

Gravedad de afeccion de incendios

4 Conclusiones

El análisis realizado ha permitido cuantificar y localizar con precisión el impacto de los incendios forestales sobre las zonas protegidas de Galicia.

Los resultados muestran que una parte significativa del territorio protegido ha sido afectado por incendios en varias ocasiones durante estos años, especialmente en provincias como Ourense y Lugo. Se han detectado ciertos años críticos con grandes superficies quemadas (como el año 2022), en general coincidiendo con condiciones meteorológicas extremas. Asimismo, se observa que algunas zonas presentan una recurrencia elevada, lo que podría indicar una mayor vulnerabilidad o falta de medidas preventivas eficaces.

Observando los resultados obtenidos, se plantean como recomendables las siguientes formas de actuación y/o compensación a los daños:

• Designar recursos públicos a la elaboración de datos de severidad de los incendios para valorar el impacto ecológico real, más allá de la mera afección superficial.
  
  
• Colaborar con administraciones locales y gestores ambientales para traducir los análisis realizados en acciones concretas de conservación y prevención.
  
  
• Incorporar imágenes satelitales (p. ej. Sentinel-2 o Landsat) para realizar análisis multitemporales más detallados y detectar variaciones en la cubierta vegetal post-incendio.
  
  
• Identificar aquellas actividades antrópicas que tengan potencial de generar incendios forestales, especialmente en aquellas zonas reincidentes.
  
  
• Desarrollar modelos predictivos de riesgo utilizando variables ambientales, topográficas y climáticas.
  
  

Por último, queremos destacar la utilidad de los SIG como una herramienta fundamental para el análisis espacial y temporal de fenómenos ambientales complejos. Su uso permite obtener una visión clara y detallada del problema, además de generar información valiosa que contribuye de manera

5 Anejos

5.1 ANEJO I. EXPLICACIÓN DE LA METODOLOGIA APLICADA

Introducción Dado que el apartado de metodología del cuerpo principal del trabajo presentaba una extensión limitada, en este anejo se expone de forma más detallada el procedimiento seguido mediante el software QGIS para la elaboración del análisis geográfico. Se describen las fuentes de datos utilizadas, los procesos de tratamiento, análisis y representación espacial aplicados, así como las herramientas específicas que permitieron obtener los resultados presentados en el informe principal.

El objetivo de este anejo es ofrecer una visión técnica y transparente del flujo de trabajo en QGIS, facilitando la comprensión y posible replicación del estudio por parte de otros usuarios o evaluadores.

1-Obtencion de datos
1.1 Obtención de la capa de masa forestal quemada
En primer lugar, se procedió a identificar y recopilar las capas necesarias para el desarrollo del análisis. La capa principal correspondía a las zonas de masa forestal quemada, elemento central del estudio.

Inicialmente, se consultó el Instituto Galego de Estatística (IGE), donde se localizaron estadísticas relacionadas con incendios forestales. Sin embargo, los datos disponibles abarcaban únicamente el período 2001–2006 y no incluían cartografía geoespacial detallada, por lo que esta fuente fue descartada.

También se revisaron los recursos del Ministerio para la Transición Ecológica y el Reto Demográfico (MITECO). Aunque se encontraron indicadores como la frecuencia de incendios forestales por municipios, no se disponía de capas espaciales con las áreas quemadas, motivo por el cual esta fuente también fue descartada.

Finalmente, se identificó como fuente adecuada la plataforma Global Forest Watch (GFW), desarrollada por el World Resources Institute (WRI). Esta herramienta proporciona acceso a datos geoespaciales actualizados sobre cobertura forestal, pérdida de masa arbórea, incendios y conservación de bosques. A través del sitio web de GFW se redirige a fuentes técnicas específicas, como el portal de la Universidad de Maryland, donde se encuentra disponible la serie de datos satelitales correspondiente.

La capa utilizada fue descargada desde el siguiente enlace oficial: https://glad.umd.edu/users/Alexandra/Fire_GFL_data/2001-23/, que contiene datos anuales de áreas afectadas por fuego a nivel global desde 2001 hasta 2023, derivados del análisis de imágenes Landsat por el GLAD Lab (Global Land Analysis & Discovery) de dicha universidad.

Se descargó la capa correspondiente a Europa anual, con resolución espacial de 30 metros en formato raster (.tif).

Cabe señalar que, en el momento de la descarga (abril de 2025), los datos correspondientes al año 2024 aún no estaban disponibles en la fuente oficial, por lo que el análisis se limitó a la información disponible hasta el año 2023.

Global Forest Watch (GFW) y el portal del GLAD Lab constituyen fuentes de datos altamente fiables y reconocidas internacionalmente. La información utilizada proviene del análisis de imágenes satelitales Landsat, procesadas por el Global Land Analysis and Discovery (GLAD) Lab, un grupo de investigación de la Universidad de Maryland especializado en monitoreo forestal. Además, GFW es una iniciativa respaldada por el World Resources Institute (WRI) y cuenta con el apoyo de socios estratégicos como Google, NASA, Esri, y organizaciones como UNEP-WCMC o Rainforest Foundation Norway, lo que garantiza la transparencia, solidez técnica y actualización periódica de los datos que ofrece.

1.2 Capa de Zonas Protegidas
El segundo objetivo fue la descarga de la capa de espacios naturales protegidos, necesaria para analizar la relación entre los incendios forestales y las áreas sujetas a algún régimen de protección ambiental.

Esta información se obtuvo desde el portal oficial del Ministerio para la Transición Ecológica y el Reto Demográfico (MITECO), específicamente a través del Centro de Descargas del Centro Nacional de Información Geográfica (CNIG). La capa fue descargada en formato vectorial Shapefile (.shp), ampliamente compatible con QGIS.

La capa incluye la delimitación de diversas figuras de protección, tales como Parques Nacionales, Parques Naturales, Reservas de la Biosfera, Lugares de Importancia Comunitaria (LIC), Zonas de Especial Conservación (ZEC) y Zonas de Especial Protección para las Aves (ZEPA), entre otras. Esta información fue clave para evaluar si las zonas afectadas por incendios coincidían con áreas de especial valor ecológico.

1.3 Capas de la comunidad Autonoma de Galicia y sus provincias
Con el fin de proporcionar un mayor contexto geográfico a los mapas elaborados, se incorporaron capas complementarias obtenidas a través del Centro de Descargas de Información Geográfica de la Xunta de Galicia .En concreto, se utilizaron las capas vectoriales correspondientes a los límites autonómicos y provinciales, lo que permitió enmarcar espacialmente los resultados del análisis de manera más precisa y comprensible.

2 Sistema de referencia de coordenadas
Todas las capas espaciales utilizadas en el proyecto se importaron inicialmente en el sistema de referencia WGS 84 (EPSG:4326), que es un sistema geodésico global basado en coordenadas geográficas (latitud y longitud). Sin embargo, para garantizar la precisión en los cálculos de superficie y otras operaciones espaciales, se procedió a reproyectar las capas al sistema de referencia UTM zona 29N (EPSG:25829), que utiliza coordenadas métricas (metros) y es más adecuado para el área geográfica de Galicia.

Esta reproyección se realizó antes de ejecutar el cálculo del área ($area) en la tabla de atributos, asegurando que las medidas resultantes fueran exactas y consistentes con la realidad territorial.

3 Metodologia
En primer lugar, importamos la capa ráster proporcionada, la cual contenía una única banda en escala de grises degradada. Para visualizar correctamente los resultados, modificamos la simbología de la capa, dividiéndola en 24 clases. La primera clase correspondía al fondo, por lo que se le asignó una transparencia total, ya que no nos interesaba que se mostrara en el análisis. Posteriormente, para evitar trabajar con la totalidad del ráster que cubre toda Europa —lo cual implicaría un procesamiento muy pesado al convertirlo a vectorial (shapefile)— decidimos recortar el ráster a una extensión aproximada de Galicia.

Para ello, utilizamos la herramienta:
Barra de herramientas > Ráster > Extracción > Cortar ráster por extensión En la ventana emergente, seleccionamos la opción "Draw on Map Canvas", lo que nos permitió delimitar de forma manual la zona de interés: Galicia. Este recorte se aplicó a la capa de incendios forestales. Una vez recortada la capa ráster a la región de Galicia, procedimos a convertirla al formato vectorial (shapefile) con el fin de poder aplicar análisis espaciales, como la intersección. Para ello, utilizamos la herramienta:
Barra de herramientas > Ráster > Conversión > Poligonizar Una vez finalizado el proceso de poligonización, obtuvimos una capa vectorial con múltiples polígonos. Como nuestro objetivo era centrarnos exclusivamente en Galicia, fue necesario recortar los polígonos de incendios y de áreas protegidas que se extendían más allá de sus límites (en zonas fronterizas con Portugal, Asturias y Castilla y León).

Para realizar estos recortes, seguimos los siguientes pasos:

1. Visualización clara de límites: Aplicamos transparencia a las capas vectoriales para poder identificar con precisión los límites geográficos de Galicia.
2. Selección manual: Seleccionamos los polígonos que se extendían fuera de Galicia.
   
   
3. Edición de geometrías: Utilizamos la herramienta: Barra de herramientas > Editar > Editar geometría > Dividir objetos espaciales, lo que nos permitió dividir los polígonos según el contorno de Galicia y eliminar las porciones no deseadas.

Este procedimiento se aplicó por separado a cada capa vectorial (incendios y áreas protegidas), asegurando que los datos analizados se limitaran exclusivamente al territorio gallego. Nuestro objetivo principal en esta fase del proyecto era obtener el área forestal que no había sido afectada por los incendios, con el fin de, posteriormente, realizar un análisis porcentual que nos permitiera identificar qué zonas habían sido más dañadas.

Para ello, utilizamos la herramienta de diferencia espacial entre la capa vectorizada de masa forestal y la capa de incendios. Esta operación nos devolvió un polígono que representaba las zonas forestales no quemadas.

Posteriormente, mediante procesos de análisis adicionales (como la comparación entre superficies originales y resultantes), pudimos calcular el área afectada por los incendios y establecer percentiles de daño, lo cual permitió determinar las regiones con mayor pérdida de masa forestal.

Durante este procedimiento, surgió un inconveniente técnico: una geometría no estaba correctamente cerrada, lo que impedía realizar la operación de intersección. Para solucionarlo, validamos las geometrías con la herramienta:
Barra de herramientas > Vectorial > Herramientas de geometría > Comprobar validez Una vez identificado el polígono problemático gracias a la herramienta de validación, procedimos a comprobar todas las geometrías de la capa vectorial. El análisis indicó que existían dos geometrías inválidas. Para resolver este problema, utilizamos la herramienta "Corregir geometrías", accesible desde el buscador de herramientas. Posteriormente, volvimos a ejecutar la validación y confirmamos que todas las geometrías eran ahora válidas. Con las geometrías corregidas, procedimos a ejecutar la operación de intersección: Barra de herramientas > Vectorial > Herramientas de geoproceso > Intersección Para continuar con el análisis, accedimos a la tabla de atributos de esta nueva capa. Utilizamos la calculadora de campo, aplicando la función $área, con el objetivo de calcular la superficie de cada polígono afectado. Esta operación nos permitió cuantificar con precisión las áreas protegidas dañadas por incendios. Para calcular los percentiles de afección, consideramos que la forma más eficaz de trabajar era exportar los datos a formato Excel, facilitando así su análisis estadístico (véase Anexo II).

Una vez obtenidos los resultados en formato Excel, procedimos a importarlos nuevamente al proyecto para poder vincularlos a la capa de zonas protegidas. Decidimos realizar la unión sobre esta capa porque consideramos que era más representativa en su totalidad, evitando trabajar con polígonos fragmentados, los cuales podrían dificultar la asignación correcta de percentiles debido a su reducido tamaño y dispersión.

Para ello, importamos el archivo Excel al proyecto y, en la caja de capas, seleccionamos la capa de zonas protegidas. Hicimos clic derecho sobre ella y accedimos a Propiedades. Dentro de la ventana emergente, nos dirigimos al apartado Uniones, donde configuramos la vinculación con el archivo Excel. Establecimos como clave de unión el campo sitecode, tanto en la capa como en el archivo Excel, para asegurar la correcta asociación de los datos entre ambas tablas. Esta operación nos permitió disponer de una única capa vectorial con toda la información integrada, lista para aplicar la simbología en función de los percentiles de afección y continuar con el análisis espacial.

Para finalizar, dentro del mismo cuadro de propiedades de la capa de zonas protegidas, accedimos al apartado Simbología, donde configuramos la representación cartográfica en función del percentil calculado previamente. De esta forma, pudimos visualizar de manera clara y diferenciada qué zonas protegidas habían sido más afectadas por los incendios según su grado de daño, completando así todo el proceso metodológico.

5.2 ANEJO II. EXPLICACIÓN EXCEL

Durante la metodología del proyecto, parte de los cálculos sobre capas vectoriales acerca del estado de pérdida de los bosques de Galicia por incendios se han realizado mediante Microsoft Excel. Para ello, se ha creado un archivo de Microsoft Excel dedicado, el cual funciona según el siguiente procedimiento:

1. Se han exportado las tablas de atributos de las capas vectoriales [Red Natura 2000] y [zonas protegidas no afectadas] desde QGIS hasta Microsoft Excel, por medio de archivos en formato “.xlsx”.
2. Mediante formulación de Microsoft Excel, se ha creado un “procesador” de cálculo tal que permite filtrar por columnas las áreas de las distintas geometrías poligonales según el nombre de su zona protegida de la Red Natura 2000 asociada y el tipo de dicha zona protegida.
3. Una vez filtradas las distintas áreas, para cada una de las distintas zonas protegidas de la Red Natura 2000, se obtiene su porcentaje de superficie no afectada por incendios respecto superficie total.
4. Se transponen los resultados de porcentajes de áreas afectadas desde fila a columna, y se calculan los valores a importar de nuevo en QGIS para elaborar la cartografía final aplicando categorización.
5. Finalmente, se importan los datos generados en columnas a QGIS, obteniendo con ello en QGIS los porcentajes de zonas de la Red Natura 200 afectados por incendios forestales en la capa [Porcentaje de área afectada]

5.3 ANEJO III. RESULTADOS

1. Cartografía temática de espacios de la Red Natura 2000 clasificados como zonas ZEPA, ZEC/LIC, y ZEPA/ZEC/LIC.
   
   
2. Cartografía temática de la distribución de incendios forestales en Galicia entre los años 2001 y 2023.
   
   
3. Cartografía específica del porcentaje de áreas afectadas por incendios forestales entre los años 2001 y 2023 de las zonas ZEPA, ZEC/LIC, y ZEPA/ZEC/LIC de la Red Natura 2000 en Galicia.
   
   
4. Tabla de identificación, caracterización y estado de las zonas ZEPA, ZEC/LIC, y ZEPA/ZEC/LIC de la Red Natura 2000 en Galicia (información complementaria al resto de cartografía).

Zonas Red Natura 2000 en Galicia

Incendios Forestales en Galicia (2001-2023)

Porcentajes de afeccion por incendios a la Red Natura 2000 en Galicia

6 Bibliografía

Instituto Galego de Estatística. (s.f.). IGE - Banco de Datos. https://www.ige.gal/igebdt/selector.jsp?COD=238

Ministerio para la Transición Ecológica y el Reto Demográfico. (s.f.). Incendios forestales - Información disponible. Gobierno de España. https://www.miteco.gob.es/es/biodiversidad/servicios/banco-datos-naturaleza/informacion-disponible/incendios-forestales.html

Global Forest Watch. (s.f.). Mapa interactivo de cobertura forestal y fuegos. https://www.globalforestwatch.org/map/?map=eyJkYXRhc2V0cyI6W3siZGF0YXNldCI6InRyZWUtY292ZXItbG9zcy1maXJlcyIsIm9wYWNpdHkiOjEsInZpc2liaWxpdHkiOnRydWUsImxheWVycyI6WyJ0cmVlLWNdmVyLWxvc3MtZmlyZXMiXX0seyJkYXRhc2V0IjoicG9saXRpY2FsLWJvdW5kYXJpZXMiLCJsYXllcnMiOlsiZGlzcHV0ZWQtcG9saXRpY2FsLWJvdW5kYXJpZXMiLCJb2xpdGljYWwtYm91bmRhcmllcyJdLCJvcGFjaXR5IjoxLCJ2aXNpYmlsaXR5Ijp0cnVlfV19&mapMenu=eyJtZW51U2VjdGlvbiI6ImRhdGFzZXRzIiwiZGF0YXNldENhdGVnb3J5IjoiZm9yZXN0Q2hhbmdlIn0%3D

Hansen, M. C., & Giglio, L. (s.f.). Global Fire Atlas Data (2001–2023). University of Maryland. https://glad.umd.edu/users/Alexandra/Fire_GFL_data/2001-23/

Ministerio para la Transición Ecológica y el Reto Demográfico. (s.f.). Red Natura 2000 - Información disponible. Gobierno de España. https://www.miteco.gob.es/es/biodiversidad/servicios/banco-datos-naturaleza/informacion-disponible/rednatura_2000_desc.html

Xunta de Galicia. (s.f.). Centro de Descargas - Información Geográfica. https://mapas.xunta.gal/es/mapas/informacion-geografica/centro-de-descargas