Introducción

En las últimas décadas, la población de los 14 municipios que conforman la microcuenca Tula ha incrementado en promedio un 2% anual, con una población de 288 196 habitantes en 2010 (^{INEGI, 2013a}). Además, se incrementó la densidad de población 205 a 278 habitantes km^-2 (1990-2010), ubicándose por arriba de la media nacional 61 habitantes km^-2 (^{INEGI, 2013a}). Lo que ha intensificado la presión sobre los recursos naturales de la región, debido a la demanda de bienes y servicios de parte de la sociedad (^{Zamora et al., 2020}).

En la zona se ha registrado un crecimiento desordenado de áreas urbanas, desplazando áreas agrícolas y forestales, además de un aumento en el área de extracción de materiales para construcción (^{INEGI, 2017a}, ^2017b, ^2017c), lo que ha provocado un aumento en el cambio de uso del suelo.

De acuerdo a ^{Zamora et al. 2020}, en el periodo del año 2000 al 2017 hubo un incremento de áreas agrícolas de 4.17% en la microcuenca Tula, siendo esta una actividad económica que deterioran el recurso suelo, además en la microcuenca existen áreas abandonadas debido a la reducción o pérdida de su productividad, que de acuerdo a ^{Zavala-Cruz et al. (2011}) implica que el recurso suelo esta degradado. Por tal motivo, es necesario buscar alternativas de especies vegetales que tengan un óptimo crecimiento en la región y permitan además la rehabilitación de suelos degradados.

La especie M. atropurpureum se localiza desde el Sur de Estados Unidos de América y México hasta Centroamérica y Sudamérica, específicamente en Brasil, Colombia, Ecuador y Perú (^{McVaugh, 1987}). Mientras que, V. sativa L. es originaria de Europa, el norte de África y el oeste de Asia y considerada una especie exótica, en México (^{Vibrans, 2009}b).

Ambas especies pertenecen a la familia Fabaceae, a la sub familia Faboideae y son utilizadas para forraje, en México (^{Vibrans, 2009}a, ^2009b). De acuerdo con ^{Lewis et al. (2005}), la familia Fabaceae tiene una extensa distribución geográfica en el mundo; es decir, 727 géneros, con 19 327 especies de distribución cosmopolita, ubicadas principalmente en regiones cálidas y templadas. Mientras que, en México, 1 851 especies se han reconocido (^{Sousa et al., 2001}). M. atropurpureum y V. sativa son especies que se caracterizan por tolerar la sequía, competir con plantas indeseadas y fijar nitrógeno (^{La O et al., 2006}; ^{Flores et al., 2007}).

Sin embargo, hasta la fecha, no existe información acerca de su distribución geográfica potencial a nivel regional a pesar de que por sus características es posible emplearlas como herbáceas de usos múltiples. En México, la especie M. atropurpureum se localiza principalmente en potreros y bordes de caminos y se ha registrado en 28 estados, mientras que V. sativa se sitúa especialmente en bosques de pino encino y tiene registro en 8 estados (^{Flores et al., 2007}; ^{Vibrans, 2009a}, ^2009b).

Por otra parte, los factores que determinan la presencia de una especie en un área geográfica son, entre otros, las condiciones abióticas (clima, topografía, condiciones edáficas, etc.) es decir, son las que determinan su nicho ecológico (^{Austin, 2002}; ^{Soberon y Nakamura, 2009}). En los últimos 15 años, el desarrollo de programas para la modelación de nichos ecológicos ha sido un campo nuevo de investigación con alto crecimiento (^{Kearney, 2006}; ^{Lobo et al., 2010}).

Estos programas realizan modelos, basándose en el concepto de nicho ecológico de ^{Hutchinson (Hutchinson,1957}; ^{Guisan y Zimmermann, 2000}; ^{Elith et al., 2010}) relacionando la información biológica de la especie (presencia-ausencia y abundancia) y la información ambiental (^{Soberón y Peterson, 2005}; ^{Soberón, 2007}). BIOCLIM es uno de los primeros algoritmos usados para la modelación de nichos ecológicos, y está sustentado en funciones envolventes (^{Nix, 1986}; ^{Busby, 1991}); es decir, identifica áreas con características climáticas similares a las áreas en donde se ha registrado la presencia de una determinada especie, generando un intervalo ecológico para la especie conforme a las variables predictivas (^{Palacios et al., 2016}).

La modelación de nichos con BIOCLIM ha sido empleada en varios estudios entre los que se encuentran el de ^{Finch et al. (2006}) con invertebrados en África, ^{Aguirre-Morales et al. (2016}) con Passiflora subgénero Astrophea (Passifloreceae) en Colombia y ^{Palacios et al. (2016}) con Prosopis laevigata (Humb. et Bonpl. ex Willd) M. C. Johnston, en cada estudio se determinó la distribución potencial de la especie. Por lo anterior, el objetivo de esta investigación fue evaluar la distribución geográfica potencial de dos especies leguminosas forrajeras: Macroptilium atropurpureum (DC.) Urb y Vicia sativa L., en la microcuenca Tula, México; a través, de la modelación de sus respectivos nichos ecológicos, usando el algoritmo BIOCLIM, para determinar si éstas pueden tener un óptimo crecimiento en esta región para emplearse como alternativa de rehabilitación de suelos.

Materiales y métodos

Evaluación de la supervivencia en campo

De acuerdo con ^{Ramírez (2011}), la evaluación de la supervivencia en campo permite obtener una medida cuantitativa del éxito de la plantación bajo la influencia de los factores del sitio. La evaluación de la supervivencia de las especies V. sativa y M. atropurpurem se llevó a cabo durante los años 2018 y 2019, en el Ejido de San Felipe Teotitlán, que abarca parte de los municipios de Axapusco y Nopaltepec, en el Estado de México, dentro de la microcuenca Tula (Figura 1), el experimento se realizó en un área de 20 x 55 m de ancho y largo, respectivamente, ubicada entre las coordenadas geográficas 19° 50’ 24 y 19° 49’ 00’’ de latitud norte y 98° 42’ 00’’ de longitud oeste.

Figura 1 Localización del área experimental de la introducción de especies en Ejido San Felipe Teotitlán, México. 

La implementación del experimento consistió en producir plántulas de M. atropurpurem, la cual es una especie perene, en el invernadero del Departamento de Zootecnia en la Universidad Autónoma Chapingo (UACH), durante el periodo de mayo a julio del 2018. También, se realizó la preparación del terreno, comenzando con la roturación del suelo a curvas de nivel con un espaciamiento de 11 m, y la apertura de zanjas de 20 x 0.3 x 0.2 m de largo, ancho y profundidad, respectivamente.

Posteriormente, las zanjas se rellenaron con suelo de la región y zeolita (0.5 kg planta^-1), debido a la presencia de piedra y tepetate en los primeros 0.15 m de profundidad. Finalmente, las plántulas de M. atropurpurem fueron trasplantadas durante la primera semana del mes de agosto a una distancia de 0.2 m entre cada planta, en los primeros 10 metros de la zanja. Cabe mencionar que V. sativa, la cual es una especie anual fue sembrada directamente en campo en los 10 metros restantes de la zanja, en el mismo periodo y a la misma distancia.

En el mes de diciembre de 2018, se evaluó la supervivencia, definida como el número de plántulas que sobrevivieron al final del control de seguimiento, respecto al número de plantas al inicio del estudio este valor es expresado en porcentaje (^{Ramos-Font et al., 2015}). En el mes de junio de 2019, se repusieron todas las plantas muertas, esta vez la siembra fue directa en ambas especies y se evaluó la supervivencia, en el mes de octubre de 2019.

Resultados y discusión

Los registros de presencia obtenidos de las bases de datos MEXU y REMIB, para la especie M. atropurpureum fueron de 219 después del proceso de depuración de los 856 registros con el software QGIS 2.14.3 y para la especie V. sativa se obtuvieron 19 registros después de la depuración de los 124 registros.

Los registros de presencia indicaron que M. atropurpurem posee una extensa distribución en México, y se le identificó principalmente en los estados de Oaxaca (68 registros), Puebla (29), Campeche (25) y Chiapas (20), lo cual coincide con lo que reportado por ^{Vibrans (2009}a). Mientras que V. sativa se registró en siete estados, siendo el Estado de México (9 registros) el que más tiene, lo cual coincide con lo que reportó ^{Vibrans (2009b)}; ^{Flores et al. (2007)} (Figura 2).

Figura 2 Estados con registros georreferenciados de V. sativa y M. atropurpurem en Méx.  

De acuerdo con los valores obtenidos de los registros para V. sativa y M. atropurpurem, se obtuvieron los rangos de las variables ambientales en que se han reportado mayor número de individuos a nivel nacional. Para V. sativa, el intervalo de temperatura media anual obtenido fue de 13.8 a 17.3 °C (Cuadro 2 ), este rango indica que no tolera el frío para un crecimiento óptimo, tal como lo menciona ^{Alfonso (2007)}.

Respecto al rango de precipitación de los sitios en los cuales obtuvo mayor frecuencia de presencia, éste fue de 821 a 1 111 mm, con un intervalo general de 531 a 1 981 mm, rango que difiere al reportado por ^{Flores et al. (2007)}, quienes indican que crece y se desarrolla de 200 a 400 mm debido a que la especie soporta el estrés hídrico.

Para la especie M. atropurpurem, el mayor número de individuos se encontraron en una temperatura media anual de 24.5 a 27.6 °C, lo que indica que tiene menor tolerancia al frío que la especie V. sativa. Mientras que, el rango de precipitación anual fue de 733.2 a 1 396.4 mm, rango que está dentro del intervalo 650-1 500 mm reportado por ^{Vibrans (2009}a) para esta especie y entre el reportado por (^{Guerrero et al., 2009}) 650 a 1 750 mm, el rango general obtenido fue de 70 a 3 386 mm, lo cual explica su extensa distribución en México y su tolerancia al estrés hídrico.

Los rangos de las variables ambientales obtenidos para México y la información climática de la microcuenca Tula, la cual presenta una temperatura media anual de 14 °C y una precipitación media anual de 600 mm (^{INEGI, 2017a}, ^2017b, ^2017c), indican que V. sativa de acuerdo con la temperatura si se puede tener un adecuado crecimiento. Sin embargo, en la microcuenca no se tiene una óptima precipitación para el crecimiento óptimo de esta especie. En el caso de M. atropurpurem en la microcuenca no se tiene una temperatura media anual ni una precipitación media anual adecuada para su óptimo crecimiento.

Los resultados obtenidos de la distribución geográfica potencial mediante la modelación con BIOCLIM, para la especie V. sativa, indica que la microcuenca tiene una superficie total de 41.51 km² con nulo potencial, que equivale al 4% de la superficie, mientras que 996.15 km² (96%) tienen potencial; de esta superficie, 21% (209.19 km²) corresponde a bajo potencial y 79% (786.96 km²) a un potencial muy bajo (Figura 3).

Figura 3 Distribución potencial de V. sativa en la microcuenca Tula, México. 

Esto indica que en la microcuenca si se desarrollará la especie, pero no alcanzará su óptimo crecimiento y por ende rendimiento, ya que la especie requiere de una mayor precipitación media anual. La distribución geográfica potencial de la especie M. atropurpureum en la microcuenca, ésta fue de 3% (31.12 km²) con nulo potencial y de 97% (1 006.53 km²) con muy bajo potencial (Figura 4), lo que indica que la especie no tendrá un crecimiento óptimo en esta región, debido a que requiere de una mayor temperatura y precipitación media anual.

Respecto a la evaluación de la supervivencia de la especie en campo, la cual permite medir cuantitativamente el éxito de la plantación bajo la influencia de los factores del sitio (^{Ramírez, 2011}), durante el año 2018, ésta fue de 100% para M. atropurpurem y de 0% V. sativa, éste último valor se adjudicó a la presencia de fauna silvestre. Sin embargo, en la evaluación del año 2019, ambas especies registraron una supervivencia regular, con 66% para la especie M. atropurpurem y 70% para V. sativa.

Figura 4 Distribución potencial de M. atropurpureum en la microcuenca Tula, México. 

Ésta última especie compitió favorablemente contra las herbáceas anuales del sitio y alcanzó la formación de vainas, mientras que M. atropurpurem tuvo un crecimiento lento con una altura menor a 0.01 m en estado vegetativo. Cabe mencionar, que es posible mejorar su adaptación al mejorar los aspectos edáficos como lo hicieron en Zacatecas con la especie V. sativa, al aplicar riegos auxiliares, fertilizar y controlar malezas (^{Flores et al., 2007}; ^{Flores et al., 2016}). Sin embargo, por los costos que implican no se recomienda sean utilizadas en proyectos para esta región.

Finalmente, los valores del AUC, obtenidos para la validación de la precisión de predecir de los modelos, para M. atropurpureum y V. sativa, fue de 0.801 y 0.967 respectivamente, valores que indicaron que los modelos tienen una buena precisión de predicción (^{Peterson et al., 2001}; ^{Araujo y Guisan, 2006}; ^{Phillips y Dudík, 2008}).

Conclusiones

Los modelos de nicho ecológico obtenidos con el algoritmo de BIOCLIM para las especies M. atropurpurem y V. sativa tienen una buena precisión de predicción para la obtención de la distribución potencial de las especies. El bajo potencial de distribución obtenido en la modelación y la regular supervivencia de las especies leguminosas M. atropurpurem y V. sativa indicó que estas especies no son las más aptas para introducirse en la microcuenca Tula; sin embargo, hace falta realizar estudios acerca de la adaptación de estas especies a la región al mejorar los aspectos edáficos del sitio de introducción.