Tree Physiology 26, 1173–1184 © 2006 Heron Publishing—Victoria, Canada Seasonal fluctuations and temperature dependence of leaf gas exchange parameters of co-occurring evergreen and deciduous trees in a temperate broad-leaved forest YOSHIKO KOSUGI 1,2 and NAOKO MATSUO1 1 Graduate School of Agriculture, Kyoto University, Kyoto 606-8502, Japan 2 Corresponding author ([email protected]) Received March 21, 2005; accepted September 22, 2005; published online June 1, 2006 Summary Seasonal fluctuations in leaf gas exchange parameters were investigated in three evergreen (Quercus glauca Thunb., Cinnamomum camphora Sieb. and Castanopsis cuspidata Schottky) and one deciduous (Quercus serrata Thunb.) co-occurring, dominant tree species in a temperate broadleaved forest. Dark respiration rate (R n ), maximum carboxylation rate (Vcmax ) and stomatal coefficient (m), the ratio of stomatal conductance to net assimilation rate after adjustment to the vapor pressure deficit and internal carbon dioxide (CO2) concentration, were derived inversely from instantaneous field gas exchange data (one-point method). The normalized values of R n and Vcmax at the reference temperature of 25 °C (R n25 , Vcmax25 ) and their temperature dependencies (ΔH a (R n ), ΔH a (Vcmax )) were analyzed. Parameter Vcmax25 ranged from 24.0–40.3 µmol m – 2 s – 1 and ΔH a (Vcmax ) ranged from 29.1– 67.0 kJ mol – 1. Parameter R n25 ranged from 0.6–1.4 µmol m – 2 s – 1 and ΔH a (R n ) ranged from 47.4–95.4 kJ mol – 1. The stomatal coefficient ranged from 7.2–8.2. For the three evergreen trees, a single set of Vcmax25 and R n25 parameters and temperature dependence curves produced satisfactory estimates of carbon uptake throughout the year, except during the period of simultaneous leaf fall and leaf expansion, which occurs in April and May. In the deciduous oak, declines in Vcmax25 were observed after summer, along with changes in Vcmax25 and R n25 during the leaf expansion period. In all species, variation in m during periods of leaf expansion and drought should be considered in modeling studies. We conclude that the changes in normalized gas exchange parameters during periods of leaf expansion and drought need to be considered when modeling carbon uptake of evergreen broad-leaved species. Keywords: dark respiration rate, maximum carboxylation rate, photosynthesis, stomatal conductance, temperate broadleaved forest. Introduction Leaf gas exchange characteristics must have an important bearing on forest carbon uptake (e.g., Baldocchi and Meyers 1998, Wilson et al. 2001, Lai et al. 2002), yet, to date, studies of leaf-scale gas exchange have yielded few results applicable to modeling long-term canopy or landscape carbon exchange (Xu and Baldocchi 2003, Kosugi et al. 2003). It is still unclear how gas exchange parameters differ among species, seasons, locations and environmental conditions in the field, although it has recently been recognized that it is essential to include seasonal changes in gas exchange characteristics to model accurately the carbon uptake of deciduous broad-leaved forests. This requirement was confirmed by Wilson et al. (2001) who showed that modeled net ecosystem exchange of carbon better matched eddy covariance measurements when seasonal variation in maximum carboxylation rate (Vcmax ; Wilson et al. 2000) was taken into account. Kosugi et al. (2003) used field gas exchange data to demonstrate that accurate modeling of gas exchange in deciduous broad-leaved trees across seasons requires information on seasonal changes in Vcmax, dark respiration rate (R n ) at a reference temperature and stomatal coefficient (m), the ratio of stomatal conductance to net assimilation rate after adjustment for the vapor pressure deficit and internal carbon dioxide (CO2) concentration in the Ball-type stomatal conductance model, during leaf expansion. Spatial variations in gas exchange parameters have also been reported in a deciduous and an evergreen broad-leaved tree (Dungan et al. 2003) and a pine tree (Han et al. 2003). Seasonal and spatial variations in gas exchange parameters can be explained, in part, by changes in leaf nitrogen content, although the relationship with leaf nitrogen is not necessarily constant. Thus, gas exchange parameters and their temporal and spatial variations cannot be described simply as a constant function of leaf nitrogen content. The influences of nitrogen allocation, acclimation, enzyme activity and other factors on leaf gas exchange parameters need, therefore, to be quantified if the relationship between gas exchange at the leaf scale and at the canopy scale, is to be understood. There is a need for information, based on long-term field observations, about which parameter set should be used for long-term and larger scale simulations of each forest type. Although several studies have clarified seasonal changes in gas exchange parameters in temperate deciduous broad-leaved forest, there is little information about seasonal fluctuations in 1174 KOSUGI AND MATSUO gas exchange characteristics in other types of forests such as temperate evergreen broad-leaved forests (Tenhunen et al. 1990). Consequently, most larger scale simulation models, which are now one of the main tools for evaluating carbon uptake by forests, use limited sets of published leaf gas exchange parameters. Detailed information on the seasonal variation in the value of these parameters in evergreen forests would allow us to determine if we can evaluate carbon uptake for all seasons with one set of leaf gas exchange parameters. We chose to study three evergreen and one deciduous tree species co-occurring and dominant in a temperate evergreen broad-leaved forest in western Japan. We investigated seasonal fluctuations in three major gas exchange parameters (Rn , Vcmax and stomatal coefficient of the improved Ball-type stomatal conductance modelt (m) (Leuning 1995)) to obtain primary information for a canopy-scale long-term estimate of gas exchange with a process-based multi-layer model based on field observations (e.g., Tanaka et al. 2002). Few studies have examined the temperature dependencies of major gas exchange parameters of canopy leaves in the forest based on in situ field observation, especially in temperate evergreen broad-leaved forests. The objectives of our study were: (1) to develop a practical method for determining the magnitude and seasonal fluctuations in leaf gas exchange parameters of the study species as a basis for modeling forest gas exchange; and (2) to determine if it is possible to evaluate the carbon uptake of a temperate evergreen broad-leaved forest in all seasons with a single set of leaf gas exchange parameters. Li-Cor, Lincoln, NE) under light-saturated conditions in three sunlit current-year or previous-year leaves in the canopy of each tree studied. Measurements were made in the morning (between 0830 and 1300 h) on April 27, May 10, 17, 24, June 7, 26, July 12, August 30 and October 16 in 2001 and on January 23 and March 6 in 2002. Measurements were made once on each leaf with an irradiance of 1000 µmol m – 2 s – 1 supplied by a red/blue LED (Li-Cor, LI-6400-02) under ambient temperature, humidity and CO2 concentration. Previous studies of gas exchange in leaves of these trees indicated that maximum gas exchange rates were usually reached at a photosynthetic photon flux (PPF) of 1000 µmol m – 2 s – 1, at which irradiance there was neither a sharp increase in leaf temperature nor a depression in photosynthesis. After the gas exchange measurements in the light, dark gas exchange rates of the same leaves were measured with the LI-6400 by placing the leaves in a dark chamber under ambient conditions. Environmental conditions were measured at a meteorological observation tower close to the gas exchange observation tower. Air temperature and relative humidity above the canopy were measured with a Visala-type hygrothermometer (HMP35C, Campbell Scientific), downward solar radiation was measured with a four-component radiometer (MR-40, Eiko, Japan), rainfall data were provided from a fire station located approximately 1 km from the plantation, soil capillary pressure was measured with tensiometers buried at depths of 10, 20, 30 and 40 cm and volumetric soil water content was measured with a water content reflectometer (CS615, Campbell Scientific) buried in the soil beside the tower at depths of 0–30 cm. Materials and methods Parameterization procedure Observations Field observations were made in a forest plantation established in 1987 at 34°44′ N, 134°22′ E in Akou, Hyogo prefecture on the east coast of Japan. The plantation includes 12 temperate evergreen trees species and four deciduous broad-leaved tree species (Kosugi et al. 2005). Three-year-old seedlings were planted in 1987 at a density of 1–2 stems m – 2. In January 2002, trees had a mean height of 5.3 m (including small shaded trees) and a mean basal area of 34 m – 2 ha – 1. Tree density was 11,150 stems ha – 1. We made leaf gas exchange measurements from April 2001 to March 2002 on one tree of each of three evergreen broad-leaved species (Quercus glauca Thunb., Cinnamomum camphora Sieb. and Castanopsis cuspidata Schottky) and one co-occurring deciduous oak species (Quercus serrata Thunb.) growing near a gas exchange observation tower in the center of the forest. In the evergreen trees, simultaneous leaf fall and leaf expansion occur in April and May, when the majority of previous-year leaves abscise and most current-year leaves expand. Expansion occurs from mid-April in Q. glauca and C. camphora and from early May in C. cuspidata. C. camphora also expands some leaves in July and August. In the deciduous oak (Q. serrata), expansion of new leaves occurs from mid-April and leaf senescence occurs in November. Gas exchange measurements were made with a portable steady-state photosynthesis measurement system (LI-6400, We derived three major gas exchange parameters: normalized dark respiration rate (R n25 ), normalized maximum carboxylation rate (Vcmax25 ) and m. The model used for the analysis was a modified Farquhar et al. (1980) biochemical model of photosynthesis for C3 plants. Net CO2 assimilation rate is described by: ⎛ p ( Γ* ) ⎞ A = Vc ⎜1 – ⎟ – Rd p ( C c )⎠ ⎝ (1) p (O) 2τ (2) p ( Γ* ) = where A is net assimilation rate (µmol m – 2 s – 1 ), Vc is rate of carboxylation in the PCR cycle (µmol m – 2 s – 1 ), R d is day (non-photorespiratory) respiration rate (µmol m – 2 s – 1 ), p(Γ* ) is the CO2 compensation point without non-photorespiraory respiration (Pa), τ the specificity factor of Rubisco and p(C c ) (Pa) and p(O) (21,000 Pa) are the partial pressures of CO2 and oxygen (O2) at the sites of carboxylation and oxygenation, respectively. The lower value between the ribulosebisphosphate (RuBP) regeneration or electron-transport limited rate of carboxylation (Wj ; µmol m – 2 s – 1) and the RuBP saturated rate of carboxylation (Wc ; µmol m – 2 s – 1 ) is taken as the velocity of carboxylation (Vc ), with Wc expressed as: TREE PHYSIOLOGY VOLUME 26, 2006 SEASONAL FLUCTUATIONS IN LEAF GAS EXCHANGE Wc = Vcmax p (C c ) ⎛ p (O) ⎞ p ( C c ) + K c ⎜1 + ⎟ Ko ⎠ ⎝ (3) where K c and K o (Pa) are the Rubisco Michaelis-Menten constants for CO2 and O2 , respectively. Parameter Vcmax (µmol m – 2 s – 1 ) was estimated by the one-point method as described by Kosugi et al. (2003). Briefly, Vcmax was derived from any point within the low intercellular carbon dioxide concentration (Ci ) range of the A–C i relationship after determining R d. Thus, after excluding the points at which Wc is not the limiting factor (under field conditions, this usually occurs at low PAR and ambient CO2 concentration), each CO2 and H2O gas exchange data point includes information on Vcmax if we know R d. This method can reveal more information than the analysis of an A–C i curve obtained with a controlled chamber in the field. Values of K c and K o , τ and their activation energies used to calculate temperature dependence were taken from Harley et al. (1992). The Arrhenius function was used to determine the temperature dependencies of parameters K c , K o , τ and R n as: ⎛ ⎡ T ⎤ ΔHa ⎞ f (T1. k ) = f (T ref ) exp ⎜ ⎢1 – ref ⎥ ⎟ ⎝ ⎣ T l. k ⎦ RT ref ⎠ (4) A simplified equation from Sharpe and DeMichele (1977) was used to determine the temperature dependence of Vcmax: ⎛ ⎡ T ⎤ ΔHa ⎞ f (T ref ) exp ⎜ ⎢1 – ref ⎥ ⎟ ⎝ ⎣ T1. k ⎦ RT ref ⎠ f ( T1. k ) = ⎛ Δ ST1. k – Δ H d ⎞ 1 + exp ⎜ ⎟ RT1. k ⎝ ⎠ (5) where f (Tl.k ) is the value of a given parameter at leaf temperature Tl.k (K), f (Tref ) is the reference value of that parameter at 25 °C (K c25 , K o25 , R n25 , τ 25 and Vcmax25 ), ΔH a is the activation energy (J mol –1 ), ΔH d is the deactivation energy (J mol –1 ) and ΔS is an entropy term (J K –1 mol –1 ). First, we determined R n at the reference temperature and its temperature dependence (Equation 4). Here, net assimilation rate measured in the dark chamber was used to derive R n25 and ΔH a (R n ) for each species. During the leaf expansion and leaf fall periods, R n25 was optimized with a nonlinear least squares optimization procedure separately for each species based on the same activation energy as used in the other periods. Day respiration rate (R d ) was scaled using the relationship with R n (Lloyd et al. 1995) based on the results of Brooks and Farquhar (1985): R d = R n , when 0 < Q < 10 µmol m – 2 s – 1 1175 Next, we obtained the instantaneous Vcmax for each value measured under light-saturated conditions by the one point method. Temperature dependence and seasonal fluctuations in Vcmax are discussed based on these data. An improved Ball et al. (1987)-type model (Leuning 1995) was used to estimate stomatal conductance: g sw = m f (D) = A f (D) + g sw min Cs – Γ 1 1+ (7) D D0 where g sw is stomatal conductance of H2O (mol H2O m – 2 s – 1 ), A is net assimilation rate (µmol m – 2 s – 1 ), f (D) is the nonlinear function of the vapor pressure deficit, m is the slope of the relationship between the stomatal index (Af (D)/(C s – Γ)) and stomatal conductance, C s is CO2 concentration at the leaf surface (µmol mol – 1 ), Γ is the CO2 compensation point ( µmol mol – 1 ) and g swmin is minimum stomatal conductance. A Lohammer-type (Lohammer et al. 1980) hyperbolic form function is used for f (D), where D is vapor pressure deficit (kPa) of air and D 0 is the empirical value at which f (D) halves. The CO2 concentration at the leaf surface (Cs ) was calculated based on the CO2 concentration of air in the chamber, a constant boundary layer conductance and net assimilation rate. An important advantage of this improved model of Leuning (1995) for long-term analyses is that it can describe stomatal response to vapor pressure deficit with one curve over a wide temperature range. If we use relative humidity for f (D), as in the original version of Ball-type model (Ball et al. 1987), stomatal response to vapor pressure deficit becomes too large in winter and also differs from the summer curve, making the analysis of the meaning of m obscure. Results Environmental conditions Seasonal fluctuations in daily total solar radiation, daily mean, maximum and minimum air temperatures, daily mean and maximum vapor pressure deficits, daily total precipitation, volumetric soil water content at 0–30 cm and soil capillary pressure at 10, 20, 30 and 40 cm are shown in Figure 1. Summer droughts at the study site varied across years and the influence of severe drought on gas exchange processes has been detected (Tanaka et al. 2002). The year 2001 was moderate in terms of precipitation (annual precipitation was 788 mm in 2000; 1078 mm in 2001; 578 mm in 2002; and 1230 mm in 2003). Dark respiration rate R d = ( 0.5 – 0.05 ln (Q ))R n , when Q > 10 μmol m – 2 s – 1 (6) where Q is the incident PPF (µmol m – 2 s – 1 ). Dark respiration rates are plotted against simultaneously measured leaf temperatures in Figure 2. Here we defined expanding leaves and old leaves as leaves having different R n characteristics compared with mature leaves of the same species. For TREE PHYSIOLOGY ONLINE at http://heronpublishing.com 1176 KOSUGI AND MATSUO Figure 1. Seasonal fluctuations in daily total solar radiation, daily mean, maximum and minimum air temperatures, daily mean and maximum vapor pressure deficits (VPD), daily total precipitation, volumetric soil water content at 0–30 cm and soil capillary pressure at depths of 10, 20, 30 and 40 cm. Figure 2. Dark respiration rate (R n ) observed in three leaves of each tree and on each observation day, plotted against leaf temperature. Circles indicate mean values for three leaves (䊊 = expanding leaves; 䊉 = mature leaves; and 䊉 = old leaves). Error bars are ± 1 standard deviation. all species, R n values were significantly higher for expanding leaves than for mature and old leaves (cf. Xu and Baldocchi 2003, Kosugi et al. 2003). In old leaves (about one month before abscission) of the evergreen trees, R n values were lower than expected based on the temperature dependence curve of mature leaves. After excluding values obtained when leaves were expanding or old, seasonal changes in R n could be explained by one temperature dependence curve for each species. Values of R n25 and ΔH a (R n ) used in the temperature dependence curve for each species (Equation 4) are shown in Table 1. The dark respiration rate at 25°C ranged between 0.6–1.4 µmol m – 2 s – 1 and ΔH a (R n ) ranged between 47.4– 95.4 kJ mol – 1. Maximum carboxylation rate Maximum carboxylation rates are plotted against leaf temperature in Figure 3. We assumed that expanding leaves and old leaves have different Vcmax characteristics compared with mature leaves of the same species. Thus, the data for expanding and old leaves were excluded from the line fitting process. The shift in observed Vcmax from expanding to mature leaves occurred earlier than the shift in R n (May 24 for Q. glauca, June 6 for C. camphora, May 10 for Q. serrata; cf. Figures 2 and 3). The post-summer decline in Vcmax that was detected in the deciduous oak (August 29 and October 16 for Q. serrata) was also considered a characteristic of old leaves. Apart from the TREE PHYSIOLOGY VOLUME 26, 2006 SEASONAL FLUCTUATIONS IN LEAF GAS EXCHANGE 1177 Table 1. Dark respiration rates at the reference temperature (R n25 ), the activation energy (ΔH a (R d )), maximum carboxylation rate at the reference temperature (Vcmax25 ), the activation energy (ΔH a (Vcmax )), the deactivation energy (ΔH d (Vcmax )) and the entropy term (ΔS) in mature leaves of each tree species. Species Rn25 (µmol m – 2 s – 1 ) ΔH a (R n ) (kJ mol – 1 ) Vcmax25 (µmol m – 2 s – 1 ) ΔH a (Vcmax ) (kJ mol – 1 ) ΔH d (Vcmax ) (kJ mol – 1) ΔS (J K – 1 mol – 1 ) Q. glauca Thunb. C. camphora Sieb. C. cuspidata Schottky Q. serrata Thunb. 1.3 1.4 1.0 0.6 47.4 48.0 60.0 95.4 24.0 40.3 33.1(27.8) 3 32.9 59.5 54.9 29.1(56.2) 3 67.0 250.0 1 250.0 1 250.0 1 204.62 600 1 600 1 600 1 650 2 1 2 3 Tenhunen et al. 1990. Kosugi et al. 2003. Values in parentheses in Castanopsis cuspidata are optimized based on summer and autumn data only (June 7–October 16, 2001). leaf expansion and leaf fall periods, seasonal fluctuations in Vcmax for the three evergreen broad-leaved trees could be explained by one temperature dependence curve for each species. A post-summer decline in normalized Vcmax needs to be incorporated into models involving deciduous trees because our results for the deciduous oak corroborate those reported for several other temperate deciduous tree species (Wilson et al. 2001, Xu and Baldocchi 2003, Dungan et al. 2003, Kosugi et al. 2003). Moreover, similar seasonal fluctuations in Vcmax have been observed in wheat (Grossman-Clarke et al. 1999), Sphagnum and Pleurozium (Williams and Flanagan 1998), although Falge et al. (1996) reported no obvious seasonality in normalized Vcmax of an evergreen conifer (Picea abies (L.) Karst.) growing in a forest. Although more field results are necessary to clarify the relationships between seasonality and plant lifestyle, we found that normalized Vcmax in our temperate evergreen broad-leaved trees changed significantly from the year-round value only during the leaf expansion and leaf fall periods. Values of Vcmax25 and ΔH a (Vcmax ) used in the temperature dependence curve (Equation 5) are shown in Table 1. The deacti- vation energy (ΔH d (Vcmax )) and the entropy term (ΔS) were fixed as constant values in the optimization procedure because the data we used provided insufficient information about the behavior of Vcmax at high temperature to evaluate these parameters, except to conclude that there was no obvious decline in Vcmax at temperatures above 30 °C. The constant values used were 250,000 for ΔH d (Vcmax ) and 600 for ΔS, following the results obtained for the evergreen oak Quercus coccifera L. by Tenhunen et al. (1990). For calculations involving the deciduous oak, we used 204, 600 for ΔH d (Vcmax ) and 650 for ΔS (cf. Kosugi et al. (2003). Parameter Vcmax25 ranged among species from 24.0 to 40.3 µmol m – 2 s – 1 and ΔH a (Vcmax ) ranged from 29.1 to 67.0 kJ mol – 1. Stomatal conductance The relationship between stomatal index (defined as Af (D)/(C s – Γ)) and stomatal conductance is shown in Figure 4. The linear relationship between these parameters and the constant m (estimated at between 7.2 and 8.2) appeared to be more or less maintained through all seasons. A more detailed examination of the data showed that there were differ- Figure 3. Maximum carboxylation rate (Vcmax ) observed in three leaves of each tree on each observation day, against leaf temperature. Circles indicate mean values for three leaves (䊊 = expanding leaves; 䊉 = mature leaves; and 䊉 = old leaves). Error bars are ± 1 standard deviation. TREE PHYSIOLOGY ONLINE at http://heronpublishing.com 1178 KOSUGI AND MATSUO ences in the slope m on each observation day. A time series for m (optimized for each tree species and each observation day) against volumetric soil water content is shown in Figure 5. The value of m was large during the leaf expansion period and then gradually settled into a narrow range following leaf maturation, with fluctuations in response to changes in soil water content. Discussion Dark respiration rate The activation energy (ΔH a (R n )) values of three evergreen species obtained in our analyses, based on the assumption of one temperature dependence curve for most of the year, were in the moderate range compared to previously published results (Table 2). Parameter ΔH a (R n ) of the deciduous oak was higher than that reported in other studies listed in Table 2. However, this difference may be unimportant, reflecting only a lack of data at low temperatures. Mixing the effects of acclimation and physiological seasonality with the temperature dependence response curve could pose problems for the analysis of field data. However, Falge et al. (1996) found no obvious differences in ΔH a (R n ) for P. abies across four seasons of the year. Xu and Baldocchi (2003) reported that the normalized R n25 of deciduous Quercus douglasii Hook. & Arn. (blue oak) Figure 4. Relationship between stomatal index (Af(D)/(C s – Γ)) and stomatal conductance observed in three leaves of each tree, on each observation day. Circles indicate mean values for three leaves (䊊 = expanding leaves; 䊉 = mature leaves; and 䊉 = old leaves). Error bars are ± 1 standard deviation. Figure 5. Time series showing changes in stomatal coefficient m (optimized for each tree and each observation day) in relation to volumetric soil water content at 0–30 cm depth. TREE PHYSIOLOGY VOLUME 26, 2006 SEASONAL FLUCTUATIONS IN LEAF GAS EXCHANGE 1179 Table 2. Activation energies of dark respiration (ΔH a (Rd )) and maximum carboxylation rate (ΔH a (Vcmax )), and deactivation energy (ΔH d (Vcmax )) and entropy (ΔS) of maximum carboxylation rate in mature leaves of each tree species investigated in this study and for other species reported in the literature. Species ΔH a (R d ) (J mol – 1 ) Evergreen trees Quercus glauca Thunb. 47400 Cinnamomum camphora Sieb. 48000 Castanopsis cuspidata Schottky 60000 Quercus coccifera L. Arbutus unedo L. Eucalyptus pauciflora Sieber ex. A. Spreng. Macadamia integrifolia Maiden & Betche and Litchi chinensis Senn. Aristotelia serrata Deciduous trees Quercus serrata Thunb. Quercus alba L. White Acer saccharum Marsh. Acer pseudoplatanus L. Betula pendula Roth Fagus sylvatica L. Quercus petraea L. ex Liebl. Quercus robur L. Fraxinus excelsior L. Juglans regia L. Juglans regia L. Fuchsia excorticata J.R. Forst & G. Forst. Conifer Picea abies (L.) Karst. Pinus sylvestris L. Pinus radiata D. ΔH a (Vcmax ) (J mol – 1 ) ΔH d (Vcmax ) (J mol – 1 ) (250000) 1 (250000) 1 (250000) 1 This study This study This study 600 – – Tenhunen et al. 1990 Harley et al. 1986 Kirschbaum and Farquhar 1984 – – Lloyd et al. 1995 – – Dungan et al. 2003 41500 22260 no data – – 53000 (53000) 2 26200–47400 43300 95400 38700 42000 67000 55000 74000 60800 69700 75400 67600 61100 88000 89000 109500 39100 62500–64500 32290–42580 9593–10315 75750 52550–53050 52762–45027 no data 65000 Others Atriplex glabriuscula Edmonston (C4) Spinach? Spinach Cotton Tobacco Mountain grassland 29 species 66405 no data 84450 46390 no data Averaged curve 59 dataset average 7 deciduous trees – – Evergreen trees (n = 7) Deciduous trees (n = 13) Conifers (n = 3) – – – 58520 53000 116300 65330 51300–149061 Source (600) 1 (600) 1 (600) 1 59500 54900 29100 56200 65000 59780 45000 21100 24950 36500 27950 41550 36500 37950 (84450) 4 32900 ΔS (J K – 1 mol – 1) 250000 (204600) 3 (200000) 203000 401000 404000 175000 144000 247000 90000 265000 199500 – (650) 4 656 645 1285 1285 559 451 778 293 851 (650) 4 – (200000) 5 202600–203100 206083–203594 (656) 5 665–669 (650) 4 This study Harley and Baldocchi 1995 Ninemets and Tenhunen 1997 Weber et al. 1985 Dreyer et al. 2001 Dreyer et al. 2001 Dreyer et al. 2001 Dreyer et al. 2001 Dreyer et al. 2001 Dreyer et al. 2001 Dreyer et al. 2001 Le Roux et al. 1999 Dungan et al. 2003 Falge et al. 1996 Wang et al. 1996 Walcroft et al. 1997 – – Badger and Collatz 1977 – – – – 194800–204000 656 Farquhar et al. 1980 Jordan and Ogren 1984 Harley et al. 1992 Bernacchi et al. 2001 Wohfahrt et al. 1999 73637 67600 149252 204600 486 (650) 4 Leuning 2002 Kosugi et al. 2003 55200 65900 56100 > 220000 204000 199500 (650)4 (650) 4 (650) 4 Estimated in this study Estimated in this study Estimated in this study 202900 – 650 – Values in parentheses are assumed to apply, but are not directly determined for the respective species. 1 Tenhunen et al. 1990. 2 Kirschbaum and Farquhar 1984. 3 Kosugi et al. 2003. 4 Harley et al. 1992. 5 Harley and Baldocchi 1995. TREE PHYSIOLOGY ONLINE at http://heronpublishing.com 1180 KOSUGI AND MATSUO remained stable through most of the growing season, except during the leaf expansion period, suggesting that ΔH a (R n ) is constant for several months. Wang et al. (1996) reported that there was no significant difference in ΔH a (R n ) between ambient and elevated temperature treatments. Conversely, Dungan et al. (2003) showed that ΔH a (R n ) in Aristotelia serrata Hook. f. (evergreen wineberry) was slightly lower in winter than in summer, indicating that dark respiration at low temperatures was higher in leaves in winter than in mature leaves in summer. We found no obvious hysteresis in the relationship between temperature and R n in three evergreen broad-leaved trees during summer and autumn (Figure 2). Dark respiration rates increased with increasing temperatures in early spring (March) and dropped below the temperature dependence curve just before leaf fall in April–May (Figure 2). Consequently, we consider that expression and parameterization based on one temperature dependence curve is acceptable for temperate evergreen species, except during the leaf expansion and leaf fall periods. Seasonal fluctuations in observed R n and simulated R n (with the parameter sets listed in Table 1) are compared in Figure 6. At the study site, R n can safely be estimated with one set of parameters, except during the period of simultaneous leaf expansion and leaf fall in April and May in the case of the evergreen trees and during the leaf expansion period in the case of the deciduous oak. by Kosugi et al. (2003) from published data for seven deciduous broad-leaved tree species (67.6 kJ mol – 1 ). These deciduous-tree values also correspond closely to results obtained for tobacco in laboratory experiments (Bernacchi et al. 2001) and are a little larger than the values we obtained for our three evergreen species. For comparison, we derived mean normalized temperature dependence curves based on the data in Table 2, for seven evergreen broad-leaved trees, 13 deciduous broad-leaved trees and three conifers (Figure 7). Leuning’s (2002) mean curve for 59 datasets and the temperature dependence curve for cotton (Harley et al. 1992) and tobacco (Bernacchi et al. 2001), which are often used as examples of typical temperature dependence of Vcmax , are also shown in Figure 7. The entropy term (ΔS) was fixed at 650 for simplicity because no dataset enabling precise optimization of ΔS exists. The result of the comparison based on these limited datasets suggests that the activation energy is lower in evergreen broad-leaved trees and conifers than in deciduous broad-leaved trees, cotton or tobacco. This might result from acclimation of leaves growing at low temperatures. In practice this means that Vcmax of evergreen trees and conifers tends to be high at low temperatures. These tendencies suggest that spe- Maximum carboxylation rate Mean Vcmax25 of the three evergreen trees (mean and SD of 32.5 ± 8.2) was slightly lower than mean Vcmax25 (39.9 ± 10.7) of the four deciduous tree species investigated by Kosugi et al. (2003). However, there were no significant differences between the Vcmax25 of the co-occurring deciduous oak and the evergreen trees. Differences in normalized Vcmax could be associated with environmental conditions at each site. The lack of a dataset for evergreen tree species grown under various environmental conditions and the lack of information on potential species-specific differences in the kinetic parameters of Rubisco make it difficult to explain the differences between our estimated Vcmax and those of previous studies. We compared our temperature dependence curve parameters for mature leaves with those reported in previous studies (Table 2). The ΔH a (Vcmax ) of the three evergreen species ranged from 29.1 to 67.0 kJ mol – 1. The ΔH a (Vcmax ) of C. cuspidata was significantly lower than the values obtained for the other two species, perhaps indicating winter acclimation. Consequently, we constructed another temperature curve based on summer and autumn data only (June 7–October 16, 2001) (Figure 3; gray line). Estimated ΔH a (Vcmax ) then closely approximated the results obtained for the other species (56.2 kJ mol –1 ). With our experimental design, we could not distinguish between the simple temperature dependence of Rubisco activity from acclimation effects, although “apparent” temperature dependence obtained from long-term field observations would improve prediction accuracy. The ΔH a (Vcmax ) of the deciduous oak was optimized to 67.0 kJ mol – 1, which is close to the value of the average curve derived Figure 6. Comparison of seasonal fluctuations in observed and simulated dark respiration rate (R n ) based on the parameter sets in Table 1. Values are means for three leaves. Symbols: 䊊 = previous year leaves; 䊉 = current year leaves; and 䊉 = simulated value using a constant R n25 . Error bars are ± 1 standard deviation. TREE PHYSIOLOGY VOLUME 26, 2006 SEASONAL FLUCTUATIONS IN LEAF GAS EXCHANGE 1181 Figure 7. Mean normalized temperature dependence curve based on the data in Table 2 for seven evergreen broad-leaved trees, 13 deciduous broad-leaved trees and three conifers. The mean curve for 59 datasets by Leuning (2002) and the temperature dependence curves for cotton (Harley et al. 1992) and tobacco (Bernacchi et al. 2001) are shown for comparison. cies of similar life type and biome have similar temperature dependence curves as a result of adaptation to the environmental conditions of their habitat. When there is insufficient information available to fully describe actual temperature dependence, we propose that parameters associated with an averaged temperature dependence curve for each life style or species (Figure 7; Table 2) can be used with relatively little error. However, the parameters in Table 2 for conifers may result in substantial errors because they are based on only three data sets. Seasonal fluctuations in observed Vcmax and simulated Vcmax based on the parameter sets listed in Table 1 are shown in Figure 8 (29,100 was used for C. cuspidata). For the three evergreen broad-leaved trees at this site, Vcmax could be estimated for a whole year based on a single set of parameters, except during the periods of leaf expansion and leaf fall in April and May. Deciduous broad-leaved trees, however, showed a postsummer decline in Vcmax at the reference temperature, as well as differences during the leaf expansion period and both of these fluctuations would need to be incorporated into models. In situ field data measured in the forest canopy under natural conditions are most reliable for revealing the exchanges in canopy leaves and hence, represent the most suitable data for the parameterization of leaf gas exchange models. However, the values, range and behavior of Vcmax detected in the field should be considered apparent rather than actual because Vcmax is determined by Rubisco activity which is directly related to the amount of Rubisco and leaf nitrogen content through complex mechanisms of nitrogen allocation (e.g., Field 1983, Evans 1987, Evans 1993, Hikosaka and Terashima 1995, Niinemets and Tenhunen 1997, Takashima et al. 2004). In addition, there are many other factors that may influence Vcmax. For example, most studies investigating gas exchange parameters ignore internal (mesophyll) conductance, even though many isotope studies have determined that Cc of tree species is around 0.1 to 0.2 CO2 mol m – 2 s – 1 or lower (e.g., von Caemmerer and Evans 1991, Loreto et al. 1992, Epron et al. 1995, Hanba et al. 1999, 2001), which means that Cc may be a major factor determining the range of apparent Vcmax for each species (Loreto 1992). The influence on apparent Vcmax and photoinhibition of patchy stomatal closure should also be considered in some sit- uations. We observed neither marked patchiness in stomatal closure, which causes severe decreases in apparent Vcmax (Takanashi et al. 2006) nor depression of electron transport rate except during the severe summer drought in 2001 (Kosugi and Matsuo, unpublished data). To evaluate the possible influences of internal conductance (g i ) and error in the estimation of R d on the estimation of apparent Vcmax , we calculated seasonal Figure 8. Comparison of seasonal fluctuations in observed and simulated maximum carboxylation rate (Vcmax ) based on the parameter sets listed in Table 1. Values are means for three leaves. Symbols: 䊊 = previous year leaves; 䊉 = current year leaves; and 䊉 = simulated value using a constant Vcmax25. Error bars are ± 1 standard deviation. TREE PHYSIOLOGY ONLINE at http://heronpublishing.com 1182 KOSUGI AND MATSUO fluctuations in Vcmax by several methods (Figure 9). The possible influence of the error in the estimation of R d is shown on the left side of Figure 9. There were differences between apparent Vcmax calculated by Equation 3 and calculated assuming that R d = R n and this assumption has been questioned in several other studies (e.g., Brooks and Farquhar 1985, Kirschbaum and Farquhar 1987, Villar et al.1994, 1995). The differences in apparent Vcmax were smaller when it was assumed that R d = 0.5R n (e.g., Niinemets and Tenhunen1997). In several other studies reporting R d derived from the A–C i curve and where R n was measured in a dark chamber (e.g., Harley et al. 1992, Dungan et al. 2003), it was suggested that a scaling ratio of about 0.5 or lower is adequate. The possible influence of gi was quite large as shown on the right side of Figure 9 (cf. Loreto et al. 1992). Values of gi of ~ 0.1 CO2 mol m – 2 s – 1 were reported for several Japanese evergreen (Hanba et al. 1999) and deciduous (Hanba et al. 2001) broad-leaved tree species. Although Vcmax is generally calculated assuming that g i = ∞. The possible influence of changing g i with leaf age might be important in our study because lower values of g i have been reported for both expanding leaves (Hanba et al. 2001) and old leaves (Loreto et al. 1994). However, Figure 9 shows that differences in g i associated with leaf age within the range reported by Hanba et al. (2001) do not explain the low Vcmax during leaf expansion and leaf fall. Stomatal conductance Figure 10 shows seasonal fluctuations in parameter m compared with the constant values optimized for each tree. Parameter m has a significant physiological meaning related to intrinsic water-use efficiency, indicating species-dependent Figure 10. Comparison of seasonal fluctuations in observed and constant stomatal coefficient (m) (䊊 = previous-year leaves; 䊉 = current-year leaves; and 䊉 = simulated values). Figure 9. Effects of daytime respiration rate (R d ) and internal conductance (g i ) on the estimation of “apparent” maximum carboxylation rate (Vcmax ). Seasonal fluctuations in Vcmax calculated by several methods are shown. Left side: 䊉 = Vcmax estimated with R d = Equation 3; 䊉 = R d = R n; and 䊊 = R d = 0.5R n. Right side: 䊉 = g i = ∞; 䊉 = g i = 0.2 CO2 mol m –2 s –1; and 䊊 = gi = 0.1 CO2 mol m –2 s –1. TREE PHYSIOLOGY VOLUME 26, 2006 SEASONAL FLUCTUATIONS IN LEAF GAS EXCHANGE regulation of water fluxes; slight fluctuations in parameter m are therefore important. Some studies have reported that m decreases in response to drought (Sala and Tenhunen 1996) and in old trees (Falge et al. 1996). The evidence suggests that different m values should be applied during periods of leaf expansion and during drought in the long-term simulation of gas exchange for all of our study species. Conclusions Seasonal fluctuations in three major leaf gas exchange parameters (R n , Vcmax and m) were investigated in co-occurring evergreen and a deciduous tree species in a temperate broad-leaved forest. For the three evergreen trees, one set of the parameters Vcmax25 and R n25 and temperature dependence curves satisfactorily estimated fluctuations in gas exchange throughout the year, except during the period of simultaneous leaf fall and leaf expansion in April and May. In the deciduous oak, post-summer declines in Vcmax25 were observed as well as changes in Vcmax25 and R n25 during the leaf expansion period. The period of leaf expansion, when leaf gas exchange characteristics differ from those of mature leaves, was shorter in the deciduous oak than in the evergreen species. The temperature dependence of the activation energy of Vcmax was investigated using both our data and data from previously published studies and was found to be smaller in evergreen broad-leaved trees than in deciduous broad-leaved trees. For all of our study species, the difference in m should be considered during periods of leaf expansion and drought. Canopy-scale flux observations at this site demonstrated changes in the photosynthetic ability and stomatal behavior characteristics of the whole canopy during periods of leaf expansion and drought (Kosugi et al. 2005). These results strongly suggest that the changes in the normalized gas exchange parameters, and thus, the impact of leaf physiology on seasonal fluctuations in gas exchange, should be incorporated in models of carbon uptake of evergreen broad-leaved forests. Acknowledgments We thank the Kansai Electric Power Co. and Kansai Environmental Engineering Co. (KANSO) for help with our field observation at Akou Power Station. We also thank Dr. Tsunahide Shidei, Dr. Makoto Tani and Dr. Shozo Shibata for their support with the project on this site and Ms. Motoko Higuchi, Dr. Satoru Takanashi, Dr. Hiroki Tanaka, Ms. Noriko Hama and Dr. Nobuhito Ohte for collecting field data. References Badger, M.R. and G.J. Collatz. 1977. Studies on the kinetic mechanisms of ribulose-1, 5-bisphosphate carboxylase and oxygenase reactions with particular reference to the effect of temperature on kinetic parameters. Carnegie Institute of Washington Year Book 76:355–361. Baldocchi, D. and T. Meyers. 1998. On using eco-physiological, micrometeorological and biogeochemical theory to evaluate carbon dioxide, water vapor and trace gas fluxes over vegetation: a perspective. Agric. For. Meteorol. 90:1–25. 1183 Ball, J.T., I.E. Woodrow and J.A. Berry. 1987. A model predicting stomatal conductance and its contribution to the control of photosynthesis under different environmental conditions. In Progress in Photosynthesis Research. Ed. I. Biggins. Martinus Nijhoff Publishers, Dordrecht, The Netherlands, pp 221–224. Bernacchi, C.J., E.L. Singsaas, C. Pimentel, Jr., A.R. Portis and S.P. Long. 2001. Improved temperature response functions for models of Rubisco-limited photosynthesis. Plant Cell Environ. 24:253–259. Brooks, A. and G.D. Farquhar. 1985. Effect of temperature on the CO 2 /O2 specificity of ribulose-1,5-bisphosphate carboxylase/oxygenase and the rate of respiration in the light. Planta 165: 397–406. von Caemmerer, S. and J.R. Evans. 1991. Determination of the average partial pressure of CO2 in chloroplast from leaves of several C3 plants. Aust. J. Plant Physiol. 18:287–305. Dreyer, E., X. Le Roux, P. Montpied, F.A. Daudet and F. Masson. 2001. Temperature response of leaf photosynthetic capacity in seedlings from seven temperate tree species. Tree Physiol. 21: 223–232. Dungan, R.J., D. Whitehead and R.P. Duncan. 2003. Seasonal and temperature dependence of photosynthesis and respiration for two co-occurring broad-leaved tree species with contrasting leaf phenology. Tree Physiol. 23:561–568. Epron, D., D. Godard, G. Cornic and B. Genty. 1995. Limitation of net CO2 assimilation rate by internal conductances to CO2 transfer in the leaves of two tree species (Fagus sylvatica L. and Castanea sativa Mill.). Plant Cell Environ. 18:43–51. Evans, J.R. 1987. The dependence of quantum yield on wavelength and growth irradiance. Aust. J. Plant Physiol. 14:69–79. Evans, J.R. 1993. Photosynthetic acclimation and nitrogen partitioning within a lucerne canopy. I. Canopy characteristics. Aust. J. Plant Physiol. 20:55–67. Falge, E., W. Graber, R. Siegwolf and J.D. Tenhunen. 1996. A model of the gas exchange response of Picea abies to habitat conditions. Trees 10:277–287. Farquhar, G.D., S. von Caemmerer and J.A. Berry. 1980. A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta 149:78–90. Field, C. 1983. Allocating leaf nitrogen for the maximization of carbon gain: lead age as a control on the allocation program. Oecologia 56:341–347. Grossman-Clarke, S., B.A. Kimball, D.J. Hunsaker et al. 1999. Effects of elevated atmospheric CO2 on canopy transpiration in senescent spring wheat. Agric. For. Meteorol. 93:95–109. Han, Q., T. Kawasaki, S. Katahara, Y. Mukai and Y. Chiba. 2003. Horizontal and vertical variations in photosynthetic capacity in a Pinus densiflora crown in relation to leaf nitrogen allocation and acclimation to irradiance. Tree Physiol. 23:851–857. Hanba, Y.T., S. Miyazawa and I. Terashima. 1999. The influence of leaf thickness on the CO2 transfer conductance and leaf stable carbon isotope ratio for some evergreen tree species in Japanese warm-temperate forests. Funct. Ecol. 13:632–639. Hanba, Y.T., S. Miyazawa, H. Kogami and I. Terashima. 2001. Effects of leaf age on internal CO2 transfer conductance and photosynthesis in tree species having different types of shoot phenology. Aust. J. Plant Physiol. 28:1075–1084. Harley, P.C. and D.D. Baldocchi. 1995. Scaling carbon dioxide and water vapor exchange from leaf to canopy in a deciduous forest. I. Leaf model parametrization. Plant Cell Environ. 18:1146–1156. Harley, P.C., J.D. Tenhunen and O.L. Lange. 1986. Use of an analytical model to study limitations on net photosynthesis in Arbutus unedo under field conditions. Oecologia 70:393–401. TREE PHYSIOLOGY ONLINE at http://heronpublishing.com 1184 KOSUGI AND MATSUO Harley, P.C., R.B. Thomas, J.F. Reynolds and B.R. Strain. 1992. Modeling photosynthesis of cotton grown in elevated CO2. Plant Cell Environ. 15:271–282. Hikosaka, K. and I. Terashima. 1995. A model of the acclimation of photosynthesis in the leaves of C3 plants to sun and shade with respect to nitrogen use. Plant Cell Environ. 18:605–618. Jordan, D.B. and W.L. Ogren. 1984. The CO2 /O2 specificity of ribulose-1,5-bisphosphate carboxylase/oxygenase. Planta 161: 308–313. Kirschbaum, M.U.F. and G.D. Farquhar. 1984. Temperature dependence of whole-leaf photosynthesis in Eucalyptus pauciflora Sieb. Ex Spreng. Aust. J. Plant Physiol. 11:519–538. Kirschbaum, M.U.F. and G.D. Farquhar. 1987. Investigation of the CO2 dependence of quantum yield and respiration in Eucalyptus pauciflora. Plant Physiol. 83: 1032–1036. Kosugi, Y., S. Shibata and S. Kobashi. 2003. Parameterization of the CO2 and H2O gas exchange of several temperate deciduous broadleaved trees at the leaf scale considering seasonal changes. Plant Cell Environ. 26:285–301. Kosugi, Y., H. Tanaka, S. Takanashi, N. Matsuo, N. Ohte, S. Shibata and M. Tani. 2005. Three years of carbon and energy fluxes from Japanese evergreen broad-leaved forest. Agric. For. Meteorol. 132: 329–343. Lai, C.T., G. Katul, J. Butnor, M. Siqueira, D. Ellsworth, C. Maier, K. Johnsen, S. Mckeand and R. Oren. 2002. Modeling the limits on the response of net carbon exchange to fertilization in a south-eastern pine forest. Plant Cell Environ. 25:1095–1119. Le Roux, X., S. Grand, E. Dreyer and F.A. Daudet. 1999. Parameterization and testing of a biochemically based photosynthesis model for walnut (Juglans regia) trees and seedlings. Tree Physiol. 19:481–492. Leuning, R. 1995. A critical appraisal of combined stomatal-photosynthesis model for C3 plants. Plant Cell Environ. 18:339–355. Leuning, R. 2002. Temperature dependence of two parameters in a photosynthesis model. Plant Cell Environ. 25:1205–1210. Lloyd, J., S.C. Wang, J.M. Styles, D. Batten, R. Riddle, C. Turnbull and C.A. McConchie. 1995. Measuring and modeling whole-tree gas exchange. Aust. J. Plant Physiol. 22:987–1000. Loreto, F., P.C. Harley, G.D. Marco and T.D. Sharkey. 1992. Estimation of mesophyll conductance to CO2 flux by three different methods. Plant Physiol. 98:1437–1443. Loreto, F., P.C. Harley, G.D. Marco and T.D. Sharkey. 1994. Measurements of mesophyll conductance, photosynthetic electron transport and alternative electron sinks of field grown wheat leaves. Photosynth. Res. 41:397–403. Lohammer, T., S. Larsson, S. Linder and S.O. Falk. 1980. FAST-Simulation models of gaseous exchange in Scots Pine. Ecol. Bull. 32:505–523. Niinemets, U. and J.D. Tenhunen. 1997. A model separating leaf structural and physiological effects on carbon gain along light gradients for the shade-tolerant species Acer saccharum. Plant Cell Environ. 20:845–866. Sala, A. and J.D. Tenhunen. 1996. Simulations of canopy net photosynthesis and transpiration in Quercus ilex L. under the influence of seasonal drought. Agric. For. Meteorol. 78:203–222. Sharpe, P.J.H. and D.W. DeMichele. 1977. Reaction kinetics of poikilotherm development. J. Theor. Biol. 64:649–670. Takanashi, S., Y. Kosugi, N. Matsuo, M. Tani and N. Ohte. 2006. Patchy stomatal behavior in broad-leaved trees grown in different habitats. Tree Physiol. In press. Takashima, T., K. Hikosaka and T. Hirose. 2004. Photosynthesis or persistence: nitrogen allocation in leaves of evergreen and deciduous Quercus species. Plant Cell Environ. 27:1047–1054. Tanaka, K., Y. Kosugi and A. Nakamura. 2002. Impact of leaf physiological characteristics on seasonal variation in CO2, latent and sensible heat exchanges over a tree plantation. Agric. For. Meteorol. 114:103–122. Tenhunen, J.D., A. Sala, P.C. Harley, R.L. Doughherty and J.F. Reynolds. 1990. Factors influencing carbon fixation and water use by Mediterranean sclerophyll shrubs during summer drought. Oecologia 82:381–393. Villar, R., A.A. Held and J. Merino. 1994. Comparison of methods to estimate dark respiration in the light in leaves of two woody species. Plant Physiol. 105: 167–172. Villar, R., A.A. Held and J. Merino. 1995. Dark leaf respiration in light and darkness of an evergreen and a deciduous plant species. Plant Physiol. 107: 421–427. Walcroft, A.S., D. Whitehead, W.B. Silvester and F.M. Kelliher. 1997. The response of photosynthetic model parameters to temperature and nitrogen concentration in Pinus radiata D. Don. Plant Cell Environ. 20:1338–1348. Wang, K.W., S. Kellomaki and K. Laitinen. 1996. Acclimation of photosynthetic parameters in Scots pine after three years exposure to elevated temperature and CO2. Agric. For. Meteorol. 82: 195–217. Weber, J.A., T.W. Jurik, J.D. Tenhunen and D.M. Gates. 1985. Analysis of gas exchange in seedlings of Acer saccharum: integration of field and laboratory studies. Oecologia 65:338–347. Williams, T.G. and L.B. Flanagan. 1998. Measuring and modeling environmental influences on photosynthetic gas exchange in Sphagnum and Pleuozium. Plant Cell Environ. 21:555–564. Wilson, K.B., D.D. Baldocchi and A.P.J. Hanson. 2000. Spatial and seasonal variability of photosynthetic parameters and their relationship to leaf nitrogen in a deciduous forest. Tree Physiol. 20: 565–578. Wilson, K.B., D.D. Baldocchi and A.P.J. Hanson. 2001. Leaf age affects the seasonal pattern of photosynthetic capacity and net ecosystem exchange of carbon in a deciduous forest. Plant Cell Environ. 24:571–583. Wohlfahrt, G., M. Bahn, E. Haubner, I. Horak, W. Michaeler, K. Rottmar, U. Tappeiner and A. Cernusca. 1999. Inter-specific variation of the biochemical limitation to photosynthesis and related leaf traits of 30 species from mountain grassland ecosystems under different land use. Plant Cell Environ. 22:1281–1296. Xu, L. and D.D. Baldocchi. 2003. Seasonal trends in photosynthetic parameters and stomatal conductance of blue oak (Quercus douglasii) under prolonged summer drought and high temperature. Tree Physiol. 23:865–877. TREE PHYSIOLOGY VOLUME 26, 2006
© Copyright 2025 Paperzz