Received 13 March 2016; accepted 23 April 2016; published 26 April 2016
Phenological sensitivity reflects the ability of a species to track climate change and indicates whether the species will persist or decline in response to climate changes  . Phenological events are easily observable and most sensitive traits to climate change  . Plants responded to changes in growing conditions and altered availability of resources induced by climate change by making phenotypic and physiological adjustments, or by migrating to a new area having similar environmental conditions that they are adapted to, or by combinations of the above responses and eventually by genetically adapting to the new conditions  . Historical evidence indicates that all the above mechanisms have occurred in plants’ responses to climate change  . During the process of migration and adaptation plants encounter multiple challenges. Since it is almost impossible to get exactly the same environmental conditions to migrate to, plants usually migrate to a relatively favorable site and phenotypic plasticity plays a critical role in the success of such a migration. Phenotypically, plastic changes have a large influence on the immediate fate of the migrating plants via their effects on survival and reproduction  .
Each environmental factor such as atmospheric carbon dioxide concentration [CO2], temperature and photoperiod has its own impacts on the growth and biomass allocations in trees  -  . Increased atmospheric [CO2] generally increases the photosynthetic capacity of plants leading to an increased carbon balance, growth and yield  -  because of increased photosynthetic carbon acquisition and changes in biomass allocation  . Elevated [CO2] can markedly alter the structure and physiology of plants, such as increased leaf expansion and cell wall extensibility and often cell turgor pressure, leading to increased leaf and root growth. Increased atmospheric [CO2] can increase total leaf area  , leaf mass  , leaf mass to area ratio  , and branching frequency  . Elevated atmospheric [CO2] enhances root growth   and affects many important root traits  . Root biomass, root length, root branching and lateral root production are reported to increase in response to elevated [CO2]    . Elevated [CO2] can affect the growth rhythm of trees by altering the timing of bud burst and growth cessation and results in an increased seedling cold tolerance in early fall   . It can significantly advance the date of bud burst and increase total foliage area as a result of increased number and length of shoots and increased individual foliage area  . Early bud burst, however, may increase the risk of frost damage from late spring frosts  .
Soil temperature influences the growth and physiology in plants, such as biomass allocation, shoot and root growth, stomatal conductance and nutrient and water uptake   -  . The soil temperature in the boreal forest is often too low for the optimal physiological activities and growth  . For example, low soil temperature can inhibit nutrient and water uptake, limit shoot growth   and reduce photosynthetic productivity  . Therefore, the increase of soil temperature by climatic warming would be beneficial to boreal plants    . Warmer soil temperatures improve water uptake, root growth and root permeability  .
Photoperiod regulates annual growth cycles in temperate and boreal trees and influences their geographic distribution  -  . Plants are generally adapted to the photoperiod regime of their habitats over hundreds and thousands of years  . Photoperiod regulates the development of cold hardness in the fall and the timing of bud flush in the spring    . However, photoperiod is not considered in climate envelop models. If species immigrate naturally or with human assistance based on the shift of their climate envelopes, it is possible that the species may not be able to acclimate to the photoperiod regimes at their new location.
Although the individual effects of elevated atmospheric [CO2], warmer soil temperature and prolonged photoperiod on plant growth are mostly positive, their combined effects may not be equal to the sum of the individual effects   . Soil temperature and light regime influence the responses of plants to elevated [CO2]   -  . However, elevated [CO2] in combination with increased temperature have no significant effect on the overall height growth of jack pine seedlings, but significantly increase the total biomass  . Cantin et al. (1997)  report that the growth of jack pine seedlings under elevated [CO2] and temperature results in reduced leaf mass ratio as more biomass is allocated to roots. The effects of individual environmental factors or bi-factor effects on boreal tree growth are well documented    -  , but the understanding of interactive effects of multiple environmental variables such as [CO2], soil temperature and photoperiod, on boreal tree growth is still limited.
Jack pine (Pinus banksiana Lamb.) is a major tree species in the boreal forests of Canada holding great ecological and commercial values and thus deserves special attention in the context of climate change. Jack pine is the second most planted tree species in Canada after black spruce   . In Ontario jack pine comprises approximately 37% of the total annual softwood harvest  . Atmosphere Ocean General Circulation Models predict a 10 degree (approximately 1000 kilometers ) northward shifts in the climate envelopes of 130 North American tree species between 2071 and 2100   . Following the predicted shift in climate envelopes jack pine might need to migrate 10˚ northward between 2071 and 2100  . In that case, the species will be exposed to a different photoperiod regime, e.g. the photoperiods will be longer in the summer and shorter in the winter with faster transition between seasons than the regimes that it has adapted to, which might affect the phenological events of the species. But, the impacts of changes in photoperiod regimes associated with migration or seed transfer of jack pine are not yet well documented. Since the impacts of elevated atmospheric [CO2] and warmer soil temperature on the growth of jack pine are mostly positive    while autumn leaf senescence and cold hardening are predominantly controlled by photoperiod   and the time period for cold hardening becomes shorter because of the faster rate of photoperiod declines in the fall at higher latitudes, we hypothesize that elevated [CO2] and soil temperature and photoperiods associated with northward migration or seed transfer will enhance the growth performance of jack pine but reduce the cold hardness. The overall objective of this study was to assess the interactive effects of photoperiod, soil temperature and CO2 elevation on the morphological traits of jack pine in the context of northward migration. The specific objectives of the study were to: 1) assess the tri-factor effects on the phenological traits of jack pine seedlings, 2) assess the impacts of elevated [CO2], soil temperature and photoperiod on the growth performance and biomass allocation patterns and 3) investigate the spatial extent of possible northward migration following the predicted shift in climate envelope.
2. Materials and Methods
2.1. Plant Materials
One-year-old jack pine (Pinus banksiana Lamb.) seedlings were obtained from a commercial tree seedling nursery. The seedlings were raised from seeds collected from Kakabeka region (48˚57'N & 90˚44'W). A total of 608 seedlings of relatively uniform size (average height 13.05; averageroot collar diameter 0.21 cm) were planted in plastic pots of size 15cm in height and 13cm in diameter, filled up with a mixture of premium grade vermiculite and peat moss (50:50, v/v).
2.2. Experimental Design
The experiment was set on November 20, 2013 and continued until May 31, 2014 mimicking the actual day lengths and temperatures for the period of April 16 to October 31 (Table 1).
The experiment was conducted in four green houses at Lakehead University’s Thunder Bay campus. The treatments were consisted of two levels of [CO2] (400 and 950 µmol∙mol−1), two soil temperatures (the average soil temperature at seed origin (Ts) and 5˚C warmer (Te)) and three photoperiod regimes (Ps, Pnm1 and Pnm2, indicating photoperiods at seed origin, at 5˚ north of the seed origin and 10˚ north of the seed origin, respectively). The 10˚ north of the seed origin represents the predicted northward shift of the climate envelopes for boreal tree species   and the 5˚ north is the midpoint between the current and the predicted future location of the species. The soil temperature of seed origin was derived from the biweekly mean soil temperature of the years 2012 and 2013 at Kakabeka region. Average temperatures at 5 cm and 20 cm soil depths were used in the experiment. Based on the biweekly average, soil temperatures were adjusted periodically. The experiment was carried out following a split plot design, with [CO2] as the whole plot, soil temperature as the sub-plot and photoperiod as the sub-sub plot. Two levels of [CO2] were randomly assigned to four independent greenhouses with two replicates of each. The two levels of soil temperature were duplicated within each level of [CO2] and then three levels of photoperiod were nested within each soil temperature.
2.3. Environmental Controls
To elevate [CO2] in the greenhouses, electronic ignition natural gas CO2 generators (model GEN-2E, Custom
Table 1. Biweekly mean day length and mean air and soil temperature for the period of April 16 to October 31. Air temperature is based on 10 years average (2004-2013). Soil temperatures are biweekly averages of soil temperatures at 5 and 20 cm depths in 2012 and 2013.
Note: Ps, Pnm1 & Pnm2 represent the seed origin, 5˚ north of seed origin & 10˚ north of seed origin, respectively.
Automated Products Inc, Riverside, CA) were used. Circulation fans were used to make the air composition even throughout the greenhouse. The [CO2] in each greenhouse was monitored and controlled automatically with an Argus Environment Control System (Argus, Vancouver, BC, Canada). The soil temperatures were controlled using a soil temperature control system consisting of a large leak-proof wooden box (196 cm long, 112 cm wide and 16 cm deep) designed by Cheng et al. (2000)  . The seedling containers were mounted within the soil temperature control box and temperature-controlled water was circulated in the space between the containers to maintain the desired soil temperature. To ensure the even distribution of temperature within the system, each control unit was equipped with circulatory pumps (model AC-2CP-MD, March Mfg. Inc., Glenview, Illinois, USA). The system was insulated to minimize heat exchange with the greenhouse air. Each day-time soil temperature was lowered by approximately 4˚C - 6˚C at night to cater for lower night temperature. The lengths of the photoperiod for active growing and cold hardening phases for each of the three locations were set periodically based on the actual biweekly average photoperiods during the growing season and cold hardening phase at corresponding locations. Each photoperiod treatment was started with summer long days followed by the corresponding short days. The length of the natural photoperiod was extended using high-pressure sodium lamps when natural day lengths were shorter than the set photoperiods. Photoperiods shorter than the natural day length were achieved through manual shading of the seedlings in the early mornings and late evenings following standard blackout techniques used in tree nurseries. All the seedlings were fertilized biweekly with N:P:K fertilizer (7.14, 0.56 and 2.12 mmol per liter of water, respectively during the rapid growth phase and 1.78, 2.20 and 4.64 mmol per liter of water, respectively during the hardening phase)  . The moisture level was maintained at 60% - 70% of the field water capacity. The moisture content of the growing medium was monitored daily using a Delta-T ML2x probe and HH2 moisture meter (Delta-T Devices, Cambridge, UK). In all the greenhouses the relative humidity was maintained at 55% - 60% during the rapid growth phase and reduced to 45% - 50% during the hardening phase. The day and night air temperatures were maintained mimicking the mean biweekly day and night air temperatures at the seed origin. Relative humidity and air temperature were controlled automatically using the Argus Environment Control System.
Height and root collar diameter (RCD) of each individual seedling were measured at the beginning of the experiment. Ten seedlings from each treatment combination were randomly chosen to investigate bud break at the beginning and bud setting during the cold hardening phase. Bud break was determined by visual observation of the terminal meristem of seedlings. Bud break was considered completed when bud scales were mostly fallen, with tips of needles protruding about 2 mm  . Days to bud break was counted from the date of starting the treatment (November 20, 2013). Bud set was also determined by visual observation and it was considered accomplished when bud scales were completely closed  . Bud setting time was counted from April 16, 2014, when the day length and day and night air temperature were reduced to 12 hours and 14˚C and 8˚C, respectively.
Heights and RCDs of three randomly chosen seedlings from each treatment combination were measured on May 28, 2014. The seedlings were then separated into foliage, stem (including branches) and roots. The needles were scanned with WinSeedle system (Regent Instruments Inc, Quebec, Canada) to determine the projected leaf area. The roots were scanned to analyze root traits using a WinRhizo system (Regent Instruments Inc, Quebec, Canada). The foliage, stem and roots were subsequently oven-dried at 70˚C for 48 hours and weighed on an analytical balance (0.001 g precision) to determine the dry mass. Stem mass ratio (ratio of stem mass to total biomass, SMR), root mass ratio (ratio of root mass to total biomass, RMR), leaf mass ratio (ratio of leaf mass to total biomass, LMR), shoot to root mass ratio (ratio of above ground to belowground mass, SRR) were used as indices of biomass allocation. Specific leaf area (SLA) was determined from the projected leaf area and leaf dry mass. Specific root length (SRL) and specific root surface area (SRA) were determined as root length per unit root dry mass and root surface area per unit root dry mass, respectively. Root length to leaf area ratio (RLA) was calculated to estimate the water supply/demand index. The stem Volume (V) was calculated from height (H) and RCD (D) using the equation  : V = (πD2/4) H/3. The relative growth rate of height (HT), root-collar diameter (RCD) and Volume (Vol) were determined by dividing the increment by the corresponding initial measurement.
2.5. Statistical Analysis
All data were analyzed using Analysis of Variance (ANOVA). Before the ANOVA test, Shapiro-Wilk and Bartlett tests were conducted to test the normality of distribution and homogeneity of variance, respectively. Whenever necessary, suitable transformation was done to normalize the data. Tukey’s HSD Post-hoc comparisons were carried out when ANOVA showed a significant interaction (P < 0.05). In the analyses, [CO2], soil temperature and photoperiod were treated as fixed factors. All analyses were done using R 3.1.3  .
On the average, budburst advanced for 5 days by elevated [CO2] (Ce) and 2 days by warmer soil temperature (Te) compared to those in current [CO2] (Ca) and soil temperature (Ts) (Figure 1(a) and Figure 1(b)). The effect of soil temperature was mostly in Ce (significant CxT interaction (P < 0.05) see Table 1 and Figure 1(d)). The total day to budburst was significantly less with the photoperiod regimes associated with northward migration: 34.63 ± 0.55, 32.63 ± 0.55 and 29.69 ± 0.49 days for Ps, Pnm1 and Pnm2, respectively (Figure 1(c)). The examination of the 3-way interaction (P < 0.05) indicated the longest days to budburst for Ca × Ts × Ps treatment combination (38.7 ± 0.874) and the shortest for Ce × Te × Pnm2 combination (25.7 ± 0.7) (Figure 1(e)). The total day to bud setincreased with photoperiod (Ps) from Ps (32.16 ± 0.2 days) to Pnm2 (34.58 ± 0.157 days) (p < 0.001, Figure 1(f)) and the variation among Ca × Te within Ps was marginally significant effect (p = 0.087).
The relative growth rate of height (RGRH) increased (P < 0.05) with the length of photoperiod (1.24 ± 0.04 for Pnm2 vs. 1.12 ± 0.03 for Ps, Figure 2(a)). Similarly, the relative growth rate of RCD (RGRD), relative growth rate of stem volume growth (RGRV), and seedling total biomass (M) became greater by elevated [CO2] (Figure 2(b), Figure 2(d), Figure 2(e)). The significant CxT interaction (P < 0.05) on RGRD indicated greater [CO2] effect under current soil temperature (see Table 2, Figure 2(c)).
The biomass allocation of jack pine seedlings was significantly affected by atmospheric [CO2] and soil temperature (P < 0.05, Table 2). The total projected leaf area was significantly greater (Figure 3(a)), and leaf mass ratio (0.34 ± 0.00 vs. and 0.29 ± 0.00), shoot to root mass ratio (1.56 ± 0.05 vs. 1.31 ± 0.04) (Figure 3(b) and Figure 3(c)), specific root length (10.57 ± 0.38 vs. 7.35 ± 0.15 m/gm, and specific root area (219.63 ± 6.02 vs. 159.29 ± 4.52 cm2/gm) (Figure 3(d) & Figure 3(e)) were smaller under elevated [CO2]. The root length to leaf
Figure 1. Least square means (+SE) of total days to budburst and bud set in jack pine seedlings grown under two levels of [CO2] (Ca = 400 & Ce = 950 µmol∙mol−1), two soil temperature (Ts, seed origin and Te, 5˚ warmer than Ts) and three photoperiods (Ps, seed origin and Pnm1 and Pnm2, values for 5˚ & 10˚ northward migration, respectively). Different lowercase letters indicate statistically significant differences at P ≤ 0.05.
Figure 2. Least square means (+SE) of relative growth rates for height, root collar diameter (RCD) and stem volume, and total seedling biomass in jack pine seedlings grown under two levels of [CO2] (Ca & Ce), two soil temperature (Ts & Te) and three photoperiods (Ps, Pnm1 & Pnm2). Different lowercase letters indicate statistically significant differences at P ≤ 0.05.
Figure 3. Least square means (+SE) of total leaf area, leaf mass ratio, shoot to root mass ratio, specific root length, specific root surface area and root length to leaf area ratio in jack pine seedlings grown under two levels of [CO2] (Ca & Ce), two soil temperature (Ts & Te) and three photoperiods (Ps, Pnm1 & Pnm2). Different lowercase letters indicate statistically significant differences at P ≤ 0.05.
The advance of budburst in jack pine seedlings with [CO2] elevation, soil temperature warming and photoperiod increase is consistent with the results of many predictions  -  , but different from others including delayed or no responses to elevated [CO2]   -  and to elevated soil temperature on various tree species  -  . As budburst is primarily controlled by temperature and longer photoperiod promotes dormancy release and budburst only in some species     -    , the advance of budburst by Pnm2 may be more due to the greater thermal accumulation under longer daytime hours in spring time. The individual effects of elevated [CO2], soil temperature and longer growing season photoperiods appeared to be additive; the budburst in the Ce × Te × Pnm2 combination was 13 days earlier than that in the Ca × Ts × Ps treatment combination. Advances in budburst will likely increase the risk of spring frosts, which may be an obstacle to the northward migration of jack pine.
The lack of treatment effect by [CO2] and soil temperature levels is supported by the general believe that bud set is induced largely by shortening photoperiod   , although a positive effect of elevated [CO2] has been reported on some tree species including black spruce  and Sitka spruce  . The earlier bud set in jack pine seedlings by the photoperiod regimes under northward migration is probably associated with their shorter photoperiod in fall time (Table 1).
As expected, the elevated [CO2] increased the relative growth rates of root collar diameter (RGRD) and stemvolume (RGRV) as well as total seedling biomass (TB). In C3 plants the rate of net CO2 assimilation is not
Table 2. ANOVA P-values for the effects of carbon-dioxide concentration (C), soil temperature (T), photoperiod (P) and their interactions on total days to budburst and bud set, relative growth rate of height, root collar diameter and volume (RGRH, RGRD and RGRV, respectively), total biomass (M), total leaf area (projected) per seedling (LA), stem mass ratio (SMR), root mass ratio (RMR), leaf mass ratio (LMR), shoot to root mass ratio (SRR), specific leaf area (SLA), specific root length (SRL), specific root surface area (SRA) and root length to leaf area ratio (RLA). Seedlings were grown under two levels of [CO2], two levels of soil temperature and three photoperiod regimes.
Note: Significant effect (P ≤ 0.05) was bold-faced.
CO2-saturated at the current level of atmospheric [CO2]  and higher [CO2] would positively affect growth if nutrients are not limiting (as would be in this study). The enhanced growth by elevated [CO2] may also be related to the reductions in dark respiration and earlier budburst and therefore longer growing season      -  . Similarly, RGRH was greater under photoperiod regimes associated with higher latitudes in jack pine seedlings, possibly due to longer hours of photosynthesis. We did not find significant increases of RGRH with elevated [CO2] and seedling growth with warmer soils, as has been reported by others     -  . We speculate that the soil temperature at seed origin is not too low to affect the physiological activities and growth of jack pine seedlings, as indicated by others that soil temperature effect is more species-specific    .
As found by others       , elevated [CO2] affected biomass allocation and functional traits of leaf and root in jack pine seedlings such as lower leaf mass ratio (LMR), shoot to root ratio (SRR), specific root length (SRL) and specific root surface area (SRA), even though not all of the treatment effects (e.g., specific leaf area (SLA)) reached the level of significance.
Though the main effects on phenological traits, growth and biomass allocation were mostly significant, interactive effects were generally statistically insignificant in this study with the exception of budburst and RCD. Rather than suggesting the lack of real significant interactions, these results indicate the difficulties in detecting interactive effects when multiple environmental variables are being investigated. As the number of treatment factors increase, the degree of freedom for the error term used to test interactions declines, causing a loss of the capability of the experiment to detect interactive effects. A remedy to the problem will be to increase the sample size. However, such a remedy often faces great logistical challenges. A more practical approach for investigating the interactive effects of multiple factors is yet to be found.
 Cleland, E.E., Allen, J.M., Crimmins, T.M., Dunne, J.A., Pau, S., Travers, S.E., Zavaleta, E.S. and Wolkovich, E.M. (2012) Phenological Tracking Enables Positive Species Responses to Climate Change. Ecology, 93, 1765-1771.
 Soolanayakanahally, R.Y., Guy, R.D., Silim, S.N. and Song, M. (2013) Timing of Photoperiodic Competency Causes Phenological Mismatch in Balsam Poplar (Populus balsamifera L.). Plant, Cell and Environment, 36, 116-127.
 Aitken, S.N., Yeaman, S., Holliday, J.A., Wang, T.L. and Curtis-McLane, S. (2008) Adaptation, Migration or Extirpation: Climate Change Outcomes for Tree Populations. Evolutionary Applications, 1, 95-111.
 Johnston, M.H., Campagna, M., Gray, P.A., Kope, H.H., Loo, J.A., Ogden, A.E., O’Neill, G.A., Price, D.T. and Williamson, T.B. (2009) Vulnerability of Canada’s Tree Species to Climate Change and Management Options for Adaptation: An Overview for Policy Makers and Practitioners. Canadian Council of Forest Ministers. Northern Forestry Centre, 44.
 Huang, J.G., Bergeron, Y., Denneler, B., Berninger, F. and Tardif, J. (2007) Response of Forest Trees to Increased Atmospheric CO2. Critical Reviews in Plant Sciences, 26, 265-283.
 Lukac, M., Calfapietra, C., Lagomarsino, A. and Loreto, F. (2010) Global Climate Change and Tree Nutrition: Effects of Elevated CO2 and Temperature. Tree Physiology, 30, 1209-1220.
 Ainsworth, E.A. and Long, S.P. (2005) What Have We Learned from 15 Years of Free-Air CO2 Enrichment (FACE)? A Meta-Analytic Review of the Responses of Photosynthesis, Canopy Properties and Plant Production to Rising CO2. New Phytologist, 165, 351-372.
 Curtis, P.S. and Wang, X. (1998) A Meta-Analysis of Elevated CO2 Effects on Woody Plant Mass, Form, and Physiology. Oecologia, 113, 299-313.
 Idso, K.E. and Edso, S.B. (1994) Plant Responses to Atmospheric CO2 Enrichment in the Face of Environmental Constraints: A Review of the Past 10 Years’ Research. Agricultural and Forest Meteorology, 69, 153-203.
 Marfo, J. and Dang, Q.L. (2009) Interactive Effects of Carbon Dioxide Concentration and Light on the Morphological and Biomass Characteristics of Black Spruce Seedlings. Botany, 87, 67-77.
 Norby, R.J., Wullschleger, S.D., Gunderson, C.A., Johnson, D.W. and Ceulemans, R. (1999) Tree Responses to Rising CO2 in Field Experiments: Implication for Future Forest. Plant, Cell and Environment, 22, 683-714.
 Pregitzer, K.S., King, J.S., Andrew, J.B. and Brown, E. (2000) Responses of Tree Fine Roots to Temperature. New Phytologist, 147, 105-115.
 Eamus, D. and Jarvis, P.G. (1989) The Direct Effects of Increase in the Global Atmospheric CO2 Concentration on Natural and Commercial Temperate Trees and Forests. Advances in Ecological Research, 19, 1-55.
 Koch, K.E., Jones, P.H., Avigne, W.T. and Allen, L.H. (1986) Growth and Dry Matter Partitioning and Diurnal Activities of RuBP Carboxylase in Citrus Seedlings Maintained at Two Levels of CO2. Physiologia Plantarum, 67, 477-484.
 Norby, R.J. and O’Neill, E.G. (1989) Growth Dynamics and Water Use of Seedlings of Quercus alba L. in CO2-Enriched Atmospheres. New Phytologist, 111, 491-500.
 Pettersson, R., McDonald, A.J.S. and Stadenberg, I. (1993) Response of Small Birch Plants (Betula pendula Roth.) to Elevated CO2 and Nitrogen Supply. Plant, Cell & Environment, 16, 1115-1121.
 Rogers, H.H., Runion, G.B. and Krupa, S.V. (1994) Plant Responses to Atmospheric CO2 Enrichment with Emphasis on Roots and the Rhizosphere. Environmental Pollution, 83, 155-189.
 Tingey, D.T., Phillips, D.L. and Johnson, M.G. (2000) Elevated CO2 and Conifer Roots: Effects on Growth, Life Span and Turnover. New Phytologist, 147, 87-103.
 Day, F.P., Weber, E.P., Hinckle, C.R. and Drake, B.G. (1996) Effects of Elevated Atmospheric CO2 on Fine Root Length and Distribution in an Oak-Palmetto Scrub Ecosystem in Central Florida. Global Change Biology, 2, 143-148.
 Janssens, I.A., Crookshanks, M., Taylor, G. and Ceulemans, R. (1998) Elevated CO2 Increases Fine Root Production, Respiration, Rhizosphere Respiration and Soil CO2 Efflux in Scots Pine Seedlings. Global Change Biology, 4, 871-878.
 Bigras, F.J. and Bertrand, A. (2006) Responses of Piceamariana to Elevated CO2 Concentration during Growth, Cold Hardening and Dehardening: Phenology, Cold Tolerance, Photosynthesis and Growth. Tree Physiology, 26, 875-888.
 Ceulemans, R., Jiang, X.N. and Shao, B.Y. (1995) Effects of Elevated Atmospheric CO2 on Growth, Biomass Production and Nitrogen Allocation of Two Populus Clones. Journal of Biogeography, 22, 261-268.
 Jach, M.E. and Ceulemans, R. (1999) Effects of Elevated Atmospheric CO2 on Phenology, Growth and Crown Structure of Scots Pine (Pinus sylvestris) Seedlings after Two Years of Exposure in the Field. Tree Physiology, 19, 289-300.
 Murray, M.B., Cannell, M.G.R. and Smith, R.I. (1989) Date of Bud Burst of Fifteen Tree Species in Britain Following Climate Warming. Journal of Applied Ecology, 26, 693-700.
 Alvarez-Uria, P. and Korner, C. (2007) Low Temperature Limits of Root Growth in Deciduous and Evergreen Temperate Tree Species. Functional Ecology, 21, 211-218.
 Cai, T.B. and Dang, Q.L. (2002) Effects of Soil Temperature on Parameters of a Coupled Photosynthesis-Stomatal Conductance Model. Tree Physiology, 22, 819-828.
 Camm, E.L. and Harper, G.J. (1991) Temporal Variations in Cold Sensitivity of Root Growth in Cold-Stored White Spruce Seedlings. Tree Physiology, 9, 425-431.
 Dang, Q.L. and Cheng, S. (2004) Effects of Soil Temperature on Ecophysiological Traits in Seedlings of Four Boreal Tree Species. Forest Ecology and Management, 194, 379-387.
 Landhausser, S.M., DesRochers, A. and Lieffers, V.J. (2001) A Comparison of Growth and Physiology in Picea glauca and Populus tremuloides at Different Soil Temperatures. Canadian Journal of Forest Research, 31, 1922-1929.
 Peng, Y.Y. and Dang, Q.L. (2003) Effects of Soil Temperature on Biomass Production and Allocation in Seedlings of Four Boreal Tree Species. Forest Ecology and Management, 180, 1-9.
 Tierney, G., Fahey, T., Groffman, P., Hardy, J., Fitzhugh, R. and Driscoll, C. (2001) Soil Freezing Alters Fine Root Dynamics in A Northern Hardwood Forest. Biogeochemistry, 56, 175-190.
 Ambebe, T.F., Dang, Q.L. and Li, J. (2010) Low Soil Temperature Inhibits the Effects of High Nutrient Supply on Photosynthetic Response to Elevated Carbon Dioxide Concentration in White Birch Seedlings. Tree Physiology, 30, 234-243.
 Domisch, T., Finer, L., Lehto, T. and Smolander, A. (2002) Effect of Soil Temperature on Nutrient Allocation and Mycorrhizas in Scots Pine Seedlings. Plant and Soil, 239, 173-185.
 Haussling, C.A., Jorn, C.A., Lehmbecker, G., Hecht-Buchholz, C. and Marschner, H. (1988) Ion and Water Uptake in Relation to Root Development in Norway Spruce (Picea abies [L.] Karst.). Journal of Plant Physiology, 133, 486-491.
 Campbell, R.K. and Sugano, A.I. (1975) Phenology of Bud Burst in Douglas-Fir Related to Provenance, Photoperiod, Chilling, and Flushing Temperature. International Journal of Plant Sciences, 136, 290-298.
 Campbell, R.K. (1979) Genecology of Douglas-Fir in a Watershed in the Oregon Cascades. Ecology, 60, 1036-1050.
 Weiser, C.J. (1970) Cold Resistance and Injury in Woody Plants. Science, 169, 1269-1278.
 Hanninen, H. (2006) Climate Warming and the Risk of Frost Damage to Boreal Forest Trees: Identification of Critical Ecophysiological Traits. Tree Physiology, 26, 889-898.
 Colombo, S.J., Zhao, S.Y. and Blumwald, E. (1995) Frost Hardiness Gradients in Shoots and Roots of Picea mariana Seedlings. Scandinavian Journal of Forest Research, 10, 32-36.
 Fracheboud, Y., Luquez, V., Bjorkén, L., Sjodin, A., Tuominen, H. and Jansson, S. (2009) The Control of Autumn Senescence in European Aspens. Plant Physiology, 149, 1982-1991.
 Curtis, P.S., Vogel, C.S., Wang, X., Pregistzer, K.S., Zak, D.R., Lussenhop, J., Kubiske, M. and Teeri, J.A. (2000) Gas Exchange, Leaf Nitrogen, and Growth Efficiency of Populus tremuloides in a CO2-Enriched Atmosphere. Ecological Applications, 10, 3-17.
 Nowak, R.S., Ellsworth, D.S. and Smith, S.D. (2004) Functional Responses of Plants to Elevated Atmospheric CO2—Do Photosynthetic and Productivity Data from FACE Experiments Support Early Predictions? New Phytologist, 162, 253-280.
 Allen, S.G., Idso, S.B. and Kimball, B.A. (1990) Interactive Effects of CO2 and Environment on Net Photosynthesis of Water-Lily. Agriculture, Ecosystems & Environment, 30, 81-88.
 Mishra, R.S., Abdin, M.Z. and Uprety, D.C. (1999) Interactive Effects of Elevated CO2 and Moisture Stress on the Photosynthesis, Water Relation and Growth of Brassica Species. Journal of Agronomy & Crop Science, 182, 223-229.
 Robredo, A., Pérez-López, U., Sainz, H., de la Maza, B., González-Moro, M.L., Mena-Petite, A. and Munoz-Rueda, A. (2007) Elevated CO2 Alleviates the Impact of Drought on Barley Improving Water Status by Lowering Stomatal Conductance and Delaying Its Effects on Photosynthesis. Environmental and Experimental Botany, 59, 252-263.
 Zebian, K.J. and Reekie, E.G. (1998) The Interactive Effects of Atmospheric Carbon Dioxide and Light on Stem Elongation in Seedlings of Four Species. Annals of Botany, 81, 185-193.
 Cantin, D., Tremblay, M.F., Lechowicz, M.J. and Potvin, C. (1997) Effects of CO2 Enrichment, Elevated Temperature, and Nitrogen Availability on the Growth and Gas Exchange of Different Families of Jack Pine Seedlings. Canadian Journal of Forest Research, 27, 510-520.
 DeLucia, E.H. and Smith, W.K. (1987) Air and Soil Temperature Limitations on Photosynthesis in Engelmann Spruce during Summer. Canadian Journal of Forest Research, 17, 527-533.
 Li, J., Dang, Q.L. and Man, R. (2015) Photoperiod and Nitrogen Supply Limit the Scope of Northward Migration and Seed Transfer of Black Spruce in a Future Climate Associated with Doubled Atmospheric CO2 Concentration. American Journal of Plant Sciences, 6, 189-200.
 Way, D.A. and Sage, R.F. (2008) Thermal Acclimation of Photosynthesis in Black Spruce [Picea mariana (Mill.) B.S.P.]. Plant, Cell and Environment, 31, 1250-1262.
 Way, D.A. and Sage, R.F. (2008) Elevated Growth Temperatures Reduce the Carbon Gain of Black Spruce [Picea mariana (Mill.) B.S.P.]. Global Change Biology, 14, 624-636.
 McKenney, D.W., Pedlar, J.H., Lawrence, K.E.V.I., Campbell, K. and Hutchinson, M.F. (2007) Potential Impacts of Climate Change on the Distribution of North American Trees. BioScience, 57, 939-948.
 McKenney, D.W., Pedlar, J.H., Rood, R.B. and Price, D.A.V.I. (2011) Revisiting Projected Shifts in the Climate Envelopes of North American Trees Using Updated General Circulation Models. Global Change Biology, 17, 2720-2730.
 Pearson, R.G. and Dawson, T.P. (2003) Predicting the Impacts of Climate Change on the Distribution of Species: Are Bioclimate Envelope Models Useful? Global Ecology and Biogeography, 12, 361-371.
 Tjoelker, M.G., Oleksyn, J. and Reich, P.B. (1998) Seedlings of Five Boreal Tree Species Differ in Acclimation of Net Photosynthesis to Elevated CO2 and Temperature. Tree Physiology, 18, 715-726.
 Cheng, S., Dang, Q.L. and Cai, T.B. (2000) A Soil Temperature Control System for Ecological Research in Greenhouses. Journal of Forest Research, 5, 205-208.
 Scarratt, J.B. (1986) An Evaluation of Some Commercial Soluble Fertilizers for Culture of Jack Pine Container Stock. Canadian Forestry Service, Great Lakes Forestry Centre, Sault Ste. Marie, Ontario, Information Report O-X-377, 21 p.
 van den Driessche, R. (1992) Absolute and Relative Growth of Douglas-Fir Seedlings of Different Sizes. Tree Physiology, 10, 141-152.
 R Core Team (2015) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria.
 Basler, D. and Korner, C. (2012) Photoperiod Sensitivity of Bud Burst in 14 Temperate Forest Tree Species. Agricultural and Forest Meteorology, 165, 73-81.
 Heide, O.M. (1993) Daylength and Thermal Time Responses of Bud Burst during Dormancy Release in Some Northern Deciduous Trees. Physiologia Plantarum, 88, 531-540.
 Heide, O.M. (1993) Dormancy Release in Beech Buds (Fagus sylvatica) Requires Both Chilling and Long Days. Physiologia Plantarum, 89, 187-191.
 Korner, C. and Basler, D. (2010) Phenology under Global Warming. Science, 327, 1461-1462.
 Laube, J., Sparks, T.H., Estrella, N., Hofler, J., Ankerst, D.P. and Menzel, A. (2014) Chilling Outweighs Photoperiod in Preventing Precocious Spring Development. Global Change Biology, 20, 170-182.
 Myking, T. and Heide, O.M. (1995) Dormancy Release and Chilling Requirement of Buds of Latitudinal Ecotypes of Betula pendula and B. pubescens. Tree Physiology, 15, 697-704.
 Partanen, J., Koski, V. and Hanninen, H. (1998) Effects of Photoperiod and Temperature on the Timing of Bud Burst in Norway Spruce (Picea abies). Tree Physiology, 18, 811-816.
 Zohner, C.M. and Renner, S.S. (2014) Common Garden Comparison of the Leaf-Out Phenology of Woody Species from Different Native Climates, Combined with Herbarium Records Forecasts Long-Term Change. Ecology Letters, 17, 1016-1025.
 Zohner, C.M. and Renner, S.S. (2015) Perception of Photoperiod in Individual Buds of Mature Trees Regulates Leaf-Out. New Phytologist, 208, 1023-1030.
 Jach, M.E., Ceulemans, R. and Murray, M.B. (2001) Impacts of Greenhouse Gases on the Phenology of Forest Trees. In: Karnosky, D.F., Ceulemans, R., Scarascia-Mugnozza, G.E. and Innes, J.L., Eds., The Impact of Carbon Dioxide and Other Greenhouse Gases on Forest Ecosystems, CABI Publishing, CAB International, Wallingford, Oxford, 193-235.
 Slaney, M., Wallin, G., Medhurst, J. and Linder, S. (2007) Impact of Elevated Carbon Dioxide Concentration and Temperature on Bud Burst and Shoot Growth of Boreal Norway Spruce. Tree Physiology, 27, 301-312.
 Apple, M.E., Lucash, M.S., Olszyk, D.M. and Tingey, D.T. (1998) Morphogenesis of Douglas Fir Buds Is Altered at Elevated Temperature but Not at Elevated CO2. Environmental and Experimental Botany, 40, 159-172.
 Olszyk, D., Wise, C., VanEss, E., Apple, M. and Tingey, D. (1998) Phenology and Growth of Shoots, Needles, and Buds of Douglas-Fir Seedlings with Elevated CO2 and (or) Temperature. Canadian Journal of Botany, 76, 1991-2001.
 Roberntz, P. (1999) Effects of Long-Term CO2 Enrichment and Nutrient Availability in Norway Spruce. I. Phenology and Morphology of Branches. Trees, 18, 188-198.
 Bergh, J. and Linder, S. (1999) Effects of Soil Warming during Spring on Photosynthetic Recovery in Boreal Norway Spruce Stands. Global Change Biology, 5, 245-253.
 Domisch, T., Finer, L. and Lehto, T. (2001) Effects of Soil Temperature on Biomass and Carbohydrate Allocation in Scots Pine (Pinus sylvestris) Seedlings at the Beginning of the Growing Season. Tree Physiology, 21, 465-472.
 Vapaavuori, E.M., Rikala, R. and Ryyppo, A. (1992) Effects of Root Temperature on Growth and Photosynthesis in Conifer Seedlings during Shoot Elongation. Tree Physiology, 10, 217-230.
 Lyr, H. and Garbe, V. (1995) Influence of Root Temperature on Growth of Pinus sylvestris, Fagus sylvatica, Tilia cordata and Quercus robur. Trees, 9, 220-223.
 Basler, D. and Korner, C. (2014) Photoperiod and Temperature Responses of Bud Swelling and Bud Burst in Four Temperate Forest Tree Species. Tree Physiology, 34, 377-388.
 Centritto, M., Lee, H.S.J. and Jarvis, P.G. (1999) Long-Term Effects of Elevated Carbon Dioxide Concentration and Provenance on Four Clones of Sitka Spruce (Picea sitchensis). I. Plant Growth, Allocation and Ontogeny. Tree Physiology, 19, 799-806.
 Murray, M.B., Smith, R.I., Leith, I.D., Fowler, D., Lee, H.S.J., Friend, A.D. and Jarvis, P.G. (1994) Effects of Elevated CO2, Nutrition and Climatic Warming on Bud Phenology in Sitka Spruce (Picea sitchensis) and Their Impact on the Risk of Frost Damage. Tree Physiology, 14, 691-706.
 Surano, K.A., Daley, P.F., Houpis, J.L.J., Shinn, J.H., Helms, J.A., Palassou, R.J. and Costella, M.P. (1986) Growth and Physiological Responses of Pinus ponderosa Dougl ex P. Laws. to Long-Term Elevated CO2 Concentrations. Tree Physiology, 2, 243-259.
 Guehl, J.M., Picon, C., Aussenac, G. and Gross, P. (1994) Interactive Effects of Elevated CO2 and Soil Drought on Growth and Transpiration Efficiency and Its Determinants in Two European Forest Tree Species. Tree Physiology, 14, 707-724.
 Johnsen, K.H. and Seiler, J.R. (1996) Growth, Shoot Phenology and Physiology of Diverse Seed Sources of Black Spruce: I. Seedling Responses to Varied Atmospheric CO2 Concentrations and Photoperiods. Tree Physiology, 16, 367-373.
 DeLucia, E.H. and Thomas, R.B. (2000) Photosynthetic Responses to CO2 Enrichment of Four Hardwood Species in a Forest Understory. Oecologia, 122, 11-19.
 Zhang, S. and Dang, Q.L. (2006) Effects of Carbon Dioxide and Nutrition on Photosynthetic Functions of White Birch Seedlings. Tree Physiology, 26, 1457-1467.
 Pushnik, J.C., Demaree, R.S., Houpis, J.L.J., Flory, W.B., Bauer, S.M. and Anderson, P.D. (1995) The Effect of Elevated Carbon Dioxide on a Sierra-Nevadan Dominant Species: Pinus ponderosa. Journal of Biogeography, 22, 249-254.
 Ambebe, T.F., Danyagri, G. and Dang, Q.L. (2013) Low Soil Temperature Inhibits the Stimulatory Effect of Elevated [CO2] on Height and Biomass Accumulation of White Birch Seedlings Grown under Three Non-Limiting Phosphorus Conditions. Nordic Journal of Botany, 31, 239-246.
 Dawes, M.A., Hagedorn, F., Zumbrunn, T., Handa, I.T., Hattenschwiler, S., Wipf, S. and Rixen, C. (2011) Growth and Community Responses of Alpine Dwarf Shrubs to in Situ CO2 Enrichment and Soil Warming. New Phytologist, 191, 806-818.
 Walker, M.D., Wahren, C.H., Hollister, R.D., Henry, G.H.R., Ahlquist, L.E., Alatalo, J.M., Bret-Harte, M.S., Calef, M.P., Callaghan, T.V., Carroll, A.B., et al. (2006) Plant Community Responses to Experimental Warming across the Tundra Biome. Proceedings of the National Academy of Sciences of the United States of America, 103, 1342-1346.
 Butler, S.M., Melillo, J.M., Johnson, J.E., Mohan, J., Steudler, P.A., Lux, H., Burrows, E., Smith, R.M., Vario, C.L., Scott, L., Hill, T.D., Aponte, N. and Bowles, F. (2012) Soil Warming Alters Nitrogen Cycling in a New England Forest: Implications for Ecosystem Function and Structure. Oecologia, 168, 819-828.
 Melillo, J.M., Butler, S., Johnson, J., Mohan, J., Steudler, P., Lux, H., Burrows, E., Bowles, F., Smith, R., Scott, L., Vario, C., Hill, T., Burton, A., Zhou, Y.M. and Tang, J. (2011) Soil Warming, Carbon-Nitrogen Interactions, and Forest Carbon Budgets. Proceedings of the National Academy of Sciences of the United States of America, 108, 9508-9512.
 Zhang, S., Dang, Q.L. and Yu, X. (2006) Nutrient and [CO2] Elevation Had Synergistic Effects on Biomass Production but Not on Biomass Allocation of White Birch Seedlings. Forest Ecology and Management, 234, 238-244.
 Cao, B., Dang, Q.L., Yu, X. and Zhang, S. (2008) Effects of [CO2] and Nitrogen on Morphological and Biomass Traits of White Birch (Betula papyrifera) Seedlings. Forest Ecology and Management, 254, 217-224.