AOBPreview originally published online on April 8, 2004
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Annals of Botany 93: 699-710, 2004
© 2004 Annals of Botany Company
Effect of 3D Nitrogen, Dry Mass per Area and Local Irradiance on Canopy Photosynthesis Within Leaves of Contrasted Heterogeneous Maize Crops
1 Unité Mixte de Recherche INRA-INAPG Environnement et Grandes Cultures, BP 01, 78850 Thiverval-Grignon, France and 2 INRA Unité Agropédoclimatique de la zone caraïbe, Domaine Duclos, 97170 Petit-Bourg, Guadeloupe, France
* For correspondence. E-mail Jean-Louis.Drouet{at}grignon.inra.fr
Received: 24 October 2003; Returned for revision: 17 November 2003; Accepted: 13 February 2004; Published electronically: 8 April 2004
| ABSTRACT |
|---|
|
|
|---|
Background and Aims Nitrogen partitioning within stands has been described fairly comprehensively, especially for C3 plants in dense stands where the horizontal heterogeneity of foliage distribution is relatively small. Nitrogen has been shown to be distributed vertically and in parallel to light, maximizing carbon assimilation and stand productivity. Conversely, row crops such as maize (C4 plants) are characterized by strong horizontal heterogeneity of foliage distribution, and a three-dimensional (3D) approach is required to investigate the combined effect of spatial distribution of nitrogen and light on canopy photosynthesis.
Model The 3D geometry of maize canopies was modelled with varying densities and at different developmental stages using plant digitizing under field conditions. For lamina parts, photosynthesis was measured and nitrogen content per unit area (Na) was described from analysis of nitrogen content per unit mass (Nm) and dry mass per unit area (Ma). Hyperbolic relationships between photosynthesis at irradiance saturation (Pmax) and Na were established as well as a linear relationship between dark respiration (Rd) and Na, whereas quantum efficiency (
) was found to be independent of Na.
Key Results and Conclusions Nm, Ma and Na were shown to change over time vertically (i.e. between laminae), which has been largely reported previously, and horizontally (i.e. within laminae), which has scarcely been described previously. Even if Ma played a major role in Na, a strong relationship between Na and Ma could not be demonstrated, whereas several previous studies have found that Na was essentially related to Ma rather than Nm. From simulations of radiative exchange using a 3D volume-based approach and lamina photosynthesis using a hyperbola, it was shown that real patterns of Na partitioning could increase daily crop photosynthesis by up to 8 % compared with uniform patterns of Na, especially for the earliest stages of stand development.
Key words: 3D plant architecture, heterogeneous crop, dry mass per unit area, irradiance, lamina, nitrogen content per unit area, nitrogen content per unit mass, maize, photosynthesis, virtual plant, Zea mays.
| INTRODUCTION |
|---|
|
|
|---|
Nitrogen supply affects plant growth and productivity by altering both leaf area and photosynthetic capacity (Novoa and Loomis, 1981). About 75 % of leaf nitrogen is involved in the photosynthetic processes (Evans, 1989) and its partitioning inside canopies has been largely described (see Grindlay, 1997, for review). Several authors have observed vertical gradients of leaf nitrogen in crops (Shiraiwa and Sinclair, 1993; Connor et al., 1995; Dreccer et al., 2000), herbaceous plants (Field, 1983; Hirose and Werger, 1987, Lemaire et al., 1991) and trees (DeJong and Doyle, 1985; Hollinger, 1996). Leaf nitrogen is distributed in parallel to light distribution (Anten et al., 1995; Kull and Jarvis, 1995), which corresponds to the optimal distribution of leaf nitrogen that maximizes carbon assimilation and crop productivity (Mooney and Gulmon, 1979; Field, 1983; Gutschick and Wiegel, 1988; Dreccer et al., 2000). Chen et al. (1993) proposed an alternative theory called coordination where canopy nitrogen is distributed so that a balance is maintained between two processes depending on nitrogen content: the Rubisco-limited rate of carboxylation and the electron transport rate of carboxylation. They found that carbon assimilation gain was similar, using either the theory of coordination or optimization (Hirose and Werger, 1987). These vertical gradients of nitrogen were essentially observed within dense stands and the combined effect of light and nitrogen on canopy photosynthesis was assessed by using a vertical description of these two variables (i.e. by canopy layer, a 1D approach). However, studies within heterogeneous canopies have been scarce, especially in row crops such as maize, sorghum, wheat and sunflower (see Grindlay, 1997), and have seldom dealt with the relationships between nitrogen and light gradients by using a 3D approach (e.g. in maize, Drouet and Bonhomme, 1999). The effect of canopy heterogeneity on light interception has been assessed previously (Drouet et al., 1999), but the combined effect of spatial distribution of nitrogen and light on canopy photosynthesis within row crops has not been investigated. Since light interception is an area-based phenomenon (Biscoe and Gallagher, 1978), leaf photosynthesis has more often than not been simulated from leaf nitrogen expressed per unit leaf area, Na (Field and Mooney, 1986; Sinclair and Horie, 1989). Na is the product of leaf nitrogen per unit mass, Nm, and dry mass per unit area, Ma, and numerous studies have dealt with relationships between Na and Ma, and between Na and Nm. Many authors have shown that changes in Na were related to changes in Ma (i.e. anatomical changes, Gulmon and Chu, 1981; DeJong et al., 1989; Hirose et al., 1989; Ellsworth and Reich, 1993; Rosati et al., 2000; Le Roux et al., 2001; Meir et al., 2002), whereas various results have been found with changes in Na related to Nm (i.e. chemical changes, Kull and Niinemets, 1993; Grassi et al., 2002). Moreover, results concerning nitrogen partitioning related to irradiance and concerning relationships between Na, Nm and Ma have been described predominantly on C3 plants and rarely on C4 plants such as maize. Those two groups of plants differ with regard to leaf nitrogen content and photosynthesis: leaf nitrogen content is lower for C4 plants than for C3 (Lemaire and Gastal, 1997) and nitrogen use efficiency is higher for C4 plants than for C3 (Gosse et al., 1986; Sage and Pearcy, 1987; Sinclair and Muchow, 1999).
In an earlier paper (Drouet and Bonhomme, 1999), it was shown that local leaf irradiance plays a major role in leaf nitrogen partitioning within maize crops. In the present study, two questions are addressed: (1) do nitrogen gradients within and between laminae in relation to local irradiance increase canopy photosynthesis compared with uniform nitrogen partitioning, and (2) what is the source of the variations in nitrogen per unit area, Na: changes in nitrogen concentration, Nm, or in dry mass per unit area, Ma, or both?
| MATERIALS AND METHODS |
|---|
|
|
|---|
Experimental design
Two field experiments were carried out in Grignon (France, 48°N, 2°E) using the maize (Zea mays L.) hybrid Déa. For the first one, maize was sown in the early summer of 1996 at two initial densities: 20 plants m2 (D-density) and 10 plants m2 (d-density). Each trial area comprised 40 rows, 50 m long. Two open-thinned plots at low density were obtained by removing plants during stem elongation from the area of initial density 10 plants m2, as follows. (1) Three plants out of four were removed from within each row, which resulted in a final density of 2·5 plants m2 (d10
2·5-density), and (2) one row out of two was removed as well as nine plants out of ten within the remaining row, which resulted in a final density of 0·5 plants m2 (d10
0·5-density). For the second experiment, maize was sown in the early summer of 1999 at two initial densities: 30 plants m2 and 10 plants m2. One plot at low density was obtained by removing plants during stem elongation from one plot of initial density 10 plants m2. Three plants out of four were removed within each row, which resulted in a final density of 2·5 plants m2. In both experiments, one amount of mineral nitrogen (5·5 g m2) was applied at sowing. The plots were weeded and plants were kept free of water stress by liberal drip irrigation.
Description of plant structure, lamina nitrogen and lamina photosynthesis
In the first experiment, the 3D geometric structure of the plants was measured with a magnetic digitizing device (Polhemus, 1993; see Sinoquet and Rivet, 1997, and Drouet, 2003, for more details). For each plant, the coordinates along the axis of the stem and the midrib of each lamina were recorded. The number of points per axis varied from ten to 30 according to the length and the curvature of the organ. To examine densities of 10 and 20 plants m2, data were obtained from 20 plants (four rows with five plants per row) at three stages of development: beginning of stem elongation (60 days after sowing, DAS), end of stem elongation (74 DAS) and post-silking (90 DAS). In both low density plots (2·5 and 0·5 plants m2), data were measured at post-silking on 12 plants (four rows with three plants per row). For each plot, the corresponding leaf area index (LAI) is shown in Table 1. Measurements were taken in the morning to minimize possible effects of wind, water stress and solar position.
|
In the first experiment, total nitrogen analyses were carried out on the six central plants of each plot (two rows with three plants per row). To study nitrogen partitioning between laminae, the even-numbered laminae from 6 to 14 were analysed. Nitrogen partitioning within laminae was assessed from analysis of nitrogen content on laminae 6 (located at the bottom of the canopy), 10 (in the middle of the canopy) and 14 (at the top of the canopy). The velum of each lamina was separated from the midrib and then subdivided into two or three parts of equal length, according to the length of the lamina. Dry mass was determined after oven drying at 70 °C for at least 2 d. After crushing, total nitrogen (mineral and organic) content per unit mass (Nm) was determined with a carbonnitrogen analyser (ANA 1500 CN; Thermoquest, Les Ulis, France) based on the Dumas method (dry process method, Dumas, 1831).
In both experiments, carbon dioxide assimilation was measured at the lamina level by using a LI-6400 portable photosynthesis system (LI-COR, Inc., NB, USA). Photo synthesisirradiance response curves were obtained by measuring photosynthesis successively at 500, 200, 150, 100, 50, 500, 1000, 1500 and 2000 µmol PAR m2 s1. Lamina photosynthesis at irradiance saturation was measured at 2000 µmol PAR m2 s1 (P2000). Samples were taken from the bottom (i.e. lamina 6) to the top (i.e. lamina 14) of the canopy, at several dates around silking within the plots at densities of 20 and 10 plants m2 in 1996, and within the plots at densities of 30, 10 and 2·5 plants m2 in 1999 (Table 2). All measurements of lamina photosynthesis were carried out between 0800 and 1100 h (see Bunce, 1990). Total nitrogen analyses were carried out on all samples where photosynthesis was measured.
|
Reconstruction of 3D shoot structure and restitution of lamina nitrogen
The reconstruction of the 3D shoot structure (Fig. 1A, B) and the restitution of lamina nitrogen partitioning (Fig. 1C, D) has been described in Drouet and Bonhomme (1999) and Drouet (2003). Dry mass per unit lamina area (Ma) was evaluated from dry mass and lamina area estimated from allometric relationships between lamina length and maximal width (Bonhomme and Varlet-Grancher, 1978). Nitrogen content per unit area (Na) was determined from Nm and Ma. Each plant was geometrically represented by a set of about 1000 triangles and individual values of Nm, Ma and Na were assigned to each triangle.
|
Simulation of light distribution and canopy photosynthesis
A 3D volume-based version of the light transfer model RIRI (Radiation Interception in Row Intercropping, Sinoquet and Bonhomme, 1992) was used to calculate irradiance distribution inside the canopies. It had previously been validated using radiation measurements for several crops, especially maize canopies (see Sinoquet and Bonhomme, 1989, 1992). This model is based on the turbid-medium analogy. In this method, the canopy structure is abstracted by an array of 3D cells (0·1 m wide), which may contain foliage or be empty. For each canopy cell, the lamina area density, the lamina angle distribution (not shown) and the lamina nitrogen content (Fig. 1C, D) are calculated from the area, the orientation and the nitrogen content of the triangles (Fig. 1A, B) (B. Andrieu, pers. comm.). The model deals with direct and diffuse incident radiation and scattered radiation, which makes it possible to obtain, within each cell of the canopy, the instantaneous lamina irradiance for sunlit and shaded lamina area (for more details, see Sinoquet and Bonhomme, 1992; Drouet, 1998). For each plot, simulations were performed in the photosynthetic active radiation waveband (PAR, 400700 nm) from six values of daily global radiation (Table 3). Instantaneous values of direct radiation and diffuse radiation were simulated for each time step (0·1 h) according to Spitters et al. (1986) and J.-M. Allirand ( pers. comm.). For each cell of the canopy, the daily average lamina irradiance (Id,a) was computed from the instantaneous values (Fig. 1E, F). The model relies upon the optical properties of the foliage. We found that reflectance and transmittance were independent of Na, at least during the period studied (data not shown). Consequently, only one average value of reflectance/transmittance was used (0·07). Soil reflectance was set equal to 0·10 (for more details, see Drouet, 1998).
|
Lamina photosynthesis in response to irradiance was simulated by using a rectangular hyperbola (Chartier, 1966, 1969; Thornley, 1976):
Pn,i,a(I) = (Pmax
I)/(Pmax +
I) Rd(1)
where Pn,i,a is net instantaneous photosynthesis per unit lamina area, I is PAR, Pmax is photosynthesis at irradiance saturation,
is quantum efficiency and Rd is dark respiration.
Pachepsky et al. (1996) indicated that such a model is quite adequate to predict quantitatively the biomass production of crops and has the advantage of requiring few parameters. These authors pointed out that it is not necessary to introduce into crop models more complicated and more sophisticated descriptions of photosynthesis. Pmax,
and Rd were fitted by non-linear regression between Pn,i,a and I.
Pmax as a function of Na was modelled by using a hyperbola (Sinclair and Horie, 1989). As photosynthesis was measured at 2000 µmol PAR m2 s1, this model was used to simulate P2000 as a function of Na:
P2000(Na) = P2000,max [2/{1 + exp[ß (Na N0)]} 1](2)
where P2000,max, ß and N0 are three parameters fitted by non-linear regression to the median P2000 value for each 0·2 g m2 range in Na (Muchow and Sinclair, 1994).
Pmax derives from P2000, and Dwyer and Stewart (1986) pointed out that the variability in P2000 is lower than that in Pmax.
and Rd were assumed to vary linearly with Na (Hirose and Werger, 1987, Muchow and Sinclair, 1994):
(Na) = a
+ b
Na(3)
Rd(Na) = aRd + bRd Na(4)
where a
, b
, aRd and bRd are parameters.
For each cell of the canopy, the net daily photosynthesis per unit lamina area (Pn,d,a, Fig. 1G, H) was computed by summing the instantaneous values (Pn,i,a). The net daily photosynthesis of the whole canopy (Pn,d,c) resulted from the sum of Pn,d,a (de Wit, 1965).
To quantify the photosynthetic benefit of real Na partitioning, we compared Pn,d,c values between real and associated hypothetical canopies. For each studied canopy, two hypothetical canopies were generated by redistributing Na from the real partitioning. In the first one, Na was redistributed uniformly along the velum by using, for each lamina, its mean Na value. In the second one, Na was redistributed uniformly within the whole canopy corresponding to the mean value of Na within the canopy.
Data were analysed (regression tests, mean comparisons, graphs) using the S-PLUS computer package (S-PLUS, 1996).
| RESULTS |
|---|
|
|
|---|
Photosynthetic capacity of laminae
The relationship between P2000 and Na was described by using the hyperbolic model proposed by Sinclair and Horie (1989) (Fig. 2A, B). Before silking, fitted values of the parameters ß and N0 (Table 4) were close to those found by Muchow and Sinclair (1994): 2·45 vs. 3·68 for ß and 0·27 vs. 0·20 for N0. Conversely, values of P2000,max strongly varied between our fitted values and those found by these authors (36·8 vs. 52·0). This difference can be explained by different climatic conditions, and consequently different photosynthetic capacity of laminae, between the two experimental sites, especially the air temperature: 2122 °C in Grignon (48°N) vs. 2830 °C in Australia (14°S). After silking, the values of the parameters ß and N0, and consequently the shape of the P2000(Na) curve, changed: 1·41 after silking vs. 2·45 before silking for ß and 0·43 vs. 0·27 for N0, while P2000,max changed slightly (40·2 vs. 36·8).
|
|
A similar behaviour was found between the pre- and post-silking photosynthetic measurements for quantum efficiency,
(Fig. 2C), on the one hand and for dark respiration Rd (Fig. 2D) on the other. No effect of Na was found on
, while Rd varied linearly with Na (Table 4). These results are in accordance with the literature (e.g. Schieving et al., 1992; Connor et al., 1993; Anten et al., 1995), except for findings by Hirose and Werger (1987) and Muchow and Sinclair (1994), who proposed a linear relationship between
and Na, and for findings by Pons et al. (1989) who proposed a curvilinear relationship between
and Na.
The model of lamina photosynthesis was globally assessed by comparing the simulated values of Pn,i,a (Pn,i,a,sim) from P2000(Na),
and Rd(Na) with the observed values (Pn,i,a,obs). We obtained Pn,i,a,sim = 1·001 ± 0·01 Pn,i,a,obs 0·251 ± 0·211 (P value = 0·0001, R2 = 0·97). Thus, the model with its parameterization is a suitable description of lamina photosynthesis.
Nm, Ma and Na partitioning within and between laminae
Results for Na have been described in a previous paper (Drouet and Bonhomme, 1999). The main trends for Na are recalled here in order to point out the links between Na, Nm and Ma.
Partitioning along the velum (Table 5).
At the beginning of stem elongation (60 DAS) under both density conditions (D- and d-density), Nm and Ma increased from the base to the tip of velums 6 and 10, especially between the base and the middle of the velum for Nm (significant at P < 0.05) and between the middle and the tip of the velum for Ma (P < 0.05). This produced a significant increase in Na from the base to the tip of velums 6 and 10. At the end of stem elongation (74 DAS), the same behaviour was observed for laminae 10 and 14 in both densities: Nm and Ma increased from the base to the tip of the velum (only between the middle and the tip of the velum for Ma), producing an increase in Na. At post-silking (90 DAS), slight variations in Nm and Ma in the opposite direction generated slight variations in Na along velum 6. This behaviour was also observed for velum 10, but only in the D-density. For velum 10 in d-, d10
2·5- and d10
0·5-densities, Nm and Ma also varied slightly in the opposite direction and no significant change was observed for Na. Nm remained constant along velum 14 in all densities, while Ma decreased significantly. A significant decrease in Na was also observed from the base to the tip of velum 14, except in the D-density.
Partitioning between laminae (Fig. 3).
Regarding changes in Nm, Ma and Na between laminae over time (D-density, Fig. 3, first line), significant differences were observed (P < 0.05) in Nm between laminae and between dates of measurements. Nm varied slightly around the mean values. Whatever the date, no strong gradient in Nm was observed between laminae. On the contrary, Ma in lower laminae, which had finished their growth (i.e. laminae 6 and 8), varied very slightly with no significant difference between mean values. Growing laminae had lower values of Ma than mature ones (e.g. for lamina 12, mean values of Ma increased from 33·2 g m2 at 60 DAS to 56·6 g m2 at 90 DAS). Strong gradients in Ma were observed between laminae (e.g. at 90 DAS, mean values of Ma varied from 40·7 g m2 for lamina 6 to 56·6 g m2 for lamina 12). Consequently, Na varied slightly over time in the laminae that had finished their growth (i.e. laminae 6 and 8) and increased in growing laminae (i.e. laminae 10, 12 and 14 between 6090 DAS). At post-silking (90 DAS), a gradient in Na was observed from the bottom to the top of the canopy. This may be the result of nitrogen remobilization from the lower laminae to the upper ones, nitrogen absorption and growth in Ma (and lamina thickness) of developing laminae. Regarding changes in Nm, Ma and Na between laminae as a function of plant density at post-silking (Fig. 3, second line), Nm increased significantly (P < 0.05) from the bottom to the top of the dense canopies (D- and d-density). Variations in Nm were not apparent in the open-thinned canopies (d10
2·5- and d10
0·5-density). Consequently, an increase in Nm from dense to open canopies was observed for lower laminae (i.e. laminae 6 and 8), while Nm remained relatively constant between dense and open canopies for upper laminae (i.e. laminae 1014). Similar vertical profiles of Ma were observed between densities, with higher values of Ma in the open canopies (d10
2·5- and d10
0·5-density) than in the dense ones (D- and d-density). Consequently, the profiles of Na partitioning were characterized by a stronger vertical gradient in the dense canopies (D-density) than in the open ones.
|
| DISCUSSION |
|---|
|
|
|---|
Contribution of a 3D approach to simulate canopy photosynthesis within heterogeneous row crops
Changes in Na along the velum (i.e. horizontal changes) were essentially observed until stem elongation (Table 5), whereas changes in Na between laminae (i.e. vertical changes) were observed after stem elongation (Fig. 3C). The effect of Na partitioning within canopies was quantified by comparing the daily net photosynthesis (Pn,d,c) between real and associated hypothetical canopies. Before stem elongation within the dense canopy (D-density), real Na partitioning (Table 5 and Fig. 3) made it possible to increase Pn,d,c from 1 to 7 % compared with uniform Na along each velum, and from 2 to 8 % compared with uniform Na within the whole canopy (Fig. 4A). After stem elongation, which corresponded to canopy closure, similar behaviour was observed but the increase in Pn,d,c was lower: from 2 to 4 % in the first case and from 4 to 6 % in the second. After silking, gradients along the velum were not significantly different in most canopies (Table 5) that involved no difference in Pn,d,c between the real D-density canopy and its associated hypothetical canopy with uniform Na along the velum. On the contrary, strong gradients in Na between laminae (Fig. 3C) involved increased Pn,d,c from 5 to 7 % compared with uniform Na within the whole canopy (Fig. 4A). This behaviour was also observed within the d-density canopy, but the difference between real and hypothetical canopies was lower: Pn,d,c increased by 2 % after silking (Fig. 4B). For the open-thinned canopies, Na gradients within and between laminae had disappeared, and the real Na partitioning gave similar Pn,d,c to uniform Na within laminae and within the whole canopy (Fig. 4B). The first hypothesis for Na partitioning (vertical Na gradients only) is thus valuable for the dense canopies after stem elongation. The second one, which does not take into account Na partitioning to simulate photosynthesis and biomass production, is valuable for the open-thinned canopies where no horizontal and vertical gradients were observed.
|
|
In heterogeneous canopies, vertical descriptions of Na partitioning can be used in dense ones to simulate canopy photosynthesis (e.g. Hirose and Werger, 1987; Anten at al., 1995; Dreccer et al., 2000), even if a vertical exponential Na partitioning underestimated Pn,d,c (Dreccer et al., 2000). But when lamina area density partitioning varies greatly along an axis perpendicular to the row direction (d-density, Fig. 1A, B), simulations indicated differences of canopy photosynthesis up to 8 % between real and associated hypothetical canopies. That involved taking into account spatial variations in Na (i.e. Na within and between laminae) by using 3D representations of the canopy variables. However, those differences are within the range of error of photosynthesis measurements in canopies.
In the open-thinned canopies (d10
2·5- and d10
0·5-density), total shoot nitrogen content increased within the plant compared with the d-density canopy: from 1·41 g per plant to 2·13 g per plant after the moderate thinning and from 1·41 g per plant to 2·30 g per plant after the strong thinning (Table 1). Thinning plants involved uniform Na partitioning, with Na tending towards a maximal value around 2 g m2 (Fig. 3F, see Drouet and Bonhomme, 1999), which might correspond to the maximal photosynthetic capacity of laminae. Simulations indicated that Na partitioning after removing plants involved an increase in Pn,d,c from 6 to 7 % within the open-thinned canopies compared with the Na partitioning within the d-density canopy. Because of foliage heterogeneity within the open-thinned canopies, 3D representations of the canopy variables were needed. But it was shown that, in our experimental open-thinned canopies, this approach was not finally necessary because of uniformity of Na and Id,a partitioning (see Drouet and Bonhomme, 1999).
Spatial changes in Na partitioning: changes in Nm or regulation by Ma?
Na is the product of Nm and Ma. Increased Ma is the result of increased lamina thickness (Maurice et al., 1997) and increased lamina density (Witkowski and Lamont, 1991). Relationships between Na and Nm and between Na and Ma for all canopies indicated that changes in Na jointly result from changes in Nm and Ma (Fig. 5). However, the relationship between Na and Ma was better (Table 6; RSE = 0·21 and R2 = 0·78) than that between Na and Nm (RSE = 0·28 and R2 = 0·61).
|
|
For each canopy, the relationship between Na and Nm was slightly stronger than that between Na and Ma at the beginning of stem elongation. At this stage, we can consider that Nm and Ma contributed equally to the variability of Na. Nm and Ma decreased simultaneously from the bottom to the top of the D-density (Fig. 3A, B) and d-density canopies (data not shown), which produced a decrease in Na. At the end of stem elongation, the tendency reversed: the relationship between Na and Nm was slightly weaker than that between Na and Ma. The vertical partitioning of Nm remained constant, whereas Na changed between laminae in the same way as Ma, which corresponded to vertical gradients of Id,a.
At post-silking, the relationship between Na and Ma was better than that between Na and Nm, except in the D-density canopy. Vertical Na distribution (Fig. 3C) corresponded with vertical Id,a distribution (see Drouet and Bonhomme, 1999), except at the top of the canopy where low Na values were observed within the most illuminated laminae (Fig. 1F). That might be due to Na remobilization from not only lower but also upper laminae towards the grains already developed at this period. This hypothesis agrees with results from several authors (e.g. Crafts-Brandner and Poneleit, 1987; Thomas and Smart, 1993; Sadras et al., 1993, 2000), who indicated an effect of reproductive growth on lamina senescence and lamina nitrogen partitioning. The latter in particular showed that the important changes in the profile of lamina nitrogen in maize and sunflower during grain filling were unrelated to the light regime. Vertical gradients were also observed in Ma (Fig. 3B), whereas vertical gradients in Nm were not apparent (Fig. 3A). In the open-thinned canopies, no correlation was found between Na and Nm (R2 = 0·09), whereas clear relationships were observed between Na and Ma (Table 6). These observations suggest that changes in Na might result from changes in Ma rather than changes in Nm, but this result was not as obvious as in previous studies (e.g. Hirose et al., 1989; Ellsworth and Reich, 1993; Rosati et al., 2000; Le Roux et al., 2001). Ma distribution therefore cannot provide a means of distributing Na within plant canopies, contrasting with results from Ellsworth and Reich (1993), for example, who found that the distribution of Ma in tree canopies explained 95 % of the distribution of Na in Acer saccharum. Furthermore, Hirose et al. (1988) showed that Ma in the canopy was controlled not only by light climate but also by nitrogen availability to each lamina. Light distribution therefore cannot provide a means of predicting Ma (and consequently Na) distribution, even if light plays a major role in Ma and Na distribution.
| CONCLUSIONS |
|---|
|
|
|---|
By using one cultivar of maize grown during a single season, it was shown that Na partitioning was closely related to the heterogeneity of irradiance within maize canopies until silking. Using a 3D approach to describe lamina nitrogen and irradiance, it was shown that the real patterns of Na partitioning could increase daily canopy photosynthesis by up to 8 % compared with uniform patterns of Na. That was especially observed for the earliest stages of development where lamina area density partitioning varies greatly along an axis perpendicular to the row direction. In row heterogeneous crops, the classical approach of vertical gradients for nitrogen and light within dense canopies therefore could not be applied and a 3D description of the canopy variables (irradiance, Nm, Ma, Na) would be required. Most studies have shown that Na is related to dry mass per unit area (Ma) rather than nitrogen concentration (Nm). Our data indicated that, even if Ma plays a major role in Na, we cannot evidence a strong relationship between Na and Ma. Na partitioning resulted mostly from Ma partitioning but also from Nm partitioning.
| ACKNOWLEDGEMENTS |
|---|
The authors are grateful to Dr V. O. Sadras and Dr N. P. R. Anten for valuable suggestions when reviewing the manuscript. We acknowledge Dr H. Sinoquet for valuable consultations concerning his RIRI model and, with P. Rivet, on the digitizing system, Polhemus. We thank Dr B. Andrieu for providing us with his surface-to-volume procedure. Thanks to F. Lafouge for the nitrogen analysis, as well as P. Bonchrétien, M. Lauransot, C. Civet, A. Fortineau and M. and J. Chartier for their technical assistance. S. Tanis-Plant gave us editorial advice.
| LITERATURE CITED |
|---|
|
|
|---|
-
Anten NPR, Schieving F, Werger MJA. 1995. Patterns of light and nitrogen distribution in relation to whole canopy carbon gain in C3 and C4 mono- and dicotyledonous species. Oecologia 101: 504513.[CrossRef][Web of Science]
Biscoe PV, Gallagher JN. 1978. A physiological analysis of cereal yield. I. Production of dry matter. Agricultural Progress 53: 3450.
Bonhomme R, Varlet-Grancher C. 1978. Estimation of the gramineous crop geometry by plant profiles including leaf width variations. Photosynthetica 12: 193196.
Bunce JA. 1990. Afternoon inhibition of photosynthesis in maize. I. Evidence, and relationship to stand density. Field Crops Research 24: 251260.[CrossRef][Web of Science]
Chartier P. 1966. Étude du microclimat lumineux dans la végétation. Annales Agronomiques 17: 571602.[Web of Science]
Chartier P. 1969. Assimilation nette dune culture couvrante. II. La réponse de lunité de surface de feuille. Annales de Physiologie végétale 11: 221263.
Chen JL, Reynolds JF, Harley PC, Tenhunen JD. 1993. Coordination theory of leaf nitrogen distribution in a canopy. Oecologia 93: 6369.[Web of Science]
Connor DJ, Hall AJ, Sadras VO. 1993. Effect of nitrogen content on the photosynthetic characteristics of sunflower leaves. Australian Journal of Plant Physiology 20: 251263.[Web of Science]
Connor DJ, Sadras VO, Hall AJ. 1995. Canopy nitrogen distribution and the photosynthetic performance of sunflower crops during grain filling a quantitative analysis. Oecologia 101: 274281.[CrossRef][Web of Science]
Crafts-Brandner SJ, Poneleit CG. 1987. Effect of ear removal on CO2 exchange and activities of ribulose bisphosphate carboxylase/oxygenase and phosphoenolpyruvate carboxylase of maize hybrids and inbred lines. Plant Physiology 84: 6169.
DeJong TM, Doyle JF. 1985. Seasonal relationships between leaf nitrogen content (photosynthetic capacity) and leaf canopy light exposure in peach (Prunus persica). Plant, Cell and Environment 8: 701706.
DeJong TM, Day KR, Johnson RS. 1989. Partitioning of leaf nitrogen with respect to within canopy light exposure and nitrogen availability in peach (Prunus persica). Trees 3: 8995.
Dreccer MF, van Oijen M, Schapendonk AHCM, Pot CS, Rabbinge R. 2000. Dynamics of vertical leaf nitrogen distribution in a vegetative wheat canopy. Impact on canopy photosynthesis. Annals of Botany 86: 821831.
Drouet J-L. 1998. Variations de la disposition spatiale et de la teneur en azote des feuilles dun peuplement de maïs en phase végétative. Étude par simulation de linfluence de ces variations sur la photosynthèse potentielle du couvert. PhD Thesis, Institut National Agronomique, Paris-Grignon, France.
Drouet J-L. 2003. MODICA and MODANCA: modelling the three-dimensional shoot structure of graminaceous crops from two methods of plant description. Field Crops Research 83: 215222.[CrossRef][Web of Science]
Drouet J-L, Bonhomme R. 1999. Do variations in local leaf irradiance explain changes to leaf nitrogen within row maize canopies? Annals of Botany 84: 6169.
Drouet J-L, Moulia B, Bonhomme R. 1999. Do changes in the azimuthal distribution of maize leaves over time affect canopy light absorption? Agronomie 19: 281294.[CrossRef][Web of Science]
Dumas JBA. 1831. Procédés de lanalyse organique. Annales de Chimie et de Physique 2: 198213.
Dwyer LM, Stewart DW. 1986. Effect of leaf age and position on net photosynthetic rates in maize (Zea mays L.). Agricultural and Forest Meteorology 37: 2946.
Ellsworth DS, Reich PB. 1993. Canopy structure and vertical patterns of photosynthesis and related leaf traits in a deciduous forest. Oecologia 96: 169178.[CrossRef][Web of Science]
Evans JR. 1989. Photosynthesis and nitrogen relationships in leaves of C3 plants. Oecologia 78: 919.[CrossRef][Web of Science]
Field C. 1983. Allocating leaf nitrogen for the maximization of carbon gain: leaf age as a control on the allocation program. Oecologia 56: 341347.[CrossRef][Web of Science]
Field C, Mooney HA. 1986. The photosynthesisnitrogen relationship in wild plants. In: Givnish TJ, ed. On the economy of plant form and function. Cambridge: Cambridge University Press, 2555.
Gosse G, Varlet-Grancher C, Bonhomme R, Chartier M, Allirand JM, Lemaire G. 1986. Maximum dry matter and solar radiation intercepted by a canopy. Agronomie 6: 4756.
Grassi G, Meir P, Cromer R, Tompkins D, Jarvis PG. 2002. Photosynthetic parameters in seedlings of Eucalyptus grandis as affected by rate of nitrogen supply. Plant, Cell and Environment 25: 16771688.[CrossRef]
Grindlay DJC. 1997. Towards an explanation of crop nitrogen demand based on the optimization of leaf nitrogen per unit leaf area. Journal of Agricultural Science 128: 377396.[CrossRef]
Gulmon SL, Chu CC. 1981. The effects of light and nitrogen on photosynthesis, leaf characteristics and dry matter allocation in the Chaparral shrub Diplacus aurantiacus. Oecologia 49: 207212.[CrossRef][Web of Science]
Gutschick VP, Wiegel FW. 1988. Optimizing the canopy photosynthetic rate by patterns of investment in specific leaf mass. American Naturalist 132: 6786.[CrossRef][Web of Science]
Hirose T, Werger MJA. 1987. Maximizing daily canopy photosynthesis with respect to the leaf nitrogen allocation pattern in the canopy. Oecologia 72: 520526.[CrossRef][Web of Science]
Hirose T, Werger MJA, Pons TL, van Rheenen JWA. 1988. Canopy structure and leaf nitrogen distribution in a stand of Lysimachia vulgaris L. as influenced by stand density. Oecologia 77: 145150.[CrossRef][Web of Science]
Hirose T, Werger MJA, van Rheenen JWA. 1989. Canopy development and leaf nitrogen distribution in a stand of Carex acutiformis. Ecology 70: 16101618.[CrossRef][Web of Science]
Hollinger DY. 1996. Optimality and nitrogen allocation in a tree canopy. Tree Physiology 16: 627634.
Kull O, Jarvis PG. 1995. The role of nitrogen in a simple scheme to scale up photosynthesis from leaf to canopy. Plant, Cell and Environment 18: 11741182.[CrossRef]
Kull O, Niinemets U. 1993. Variation in leaf morphometry and nitrogen concentration in Betula pendula Roth. Corylus avalana L. Lonicera xylosteum L. Tree Physiology 12: 311318.
Lemaire G, Gastal F. 1997. N uptake and distribution in plant canopies. In: Lemaire G, ed. Diagnosis of the nitrogen status in crops. Berlin: Springer-Verlag, 343.
Lemaire G, Onillon B, Gosse G, Chartier M, Allirand J-M. 1991. Nitrogen distribution within a lucerne canopy during regrowth: relation with light distribution. Annals of Botany 68: 483488.
Le Roux X, Walcroft AS, Daudet FA, Sinoquet H, Chaves MM, Rodrigues A, Osorio L. 2001. Photosynthetic light acclimation in peach leaves: importance of changes in mass:area ratio, nitrogen concentration, and leaf nitrogen partitioning. Tree Physiology 21: 377386.
Maurice I, Gastal F, Durand JL. 1997. Generation of form and associated mass deposition during leaf development in grasses: a kinematic approach for non-steady growth. Annals of Botany 80: 673683.
Meir P, Kruijt B, Broadmeadow M, Barbosa E, Kull O, Carswell F, Nobre A, Jarvis PG. 2002. Acclimation of photosynthetic capacity to irradiance in tree canopies in relation to leaf nitrogen concentration and leaf mass per unit area. Plant, Cell and Environment 25: 343357.
Mooney HA, Gulmon SL. 1979. Environmental and evolutionary constraints on the photosynthetic characteristics of higher plants. In: Solbrig OT, Jain S, Johnson GB, Raven PH, eds. Topics in plant population biology. New York: Columbia University Press, 316337.
Muchow RC, Sinclair TR. 1994. Nitrogen response of leaf photosynthesis and canopy radiation use efficiency in field-grown maize and sorghum. Crop Science 34: 721727.
Novoa R, Loomis RS. 1981. Nitrogen and plant production. Plant and Soil 58: 177204.[CrossRef][Web of Science]
Pachepsky LB, Haskett JD, Acock B. 1996. An adequate model of photosynthesis. I. Parametrization, validation and comparison of models. Agricultural Systems 50: 209225.
Polhemus. 1993. 3SPACE Fastrak Users Manual, Revision F. Colchester: Polhemus Vt, USA.
Pons TL, Schieving F, Hirose T, Werger MJA. 1989. Optimization of leaf nitrogen allocation for canopy photosynthesis in Lysimachia vulgaris (L.). In: Lambers H, Cambridge ML, Konings H, Pons TL, eds. Causes and consequences of variation in growth rate and productivity of higher plants. The Hague: SPB Academic Publishing, 175186.
Rosati A, Day KR, DeJong TM. 2000. Distribution of leaf mass per unit area and leaf nitrogen concentration determine partitioning of leaf nitrogen within tree canopies. Tree Physiology 20: 271276.
Sadras VO, Hall AJ, Connor DJ. 1993. Light-associated nitrogen distribution profile in flowering canopies of sunflower (Helianthus annuus L.) altered during grain growth. Oecologia 95: 488494.[CrossRef][Web of Science]
Sadras VO, Echarte L, Andrade FH. 2000. Profiles of leaf senescence during reproductive growth of sunflower and maize. Annals of Botany 85: 187195.
Sage RF, Pearcy RW. 1987. The nitrogen use efficiency of C3 and C4 plants. II. Leaf nitrogen effects on the gas exchange characteristics of Chenopodium album L. and Amaranthus retroflexus L. Plant Physiology 84: 959963.
Schieving F, Werger MJA, Hirose T. 1992. Canopy structure, nitrogen distribution and whole canopy photosynthetic carbon gain in growing and flowering stands of tall herbs. Vegetatio 102: 173181.[CrossRef][Web of Science]
Shiraiwa T, Sinclair TR. 1993. Distribution of nitrogen among leaves in soybean canopies. Crop Science 33: 804808.
Sinclair TR, Horie T. 1989. Leaf nitrogen, photosynthesis, and crop radiation use efficiency: a review. Crop Science 29: 9098.
Sinclair TR, Muchow RC. 1999. Radiation use efficiency. Advances in Agronomy 65: 215265.[CrossRef][Web of Science]
Sinoquet H, Bonhomme R. 1989. Modelling of solar radiation interception in row crops, II. Crop geometry and validation of the model. Agronomie 9: 619628.[CrossRef][Web of Science]
Sinoquet H, Bonhomme R. 1992. Modeling radiative transfer in mixed and row intercropping systems. Agricultural and Forest Meteorology 62: 219240.[CrossRef][Web of Science]
Sinoquet H, Rivet P. 1997. Measurement and visualization of the architecture of an adult tree based on a three-dimensional digitising device. Trees 11: 265270.[CrossRef]
Spitters CJT, Toussaint HAJM, Goudriaan J. 1986. Separating the diffuse and direct component of global radiation and its implication for modeling canopy photosynthesis, I. Components of incoming radiation. Agricultural and Forest Meteorology 38: 217229.[CrossRef][Web of Science]
S-PLUS. 1996. S-Plus users guide. Version 3.4. Seattle: MathSoft, WA.
Thornley JHM. 1976. Mathematical models in plant physiology. A quantitative approach to problems in crop and plant physiology. New York: Academic Press, USA.
Thomas H, Smart CM. 1993. Crops that stay green. Annals of Applied Biology 123: 193219.[CrossRef][Web of Science]
de Wit CT. 1965. Photosynthesis of Leaf Canopies. Agricultural Research Report No. 663. Center for Agricultural Publication and Document ation, Wageningen, The Netherlands.
Witkowski ETF, Lamont BB. 1991. Leaf specific mass confounds leaf density and thickness. Oecologia 88: 486493.[Web of Science]
This article has been cited by other articles:
![]() |
J. Bertheloot, P. Martre, and B. Andrieu Dynamics of Light and Nitrogen Distribution during Grain Filling within Wheat Canopy Plant Physiology, November 1, 2008; 148(3): 1707 - 1720. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||





