Elevated CO2 Did Not Stimulate Stem Growth in 11 Provenances of a Globally Important Hardwood Plantation Species

Wesolowski, Anita; Blackman, Chris J.; Smith, Renee A.; Tissue, David T.; Pfautsch, Sebastian

doi:10.3389/ffgc.2020.00066

ORIGINAL RESEARCH article

Front. For. Glob. Change, 04 June 2020

Sec. Forest Ecophysiology

Volume 3 - 2020 | https://doi.org/10.3389/ffgc.2020.00066

Elevated CO₂ Did Not Stimulate Stem Growth in 11 Provenances of a Globally Important Hardwood Plantation Species

$\nAnita Wesolowski$ Anita Wesolowski¹

Chris J. Blackman^1,2

Renee A. Smith¹

David T. Tissue¹

Sebastian Pfautsch^1,3^*

¹Hawkesbury Institute for The Environment, Western Sydney University, Richmond, NSW, Australia
²Université Clermont-Auvergne, INRA, PIAF, Clermont-Ferrand, France
³Urban Studies, School of Social Science, Western Sydney University, Parramatta, NSW, Australia

Elevated atmospheric carbon dioxide (eCO₂) often enhances rates of photosynthesis leading to increased productivity in trees. In their native habitats in Australia, eucalypts display considerable phenotypic plasticity in response to changes in environmental conditions. Little is known whether this plasticity can be harnessed effectively under future atmospheric eCO₂ conditions and be used to identify provenances with superior growth. Here, we report two experiments that assessed the physiological and growth responses of Eucalyptus grandis—one of the world's most important hardwood plantation species—to eCO₂. We used 11 provenances from contrasting climates. Our selection was based on site-specific information of long-term temperature and water availability. In Experiment 1, four provenances exhibited significant variation in light-saturated photosynthetic rates (A_sat), stomatal conductance (g_s), and concentrations of non-structural carbohydrates in leaves, stems and roots when grown under ambient CO₂ (aCO₂). Biomass of leaves, stems and roots varied significantly and were negatively correlated with mean annual temperature (MAT) at seed origin, indicating that provenances from cooler, wetter climates generally produced greater biomass. Yet, stem growth of these provenances was not stimulated by eCO₂. Given the vast environmental gradient covered by provenances of E. grandis, we expanded the selection from four to nine provenances in Experiment 2. This allowed us to validate results from Experiment 1 with its small selection and detailed measurements of various physiological parameters by focusing on growth responses to eCO₂ across a wider environmental gradient in Experiment 2. In Experiment 2, nine provenances also exhibited intraspecific differences in growth, but these were not related to climate of origin, and eCO₂ had little effect on growth traits. Growth responses under eCO₂ varied widely across provenances in both experiments, confirming phenotypic plasticity in E. grandis, though responses were not systematically correlated with climate of origin. These results indicate that selection of provenances for improved stem growth of E. grandis under future eCO₂ cannot be based solely on climate of origin, as is common practice for other planted tree species.

Introduction

Elevated atmospheric CO₂ (eCO₂) may affect tree growth by increasing photosynthesis (Way et al., 2015). Empirical and theoretical work for a range of eucalypt species has shown that eCO₂ will increase biomass production if access to water and nutrients is not limited (e.g., Wong et al., 1992; Ghannoum et al., 2010; Battaglia and Bruce, 2017; Ellsworth et al., 2017). This response is commonly observed due to the expansion of leaf area (LA), improved photosynthesis (A) and water use efficiency (WUE) under eCO₂ (McCarthy et al., 2007; Drake et al., 2011; Norby and Zak, 2011; Resco de Dios et al., 2016). However, it is also well-known that genotypes and provenances of a single tree species may differ in their response to eCO₂ (Aspinwall et al., 2015). High genetic variability across (and within) populations may result in intraspecific differences of phenotypic responses to environmental cues, commonly termed “genotype-by-environment interaction,” or G × E (Matheson and Raymond, 1986; Josephs, 2018). This intraspecific variation may be due to heritability of specific trait characteristics and due to phenotypic plasticity (Moran et al., 2016).

For decades tree breeding programs have been used to identify species and provenances that display fast stem growth (e.g., Melville, 1940; Lavigne, 1996; Pan et al., 2020) and increase productivity of commercial plantations (e.g., Cornelius, 1994; Erskine et al., 2006; Calderia et al., 2020). However, provenances selected for high rates of stem growth today may not necessarily display these favorable traits under future environmental conditions. Higher rates of growth under eCO₂ can lead to higher respiratory C loss (Dusenge et al., 2019) and metabolic activity due to higher rates of photosynthetic C assimilation and subsequent increases in non-structural carbohydrates in leaves (e.g., Moore et al., 1999; Poorter et al., 2009; Smith et al., 2012; Schmidt et al., 2017). Under such conditions, the carbon use efficiency (CUE, ratio of net primary production to gross primary production, sensu Gifford, 2003) may result in zero growth stimulation under eCO₂. If CUE is improved under eCO₂, carbon allocation patterns may change such that increased rates of C assimilation stimulate growth of roots rather than stems (Finzi et al., 2007; De Kauwe et al., 2014; Zaehle et al., 2014), which would be undesirable for the timber growing industry, yet potentially improve resilience during dry conditions (Grote et al., 2016). However, identifying provenances that show net-positive stem growth under eCO₂ is critical to meeting a rising global demand for wood and wood products (Wang et al., 2010; Booth et al., 2015; Loik et al., 2017).

Nine Eucalyptus species dominate hardwood plantations around the world due to their fast growth under a wide range of environmental conditions (Harwood, 2011; Müller da Silva et al., 2019), providing considerable capacity as sinks for atmospheric C (e.g., Binkley et al., 2017; Viera and Rodríguez-Saolliero, 2019). One of the most important eucalypt species used in plantations is E. grandis Hill ex Maid (comm. Flooded Gum or Rose Gum). Across its natural habitat this tall, fast-growing species displays considerable plasticity in tolerating a range of climates, from temperate to tropical along the coast of eastern Australia (Figure 1). Mean annual temperature (MAT) can vary from 12 to 25°C and mean annual precipitation (MAP) from 750 to 3,500 mm across its distributional range (Aspinwall et al., 2017). At least 40 seed provenances are listed for E. grandis by the Australian Tree Seed Centre and progeny trials have been operated for some time to identify seed sources of E. grandis and other eucalypt species that produce trees with favorable timber properties (Bamber and Humphreys, 1963; Wang et al., 1984; Booth, 2013).

FIGURE 1

Figure 1. Climate parameters and location for the 11 selected E. grandis sites (each provenance with a unique color). (A) Relationship between evapotranspiration/mean annual precipitation (ET/MAP) and mean annual air temperature (MAT). (B) Records of E. grandis from Atlass of Living Australia (2020) (light gray dots) and provenances location on a map of Eastern Australia (area of interested enlarged in inlet). Note that higher ET/MAP signifies greater water limitation.

High levels of plasticity in E. grandis have been demonstrated in response to nutrient addition, and water deficit (Grassi et al., 2002; Battie-Laclau et al., 2013, 2014). Drake et al. (2015) demonstrated high growth plasticity in E. grandis in response to increasing temperatures. Another recent study showed that eCO₂ increased leaf dry mass per unit area, but not respiration, in 15 provenances of E. grandis (Aspinwall et al., 2017). Yet, these studies did not report how stem growth was affected by eCO₂. In absence of this information, valuable insight can be drawn from another study that used a widely distributed eucalypt species. Eucalyptus camaldulensis was exposed to eCO₂, which stimulated biomass production, although stem biomass increased the least and root biomass the most (Aspinwall et al., 2018). Moreover, the study found no clear trend indicating that the positive growth response to eCO₂ was related to climate of origin of the different provenances tested.

Here we report two experiments that assessed phenotypic plasticity in response to growth under eCO₂ of E. grandis provenances originating from different climates. Our selection was informed by results from Drake et al. (2015) who showed that productivity of E. grandis provenances generally increased from cool to warm temperate climates when grown under aCO₂. Based on this observation, we tested if superior growth performance in provenances from warmer climates would persist under eCO₂, using two individual experiments. In Experiment 1, we examined the response to eCO₂ for a range of physiological and growth traits in four provenances of E. grandis. The provenances originated from the southern distribution cluster of the species (Figure 1) and varied in climate of origin (precipitation and temperature). Given the vast environmental gradient covered by provenances of E. grandis, we expanded the geographical range of provenances in Experiment 2, where we used nine provenances to exclusively assess the effect of eCO₂ on growth traits. Two provenances of Experiment 1 were included in Experiment 2 to allow cross-validation of results in both experiments and their relevance across wider environmental gradients.

Materials and Methods

This study reports two experiments on Eucalyptus grandis seedlings, grown in pots in a glasshouse, originating from a variety of climates at seed origin. Experiment 1 included four provenances (four climates) with five individuals per provenance in order to examine the intraspecific response of leaf physiology traits, non-structural carbohydrates and growth traits to aCO₂ and eCO₂. Following observation of intraspecific variation, we extended the selection to nine provenances in Experiment 2 to assess phenotypic plasticity in growth enhancement and correlations with climate of origin using 104 seedlings. Experiment 2 allowed validation of results from Experiment 1 by focusing on growth responses to eCO₂ across wider environmental gradients compared to more detailed physiological measurements conducted in Experiment 1. For both experiments, seed from E. grandis mother trees were collected at the 11 locations by personnel of the Australian Tree Seed Centre (Commonwealth Scientific and Industrial Research Organization, Government of Australia).

Experiment 1

Plant Material

The objective for Experiment 1 was to (a) identify if temperature or water availability at seed origin affected plant physiology and growth in aCO₂ and (b) if the observed trends remain under eCO₂. Experiment 1 lasted from March to June in 2015. Four provenances of E. grandis were selected (Table 1, Figure 1). Detailed climate-of-origin data were obtained from eMAST (www.emast.org.au) to calculate MAT, MAP and evapotranspiration (ET) between 1971 and 2010 for each location (see Aspinwall et al., 2017). Final selection for Experiment 1 was based on forming three pairs of MAT and water availability regimes (ratio of ET and MAP). Pair 1: Mt Lindsay and Coffs Harbour had the same low ET/MAP (wet) but different MAT (cool/warm). Pair 2: Coffs Harbour and Kilcoy Creek had similar MAT (warm) but different ET/MAP (wet/dry). Pair 3: Kilcoy Creek and Woondum had similar ET/MAP (dry) but different MAT (warm/hot). Seed was obtained from five mother trees (family) within each provenance of Experiment 1. Seeds were germinated under ambient atmospheric conditions and individuals grown at a plant nursery to enable collection of vegetative cuttings. Randomly selected seedlings from each provenance were coppiced to produce hedges for rooted cutting propagation. After development of sufficiently strong roots, cuttings were transported into temperature and CO₂ controlled glasshouse rooms (12 March 2015) at the Hawkesbury campus of Western Sydney University (Richmond, NSW, Australia). There, seedlings were transplanted into pots made of PVC pipe (15 cm in diameter, 40 cm high; 7 L volume) and filled with a customized soil, which consisted of 90% coarse washed river sand, 10% diatomaceous earth (intermediate grade: 0.9–2 mm from Mt Sylvia, QLD) and mineral soil (Zeogreen, premium turf grade 0.5–1.6 mm from ZeoGreen, VIC, Australia). Pots were watered daily to saturation. Liquid fertilizer (500 mL Aquasol, at 1.6 g l⁻¹; 23% N, 4% P, 18% K, 0.05% Zn, 0.06% Cu, 0.013% Mo, 0.15% Mn, 0.06% Fe, 0.011% B; Yates Australia, NSW, Australia) was added to pots every 14 days to minimize potential nutrient limitation.

TABLE 1

Table 1. Location and climate variables characterizing the origin of E. grandis seed provenance for Experiment 1 (n = 40) and Experiment 2 (n = 104); see superscript numbers.

Seedlings were grown for 14 weeks (12 March to 15 Jun 2015) under natural light in four adjacent temperature and CO₂ controlled glasshouse rooms. Temperature (T) in all glasshouse rooms was programmed to represent the 30 year average of November to May conditions in Richmond NSW and diel variation (0600–1,000 h: 22°C; 1,000–1,600 h: 26°C; 1,600–2,000 h: 24°C; 2,000–0600 h: 15°C). Ambient and elevated CO₂ concentrations were set at 400 and 640 μmol mol⁻¹ in two rooms per treatment, respectively, and controlled automatically (Lambda T, ADC BioScientific Ltd., Herts, UK). Environmental conditions inside the glasshouse bays were logged every 15 min. Ten plants of each provenance (40 seedlings total) were randomly assigned to the four glasshouse rooms and rotated between and within rooms every 14 days. Initially, 20 seedlings were grown in aCO₂ and 20 seedlings in eCO₂. At the end of the experiment, seedlings were destructively harvested. Leaf number was determined for each plant and leaf area (LA) of the canopy of each seedling was measured by placing all leaves into a leaf area meter (LI-3100C, Li-Cor Inc., NE, USA). Roots were carefully washed free of soil medium. Leaves, stems and roots were dried at 70°C to constant mass to determine dry weight.

Leaf Gas Exchange

Light-saturated photosynthetic rates of CO₂ assimilation (A_sat, μmol m⁻² s⁻¹) and stomatal conductance of water vapor (g_s, mol m⁻² s⁻¹) of fully expanded leaves were measured during the day (0930–1,530 h, sunny days) using an infrared gas analyzer (LI-6400XT, Li-Cor Inc., NE, USA). Measurement conditions for A_sat were set as: saturating light = 1600 μmol m⁻² s⁻¹; block T = 26°C; CO₂ concentration = 400 or 640 μmol mol⁻¹; and flow rate = 500 μmol s⁻¹. Leaf vapor pressure deficit (VPD) was regulated by manipulating relative humidity inside the cuvette to maintain VPD between 1.2–1.6 kPa. Leaves were allowed to equilibrate to environmental conditions inside the cuvette for at least 5 min prior to measurement. Gas analysers were always matched before logging data. A minimum of five measurements per seedling were logged at 5 s intervals and later averaged. Care was taken that leaves covered the full 2 × 3 cm measurement chamber.

Leaf δ¹³C

Relative abundance of ¹³C to ¹²C in leaf material can be used as an integrated measurement of water use efficiency (WUE) over a longer timespan (Farquhar et al., 1989). For this purpose, one fully developed leaf per seedling was collected on days of leaf gas exchange measurements during Experiment 1. Leaves were snap-frozen in liquid nitrogen, dried at 60°C for 3 days and ground to a fine powder. Only leaves that were used for gas exchange measurements were collected. Analyses of leaf bulk tissue (δ¹³C_leaf) were completed at the Australian National University using a continuous-flow stable isotope ratio mass spectrometer (Fison Isochrom CF-IRMS, former Fisons Instruments, UK). The photosynthetic discrimination rate (Δ¹³C), which relates to WUE for the hours of 0900–1600, was calculated from leaf material (δ¹³C_leaf) and known atmospheric values (δ¹³C_air), expressed in ‰ (Farquhar et al., 1982):

\begin{array}{l} Δ^{13} C = \frac{δ^{13} C_{a i r} - δ^{13} C_{l e a f}}{1 + {\frac{δ^{13}}{C}}_{leaf} 1000} & (1) \end{array}

Since glasshouse air was artificially enriched with ¹³C depleted CO₂, the δ¹³C fraction of glasshouse room air (δ¹³C_air) was also determined to calculate Δ¹³C. This was achieved by collecting measurements of air of all four glasshouse rooms separately using a Tuneable Diode Laser (TGA100; Campbell Scientific, UT, USA), repeated for each room three times during a diel cycle to account for diurnal variation of CO₂. Mean Δ¹³C was calculated for each glasshouse room. After determining Δ¹³C, WUE was calculated as:

\begin{array}{l} WUE = \frac{A}{g_{s}} = \frac{c_{a}}{1.6} (\frac{b^{′} - Δ^{13} C}{b^{′} - a}) & (2) \end{array}

where g_s is stomatal conductance of water vapor, a is the average fractionation constant (4.4‰, O'Leary, 1981) for the diffusion of gases through stomata and b' represents the carboxylation fractionation by Rubisco and PEP carboxylase (27‰, Farquhar and Richards, 1984). The factor 1.6 accounts for the greater diffusivity of water vapor vs. CO₂ in air and c_a represents the partial pressure of CO₂ in the outer atmosphere.

Mean response ratios (RR) for A_sat, g_s, and Δ¹³C between eCO₂ and aCO₂ were calculated by dividing the mean value of the parameter in eCO₂ by the mean value in aCO₂. Variance of these RR was calculated according to Hedges et al. (1999).

Non-structural Carbohydrates

Plant tissues (leaves, stem, and roots) were sampled during morning hours of the final harvest of Experiment 1 to assess concentration of non-structural carbohydrates (NSC). Tissues were microwaved for 10 s to stop metabolic activity and then dried at 70°C for 3 days. A detailed protocol for NSC extraction can be found in Tissue and Wright (1995). In short, subsamples of plant material were dried and ground to a fine powder before extracting NSC. Starch was separated from soluble sugars, lipids and amino acids using a methanol/chloroform/water solution. Water soluble sugars were partitioned from remaining cell components by phase separation. Starch was digested into ss with perchloric acid. Sugar content of both fractions was quantified colorimetrically in a spectrophotometer (DU800, Beckman Coulter Australia Pty Ltd, NSW, Australia) via a phenol-sulfuric acid reaction (Tissue and Wright, 1995).

Biomass

At harvest (15–18 June 2015, 14 weeks since repotting), leaf area (LA) of the whole canopy was measured by placing all leaves into a leaf area meter (LI-3100C, Li-Cor Inc., NE, USA); leaf size was estimated by dividing LA by leaf number. Leaf mass per area (LMA) was calculated by dividing dry weight of its leaves (DW_leaf) by LA. Measurements of height (h, cm) and diameter (d, mm; measured 2 cm above soil surface) commenced at week 2 since repotting into the glasshouse and were subsequently conducted on a weekly basis. Stem volume (v_stem, cm³) was estimated using h and d, assuming the shape of a cone using the following formula:

\begin{array}{l} v_{s t e m} = \frac{1}{3} π r^{2} h & (3) \end{array}

Total plant dry mass (DW_plant) was calculated by summing DW_leaf, DW_stem, and DW_root. Biomass fraction was calculated by dividing DW of plant organs by DW_plant as leaf-mass-fraction (LMF), stem-mass-fraction (SMF) and root-mass-fraction (RMF).

Data Analyses

Data were analyzed using R-3.6.1 software (R Core Team, 2019). Analysis of variance (ANOVA) was used to test if the [CO₂] treatment as well as the origin of the four provenances had a significant effect on leaf gas exchange, growth-related traits (DW, v_stem, LMF, SMF, RMF) leaf specific traits (leaf number, leaf size, LA, LMA) and NSC. Statistical models were tested with all interactions between factors ([CO₂] and provenance). Assumptions of data normality and equal variance were tested using Shapiro-Wilk and Levene's test, as well as a quantile-quantile plot for visual assessment. The “car” package (Fox and Weisberg, 2019) was used for Anova() and leveneTest() functions with Type III sums of squares and Levene's test for equal variance among test groups, respectively. Post-hoc tests were conducted when plant traits were significantly different across provenances. To establish whether A_sat or g_s was the driver for changes in Δ¹³C correlations were analyzed with Pearson's correlation coefficient using the cor.test() function of the base package. All means are presented with one standard error (SE).

Experiment 2

Plant Material

The objective of Experiment 2 was to validate results of Experiment 1 recorded for a narrow climate of seed origin gradient and be able to identify a generalized growth response strategy of E. grandis provenances. Experiment 2 lasted from February to May 2017. The range of climate at seed origin was extended to nine provenances in Experiment 2 including two provenances from Experiment 1 (Table 1, Figure 1). Detailed selection criteria for the provenances used in Experiment 2 can be found in Aspinwall et al. (2017). Seed was obtained from 2 to 5 mother trees within each provenance. Seeds were germinated as described for Experiment 1 and were transported into temperature and CO₂ controlled glasshouse rooms on 27 February 2017 using similar pots and soil as in Experiment 1.

Seedlings of Experiment 2 were grown at the same location as Experiment 1 for 10 weeks (27 February to 11 May 2017). Glasshouse environmental conditions were comparable to Experiment 1. Experiment 2 included 3–7 plants per provenance, of which 52 were randomly assigned and grown in aCO₂ and 52 plants in eCO₂. Seedlings were destructively harvested in the same way as Experiment 1.

Biomass

Harvest was conducted on 10 and 11 May 2017, 10 weeks since repotting. The same measurements were taken as in Experiment 1. Measurements of h and d commenced during the first week following replanting into the glasshouse.

Data Analyses

Data were analyzed using R-3.6.1 software (R Core Team, 2019). Analysis of variance (ANOVA) was used to test if the [CO₂] treatment as well as climate origin of the nine provenances had a significant effect on growth-related traits (DW, v_stem) and leaf specific traits (LA, LMA). See Experiment 1 for further details. Correlations between DW_shoot and DW_stem, respectively, and climate variables (MAT, MAP, ET/MAP) were analyzed with Pearson's correlation coefficient using the cor.test() function of the base package. The same approach was used for determining the correlation between DW_stem and LMA. All figures were drawn using Python (Phyton Software Foundation)¹.

All figures were plotted using Python (Phyton Software Foundation).

Results

Experiment 1

Environmental Conditions

Throughout Experiment 1, relative humidity ranged from 57–85%, T averaged 26.9 ± 0.1(±1SE) °C during the photoperiod, and excellent [CO₂] control was achieved (aCO₂: 401 ± 0.2 μmol mol⁻¹; eCO₂: 645 ± 0.1 μmol mol⁻¹). From March to June 2015, mean photon flux density during mid-day was 950 μmol m⁻² s⁻¹, ranging from 115 μmol m⁻² s⁻¹ on a cloudy day to 1,900 μmol m⁻² s⁻¹ on a sunny day.

Physiological Responses to eCO₂

A_sat differed significantly by provenance and [CO₂] treatment. However, there were no provenance by [CO₂] interactions statistically confirming these observed trends. In aCO₂, A_sat of seedlings originating from the warmest climate (provenance Woondum) was significantly lower (18.2 ± 1.4 μmol m⁻² s⁻¹, p = 0.002; Figure 2) than A_sat of seedlings in the other three provenances, which were similar to each other (22.1 ± 0.6 μmol m⁻² s⁻¹). Across all provenances, A_sat was on average 23% greater in eCO₂ (23.3 ± 0.9 μmol m⁻² s⁻¹) compared to those grown in aCO₂ (19.0 ± 0.6 μmol m⁻² s⁻¹, p = 0.003). In eCO₂, A_sat of seedlings from Woondum was 26% lower than seedlings originating from cooler climates (Figure 2). Notably, seedlings from the coolest climate (Mt Lindsay) showed on average the largest response to eCO₂, increasing A_sat by 40% in eCO₂ compared to aCO₂, yet also displayed high variability within this provenance (Figure 3).

FIGURE 2

Figure 2. Light-saturated photosynthesis (A_sat; A,B) stomatal conductance (g_s; C,D) and carbon isotope discrimination (Δ¹³C; E,F) of four provenances of E. grandis selected for Experiment 1. Panels on the left show seedlings grown in ambient [CO₂] (aCO₂), panels on the right show seedlings grown in elevated [CO₂] (eCO₂; n = 5 per provenance). Error bars represent 1SE. Different capital case letters represent statistical significance of traits among provenances within the same [CO₂] treatment. Different lower-case letters represent statistical significance of traits between [CO₂] treatments per provenance.

FIGURE 3

Figure 3. Response ratios of leaf physiological traits calculated for Experiment 1. Response ratios represent differences in four different provenances of E. grandis grown under eCO₂ compared to aCO₂ (n = 5 per provenance and [CO₂] treatment). Traits: light-saturated photosynthetic carbon assimilation (A_sat; A), stomatal conductance (g_s; B), carbon isotope discrimination (Δ¹³C; C). Error bars represent 1SE.

[CO₂] treatment did not affect g_s, but values differed across provenances. In May 2015, seedlings growing in aCO₂ and originating from provenance Woondum had 49% lower g_s (0.45 ± 0.11 mol m⁻² s⁻¹) than seedlings from provenance Coffs Harbour, averaging 0.88 ± 0.07 mol m⁻² s⁻¹ (p = 0.001; Figures 2, 3). Seedlings from provenances Mt Lindsay and Kilcoy Creek had intermediate mean values of 0.69 ± 0.05 mol m⁻² s⁻¹ and 0.75 ± 0.09 mol m⁻² s⁻¹, respectively. In eCO₂, seedlings from provenance Woondum had similar low g_s averaging 0.44 ± 0.12 mol m⁻² s⁻¹ while seedlings from provenance Mt Lindsay had the highest values with 0.93 ± 0.06 mol m⁻² s⁻¹ (35% increase from aCO₂ to eCO₂).

Provenances did not differ in Δ¹³C when grown in aCO₂. Seedlings grown in eCO₂ had an improved WUE, indicated by Δ¹³C being significantly lower compared to seedlings grown in aCO₂ (22.5 ± 0.3o vs. 25.2 ± 0.3o, p < 0.001; Figures 2, 3). A_sat was negatively correlated with Δ¹³C (Pearson's r = −0.59, p > 0.001) whereas g_s was not (p = 0.99).

Non-structural Carbohydrates

eCO₂ did not change concentrations of soluble sugars or starch in any provenance or tissue. For instance, soluble sugars in leaves averaged 55.0 ± 2.1 mg g⁻¹ dry weight in aCO₂ and 53.8 ± 1.6 mg g⁻¹ dry weight in eCO₂. Leaf starch averaged 68.0 ± 1.7 mg g⁻¹ dry weight in aCO₂ and 64.5 ± 1.4 mg g⁻¹ dry weight in eCO₂. Data were pooled for [CO₂] treatments to assess variation of NSC concentration in plant organs among provenances. Concentrations of soluble sugars at harvest differed significantly in leaves, stem and roots among provenances. Seedlings originating from Woondum had the greatest concentrations of soluble sugars in leaves (60.1 ± 2.5 mg g⁻¹ dry weight, p = 0.02) and stem tissue (11.4 ± 1.0 mg g⁻¹ dry weight, p < 0.001), respectively. Concentrations in leaves were 21% greater than soluble sugars in seedlings from Coffs Harbour and 64% greater than soluble sugars in stem tissue of Kilcoy Creek. In roots, concentration of soluble sugars differed significantly among provenances (p = 0.02), where seedlings originating from Coffs Harbour contained the greatest amount of soluble sugars (18.3 ± 0.9 mg g⁻¹ dry weight) followed by those from Woondum (16.7 ± 1.8 mg g⁻¹ dry weight). Concentrations of starch at harvest only differed among provenances in leaves (p = 0.05; Figure 4). Seedlings from Kilcoy Creek had the greatest starch concentrations in leaves (69.7 ± 2.0 mg g⁻¹ dry weight), whereas those from Woondum had the lowest (61.3 ± 2.0 mg g⁻¹ dry weight; Figure 4).

FIGURE 4

Figure 4. Non-structural carbohydrates (NSC) in three tissues of E. grandis. Samples were collected at the final harvest. Data are shown for ambient [CO₂] treatment (aCO₂; left; n = 5) and elevated [CO₂] treatment (eCO₂; right; n = 5). NSC are presented for two fractions: soluble sugars (ss; light colored area) and starch (st; full colored area) in leaves (NSC_leaf; A,B), stem (NSC_stem; C,D) and roots (NSC_root; E,F). Error bars represent 1SE. Colors and abbreviations represent different provenances. Different capital case letters represent statistical significance of ss or st, respectively, among provenances within the same [CO₂] treatment.

Growth

At harvest, seedlings were 62–118 cm tall. Throughout the experiment, stem volume (v_stem) was not significant across [CO₂] treatments (p = 0.15) nor [CO₂] x provenance interaction (p = 0.91) but varied across provenances (p < 0.001; Figures 5A,B). Specifically, seedlings originating from the warmest climate (Woondum) had continuously lower stem volume than seedlings from the coolest climate (Mt Lindsay). Trends were detected for relative responses of v_stem from aCO₂ to eCO₂ with seedlings from Coffs Harbour indicated the greatest increase (+33%). At harvest, DW_plant of individual seedlings varied from 7.3 to 59.5 g. Leaf, stem and root biomass was only significant across provenances, but not [CO₂] nor their interactions (Table 7). Under aCO₂, seedlings from Woondum had the lowest DW_leaf (p = 0.004), DW_stem (p < 0.001; Figure 6A), and DW_root (p < 0.001) followed by Coffs Harbour. Seedlings from Mt Lindsay had the greatest DW_plant (39.8 ± 3 g; Table 3; Figure 6A). DW_root was essentially statistically significant (p = 0.07) increasing from 4.9 ± 0.6 g in aCO₂ to 6.3 ± 0.8 g in eCO₂ (Table 3). Responses of DW_plant from aCO₂ to eCO₂ widely differed across provenance for biomass of all plant tissues (DW_stem as example shown in Figures 6B, 7A). As a trend, seedlings from Coffs Harbour increased mass in all tissues resulting in a mean DW_plant of 34.7 ± 4.4 g (+44%). DW_stem was correlated to climate parameters MAT, MAP, ET/MAP, mean annual minimum and maximum temperature. The best fit was found for MAT (r² = 0.38, p = 0.004; Figures 8A,B) where biomass decreased with greater MAT. In aCO₂, this linear fit had a negative slope of −1.54, whereas the slope in eCO₂ was slightly steeper with −1.74. Biomass partitioning calculated as LMF, SMF and RMF did not differ among provenances, but LMF was statistically significant and RMF essentially significant, respectively between [CO₂] treatments (Table 7). In aCO₂, LMF averaged 0.40 ± 0.008 compared to 0.37 ± 0.009 in eCO₂ (p = 0.01; Table 3). RMF showed an opposite trend increasing from 0.17 ± 0.006 in aCO₂ to 0.19 ± 0.01 in eCO₂ (p = 0.09; Table 3).

FIGURE 5

Figure 5. Stem volume of E. grandis during Experiment 1 and 2. Seedlings (n = 3–7 per provenance) were grown in ambient (aCO₂; left) and ambient [CO₂] (eCO₂; right). Colors in (A,B) represent the four selected provenances of Experiment 1. The nine selected provenances of Experiment 2 are shown in (C,D). Shaded bands represent 1SE.

FIGURE 6

Figure 6. Mean stem dry weight (g) of E. grandis at the end of Experiment 1 and 2. Seedlings (n = 3–7 per provenance) were grown in ambient (aCO₂; left) and elevated CO₂ (eCO₂; right). Colors in (A,B) represent the four selected provenances of Experiment 1. The nine selected provenances of Experiment 2 are shown in (C,D). Error bars represent 1SE.

FIGURE 7

Figure 7. Response ratios of stem dry weight at harvest (DW_stem) for Experiment 1 and 2. Response ratios represent differences in the different provenances of E. grandis grown under elevated [CO₂] (B) compared to ambient [CO₂] (A; n = 3–7 per provenance and [CO₂] treatment). Colors in (A) represent the four selected provenances of Experiment 1 and colors in (B) represent the nine selected provenances of Experiment 2. Error bars represent variance.

FIGURE 8

Figure 8. Correlation between stem dry weight at harvest (DW_stem) and mean annual temperature (MAT) at seed origin for Experiment 1 and 2 for ambient (aCO₂; A,B) and elevated [CO₂] (eCO₂; n = 3–7 per provenance and [CO₂] treatment; C,D). Circles in (A,B) represent data of the four provenances in Experiment 1, squares in (C,D) represent data of the nine provenances in Experiment 2. Dark blue line shows linear fit of the data with light blue areas representing confidence intervals.

At final harvest, leaf number and LA differed significantly among the four provenances, while leaf size and LMA did not. Leaf traits did not significantly respond to [CO₂] nor [CO₂] x provenance interaction. Under aCO₂, seedlings from Mt Lindsay had developed the most leaves (236 ± 24 plant⁻¹, p = 0.01) and had also the largest leaf area (4249 ± 143 cm² plant⁻¹, p = 0.003; Table 2). Seedlings from provenance Woondum had developed about half of the number and area of leaves compared to those from Mt Lindsay. Seedlings from Kilcoy Creek had the largest leaf size (19.0 ± 3.0 cm²) and those from Coffs Harbour the smallest (15.1 ± 1.9 cm²; Table 2). LMA under aCO₂ was similar among provenances, averaging 36.5 to 40.0 g m⁻². Although a trend toward increased leaf size and area was apparent in all seedlings when grown in eCO₂, large intra-specific variation resulted in non-significant changes compared to seedlings grown in aCO₂.

TABLE 2

Table 2. Leaf characteristics of E. grandis seedlings of Experiment 1 grown in ambient [aCO₂; (a)] and elevated [CO₂] [eCO₂; (b)].

TABLE 3

Table 3. Dry mass (g) and biomass fractions of E. grandis at the end of Experiment 1.

Experiment 2

Environmental Conditions

Relative humidity during Experiment 2 ranged from 62 to 88% between the hours of 10:00 and 16:00. T averaged 27.6 ± 0.2°C and mean [CO₂] were for aCO₂ 396 ± 0.2 μmol mol⁻¹ and for eCO₂ 643 ± 0.2 μmol mol⁻¹.

Growth

At harvest, seedlings were 58–109 cm tall. Throughout Experiment 2, v_stem was not significant across [CO₂] treatments nor [CO₂] x provenance interaction yet varied across provenances (p < 0.001; Figures 5C,D). Seedlings originating from provenance Bellthorpe had continuously greater stem volume than seedlings from provenance Mt Mee (Figure 5). Leaf and stem biomass were not significant across [CO₂] treatments (p =) nor [CO₂] x provenance interaction but differed among the nine provenances (Table 5). Under aCO₂, seedlings from provenance Mt Mee of the warmest and wettest climate had the lowest DW_leaf (p = 0.008) and DW_stem (p = 0.01; Figure 6C). Seedlings from provenance Connondale had the greatest DW_shoot (67.3 ± 6.9 g; Table 5), 91% greater than seedlings from Mt Mee. Similar to Experiment 1, trends for intraspecific variability in growth response to eCO₂ were noticeable (Table 5; Figures 6D, 7B). Seedlings from provenance Wedding Bells SF, similar to seedlings from Coffs Harbour in Experiment 1, had the greatest increase of DW_shoot resulting in a mean of 49.1 ± 7.1 g (+31% from aCO₂ to eCO₂). Yet, dry weight of shoot and stem was, contrary to Experiment 1, not correlated to any climate parameters of seed origin in aCO₂ or eCO₂ (Table 6; Figures 8C,D).

LA and LMA of the nine selected provenances only differed significantly among provenances at the end of Experiment 2 (Table 4). Under aCO₂, seedlings from Connondale originating from warmer and wetter climate had the largest LA (8212 ± 779 cm², p = 0.002; Table 4). Seedlings from provenance Mt Mee, originating from the wettest climate, had only half the LA compared to Connondale. LMA under aCO₂ differed among the nine provenances. Seedlings from provenance Richmond Range with relatively wet climate had the greatest LMA averaging 55.9 ± 3.0 g m⁻². Similar to Experiment 1, seedlings from provenance Mt Lindsay had the lowest LMA (41.9 ± 1.2 g m⁻²). Trends for intraspecific variation in the response of leaf traits to eCO₂ were noticeable. Six of the nine provenances had increased LA under eCO₂ with the largest increase in seedlings from provenance Wedding Bells SF (+28%; Table 4), resembling similar climate to seedlings from provenance Coffs Harbour in Experiment 1. In aCO₂, LA was negatively correlated to MAT (Pearson's r = −0.27, p = 0.05) and ET of seed origin (Pearson's r = −0.27, p = 0.05; Table 6). LMA was positively correlated to DW_stem (Pearson's r = 0.29, p = 0.04 in aCO₂ and Pearson's r = 0.44, p = 0.001 in eCO₂, respectively).

TABLE 4

Table 4. Leaf characteristics of E. grandis seedlings of Experiment 2 grown in ambient [aCO₂; (a)] and elevated [CO₂] [eCO₂; (b)].

TABLE 5

Table 5. Dry mass (g) of E. grandis at the end of Experiment 2.

TABLE 6

Table 6. Pearson's correlation coefficient (Pearson's r) between leaf area (LA), shoot dry weight (DW_shoot) and stem dry weight (DW_stem), and climate variables of seed origin of the nine provenances selected for Experiment 2: mean annual temperature (MAT, °C), difference between mean annual maximum and minimum temperature (ΔT, °C), mean annual precipitation (MAP, mm), mean annual evapotranspiration (ET, mm), and evapotranspiration/MAP (ET/MAP).

TABLE 7

Table 7. Statistical results with degrees of freedom (df), f-values (f), and p-values (p) for Analysis of Variance (ANOVA) for E. grandis.

Discussion

Two experiments on E. grandis were used to investigate phenotypic plasticity in this important hardwood plantation species. These experiments revealed significant variability regarding leaf physiology traits, concentrations of non-structural carbohydrates and growth parameters among provenances grown in aCO₂. Yet, contrary to our expectation, variation in biomass was not systematically related to climate at seed origin in E. grandis seedlings.

Although eCO₂ increased A_sat and WUE (calculated from leaf carbon isotope discrimination and reflecting the lifespan of the leaf), this did not result in a systematic stimulation of growth and increased biomass among provenances. Systematic stimulation of growth by eCO₂ was not apparent, possibly due to the limited number of selected provenances and high intraspecific plant-to-plant variation. However, even by increasing the number of tested provenances to improve statistical power in Experiment 2, we did not observe a uniform growth stimulation and biomass (stem, leaves, roots) increase in eCO₂ nor a significant provenance by CO₂ effect. However, some provenances did show greater responses of physiology (A_sat and g_s) or growth in eCO₂ than others, potentially related to environmental conditions of seed origin—as seen with provenances Coffs Harbour and Wedding Bells SF with strong positive biomass increases under eCO₂. Yet, unidirectional phenotypic plasticity among provenances as result of eCO₂ was not observed due to great variability within groups.

Physiological Differences Among Provenances

Greater carbon assimilation, due to greater A_sat and a larger canopy, enabled seedlings from regions with low ET/MAP to also produce more stem biomass (i.e., provenance Coffs Harbour) relative to seedlings from drier climates (i.e., provenance Woondum and Mt Mee). Variability in LMA among provenances can also influence plant growth (Poorter et al., 2009) as shown with our LMA data correlating with stem biomass. As LMA has a strong positive relationship with leaf density, our results suggest that provenances MT and RR have greater leaf density than provenances KC and ML. Higher leaf density may be caused by a greater proportion of mesophyll tissue or lignified tissue resulting in greater leaf toughness and consequently plant survival (Poorter et al., 2009). Despite these observed trends, stem biomass of E. grandis was not correlated to climate of origin in Experiment 2. We thus had to reject the notion that results from Experiment 1 could be used to extrapolate growth responses to eCO₂ across wider environmental gradients. However, intraspecific differences across provenances or families are frequently reported for trees in response to their local climate conditions. European pine and beech species (Correia et al., 2008; Rose et al., 2009; Bachofen et al., 2018) as well as North American poplars (Kaluthota et al., 2015) as well as eucalypts (Blackman et al., 2016; Aspinwall et al., 2018) have been shown to display considerable intraspecific variation in physiological and growth-related traits when grown as potted plants or in common gardens. These variations were correlated to temperature, water availability and the length of the growing season at seed origin. Seedlings from provenances originating either from areas with wet and cool or wet and warm climate produced the greatest biomass when grown at aCO₂. In support of this trend of adaptation to local climate conditions, biomass production in E. grandis was greater in trees from cooler and wetter climates (Drake et al., 2015). Thus, E. grandis does show signs of a systematic adaption to environmental conditions at seed origin, yet climate conditions alone are insufficient to explain observed responses for every provenance.

Physiological Responses of Provenances to eCO₂

Intraspecific variation across provenances was more prominent than responses to eCO_2. Hence, ranking among provenances for physiology and growth traits persisted, supporting Hypothesis 2. A_sat increased by 24% on average, similar to the stimulation of A_sat under eCO₂ in other trees (e.g., Ainsworth and Rogers, 2007; Ghannoum et al., 2010; Leakey et al., 2012; Resco de Dios et al., 2016), which predominately results from increased rates of carboxylation by Rubisco and reduced rates of photorespiration (Drake et al., 1997; Long et al., 2004). Additionally, eCO₂ increased WUE in all provenances (i.e., Δ¹³C declined), an effect that has been attributed to decreases in g_s (Saxe et al., 1998; Medlyn et al., 2001) and subsequent decline in conductance of water vapor by leaves (e.g., Ainsworth and Rogers, 2007; Adair et al., 2011; Keenan et al., 2013; Reich et al., 2014; Frank et al., 2015; Dekker et al., 2016). Yet, in E. grandis, improved WUE was due to increasing carbon uptake (increases in A_sat) under non-limiting water supply rather than decreases in g_s, a strategy also confirmed for Fagus sylvatica (Aranda et al., 2017). Under eCO₂, g_s is commonly reported to decline (e.g., Ainsworth and Long, 2005). However, a lack of downregulation of g_s under eCO₂ was also found in Larix and Pinus trees (Streit et al., 2014). Furthermore, the strong influence of leaf tissue structure on gas exchange in E. camaldulensis (Blackman et al., 2016) may explain why g_s response to eCO₂ was not significant given that LMA was also not affected in E. grandis. Regrettably, we did not investigate stomatal traits, which might have helped differentiate anatomical from biochemical constraints on WUE (Blackman et al., 2016). Limited information suggests changes in stomatal density of broadleaved trees may occur as an interactive response to [CO₂], water availability and origin of provenance (Pyakurel and Wang, 2014). In agreement with the current study, most leaf physiological traits in 14 genotypes within six provenances of E. camaldulensis showed extensive intraspecific variation, but a limited growth response to eCO₂ (Blackman et al., 2016).

When trees are grown under eCO₂, they often exhibit increased carbon assimilation rates. However, these increases sometimes are not reflected in simultaneous increases in stem biomass, indicating that additional carbon gained is not necessarily allocated to stem growth (Resco de Dios et al., 2016). For example, net carbon assimilation of E. camaldulensis was greatly stimulated by eCO₂, yet no increase in stem biomass or LMA (leaf structure) was detected (Blackman et al., 2016). In the present study, concentration of starch in leaves of E. grandis was greater in eCO₂ compared to aCO₂, suggesting that leaves may serve as temporary storage organs under eCO₂ rather than building new leaf structure (e.g., larger or thicker leaves).

Another possible fate of additional carbon gained in eCO₂ can be its loss through elevated rates of cell respiration and/or root exudation. We did not examine either of the two processes, but studies have shown that large variability exists in the effects of eCO₂ on plant mitochondrial respiration rates (see Dusenge et al., 2019 for review). Studies have found respiration to remain stable, increase or decline when plants are exposed to eCO₂ (Curtis, 1996; Wang et al., 2001; Tissue et al., 2002; Ayub et al., 2011; Crous et al., 2012; Li et al., 2013; Gauthier et al., 2014). Increased rates of respiration would explain the absence of greater allocation of assimilated carbon into aboveground biomass in E. grandis. However, higher respiration rates have also been linked to greater carbohydrate accumulation in pine (Tjoelker et al., 2009; Li et al., 2013), but it was not observed here in E. grandis.

Growth is often enhanced in other species of Eucalyptus grown under eCO₂ (Ghannoum et al., 2010; Ayub et al., 2011; Smith et al., 2012; Duan et al., 2013; Aspinwall et al., 2018), but was not statistically significant here for stem and leaf biomass of E. grandis. However, the lack of significant aboveground growth stimulation in eCO₂ has been observed in other species of Eucalpytus (Duff et al., 1994) and Populus (Liberloo et al., 2005) despite increased carbon assimilation (Cantin et al., 1997; Atwell et al., 2007; Lawson et al., 2017; Bachofen et al., 2018; Killi et al., 2018). It has also been observed that carbon may be allocated belowground to facilitate increased nutrient uptake when photosynthesis is increased in eCO₂ (Hättenschwiler and Körner, 1997; Iverson et al., 2008; Pritchard et al., 2008; Souza et al., 2016; Aspinwall et al., 2018). Here, 13% more biomass was on average allocated to roots rather than leaves (−8%) in Experiment 1. Thus, biomass responses of E. grandis to eCO₂ may be due to carbon allocation to roots rather than aboveground biomass allocation.

Environmental Drivers

Cooler climate and low ET/MAP at seed origin of trees from Mt Linsday may represent contributing factors to the pronounced increase of A_sat in the different provenances when grown in eCO₂. A common garden study on European Picea abies found that populations from cooler climates at high-altitudes had higher leaf nitrogen concentrations, greater A_max and higher concentrations of chlorophyll than those from lower elevations (Oleksyn et al., 1998). In E. camaldulensis, genotypes with the largest growth response to eCO₂ had large increases in root mass fraction and photosynthetic nitrogen-use efficiency (Aspinwall et al., 2018). The Coffs Harbour (Experiment 1) and Wedding Bells (Experiment 2) provenances had the largest relative increase of leaf area and biomass when grown in eCO₂. Both locations are in very close proximity to each other and experience a similarly low ET/MAP compared to Mt Lindsey. Shifting such species from cooler to warmer growth conditions could influence temperature optima for photosynthesis. This physiological optimization due to warming while maintaining optimal access to water is in agreement with previous temperature-shift experiments using E. grandis (Drake et al., 2015).

In addition to climate, soil characteristics are known to influence productivity of eucalypt trees (Cavalli et al., 2020). Eucalyptus grandis, however, tolerates a wide range of soil conditions, reflected in high growth rates of provenances in their native habitats in Australia as well as plantations around the world. Here we controlled the influence of soil on growth by providing all provenances with the same customized, non-nutrient limited soil. The impact of eCO₂ on stem growth of individual provenances, when grown in soil from their native habitat, should be investigated in future glasshouse experiments.

Conclusions

It now seems well-established that tree genotypes and provenances respond differentially to eCO₂ (Ceulemans et al., 1996; Dickson et al., 1998; Isebrands et al., 2001; Mohan et al., 2004). While individual species may show increased stem growth and biomass accumulation under eCO₂, this is not a widely observed trend (Resco de Dios et al., 2016). Here we added one of the most important hardwood plantation species to the growing list of tree species that seem to accelerate photosynthetic carbon assimilation under eCO₂ yet do not use surplus photosynthates to increase stem growth. Selection of high performing tree provenances for the next generation of hardwood plantations in a higher CO₂ world must account for this effect to remain economically profitable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

AW, SP, and DT designed Experiment. AW, CB, DT, and SP designed Experiment. RS managed the controlled environment facilities, including CO₂ treatments. AW analyzed the data and wrote the complete first draft. All authors collected plant physiological measurements, harvested data, and assisted in developing the final version of the manuscript.

Funding

This study was supported by the Science and Industry Endowment Fund (SIEF grant RP04-122) and the Higher Degree Research Scholarship from the Hawkesbury Institute for the Environment.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank Carrie Drake, Wen Shi, and Craig Barton for technical support and helpful advice for development of this manuscript by Mark Tjoelker and Victor Resco de Dios is acknowledged.

Footnotes

1. ^Python Software Foundation. Python Language Reference, version 3.6. Available online at http://www.python.org.

References

Adair, E. C., Reich, P. B., Trost, J., and Hobbie, S. E. (2011). Elevated CO₂ stimulates grassland soil respiration by increasing carbon inputs rather than by enhancing soil moisture. Glob. Change Biol. 17, 3546–3563. doi: 10.1111/j.1365-2486.2011.02484.x

ORIGINAL RESEARCH article

Elevated CO2 Did Not Stimulate Stem Growth in 11 Provenances of a Globally Important Hardwood Plantation Species

Introduction

Materials and Methods

Experiment 1

Plant Material

Leaf Gas Exchange

Leaf δ13C

Non-structural Carbohydrates

Biomass

Data Analyses

Experiment 2

Plant Material

Biomass

Data Analyses

Results

Experiment 1

Environmental Conditions

Physiological Responses to eCO2

Non-structural Carbohydrates

Growth

Experiment 2

Environmental Conditions

Growth

Discussion

Physiological Differences Among Provenances

Physiological Responses of Provenances to eCO2

Environmental Drivers

Conclusions

Data Availability Statement

Author Contributions

Funding

Conflict of Interest

Acknowledgments

Footnotes

References

Elevated CO₂ Did Not Stimulate Stem Growth in 11 Provenances of a Globally Important Hardwood Plantation Species

Leaf δ¹³C

Physiological Responses to eCO₂

Physiological Responses of Provenances to eCO₂