Introducción
La producción de alimentos se está convirtiendo en un verdadero reto en un clima que está cambiando. Bajo este contexto, la diversificación de la alimentación animal es un aspecto clave para adaptarse mejor al cambio climático. Este es el motivo para evaluar el impacto del cambio climático sobre la presencia y distribución potencial de las especies forrajeras. Este es el caso específico de T. diversifolia, nativa y con distribución actual en las tierras bajas del sudeste de México, Centro y Sudamérica, especie que se considera un importante recurso genético y una planta excepcional como alternativa en la alimentación animal1,2. El riesgo futuro para muchas especies de plantas debido al cambio climático no ha sido bien establecido aún.
La temperatura global ha aumentado 0.85 °C en el periodo de 1901 al 20113. Se ha establecido sobre la base de los modelos de predicción del clima, que la presión de saturación del vapor ha sido muy sensible a los cambios de temperatura, por lo que se espera que se manifiesten futuras rupturas en el ciclo global del agua4. Por otro lado, se ha mencionado que la disponibilidad de agua presentará una marcada reducción anual en el sureste de Estados Unidos de Norteamérica, el Caribe y diversas partes de México5. Además, la agricultura se expande e intensifica principalmente en los trópicos6; se estima que alrededor del 80 % de la nueva tierra de cultivo será a base de sustitución de bosques7. Por lo tanto, al retrasar o retener la expansión de la agricultura en los trópicos se reducirán las emisiones de carbono y la pérdida de la biodiversidad8.
Con el cambio climático, la diversidad de climas puede ser modificada, así como la composición de los ecosistemas9, que incluirían alteraciones en la distribución, fenología y un aumento de riesgo de especies en peligro de extinción10. Se considera que varios componentes del cambio climático afectarán a todos los niveles de la biodiversidad; es decir, desde los organismos y poblaciones hasta las áreas bióticas11. En los niveles básicos, el cambio climático es capaz de disminuir la diversidad genética en las poblaciones debido a una selección direccional y una migración dinámica, que podría interferir en la capacidad de resiliencia y la funcionalidad de los ecosistemas12; por lo tanto, podría causar la modificación de la red de interacciones a nivel de la comunidad13. Además, un efecto del cambio climático podría inducir la invasión de especies potencialmente peligrosas para un ecosistema14; afectando la fisiología, fenología y el comportamiento de especies15.
De acuerdo con lo anterior, México está interesado en la construcción de un inventario en lo que respecta a los efectos previstos del cambio climático en las especies silvestres deseables. Debido a eso, el objetivo fue estimar el impacto del cambio climático en la distribución potencial de T. diversifolia.
De acuerdo con diferentes autores los usos reportados de Tithonia incluyen producción de forrajes, control de erosión del suelo, como material de construcción y un refugio para las aves16. Como contrapartida T. diversifolia es considerada como una hierba alelopática con aleloquímicos solubles en agua en partes de la planta y con tal potencia fitotóxica, que podría suprimir el crecimiento y la acumulación de nutrientes en cultivos asociados17.
Material y métodos
Base de datos
Se consideraron los datos de 52 sitios de población de T. diversifolia, con distribución actual en las tierras bajas del sudeste de México1,2. Las bases de datos se obtuvieron a partir de dos fuentes: Inventario Nacional Forestal y de Suelos (INFyS) de la Comisión Nacional Forestal de México (CONAFOR) y del sitio http://www.tropicos.org/.
Áreas potenciales de distribución
En este estudio se utilizó el modelo de MaxEnt (máxima entropía)18 para modelar el nicho ecológico y predecir la distribución geográfica más probable de T. diversifolia. MaxEnt ha sido ampliamente utilizado para estimar distribuciones potenciales de plantas y animales con alta precisión; algunos grupos terrestres son excelentes ejemplos19.
En el modelo MaxEnt la distribución de una especie se representa mediante una función de probabilidad P sobre un conjunto X de sitios en el área de estudio. Se construye un modelo de P mediante un conjunto de restricciones derivado de datos empíricos de presencia de la especie20. Las restricciones se expresan como funciones simples de variables ambientales conocidas. El algoritmo de MaxEnt obliga a que la media de cada función de cada variable se aproxime a la media real en las áreas de presencia de la especie. De todas las posibles combinaciones de funciones, se selecciona la que logra minimizar la función de entropía; la cual es medida con el índice de Shannon. La expresión general de la función de probabilidad para i variables ambientales es21:
Donde:
P(x)= función de probabilidad;
λ= vector de coeficientes de ponderación;
f= vector correspondiente de funciones de variables ambientales;
Z= constante de normalización para asegurar que P(x) = 1.
Los valores P(x) así obtenidos, representan valores de idoneidad relativa para la presencia de la especie, constituyendo así la base de un modelo de distribución potencial. Un alto valor de la función de distribución en cada píxel indica que tiene condiciones favorables para la especie22.
En este estudio con el uso de capas de parámetros ambientales y los datos de presencia de especies, se generó un modelo. El 75 % de los registros de presencia se utilizaron como puntos de entrenamiento y 25 % como puntos de validación, además se utilizó el índice AUC (área bajo la curva), para evaluar el modelo estadístico, ya que este índice es de los más utilizados para medir la calidad de los modelos20. Los ajustes que se pueden determinar en el modelo de MaxEnt, afectan su precisión, ya que el modelo puede trabajar con variables ambientales simples o complejas.
Bases de datos y parámetros ambientales
Los datos mensuales, estacionales y anuales de precipitación, temperatura máxima, temperatura mínima, temperatura media y oscilación térmica de los períodos 1950-2000 (climatología de referencia) y 2041-2060 y 2061-2080 (climatología futura), se utilizaron para una caracterización ambiental de los sitios de presencia de las especies. Estos datos climáticos se obtuvieron del portal Global Climate Data de WorldClim y las imágenes se trabajaron en formato ascii y raster, con una resolución de 2.5 min de arco. Para los períodos 2041-2060 y 2061-2080, el MCG HadGEM2-ES fue considerado con dos vías de concentración representativas de gases de efecto invernadero: Rcp 4.5 y Rcp 8.5. Este modelo ha sido utilizado en los nuevos escenarios de cambio climático para México; estudio presentado en el quinto informe de evaluación sobre el cambio climático IPCC (Grupo Intergubernamental de Expertos sobre el Cambio Climático)23. Otras variables que se incluyeron en la modelación de la distribución de T. diversifolia, fueron la altitud y la pendiente del suelo, imágenes que se obtuvieron del Sistema Nacional de Información Ambiental (SIAN) del Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias24. La lista de variables utilizadas para el modelo de nicho espacial fueron: pendiente del suelo (%), altitud (m), temperatura máxima media anual (Txan), temperatura máxima media de mayo a octubre (tx5-10), temperatura máxima media noviembre-abril (tx11-4), temperatura máxima del mes más cálido (txmc), temperatura media anual (tman), temperatura media de mayo a octubre (tm5-10), temperatura media noviembre-abril (tm11-4), temperatura mínima del mes más frío (timf), temperatura media del mes más frío (tmmf), temperatura mínima promedio anual (tian), temperatura media mínima de mayo a octubre (ti5-10), temperatura media mínima noviembre-abril (ti11-4), precipitación acumulada anual (pa), precipitación acumulada mayo a octubre (pa5-10) precipitación acumulada noviembre-abril (pa11-4), precipitación acumulada del mes más seco (pms), precipitación del mes más húmedo (pmh), oscilación térmica media anual (otan), oscilación térmica de mayo a octubre (ot5-10) y oscilación térmica noviembre a abril (ot11-4).
Áreas con probabilidad de aptitud ambiental
El modelo obtenido de predicción MaxEnt para la distribución potencial de T. diversifolia fue examinado con el sistema de Idrisi Selva25, y se hizo un mapa con los valores del umbral de los píxeles correspondientes al percentil 1026. Para el cálculo de la superficie potencial de distribución de la especie, las áreas ocupadas por las masas de agua y los centros urbanos no fueron consideradas. Estas capas temáticas se obtuvieron mediante la manipulación del uso del suelo y la capa de vegetación27.
Descriptores ecológicos
Los descriptores ecológicos para T. diversifolia se determinaron sobre la base de los intervalos ambientales derivados de la caracterización de los parámetros de sitios de presencia de la especie. Esto se hizo utilizando el sistema IDRISI y con las imágenes raster de cada variable ambiental y las coordenadas geográficas de cada sitio de presencia.
Resultados y discusión
En el modelo del nicho de distribución potencial de T. diversifolia, la curva operacional mostró resultados satisfactorios (Figura 1). Al considerar los datos de entrenamiento (75 %), el área bajo la curva (AUC) de T. diversifolia fue de 0.979; mientras que cuando se utilizan los datos de prueba (25 %), la validación del modelo reportó una AUC de 0.966, lo que indica que la capacidad del modelo para representar la presencia de especies fue satisfactoria28.
La Figura 2 muestra la distribución geográfica actual de T. diversifolia. Su presencia se encuentra principalmente en la región Centro-Sur, Golfo de México, las zonas costeras del Pacífico y la Península de Yucatán; es decir en la parte tropical del país (<23° 26' 14'' Norte). La mayor parte de las poblaciones de T. diversifolia se concentran en la parte sureste de México, que está en correspondencia con la aseveración acerca de su origen geográfico (México y América Central y del Sur)2. Zonas tropicales y subtropicales de todo el mundo, lo que le ha permitido, crecer bajo una amplia gama de condiciones ambientales y, por tanto, desarrollar un amplio rango de adaptación29.
Los descriptores ecológicos de T. diversifolia y las variables ambientales que contribuyen a la construcción del modelo se pueden observar en el Cuadro 1. Según estos resultados, la distribución de T. diversifolia se determina principalmente por la variable pa5-10 (78.2 %), esto es, la cantidad de precipitación durante el período mayo-octubre, la cual es un elemento clave para la presencia de esta especie, así como para otras especies, en donde la precipitación en pequeños pero frecuentes eventos durante el verano es clave para el buen desarrollo30. El descriptor de T. diversifolia para esta variable indica un rango de 297 a 3,404 mm (para la precipitación anual de 356 a 3,828 mm), lo que da una amplia variación de la precipitación espacial. Estos resultados son consistentes en el sentido que esta especie crece adecuadamente bajo una amplia gama de precipitación. Sin embargo, se han establecido límites de precipitación anual de 600 y 5,000 mm30. Por lo anterior, se puede concluir que T. diversifolia puede crecer incluso con menos precipitación estacional y anual de lo que indica la literatura actual.
Variable | Valor mínimo | Valor máximo | Contribución (%) |
---|---|---|---|
pa5-10, mm | 297 | 3404 | 78.2 |
Txmc, °C | 21.9 | 34.4 | 5.9 |
Pendiente, % | 0 | 14 | 3 |
Tmmf, °C | 11 | 26.15 | 2.8 |
Pmh, mm | 67 | 713 | 2.1 |
El segundo descriptor ecológico más contribuyente fue txmc (temperatura máxima media del mes más cálido) y T. diversifolia también mostró tener un amplio rango térmico; 21.9 a 34.4 °C. Teniendo en cuenta que T. diversifolia también tiene un amplio rango para el descriptor timf (temperatura media mínima del mes más frío, 1.8 a 19.8 ° C, Cuadro 1) se puede concluir que esta especie es capaz de explorar y colonizar ambientes térmicamente extremos. Otros autores2 habían indicado previamente que esta especie está adaptada a diferentes climas y altitudes.
Los resultados de los nichos potenciales para T. diversifolia bajo la climatología de referencia (1951-2000) así como las futuras climatologías, también se muestran en la Figura 2. Los mapas de esta figura muestran la presencia de áreas adecuadas, inadecuadas, de retracción y expansión para Tithonia diversifica en comparación con los escenarios de la climatología de referencia y las climatologías futuras. En todos los escenarios futuros aparecerán áreas de retracción, así como las áreas de expansión; mostrando la posible dinámica de las áreas con aptitud ambiental en función de los cambios climáticos futuros.
El umbral logístico de entrenamiento de presencia del percentil 10 estima un área con aptitud ambiental (para la climatología de referencia) que representa el 30.71 % de la extensión territorial de México. Después de un balance entre las áreas de contracción y expansión de las climatologías de cambio climático, las áreas óptimas se incrementarán al 31.62 %, 31.83 %, 32.45 % y 32.45 % del territorio mexicano, en los escenarios 2041-2060 Rcp8.5, 2041-2060 Rcp4.5 y 2061-2080 Rcp4.5, 2061-2080 Rcp8.5, respectivamente. Estos resultados muestran que los cambios climáticos futuros serán aparentemente benéficos para T. diversifolia. Al considerar que las variables más significativas en los modelos de los nichos de distribución son la precipitación acumulada de mayo a octubre y la temperatura máxima media del mes más cálido, se podría inferir que las combinaciones de estos parámetros bajo escenarios de cambio climático tendrán un efecto positivo en las áreas con aptitud ambiental para T. diversifolia.
En algunos estudios se ha señalado que el efecto del impacto del cambio climático sobre la distribución de las especies, pueden resultar en especies ganadoras y perdedoras31. El resultado final va a depender de su adaptación evolutiva. Por ejemplo, se han reportado expansiones de áreas adecuadas debido al cambio climático para otras especies como Leucospermum hypophyllocarpodendron subsp. hypophyllocarpodendron11.
La mayoría de los estudios sobre los efectos del cambio climático en la distribución de especies han reportado impactos negativos, siendo las áreas de retracción mayores que las áreas de expansión32. En el presente estudio se obtuvo lo contrario; es decir, áreas de expansión mayores que las áreas de retracción. Como se puede visualizar en la Figura 2, estos dos tipos de áreas tienen un patrón geográfico; las superficies de contracción se ubican principalmente en la región occidental del país y cerca de la costa del Pacífico, mientras que las de expansión se localizan principalmente en el lado oriental y cerca de la costa del Golfo de México. Este hecho apunta a dos regiones del país donde los futuros patrones climáticos difieren en las condiciones ambientales y, por lo tanto, difieren en el nivel de idoneidad para T. diversifolia. Esta situación en particular ya se ha manifestado, ya que la región correspondiente a las áreas de retracción ha sido reportada previamente con cambios en los patrones de cultivos como consecuencia de los cambios climáticos regionales33,34.
Aparentemente en los mapas que se ejemplifican en la Figura 2 ninguna población actual de T. diversifolia se ve afectada por la retracción de la aptitud ambiental. Sin embargo, en el caso de que algunas poblaciones de T. diversifolia se ubiquen en áreas de retracción en los próximos años, se debe considerar que estas poblaciones tendrán sólo dos opciones para sobrevivir: a) migrar a ambientes más favorables o; b) adaptarse a las nuevas condiciones climáticas35, que dependerá principalmente de la diversidad genética de las especies36,37. En el caso del presente estudio, se puede hipotetizar que la idoneidad medioambiental ganará un poco más de superficie debido al cambio climático, y esto podría aparentemente promover la futura dispersión de T. diversifolia a nuevas áreas. No obstante, es importante considerar las posibles capacidades de esta especie para competir contra otras del ecosistema o con especies invasoras que constituyen una amenaza para la estabilidad del ecosistema38,39,40. Además, es importante considerar que T. Diversifolia en el futuro también dependerá de su capacidad de adaptación al cambio climático41, lo que está en función de su capacidad para colonizar nuevas áreas o (cuando sea necesario) su capacidad de implementar modificaciones fisiológicas para adaptarse al nuevo entorno42. Al adaptarse al cambio climático y a los nuevos climas, la plasticidad ecológica de T. diversifolia1,30 constituye un aspecto clave43).
De acuerdo con los resultados obtenidos, el cambio climático previsto para ambos períodos tendrá un efecto más positivo que negativo sobre la idoneidad ambiental de T. diversifolia (Figura 3), permitiendo así en el futuro su posible expansión territorial. Sin embargo, está claro que esta ventaja ambiental que será causada por el cambio climático futuro, resultará también favorable para otras especies tal vez en un contexto de competencia con T. diversifolia. Bajo este tipo de escenario, el amplio intervalo para la mayoría de los descriptores ecológicos de T. diversifolia (Cuadro 1) podría ser una característica ventajosa, aún más si se considera la diversidad genética de esta especie en otras partes del mundo, lo cual podría aumentar las capacidades de adaptación y colonización de T. diversifolia36,44.
No obstante que los escenarios climáticos actual y futuros parecen ser favorables para T. diversifolia; deben adoptarse medidas para preservar las poblaciones actuales de T. diversifolia a fin de asegurar su presencia en su medio natural. Por otro lado, dado que se ha encontrado que T. diversifolia es una maleza alelopática con aleloquímicos solubles en agua en varias partes de la planta17, se debe prestar atención a supervisar su posible expansión territorial debido a posibles afectaciones a los cultivos asociados. También debe monitorearse la lucha contra la presencia de T. diversifolia como maleza en los cultivos, lo cual podría amenazar las poblaciones naturales de esta especie.
Conclusiones e implicaciones
Actualmente, las áreas con aptitud ambiental para T. diversifolia representan cerca de un tercio del territorio mexicano; adaptándose a una amplia gama de condiciones de temperatura y humedad. Su presencia y distribución espacial está determinada por la cantidad de precipitación acumulada durante el período mayo-octubre, con un rango de 297 a 3404 mm. Los cambios climáticos futuros provocarán tanto retracción como expansión de las áreas con aptitud ambiental en diferentes partes del país; haciendo un análisis entre estos dos efectos se concluye que en el presente siglo se incrementarán ligeramente estas áreas, entre 0.91 y 1.74 % según los modelos utilizados. Por lo anterior, se puede considerar que esta especie puede ser una buena alternativa de producción de forraje para el futuro. Por otro lado, dado que se ha encontrado que T. diversifolia es una hierba alelopática, se debe poner atención en la supervisión de su posible expansión territorial debido a las posibles afectaciones a los cultivos asociados. La lucha contra la presencia de T. diversifolia en los campos de cultivo también debe supervisarse para que no constituya una amenaza para las poblaciones naturales de esta especie.