Our objectives were to predict and map the climatic niche for Pinus leiophylla for a period of normalization (years 1961-1990) and future (2030, 2060 and 2090) climates, and to suggest management strategies to accommodate climate changes, and discuss implications for conservation. A bioclimate model predicting the presence or absence of P. leiophylla (lumped with its putative variety P. leiophylla var. chihuahuana ) was developed by using the Random Forests classification tree on Mexican and Unites States of America forest inventory data. The bioclimatic model had an average error of prediction of 4.6 %. The model used six predictor variables, dominated by precipitation variables. Projecting the 1961-1990 climate niche into future climates provided by three general circulation models and two greenhouse-effect gas emission scenarios, suggested that the area occupied by the niche should diminish rapidly over the course of the century: a decrease of 35 % by the decade surrounding 2030, 50 % for 2060, and 76 % for 2090. The most serious habitat reduction occurs at both latitudinal extremes of the species distribution: Chiricagua Mountains, Arizona, Unites States of America in the northern extreme, and at Oaxaca State, Mexico, in the southernmost extreme. There is no indication at all of expansion of suitable climatic habitat northwards. We urge establishing seed banks encompassing seed from provenances sampled from the largest part possible of the natural distribution, and start assisted migration tests, to realign the natural populations with the climate for which they are adapted and that will occur at higher altitudes.

Nuestros objetivos fueron predecir y mapear el nicho climático para Pinus leiophylla, para los climas del período 1961-1990 y del futuro (2030, 2060 y 2090), sugerir estrategias de manejo para adaptarse al cambio climático y discutir implicaciones para su conservación. Se desarrolló un modelo bioclimático que predice la presencia o ausencia de P. leiophylla (agrupado con su supuesta variedad P. leiophylla var. chihuahuana) utilizando la técnica de árboles de clasificación Random Forests, con datos del inventario forestal de México y Estados Unidos. El modelo bioclimático tuvo un error promedio de predicción de 4.6 %. El modelo utilizó seis variables de predicción, dominadas por variables de precipitación. La proyección del nicho climático 1961-1990 en climas futuros, a partir de tres modelos de circulación general y dos escenarios de emisiones de gases de efecto invernadero, sugieren que el área ocupada por el nicho disminuirá rápidamente durante el siglo: una disminución del 35 % en la década alrededor del 2030, del 50 % para 2060, y 76 % para 2090. La reducción más grave de hábitat se produce en los dos extremos latitudinales de la distribución de las especies: Montañas Chiricagua, Arizona, en Estados Unidos en el extremo norte y en el estado de Oaxaca, México, en el extremo sur; no hay indicio alguno de la expansión del hábitat climático propicio hacia el norte. Se sugiere el urgente establecimiento de bancos de semillas de procedencias colectadas en la mayor parte posible de su distribución natural, e iniciar ensayos de migración asistida, para realinear las poblaciones al clima para el cual están adaptadas, mismo que ocurrirá a mayores altitudes.

Understanding the patterns of morphological variation along environmental gradients is useful to distinguish taxonomical entities and take management decisions. Morphological variation of cones, seeds, fascicles and buds were analyzed on individuals from 16 natural Pinus devoniana (ex P. michoacana) populations collected along an altitudinal transect (from 1,600 to 2,450 m, one population every 50 m of altitudinal difference) near Morelia city, in the state of Michoacán (Mexico). Univariate (ANOVA) and multivariate (principal components and cluster) statistical analyses were conducted, in order to clarify if putative varieties exist. There is a significant pattern of morphological altitudinal variation, where populations from lower altitudes have larger cones, needles and seeds than populations from middle and higher altitudes. The multivariate analysis discriminated individuals and populations in two clearly distinctive groups: one group made of individuals from populations originated between 1,600 m to 1,850 m of altitude, which would correspond to Pinus devoniana var. devoniana (equivalent to P. michoacana typical); a second group is composed of individuals from populations originated between 1,950 m to 2,450 m, which would correspond to Pinus devoniana var. cornuta (ex P. michoacana var. cornuta). Thus, we suggest the separation of P. devoniana in at least two varieties: P. devoniana var. devoniana and P. devoniana var. cornuta.

La comprensión de los patrones de variación morfológica a lo largo de gradientes altitudinales es útil para distinguir entidades taxonómicas y tomar decisiones de manejo. Conos, semillas, fascículos y yemas fueron analizados en individuos de 16 poblaciones naturales de Pinus devoniana (ex P. michoacana). El material se recolectó a lo largo de un transecto de 1,600 a 2,450 msnm, a razón de una población cada 50 m, cerca de la ciudad de Morelia, en el estado de Michoacán. Se realizaron análisis univariados (ANOVA) y multivariados (componentes principales y de conglomerados) para determinar las variedades putativas presentes. Los resultados evidenciaron la existencia de un patrón altitudinal significativo de variación morfológica, en el cual las poblaciones de localidades bajas tienen hojas, semillas y conos más largos que las desarrolladas en lugares con altitudes medias y altas. El análisis multivariado discriminó dos grupos evidentemente distintivos: uno compuesto por individuos de poblaciones originadas entre 1,600 y 1,850 msnm: Pinus devoniana var. devoniana (equivalente a P. michoacana típica), y otro en el que se integraron ejemplares de poblaciones procedentes de sitios ubicados en un intervalo de 1,950 a 2,450 msnm: Pinus devoniana var. cornuta (ex P. michoacana var. cornuta). Así, se sugiere la separación de las poblaciones estudiadas de P. devoniana en al menos dos variedades: P. devoniana var. devoniana y P. devoniana var. cornuta.