Revista Cubana de Ciencias Forestales. 2021; mayo-agosto 9(2): 158-178

Proyección del hábitat potencial de Magnolia cubensis urb. subsp. cubensis en el oriente de Cuba

Proyection of potencial habitats of Magnolia cubensis urb. subsp. cubensis in the eastern Cuba

Projeção de habitats potenciais da Magnólia cubensis urb. subsp. cubensis no oriente de Cuba

Yenia Molina-Pelegrín^1* https://orcid.org/0000-0002-9037-3817

Gretel Geada López² https://orcid.org/0000-0002-8421-0624

Adonis Sosa-López¹ https://orcid.org/0000-0002-3492-4514

Alain Puig-Pérez¹ https://orcid.org/0000-0002-5261-2709

José Luis Rodríguez-Fonseca¹ https://orcid.org/0000-0003-1395-1588

Adonis Ramón-Puebla³ https://orcid.org/0000-0002-2515-2508

¹Instituto de Investigaciones Agroforestales. UCTB Estación Experimental Agroforestal Guisa. Cuba.
²Universidad de Pinar del Río "Hermanos Saíz Montes de Oca". Pinar del Río, Cuba.
³Ministerio de Ciencia, Tecnología y Medio Ambiente. Delegación Provincial del CITMA en Granma. Órgano de Montaña Sierra Maestra. Cuba.

^*Autor para la correspondencia: yenia@guisa.inaf.co.cu

Recibido: 25/03/2021
Aprobado: 29/03/2021

RESUMEN

Los modelos de distribución de especies son una herramienta importante en estudios de conservación de especies amenazadas. Magnolia cubensis subsp. cubensis es una especie endémica categorizada como vulnerable de bosques pluvisilvas de montaña del oriente de Cuba. El presente estudio tuvo como objetivo determinar la distribución potencial actual de M. cubensis subsp cubensis en el oriente de Cuba, como base para estudios de impacto de cambio climático e identificación de sitios para la conservación. La distribución real de la subespecie mostró un patrón de distribución fragmentado para el macizo montañoso Sierra Maestra que se ubica entre los 702 a 1 396 m s.n.m y ocupa un área de 0,02 km². La modelación de los hábitats potencial muestra un comportamiento de fragmentado con un área de idoneidad de 12 428,62 ha. Las variables que más influyeron en la distribución potencial de la especie fueron las relacionadas el régimen térmico, especialmente temperatura mínima del período más frío y la altitud.

Palabras clave: Magnolia cubensis subsp. cubensis; Modelación; Distribución espacial.

ABSTRACT

Species distribution models are an important tool in conservation studies of threatened species. Magnolia cubensis subsp. cubensis is an endemic species categorized as vulnerable in mountain rainforests of eastern Cuba. The present study aimed to determine the current potential distribution of M. cubensis subsp. cubensis in eastern Cuba, as a basis for climate change impact studies and identification of sites for conservation. The actual distribution of the subspecies showed a fragmented distribution pattern for the Sierra Maestra mountain range, which is located between 702 and 1396 m a.s.l. and occupies an area of 0.02 km². The modeling of potential habitats shows a fragmented behavior with a suitability area of 12,428.62 ha. The variables that most influenced the potential distribution of the species were those related to the thermal regime, especially the minimum temperature of the coldest period and altitude.

Keywords: Magnolia cubensis subsp. cubensis; Modeling; Spatial distribution.

RESUMO

Os modelos de distribuição das espécies são um instrumento importante nos estudos de conservação de espécies ameaçadas. Magnólia cubensis subsp. cubensis é uma espécie endémica classificada como vulnerável nas florestas tropicais de montanha do oriente de Cuba. O presente estudo visava determinar a atual distribuição potencial de M. cubensis subsp. cubensis no oriente de Cuba, como base para estudos de impacto das alterações climáticas e identificação de sítios para conservação. A distribuição real da subespécie mostrou um padrão de distribuição fragmentado para a Serra Maestra, que varia de 702 a 1 396 m a.s.l. e ocupa uma área de 0,02 km². A modelação dos habitats potenciais mostra um comportamento fragmentado com uma área de aptidão de 12 428,62 ha. As variáveis que mais influenciaram a distribuição potencial da espécie foram as relacionadas com o regime térmico, especialmente a temperatura mínima do período mais frio e a altitude.

Palavras chave: Magnólia cubensis subsp. cubensis; Modelação; Distribuição espacial.

INTRODUCCIÓN

En la actualidad, uno de los mayores retos para la conservación y el manejo sostenible de la biodiversidad es mantener poblaciones viables en sus entornos naturales (Aitken et al., 2008). La conservación de la biodiversidad requiere igualmente un conocimiento razonable sobre la distribución de las especies y sus poblaciones (Margules y Sarkar 2007; Rehfeldt et al., 2015). Sin embargo, la falta de información detallada sobre la distribución geográfica de estas constituye una limitante ante tal demanda (Pinkard et al., 2015).

La información disponible para realizar iniciativas de conservación en sitios degradados o en aquellos compuestos por especies raras o escasamente estudiadas, muchas veces puede estar constituida por antecedentes incompletos (Rehfeldt et al., 2015). Esta dificultad para obtener información del sitio puede ser particularmente relevante cuando se intenta trabajar con especies o ecosistemas circunscritos a pequeños fragmentos inmersos dentro de una matriz alterada por el ser humano. La falta de información limita la capacidad para generar programas a gran escala, y como consecuencia las iniciativas propuestas, en muchos casos, tienen un éxito limitado. Por ende, un tema relevante en esta materia, es poder predecir cuáles especies o grupos de especies podrían obtener mejores resultados para ser utilizadas en proyectos de conservación para sitios específicos (Contreras-Medina et al., 2010; Morales 2012).

Para sobrepasar estas barreras de información, los modelos de distribución de especies son una valiosa herramienta para determinar la distribución de especies poco estudiadas y pueden ser de gran ayuda en la generación de información biológica base no disponible. Los modelos de distribución de especies son una herramienta empírica, robusta, repetible y fácil de utilizar, que pueden ayudar a identificar zonas potenciales a conservar y/o restaurar (Phillips y Dunik, 2008 y Arribas et al., 2012). El valor de estos estudios en la conservación está basado en que 1) indica potenciales sitios donde intervenir en el caso de necesitar hacer exploraciones directas, 2) entrega lineamientos para determinar el comportamiento de las especies frente a alteraciones exógenas, principalmente antrópicas y 3) entrega información sobre los sitios a conservar (Hernández et al., 2008 y Arribas et al., 2012).

Los algoritmos de modelado de distribución de especies toman las ubicaciones de recolección de una especie e identifican los valores de un conjunto de variables ambientales explicativas que inciden en esos sitios (Eliht et al., 2011) y devuelven las probabilidades proyectadas de presencia para esa especie a través de un área de estudio predefinida, dependiente de las variables de conocimiento. Siempre que la distribución geográfica de la especie esté bien representada por los sitios muestreados, los modelos deben identificar con precisión el nicho ecológico para la especie de interés.

Magnolia cubensis Urb. subsp. cubensis (marañón de Sierra Alta), es un taxón endémico del bosque pluvisilva de montaña en el oriente de Cuba, se desarrolla sobre los 800 m s.n.m de la Sierra Maestra y Sierra de Gran Piedra (Testé et al., 2019). Dentro del género Magnolia es uno de los taxones con más amplia distribución (Palmarola et al., 2015; 2016), pero con poblaciones fragmentadas dentro del macizo montañoso, es considerada dentro de las especies con madera preciosas cubanas (Molina-Pelegrín et al., 2014). En los últimos años, sus poblaciones naturales se han visto seriamente afectadas por el aprovechamiento maderero, otras actividades antrópicas y los eventos meteorológicos (Palmarola et al. 2015) por lo que la especie está categorizada como vulnerable (González-Torres et al., 2016). En la actualidad, se han iniciado una serie de actividades encaminadas a fortalecer su conservación y utilización sostenible, así como estudios básicos de su estructura poblacional en la Sierra de Gran Piedra (Testé et al., 2019).

Esta subespecie se encuentra sujeta a niveles de riesgo por experimentar reducciones potenciales en sus poblaciones, a causa de la falta de espacios naturales y escaso manejo. De ahí, el importante papel que juegan las áreas protegidas ya que gran parte de sus poblaciones se encuentran en estas, las acciones a realizar abarcan desde el conocimiento de los recursos del área y el manejo de especies amenazadas, hasta la educación ambiental e integración comunitaria a las tareas de conservación.

Todo lo anterior indica la importancia de realizar estudios de modelación de la distribución espacial de M. cubensis subsp. cubensis que contribuyan a la conservación de este taxón en la región oriental de Cuba. La presente investigación tiene por objetivo determinar la distribución potencial actual de M. cubensis subsp cubensis en el oriente de Cuba, como base para estudios de impacto de cambio climático e identificación de sitios para la conservación.

MATERIALES Y MÉTODOS

Caracterización del área de estudio

La subespecie tiene su distribución natural en la Sierra Maestra comprende cinco áreas protegidas: Parque Nacional La Bayamesa (PNLB) y Parque Nacional Turquino (PNT), áreas protegidas de significación nacional, Paisaje Natural Protegido Gran Piedra, así como, las Reservas Ecológicas Pico Caracas y El Gigante. Conforman el bloque montañoso por encima de los 1 200 m s.n.m. con mayor área en Cuba.

El clima del territorio es poco variado, aunque existen diferencias entre la parte norte y el borde meridional. La precipitación es elevada, pues oscila entre 1 800 y 2 300 mm anuales. La evaporación media anual es baja. La temperatura media anual varía entre 16°C y 20°C. Las temperaturas mínimas absolutas varían entre 12°C y 16°C, la media de las temperaturas máximas absolutas es de 26°C a 30°C, mientras que la media de las mínimas absolutas es de 4°C a 8°C (Montenegro, 1991).

En la Sierra Maestra, se encuentra una gran diversidad de formaciones vegetales: bosque nublado, matorral nublado, pluvisilva de montaña, pinares naturales, bosque siempreverde mesófilo, herbazal de galería y vegetación antrópica. Estas dependen fundamentalmente de las diferencias ecológicas, las que están condicionadas principalmente por la altitud y en ocasiones el edátopo (Maceira et al., 2005).

Magnolia cubensis subsp. cubensis, habita en la formación vegetal de pluvisilva de montaña, que predomina en sitios entre los 800 y los 1 400 m s.n.m. En las áreas superiores transiciona con el bosque nublado, y en las inferiores con el bosque siempreverde mesófilo. El relieve es abrupto, frecuentemente con pendientes de más de 30 grados sobre suelos ferralíticos rojos lixiviados muy pobres y ácidos (Testé et al., 2019).

Distribución geográfica

Distribución real

El área de distribución real de M. cubensis subsp. cubensis se determinó a partir de las coordenadas de cada individuo registrados con GPS Gramin. Los individuos adultos censados cuentan con una edad que oscila entre los 12 y los 60 años de edad estimada de acuerdo al fenotipo. El período de recopilación de los puntos de presencia comprende los años entre 2008 y 2017. Con estos datos se generó el mapa de distribución real de la subespecie mediante el empleo del programa ArcGis versión 10.3. Se calculó el área de ocupación real de la población, multiplicando la cantidad de celdas con presencia de al menos un individuo por el área de la celda (400 m²).

Transformación y conversión de variables bioclimáticas y topográficas

Para el estudio se emplearon 19 variables bioclimáticas: (bio01) temperatura media anual, (bio02) rango medio diurno de temperatura, (bio03) isotermalidad, (bio04) estacionalidad de la temperatura, (bio05) temperatura máxima del período más cálido, (bio6) temperatura mínima del período más frío, (bio07) variación anual de temperatura (bio05-bio06), (bio08) temperatura media del trimestre más lluvioso, (bio09) temperatura media del trimestre más seco, (bio10) temperatura media del trimestre más cálido, (bio11) temperatura media del trimestre más frío, (bio12) precipitación anual, (bio13) precipitación del período más lluvioso, (bio14) precipitación del período más seco, (bio15) estacionalidad de las precipitaciones, (bio16) precipitación en el trimestre más húmedo, (bio17) precipitación del trimestre más seco, (bio18) precipitación del trimestre más cálido y (bio19) precipitación del trimestre más frío, obtenidas de la base de datos global de superficies climáticas WorldClim, (http://www.worldclim.org). Además, se utilizó la altitud como variable, a partir del modelo digital de elevación disponible en la colección WordClim basada en los registros del radar SRTM (Shuttle Radar Topography Mission).

Primeramente, los archivos raster de las variables bioclimáticas fueron recortados para tomar el área de estudio. Para esto se utilizó una máscara proveniente de un shapefile de la región oriental de Cuba, el shapefile tuvo que reproyectarse al sistema de coordenadas WGS 84 UTM Zone 18N, una vez reproyectado, constituyó el soporte para la transformación del resto de las variables bioclimáticas y altitud.

La configuración de los raster ambientales se llevó a cabo con la aplicación de la herramienta Entornos Ambientales, que permite que todas las capas ambientales tengan la misma extensión, tamaño de celda y sistema de coordenadas; estos son requisitos indispensables para que los archivos sean procesados por el programa MaxEnt. Una vez configurados los raster fueron convertidos a formato ASCII para poder ser ejecutados por MaxEnt, el proceso de conversión se realizó con el empleo del programa ArcGis 10.3.

Modelo de la distribución

Modelo de la distribución potencial actual

Se modeló la distribución de M. cubensis subsp. cubensis utilizando el programa MaxEnt v. 3.3.3 (Phillips et al., 2006). El resultado de la modelación con MaxEnt, revela la probabilidad relativa de distribución de una especie en todas las cuadrículas o celdas en el espacio geográfico definido, en el cual un valor de alta probabilidad asociado a una cuadrícula particular indica la probabilidad de la misma, de tener condiciones ambientales favorables para la especie modelada (Elith y Leathwick, 2009).

Para la modelación se partió de las coordenadas de los 589 puntos de presencia, puntos generados durante el censo a los individuos de M. cubensis subsp. cubensis, que comprenden todo el área de distribución natural. Las variables utilizadas en esta primera modelación fueron las 19 variables ambientales y la altitud como variable topográfica.

La robustez de predicción del modelo se determinó a través de la técnica de evaluación de la curva operacional (ROC) a través del área bajo la curva ROC (AUC) generando un análisis de sensibilidad vs. Especificidad (Phillips et al., 2006), para esto se utilizó un 75 % de las accesiones de entrenamiento y un 25 % de prueba en la opción random test percentage, la cual le indica al programa que debe escoger de manera aleatoria y apartar el 25 % de los registros de muestras para usarlos en la prueba para la comprobación de los modelos resultantes. El AUC se puede entonces interpretar como la mayor probabilidad de que un punto de presencia seleccionado aleatoriamente esté ubicado en una celda del raster con un alto valor de probabilidad para la presencia de la especie en comparación al de un punto generado aleatoriamente (Elith et al., 2011 y Ávila-Coria et al., 2014). La mayor capacidad de predecir de un modelo generado por MaxEnt se alcanza cuando el AUC tiene un valor de 1 (Araújo y Guisan, 2006).

En la configuración, se especificó un máximo de 500 interacciones con un límite de convergencia prefijado en 0,00001 que garantiza la convergencia del algoritmo, es decir que mantiene una estabilidad en sus análisis (Elith et al., 2011). La definición de estos parámetros del límite de convergencia son usualmente los empleados (Ávila et al., 2014). Con el fin de evitar datos de presencia con registros duplicados se activó borrar duplicados, lo que permitió reducir los sesgos del muestreo.

Para validar y entrenar al mismo tiempo sin que se produzca sobreajuste en la red neural se utilizó la opción crossvalidate. Se estableció como formato de salida la función logística, por ser la más sencilla de conceptualizar, ya que proporciona un estimado entre cero y uno de probabilidad de presencia (Phillips y Dudik, 2010). Esta representa la probabilidad estimada de presencia de la especie, dadas las restricciones impuestas por las variables predictivas. Con este formato, las cuadrículas con valores bajos (cercanos a cero) se consideran inapropiadas para la especie dado su nicho ecológico (Phillips y Dudik, 2010; Elith et al., 2011). El modelo final de la distribución potencial se tomó a partir del mapa en formato ASCII generado por el programa, este mapa fue convertido a formato raster mediante el empleo del programa ArcGis 10.3. Se determinó a su vez, el área que abarca la distribución potencial actual, mediante la creación de polígonos y determinando el área que ocupa dentro del territorio, para ello se utilizaron las herramientas de ArcMap.

Para evaluar la contribución de las variables bioclimáticas dentro del modelo se utilizó el método de remuestreo Jackknife. Esta prueba analiza la contribución de cada variable de manera aislada y adicionalmente genera otro con todas las variables, en donde indica la ganancia de cada una de estas. Para realizar el análisis de la influencia de cada variable climática en la presencia de la especie fueron generadas las curvas de respuesta, ya que estas expresan cómo afectan cada una de estas la predicción de MaxEnt.

Evaluación de zonas de conservación de la subespecie

La determinación de las zonas potenciales y críticas para la conservación de la subespecie, se propusieron mediante la unión de los mapas de distribución geográfica generados por el programa MaxEnt (puntos de presencia y distribución potencial) con la zonificación del Sistema Nacional de Áreas Protegidas (SNAP) de acuerdo al CNAP (2014). Para ello, se utilizó la versión digital del mapa del Centro Nacional de Áreas Protegidas (CNAP) para el período 2014-2020. Se cuantificó el área de las zonas de alta idoneidad para la presencia de la subespecie en cada una de las áreas protegidas objeto de estudio con el programa ArGis 10.3.

RESULTADOS Y DISCUSIÓN

Distribución geográfica

Distribución real

Magnolia cubensis subsp. cubensis ocupa un área de 0,02 km², de ellos 0,0184 km² corresponden a la subpoblación localizada en la provincia Granma y 0,0016 km² a la subpoblación de la provincia Santiago de Cuba (Figura 1). Desde el punto de vista biológico, la distribución pudiera estar relacionada con condiciones ecológicas favorables para su desarrollo solo en las montañas de la Sierra Maestra (Testé et al., 2019).

La disposición actual de su distribución en áreas pequeñas y tan separadas, desde punto de vista geográfico, puede deberse a la fragmentación del hábitat como consecuencia del aprovechamiento del taxón, la poca disponibilidad de individuos reproductores y problemas relacionados con el banco de semillas; sin embargo, estudios que se concentren en estos aspectos son muy escasos (Molina-Pelegrín, 2014). Hasta el momento solo se dispone de un estudio del impacto de la fragmentación del hábitat en la estructura poblacional y genética de M. cubensis subsp. acunae a través de análisis espacial y molecular, que permitió concluir que esta es la causa principal de la reducción de su distribución natural y que, además, incrementa el riesgo ante los escenarios de cambio climático (Hernández et al., 2020). Por lo que estudios moleculares y de autoecología en la subespecie cubensis son necesarios.

Este patrón de distribución en M. cubensis subsp. cubensis y la estructura poblacional de clases diamétricas y de altura evaluada por Molina-Pelegrín et al., (2014) y por Testé et al., (2019) para Gran Piedra, revelaron el envejecimiento de estas poblaciones, lo cual podría explicar la reducida área de ocupación actual.

La distribución geográfica de los individuos de la subespecie comprende desde los 19°58' 10,2'' N y 77°00'29,6'' W hasta los 20°0'46,0'' N y 75°37'50,1'' W. El rango altitudinal se encuentra entre los 702 y los 1 396 m s.n.m, aunque mayor parte de los individuos se encuentran en el rango de 1 000 hasta 1 200 m s.n.m.

Figura 1. - Distribución del área real ocupada, a partir de la georreferenciación de los 587 individuos, de Magnolia cubensis subsp. cubensis

Modelo de distribución potencial actual

Existe un buen ajuste para predecir la presencia de la especie, ya que el valor del área bajo la curva ROC fue de 0,996 (Figura 2), valores por encima de 0,9 son considerados como indicadores de buena robustez de predicción del modelo (Phillips et al., 2008 y Elith et al., 2011). Este resultado fue similar en estudios de predicción de hábitat potencial de coníferas (Contreras-Medina et al., 2010; García-Aranda et al., 2012, Ávila-Coria et al., 2014 y Miranda-Sierra 2017) y otras especies de interés forestal (Smith et al., 2012). Sin embargo, en especies de animales se reportan valores inferiores de AUC en la modelación con variables climáticas, al parecer por la influencia de otras variables ambientales en la definición de su nicho climático (Azor-Hernández y Barro-Cañamero 2014; Cobos-Cobos, 2016).

Figura 2. - Curva de respuesta operativa (ROC) del modelo de distribución potencial actual de Magnolia cubensis subsp. cubensis en el oriente de Cuba
Nota: en rojo la curva ROC promedio para los datos de entrenamiento. La línea negra representa la predicción aleatoria (AUC=0,5)

El mapa de probabilidad de presencia potencial, según el umbral elegido, indicó ausencias de la subespecie por debajo de 0,15. El modelo destacó dos núcleos geográficos principales: Sierra Maestra (provincia Granma) y Gran Piedra, también perteneciente al mismo macizo montañoso, pero este último ubicado en la provincia de Santiago de Cuba. Aunque se muestran seis pequeños núcleos poblacionales en el macizo montañoso Sagua Nipe Baracoa, con posibilidades de presencia del taxón. De forma general, se manifiesta un área con probabilidad media (0,38-0,69) de presencia y se reflejó un patrón disyunto en la distribución de la subespecie (Figura 3). La región de Sierra Maestra, en la parte oeste (provincia Granma), comprende la mayor área de ocupación potencial con 54 081,48 ha (92,80 %), mientras Gran Piedra (este) comprende 4 168, 36 ha (7,15 %), para un total de 58 249, 83 ha.

Figura 3. - Mapa de distribución potencial actual de Magnolia cubensis subsp. cubensis en el oriente de Cuba

El análisis de la importancia de cada variable mediante la prueba de jackknife (Figura 4), muestra que las variables que más influyen en la presencia de la especie, son las relacionadas con las temperaturas (bio01, bio05, bio06, bio08, bio09 y bio10) y la altitud. Estos resultados permiten identificar las exigencias climáticas de las especies y sus rangos que determinan la probabilidad de presencia. Por lo tanto, este constituye un primer acercamiento a la evaluación de la influencia de las variables climáticas en la distribución de la especie. Asimismo, esta ha sido una alternativa empleada por varios autores para algunas especies endémicas o en peligro de extinción y de esta forma facilitar la toma de decisiones en la conservación como: Cedrela odorata (Gómez-Díaz et al., 2007); Taxus globosa (García-Aranda et al., 2012); Pinus strobus (Joyce y Rehfeldt, 2013); Pinus herrerae (Ávila-Coria et al., 2014); Magnolia schiedeana (Vásquez-Morales et al., 2014) y otras (ver revisión Gray y Hamann 2012).

En este caso, se explica que las variables relacionadas con el régimen térmico son las que afectan, en primer orden, la distribución potencial del taxón. Por tanto, si hay cambios bruscos y prolongados en los rangos para una o varias variables climáticas entonces se vería afectada las condiciones climáticas favorables para la presencia del taxón. Al respecto, McLean (2015) proyecta un clima más seco e incrementos de las temperaturas en el área del Caribe, en particular al oriente de Cuba; esto podría suponer afectaciones en la especie relacionados con sus procesos fisiológicos, la producción de biomasa e inclusive la floración y producción de semilla (Álvarez-Brito, 2017 y Sáenz-Romero, 2014). Esto hará que estas sean sensibles como resultado de alteraciones en su nicho climático.

Estudios de la influencia de los cambios climáticos en procesos fisiológicos y reproductivos ha sido estudiado en varias coníferas y especies latifolias con amplia distribución (Aitken et al., 2008), pero en especies tropicales e insulares son muy escasos (Álvarez-Brito et al., 2014 y Álvarez-Brito, 2018) y estos apoyado además en estudios genéticos de la especie y sus poblaciones son la base para promover la resiliencia y el manejo adaptativo. En Cuba, dichos estudios, aunque incipientes, se dispone de información de alteraciones de las fenofases vegetativas y reproductivas en tres especies forestales (Hechavarría-Kindelán et al., 2008; Hechavarría-Kindelán 2012 y Álvarez-Brito y Mercadet-Portillo, 2014). Por lo que, investigaciones sobre la fenología y biología reproductivas en especies de árboles y sus cambios por efecto de las variaciones en el clima, junto a estos, pero a partir de la simulación de cambios en las variables climáticas en condiciones de viveros, áreas de conservación in situ son cruciales en la comprensión del impacto de los cambios en la estructura demográfica de las especies.

Figura 4. - Contribución de las variables ambientales al modelo de distribución potencial actual de Magnolia cubensis subsp. cubensis en el oriente de Cuba. Las barras azules representan la ganancia promedio del modelo usando solo una variable y las verdes la ganancia cuando cada variable fue excluida

En el análisis de la respuesta de las variables que más contribuyen al modelo reveló que, en especial, la curva de respuesta de la variable temperatura mínima del período más frío (°C) tuvo una forma exponencial inversa, a menores de temperatura correspondió a aumentos en la probabilidad de presencia de la especie, es aumentos de la idoneidad para la presencia de la subespecie (Figura 5A). Para el caso de la variable altitud, se muestra que la probabilidad de presencia del taxón aumentó desde los 700 hasta alrededor de los 1 400 m s.n.m, a partir del cual comenzó a perder relevancia (Figura 5B). La implicación de la altitud como una de las variables más explicativas en los modelos concuerda con lo reportado por González-Torres et al., (2013) y Palmarola et al., (2018) para la subespecie M. cubensis subsp. acunae en el macizo de Guamuhaya y con los estudios de M. schiedeana en México (Vásquez-Morales et al., 2014).

Figura 5. - Curvas de respuestas promedio para (A) temperatura mínima del período frío (BIO 6) y (B) elevación sobre el nivel del mar (Altitud) en el modelo de distribución potencial actual de Magnolia cubensis subsp. cubensis en el oriente de Cuba

Zonas potenciales y críticas para la conservación de la subespecie

El 78,58 % de la distribución potencial actual de M. cubensis subsp. cubensis se encuentra bajo objeto de conservación de áreas pertenecientes al Sistema Nacional de Áreas Protegidas de Cuba (SNAP), con una superficie de 45 773, 38 ha. A continuación, se muestran los mapas de las zonas potenciales para la conservación in situ del taxón (Figura 6, 7, 8, 9 y 10) y que corresponden a cada una de las áreas protegidas.

Figura 6. - Mapa de zonas potenciales para la conservación in situ de Magnolia cubensis subsp. cubensis en el Parque Nacional La Bayamesa, basado en la distribución potencial

Figura 7. - Mapa de zonas potenciales para la conservación in situ de Magnolia cubensis subsp. cubensis en el Parque Nacional Turquino, basado en la distribución potencial

Figura 8. - Mapa de zonas potenciales para la conservación in situ de Magnolia cubensis subsp. cubensis en la Reserva Ecológica El Gigante, basado en la distribución potencial

Figura 9. - Mapa de zonas potenciales para la conservación in situ de Magnolia cubensis subsp. cubensis en el Paisaje Natural Protegido Gran Piedra, basado en la distribución potencial

Figura 10. - Mapa de zonas potenciales para la conservación in situ de Magnolia cubensis subsp. cubensis en la Reserva Ecológica Pico Caracas, basado en la distribución potencial

Se muestra una baja correspondencia de la distribución potencial de M. cubensis subsp. cubensis en el SNAP que contrasta con una alta representatividad real de la subespecie en el estudio desarrollado.

Se muestra una baja correspondencia en la distribución de taxón censada por el SNAP y la presentada en este estudio como distribución actual real, lo que demuestra que no toda la especie está únicamente dentro de las áreas protegidas, y que, por tanto, quedan aún individuos formando parte de espacios donde no existe protección. En este sentido, Hernández y Estrada (2004) plantean que la delimitación de las áreas protegidas se ha establecido en base a razones logísticas, a partir de criterios antropocéntricos o geográficos y sin tomar en cuenta los requerimientos de las especies a proteger. De igual manera, ocurre con relación al área potencial y su correspondiente distribución real, en casi todas las áreas protegidas hay hábitats potenciales donde la especie podría desarrollarse de hacerse acciones para su propagación.

Las zonas críticas de conservación de la subespecie, principalmente se encuentran en la zona de amortiguamiento de las áreas protegidas; esto refuerza la importancia de estas zonas dentro de la estrategia de conservación del taxón y la necesidad de aunar acciones más interdisciplinarias e interadministrativas de las áreas protegidas. Se determinó la existencia de 12 428,62 ha que no forman parte del objeto de conservación de la subespecie y que muestran sitios con probabilidades de presencia del taxón (Figura 11). Algunos de estos sitios están contenidos dentro de la zona de amortiguamiento del Parque Nacional Pico La Bayamesa, lugar donde además se encuentran los mayores asentamientos poblacionales rurales, factor que constituye un riesgo para la conservación in situ de la subespecie.

Figura 11. - Mapa de sitios con probabilidades para hábitat potencial Magnolia cubensis subsp. cubensis, no incluidas dentro de las áreas protegidas en el oriente de Cuba

CONCLUSIONES

Magnolia cubensis subsp. cubensis se distribuye de forma fragmentada en el macizo montañoso Sierra Maestra. Los individuos mayoritariamente se distribuyen en el rango altitudinal comprendido desde los 1 000 hasta los 1 200 m s.n.m.

El área de distribución potencial actual, en el oriente de Cuba, presenta un patrón disyunto con un área de idoneidad total en el macizo montañoso Sierra Maestra, condicionada fundamentalmente por las variables: temperatura mínima del período más frío y elevación sobre el nivel del mar.

El 78,58 % de la distribución potencial actual de Magnolia cubensis subsp. cubensis se encuentra bajo objeto de conservación por parte del Sistema Nacional de Áreas Protegidas de Cuba (SNAP).

REFERENCIAS BIBLIOGRÁFICAS

AITKEN, S.N., YEAMAN, S., HOLLIDAY, J.A., WANG, T. y CURTIS-MCLANE, S., 2008. Adaptation, migration or extirpation: climate change outcomes for tree populations. Evolutionary Applications, vol. 1, no. 1, pp. 95-111. ISSN 1752-4571. DOI 10.1111/j.1752-4571.2007.00013.x.

ÁLVAREZ BRITO, A., 2017. Impactos y adaptación al cambio climático en el sector forestal cubano: Sexta aproximación. La Habana: Instituto de Investigaciones Agroforestales. ISBN 978-959-7215-32-5. Disponible en: https://rc.upr.edu.cu/jspui/handle/DICT/2876

ÁLVAREZ BRITO, A., ORTIZ, O., CORDERO, E., HECHAVARRÍA, O., SUÁREZ, T. y ESCARRÉ, A., 2014. El sector forestal cubano y el cambio climático. Anales de la Academia de Ciencias de Cuba [en línea], vol. 4, no. 2, pp. 1-11. [Consulta: 30 marzo 2021]. ISSN 2304-0106. Disponible en: http://www.revistaccuba.cu/index.php/revacc/article/view/135.

ARAÚJO, M.B. y GUISAN, A., 2006. Five (or so) challenges for species distribution modelling. Journal of Biogeography [en línea], vol. 33, no. 10, pp. 1677-1688. [Consulta: 30 marzo 2021]. ISSN 1365-2699. DOI https://doi.org/10.1111/j.1365-2699.2006.01584.x. Disponible en: https://onlinelibrary.wiley.com/doi/abs/10.1111/j.1365-2699.2006.01584.x.

ARRIBAS, P., ABELLÁN, P., VELASCO, J., BILTON, D.T., LOBO, J.M., MILLÁN, A. y SÁNCHEZ FERNÁNDEZ, D., 2012. La vulnerabilidad de las especies frente al cambio climático, un reto urgente para la conservación de la biodiversidad. Ecosistemas [en línea], vol. 21, no. 3, pp. 79-84. [Consulta: 30 marzo 2021]. ISSN 1697-2473. DOI 10.7818/ECOS.701. Disponible en: https://www.revistaecosistemas.net/index.php/ecosistemas/article/view/701.

ÁVILA CORIA, R., VILLAVICENCIO GARCÍA, R. y RUIZ CORRAL, J.A., 2014. Distribución potencial de Pinus herrerae Martínez en el occidente del estado de Jalisco. Revista mexicana de ciencias forestales [en línea], vol. 5, no. 24, pp. 92-109. [Consulta: 15 diciembre 2020]. ISSN 2007-1132. Disponible en: http://www.scielo.org.mx/scielo.php?script=sci_abstract&pid=S2007-11322014000400009&lng=es&nrm=iso&tlng=es.

COBOS-COBOS, M.E., 2016. Posibles implicaciones del cambio climático sobre la distribución de las especies del género Peltophryne (Anura: Bufonidae) en Cuba. La Habana, 2016. Tesis de Maestría. La Habana: Universidad de La Habana. Disponible en: https://www.semanticscholar.org/paper/Posibles-implicaciones-del-cambio-clim%C3%A1tico-sobre-Cobos-Emanuel/9b268065d35af5208fe6c45f91b41aeffff17579

CONTRERAS MEDINA, R., LUNA VEGA, I. y RÍOS MUÑOZ, C., 2010. Distribución de Taxus globosa (Taxaceae) en México: Modelos ecológicos de nicho, efectos del cambio del uso de suelo y conservación. Revista Chilena de Historia Natural [en línea], vol. 83, pp. 421-433. DOI 10.4067/S0716-078X2010000300009. Disponible en: https://www.researchgate.net/publication/233785982_Distribucion_de_Taxus_globosa_Taxaceae_en_Mexico_Modelos_ecologicos_de_nicho_efectos_del_cambio_del_uso_de_suelo_y_conservacion

ELITH, J. y LEATHWICK, J.R., 2009. Species distribution models: ecological explanation and prediction across space and time. Annual Review of Ecology, Evolution and Systematics [en línea], vol. 40. [Consulta: 26 abril 2021]. Disponible en: https://findanexpert.unimelb.edu.au/scholarlywork/318062-species-distribution-models?ecological-explanation-and-prediction-across-space-and-time.

ELITH, J., PHILLIPS, S.J., HASTIE, T., DUDÍK, M., CHEE, Y.E. y YATES, C.J., 2011. A statistical explanation of MaxEnt for ecologists. Diversity and Distributions [en línea], vol. 17, no. 1, pp. 43-57. [Consulta: 26 abril 2021]. ISSN 1472-4642. Disponible en: https://onlinelibrary.wiley.com/doi/abs/10.1111/j.1472-4642.2010.00725.x.

GARCÍA ARANDA, M., AYALA, C., ESTRADA, E., PANDO-MORENO, M. y MORENO TALAMANTES, A., 2012. Distribución actual y potencial de Taxus globosa Taxaceae en México. Journal of the Botanical Research Institute of Texas [en línea], vol. 6, no. 2, pp. 587-598. Disponible en: https://www.researchgate.net/publication/246548792_DISTRIBUCION_ACTUAL_Y_POTENCIAL_DE_TAXUS_GLOBOSA_TAXACEAE_EN_MEXICO.

GÓMEZ DÍAZ, J.D., MONTERROSO RIVAS, A.I., TINOCO RUEDA, J.Á., GÓMEZ DÍAZ, J.D., MONTERROSO RIVAS, A.I. y TINOCO RUEDA, J.Á., 2007. Distribución del cedro rojo (Cedrela odorata L.) en el estado de Hidalgo, bajo condiciones actuales y escenarios de cambio climático. Madera y bosques [en línea], vol. 13, no. 2, pp. 29-49. [Consulta: 27 abril 2021]. ISSN 1405-0471. DOI 10.21829/myb.2007.1321227. Disponible en: http://www.scielo.org.mx/scielo.php?script=sci_abstract&pid=S1405-04712007000200029&lng=es&nrm=iso&tlng=es.

GONZÁLEZ TORRES, L.R., BEJERANO, A., BÉCQUER, E.R., BERAZAÍN, R., BARRIOS, D. y ECHEVERRÍA, J., 2013. Top 50 - las 50 plantas más amenazadas de Cuba (Bissea 7 NE1) [en línea]. La Habana: Jardín Botánico Nacional. Disponible en: https://www.researchgate.net/publication/281242914_Top_50_-_las_50_plantas_mas_amenazadas_de_Cuba_Bissea_7_NE1.

GONZÁLEZ TORRES, L.R., PALMAROLA BEJERANO, A., GONZÁLEZ OLIVA, L., BÉCQUER, E.R., TESTÉ, E., BARRIOS VALDÉS, D., ACOSTA RAMOS, Z., ALOMÁ MORENO, O., ÁLVAREZ MONTES DE OCA, J.C., BERAZAÍN ITURRALDE, R.C., BONET MAYEDO, W.E., CABALLERO TIHERT, L., CAPOTE LÓPEZ, R.P., CARMENATE REYES, W., CASTAÑEDA NOA, I., CASTAÑEIRA COLOMÉ, M.A., CATASÚS GUERRA, L.J., CEJAS RODRÍGUEZ, F., FAGILDE ESPINOSA, M. del C., FALCÓN HIDALGO, B., FERNÁNDEZ GRANDA, L. y FERNÁNDEZ ZEQUEIRA, M.D., 2016. Lista Roja de la Flora de Cuba 2016 [en línea]. La Habana: Jardín Botánico Nacional. [Consulta: 14 septiembre 2020]. ISBN 978-959-300-113-7. Disponible en: http://repositorio.geotech.cu/jspui/handle/1234/1054.

GRAY, L. y HAMANN, A., 2012. Tracking suitable habitat for tree populations under climate change in western North America. Climatic Change [en línea], vol. 117, pp. 1-2. DOI 10.1007/s10584-012-0548-8. Disponible en: https://www.researchgate.net/publication/257548020_Tracking_suitable_habitat_for_tree_populations_under_climate_change_in_western_North_America.

HECHAVARRÍA KINDELAN, O., 2012. Alteraciones en la época de recolección de frutos de especies forestales por aumento de la temperatura en zonas montañosas de Topes de Collantes. Revista Forestal Baracoa [en línea], vol. 31, no. 1, pp. 73-78. ISSN 0138-6441. Disponible en: http://www.actaf.co.cu/revistas/rev_forestal/Baracoa-2012-1/FAO1%202012/ALTERACIONES%20EN%20LA%20%C3%89POCA%20DE%20RECOLECCI%C3%93N.pdf

HECHAVARRÍA KINDELAN, O., ÁLVAREZ BRITO, A. y MONTALVO GUERREO, J.M., 2008. Respuesta fenológica de Juglans jamaicensis subsp. jamaicensis al aumento de la temperatura en bosque pluvial montano. Revista Forestal Baracoa [en línea], vol. 27, no. 2, pp. 81-89. Disponible en: http://www.actaf.co.cu/revistas/rev_forestal/Baracoa-2008-2/FAO%202%202008/RESPUESTA%20JUGLANS.pdf

HERNÁNDEZ, L.A. y CAÑAMERO, A.B., 2014. Modelación de la distribución potencial de mariposas endémicas cubanas (Lepidoptera: Papilionoidea). Revista Cubana de Ciencias Biológicas [en línea], vol. 3, no. 3, pp. 18-30. [Consulta: 30 marzo 2021]. ISSN 2307-695X. Disponible en: http://www.rccb.uh.cu/index.php/RCCB/article/view/75.

HERNÁNDEZ, M., BEJERANO, A., VELTJEN, E., ASSELMAN, P., TESTÉ, E., LARRIDON, I., SAMAIN, M.-S. y LUIS ROBERTO, G.-T., 2020. Population structure and genetic diversity of Magnolia cubensis subsp. acunae (Magnoliaceae): effects of habitat fragmentation and implications for conservation. Oryx [en línea], vol. 54, no. 4, pp. 1-9. DOI 10.1017/S003060531900053X. Disponible en: https://www.researchgate.net/publication/339364440_Population_structure_and_genetic_diversity_of_Magnolia_cubensis_subsp_acunae_Magnoliaceae_effects_of_habitat_fragmentation_and_implications_for_conservation

HERNANDEZ, P.A., FRANKE, I., HERZOG, S., PACHECO, V., PANIAGUA, L., HEIDI, Q., SOTO, A., SWENSON, J., TOVAR INGAR, C., VALQUI, T., VARGAS, J. y YOUNG, B., 2008. Predicting species distributions in poorly-studied landscapes. Biodiversity and Conservation [en línea], vol. 17, no. 6, pp. 1353-1366. DOI 10.1007/s10531-007-9314-z. Disponible en: https://www.researchgate.net/publication/225117738_Predicting_species_distributions_in_poorly-studied_landscapes.

JOYCE, D.G. y REHFELDT, G.E., 2013. Climatic niche, ecological genetics, and impact of climate change on eastern white pine (Pinus strobus L.): Guidelines for land managers. Forest Ecology and Management [en línea], vol. 295, pp. 173-192. [Consulta: 27 abril 2021]. ISSN 0378-1127. DOI 10.1016/j.foreco.2012.12.024. Disponible en: https://www.sciencedirect.com/science/article/pii/S0378112712007475.

MACEIRA, D., FONG, A., ALVERSON, W.S. y WACHTER, T., 2005. Cuba: Parque Nacional La Bayamesa [en línea]. Chicago: The Field Museum. ISBN 0-914868-62-4. Disponible en: https://www.worldcat.org/title/cuba-parque-nacional-la-bayamesa/oclc/67786744

MARGULES, C. y SARKAR, S., 2007. Systematic Conservation Planning [en línea]. Cambridge, UK: Cambridge University Press. [Consulta: 27 abril 2021]. ISBN 978-0-521-87875-3. Disponible en: http://www.cambridge.org/au/academic/subjects/life-sciences/ecology-and-conservation/systematic-conservation-planning?format=PB&isbn=9780521703444.

MCLEAN, N., STEPHENSON, T., TAYLOR, M. y CAMPBELL, J., 2015. Characterization of Future Caribbean Rainfall and Temperature Extremes across Rainfall Zones. Advances in Meteorology [en línea], vol. 2015, no. 9, pp. 1-18. DOI 10.1155/2015/425987. Disponible en: https://www.researchgate.net/publication/276091258_Characterization_of_Future_Caribbean_Rainfall_and_Temperature_Extremes_across_Rainfall_Zones.

MIRANDA SIERRA, C.A., 2017. Cambios en la distribución potencial de pinares del occidente de Cuba en diferentes escenarios de cambio climático [en línea]. Thesis. Pinar del Río: Universidad de Pinar del Río «Hermanos Saíz Montes de Oca». Facultad de Ciencias Forestales y Agropecuarias. Departamento Forestal. [Consulta: 27 abril 2021]. Disponible en: https://rc.upr.edu.cu/jspui/handle/DICT/2459.

MOLINA PELEGRÍN, Y., SANTOS CHACÓN, W., SOSA LÓPEZ, A., ARCIA CHÁVEZ, M., HECHAVARRÍA KINDELAN, O. y ROSALES RODRÍGUEZ, M., 2014. Estructura poblacional de Magnolia cubensis subsp. cubensis en la Reserva Ecológica El Gigante. Revista Forestal Baracoa [en línea], vol. 33, no. 1, pp. 15-23. [Consulta: 27 abril 2021]. DOI 10.5281/zenodo.4432698. Disponible en: https://zenodo.org/record/4432698.

MONTENEGRO, V., 1991. Atlas de Santiago de Cuba. [carte]. Santiago de Cuba: Academia de Ciencias de Cuba.

MORALES, N., 2012. Modelos de distribución de especies: Software MaxEnt y sus aplicaciones en Conservación. Revista Conservación Ambiental [en línea], vol. 2, no. 1, pp. 1-5. [Consulta: 27 abril 2021]. Disponible en: https://issuu.com/fundacionecomabi/docs/revista_conservaci__n_ambiental_mas.

PALMAROLA, A., GRANADO, L., TESTÉ, E., RODRÍGUEZ, M.H., ALBELO, N. y GONZÁLEZ-TORRES, L.R., 2018. Estructura poblacional y distribución de Magnolia cubensis subsp. acunae (Magnoliaceae). Revista del Jardín Botánico Nacional [en línea], vol. 39, pp. 103-111. [Consulta: 27 abril 2021]. ISSN 2410-5546. Disponible en: http://www.rjbn.uh.cu/index.php/RJBN/article/view/429.

PALMAROLA, A., ROMANOV, M.S., BOBROV, A.V.F.C. y GONZÁLEZ-TORRES, L.R., 2016. Las magnolias de Cuba: Talauma taxonomía y nomenclatura: Revista del Jardín Botánico Nacional [en línea], vol. 37, pp. 1-10. [Consulta: 27 abril 2021]. ISSN 0253-5696. Disponible en: https://arcabc.ca/islandora/object/dc:43832.

PHILLIPS, S. y DUDÍK, M., 2008. Modeling of species distributions with MAXENT: new extensions and a comprehensive evaluation. Ecography [en línea], vol. 31, no. 2, pp. 161-175. DOI 10.1111/j.0906-7590.2008.5203.x. Disponible en: https://www.researchgate.net/publication/227635331_Modeling_of_species_distributions_with_MAXENT_new_extensions_and_a_comprehensive_evaluation.

PHILLIPS, S.J., ANDERSON, R.P. y SCHAPIRE, R.E., 2006. Maximum entropy modeling of species geographic distributions. Ecological Modelling [en línea], vol. 190, no. 3, pp. 231-259. [Consulta: 27 abril 2021]. ISSN 0304-3800. DOI 10.1016/j.ecolmodel.2005.03.026. Disponible en: https://www.sciencedirect.com/science/article/pii/S030438000500267X.

PINKARD, E., BATTAGLIA, M., BRUCE, J., MATTHEWS, S., CALLISTER, A., HETHERINGTON, S., LAST, I., MATHIESON, S., MITCHELL, C., MOHAMMED, C., MUSK, R., RAVENWOOD, I., ROMBOUTS, J., STONE, C. y WARDLAW, T., 2014. A history of forestry management responses to climatic variability and their current relevance for developing climate change adaptation strategies. Forestry [en línea], vol. 88, no. 2, pp. 155-171. DOI 10.1093/forestry/cpu040. Disponible en: https://www.researchgate.net/publication/276835826_A_history_of_forestry_management_responses_to_climatic_variability_and_their_current_relevance_for_developing_climate_change_adaptation_strategies

REHFELDT, G., WORRALL, J., MARCHETTI, S. y CROOKSTON, N., 2015. Adapting forest management to climate change using bioclimate models with topographic drivers. Forestry [en línea], vol. 88, no. 5. DOI 10.1093/forestry/cpv019. Disponible en: https://www.researchgate.net/publication/279222422_Adapting_forest_management_to_climate_change_using_bioclimate_models_with_topographic_drivers.

SÁENZ ROMERO, C., 2014. Guía técnica para la planeación de la reforestación adaptada al cambio climático [en línea]. México: Comisión Nacional Forestal. Disponible en: http://www.conafor.gob.mx:8080/documentos/docs/19/6688Gu%C3%ADa%20T%C3%A9cnica%20para%20la%20Planeaci%C3%B3n%20de%20la%20Reforestaci%C3%B3n.pdf

SMITH, A., PAGE, B., DUFFY, K. y SLOTOW, R., 2012. Using Maximum Entropy modeling to predict the potential distributions of large trees for conservation planning. Ecosphere [en línea], vol. 3, no. 6, pp. 56. DOI 10.1890/ES12-00053.1. Disponible en: https://www.researchgate.net/publication/271200752_Using_Maximum_Entropy_modeling_to_predict_the_potential_distributions_of_large_trees_for_conservation_planning.

TESTÉ, E., GORDILLO, M., PALMAROLA, A., HERNÁNDEZ, M. y GONZÁLEZ TORRES, L.R., 2019. Estructura poblacional de Magnolia cubensis subsp. cubensis (Magnoliaceae) en el Paisaje Natural Protegido Gran Piedra. Revista del Jardín Botánico Nacional [en línea], vol. 40, pp. 19-21. [Consulta: 27 abril 2021]. ISSN 2410-5546. Disponible en: http://www.rjbn.uh.cu/index.php/RJBN/article/view/439.

Conflicto de intereses:
Los autores declaran no tener conflictos de intereses.

Contribución de los autores:
Yenia Molina-Pelegrín: Concepción de la idea, búsqueda y revisión de literatura, confección de instrumentos, aplicación de instrumentos, confección de tablas, gráficos e imágenes, redacción del original (primera versión), traducción de términos o información obtenida, revisión de la aplicación de la norma bibliográfica aplicada.

Gretel Geada López: Búsqueda y revisión de literatura, asesoramiento general por la temática abordada, revisión y versión final del artículo, corrección del artículo, traducción de términos o información obtenida.

Adonis Sosa-López: Concepción de la idea, confección de instrumentos, aplicación de instrumentos, recopilación de la información resultado de los instrumentos aplicados, análisis estadístico, confección de tablas, gráficos e imágenes.

Alain Puig-Pérez: Recopilación de la información resultado de los instrumentos aplicados, confección de base de datos.

José Luis Rodríguez-Fonseca: Confección de base de datos.

Adonis Ramón-Puebla: Confección de instrumentos, recopilación de la información resultado de los instrumentos aplicados.

Esta obra está bajo una licencia de Creative Commons Reconocimiento-NoComercial 4.0 Internacional.
Copyright (c) 2021 Yenia Molina-Pelegrín, Gretel Geada López, Adonis Sosa-López, Alain Puig-Pérez, José Luis Rodríguez-Fonseca, Adonis Ramón-Puebla