1. IntroducciónTop
El pinsapo (Abies pinsapo. Boiss) es un endemismo Ibérico (Bético) cuyas poblaciones se localizan en la Serranía de Ronda[1]. Forma parte de la estirpe de los abetos circum-mediterráneos, especies con un antepasado común que evolucionaron de manera homóloga desde el Terciario en función de las pulsaciones climáticas del Cuaternario, aprovechando microambientes favorables en las principales cadenas montañosas que bordean la Cuenca Mediterránea (Alizoti et al., 2011; Linares, 2011).
En al-Ándalus se empleaba el fitónimo «Juayha» para referirse al pinsapo, un término que también se utilizaba para designar a los abetos rifeños (Martínez Enamorado, 2013). También reconocido como abeto –concretamente «abete»– a finales del siglo XVIII (Bowles, 1782), en la zona de Málaga adoptó la denominación de pinsapo a partir del fitónimo «pinabeto», que parece derivar de pinus (pino) y sapinus (abeto), mientras que en la zona de Cádiz se denominaba pino, de ahí el topónimo de Sierra del Pinar (Grazalema). Aunque ya fue descrito por Simón de Rojas Clemente y Rubio, su identificación científica corresponde a Edmund Boissier, muy probablemente influido por el trabajo del botánico español (Cabezudo et al., 2010).
Es un árbol de porte piramidal que puede alcanzar hasta los 30 m de altura, con un tronco de hasta 1 m de diámetro, si bien en condiciones naturales no suele sobrepasar los 20 m (López, 2007). Florece entre abril y mayo y forma las piñas entre junio y octubre, un periodo relativamente largo, quizás como respuesta a la sequía estival, diseminando las semillas a continuación (Pérez-Latorre & Cabezudo, 2012). Desde el punto de vista ecofisiológico, las principales diferencias con respecto al otro abeto que habita en la península Ibérica, Abies alba, son su mejor adaptación a la sequía y su menor tolerancia al frío (Peguero-Pina et al., 2011; Sancho-Knapik et al., 2014). La última revisión sobre la distribución actual de los abetos circum-mediterráneos (San-Miguel-Ayanz et al., 2016), muestra que la mayor parte de ellos persisten en poblaciones pequeñas y medianas en las principales alineaciones montañosas, localizadas en lugares que presentan unas elevadas precipitaciones invernales y una prolongada sequía estival. En general, los abetos tienden a formar rodales puros cuando se encuentran en hábitats óptimos. Hoy en día, con la excepción del abeto blanco (Abies alba) y el abeto del Cáucaso (Abies nordmanniana), los abetos circum-mediterráneos no presentan un gran interés comercial. Sin embargo, la conservación de los abetos circum-mediterráneos ha cobrado una importancia capital, especialmente en aquellas especies con áreas de distribución muy restringidas, con la creación de reservas protegidas y programas de conservación.
Según la Unión Internacional para la Conservación de la Naturaleza (IUCN, International Union for Conservation of Nature), el pinsapo es considerado una especie en peligro de extinción por la existencia de poblaciones fragmentadas en las que se han observado, inferido y / o proyectado diversos procesos de decaimiento. Según los mismos criterios, la especie se encuentra incluida en la Lista Roja de la Flora Vascular Española y, conforme la Directiva 92/43/CEE, dentro del Hábitat de Interés Comunitario 9520 «Abetales de Abies pinsapo Boiss». El pinsapo está catalogado «En Peligro de Extinción» (EN) dentro del Catálogo Andaluz de Especies Amenazadas del Listado Andaluz de Especies Silvestres en Régimen de Protección Especial (LAESPE) emanado de la Ley 8/2003 de la Flora y Fauna Silvestre de Andalucía.
La Consejería de Medio Ambiente de la Junta de Andalucía viene realizando acciones de manejo del pinsapar orientadas a garantizar la conservación del hábitat y reducir los factores de riesgo. Entre ellas destacan: reforestaciones experimentales en áreas antes habitadas por la especie o afectadas por incendios, resguardo de algunas zonas con cerramiento ganadero, protección individual de aquellos individuos que se hallaban aislados, empleo de tratamientos fitosanitarios para fortalecer las poblaciones contra las plagas, eliminación de pies muertos de pinsapo también para el control de incendios y plagas, etc. También se han desarrollado actuaciones de conservación «ex situ» dentro de la Red Andaluza de Jardines Botánicos y más concretamente a través de la Red de Viveros mediante el Laboratorio de Propagación Vegetal y a través del Banco de Germoplasma Vegetal Andaluz.
La mayor parte de estas acciones se iniciaron a mediados de la década de los noventa del siglo XX, se incrementaron en la primera década del siglo XXI con la puesta en marcha de la Red de Equilibrios Biológicos en ecosistemas con presencia de Abies pinsapo (Navarro-Cerrillo & Calzado Martínez, 2004), y finalmente se han formalizado en el Plan de Recuperación del Pinsapo aprobado por la Junta de Andalucía. Se han diagnosticado diferentes amenazas que afectan a los pinsapares Béticos: los incendios forestales, el sobrepastoreo, el aislamiento de las poblaciones, la modificación y fragmentación del hábitat, la incidencia de los patógenos, el calentamiento global y sus efectos en el clima de la zona, la hibridación del pinsapo, la contaminación atmosférica y la propia estructura demográfica de las poblaciones (López Quintanilla, 2013; Caudullo & Tinner, 2016).
El proyecto de Ley de Medidas Frente al Cambio Climático de la Comunidad Autónoma de Andalucía (Consejería de Medio Ambiente y Ordenación del Territorio, 2017) recoge los principales impactos derivados del calentamiento global en la región, así como un programa de mitigación y adaptación ante los nuevos escenarios climáticos previstos. Como hemos comentado anteriormente en relación con el Plan de Recuperación del Pinsapo, buena parte de estas medidas ya se han puesto en práctica en el campo del manejo de la biodiversidad.
De acuerdo con los escenarios de cambio climático previstos por el 5º Informe del IPCC (IPCC, 2013), se auguran unas expectativas poco favorables para la supervivencia de Abies pinsapo hacia el horizonte del año 2100 (Gutiérrez-Hernández, 2018). Sin embargo, la estrategia de conservación y recuperación del pinsapo impulsada por la Consejería de Medio Ambiente de la Junta de Andalucía pretende, a través de una gestión activa, la adaptación de la especie a este escenario cambiante.
Dado que los efectos del cambio climático son cada vez más evidentes a escala global (Mauritsen & Pincus, 2017), y ya se están manifestando tendencias a escala local, incluso en los que respecta a las precipitaciones (Ruiz Sinoga et al., 2011), cabría esperar un posible efecto negativo sobre la dinámica de las formaciones forestales dominadas por Abies pinsapo. Sin embargo, también cabría esperar que todo el esfuerzo realizado por la Consejería de Medio Ambiente hubiera arrojado resultados esperanzadores en relación con la recuperación de los ecosistemas del pinsapo. Nuestra hipótesis apunta en esta segunda dirección, por lo menos a corto plazo.
Ante este complejo panorama, la monitorización de la vegetación por medio de la teledetección se presenta como uno de los sistemas más adecuados para registrar objetivamente la naturaleza de los cambios desde un punto de vista global y sistemático (Jones & Vaughan, 2010).
El objetivo general de este trabajo es estudiar la dinámica temporal, interanual y estacional, del vigor de la vegetación y fenología en las superficies forestales dominadas por Abies pinsapo. Se trata de identificar los cambios ocurridos a lo largo de los últimos quince años en los pinsapares de Andalucía y detectar las posibles tendencias a partir de información obtenida por el sensor MODIS, tomando como unidad de medida el NDVI, con frecuencia mensual, calculado desde el año 2002 hasta el año 2016, a una resolución espacial de 250 m, para componer una serie temporal de 180 imágenes de satélite. Serie que consideramos lo suficientemente larga y consistente para realizar un primer diagnóstico sobre la evolución del pinsapar en los albores del siglo XXI.
2. MetodologíaTop
2.1. Área de estudio
Las poblaciones de Abies pinsapo se distribuyen en su mayor parte por las sierras de las Nieves, Bermeja y Blanca en la provincia de Málaga, y la sierra de Grazalema en Cádiz, conjunto que se corresponde con el ámbito geográfico-histórico de la Serranía de Ronda (Figura 1).
|
Figure 1: Study area. Source: own elaboration. Based on DERA-IECA. |
La especie habita entre los 700 y 1700 m de altitud; excepcionalmente fuera de este rango, siendo más abundante en torno a los 1200 m de altitud, en laderas generalmente empinadas (35% de pendiente media) y orientadas al norte (>70% de la superficie), y en lugares donde las lluvias anuales superan los 800 mm (Gutiérrez-Hernández et al., 2017). Aunque la especie asoma frecuentemente sobre terrenos calizos, realmente es indiferente del sustrato, y es posible encontrarla sobre peridotitas, pizarras o gneises. En el límite inferior de su distribución, los pinsapares se entremezclan con formaciones mediterráneas de Quercus rotundifolia, Quercus faginea, Quercus suber, Pinus halepensis o Pinus pinaster. Por encima de los 1.200 m de altitud, el pinsapo tiende a formar bosques con estrato arbóreo monoespecífico (Valladares, 2009).
Por tanto, el pinsapo es un abeto claramente vinculado a los ecosistemas de montaña del suroeste de la península Ibérica, que aparece puntualmente en diferentes puntos de la Serranía de Ronda como consecuencia de una serie de vicisitudes biogeográficas; habita dentro de un rango de tolerancias bioclimáticas relativamente restringido, entre la media-alta montaña y la media montaña mediterránea, donde predominan los inviernos frescos y húmedos, manifestando algunos registros propios de los bosques de coníferas templado boreales en los pisos superiores, pero dentro del contexto de un clima mediterráneo que impone la severidad de la sequía y una intensa evapotranspiración en el periodo estival y determina unas transiciones forestales que tienden hacia la mediterraneidad más propia de los pisos inferiores, donde el pinsapo se entremezcla con especies más termófilas (Figura 2).
|
Figure 2: (1) Centenarian forest of Spanish fir in high places; (2) Spanish fir branches; (3) Spanish fir forest on peridotite rocks; (4) Spanish fir forest on limestone rocks; (5) Snow-covered forest of Spanish fir in high places (6) Mixed forest: cork oak, Portuguese oak, Pinus haplependis, and Spanish fir. Photos: Paco G. Portillo (2,4,6), Ignacio Diaz Triviño (1,5), Javier Recio-Moreno (3). |
2.2. Pre-procesamiento de imágenes de satélite del sensor MODIS
Empleamos un total de 180 imágenes MODIS para representar una serie con frecuencia mensual de 15 años de duración, entre enero de 2002 hasta diciembre de 2016, donde cada imagen registró el NDVI mensual calculado con las bandas 1 (rojo) y 2 (infrarrojo cercano) del sensor MODIS a bordo de los satélites Terra y Aqua (Rouse et al., 1974), de acuerdo con las especificaciones técnicas del MODIS Science Team (NASA).
La resolución de los datos de entrada se estableció a través de una malla regular de 250 metros, de acuerdo con la resolución espacial de las imágenes de satélite empleadas. Los datos georreferenciados se integraron en un Sistema de Información Geográfica (Grass GIS), donde definimos la misma región, máscara, resolución, alineación y codificación para todas las imágenes, y asimismo proyectamos los geodatos usando el sistema de coordenadas UTM ETRS89 HUSO 30.
Corregimos los píxeles con valores erróneos (nieve, nubes, errores de registro) a través de una adaptación del procedimiento de HANTS - Harmonic Analysis of Time Series (Roerink et al., 2000) usando Terrset (Eastman, 2015). Aunque es posible generar una nueva serie en la que todos los píxeles son reemplazados por sus valores predichos basados en una regresión armónica, sólo aplicamos este procedimiento para rellenar los valores ausentes ajustando una regresión armónica a una ventana de tiempo y utilizando el valor de cada píxel predicho para reemplazar los píxeles sin información. Como resultado del pre-procesamiento, obtuvimos una serie temporal completa, uniforme y coherente, válida para todo tipo de análisis estadístico.
2.3. Estrategia de muestreo
Adoptamos como referencia acerca de la distribución y delimitación de las superficies cubiertas parcial o totalmente por Abies pinsapo, la cartografía oficial del Hábitat de Interés Comunitario 9520 «Abetales de Abies pinsapo Boiss» (Valladares, 2009) producida por la REDIAM (Consejería de Medioambiente, Junta de Andalucía) en cumplimiento de la Directiva 92/43/CEE.
Dadas las limitaciones establecidas por la resolución espacial de las imágenes MODIS, que fijan una escala espacial de captura del fenómeno de 250 m, sobre la base de imágenes de satélite superpusimos una retícula de polígonos de 250 m de lado. En total, 573 polígonos incluyeron presencia del pinsapo. Utilizamos la ortofotografía más actualizada del PNOA, con fecha de captura julio de 2016, para proyectar (y validar, en su caso) los datos de la cartografía de referencia dentro de la malla. Hay que tener en cuenta que el NDVI registrado en cada rejilla de 250 m con presencia de pinsapo, puede informar sobre copas de pinsapos, copas de otros árboles, arbustos o herbáceas e incluso la señal del roquedo. Es por este motivo por el que recurrimos a la fotointerpretación para caracterizar la información contenida en cada polígono con presencia de pinsapo. Así, dentro de cada polígono de 250 m de lado, generamos cuatro teselas de 125 m de lado para caracterizar la información contenida en el conjunto del polígono. Mediante la fotointerpretación también distinguimos la presencia de formaciones arboladas mixtas y la presencia de otras coberturas arbustivas, herbáceas o roquedo. Este trabajo se hizo para determinar con mayor exactitud qué dato aporta realmente el análisis de tendencia estacional sobre cada unidad espacial. Así pudimos discriminar cuándo un píxel informaba exclusivamente sobre copas de pinsapos o sobre otros elementos, incluidas las copas de pinsapo u otras especies. En la Figura 3 se ilustran las zonas muestreadas y un ejemplo de la fotointerpretación basada en las teselas de 125 m de lado incluidas en la malla de 250 m de lado que usamos para caracterizar las formaciones con presencia de pinsapo.
|
Figure 3: Sample design and field data validation based on photo-interpretation. 1. Very dense forest. 2. Dense forest 3. Open forest. 4. Isolated trees. Source: own elaboration. Based on REDIAM and PNOA. |
En definitiva, el muestreo precisó la máscara de análisis de los datos satelitales, mostrando las zonas (píxeles) con presencia de pinsapo: discriminamos áreas, poblaciones, formaciones, etc.
2.4. Análisis estadístico
Una vez que las imágenes de satélite fueron pre-procesadas, analizamos los píxeles correspondientes a la zona de muestreo.
Primero, calculamos medidas de tendencia central y de dispersión del NDVI para caracterizar el estado general de las masas en las principales zonas (resultados en apartado 3.1). Después, realizamos un análisis de tendencia interanual, desestacionalizando la serie y estimando los siguientes parámetros que se indican seguidamente:
- Coeficiente de producto-momento de Pearson (R), entre los valores de cada píxel a través del tiempo y una serie perfectamente lineal (resultados en apartado 3.2). El valor R2 informa sobre la linealidad de la serie, si bien es sensible al ruido en series cortas;
- Prueba de Mann-Kendall o tendencia monótona (resultados en apartado 3.2), más robusta que la anterior, que mide el grado en que una tendencia aumenta o disminuye consistentemente en el tiempo. Un valor de +1 indica que la tendencia aumenta consistentemente y nunca disminuye. Un valor -1 indica todo lo contrario. Y un valor 0 indica que no hay una tendencia consistente. Estimamos la significación estadística (p) del valor de Mann-Kendall obtenido (Neeti & Eastman, 2011).
- Finalmente, aplicamos la técnica de Análisis de Tendencia Estacional (STA, del inglés Seasonal Trend Analysis) desarrollada por ClarkLabs (Eastman et al., 2009, 2013).
El STA calcula una regresión armónica de imágenes anuales seguida de un análisis de Kendall de la amplitud y las fases de cada año generadas por la regresión armónica. La regresión armónica es similar al análisis de Fourier, pero permite la especificación de fechas julianas. Se utilizan dos o cuatro frecuencias para la regresión armónica con el fin de capturar las tendencias dominantes en los datos sin incluir el ruido de alta frecuencia (resultados en apartado 3.3).
Los resultados producen nuevas imágenes con pendiente, intercepto y significación estadística (Mann-Kendall), así como compuestos para fases y amplitudes donde especifican las imágenes medianas desde el primer año hasta el último. De este modo, estimamos la curva tendencia estacional para cada periodo de referencia, basada en la mediana de dos años consecutivos, y calculamos los parámetros fenológicos básicos derivados de los periodos de reverdecimiento o Green up (Green up onset y Green up phase) y senescencia o Green down (Green down onset y Green down phase) (Zhang et al., 2003, 2012). Para el análisis de la fase del Green up y la fase del Green down, el umbral de la curva se estableció en el punto en que la cantidad de verde excede el 40% de toda la trayectoria desde el valor mínimo hasta el máximo (Eastman, 2015). Finalmente, cartografiamos los resultados del análisis de tendencia interanual y dibujamos los resultados del análisis de tendencia estacional.
3. ResultadosTop
3.1. Estado general de las masas
Los resultados sobre el estado de los pinsapares en las zonas se ilustran en la Tabla 1 y Figura 4. La mediana del NDVI para el conjunto de las áreas con presencia de Abies pinsapo fue de 0,68, mínimo en Sierra Blanca-Canucha (0,57) y máximo en Sierra de Grazalema (0,77). Observamos un aumento del NDVI de este a oeste, que fue consistente también cuando sólo consideramos los píxeles con dominancia de pinsapo. No obstante, los valores promedios del NDVI también estuvieron influidos por las características de los recintos o fragmentos de paisaje con presencia de pinsapo, de modo que fragmentos más pequeños y con mayor densidad de bordes, mostraron un NDVI menor por la influencia de la señal del suelo (ruido). En esta dirección, encontramos que en las cotas más elevadas de la Sierra de las Nieves y Sierra de Grazalema el NDVI fue menor.
Tabla 1: Estadísticas básicas del NDVI mensual en zonas con presencia de pinsapo. Fuente: Elaboración propia.
Table 1: Basic statistics in NDVI of the Spanish fir forests. Source: own elaboration.
| Estadístico | S. de las Nieves | S. de Grazalema | S. Bermeja | S. Blanca-Canucha | Conjunto |
|---|---|---|---|---|---|
| n (malla 250 m) | 401 | 125 | 37 | 10 | 573 |
| Mínimo | 0.338 | 0.442 | 0.524 | 0.553 | 0.338 |
| Máximo | 0.798 | 0.859 | 0.776 | 0.688 | 0.859 |
| Mediana | 0.633 | 0.776 | 0.648 | 0.577 | 0.685 |
| Media | 0.602 | 0.735 | 0.657 | 0.602 | 0.664 |
| Desviación estándar | 0.115 | 0.108 | 0.052 | 0.049 | 0.091 |
|
Figure 4: Mean and standard deviation in NDVI of the Spanish fir forests. Source: own elaboration. |
3.2. Cartografía de tendencia interanual
A grandes rasgos (Figura 5), encontramos una correlación lineal temporal positiva, máxima en los pinsapares de la Sierra de las Nieves y Sierra Bermeja (Los Reales, Cerro Abanto y Sierra del Real) y mínima en los pinsapares de la Sierra de Grazalema, donde se dieron incluso valores negativos. En la misma dirección, según el test de tendencia de Mann-Kendall y su significación, la tendencia positiva fue consistente en buena parte de los pinsapares de la Sierra de las Nieves y Sierra Bermeja. Por su parte, en la mayor parte del pinsapar de la Sierra de Grazalema y Sierra Blanca-Canucha, encontramos valores no significativos.
|
Figure 5: Interannual trend in NDVI of the Spanish fir forests. Source: own elaboration. |
3.3. Análisis de tendencia estacional
Los resultados del análisis de tendencia estacional refrendaron y complementaron los análisis previos. En general, confirmamos la tendencia hacia el incremento del NDVI de los pinsapares. De acuerdo a las principales zonas (Figura 6), esta tendencia observada fue mayor en Sierra Bermeja y Sierra de las Nieves, por este orden, y en ambos casos resultó consistente en todos los meses y estaciones. En Sierra Blanca-Canucha también se observó un incremento del NDVI, especialmente en primavera. Y en Sierra de Grazalema, de nuevo, no observamos una tendencia clara. También es cierto que esta zona tiene una mayor proporción de cobertura muy densa.
|
Figure 6: Seasonal trend in the main areas. Dashed line: 2002; Continuous line: 2016. Source: own elaboration. |
Caracterizamos dos periodos de reverdecimiento (Green up): una fase más intensa en otoño-invierno y otra secundaria en primavera-verano; y dos periodos de senescencia (Green down): invierno y verano. En general, los máximos del NDVI se mantuvieron en invierno, entre diciembre y enero, y los mínimos al final del verano. Hay que tener en cuenta que estos registros tienen en cuenta el conjunto de las masas, incluyendo presencia de otras cubiertas forestales.
En el conjunto del área de estudio, la estructura demográfica y diversidad de las poblaciones de Abies pinsapo está muy influida por la desigual presión antrópica ejercida en los distintos montes sujetos a diferentes sistemas de explotación en el pasado (Linares et al., 2011a). Por este motivo, también analizamos la tendencia estacional del NDVI en distintas formaciones de Abies pinsapo localizadas en los diferentes municipios del área de estudio. Con la ayuda de la fotointerpretación, discriminamos las poblaciones de Abies pinsapo según sus niveles de cobertura y las especies o formaciones acompañantes en las coberturas representadas.
Cuando Abies pinsapo se halla próximo a su óptimo habita en formaciones densas, a menudo monoespecíficas, pero también formando bosques mixtos junto con otras especies de los géneros Quercus y Pinus en el estrato arbóreo. En la Figura 7 se ilustra la tendencia estacional del NDVI para el periodo 2002 a 2016 en los pinsapares más densos.
|
Figure 7: Seasonal trend of NDVI in dense forests. Dashed line: 2002; Continuous line: 2016. Source: own elaboration. |
Incluso en los pinsapares más densos, observamos un aumento del NDVI consistente en todas las poblaciones, con la excepción del pinsapar de la Sierra del Pinar (Grazalema). No obstante, este pinsapar presentó los valores de NDVI más elevados. Todos los pinsapares densos de la Sierra de las Nieves, incluyendo los pinsapares monoespecíficos, los bosques mixtos y los pinsapares con matorral denso, mostraron una tendencia positiva del NDVI en todas las estaciones, con mayores tasas de incremento en las estaciones equinocciales.
Por último, analizamos algunos pinsapares puros situados en la Sierra de las Nieves, ubicados en las cotas más elevadas (Figura 8). Desde el punto de vista ecológico, se trata de poblaciones que representan auténticos bioindicadores para la especie, ya que podrían evidenciar un prediagnóstico de la dinámica futura del pinsapar de acuerdo con la evolución prevista de las temperaturas. El pinsapar de Ronda (Cañada del Cuerno, Cañada de Enmedio, Cañada de las Ánimas) es uno de los más singulares, tanto por la longevidad de sus ejemplares centenarios, como por el rango altitudinal que presenta, alcanzando cotas por encima de los 1.700 metros de altitud sobre el nivel del mar. Por encima de los 1.700 m, el pinsapar de la Cañada del Cuerno mostró una dinámica positiva del NDVI, pero un destacable adelantamiento de la fase de reverdecimiento (Green up phase), de entre dos y tres semanas. Esta tendencia también se observó en los pinsapares de la Cañada de Enmedio y la Cañada de las Ánimas, con un adelantamiento general del reverdecimiento (Green up). Finalmente, en los pinsapares más extensos, los de Yunquera, también observamos un incremento del NDVI en todas las estaciones. Nuevamente, observamos un adelanto del Green up (Green up onset y Green up phase). Esta tendencia también se observó en los pinsapares más elevados de Yunquera (Puerto del Hornillo) y Tolox (Cañada de Froncaile, Cañada de los Aserradores). En general, observamos un incremento generalizado del NDVI, más intenso en los máximos; un adelanto del inicio del reverdecimiento (Green up onset) a febrero-marzo, junto con un incremento de la pendiente de la curva del Green up phase; y un retraso de la fase senescencia estival (Green down phase), pero esta última precedida de un adelanto del inicio de la senescencia (Green down onset).
|
Figure 8: Seasonal trend of NDVI in the Spanish fir forest of Sierra de las Nieves. Source: own elaboration. |
4. DiscusiónTop
Inicialmente, esta investigación se planteó como un pre-diagnóstico para verificar cuál de las hipótesis sobre la evolución futura de los pinsapares podría ser más factible a corto y medio plazo, a la vista de los datos aportados por el sensor MODIS.
Por un lado, partíamos de un conjunto de hipótesis basadas en investigaciones con base ecológica sobre los efectos del calentamiento global, que auguran unos escenarios negativos para la supervivencia de especies relictas con escasa plasticidad ecológica como Abies pinsapo. Por otro lado, consideramos la eficacia de las distintas medidas de manejo y gestión de los pinsapares llevadas a cabo por la administración competente en materia de medio ambiente en la región. Esta segunda se planteó como hipótesis alternativa.
En principio, nuestros resultados apoyan la hipótesis alternativa y sugieren un éxito de la estrategia de conservación y regeneración de las formaciones con presencia de Abies pinsapo, por lo menos en la escala temporal considerada. Los resultados confirman que el aumento del NDVI es consistente en el tiempo y puede corresponder con un proceso de regeneración de los pinsapares situados en la Sierra de las Nieves y Sierra Bermeja. En efecto, desde enero de 2002 hasta diciembre de 2016 se ha producido un proceso de reverdecimiento en todos los tipos de superficies forestales cubiertas con Abies pinsapo, incluso sobre sustratos ultramáficos en poblaciones situadas en cotas no muy elevadas, lo que avala la viabilidad de la especie en las sierras peritotíticas donde la especie estuvo presente en el Holoceno Medio y hoy es escasa por la recurrencia de los incendios forestales (Gómez-Zotano et al., 2017; Olmedo-Cobo et al., 2017). En este mismo periodo, se han sucedido años húmedos y años secos, y, de hecho, los últimos años se corresponden con un periodo seco. En el caso de la Sierra de Grazalema, el hecho de que la tendencia no fuera significativa permite inferir que estamos ante un pinsapar muy denso, maduro y estable (Arista, 1995; Abellanas et al., 2016).
Cabría pensar que el tamiz que impone la resolución espacial de 250 m de las imágenes del sensor MODIS ha supuesto una limitación insalvable para conocer el verdadero comportamiento espectral de la especie, pues la reflectancia de otras cubiertas forestales y no forestales también está contenida en el registro del NDVI de cada píxel. Pero descubrimos que la tendencia es consistente en casi todas las sierras, durante todos los meses y estaciones fenológicas, con distintos grados de cobertura y en las diferentes formaciones en las que se manifiestan los pinsapares, entre las cuales se incluyen pinsapares monoespecíficos densos, cuyas copas casi monopolizan la reflectancia en numerosos píxeles muestreados. Por otra parte, esta resolución permite trabajar a escala de bosque monoespecífico y a escala de formaciones multi-específicas, porque las especies no son islas, conviven e interaccionan entre sí, dependen de una comunidad y pueden verse favorecidas o perjudicadas por la convivencia (Morelli & Tryjanowski, 2015).
En este trabajo, también hemos puesto de manifiesto un adelantamiento de la fase de reverdecimiento (Green up phase) y un retraso de la fase de senescencia (Green down phase), hecho que podría estar relacionado con un incremento progresivo de la temperatura y la aridez, y que ya ha sido constatado a nivel global por otros autores (Badeck et al., 2004; Wang et al., 2015; Ma et al., 2016). Esta tendencia se corresponde con la Clase 1 señalada por Eastman et al. (2009), que incluye -entre otros- bosques mediterráneos y bosques templados de coníferas, donde se ha detectado un aumento del NDVI promedio (Amplitud 0) desde 1982, un adelanto de la fase de reverdecimiento y un retraso de la fase de senescencia por extensión de la estación de crecimiento. Este fenómeno ha sido detectado por los citados autores en áreas boscosas, y el adelanto del Green up es interpretado como una respuesta al incremento de la temperatura (Richardson et al., 2013). Las estaciones de crecimiento más largas y más cálidas aumentan el secuestro de carbono (White et al., 2009), pero también aumentan la evapotranspiración y el estrés por la sequía (Zhang et al., 2009) y pueden aumentar la incidencia de incendios forestales (Westerling et al., 2006). A nivel especifico, el adelanto de la estación de crecimiento del pinsapo ha sido detectado por otros autores trabajando con técnicas de fenomorfología (Pérez-Latorre & Cabezudo, 2012). Se ha observado que el pinsapo es especialmente sensible a la falta de agua durante la primavera (Linares et al., 2011b). No resulta tan clara la explicación del retraso del periodo de senescencia (Green down phase), pero es posible que se deba al incremento de productividad neta de las formaciones boscosas, que produce un incremento de los valores máximos y una elongación de la curva del NDVI. En cualquier caso, este retraso del periodo de senescencia (Green down phase) viene acompañado por un sintomático adelanto del inicio del periodo de senescencia (Green down onset), hecho que está relacionado con el incremento de la temperatura. El pinsapo es una especie sensible al estrés hídrico que cierra los estomas rápidamente cuando existe poca disponibilidad de agua en el suelo (Sánchez-Salguero et al., 2015). En cualquier caso, hay que considerar las limitaciones que impone la escala espacial para determinar el fenómeno con exactitud, pues sería necesario contrastar los resultados estudiando la variabilidad a nivel de individuos (Laskurain et al., 2018).
Desde el punto de vista de la ecología funcional, el pinsapar «parece un bosque templado de coníferas pero que está sujeto a constricciones climáticas de tipo mediterráneo» (Linares & Carreira, 2006). Nuestros resultados corroboran esta afirmación que sugiere que la singularidad del pinsapar deriva de su carácter mixto, entre bosque templado boreal de coníferas y bosque mediterráneo. El pinsapo es una especie vulnerable a la sequía (Linares et al., 2010), pero el manejo adecuado de las poblaciones puede contribuir a mitigar su efecto negativo sobre la especie (Lechuga et al., 2017). En este sentido, se ha puesto de manifiesto que la alta densidad de pies produce un efecto negativo en el estado fitosanitario de las masas (Navarro-Cerrillo & Calzado Martínez, 2004; Linares et al., 2010).
Nuestros resultados parecen poner de manifiesto un efecto positivo de las estrategias de gestión activa y adaptativa de las poblaciones puestas en práctica hasta el momento, a pesar de que también se están notando los efectos del calentamiento global en los ciclos fenológicos. Dado que se estima que los peores efectos del calentamiento global están por llegar (Ripple et al., 2017), y los últimos modelos basados en los datos aportados por el 5º Informe del IPCC (IPCC, 2013) proyectan una reducción significativa del hábitat del pinsapo en la segunda mitad de la presente centuria (Gutiérrez-Hernández, 2018), es muy posible que todavía no se hayan expresado los efectos más negativos del calentamiento global. Cuanto menos, puede admitirse que el manejo de las poblaciones ha conseguido demorar los efectos más negativos. En este sentido, trabajos basados en el estudio de anillos de crecimiento de distintas especies del género Abies proyectan, siguiendo también los escenarios de emisiones del 5º Informe del IPCC, una reducción del crecimiento y un acortamiento de la estación óptima (Camarero & Gutiérrez, 2017; Sánchez-Salguero et al., 2017).
Finalmente, en lo que se refiere a los factores que pueden limitar los resultados de esta investigación, hay que citar nuevamente la resolución espacial empleada y las limitaciones del NDVI. La resolución de 250 metros puede ser suficiente o no, dependiendo de las características de la población estudiada. En los casos en que pudiera resultar insuficiente lo sería por la influencia de la señal espectral de otras cubiertas forestales o no forestales. En este sentido, no puede descartarse que ciertos procesos de reverdecimiento puedan estar camuflando el avance de unas especies (pino, encina) con respecto al pinsapo, aunque en el presente estudio hemos discriminado la mayor parte de estas interferencias, con ayuda de la fotointerpretación. Por otra parte, sería deseable comparar los resultados obtenidos con diferentes índices de vegetación, ya que algunos son más o menos sensibles a la señal del suelo o más adecuados cuando predominan niveles de cobertura forestal muy elevados (Jones & Vaughan, 2010).
5. ConclusionesTop
Los resultados obtenidos confirman la hipótesis de que la gestión adaptativa de los pinsapares está produciendo un incremento del NDVI de las superficies cubiertas con pinsapos, por lo menos a la resolución espacial (250 m) y escala temporal (2002-2016) considerada en este trabajo. Se detecta un proceso de reverdecimiento generalizado que sugiere una regeneración de las principales masas con presencia del pinsapo. La tendencia fue significativa en casi todas formaciones y grados de cobertura. Por otro lado, revelamos un adelantamiento general de la fase de reverdecimiento (Green up phase) y un retraso de la fase de senescencia (Green down phase), aunque éste último va acompañado de un adelanto del inicio de la senescencia (Green down onset). Estas observaciones podrían estar relacionadas con el calentamiento global, y también han sido señaladas por otros autores, si bien habría que considerar la existencia de cambios en la estructura poblacional de las masas.
Futuros trabajos sobre monitorización de la actividad vegetal y fenología de los pinsapares a través de imágenes de satélite y análisis estadístico de series temporales, deberían contemplar los siguientes aspectos: una extensión de la serie temporal para abordar los cambios fenológicos de la especie desde una perspectiva más amplia; un incremento de la frecuencia temporal de captura de los datos, para determinar con mayor precisión la datación de los eventos fenológicos en el tiempo; una comparación entre distintos índices de vegetación; el análisis de series temporales superpuestas sobre variables independientes tales como temperaturas, precipitaciones, etc., para avanzar en el estudio de relaciones causales; la inclusión de imágenes de mayor resolución espacial, para obtener una visión más precisa del comportamiento espectral de la especie a escala de rodal; y, finalmente, la realización de muestreos a pie de campo para correlacionar los datos reales tomados en campo con los datos obtenidos a partir de las imágenes de satélite.