Thermal acclimation of leaf photosynthetic traits in an evergreen woodland, consistent with the coordination hypothesis

. Ecosystem models commonly assume that key photosynthetic traits, such as carboxylation capacity measured at a standard temperature, are constant in time. The temperature responses of modelled photosynthetic or respiratory rates then depend entirely on enzyme kinetics. Optimality considerations, however, suggest this assumption may be incorrect. The “coordination hypothesis” (that Rubisco- and electron-transport-limited rates of photosynthesis are co-limiting under typical daytime conditions) predicts, instead, that carboxylation ( V cmax ) capacity should acclimate so that it increases somewhat with growth temperature but less steeply than its instantaneous response, implying that V cmax when normalized to a standard temperature (e.g. 25 ◦ C) should decline with growth temperature. With additional assumptions, similar predictions can be made for electron-transport capacity ( J max ) and mitochondrial respiration in the dark ( R dark ) . To explore these hypotheses, photosynthetic measurements were carried out on woody species during the and the cool in the semi-arid Great Western Woodlands, Australia, under broadly similar light environments. A consistent proportionality between max was ,


Introduction
Net photosynthetic CO 2 uptake (A net ) depends on temperature in all vegetation models, but models commonly disregard possible acclimation of the parameters determining A net to temporal variations in the growth environment. There are plenty of data on the instantaneous (minutes to hours) temperature responses of photosynthetic uptake (Hikosaka et al., 2006;Sage and Kubien, 2007;Way and Sage, 2008), but data on the responses of photosynthetic traits on ecologically relevant timescales (days to years) are scarce (Wilson et al., 2000;Lin et al., 2013). Nonetheless there is evidence that temperature responses of biochemical processes are a function of plant growth temperature, and not just instantaneous temperature, from studies comparing species (Miyazawa and Kikuzawa, 2006;Kattge and Knorr, 2007), to experiments or serial measurements on single species (Medlyn et al., 2002b;Onoda et al., 2005). Neglecting acclimation to the growth temperature could lead to incorrect model estimates of the responses of primary production and ecosystem carbon budgets to climate change. Farquhar et al. (1980) provided the standard model to predict photosynthetic responses to the environment in C 3 plants. The model describes photosynthesis as instantaneously determined by the slower of two biochemical rates: the carboxylation of RuBP (ribulose-1,5-bisphosphate), which is dependent on Rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase) activity (V cmax ), and electron transport for RuBP regeneration in the Calvin cycle, which is dependent on light intensity and (at high light intensity) limited by the capacity of the electron transport chain (J max ). Both rates are influenced by intercellular CO 2 concentration (c i ), which in turn is partially regulated by stomatal conductance (g s ). Unlimited mesophyll conductance (Miyazawa and Kikuzawa, 2006;Lin et al., 2013) remains the most common implementation of the Farquhar model: apparent V cmax and J max values can be calculated from A-c i curves under this assumption provided K C and K O (the Michaelis-Menten coefficients for carboxylation and oxygenation, respectively) are appropriately specified. V cmax and J max are highly variable both within and between species (Wullschleger, 1993), ranging through 2 to 3 orders of magnitude. Despite the accepted importance of these parameters for predicting rates of net CO 2 exchange in natural ecosystems, a full understanding of how changes in the environment affect these quantities is lacking.
The "coordination hypothesis" (Chen et al., 1993;Field and Mooney, 1986;Maire et al., 2012) predicts that values of V cmax and J max should acclimate, in time as well as in space, in such a way that carboxylation and RuBP regeneration are co-limiting under average daytime conditions. This has been suggested often in the literature, both as a prediction based on optimality considerations (e.g. Von Caemmerer and Farquhar, 1981;Dewar, 1996;Haxeltine and Prentice, 1996) and as an empirical observation, consistent with the finding that carboxylation is a limiting factor at saturating light intensity as has been shown to apply in the great majority of cases (De Kauwe et al., 2015). The coordination hypothesis makes a number of testable predictions regarding photosynthetic properties of woody plants experiencing large seasonal variations in growth temperature, including an increase in V cmax with rising growth temperature because at higher daily temperatures, somewhat greater Rubisco activity is required to match any given rate of photosynthesis. On the other hand, it predicts that the value of V cmax normalized to 25 • C -and with it, the leaf N per unit area (Maire et al., 2012;Dong et al., 2017) -should decline with increasing daily average temperature, because the quantities of proteins needed to maintain a given level of photosynthetic activity decline more steeply with temperature than does the predicted increase in V cmax . These predictions were strongly supported by experiments on tree species grown at different temperatures (Scafaro et al., 2017). Scafaro et al. (2017) also demonstrated that the acclimation of V cmax involves changes in the Rubisco amount and the relative allocation of leaf N to Rubisco. If it is further assumed that J max covaries closely with V cmax , as is commonly observed, then similar predictions apply to the acclimation of J max . Wang et al. (2017) have provided a fuller treatment of J max acclimation, based on a cost-minimizing principle that also correctly predicts the dependence of the ratio J max /V cmax on growth temperature (Kattge and Knorr, 2007;Wang et al., 2017). Note that the coordination hypotheses as presented here is the strong form of the hypothesis, which allows the metabolic component of leaf N to be optimized -thus differing from an alternative interpretation, in which leaf N remains constant and only the J max /V cmax ratio is allowed to vary (Medlyn, 1996;Ali et al., 2016;Quebbeman and Ramirez, 2016).
R dark (leaf respiratory CO 2 release in darkness) is known to acclimate to temperature, although little is known about how this might be linked to the coordination hypothesis. Acclimation results in R dark being of similar magnitude in plants grown at different temperatures when measured at their respective growth temperatures (Larigauderie and Körner, 1995;Atkin and Tjoelker, 2003) and also results in R dark at 25 • C, increasing upon cold acclimation and declining upon acclimation to warmer temperature . Growth temperature-dependent changes in R dark at a standard temperature can occur over periods of 1-3 days (Atkin et al., 2000;Bolstad et al., 2003;Lee et al., 2005;Zaragoza-Castells et al., 2007;Armstrong et al., 2008). A data synthesis of global patterns in R dark  showed that geographic variation in R dark at growth temperature from the tropics to the tundra is much smaller than would be expected on the basis of enzyme kinetics. In that study, leaf R dark at standard temperature was shown to increase with decreasing growth temperature at a global scale . Observed global patterns in R dark are consistent with the acclimation of R dark to global patterns in growth temperature (Slot and Kitajima, 2015;Vanderwel et al., 2015).
Many ecosystem models assume optimality criteria for stomatal behaviour, in which carbon assimilation is traded off against water loss. Prentice et al. (2014) provided field evidence supporting the least-cost hypothesis, stating that plants adopt an optimal c i : c a ratio that minimizes the combined costs per unit carbon assimilation of maintaining the capacities for carboxylation (V cmax ) and water transport. This hypothesis predicts that the c i : c a ratio should increase with temperature due to lower water viscosity (reducing water costs) and higher photorespiration (increasing carboxylation costs)  while declining with vapour pressure deficit (VPD) due to increasing water costs (Prentice et al., 2011).
In the current study, we present leaf-level measurements carried out during the warm and the cool seasons in the semi-arid environment of the Great Western Woodlands of south-western Australia. By sampling during both seasons, we were able to compare observations under broadly overlapping light conditions but at contrasting temperatures. We are not aware of any previous study that has tested whether seasonal temperature acclimation is consistent with the coordination hypothesis. We explore the idea by comparing the field-observed relationships of each trait to temperature with the theoretical acclimation of photosynthetic traits (as predicted by the coordination hypothesis) and with the alternative, i.e. the relationship of each trait to temperature to be expected if it were controlled only by enzyme kinetics (no acclimation).

Site
Eight representative woody species were studied at the Great Western Woodlands SuperSite (17 • 07 S, 145 • 37 E) approximately 70 km north-west of Kalgoorlie, Western Australia. The area has a semi-arid climate and the vegetation is a wellpreserved mosaic of temperate woodland, shrubland and mallee. The mean annual precipitation (MAP) during 1970-2013 was 380 mm (Hutchinson, 2014c). The average precipitation is slightly higher during summer months, but this rain often falls during short periods as intense storms. Mean annual temperature (MAT) is 20 • C (Hutchinson, 2014a, b). Mean monthly daily temperature minima range between 6 and 18 • C and maxima between 17 and 35 • C. The area is not prone to large shifts in temperature or VPD within days (e.g. cold fronts). Data for daily temperature and shortwave radiation were obtained from the flux tower (TERN Ozflux, http://www.ozflux.org.au, last access: 21 May 2018). All trees were sampled within a 5 km radius from the tower. The species studied were the angiosperm trees Eucalyptus clelandii, E. salmonophloia, E. salubris and E. transcontinentalis and the shrub Eremophila scoparia; the nitrogenfixing leguminous tree Acacia aneura and shrub A. hemiteles; and one gymnosperm tree, Callitris columellaris. All are evergreen.

Gas exchange measurements
We measured 109 A-c i curves altogether during the warm season (late March-early April) and the cool season (late August-early September). The same individual plants were sampled in both seasons. A portable infrared gas analyser system (LI-6400; Li-Cor, Inc., Lincoln, NB, USA) was used. Sunlit terminal branches from the top one-third of the canopy were collected and immediately re-cut under water (Domingues et al., 2010). One of the youngest fully expanded leaves, attached to the cut branch, was placed in the leaf chamber. Measurements in the field were taken with the chamber block temperature close to the air ambient temperature. The CO 2 partial pressure in the chamber for the A-c i curves proceeded stepwise down from 400 to 35, back to 400 and then up to 2000 µmol mol −1 . Prior to the measurements, we tested plants to determine appropriate light-saturation levels. The photosynthetic photon flux density (PPFD) adopted for measurement ranged between 1500 and 1800 µmol m −2 s −1 . After measuring the A-c i curves over about 35 min, light was set to zero for 5 min before measuring R dark . This was a time-saving compromise to allow four or five replicate curves per machine per day, based on our experience that stable results are commonly obtained after 5 min. Moreover, this quick estimate should still be superior to the common practice of deriving R dark as one of the parameters in a curve-fitting routine. Following Domingues et al. (2010), we discarded 23 A-c i curves in which stomatal conductance declined to very low levels, resulting in 86 curves being used in further analyses. TPU (triose phosphate utilization) limitation (Sharkey et al., 2007) was not considered, as it would be unlikely to occur at our field temperatures of above 17 • C.
In theory optimal growth temperature (T opt ) could affect the calculation of V cmax and J max . However, Medlyn et al. (2002b) and Kattge and Knorr (2007) have found very good correlations between T opt and mean daily temperature; and we measured A-c i curves between morning and early afternoon, avoiding the hottest part of the day. Therefore, it is unlikely that any measurements were carried out above the optimum temperature.
The primary data are available via the TERN Su-perSites Data Portal (Prober et al., 2015). Reported ratios of c i : c a relate to chamber conditions, with ambient CO 2 ≈ 400 µmol mol −1 .

Photosynthetic parameters and their temperature responses
Apparent values of V cmax and J max were fitted using the Farquhar et al. (1980) model. Values were standardized to 25 • C (V cmax 25 and J max25 ) using the in vivo temperature dependencies given in Bernacchi et al. (2001Bernacchi et al. ( , 2003. Following Bernacchi et al. (2009), we used the Arrhenius equation to describe the kinetic temperature responses of V cmax and J max : where H a is the activation energy (J mol −1 ), R is the universal gas constant (8.314 J mol −1 K −1 ), T ref is a reference temperature and T is the leaf temperature (K). To derive R dark25 we applied a temperature-dependent Q 10 (fractional change in respiration with a 10 • C increase in temperature) equation in which Q 10 declines with increasing leaf temperature (Atkin and Tjoelker, 2003): where 3.09 and 0.04 are empirical constants, T Cref = 25 • C and T C is leaf temperature ( • C).

Nutrient analyses
After completion of each A-c i curve, leaves were retained to determine leaf area, dry mass, and mass-based nitrogen (N) and phosphorus (P) concentrations (mg g −1 ). Leaves were sealed in plastic bags containing moist tissue paper to prevent wilting. Leaf area was determined using a 600 dots/inch flatbed top-illuminated optical scanner and Image J software (http://imagej.nih.gov/ij/, last access: 1 April 2016). Leaves were dried in a portable desiccator for 48 h, to be preserved until the end of the campaign. Subsequently, in the laboratory, leaves were oven-dried for 24 h at 70 • C and the dry weight was determined (Mettler-Toledo Ltd, Port Melbourne, Victoria, Australia). Leaf mass per unit area (LMA; g m −2 ) was calculated from leaf area and dry mass. N mass and P mass were obtained by Kjeldahl acid digestion of the same leaves (Allen et al., 1974). The leaf material was digested using 98 % sulfuric acid and 30 % hydrogen peroxide. Digested material was analysed for N and P using a flow injection analyser system (LaChat QuikChem 8500 Series 2, Lachat Instruments, Milwaukee, WI, USA). Area-based N and P values (N area , P area , mg m −2 ) were calculated as the products of LMA and N mass or P mass .

Statistical analyses
All statistics were performed in R (R Core Team, 2012). For graphing we used the ggplot2 package (Wickham, 2010). V cmax , V cmax 25 , J max , J max25 , R dark , R dark25 , N area and P area data were naturally log-transformed to approximate a normal distribution of values. Log-transformation of the photosynthetic parameters has the added advantage that the slopes of log-transformed values with respect to temperature are expected to be approximately linear due to the exponential nature of biochemical temperature dependencies generally. The ratio c i : c a , however, was not transformed because of its small variance and approximately normal distribution in this study. Linear regression was used to test dependencies among the parameters. Slopes and elevations of regressions were compared using standardized major axis regression with the smatr package (Warton et al., 2006). The Welch two-sample t test was used to test pairwise differences in traits (e.g. differences between the warm season and the cool season measurements). Generalized linear models (GLMs) were used to test acclimation to temperature across species, with temperature as the single quantitative predictor and species as a factor.

Comparison of observed and predicted responses to temperature
Regression slopes of each naturally log-transformed variable to temperature (obtained using species as a factor, to control for differences in the magnitudes of values for different species) were compared with theoretically derived values The three solid lines indicate the instantaneous response of V cmax to temperature for plants grown at 20, 25 and 30 • C under identical light conditions. The dashed line indicates the predicted response of field-measured (acclimated) V cmax to the growth temperature. The instantaneous responses are due to enzyme kinetics, whereas the acclimated response is due to the slightly higher V cmax required to compensate for increased photorespiration at higher temperatures. Because the acclimated response is shallower than the instantaneous response, V cmax at a common temperature (e.g. 25 • C) declines with growth temperature.
based on alternative hypotheses: (a) based on enzyme kinetics (without acclimation), and (b) based on the coordination hypothesis for V cmax (Fig. 1), J max and R dark and the leastcost hypothesis for c i : c a . "Kinetic" temperature responses were predicted from Eq. (1), which can also be written as follows: which in turn is well approximated by linearizing around T ref : where T = T − T ref . Thus, the kinetic slopes of ln param (T ) vs. T can be predicted from Eq. (4) and compared directly with the fitted slopes obtained with a GLM. We set T ref = 298 K (which is both conventional and close to the median measurement temperature in our data set). We applied H a values based on in vivo measurement at 25 • C by Bernacchi et al. (2001) for V cmax (65 330 J mol −1 K −1 ).
For J max , we used the H a value based on in vivo measurements by Bernacchi et al. (2003) on plants that had been grown at 25 • C (43 900 J mol −1 K −1 ). For R dark , we used a value of H a (50 230 J mol −1 K −1 ) calculated from Eq. (2).
Although H a values can vary among species and with temperature (Von Caemmerer, 2000), for simplicity we adopted the same H a for all species and temperatures. This approximation could affect interspecies comparisons of H adependent parameters but should not interfere with the comparison of theoretical and fitted slopes. The "acclimated" response of V cmax according to the coordination hypothesis was predicted by setting the Rubiscoand electron-transport-limited rates of photosynthesis to be equal. To simplify matters, we disregarded the curvature of the response of electron transport to PPFD at values below the point of Rubisco limitation, giving where ϕ 0 is the intrinsic quantum efficiency of photosynthesis, I abs is the absorbed PPFD, K is the effective Michaelis-Menten coefficient for carbon fixation and * is the photorespiratory compensation point. The theoretical (acclimated) sensitivity of V cmax to temperature was then calculated from the derivative of Eq. (5): We evaluated Eq.
(6) at T = 298 K and c i = 200 µmol mol −1 (approximately the median of our observed values of c i ), using the temperature dependencies of K and * from Bernacchi et al. (2001). The temperature dependency of K was determined from those of the constituent terms K C and K O (the Michaelis-Menten coefficients for carboxylation and oxygenation, respectively) (Wang et al., 2017). We estimated the theoretical acclimated slope for J max as being the acclimated slope of V cmax , minus the difference between the kinetic slopes of V cmax and J max ; this is equivalent to assuming that the ratio of J max to V cmax at 25 • C is constant. We further assumed that R dark (on acclimation) should be an invariant fraction of V cmax , implying the same acclimated temperature response for R dark as for V cmax .
The least-cost hypothesis  provides an optimal value for c i : c a , denoted by χ o , such that where b is the (assumed constant) ratio of R dark to V cmax and a is the cost of maintaining the water-transport pathway. Both K and a are temperature-dependent (a because it is proportional to the viscosity of water, η). Holding vapour pressure deficit (D) constant, we obtained an expression for ∂χ * /∂T , and from Eqs. (7) and (8), We evaluated Eq. (9) at T = 298 K, χ o = 0.5 and c a = 400 ppm using the known temperature dependencies of K and η.

Relationships among photosynthetic parameters
When measured at near-ambient air temperature, speciesaverage V cmax values ranged across seasons from 44.4 to 105 µmol m −2 s −1 , J max from 77.4 to 160 µmol m −2 s −1 , R dark from 1.16 to 3.14 µmol m −2 s −1 and c i : c a from 0.39 to 0.60 (at ambient CO 2 ≈ 400 µmol mol −1 ). At the prevailing air temperatures, V cmax , J max and R dark were systematically higher in the warm season, while their values standardized to 25 • C were higher in the cool season (Fig. 2). The c i : c a ratio also exhibited significantly higher average values in the warm season for six out of eight species (not shown). V cmax and J max were closely correlated across species within and across seasons (Fig. 3). The warm and the cool season regression equations relating V cmax and J max were significant and statistically indistinguishable. The warm and the cool season slopes of regressions forced through the origin are shown in Fig. 3. (We show regressions through the origin here because the intercept (a) is extrapolated well beyond the range of the data, and (b) has no clear biological meaning.) Regressions between V cmax and J max for the warm and the cool seasons together yielded J max = 0.84 V cmax + 55.2 (p < 0.05, slope standard error = ±0.2). There were positive correlations between R dark25 and both V cmax and J max across seasons, as well as strong negative correlations between species-average c i : c a ratios and both V cmax and J max across seasons (Fig. 4). Log-transformed R dark was only correlated to V cmax and J max (p < 0.05) for the individual species E. salmonophloia and C. columellaris.

Leaf gas exchange trait responses to temperature
Based on data from all species together, V cmax , J max and R dark increased with leaf temperature, while the corresponding normalized (to 25 • C) values declined with leaf temperature (p < 0.05, Fig. 5). The c i : c a ratio also increased slightly but significantly with leaf temperature, (p < 0.05; Fig. 6a). The ratio J max : V cmax did not correlate with temperature based on the data from all species together, but it was negatively correlated with temperature for E. salmonophloia, E. scoparia and C. columellaris (not shown). Excluding the two N-fixing species, and/or the one gymnosperm, from the data set had no effect on these results.
The relationship between ambient air temperature (T air , • C) and leaf temperature (T leaf , • C) as measured in the gas-exchange system was T leaf = 1.01 T air + 0.35 (p < 0.05,  R 2 = 0.96). Regression slopes between photosynthetic parameters and T air showed no significant differences from those calculated using T leaf , but the goodness of fit was weaker with T air than with T leaf . We also fitted regressions using T day (the daily mean temperature). Again, the slopes did not change, but the goodness of fit was further reduced. The factor of season (included as a predictor in a generalized linear model, in addition to T leaf ) did not improve model fit.
Within individual species, we also found positive responses of V cmax and J max to temperature, and negative responses when the parameters were standardized to 25 • C (Fig. 7). The response of the c i : c a ratio to leaf temperature was similar in most species (Fig. 6b). Within-species responses of R dark to leaf temperature were weaker and less consistent (Fig. 7), suggesting that respiration had acclimated to a greater extent than was the case for V cmax and J max .
The average VPD value during the warm season was 1.13 kPa, and during the cool season it was 0.45 kPa. The average leaf-to-air VPD (i.e. corrected to leaf temperature) during the warm season was 2.5 kPa, and during the cool season it was 1.44 kPa. However, there was very little difference in stomatal conductance at light saturation: 0.064, 0.082 and 0.101 (lower quartile, median and upper quartile) mol m −2 s −1 for the warm season, and 0.057, 0.078 and 0.085 mol m −2 s −1 for the cool season.
Incoming short-wave radiation at the surface was used as a proxy for PPFD. Daily values ranged from 90 to 256 W m −2 and averaged 193 W m −2 . Averages for the warm and the cool-season sampling periods were not significantly different. None of the photosynthetic parameters showed any relationship with short-wave radiation. and species (n = 8) for photosynthetic capacity, V cmax and J max with R dark (µmol m −2 s −1 ) and the c i : c a ratio at ambient CO 2 ≈ 400 µmol mol −1 (p < 0.05). V cmax , J max and R dark were log ln transformed; c i : c a was not. Each point represents one A-c i curve (n = 86).

Leaf N and P relationships to photosynthetic traits and temperature
Area-based rates of leaf gas exchange traits were not systematically related to total leaf N area or P area (not shown). There was a positive relationship between N and P (by mass) when taking all species together and within three of the species (p < 0.05; Fig. 8). High values of foliar N : P ratios (> 16) in seven out of eight species (Fig. 9) may suggest that P in this ecosystem is more limiting to growth than N . N : P ratios declined with increasing temperature (p < 0.05). A. aneura and A. hemiteles presented the highest N : P, as expected for N-fixing species. We investigated how the ratio of V cmax (and V cmax25 ) to N area varied with ambient temperature. The results were generally very similar to the analysis of V cmax (and V cmax25 ) to temperature (Figs. 5 and 7), but they were significant (p < 0.05) for seven out of eight species (all except the gymnosperm Callitris columelaris) for both V cmax and V cmax25 (not shown).

Quantitative temperature responses
For V cmax and J max , the fitted slopes with leaf temperature were shallower than the predicted kinetic slopes by a margin that greatly exceeded their 95 % confidence limits (Table 1). The kinetic values are what would be expected if the activities of the relevant enzyme complexes remained constant with changing growth temperature. The coordination hypothesis predicts much shallower acclimated rate temperature slopes (Fig. 1). The acclimated slope for V cmax can be predicted from the temperature dependencies of K and the photorespiratory compensation point ( * ). The acclimated slope for V cmax falls just marginally above the 95 % confidence interval for the fitted slope. The acclimated slope of J max falls well within the 95 % confidence interval for the fitted slope. For R dark the acclimated and kinetic slopes are closer together, and both fall within the 95 % confidence interval of the fitted slope.
There is no kinetic response of c i : c a , but the least-cost hypothesis predicts a positive response to temperature. This was observed, although the fitted slope of the response to temperature was shallower than predicted. Figure 5. Bivariate linear regressions (p < 0.05) of naturally log-transformed V cmax , V cmax25 , J max , J max25 , R dark and R dark25 (µmol m −2 s −1 ) vs. leaf temperature (T leaf , • C). Each point represents one A-c i curve (n = 86).

Quantitative ranges of photosynthesis traits
Values of V cmax , J max , R dark and c i : c a for the eight species measured here were within commonly reported ranges. V cmax and J max were generally lower than expected for desert species (153 ± 54 µmol m −2 s −1 ) but higher than typical values for mesic perennial species (53 ± 15 µmol m −2 s −1 ) (Wullschleger, 1993). The values were also high compared with trees from savannas with twice the annual precipitation (range 20-90 µmol m −2 s −1 ) (Domingues et al., 2010). The c i : c a ratios fall within the range typical for dry environments .

Comparison between seasons
Our results were consistent with acclimation of photosynthetic traits to temperature as predicted by the coordination hypothesis. When measured at the prevailing ambient temperature, V cmax , J max and R dark were all generally higher in the warm season than in the cool season, whereas values standardized to 25 • C were generally lower in the warm season than in the cool season (Fig. 2). This is prima facie evidence for active seasonal acclimation, in agreement with the coordination hypothesis that predicts lower allocation of N to Rubisco and other photosynthetic enzymes at higher temperatures, offsetting the increase in enzyme activity with elevated temperature. Moreover, the absence of acclimation (as found in Way and Yamori, 2014) should result in no relationship between V cmax25 and leaf temperature, while the negative relationship found here is evidence of acclimation (Figs. 5 and 7).
The comparison of fitted and theoretically predicted slopes (Table 1) reveals not only that the responses of V cmax , J max and R dark to ambient temperature were smaller than would be predicted from enzyme kinetics alone but also that the observed responses were close to, or (in the case of J max ) statistically indistinguishable from, the responses predicted by the coordination hypothesis. The response of c i : c a is in the same direction (positive) as the response predicted by the least-cost hypothesis but only about half as large, probably due to the opposing effect (reduction in c i : c a ) of larger VPD in the warm season than in the cool season .
V cmax and J max were strongly and positively correlated across species (Fig. 3), and the relationship did not shift significantly between seasons. Some studies have reported a lower J max : V cmax ratio in warmer seasons compared to cooler seasons (Medlyn et al., 2002a;Lin et al., 2013). Our data show a tendency in this direction but not enough to be significant (Fig. 3). V cmax and J max have previously been reported to increase seasonally with leaf temperature. In one study on six Eucalyptus species, measurements were taken at six temperature levels in winter, spring and summer; there was an increase in V cmax and J max with air temperature in seasons with overlapping temperatures, and V cmax25 was significantly higher in the winter than in the summer (Lin et al., 2013). Miyazawa and Kikuzawa (2006) obtained similar results in five evergreen broadleaved species. Our measurements also yielded similar results (Figs. 5 and 7).

Seasonal variability in leaf nutrients
Levels of leaf N and P have been reported to change seasonally (Medlyn et al., 2002a, and Fig. 9). We found a reduction in the N : P ratio in the warm season, consistent with a reduced allocation of N to photosynthetic functions (Way and Sage, 2008). A reduction in total leaf N does not necessarily indicate changes in N allocation; however, a strong coupling between N and photosynthesis has been widely observed, even though Rubisco accounts for only 10-30 % of the total leaf N (Evans, 1989). This coupling can be verified in our data, represented by the relationship between V cmax / N area and temperature across the two seasons, which was significant for a greater number of species in comparison to the relationship between V cmax and temperature.
A reduction in leaf N in the warm season is unlikely to be caused by general growth dilution during an actively growing part of the year because (a) this is an environment with a year-round growing season, (b) PPFD during the periods of the field campaigns was similar, and (c) the measured V cmax values were shown to be consistent with the coordination hypothesis, implying similar assimilation rates in the two seasons. We did not find significant relationships of photosynthetic traits to foliar N or P despite the study area being extremely limited in supplies of both nutrients (Prober et al., 2012).

Links between photosynthetic activity, R dark and
c i : c a V cmax , J max and R dark were positively correlated with leaf temperature across a wide range (cool season 17 to 27 • C; Figure 7. Bivariate linear regressions of naturally log-transformed V cmax , V cmax25 , J max , J max25 , R dark and R dark25 (µmol m −2 s −1 ) vs. leaf temperature (T leaf , • C) within species (p < 0.05). Only significant regressions (p < 0.05) are shown. Each point represents one A-c i curve (n = 86).
warm season 26 to 37 • C), both for the data set as a whole and within individual species (Figs. 5 and 7). Photosynthetic capacity and respiratory flux are linked via the ATP (adenosine triphosphate) demand of sucrose synthesis and transport, leading to the interdependence of chloroplast and mitochondrial metabolism (Krömer, 1995;Ghashghaie et al., 2003). The parallel temperature acclimation of R dark and V cmax illustrates this close relationship. Across all species, there was a strong negative relationship between c i : c a and photosynthetic capacity. A simple interpretation of this (common) finding is that a higher photosynthetic capacity allows a stronger drawdown of CO 2 between c a and c i . However, this interpretation assumes that a higher V cmax means a higher assimilation rate, which cannot be true if V cmax is already optimized to the current light environment. Instead, the coordination hypothesis implies that the lower the c i : c a ratio, the greater the photosynthetic capacity that is required in order to achieve a given assimilation rate. Measured c i : c a ratios in this study increased with temperature as predicted, but the slope of 0.006 (Fig. 6) based on all data was shallower than the predicted slope of 0.013. This difference may reflect higher VPD in the warm season, which would be expected to close stomata and therefore act in the opposite way to the effect of temperature alone .

Implications for modelling
Data that test predictions of the coordination hypothesis regarding the seasonal acclimation of V cmax , J max and R dark are scarce. This study has presented evidence suggesting that such acclimation (at least in evergreen woody plants) is a real phenomenon, and is quantitatively consistent with the coordination hypothesis. Terrestrial models that do not allow seasonal acclimation therefore may incorrectly represent the seasonal time course of carbon exchange at the plant and ecosystem levels.
One dynamic global vegetation model, the Lund-Potsdam-Jena (LPJ) model (Sitch et al., 2003), together with several later models based on LPJ, formally assumes the coordination hypothesis (as well as the coupling between R dark and V cmax ) and thus implicitly allows photosynthetic parameters and leaf respiration to acclimate to the seasonal variation of climate. The coordination hypothesis is also invoked by the new "first-principles" global primary production model developed by Wang et al. (2014) and further elaborated and tested by Wang et al. (2017). However, there Table 1. Linear regression slopes (K −1 ) and their 95 % confidence intervals for naturally log-transformed photosynthetic traits, with species included as a factor. The values are compared to kinetic slopes (as expected in the absence of acclimation) and acclimated slopes, as predicted by the coordination hypothesis for V cmax , J max and R dark and the least-cost hypothesis for c i : c a ; n/a means not applicable. are several steps between these models and the potentially wider application of the coordination hypothesis in physical land surface and Earth System modelling. Timescales are a key issue. The timescale of acclimation of V cmax and J max is unclear, and could not be resolved empirically without frequently repeated measurements. It is generally accepted, however, that acclimation occurs over multiple days. The LPJ model requires daily or monthly inputs and assumes complete acclimation on monthly timescales. Wang et al. (2014Wang et al. ( , 2017 implemented their model on monthly time steps. Modelling diurnal cycles would require a separation of timescales, such that photosynthetic capacities would be near-constant over a diurnal cycle but would vary gradually in response to the seasonal cycle. Thus, whereas Wang et al. (2014Wang et al. ( , 2017 could make V cmax and J max "disappear" from equations describing monthly primary production (because on this timescale they were considered to depend only on the environment), multitemporal applications will have to explicitly predict time-varying values of V cmax and J max depending on antecedent environmental conditions. Field measurements of photosynthetic parameters in different seasons are required to test these predictions .
Author contributions. HFT, ICP, OKA, CM and SMP planned and designed the research. HFT performed all research and data analysis. HFT measured all the field data used in the experiment. ICP developed the theoretical predictions. HFT wrote the first draft; all authors commented on subsequent versions and assisted with data interpretation.
Competing interests. The authors declare that they have no conflict of interest.