Reciprocal Common Garden Altitudinal Transplants Reveal Potential Negative Impacts of Climate Change on Abies religiosa Populations in the Monarch Butterﬂy Biosphere Reserve Overwintering Sites

: Research Highlights: Reciprocal altitudinal transplants of Abies religiosa seedlings within the Monarch Butterﬂy Biosphere Reserve (MBBR) allow prediction of the impacts of climatic change, because they grow in sites with a climate that differs from that of their origin. Background and Objectives: Climatic change is generating a mismatch between the sites currently occupied by forest populations and the climate to which they have adapted. This study determined the effect on the survival and growth of A. religiosa seedlings of transfer to sites that were warmer or colder than that of the origin of their seeds. Materials and Methods: Eleven provenances of A. religiosa , collected along an altitudinal gradient (3000 to 3550 m a.s.l.), were assayed in common gardens in three sites of contrasting altitude: 3400, 3000 and 2600 m a.s.l. The results were evaluated by ﬁtting a response curve with a mixed model. Results: The climate transfer distance for the seasonal balance between the temperature conducive to growth (degree days above 5 ◦ C) and the available precipitation (a ratio expressed as dryness index) dominated the shape of the response function curve. The rainy season (June–October) dryness index transfer distance was critical for survival, while that of the cold and dry season (November–February) was critical for aerial biomass, and the annual index was critical for the increase in basal diameter. The effect of climatic transfer distance is much more negative (triggering about 45% mortality) when transfer is toward warmer and dryer sites (at 400 m lower in altitude, +1.9 ◦ C warmer and 16% less precipitation), than when shifting toward colder and wetter sites (400 m higher in altitude, resulting in 95% survival). Conclusions: The projected higher temperatures and lower precipitation due to climatic change will undoubtedly cause severe mortality in young A. religiosa seedlings. A 400 m shift upwards in altitude to compensate for climatic change (assisted migration) appears to be a feasible management action.


Introduction
In the near future, climatic change will generate warmer climates and displace the climatic zones towards the poles and towards higher altitudes in mountainous regions [1].
Given that, as the result of a long process of evolutionary adaptation, species and their populations are distributed in localities that are generally coupled to a climate suitable for them, climatic change will generate a decoupling between the populations and the climate to which they are adapted [2,3], influencing phenological changes, population and growth dynamics [4] and reproduction [5]. Stress induced by this decoupling favors the frequency and intensity of the development of pests and diseases [6][7][8][9][10] and defoliation through drought stress [11], inducing a severe process of forest decline, which is already underway globally [12], especially at the xeric limit (lower altitudinal or south in the northern hemisphere) [13,14].
Forest tree species that present a wide geographic distribution can present an important genetic differentiation among populations in terms of quantitative characters, such as adaptation to the diversity of environmental conditions [3,15,16]. This differentiation is expressed in different variables such as growth rate, seasonality of elongation of the apical bud and resistance to frost damage [17][18][19][20]. Provenance assays allow the expression of these differences by cultivating plants from different geographic origins (provenances) under the same environmental conditions [21]. In general, for forest populations, a pattern of variation of growth has been reported in which the provenances originating from the colder (higher altitude) part of the natural distribution express a lower potential for growth in height, but a greater resistance to frost damage while, in contrast, those from the low part of the altitudinal distribution (warmer) have greater potential for growth but are more susceptible to frost damage [22][23][24].
In Mexico, one of the most important species for biological conservation is Abies religiosa (sacred fir or oyamel), because it hosts the overwintering sites of the monarch butterfly (Danaus plexippus L.) for five months (November-March) [25]. The areas of the greatest extension of A. religiosa are found within the Trans-Mexican Volcanic Belt, at between 2400 and 3600 m a.s.l. and mainly between 19 • and 20 • in latitude North [26]. The Monarch Butterfly Biosphere Reserve (MBBR) is found within the Trans-Mexican Volcanic Belt, on the border between the states of Michoacán and Estado de México. The forests of A. religiosa can be strongly affected by climatic change-it will induce a decoupling between the populations of A. religiosa and the environment suitable for them, such that by the decade 2090, the models project that suitable climatic habitat for A. religiosa will have disappeared completely from the MBBR [26]. Previous studies have reported that A. religiosa presents a high percentage of non-viable seeds (63-79%) [27,28], in addition to very low recruitment and high mortality of sacred fir seedlings, necessitating the implementation of an active restoration program or other management actions [29]. Additionally, so that the populations of this species can remain coupled to their favorable climatic habitat, their upwards altitudinal displacement of between 300 and 500 m will be necessary by 2060 [30], due to the projected increases in temperature and diminution in precipitation. This could have a considerable effect on the overwintering of the monarch butterfly, in addition to severe environmental impacts. For these reasons, human-assisted displacement of seed sources towards higher altitudes (through reforestation programs) appears to be an option for consideration as a strategy of adaptive management to climatic change. This is because, while the populations of forest species can, and in fact currently are, migrating altitudinally towards higher altitudes through natural means, the speed at which they can accomplish this is lower than that necessary to remain coupled to the climate that favors them [31,32].
The results of two previous field tests indicate that it is feasible to conduct assisted altitudinal migration to above the provenances of A. religiosa, translocating the seed sources up to 400 m higher in altitude than their origin, but only when planted beneath the shade protection of local shrubs (such as Baccharis conferta) serving as nurse plants [33]. However, to date, no experiments have determined what would happen if young A. religiosa seedlings are exposed to greater temperatures and lower precipitation, as projected in scenarios of climatic change [26,34]. There is also evidence that there are serious difficulties in terms of the natural regeneration of young seedlings, because the warm and dry season of March-May is now more dry and hot than usual, causing a severe mortality of seedlings in April [29]. Although an important genetic differentiation has been demonstrated among populations of A. religiosa for resistance to frost damage along altitudinal gradients [20], a significant genetic differentiation has not been demonstrated among populations for resistance to drought, for example, in the differential accumulation of solutes, such as the amino acid proline, that confer physiological resistance to drought [35,36].
Considering the above, the objective of this study was to evaluate the effects of climatic transfer (climatic difference between the site of origin and the site of plantation) on the survival and growth of young A. religiosa seedlings, through reciprocal transplants of eleven provenances of A. religiosa to common gardens at three contrasting altitudes. The climatic transfer was conducted by simulating three possible scenarios: (i) towards altitudes lower than the seed source origin (simulating a future climate of higher temperature and lower precipitation); (ii) towards higher altitudes (to confirm the viability of assisted migration); and (iii) to a site of intermediate altitude as a reference. The latter two sites were located in the core zone and on the edge of the buffer zone of the MBBR, respectively, in Estado de México, in central-western Mexico. The intention was to generate guidelines for decisionmaking in terms of the transfer of seed sources in a program of assisted migration, as well as to quantify the impacts of climatic change and the risks of failing to practice adaptive management in the face of climatic change.

Study Site
The study was conducted within the Monarch Butterfly Biosphere Reserve (MBBR), which forms part of the Trans-Mexican Volcanic Belt, and in Tlalpujahua de Rayón, Michoacán. The altitudes of the MBBR range from 2200 to 3640 m a.s.l. and the annual average temperatures range from 8 to 22 • C [37]. The vegetation is mainly comprised by forests of fir (A. religiosa), pine (Pinus spp.), oak (Quercus spp.), and cedar (Cupressus spp.) [38].

Seed Collection Sites
Cones were collected (in December 2017) from eleven provenances of A. religiosa, along an altitudinal transect (3000 to 3550 m a.s.l. in altitude; altitudinal difference of 50 m between sites), within the MBBR (Table 1). At each site, cones were collected at random from ten trees, at a distance of at least 30 m apart [39], in order to reduce the possibility of endogamy.

Provenance Test Sites in Common Gardens
Three field provenance tests were established in common gardens (Table 3): (a) high altitude site at 3400 m a.s.l., to confirm the viability of altitudinally upwards assisted migration, within the core zone of the MBBR (at an approximate linear distance of 300 m from a monarch butterfly overwintering site); (b) intermediate altitude site at 3000 m a.s.l., on the edge of the buffer zone of the MBBR; and (c) low altitude site at 2600 m a.s.l., to quantify the effects of a warmer, drier climate than that of the seed origin (considering that the temperature changes by approximately 0.5 • C for each 100 m of altitudinal difference [34]). Table 3. Location and some climatic variables (their acronyms as in Table 2) of sites of tests of provenances in common gardens.

Production of Plants in Nursery
Once the seed was obtained from the mother tree, it was stratified at 4 • C for 14 days [39,41]. Each lot of seeds subjected to stratification comprised a similar quantity of seeds originating from each mother tree, so that the progenies were represented equally in the provenances.  soil, agrolite, vermiculite and peat moss (volume 60:15:15:10, respectively). Transplantation was conducted with an equal quantity of seeds germinated from each mother tree, so that they were represented equally in the provenances. To prevent differential growth among the provenances in the nursery as a result of possible environmental differences within the shade house, the plants were positioned in a completely random block design (6 blocks, 11 provenances, 54 plants per plot).

Test of Provenances in Common Gardens in the Field
In each field site (as in Table 3), three raised beds were constructed, with wooden frames of 5 × 1.5 × 0.4 m in length:width:height. The internal base of each raised bed was covered by a metal mesh (aperture size 3.5 cm) to prevent access by gophers. To ensure that the sites differed only in terms of the climate as a result of their altitude, and thus to evaluate the impact of the climatic transfer in each provenance, the raised beds were filled with the same forest soil, transported from a forest site dominated by A. religiosa at the intermediate altitude (3000 m a.s.l.), thus avoiding any confounding effect of the local soil. Finally, the raised beds were covered by a shade mesh (35%) at a height of 2.5 m, to simulate the protective effect of nurse plants [33] (Figure 1). originating from each mother tree, so that the progenies were represented equally in the provenances.
Plant production was conducted for 14 months (May 2018-June 2019) in a shadehouse (35% shade mesh) in the Instituto de Investigaciones sobre los Recursos Naturales (INIRENA), in Morelia city in Michoacán, Mexico. Plants were germinated in a seedbed and later transplanted into containers of 380 cm 3 , with a mixed substrate comprising loamy soil, agrolite, vermiculite and peat moss (volume 60:15:15:10, respectively). Transplantation was conducted with an equal quantity of seeds germinated from each mother tree, so that they were represented equally in the provenances. To prevent differential growth among the provenances in the nursery as a result of possible environmental differences within the shade house, the plants were positioned in a completely random block design (6 blocks, 11 provenances, 54 plants per plot).

Test of Provenances in Common Gardens in the Field
In each field site (as in Table 3), three raised beds were constructed, with wooden frames of 5 × 1.5 × 0.4 m in length:width:height. The internal base of each raised bed was covered by a metal mesh (aperture size 3.5 cm) to prevent access by gophers. To ensure that the sites differed only in terms of the climate as a result of their altitude, and thus to evaluate the impact of the climatic transfer in each provenance, the raised beds were filled with the same forest soil, transported from a forest site dominated by A. religiosa at the intermediate altitude (3000 m a.s.l.), thus avoiding any confounding effect of the local soil. Finally, the raised beds were covered by a shade mesh (35%) at a height of 2.5 m, to simulate the protective effect of nurse plants [33] (Figure 1). The seedlings produced in the nurseries were transplanted to the raised beds in July 2019 (beginning of the rainy season), at age 14 months. In each site, six completely random blocks were established (11 provenances per block, five seedlings per provenance per plot in line, at a spacing of 0.2 × 0.2 m and with two blocks per raised bed. At the extremes of the raised bed, a protection strip was planted, comprising a mixture of seedlings of all of the provenances, in order to control the edge effect ( Figure 1).
The trees (the seeds of which were collected at different altitudes) represented in the common garden tests by their seedlings are herein termed populations, while the location of origin of a population is called the provenance.

Measurements
For sixteen months, monthly or bimonthly evaluations of plant survival and height were conducted, with semestral evaluations of the basal diameter. The first evaluations of The seedlings produced in the nurseries were transplanted to the raised beds in July 2019 (beginning of the rainy season), at age 14 months. In each site, six completely random blocks were established (11 provenances per block, five seedlings per provenance per plot in line, at a spacing of 0.2 × 0.2 m and with two blocks per raised bed. At the extremes of the raised bed, a protection strip was planted, comprising a mixture of seedlings of all of the provenances, in order to control the edge effect ( Figure 1).
The trees (the seeds of which were collected at different altitudes) represented in the common garden tests by their seedlings are herein termed populations, while the location of origin of a population is called the provenance.

Measurements
For sixteen months, monthly or bimonthly evaluations of plant survival and height were conducted, with semestral evaluations of the basal diameter. The first evaluations of plant height and basal diameter were conducted six months after plantation, in the expectation that the plants would be well established by that time. The final measurement was taken at the end of the growing season of the following year (November 2020). Plant height was measured with a 1 mm precision up to the tip of the apical bud, supporting the ruler on a bar of aluminum, which was in turn supported on the edge of the wooden frame (to avoid measurement errors due to the irregular surface of the substrate). Basal diameter was measured at the base of the stem with a digital vernier (KNOVA ® , Cupertino, CA, USA) at a 0.1 mm precision. The increase in seedling height and basal diameter was obtained by subtracting the first measurement (December 2019) from the final measurement (November 2020). Survival was estimated as the percentage of live plants per plot. Substrate moisture content was measured with a W.E.T (water content, electrical conductivity and temperature) sensor of type W.E.T-2 (Delta-T-Devices, Cambridge, UK) as well as gravimetrically, by measuring the wet and dry weights of soil samples.

Estimation of Aerial Biomass
Following the methodology of Martinez-Luna, et al. [42], the aerial biomass of the seedlings in the common garden field assay was estimated using a regression formula from the measurements of total height and basal diameter. The formula was fitted by measuring (through destructive harvest) seedling height, basal diameter and dry weight of the stem, branches, and foliage (dried for 48 h at 65 • C in a Shel Lab drying oven (model FX4) and weighed on a Mettler Toledo analytical balance, model AB204-S). Five seedlings of diverse sizes were selected for this harvest from each of the 11 assayed provenances (harvested from a surplus lot of plants). The regression was fitted using Proc Reg of SAS [43,44], from which the following parameters were obtained: where Ba is the aerial biomass (g), Ht is the total height (mm), and Db is the basal diameter (mm 0.1).

Climatic Data
Climatic data from the provenance of origin of the seeds were obtained for the reference period 1961-1990 through climatic spline models (based on [34]), available at the website [40]. The climate of the provenance is considered the recent historic climate to which the provenances have evolved and adapted [45]. The climate of the field assay sites, which was that experienced by the seedlings while growing in the field common gardens, was estimated using three data loggers HOBOs (temperature/external channel data logger; Onset Computer Corporation ® , Bourne, MA, USA) per site, set to record temperature every 30 min. Monthly precipitation was measured by capturing rainwater, averaging the results from two rainwater traps per site.
The climatic variables estimated for both the provenances and the sites are described in Table 2, and are considered of physiological importance for representing the seasonal or annual balance between the available humidity and the temperature in terms of seedling growth [17,46]. The seasonal indices warm and dry season dryness index (WDSDI), rainy season dryness index (RSDI) and cold and dry season dryness index (CDSDI) aim to represent the seasonality of each of the A. religiosa study sites and are based on the seasonality observed in [33].
From the climate of the provenances and of the field test sites, the climatic or altitudinal transfer distance (CTD or ATD = climate or altitude of the test site − climate or altitude of the source of the provenance) was estimated in order to represent the impact of growing in a site with a climate (or altitude) that differed from that of the seed origin.

Statistical Analysis
A curve was estimated for each of the four response variables studied (survival, increase in basal diameter, increase in height, and aerial biomass of A. religiosa seedlings), to measure the effect of the climatic or altitudinal transfer (growing in a site with a climate or altitude that differed from that of the seed origin), using a mixed model in Proc Mixed of SAS [43,44], based on refs [16,46].
The fixed effects considered three components of the effect of climate as a selective force that produces a genetic differentiation among populations and/or a response of phenotypic plasticity (also considered as a genotype × environment interaction [15]): (a) the climate of the provenance (C); (b) the climatic transfer distance (D = climate of the test site − climate of the provenance); and (c) the interaction between climate of the provenance and climatic transfer distance (C × D). Random effects provided the variance that was not explained by the fixed effects; i.e., those derived from the experimental design: site, provenance, block and the interaction between the site and the provenance. The climatic transfer distance is a concept originally known as "ecological distance" [47]; more details regarding the reasoning of this model can be found in refs [16,48].

Mixed Model
The following mixed model was fitted, based on refs [16,46]: where Y ijkl is the survival or aerial biomass or increase in basal diameter or increase in plant height, corresponding to the lth tree for the jth provenance in the kth block in the ith test site, and β 0 is the intercept. D ij is the climatic transfer distance (difference between the value of a climatic variable in the test site and that of the provenance) for the jth provenance in the ith test site, C j is the value of the climatic variable in the jth provenance, and D ij × C j is the interaction between the climatic transfer distance for the jth provenance in the ith test site and the climatic variable in the jth provenance. S i is the effect of the ith test site, P j is the effect of the jth provenance, B k (S i ) is the effect of the kth block nested within the ith test site, and e ijkl is the error term.

Selection of Variables to Fit the Best Mixed Model
The climatic variables for use in the mixed model were selected using the following steps, independently for each response variable, based on refs [16,46]: (1) Five climatic variables were selected that best described the climate of the provenance, estimating the Spearman correlations between the values of the climatic variables and the average value per provenance across sites of the response variables, and selecting those with the highest |r| value. (2) The climatic variables for the climatic transfer distance were selected by fitting a reduced mixed model, eliminating from the model (Equation (2)) the term of climate of the provenance (C j ) and its respective interaction (D ij × C j ). The five climatic variables selected were those for which the model obtained the lowest (and thus the best) value of the Akaike information criterion (AIC) and, in addition, that necessarily presented the estimated parameter of the quadratic term both negative and significant, in order to ensure that it was biologically sound [15,16]. (3) Subsequently, 5 × 5 = 25 full "competing" models were run, which included all of the possible combinations of the five variables of the climate of provenance and the five variables of climatic transfer distance preselected in the previous two steps. The best model was selected based on the AIC value.
In total, 100 models (25 "competing" models for each of the four response variables) were run, in order to determine the best model for understanding the response of survival, increase in basal diameter, increase in plant height and aerial biomass.
Finally, from the regression parameters estimated for the fixed effects, quadratic curves of response to the climatic transfer were estimated for each response variable, in order to visualize the effect on survival and growth of growing in a site with a climate that differs to that of the origin of the populations.

Climatic Variables That Best Explained the Climatic Transfer Distance Response Function
The survival response function was best explained (lowest AIC value = 1651.1) by the quadratic term of the rainy season dryness index transfer distance (RSDI_td 2 ; p = 0.0112), although the linear term was not significant (p = 0.3542). The other fixed effect terms, the climate of the provenance term selected (warm and dry season dryness index, WDSDI; p = 0.0941) and the interaction between climatic transfer distance and the climate of the provenance (RSDI_td × WDSDI; p = 0.1273), were also not significant. Regarding the random effect terms, only the block nested within the site was significant (p = 0.0252), contributing 21% of the total variance, where 100% is the sum of the contribution of the random terms only (Table 4). Table 4. Analysis of the mixed model for seedling survival, biomass and increase in basal diameter. Akaike information criterion (AIC), estimated parameters (coefficients of regression for fixed effects and variances estimated for random effects), contribution to the total variance (of the random effects) and significance (p). The biomass response function was best explained (lowest AIC value = 2397.2) by the significant linear and quadratic term of the cold and dry season dryness index transfer distance (CDSDI_td, p = 0.0002; CDSDI_td 2 , p < 0.0001; respectively). The climate of the provenance term selected (warm and dry season dryness index, WDSDI) was also significant (p = 0.0198), as was the interaction CDSDI_td × CDSDI (p = 0.0395). The random effect term that contributed most to the total variance (among the random effects) was block nested within site (p = 0.0196), contributing 7% of the total variance ( Table 4).

Parameter or Source of Variation
The increase in basal diameter response function was best explained (lowest AIC value = 1917.3) by the significant linear and quadratic term of the annual dryness index (ADI) transfer distance (p < 0.0001). The climate of the provenance term selected (mean temperature of the coldest month, MTCM) was also significant (p = 0.0036), as was the interaction ADI_td × MTCM (p < 0.0001). The random effect term that most contributed to the total variance was block nested within site (p = 0.0095), contributing 15% of the total variance ( Table 4).
The increase in seedling height response function was best explained (lowest AIC value = 5122.4) by RSDI, although, as with transfer distance, it was not significant (linear: p = 0.23; quadratic: p = 0.3517), as well as the climate of the provenance term selected (MMIN, p = 0.1605). Again, only block nested within site was significant (p = 0.0096) among the random terms, contributing 11% of the total variance (results not shown in Table 4 for simplicity). In all of the four response variables analyzed, site was not significant, and the error term contributed the most among all the random terms (from 79% for survival to 88% for aerial biomass; Table 4).

Curves of Response to the Climatic and Altitudinal Transfer Distance
The predicted response curve for survival indicated that when the seed sources were moved toward warmer and dryer sites (toward lower altitude, right side of the curve in Figure 2a), there is a severe induced mortality. In sharp contrast, when the shift is conducted toward a colder and moister site (to higher altitudes, sites at 3400 m a.s.l.), survival is maintained mostly above 95%. This is even more evident when fitting a simple quadratic regression of survival against altitudinal transfer distance (Figure 3a), where exceeding the 400 m shift to lower altitudes (grossly equivalent to +2 • C, based on a lapse rate of 0.5 • C per 100 m [34]) causes significantly increased mortality. It should also be noted in Figure 2a that the maximum value of the response curve is displaced slightly towards the left side (colder and moister sites), either for the average species curve or for the extreme lower altitudinal population curve. This suggests that the optimum of survival for the species in general, and for the extreme lower population, would be a colder site than that occupied at present.    The predicted response curve for aerial biomass indicated that, when seed sources are moved toward warmer and dryer sites (to a lower altitude, right side of curve the in Figure 2b), there is little negative effect on the growth expressed as aerial biomass: the observed values fluctuated around the maximum predicted value of the curve (about 4.2 g of aerial biomass in dry weight). In contrast, when seed sources are shifted toward colder and moister sites, there is an important loss of accumulated aerial biomass. This becomes more evident when plotting the average biomass per population and per site against the altitudinal transfer distance (Figure 3b).
The predicted response curve for the increase in basal diameter indicated that the greatest growth occurs at a site that is warmer and dryer than that occupied by the populations at present, given that the curves are displaced to the right (Figure 2c). It is also clear that when they are moved toward colder and moister sites (left of the x-axis) or warmer and dryer sites (right side of the x-axis), relative to the optimum, they present a lower growth in basal diameter. The resulting curve regarding altitudinal transfer distance indicated clearly that the optimum growth, for all the provenances in general, was achieved at the site of intermediate elevation (3000 m a.s.l.; Figure 3c).
When comparing the average value across sites for the mean annual temperature and annual precipitation to the average value across sites of the reference period 1961-1990, it is evident that the sites have already increased in temperature (for the observed period) by an average of 0.58 • C, while precipitation has decreased by 27 mm (around 3%) ( Table 5). When plotting seedling mortality, monthly precipitation and mean temperature per site, and comparing among the sites (Figure 4), it is evident that precipitation differences among sites are not particularly large, but seem to have a meaningful effect, nevertheless. The low altitude site (2400 m a.s.l., 1407.4 mm total precipitation in 15 months; Figure 4a), had 16.5% lower total precipitation over the observed period than the site of intermediate altitude (3000 m a.s.l., 1685.14 mm; Figure 4b), while the site of elevated altitude (3400 m a.s.l., 1853.6 mm; Figure 4c) had 9.1% higher total precipitation than the intermediate site.
The temporal patterning of the dry/rainy season seems to be similar across the three sites. Surprisingly, the apparently relative minor differences in mean temperature between the sites at 2600 m a.s.l. (14.9 • C averaging across all the 15 months observed) and 3000 m a.s.l. (13.02 • C; 1.9 • C of difference among those two sites), combined with the 16.5% difference in precipitation, seemed to be enough to trigger much higher mortality at the site at 2600 m a.s.l., particularly during the dry months, whether cold and dry (December-February) or warm and dry (March-May; Figure 4a). Moreover, an unexpected and relatively high mortality occurred just after the rainy season abruptly ended in October 2020 (Figure 4a), suggesting that several seedlings were already debilitated before the rainy season, perhaps subsequently expressing a drought stress legacy effect after the rainy season. In sharp contrast, seedling mortality at the site at 3400 m a.s.l., which was much colder and only 9.1% more rainy than the intermediate site, was nearly zero (Figure 4c).

Discussion
The most statistically significant climatic variables in terms of the response of survival and growth in A. religiosa seedlings to the climatic transfer distance were, in all cases, the balance between the temperature useful to growth (degree days above 5 °C) and available precipitation. That balance, expressed as dryness indices in our analysis, varied only in different seasons for the response variables: the rainy season (June-October) dryness index transfer distance was critical for survival, while that of the cold and dry season (November-February) was critical for aerial biomass, and the annual index was critical for the increase in basal diameter.
The effect of climatic transfer distance was much more negative (in terms of seedling survival) when transferring toward warmer and dryer sites, than when shifting to colder and moister sites. Apparently, there is a critical threshold that, when crossed, triggers mortality: a shift to a site 400 m lower in altitude, equivalent to +2 °C considering an overall lapse rate [34], or a difference of +1.9 °C based on the observed months, or a reduction in precipitation of 16%.
Most of the seedling mortality at the low altitude site (2400 m a.s.l.) occurred during the dry season (December-May). This is consistent with previous observations of a high mortality of emergent seedlings of natural regeneration in the core zone of the Monarch Butterfly Biosphere Reserve (MBBR) [29], and seems to support the projected decrease in climatic habitat suitable for A. religiosa inside the MBBR as a result of climatic change [26,30,[49][50][51]. The apparent reason for the observed timing of this mortality is that the dry season was as dry as always, perhaps a little more so, but was warmer due to the ongoing climatic change. It was surprising to us that the relatively small differences in temperatures and precipitation between the sites at 2400 m a.s.l. and 3000 m a.s.l. made such a large difference to survival. Thus, the amount of residual moisture in the soil during the dry season, following the rainy season of the previous year, seems to be critical for survival of the young seedlings.
The high survival (around 95%) of the seed sources shifted to higher altitudes than their origin (to the site at 3400 m a.s.l.) confirms the viability of assisted migration to higher altitudes in the Trans-Mexican Volcanic Belt as a management measure for consideration,

Discussion
The most statistically significant climatic variables in terms of the response of survival and growth in A. religiosa seedlings to the climatic transfer distance were, in all cases, the balance between the temperature useful to growth (degree days above 5 • C) and available precipitation. That balance, expressed as dryness indices in our analysis, varied only in different seasons for the response variables: the rainy season (June-October) dryness index transfer distance was critical for survival, while that of the cold and dry season (November-February) was critical for aerial biomass, and the annual index was critical for the increase in basal diameter.
The effect of climatic transfer distance was much more negative (in terms of seedling survival) when transferring toward warmer and dryer sites, than when shifting to colder and moister sites. Apparently, there is a critical threshold that, when crossed, triggers mortality: a shift to a site 400 m lower in altitude, equivalent to +2 • C considering an overall lapse rate [34], or a difference of +1.9 • C based on the observed months, or a reduction in precipitation of 16%.
Most of the seedling mortality at the low altitude site (2400 m a.s.l.) occurred during the dry season (December-May). This is consistent with previous observations of a high mortality of emergent seedlings of natural regeneration in the core zone of the Monarch Butterfly Biosphere Reserve (MBBR) [29], and seems to support the projected decrease in climatic habitat suitable for A. religiosa inside the MBBR as a result of climatic change [26,30,[49][50][51]. The apparent reason for the observed timing of this mortality is that the dry season was as dry as always, perhaps a little more so, but was warmer due to the ongoing climatic change. It was surprising to us that the relatively small differences in temperatures and precipitation between the sites at 2400 m a.s.l. and 3000 m a.s.l. made such a large difference to survival. Thus, the amount of residual moisture in the soil during the dry season, following the rainy season of the previous year, seems to be critical for survival of the young seedlings.
The high survival (around 95%) of the seed sources shifted to higher altitudes than their origin (to the site at 3400 m a.s.l.) confirms the viability of assisted migration to higher altitudes in the Trans-Mexican Volcanic Belt as a management measure for consideration, in order to reduce the projected impacts of climatic change, as shown in previous field tests [33].
Suitable climatic habitat projections for the 2060 decade, under intermediate greenhouse emission scenarios (RCP 6.0 watts/m 2 ), indicate that it is required to shift the A. religiosa seed sources upwards in altitude by between 300 and 500 m in order to realign the populations with the climate to which they are adapted [30]. Our evidence that is safe to move these sources upwards in altitude by 400 m therefore supports such a management proposal.
The mean temperature of the coldest month (MTCM) was the most important climate variable at the seed source for the increase in basal diameter. This indicates that genetic differentiation among populations for growth potential is shaped by the cold temperatures, at least in our case for those seedlings that survived the dry season. This is consistent with previous patterns of genetic variation among A. religiosa populations along altitudinal clines demonstrated in a common garden test under optimal conditions, where MTCM was also found to be the climatic variable that best explained the differences among populations [20].
Aerial biomass necessarily had to be measured, at the end of the observation period, in live plants that had survived the transfer toward warmer and dryer sites. It is therefore likely that those seedlings might have had some resistance or tolerance to the drought stress, and growing in a site with higher temperature (2400 m a.s.l.), they were able to grow more. In contrast, although the proportion of seedlings that survived in those that were transferred to a colder and moister site was higher (as in Figure 2a), growth was likely limited by the colder temperatures.

Implications for Management of Abies religiosa Inside the MBBR
The results indicate that increased temperatures and reduced precipitation due to climatic change would undoubtedly have a negative impact on the recruitment of A. religiosa seedlings at the lower altitudinal limits of the natural distribution range of the species, through increased mortality during the dry season (December-May). Active ecological restoration efforts in the MBBR must therefore consider how to increase the survival of planted A. religiosa seedlings in order to restore perturbed sites and ensure the maintenance of the monarch butterfly overwintering sites. In this context, perturbed sites at the low altitude of the natural range distribution of A. religiosa should have priority. Actions should include: (a) conducting reforestation using local shrubs (such as Baccharis conferta or other native understory species) as nurse plants for the A. religiosa seedlings [33]; and (b) conducting assisted migration to higher altitudes to realign A. religiosa populations with the climate to which they are adapted. Previous estimates regarding such realignment, targeting the climate projected for the 2030 decade, have suggested an upwards shift in altitude of 350 m [20]. The results of the present study suggest that it is safe to shift up to 400 m upwards in altitude, in terms of good seedling survival. Thus, an upward shift of between 350 and 400 m in altitude seems to be advisable. If such action is not undertaken, the lack of recruitment of natural regeneration seedlings is likely to continue as a result of the high mortality of young seedlings at lower altitudes (as demonstrated by ref [29]), impeding the natural renewal of adult A. religiosa trees within the MBBR, which will endanger the monarch butterfly overwintering sites in the long term.

Conclusions
The balance between the seasonal temperature available for growth (degree days > 5 • C) and the available precipitation is a determining factor for the permanence of A. religiosa, because the survival of young seedlings depends on the residual humidity of the soil during the dry season, following the rainy season of the previous year. A 2 • C increase in temperature and a 16% decrease in precipitation would imply a high mortality rate (~45%), negatively impacting the recruitment of naturally regenerating seedlings, mainly at the lower altitudinal limit of the natural distribution of the species. These results suggest that an upward movement of 400 m in altitude, through assisted migration, will have positive effects on the survival of A. religiosa seedlings, re-coupling the populations to the climate to which they are adapted and increasing the possibility of maintaining the overwintering sites of the monarch butterfly.