Introducción: La complejidad del ecotono del límite superior del bosque de Pinus hartwegii puede funcionar como una barrera para la migración altitudinal.Objetivos: La estructura del bosque de P. hartwegii se caracterizó a través del ecotono de su límite superior en el Nevado de Toluca; además, se determinaron los posibles cambios en los límites altitudinales en los últimos 25 años.Materiales y métodos: Se determinaron la densidad arbórea, altura y diámetro normal de cada árbol, y la regeneración natural del bosque. Los límites altitudinales se determinaron mediante el método de la ventana móvil en imágenes de satélite de 1989 y 2014.Resultados y discusión: La densidad del arbolado se redujo entre 8 y 70 % en el gradiente altitudinal. Los límites altitudinales (3,980-4,130 m) no presentaron diferencias significativas (P = 0.07) entre los años evaluados. El ecotono es de forma difusa, regulado por limitaciones al crecimiento. Esto implica que el bosque modificará su distribución altitudinal en función de los cambios ambientales, de acuerdo con su amplitud ecológica y la rapidez con que puedan establecerse nuevos individuos a mayores altitudes.Conclusión: Es importante establecer monitoreos, a largo plazo, para saber si la especie es realmente capaz de migrar altitudinalmente como resultado del incremento de temperatura ambiental.

Geographical distribution of forest species is closely regulated by environmental conditions, particularly temperature and precipitation. Climate change predicted by general circulation models is expected to modify the distribution of many species’ distribution, especially those adapted to extreme environmental conditions, leading to large-scale migrations or local extinctions. The aim of this research was to determine the potential impact of climatic change on Pinus hartwegii geographic distribution and the niche breadth of its populations. Ecological niche models were used by generated with four different algorithms based on 19 bioclimatic variables in addition to altitude. Climatic niche breadth was delimited by the dispersion of species occurrence records within the intervals of the bioclimatic variables. We modelled future distribution based on three general circulation models, MIROC-ESM-CHEM, CCSM4 and HadGEM2-ES, using two representative concentration pathways (RCP) 2.6 and 8.5, for two-time horizons 2050 and 2070. Niche breadth analysis showed narrow ranges of suitability, indicating a strong relationship between the presence of P. hartwegii with the temperature of the warmest quarter and precipitation of the coldest quarter. In addition, the suitability area of P. hartwegii is predicted to be reduced up to 70% by 2070; the populations of the extreme northern and southern latitudes will be reduced in greater proportion than those of central Mexico. This suggest that environmental suitability area of P. hartwegii are reduced by the effect of the increase in environmental temperature. Therefore, it is necessary to monitor extreme populations of this species in the long term in order to establish efficient conservation strategies and well adaptive management facing climate change.