Can we estimate forest gross primary production from leaf lifespan? A test in a young Fagus crenata forest
 Author: Koyama Kohei, Kikuzawa Kihachiro
 Organization: Koyama Kohei; Kikuzawa Kihachiro
 Publish: Journal of Ecology and Environment Volume 33, Issue3, p253~260, 01 Sep 2010

ABSTRACT
It has been well established that leaf longevity is linked to the carbon economy of plants. We used this relationship to predict leaf lifetime carbon gains from leaf lifespan, and estimated the gross primary production (GPP) of a young deciduous forest of Japanese beech (
Fagus crenata ) located in central Japan. The lightsaturated photosynthetic rates of the leaves were measured repeatedly during the growing season. We used the leaf lifespan to calculate the conversion coefficient from the lightsaturated photosynthetic rate into the realized leaf lifetime carbon gain under field conditions. The leaf turnover rate was estimated using litter traps. GPP was estimated as the product of lifetime carbon gain per unit of leaf mass, and the annual leaf turnover rate. The GPP of the forest in 2007 was estimated to be 1.2 × 10^{3} g C m^{2} y^{1}, which was within the range of previously reported GPP values of beech forests in Japan, and was close to the GPP of a European beech forest, as estimated by eddy flux measurements.

KEYWORD
Fagus crenata , functional leaf longevity , gross primary production , mean labor time , temperate forest

INTRODUCTION
The gross primary production (GPP) of forests may function as a terrestrial carbon sink against climate change (Suwa et al. 2006, Canadell et al. 2007, Luyssaert et al. 2008, Saigusa et al. 2008, Lewis et al. 2009). GPP is also a fundamental energy and material source that occurs in forest ecosystems as the result of litter input (Kurokawa and Nakashizuka 2008, Litton and Giardina 2008), and the activity of ectomycorrhiza (Druebert et al. 2009). Recent "ecological scaling" theories have outlined simple relationships predicting ecosystem carbon gain from single leaf photosynthesis (Sellers et al. 1992, Koyama and Kikuzawa 2009, 2010). However, these models are useful for the spatial integration of canopy photosynthesis. To estimate GPP, temporal integration over growing seasons is required (Thornley 2002). In service of this objective, numerical calculations (e.g. Hikosaka 2003) have been applied. However, a simple model for timeintegration over the growing period has yet to be developed.
Leaf demography and longevity has been studied extensively in association with wholeplant carbon economy (Chabot and Hicks 1982, Ackerly and Bazzaz 1995, Franklin and Agren 2002, Oikawa et al. 2004, 2008, Reich et al. 2004, Hikosaka 2005, Koyama and Kikuzawa 2008, 2009, Mediavilla and Escudero 2009, Suarez 2010). Mathematical models predicted that higher photosynthetic rates in leaves would result in higher leaf turnover rates, accompanied by reduced leaf longevity (Kikuzawa 1991, Ackerly 1999, Hikosaka 2003). In support of this notion, a negative correlation was observed to exist between leaf lifespan and the instantaneous lightsaturated leaf photosynthetic rate (Chabot and Hicks 1982, Koike 1988, Reich et al. 1991, 1997, Wright et al. 2004, Vincent 2006, Ishida et al. 2008, Nagano et al. 2009). Hence, it has been suggested that this lifespancarbon gain relation can be employed in the development of a quantitative model of ecosystem productivity (SauraMas et al. 2009). However, to the best of our knowledge, no study has yet successfully estimated GPP via this relationship.
Kikuzawa and Lechowicz (2006) recently proposed a scheme that predicts GPP over the growing period, using a simple model. In essence, they demonstrated that if the leaf lifespan is a function of the daily carbon gain of that leaf, then the leaf lifetime carbon gain can be viewed as an inverse function of the leaf lifespan. In this paper, we applied this method to estimate the GPP of a temperate deciduous forest composed of Japanese beech (
Fagus crenata ). We combined leaf demography census, periodical measurements of lightsaturated photosynthetic rate of sample leaves, and litter measurements, and converted the lightsaturated photosynthetic rate into the leaf lifetime carbon gain value, by using the leaf lifespan.> Model
Our model is based on the model described by Kikuzawa and Lechowicz (2006), with some slight modifications. The relevant symbols are listed in Table 1. The two principal assumptions here are as follows: 1) leaf lifespan is determined so as to maximize wholeplant carbon gain, in concert with the numerical model of Kikuzawa (1991) (for a more elaborate version of this model, see Takada et al. 2006). This model was also supported by the empirical data (Kikuzawa and Ackerly 1999). 2) Stand leaf biomass is constant for closed steadystate canopies (Tadaki 1986, 1991). GPP (g C m^{2} y^{1}) is the sum of leaf photosynthesis of all the leaves within a stand:
B is the standing leaf biomass (g C/m^{2}) andA _{gross} (g C g C^{1} d^{1}) is the daily gross leaf photosynthetic rate averaged over the leaf lifespan and over the stand.f (d/y) is the number of days on which photosynthesis was conducted per year (i.e. growing season length). Eq. 1 can be rewritten by multiplying the right side byL _{f}/L _{f} (= 1), whereL _{f} denotes the functional leaf longevity (Kikuzawa and Lechowicz 2006):Functional leaf longevity (
L _{f}) is the number of days on which the leaf actually carries out photosynthesis within its lifespan. For evergreen leaves in seasonal environments,L _{f} is the difference between the observed leaf lifespan (total number of days, for which one leaf exists) and the length of the unfavorable period for photosynthesis (i.e. winter, dry season etc.). In the case of wet tropics, as well as the case of deciduous leaves in a seasonal climate,L _{f} equals the observed leaf lifespan. Based on the assumption of constant stand leaf biomass, the ratio of stand leaf biomassB (g C/m^{2}) andL _{f} (days) is assumed to be equal to the daily leaf turnover rate during the favorable season, which is denoted byp (g C m^{2} d^{1}):The product of functional leaf longevity (
L _{f}) and daily leaf photosynthetic rate averaged over the leaf lifespan (A _{gross}) is equivalent to the leaf lifetime gross carbon gain, which is denoted by G_{g} (g C/g C):By substituting Eq. 3 and Eq. 4 into Eq. 2, GPP can be rewritten as:
The product of
p andf is an annual leaf carbon turnover rate (i.e. annual leaf death rate is assumed to be equal to the production rate). In essence, Eq. 5 expresses GPP as the product of the "amount of carbon produced per unit leaf carbon through its lifecycle (G _{g})" and the "annual leafcarbon lifecycle turnover rate (p f )". Those equations are generic, in the sense that they may be applicable to ecosystems with different lengths of growing season (wet tropics, dry tropics, temperate, etc.), as well as with different leaf lifespans (deciduous, evergreens which has leaf lifespan longer than one year, etc.).G _{g} is the sum of leaflifetime net carbon gainG _{n} (g C/g C) and leaflifetime respirationR (g C/g C):G _{n} can be expressed as the product of the three following parameters:A _{max} (g C g C^{1} s^{1}) is an instantaneous lightsaturated net photosynthetic rate averaged over the leaf lifespan, andm (s/d) is a conversion coefficient from lightsaturated rate into the actual net photosynthetic rate averaged over its lifespan realized at the site, and referred to as the "mean labor time" (Kikuzawa et al. 2004).m is a combination of all the factors that decide the actual photosynthetic rates of leaves under each field condition, including the diurnal solar cycle, weather, selfshading, and midday depression. In a previous study (Kikuzawa et al. 2004), all of those factors were measured in order to calculate m. In this study, we employed an optimal leaf longevity theory (Kikuzawa 1991) to estimatem . This is based on the notion that, when an optimal leaf lifespan is a function of the leaf carbon budget, the leaf carbon gain should be an inverse function of the leaf lifespan. Kikuzawa and Lechowicz (2006) then proposed an equation to calculatem from the leaf lifespan:m is conventionally expressed as seconds per day (Kikuzawa et al. 2004).A _{max (0)} (g C g C^{1} s^{1}) is the massbased lightsaturated net photosynthetic rate of a fully expanded leaf. b (days) is the time at which the lightsaturated photosynthetic rate of each leaf becomes zero, assuming the linear decline of lightsaturated photosynthetic rate with respect to time.C (g C/g C) is the construction cost of a leaf, defined as the amount of carbon material required to construct one unit of leaf carbon mass. In this study,C was assumed to be 1.53 g C/g C (Griffin 1994). The substitution of Eqs. 6, 7 and 8 into Eq. 5 resulted in the following:m is conventionally expressed as seconds per day (Kikuzawa et al. 2004).A _{max (0)} (g C g C^{1} s^{1}) is the massbased lightsaturated net photosynthetic rate of a fully expanded leaf. b (days) is the time at which the lightsaturated photosynthetic rate of each leaf becomes zero, assuming the linear decline of lightsaturated photosynthetic rate with respect to time.C (g C/g C) is the construction cost of a leaf, defined as the amount of carbon material required to construct one unit of leaf carbon mass. In this study,C was assumed to be 1.53 g C/g C (Griffin 1994). The substitution of Eqs. 6, 7 and 8 into Eq. 5 resulted in the following:MATERIALS AND METHODS
> Species and site
Fagus crenata Blume is a latesuccessional deciduous canopy tree found in mountainous forests throughout Japan, and is a major component of those forests (Okaura and Harada 2002). The study site employed herein was the Ishikawa Prefectural Forest Experiment Station at the foot of Mt. Hakusan in central Japan (36^{o}25' N, 136^{o}38' E, 200 m a.s.l.). The mean annual temperature and mean annual precipitation were 13.0℃ and 2,438 mm, respectively (Ishikawa Prefectural Forest Experiment Station 2008). The stand was a 14 yearold plantation of Japanese beech (Fagus crenata Blume). There were 130 trees in an area of 171 m^{2}, and the heights of the trees were 45 m and the mean diameter at breast height was 4 cm in 2006. The canopy was fairly closed, and there was almost no understory vegetation.> Measurement of photosynthesis
In 2007, a scaffolding tower that reached the canopy was constructed. A total of 16 leaves were selected from four trees from the entire position of the canopy (0.84.7 m from the ground). Sample leaves were selected from different heights, so as to represent the entire leaf population (Fig. 1). Photosynthesis measurements were repeated for the same leaves 613 times, from May to November of 2007. The number of measurements was smaller for the leaves that fell earlier. The photosynthetic rate of each leaf was measured using a portable infrared gas analyzer (LI6400; LICOR, Lincoln, NE, USA). Photosynthetic photon flux density (PPFD) 1,500 μmol m^{2} s^{1} was supplied until the equilibration of the leaf with an LED light source (LI640002B; LICOR) within the chamber. The equilibration required 1020 minutes for the uppercanopy leaves, and 2030 minutes for the lowercanopy leaves. CO_{2} concentrations inside the chamber were controlled at 350 ppm. We did not control the air temperature inside the chamber, which ranged from 2230℃. Measurements were conducted from 8 to 11 a.m. each day. After each measurement, PPFD was changed to zero, and after equilibration (1020 minutes), the dark respiration rate was measured. The gross lightsaturated photosynthetic rate of the leaves on each day is acquired as the sum of the lightsaturated net photosynthetic rate and the dark respiration rate.
For each of the monitored leaves, we sampled three adjacent leaves on the same horizontal branch in the summer. The areas of those leaves were determined via the analysis of digitallyscanned images with Adobe Photoshop Elements (Adobe Systems Inc., San Jose, CA, USA). The leaf samples were then ovendried at 70℃ to a constant weight and weighed. The mean leaf mass per area of each of those three leaves was used as surrogate values for each of the monitored leaves. Using them, areabased photosynthetic and respiration rates for each leaf were converted into massbased values. Those leaf masses and the leaf litter were converted as follows: 1 g dry leaf = 0.5 g C (Yamasaki and Kikuzawa 2003), in order to obtain the carbonbased values.
> Leaf demography
Sample leaves for the photosynthesis measurements were monitored from leaf emergence until fall at 723day intervals in 2007. Functional leaf longevity (
L _{f}) was calculated as the difference between leaf emergence and fall, as those leaves were detected only during the growing season. Ten litter traps (0.17 m^{2} for each) were set on the stand floor such that they were evenly distributed throughout the entire forest. The litters were collected every month from April to November, and separated into leaves and branches. The leaves were ovendried at 70℃ to a constant weight and weighed. The annual leaf production rate (p f ) is assumed to be equal to the annual leaf fall from the beech trees, assuming that the loss by herbivory was negligible.> Data analysis for photosynthesis
The lifetime maximum instantaneous lightsaturated net photosynthetic rate is used as the
A _{max (0)} for each leaf (Kikuzawa 1991). TheA _{max} was calculated as the lifetime average of the net lightsaturated photosynthetic rate for each leaf. The slope of the linear regression of the lightsaturated photosynthetic rate with time was calculated and tested via Microsoft Excel (t test). The time at which the lightsaturated photosynthetic rate of each leaf becomes zero (b ) was estimated via the regression line.RESULTS AND DISCUSSION
The results were summarized in Table 2.
F. crenata is a representative species of simultaneous leaf emergence and fall (Kikuzawa 1983, 2003). At the study site, all of the leaves emerged almost simultaneously in late April, and the leaf fall was concentrated in October or November. The mean functional leaf lifespan (L _{f}) of the sample leaves was 183 days. In addition to those sampled leaves, we also measured a total of 1,066 leaves from the entire position of the canopy, to determine the lifespans of the sample leaves. The mean lifespan of the additional 1,066 leaves was 184 days, which was similar to those of the monitored leaves, thereby indicating that the lifespans of those sample leaves represented those of the entire canopy.The lightsaturated net photosynthetic rate of each leaf reduced in an almost linear fashion after full expansion (Fig. 2). The slope of the linear regression of lightsaturated photosynthetic rate with time was found to be significant for all the leaves (
P < 0.05,r ^{2} = 0.520.96,N = 16). The decline of photosynthetic capacity with increasing leaf age has been observed commonly among woody species (Koike 1990, Kitajima et al. 2002, Kikuzawa 2003, Mediavilla and Escudero 2003, Ito et al. 2006, Han et al. 2008, Kikuzawa et al. 2009, Reich et al. 2009). However, the decline of photosynthetic capacity was frequently not precisely linear, but rather curvilinear (cf. Koike 1990). Hence, the linear regression was a simplified approximation, and more elaborated models should incorporate the curvilinear change characteristics of photosynthetic capacity.The estimated GPP of the forest (1.2 × 10^{3} g C m^{2} y^{1}) was within the GPP range of beech forests in central Japan (0.57 × 10^{3}1.56 × 10^{3} g C m^{2} y^{1}, assuming that 1 g dry matter = 0.44 g C) previously reported by Kakubari (1991), which employs the conventional allometric method. This was also similar to the GPP range for a 30year old European beech forest (1,0001,300 g C m^{2} y^{1}), which was estimated via eddy flux measurements (Granier et al. 2000). To the best of our knowledge, this is the first application of the leaf lifespan theory to the calculation of GPP. Although our method may be less accurate than the present eddy flux measurements, our theoretical approach produced a value proximate to the actual one.
The principal assumption of this model is that leaf longevity is determined so as to maximize wholeplant carbon gains in concert with the numerical model developed by Kikuzawa (1991). Many previous reports conducted in the tropics have corroborated the carboneconomy model (e.g. Reich et al. 1991, 1999, 2004, Vincent 2006). However, leaf longevity should be susceptible to the effects of strong seasonality in temperate regions (Reich et al. 2004, Vincent 2006, Koyama and Kikuzawa 2008). Nonetheless, two lines of evidence appear to suggest that leaf longevity in temperate regions functions as an indicator of the carbon budget of a plant. Firstly, temperate species were included in the worldwide leaf economic spectrum (Reich et al. 1997, 1999, Ackerly and Reich 1999, Wright et al. 2004, He et al. 2009). If the effect of seasonality is predominant, species in the temperate region should be outliers within those spectra; however, this is not the case. Secondly, variations in leaf longevity in temperate regions have also been predicted by carbon economy models (Hikosaka 2005, Oikawa et al. 2006, 2008). Even within the same temperate forest, deciduous and evergreen leaf habits coexist in microclimates, with deciduous trees tending to be found in more productive environments (Monk 1966, Chabot and Hicks 1982, Kikuzawa 1984, Koyama and Kikuzawa 2008). Leaf longevity in temperate regions should, assumedly, reflect the carbon economy of plants. Hence, the method presented herein should be globally applicable, based on the generality of the leaf economic spectra.

4. Canadell J.G, Le Quere G, Raupach M.R, Field C.B, Buitenhuis E.T, Ciais P, Conway T.J, Gillett N.P, Houghton R.A, Marland G 2007 Contributions to accelerating atmospheric CO2 growth from economic activity carbon intensity and efficiency of natural sinks [Proc Natl Acad Sci U S A] Vol.104 P.1886618870

[Table 1.] List of the parameters

[Fig. 1.] Vertical distribution of the monitored leaves within the canopy. Each closed circle indicates one monitored leaf (N = 16). The broken line depicts the vertical distribution of all the leaves in the canopy (estimated by spotmeasurements using a vertical measuring pole at 25 places within the stand, N = 235 leaves). The vertical axis shows the relative position (i.e. cumulative leaf numbers of each kind of the leaves from top to that height, divided by the total number of leaves of that type. The circles were distributed along the broken line. This indicates that each monitored leaf (e.g. the middle leaf among the monitored leaf) exists at a similar height as that of the middle leaf throughout the entire canopy.

[Table 2.] The estimated parameters in 2007

[Fig. 2.] Decline of lightsaturated photosynthetic rate. Each panel shows one leaf. Each open circle indicates one measurement. Two representative leaves from the upper and lower canopy are shown. Solid lines show significant linear regressions (P < 0.05).