- Versión: v.2021-2022
- Asignatura (grado): Ecología (CCAA)
- Autor: Curro Bonet-García (fjbonet@uco.es)
- Duración: Dos sesiones de tres horas cada una.
Esta práctica tiene los siguientes objetivos:
Disciplinares: Están relacionados con la adquisición de competencias y saberes propios de la ecología.
Instrumentales: Están relacionados con la adquisición de competencias en el manejo de herramientas potencialmente útiles en ecología y en otros ámbitos.
El abordaje de los objetivos anteriores se realizará a través del estudio y análisis de una pregunta concreta. Así pues, esta práctica sigue la técnica docente denominada "aprendizaje basado en problemas". Aprenderemos todo lo enumerado anteriormente mediante la resolución de una pregunta concreta que se formulará después de contextualizar ecológicamente el objeto de esta práctica: los pinares de repoblación. Esta presentación muestra con detalle lo enumerado a continuación.
Los pinares de repoblación son masas forestales densas, monoespecíficas y coetáneas que se implantan en el territorio por una serie de razones:
La implantación de los pinares de repoblación en Andalucía consiguió frenar la intensa erosión que habí en nuestros montes a principios del siglo XX. Sin embargo, las repoblaciones se hicieron con una densidad excesiva, lo cual provoca una intensa competencia intraespecífica. Esto, junto con la escasez de tratamientos selvícolas, hace que ahora estas masas generan ciertos problemas socioecológicos.
Durante esta práctica trataremos de analizar los problemas anteriores y, sobre todo, intentaremos poner en práctica una metodología para minimizarlos: la naturalización de los pinares. Esta técnica de manejo consiste en reducir la densidad de los pinares minimizando la competencia intraespecífica. Esto provoca que dichos pinares sean susceptibles de ser "invadidos" por vegetación natural propia de la zona, lo cual aumenta la biodiversidad, reduce el riesgo de incencio, etc. De manera más específica, en esta práctica trataremos de abordar la siguiente pregunta:
¿Qué lugares de Sierra Nevada ocupados actualmente por pinares de repoblación pueden recibir tratamientos selvícolas (y reducir la competencia intraespecífica) para promover el crecimiento de la vegetación natural?
El abordaje de la pregunta anterior se sustenta en los siguientes conceptos ecológicos:
Esta práctica se desarrolla en Sierra Nevada. La razón fundamental es que en ese macizo montañoso abundan los pinares de repoblación. En concreto hay unas 40.000 Has de pinares de repoblación, que suponen aproximadamente la mitad de toda la superficie boscosa de este espacio protegido. Otra razón importante es que, con un poco de suerte, haremos una salida de campo a Sierra Nevada ;)
Este mapa muestra la ubicación de Sierra Nevada. Los pinares de repoblación son las formaciones de color verde oscuro que hay en la vertiente norte de la montaña.
Para abordar la pregunta planteada más arriba es necesario tener en cuenta una serie de variables ambientales que se describen a continuación:
La densidad del pinar es el principal factor a considerar, ya que esta variable es la que condiciona la intensidad de la competencia intraespecífica. Dado que nuestro objetivo es reducir dicha competencia, consideraremos que los sitios más adecuados para recibir tratamientos selvícola son aquellos que tienen mayor densidad.
Todo lo que vamos a trabajar en esta práctica tiene una componente espacial evidente. Así que trabajaremos con mapas. En este caso, el mapa de densidad procede de una imagen de satélite. Concretamente ha sido generada por una infraestructura de investigación denominada Copernicus, que suministra multitud de información ambiental procedente de sensores remotos. Concretamente vamos a usar una imagen denominada TCD (Tree cover density), que muestra la densidad de árboles expresada en porcentaje. La imagen cubre buena parte de Europa, así que yo la he recortado y reproyectado para que se ajuste a lo que necesitamos.
La variable "densidad del pinar" ha de ser transformada en un criterio espacializable. En este caso asumimos que a más densidad mayor idoneidad para recibir tratamientos forestales. Es decir, hay una relación directa entre la distribución de la variable (densidad) y el criterio (aptitud desde el punto de vista de la densidad). Dado que la densidad se expresa en porcentaje (de 0 a 100%), es muy fácil transformarla en un mapa de aptitud desde el punto de vista de la densidad. Lo veremos más adelante.
El proceso mediante el cual los pinares de repoblación son "invadidos" por vegetación natural depende en buena medida de la disponibilidad de propágulos de otras especies. Hay varios trabajos que justifican esta afirmación (1, 2). Un pinar que esté lejos de una mancha de vegetación natural que suministre semillas, será invadido con menor probabilidad. Es decir, la distancia de cada punto del territorio ocupado por pinares a las manchas de vegetación natural más cercana, es importante para nuestro objetivo. En este caso la variable es "distancia de cada pinar a la mancha de vegetación natural más cercana". Su transformación a criterio implica tener en cuenta que cuanto más distancia haya, menos adecuado será un punto del territorio para albergar tratamientos selvícolas. Es decir, en este caso, la transformación de variable a criterio requiere usar una función invertida (a valores más altos de la variable corresponden valores más bajos de idoneidad o aptidud).
Como veremos en las siguientes secciones, obtendremos el mapa de distancias a partir de un mapa de usos y coberturas vegetales del suelo.
Las plantas toman del suelo el agua que necesitan para vivir. Por ello es muy importante que el suelo tenga cierta cantidad de humedad. Esto es especialmente importante en ambientes mediterráneos como el nuestro, donde la sequía estival reduce notablemente la disponibilidad de agua en el suelo. Tendremos en cuenta esta variable de una manera muy rudimentaria que tiene en cuenta únicamente la topografía. Los lugares cóncavos tienden a acumular más agua que los convexos. Así, el grado de convexidad-concavidad del terreno nos ayuda a simular cómo se almacenaría el agua en el suelo. En estos casos se dice que usamos una variable fácil de obtener (concavidad-convexidad del territorio) como subrrogado de otra más compleja de calcular (humedad del suelo). Afortunadamente la REDIAM (Red de Información Ambiental de Andalucía) suministra una capa que muestra la distribución de la curvatura de toda la región. Puedes descargarla aquí. He recortado esta imagen para que se ajuste a nuestra zona de estudio. Esta capa contiene valores que van desde -0.37 hasta 0.32. Los valores más altos corresponden con lugares más convexos, es decir, sitios en los que se acumula menos agua (cumbres de montañas, divisorias de aguas). Por contra, los valores más pequeños corresponden con lugares cóncavos (cauces, barrancos, etc.). Al igual que en el caso anterior, a valores más altos de la variable corresponden valores más bajos de aptitud. Es decir, a menor índice de humedad del suelo mayor idoneidad desde el punto de vista de los objetivos de este trabajo.
La siguiente figura muestra una versión simplificada del flujo de trabajo que seguiremos durante el desarrollo de la actividad. Un flujo de trabajo es una representación gráfica de una serie de procedimientos ordenados secuencialmente y que van encaminados a lograr un objetivo concreto. Esta es nuestra primera versión del flujo de trabajo

El siguiente esquema es un documento dinámico creado con una herramienta llamada draw.io . Se trata de una versión detallada del flujo de trabajo que seguiremos. Puedes utilizarlo para realizar la práctica en casa o en clase. Aquí puedes descargar el mismo esquema que ves abajo. En caso de que lo hayas descargado, lo podrás abrir desde draw.io.
Para empezar, recordemos la pregunta planteada: Queremos identificar zonas del territorio que estén ocupadas por pinares de repoblación que sean susceptibles de recibir tratamientos forestales de reducción de la competencia intraespecífica (aclarado de pinares) con objeto de maximizar la regeneración de la vegetación natural. Es decir, de alguna manera queremos identificar las áreas del territorio que reúnan unas condiciones más adecuadas para soportar dichos tratamientos. O dicho de otra manera más, vamos a asignar una "puntuación" a los distintos pinares en función de que satisfagan o no una serie de criterios que consideramos importantes.
Así que el primer paso es identificar las variables ambientales y los criterios de manejo con los que vamos a trabajar. En nuestro caso serán los siguientes:
Aunque se irán detallando y describiendo a lo largo del guión, aquí tienes la lista del material que necesitas para ejecutar esta práctica:
Mapa de usos y coberturas vegetales del suelo de Sierra Nevada. Usaremos esta capa para generar el mapa de distancias de los pinares a las manchas de vegetación natural más cercana. Este video describe su estructura de datos: MUCVA_25_multi_snevada.zip
Mapa de densidad de los pinares, expreasada en porcentaje: TCD_pinares_23030.tif
Mapa de humedad potencial del suelo basada en la topografía: cti_pinares_23030.tif
Mapa de distribución de los pinares de repoblación en Sierra Nevada: pinares_repoblacion_23030.tif
Las siguientes acciones se describen con la misma numeración que se muestra en el flujo de trabajo descrito más arriba.
Para generar este mapa necesitamos contar con la distribución de los pinares de repoblación y con la de las formaciones vegetales que actuarán como donadoras de propágulos. Obtendremos ambos mapas a partir del mapa de usos y coberturas vegetales de Andalucía, generado por la REDIAM. En este enlace puedes descargar el mapa completo (Andalucía). Y aquí tienes un video en el que te cuento cómo está estructurado. Para esta práctica mejor usa la capa que se muestra más arriba.
De manera resumida haremos lo siguiente: Crearemos un campo nuevo en el mapa de usos y coberturas anterior y asignaremos los valores 1 a todos los polígonos que tengan pinares, mientras que pondremos el valor 2 a todos los que sean considerados como fuentes de semillas.
El algoritmo de QGIS que permite calcular el mapa de distancias necesita un raster como capa de entrada. Así que rasterizaremos el fichero de formas anterior. Una vez hecho esto podremos calcular la distancia entre todos los píxeles ocupados por pinares (1) y su mancha de vegetación natural más cercana (2).
Vamos con el paso a paso:
x"DES_UC07" = 'FOR. ARBOL. DENSA: CONIFERAS' OR "DES_UC07" = 'FOR. ARBOL. DENSA: QUERCINEAS+CONIFERAS' or "DES_UC07" = 'MATORRAL DENSO ARBOLADO: CONIFERAS DENSAS'
xxxxxxxxxx"DES_UC07" = 'FOR. ARBOL. DENSA: QUERCINEAS' or "DES_UC07" ='MATORRAL DENSO' or "DES_UC07" ='MATORRAL DENSO ARBOLADO: OTRAS FRONDOSAS' or"DES_UC07" = 'MATORRAL DENSO ARBOLADO: QUERCINEAS DENSAS' or"DES_UC07" = 'MATORRAL DENSO ARBOLADO: QUERCINEAS DISPERSAS' or"DES_UC07" = 'MATORRAL DISP. ARBOLADO: QUERCINEAS+CONIFERAS' or "DES_UC07" = 'MATORRAL DISP. ARBOLADO: QUERCINEAS. DENSO' or"DES_UC07" = 'MATORRAL DISP. ARBOLADO: QUERCINEAS. DISPERSO'
(3) Al igual que antes, haz que estos polígonos seleccionados tomen el valor 2 en el campo dist_targe. Guarda los cambios.
(4) Ahora rasterizamos el fichero de formas anterior con la opción rasterizar de QGIS (menú raster -> conversion -> rasterizar). Aplicamos los siguientes valores a los parámetros necesarios:
(5) Creamos el mapa de distancia usando el algoritmo llamado proximity raster (GDAL) en QGIS. Necesitamos añadir los siguientes parámetros:
O expresado en el lenguaje de QGIS:
xxxxxxxxxx gdal_proximity.py -srcband 1 -distunits GEO -values 2 -nodata 0.0 -ot Float32 -of GTiff "/tu ruta/dist_target.tif" "/tu ruta/distancia1.tif"
xxxxxxxxxx 1/(100-7300.68505859375)*"dist_nat@1"- (7300.68505859375/(100-7300.68505859375)) 
xxxxxxxxxx 1/(-0.370579987764359-0.329879999160767)*"cti_pinares_23030@1"- (0.329879999160767/(-0.370579987764359-0.329879999160767)) El análisis multicriterio es una técnica muy sencilla que permite conciliar en un mismo mapa criterios que pueden ser contradictorios. En nuestro ejemplo tenemos tres criterios y se trata de unificarlos en un único mapa que asigne un valor de aptitud global a cada punto ocupado por pinares de repoblación. Para agregar los tres criterios empezaremos asignando un peso a cada uno de ellos. La suma de los pesos ha de ser 1. El criterio que tenga más peso contribuirá en mayor medida al mapa final. Procedemos en dos pasos:
(10) Asignación de pesos a cada variable. Por ejemplo:
(11) Integración de las variables y sus pesos. El proceso se hace fácilmente con la calculadora de mapas. El resultado final es la suma del producto de cada criterio (capa apt) por su peso. Introducimos la siguiente ecuación en la calculadora raster:
xxxxxxxxxx ("apt_densidad@1"*0.6)+("apt_cti@1"*0.1)+("apt_distancia@1"*0.3) Llamamos apt_final.tif a la capa resultante.
Uno de los problemas del análisis multicriterio es que ocurre una compensación de criterios. Si una variable tiene un valor muy alto en un lugar determinado, puede que el resultado final en ese punto sea algo aunque su peso sea bajo. Esto puede hacer que lugares no adecuados sean etiquetados como sí adecuados.
(12) El mapa resultante apt_final.tif tiene valores númericos que van de 0 a 1. A partir de este mapa debermos seleccionar aquellos píxeles que tengan una aptitud mayor. El resto los descartaremos porque no reunen los requisitos que hemos impuesto. Para hacer esta selección aplicaremos una operación muy común en análisis raster: reclasificación. Consiste en reducir la diversidad de valores de un raster asigando nuevos en función de un rango. Por ejemplo, asignaremos el valor de 1 a todos los píxeles que tengan una aptitud igual o mayor de 0.8. Para hacer esto, construimos una tabla de reclasificación.
Buscamos el algoritmo "reclassify by table" en el buscador de procesos de QGIS. Añadimos los siguientes parámetros:
| Desde (Mínimo) | hasta (Máximo) | Valor nuevo |
|---|---|---|
| 0 | 0.9 | 0 |
| 0.9 | 1 | 1 |
El siguiente video muestra (en inglés) la secuencia de acciones necesarias para completar la práctica.
Dado que esta práctica tiene dos sesiones, en el último día deberás entregar un breve informe que contenga la siguiente información:
Lista ordenada de los pasos que has ido ejecutando con el ordenador para obtener el mapa final con los lugares más adecuados para reducir la competencia intraespecífica. Entre estos pasos debe de estar la combinación de pesos que consideres más adecuada para satisfacer nuestros objetivos. Deberás ir completando esta parte del informe conforme vayamos trabajando juntos en clase. Puedes describir cada paso como quieras. Te sugiero algo así:
Paso 1: [nombre del paso]
Resultado obtenido. Deberás incluir en el informe un mapa de una zona que elijas de Sierra Nevada. Además, deberás describir dicha zona usando términos ecológicos relacionados con la competencia intraespecífica. Por ejemplo: "la zona seleccionada se caracteriza por estar en una ladera de orientación norte. Hay un pinar de repoblación denso que no tiene cerca ninguna mancha donadora de semillas. Sin embargo, el suelo parece adecuado para almacenar agua, por lo que se ha considerado que es un lugar adecuado para realizar tratamientos forestales".
El informe se evaluará según la siguiente rúbrica. Sube el material a esta página de Turnitin.
| Criterio de evaluación | 0 puntos | 1 punto | 2 puntos | 3 puntos | 4 puntos | 5 puntos |
|---|---|---|---|---|---|---|
| Selección de variables: Se refiere a si las variables elegidas para tu análisis se ajustan a la lógica del problema planteado. | No entrega nada o es incomprensible. | El conjunto de variables elegido no se alinea con el objetivo del trabajo. | No se tienen en cuenta los grupos de variables propuestas. | Incluye todas las variables sin justificar la decisión. | Excelente selección de variables. | Además de lo anterior, incorpora variables nuevas no propuestas inicialmente en el ejercicio. |
| Descripción de los procesos realizados: Se refiere a lo bien que ha sido descrito el proceso que permite generar el resultado final. | No entrega nada o es incomprensible. | Es imposible repetir el proceso con la información suministrada. | Se entiende el proceso y se podrían repetir algunas partes del mismo, pero el conjunto no contiene la secuencia completa. | Buena descripción. Sería posible reproduccir el trabajo realizado. | Excelente descripción del flujo de trabajo. | Además, la presentación es brillante. |
| Descripción de la zona seleccionada: Este criterio evalúa en qué medida describes bien la zona elegida después del trabajo. | No entrega nada o es incomprensible. | La descripción no contiene menciones a la competencia ni a ningún concepto ecológico relevante. | Los conceptos teóricos en los que se basa el ejercicio aparecen tangencialmente. | Incorpora adecuadamente la competencia intraespecífica y las demás variables consideradas. | Además de lo anterior, mencionas adecuadamente las relaciones de la ecología con aspectos socioeconómicos. | Excelente integración de marcos conceptuales diferentes. Además de los requeridos has incluido otros. |
| Legibilidad: hace referencia a lo bien escrito que está el texto y a su legibilidad | No entrega nada o es incomprensible. | Apenas entiendo lo que has escrito | He tenido que reinterpretar casi cada frase para entenderlo | Se entiende bien todo, pero el texto no es fluido | Muy buena redacción. La lectura fluye fácilmente, cual novela. | Impecable estilo de escritura. |
Aquí puedes ver los resultados que, por ahora, arroja esta encuesta sobre grado de satisfacción de los objetivos planteados. Por favor, rellénala si no lo has hecho ya. Se tardan 10 segundos y sus resultados son de gran utilidad para mejorar los contenidos cada año. Gracias.


Además, algunos estudiantes hicieron los siguientes comentarios: