El pinsapo (Abies pinsapo. Boiss) es un endemismo Ibérico (Bético) cuyas poblaciones se localizan en la Serranía de Ronda1. 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.

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).

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.