Introducción

La Tierra absorbe radiación solar y esta es redistribuida a través de la circulación de la atmósfera y el océano de la zona ecuatorial hacia los polos, la energía es irradiada de la superficie al espacio hasta alcanzar el equilibrio entre la radiación absorbida y emitida, sin embargo, la variación de estos factores puede derivar en alteraciones en el clima (Magaña, 2004).

Las emisiones de gases de efecto invernadero derivadas de las actividades humanas han aumentado desde la era preindustrial debido al crecimiento demográfico y económico. De acuerdo con el Panel Intergubernamental de Expertos sobre Cambio Climático, IPCC, por sus siglas en inglés (2014) entre el año 2000 y 2010 las emisiones de gases de efecto invernadero alcanzaron un máximo histórico, las concentraciones de bióxido de carbono (CO₂), metano (CH₄) y óxido nitroso (NO) lograron los niveles más altos en 800 mil años, a partir de la revolución industrial y como consecuencia, se ha incrementado la temperatura global.

Los efectos del cambio climático se han reflejado, en el incremento de la temperatura media anual entre 1880 y 2012 de 0.85°C; el nivel medio del mar ha aumentado, en promedio, 3.2 mm./año en las últimas dos décadas, si se considera el escenario de mayores emisiones de gases de efecto invernadero, al año 2100, la temperatura global podría elevarse hasta 4°C y un incremento en el nivel del mar de casi un metro con respecto a lo observado entre 1980 y 1999. a fines del S. XXI (Organización de las Naciones Unidas, 2017).

Para estudiar las posibles consecuencias de la alteración de la composición de la atmósfera, el Programa de Naciones Unidas para el Medio Ambiente (PNUMA) en conjunto con la Organización Meteorológica Mundial (OMM), fue constituido el IPCC que a la fecha ha publicado cinco reportes cuyo objetivo es reunir el conocimiento científico, socioeconómico y técnico sobre el cambio climático (IPCC, 2015)

En México, el Instituto Nacional de Ecología y Cambio Climático (INECC) es el órgano gubernamental creado a partir de la Ley General de Cambio Climático para coordinar y llevar a cabo proyectos de investigación científica y tecnológica en conjunto con instituciones académicas públicas y privadas. A partir de ello surge el proyecto Actualización de escenarios de cambio climático para estudios de impactos, vulnerabilidad y adaptación en México y Centroamérica, financiado por el Fondo de Medio Ambiente Mundial (FMAM) y administrado por el Programa de Naciones Unidas para el Desarrollo (PNUD). (Cavazos, Salinas et al., 2013; Quinta Comunicación Nacional, 2012 en Fernández et al., 2015).

Dicho proyecto incluye 15 modelos generales de circulación atmosférica elaborados por instituciones de diversas partes del mundo e integran a la Quinta Comunicación Nacional y que a su vez fueron utilizados para la elaboración del modelo Reliability Ensemble Averaging, REA, así mismo, se proyectaron datos de temperatura y precipitación a tres horizontes temporales: cercano (2015-2039); intermedio (2045-2069) y lejano (2075-2099), además se generó un sistema interactivo de consulta y descarga de bases de datos de código abierto. (Fernández et al. 2015; INECC, 2016).

El modelo REA presenta dos trayectorias representativas de concentraciones o forzamientos radiativos; la 4.5 (bajas emisiones) y la 8.5 (altas emisiones) y una resolución espacial de 30”x30” (Fernández et al., 2015).

El bosque de oyamel

El ecosistema del bosque de oyamel brinda refugio a gran cantidad de biodiversidad, en él se han logrado identificar hasta 510 especies de plantas, sin embargo, solo el 8 % de ellas se encuentran de manera común en estas áreas (Sánchez et al., 2005).

El oyamel, además, sirve como sumidero de CO₂, de acuerdo con la Procuraduría de Ordenamiento Territorial (2010) el promedio de almacenamiento de carbono registrado en el suelo de conservación de la Ciudad de México es de 93.4 ton/ha, en tanto que, Razo et al. (2015) encontraron que los bosques conservados presentaron una menor capacidad de secuestro de carbono (62.6 ton/ha.) en contraste, los bosques alterados con árboles jóvenes alcanzan hasta 138.6 ton/ha. en el PN El Chico, Hidalgo.

El hábitat del oyamel se encuentra restringido y disperso debido principalmente a que requiere de condiciones específicas de temperatura y humedad, se le considera como un bosque relicto de periodos glaciales (Rzedowski & McVaugh en Guerrero et al., 2014).

Los bosques de Abies religiosa se localizan en altitudes superiores a los 2400 hasta los 3600 msnm o incluso 3800 msnm, en laderas protegidas de la radiación solar y de los fuertes vientos, así mismo, donde la humedad es alta y la lluvia apreciable se presenta en al menos 100 días a lo largo del año (Challenger 1998; Rzedowski, 2006; Villers et al., 2006) (Figura 1).

Díaz et al. (2011) indican los requerimientos ambientales del bosque de oyamel en el cual consideran cinco variables distintas y establecen cuatro categorías de acuerdo con lo favorable de las condiciones (Tabla 1).

El área de estudio se delimita al Norte con el Parque Nacional El Chico (20°13’ 25”N); al Sur Parque Nacional Izta-Popo (18°59’00”N); al Este el lago de Cuitzeo (100°59’10”) y al Oeste el Parque Nacional Pico de Orizaba (97°12’00”) (Figura 1).

Metodología

La línea base de las condiciones climáticas se construyó a partir de los datos mensuales de las estaciones climatológicas de CLICOM del Servicio Meteorológico Nacional. El periodo de tiempo considera una dimensión temporal de 1950 hasta 2010, así como al menos 20 años de datos continuos y 90% de disponibilidad de los mismos.

La validez de los datos depende principalmente de tres factores; 1) tamaño de la muestra; 2) homogeneidad; 3) coincidencia temporal a fin de contrastarlos (Bello, 2014).

Se utilizó el método estadístico de prueba de homogeneidad por rachas bajo el criterio de Doorembos. Se aplica a los datos de precipitación y consiste en dos etapas: 1) obtención de media y median de la serie; 2) cálculo de desvíos o cambios de signo. Si el número de desvíos se encuentra dentro del rango establecido por Doorembos, la serie se considera homogénea, de lo contrario la estación se descarta (Cruz et al., 2012).

Se calcularon datos de temperatura media, máxima, mínima y precipitación a partir de la determinación del gradiente altitudinal mediante ecuaciones de regresión lineal simple para puntos aleatorios cercanos a la cima del volcán Tláloc, Iztaccíhuatl, Popocatépetl y La Malinche.

El total de estaciones climatológicas, incluyendo los puntos aleatorios calculados fueron 168; Ciudad de México, 20; Estado de México, 56; Morelos, 8; Tlaxcala, 15; Puebla, 21; Hidalgo, 25; Veracruz, 12; Michoacán, 11 (Figura 2).

La base de datos quedó integrada por un total de cinco variables climatológicas:

Fuente: CLICOM, INEGI

Se trabajó con el modelo de circulación regional mexicano REA con horizonte lejano (2075-2099) con los forzamientos radiativos disponibles (4.5 y 8.5) y resolución de 30” x 30” por píxel. Los datos mensuales de las variables climatológicas se integraron a un entorno SIG, de tal manera que se obtuvieron los datos anualizados correspondientes a las 168 estaciones.

Además, se seleccionaron tres de las estaciones ubicadas en inmediaciones del bosque. Se analizó el comportamiento de las condiciones de temperatura y precipitación a lo largo del año para así contrastar las normales climatológicas con las proyecciones del modelo REA mediante climogramas de Ball.

Para la distribución potencial se analizó mediante el software MaxEnt 3.3.3. el cual basa su funcionamiento en el principio de máxima entropía, que consiste en la identificación de las relaciones existentes entre la localización real de las especies y las características biofísicas y ambientales, de tal manera que es posible estimar la probabilidad de la distribución y cuya ecuación básica se expresa del siguiente modo:

Donde π es la distribución probable desconocida, x como el conjunto de píxeles en el área de estudio. La distribución π asigna una probabilidad no negativa π(x) para cada píxel x sumando 1. π también se define como una distribución probable la cual queda definida como ; ln indica logaritmo natural. La entropía es no negativa y es como máximo, el logaritmo natural de los elementos contenidos en X (Kumar & Stohlgren, 2009; Martínez, 2010; Phillips et al., 2006).

Las tablas de puntos georreferenciados se obtuvieron a partir de la capa vectorial de Uso de Suelo y Vegetación Serie V escala 1: 250 000 de INEGI y por el otro, archivos ráster en formato ASCII los cuales contienen todas las variables climatológicas a trabajar.

Estos archivos ráster son el resultado de la interpolación de las variables con la ayuda de Sistemas de Información Geográfica (SIG) ya que esta herramienta permite trasladar al espacio geográfico los datos aportados por las estaciones.

El método de interpolación trabajado fue Kriging simple con superficie predictiva y distribución base de t-student, esto debido a que presentó los mejores niveles de desviación estándar a comparación de otros métodos analizados en cada una de las variables (Tabla 2).

Los productos obtenidos por el procesamiento incluyen entre otros, archivos de formato ASCII cuyos valores van de 0 a 1, donde los valores cercanos a cero representan la menor probabilidad de distribución, mientras que los valores cercanos a 1 son aquellos cuya potencialidad es mayor. El formato de salida es compatible en el ambiente SIG, con lo que es posible realizar el análisis espacial de los datos.

La categorización de los valores de cada escenario se realizó tomando como referencia de la propuesta por Díaz et al. (2011) a la que se añadió una quinta; a) muy adecuado; b) adecuado; c) moderadamente adecuado; d) poco adecuado; e) no adecuado.

Resultados

Análisis puntual de estaciones climatológicas.

Fuente: Elaboración propia.

La estación climatológica “Desierto de los Leones” (Álvaro Obregón, Ciudad de México), se encuentra a 2995 msnm, reporta un periodo de funcionamiento de 26. 7 años (1961-1987) y un porcentaje de datos total del 99 % (Figura 2 (A)).

La temperatura media anual de la línea base es de 10.6 °C, mayo representa la temperatura más alta del año (12.5 °C) en contraste con enero que es el mes más frío (8.2 °C); ello significa una amplitud térmica de 4.3 °C.

La estación lluviosa se establece a partir del mes de abril y se extiende hasta octubre. La precipitación acumulada anual es de 1316.9 mm y las lluvias más copiosas ocurren durante el mes de agosto (284.1 mm) y el de menor cantidad es el de diciembre (12.4 mm); la temporada de estiaje cubre un periodo de 4 meses.

En el escenario propuesto REA 4.5 la temperatura media anual muestra un aumento de 3.4 °C respecto al promedio; mayo se mantiene como el mes más cálido del año con 16.5 °C, en contraste, el más frío es enero con 11.1 °C lo que representa un incremento de 4 °C y 2.9 °C respectivamente; con una amplitud térmica de 5.5 °C a lo largo del año

En cuanto a la precipitación anual muestra una disminución de 9.5% al ubicarse en 1191.2 mm. El máximo de lluvias se presenta durante agosto con una disminución de del 9% con respecto a la línea base; el mes más seco se traslada a febrero, al reducirse un 60% en relación con la normal actual. La temporada seca del año muestra un alargamiento que pasa de cuatro meses en la actualidad a seis meses. En tanto que la extensión de la estación lluviosa se mantiene desde mayo hasta octubre.

El escenario REA 8.5 exhibe condiciones más extremas en general; la temperatura media anual se ubica en los 15.7 °C, es decir, 4.9 °C arriba del escenario base; mayo alcanza los 18.6 °C lo que representa 6.1 °C por arriba del mismo mes en la actualidad. Por su parte, el mes más frío se mantiene en enero con 12.6 °C, 4.4 °C más. La amplitud térmica alcanza los 6 °C.

La lluvia acumulada anual exhibe una reducción de 15.6% al ubicarse en 1191.2 mm; el mes más lluvioso se traslada al mes de septiembre donde la disminución es del 1 % con 237.6 mm, en contraste con el mes de agosto que muestra una disminución de 17.9 %; el mes más seco es febrero donde la precipitación se reduce al 100 %. La sequía se presenta durante seis meses a lo largo del año y la temporada de lluvias se mantiene desde mayo hasta octubre.

La segunda estación utilizada para el análisis puntual es “Presa Iturbide” (Isidro Fabela, Estado de México) se ubica a 3290 msnm, posee 33 años de funcionamiento (1977-2010) y un porcentaje de datos del 95.1 % (Figura 2 (B)).

La estación reporta una temperatura media anual de 8.9° C, la temperatura máxima se presenta durante mayo con un valor de 10.9 °C, mientras que la temperatura mínima se observa en enero alcanzando 6.3 °C. La amplitud térmica se observa de 4.6 °C.

La lluvia acumulada anual se ubica en los 1529.2 mm y la mayor precipitación se registra en durante julio con 304.8 mm, sin embargo, el más seco es diciembre con 14.2 mm. La mayor parte del año puede considerarse dada la relación temperatura-precipitación como húmedo, exceptuando diciembre, aun así, es posible distinguir un incremento significativo en las precipitaciones durante mayo, en lo que puede considerarse el establecimiento de la temporada de lluvias la cual se extiende hasta octubre, momento en el que se percibe el importante descenso de las mismas.

El escenario REA 4.5 da como resultado el incremento de la temperatura media anual hasta 12.2 °C, un incremento de 3.3 °C; el mes más cálido se presenta en mayo con 14.1 °C, 3.2 °C más; la temperatura mínima se aprecia durante enero con un valor de 9.6 °C, 3.3 °C por arriba de la cifra actual; la magnitud de la amplitud térmica es de 4.5 °C.

La lluvia acumulada presenta un descenso en términos anuales de 7.9 % con respecto a la línea base, con 1408.1 mm; agosto, por su parte, se coloca con la mayor cantidad de precipitación con 285.1 mm; julio que en la línea base es la mayor cantidad de precipitación disminuye 10 % hasta 274.3 mm; en tanto, el mes con la menor cantidad de lluvia es diciembre con 5.1 mm, 65 % menos.

En el escenario REA 8.5 la temperatura media anual es de 14.4 °C, 5.5 °C más que en la línea base; la temperatura máxima se presenta durante abril y mayo con 16.7 °C, por su parte, la mínima se observa en enero con 11.6°C, incrementos de hasta 6.2 °C y 5.3 °C respectivamente. Ello representa que la amplitud térmica se ubica en 5.1 °C.

La precipitación advierte un porcentaje de disminución de 14.6 % en el acumulado anual; el mes que concentra la mayor cantidad de lluvia se desplaza de julio, donde se ubica en la línea base, a septiembre donde se observa incluso un ligero incremento con respecto al valor inicial (0.9 %), en contraste con julio que se redujo 18.9 %. Diciembre se presenta como el más seco en donde se observa una reducción de 81.5 % con respecto a la normal actual; los meses secos sumen un total de cinco, mayo marca el inicio de la temporada húmeda que se extiende hasta octubre. Destacan los meses de septiembre y octubre en los que se advierten leves incrementos en los acumulados de lluvia.

En tanto la estación “San Antonio Arroyo Prieto” (Ixtacamaxtitlán, Puebla), se localiza a 3065 msnm tiene un periodo de funcionamiento 54.9 años (1954-2009) y concentra el 92.4 % de los datos (figura 2(C)).

La temperatura media anual actual registra 10 °C; abril y mayo son los meses donde se presenta la temperatura más alta con 12.1 °C mientras que la mínima se encuentra en enero con 8.1 °C, de tal manera que la amplitud térmica es de 4 °C.

La precipitación de la estación registrada fue de 998.4 mm anuales, en donde el mes más lluvioso es junio con 176.4 mm y el menos lluvioso es febrero con 16.3 mm. En total se observan 10 meses húmedos por dos secos; las altas precipitaciones comienzan a observarse a partir de mayo hasta octubre. Además, se presenta canícula durante en el mes de julio para alcanzar el segundo máximo en septiembre.

Las condiciones en el modelo REA 4.5 presentan una temperatura media anual de 13.2 °C, 3.2 °C más en comparación a la base actual; la temperatura máxima se advierte durante abril, con 15.5 °C, lo que significa una diferencia de 3.4 °C a la normal existente; la temperatura mínima ocurre en enero con 10.3 °C, es decir, un aumento de 2.2 °C. La amplitud térmica se ubica en 5.2 °C.

La lluvia experimenta un decrecimiento de 8.8 % al alcanzar la cantidad de 910 mm anuales; el mes más lluvioso se presenta en junio, donde la variación es de 3.9 % negativo, equivalente a 169.4 mm; en contraste el más seco es diciembre donde se presenta una disminución de 50 % de la precipitación al quedar en 8.4 mm. El periodo de aridez se extiende a lo largo de cuatro meses por ocho de humedad; la temporada de lluvias copiosas abarcan de mayo hasta octubre, este último mes muestra un incremento de 3.9 % para alcanzar 93.2 mm. Se presenta un periodo de sequía intra estival durante julio, pero más acusada que en la actualidad mientras que el segundo punto máximo de lluvias se mantiene en septiembre.

En el modelo REA 8.5 establece una temperatura media anual de 14.9 °C, lo que representa una diferencia de 4.9 °C con respecto a la normal actual; la temperatura máxima se ubica durante mayo con 17.2 °C, 5.1 °C por arriba; la temperatura mínima se presenta en enero con 11.6 °C lo que significa un incremento de 3.5 °C. La amplitud térmica es de 5.6 °C.

El caso de la precipitación tiene una variación negativa de 16.1 % con una cantidad de 837.8 mm de lluvia anual; la mayor cantidad de lluvia se presenta durante junio, sin embargo, tiene una disminución de 6 % al presente; por su parte, el más seco es enero con solo 1.2 mm, es decir, disminuye 93.2 %. En este escenario se observan un total de cinco meses de aridez, por siete húmedos; el establecimiento de la temporada de lluvias se aprecia desde mayo hasta octubre, este último presenta también un incremento en la cantidad acumulada, de 21.2 % para registrar 113 mm. Al igual que en los otros escenarios, se presenta canícula durante julio, sin embargo, ésta se presenta de manera más acusada que en los anteriores.

Análisis espacial

Fuente: Elaboración propia.

Las condiciones ambientales que se presentan propicias para la distribución potencial bajo las condiciones actuales son amplias, principalmente en las sierras circundantes a la Cuenca de México, Valle de Toluca y Valle de Puebla-Tlaxcala. Las distintas categorías se presentan a manera de anillos concéntrico conforme las condiciones del relieve.

En la línea base (Figura 3 (A)) la franja que se clasifican como muy adecuado y adecuado se presentan a partir de los 2800 msnm hasta los 3800 msnm en sitios como la RBB Mariposa Monarca, Nevado de Toluca; en la Sierra de las Cruces y Ajusco-Chichinautzin se ubica alrededor de los 2900 msnm hasta los 3800 msnm aproximadamente; en tanto que en la Sierra Nevada a partir de los 3000 y hasta los 3800 msnm; para el volcán La Malinche va de 3000 a 3600 msnm; en el área del Pico de Orizaba y Cofre de Perote desde los 3000 hasta los 3800 msnm; en la zona de Mineral del Chico inicia alrededor de 2800 y llega a 3000 msnm; en tanto que en los límites de Tlaxcala y Puebla se aprecia a partir de 2800 y por encima de 3200 msnm. Fuera de esta franja altitudinal la aptitud de albergar al oyamel disminuye.

El escenario REA 4.5 (Figura 3 (B)) muestra una importante reducción en las áreas clasificadas como muy adecuado y adecuado y un incremento sustancial en las clasificadas como poco adecuado.

El área de la RBMM pasa de predominantemente muy adecuado y adecuado a poco adecuado; por su parte se observa la fragmentación del bosque en el Nevado de Toluca, Sierra de las Cruces y Ajusco-Chichinautzin; en la Sierra Nevada destaca el ascenso de la clase moderadamente adecuado a las cercanías de las cimas de los volcanes Iztaccíhuatl y Popocatépetl; en el volcán La Malinche se advierte una cierta condición estable; el Pico de Orizaba pasa a una condición poco adecuada mientras que en el Cofre de Perote se mantiene, sin embargo el corredor que unía a estos se ve deteriorado al clasificarse como poco adecuado; en los límites de los estados de Puebla y Tlaxcala las condiciones pasan a ser poco adecuadas, a pesar de ello se mantiene una porción en adecuado en la zona circundante a la Presa Iturbide; al norte, el Parque Nacional El Chico predominan las condiciones moderadamente adecuadas, en detrimento tanto de la muy adecuada como la adecuada.

En el escenario REA 8.5 (Figura 3(C)) las condiciones adecuadas se restringen solo a altitudes superiores a 3000 msnm, el otro sitio donde se mantienen condiciones adecuadas es en el Cofre de Perote en las laderas a barlovento; mientras que las muy adecuadas de 3400 hasta 3800 msnm pero solo en las laderas del Iztaccíhuatl y Popocatépetl; las moderadamente adecuadas se localizan entorno al cono del volcán La Malinche, debajo de los 3000 metros de la Sierra Nevada y por encima de los 3800, en el Cofre de Perote en las laderas a sotavento.

El resto de las áreas que anteriormente se encontraban en categorías de moderadamente adecuado, adecuado y muy adecuado pasan a poco adecuado o se integran a la clase de no adecuado.

De las variables climatológicas utilizadas para la construcción del modelo, se encontró que la precipitación fue la que aportó mayor cantidad, en tanto las unidades fototérmicas tuvieron nula aportación al mismo; cabe resaltar que se excluye la componente del relieve por no ser de interés para la investigación. La tabla 3 muestra el porcentaje relativo de aporte al modelo de cada elemento climático.

A partir del análisis espacial se pudo obtener el área aproximada de distribución de cada clase para el oyamel, así como la variación de la superficie con respecto a la línea base (Tabla 4).

Conclusiones

Los bosques de oyamel son importantes ecosistemas que sirven de refugio de cientos de especies de animales, plantas y hongos, algunas de ellas tan representativas como la mariposa monarca; es de los biomas con mayor capacidad de captura de carbono, almacena e infiltra grandes cantidades de agua; provee diversas materias primas como madera usada en ebanistería, leña y resinas de uso industrial, así como ornamental como árboles de navidad.

La presión directa de las actividades humanas como la extracción de la madera de forma clandestina con fines comerciales, el cambio de uso de suelo para abrir nuevos terrenos agrícolas son dos de los principales motivos por los que la superficie ocupada por el oyamel se ha visto reducida, a esto se le suma la del cambio climático.

Las actuales condiciones climáticas lo colocan como sitios con características de alta humedad, escasos periodos de estiaje y con oscilaciones térmicas bajas a lo largo del año.

En términos generales, tanto el REA 4.5 y 8.5 presentan aumentos significativos de la temperatura y la disminución de las precipitaciones acumuladas. Se observa que la amplitud térmica tiende a hacerse más acusada, en tanto que, el periodo de estiaje muestra signos de volverse más aguda.

El corrimiento de la estación lluviosa, así como de la distribución de la cantidad de lluvia puede generar la alteración del ciclo fenológico del oyamel. La sequía intra estival presenta valles más bajos que en la actualidad, lo que generaría estrés adicional sobre el ecosistema.

Solo el modelo con forzamiento radiativo 8.5 muestra incrementos en la precipitación, coincidente en el mes de octubre para las tres estaciones analizadas.

Las variaciones en temperatura y precipitación como factores tensionantes representarían el aumento de la vulnerabilidad del oyamel al ataque de plagas y enfermedades, así como el potencial de incremento de ocurrencia de incendios forestales. Los modelados de ambas proyecciones muestran la tendencia de migración en la vertical, sin embargo, en altitudes que resulten insuficientes para cubrir sus necesidades puede darse el caso del eventual desplazamiento por otras especies mejor adaptadas.

A pesar de que actualmente la especie no se encuentra enlistada bajo la protección de la NOM-059-SEMARNAT-2010, muestra una tendencia generalizada a reducir su población. La progresiva pérdida de cobertura boscosa debido al cambio climático pondría en peligro de extinción no solo al oyamel sino a las especies que conforman el ecosistema que a su vez significaría el colapso de las cadenas tróficas.

Frente a estos escenarios es imperativo buscar alternativas que busquen proteger al bosque de oyamel y con ello al ecosistema completo mediante la elaboración de planes y programas de manejo especialmente dirigidos y que permitan la participación en la gestión por parte de las comunidades forestales, en su mayoría de origen indígena, en la protección y aprovechamiento sustentable de los recursos. Actualmente se desarrollan distintos proyectos de conservación, ejemplo de ello es la Universidad Michoacana de San Nicolás Hidalgo cuya línea de investigación se enfoca en la migración asistida, la creación de bancos de germoplasma y la instalación de viveros para la producción de árboles con fines de reforestación.

Es recomendable el fortalecimiento y consolidación de la red de estaciones climatológicas, así como el registro y sistematización de datos, a fin de lograr investigaciones con mayores niveles de confiabilidad, y que sirvan para elaborar estudios detallados con la finalidad de servir como una base sólida en la que se sustenten políticas públicas de planeación y conservación a largo plazo.

Referencias

Fechas de Publicación

Histórico