Skip Navigation


AOBPreview originally published online on August 4, 2004
Annals of Botany 2004 94(3):393-404; doi:10.1093/aob/mch155
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
94/3/393    most recent
mch155v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (3)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by SEGINER, I.
Right arrow Articles by BREUGELMANS, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by SEGINER, I.
Right arrow Articles by BREUGELMANS, M.
Agricola
Right arrow Articles by SEGINER, I.
Right arrow Articles by BREUGELMANS, M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?


Annals of Botany 94/3, © Annals of Botany Company 2004; all rights reserved

Modelling Ontogenetic Changes of Nitrogen and Water Content in Lettuce

IDO SEGINER1,*, PETER BLEYAERT2 and MAAIKE BREUGELMANS3

1 Faculty of Civil and Environmental Engineering, Division of Environmental Engineering, Water and Agriculture, Technion, Haifa 32000, Israel, 2 Provincial Research and Advisory Centre for Agriculture and Horticulture, B-8800 Rumbeke, Belgium and 3 Laboratory of Plant Ecology, Ghent University, B-9000, Ghent, Belgium

* For correspondence. E-mail segineri{at}tx.technion.ac.il

Received: 10 November 2003    Returned for revision: 23 March 2004    Accepted: 18 May 2004    Published electronically: 4 August 2004


   ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 

Background and Aims It is well established that the nitrogen content of plants, including lettuce, decreases with time. It has also been observed that water content of lettuce increases between planting and harvest. This paper is an attempt at modelling these observations.

Methods An existing dynamic model (NICOLET), designed to predict growth and nitrate content of glasshouse lettuce, is modified to accommodate the ontogenetic changes of reduced-nitrogen and water contents (on a dry matter basis). The decreasing reduced-N content and the increasing water content are mimicked by dividing the originally uniform plant into ‘metabolically active’ tissue and ‘support’ tissue. The ‘metabolic’ tissue is assumed to contain a higher nitrogen content and a lower water content than the ‘support’ tissue. As the plants grow, the ratio of ‘support’ to ‘metabolic’ tissue increases, resulting in an increased mean water content and a decreased reduced-N content. Simulations with the new model are compared with experimental glasshouse data over four seasons.

Key Results The empirical linear relationship between water and reduced-N contents, matches, to a good approximation, the corresponding relationship based on the model. The agreement between the two makes it possible to effectively uncouple the estimation of the ‘ontogenetic’ parameters from the estimation of the other parameters. The growth and nitrate simulation results match the data rather well and are hardly affected by the new refinement. The reduced-N and water contents are predicted much better with the new model.

Conclusion Prediction of nitrogen uptake for the substantial nitrate pool of lettuce depends on the water content. Hence, the modified model may assist in making better fertilization decisions and better estimates of nitrogen leaching.

Key words: Lactuca sativa L, lettuce, ontogenetic changes, nitrogen uptake, chemical composition, nitrogen content, water content, dynamic model, metabolic and support compartments


   INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
When plants are grown to maturity in agricultural stands, an ontogenetic decline of nitrogen content (so called ‘N-dilution’) is generally observed [Greenwood et al. (1990)Go – several C3 and C4 crops; Justes et al. (1994)Go – winter wheat; Sheehy et al. (1998)Go – rice; Colnenne et al. (1998)Go winter oilseed rape; Plénet and Lemaire (2000)Go maize; Greenwood and Draycott (1989)Go – various vegetable crops]. This phenomenon is generally modelled by viewing the plants as comprising two types of tissue. One, ‘metabolically active’, responsible for photosynthesis and nutrients uptake, with a high content of nitrogenous compounds; and the other, consisting of mechanical and vascular ‘support’ elements, with a low content of N-compounds (Caloin and Yu, 1984Go). Hirose and Werger (1987)Go and Chen et al. (1993)Go associated the metabolically active tissue with the sunlit leaves, and the support tissue with the shaded leaves and stems. As plants grow, the mass ratio of sunlit leaves to total plant material decreases. Nitrogen-rich compounds, no longer required in the shaded parts, are re-allocated to the younger leaves (Seligman, 1993Go), resulting in a decline of nitrogen content with depth in the canopy (Charles-Edwards et al., 1987Go; Grindlay et al. 1995Go; Alt et al., 2000Go) and in the canopy as a whole.

This description is inspired by the morphology of ordinary crops, where the young leaves are exposed to solar radiation more than the older leaves. In head-lettuce the young leaves are shaded and yet the same ‘dilution’ curve is evident. One way to treat this apparent inconsistency is to associate, as before, ‘metabolic’ with ‘sunlit’, despite the fact that the sunlit leaves are no longer the young ones. Another option is to assume that the composition of young and old leaves is, in evolutionary terms, a slow-changing trait, compared with leaf morphology (head formation). Finally, it is possible to abandon any attempt at an essential interpretation (for lack of suitable data) and view the model as an instrument, whose task is to produce correct predictions of crop nitrogen content. Our approach here is the last of these, retaining for convenience the description of an ‘ordinary’ crop to develop the lettuce model.

The total nitrogen content of plants is conveniently divided into nitrate (inorganic) and reduced (organic) nitrogen. The ratio nitrate-N : reduced-N has a wide range, from essentially zero (N-stressed plants; Broadley et al., 2003Go, Linker et al., 2004Go) to one or more (Table 1). The ‘critical’ nitrogen content, as traditionally plotted in N-dilution curves and believed to be required for maximum growth (Ulrich, 1952Go), seems to include all the reduced-N and possibly part of the nitrate-N. The extra nitrate, if any, is often considered to result from ‘luxury consumption’ (Justes et al., 1994Go; Grindlay, 1997Go). Whatever the definition of the components, excessive nitrate concentration is considered as posing a health hazard (European Commission, 1999Go).


View this table:
[in this window]
[in a new window]
 
TABLE 1. Summary of lettuce experiments carried out at Beitem during 2000

 


View larger version (37K):
[in this window]
[in a new window]
 
FIG. 4. Comparison between measured (asterisks) and simulated trajectories (lines) of the four data sets (columns), using variable compositional ratios. Seasons: 1, spring; 2, summer; 3, autumn; 4, winter. Cultivars: 1 and 2, Flandria; 3 and 4, Troubadour. DAP, days after (trans)planting. Five variables are presented (top to bottom): dry matter accumulation (linear), fresh matter accumulation (logarithmic), dry matter content, nitrate-N content and reduced-N content.

 


View larger version (38K):
[in this window]
[in a new window]
 
FIG. 6. Comparison between measured (asterisks) and simulated trajectories (lines) of the four data sets (columns), using fixed compositional ratios. Seasons: 1, spring; 2, summer; 3, autumn; 4, winter. Cultivars: 1 and 2, Flandria; 3 and 4, Troubadour. DAP, days after (trans)planting. Five variables are presented (top to bottom): dry matter accumulation (linear), fresh matter accumulation (logarithmic), dry matter content, nitrate-N content and reduced-N content.

 
Information regarding the ontogenetic changes of water content of plants is scant compared with information about nitrogen content. All annual plants end up dry, but during the vegetative phase of lettuce (and possibly of other crops), the water content seems to increase. Heinen et al. (1991)Go reported a significant increase, from 14·2 to 24·7 kg[water] kg–1[DM], and quoted Hansen (1976)Go as having observed a similar increase, from about 14·5 to about 25·3 kg[water] kg–1[DM]. With the ‘metabolic’/‘support’ model in mind, this implies that (during vegetative growth) the ‘support’ tissue contains more water than the ‘metabolic’ tissue, possibly due to relative abundance of vascular tissue. Results in support of the findings of Heinen and Hansen were also obtained in an experiment conducted recently in Beitem, Belgium (Bleyaert and Breugelmans; unpubl. res., 2001Go; next section). These new data clearly show that the ontogenetic decrease in reduced-N content is associated with a parallel increase in water content (Fig. 1).



View larger version (18K):
[in this window]
[in a new window]
 
FIG. 1. Relationship between water and reduced nitrogen, resulting from ontogenetic decrease of N content and increase of water content in soil-grown butterhead lettuce plants. Squares, spring; triangles, summer; diamonds, autumn; circles, winter. Spring and summer data for cultivar ‘Flandria’ and autumn and winter data for ‘Troubadour’. Data provided by P. Bleyaert and M. Breugelmans, Beitem, Belgium. The solid line is the regression through all the points. The dashed lines indicate the range between young (lower right) and mature plants (upper left).

 
Note, in passing, that we follow here the current convention (see list of references), where ‘content’ is preferred over ‘concentration’ to express mass of substance on dry mass basis. One exception is dry matter content, which is expressed on fresh mass basis.

Before moving on, it seems instructive to contrast the negative correlation of Fig. 1 with the positive correlation of Fig. 2, where water content is proportional to reduced-N content. This figure, based on recent laboratory data for relatively small isolated lettuce plants, shows that both reduced-N content and water content of N-stressed plants (bottom left) are considerably lower than those of normal plants (top right). In fact, the dry matter content (DMC) of highly N-stressed lettuce may become three to four times higher than the normal level of 4–6 % (on fresh-mass basis). These results are in agreement with earlier data, obtained during the growth of young seedlings, a stage termed ‘steady-state’ by Ingestad and Ågren (1988)Go, and ‘balanced exponential growth’ by Thornley (1997)Go. During this early stage of growth the composition of the seedlings is time-invariant and growth is exponential. If the plants are N-stressed, the N : C ratio decreases (Oscarson et al., 1989Go; Ingestad and Ågren, 1992Go) as does the water content (Oscarson et al., 1989Go), resulting in the pattern of Fig. 2.



View larger version (17K):
[in this window]
[in a new window]
 
FIG. 2. Decline of reduced-N content and of water content of hydroponically grown lettuce (various cultivars), resulting from nitrogen stress. Unstressed plants are represented by the top right data points (DMC = 4–6 %). Open symbols, data provided by F. Buwalda, Aalsmeer, The Netherlands; closed circles, Burns (1992)Go; triangles, de Pinheiro Henriques and Marcelis (2000)Go (after Seginer, 2003aGo).

 
The data in Figs 1 and 2, for the young, unstressed plants, should coincide. However, in Fig. 1 these are clustered around [2·5, 15], while in Fig. 2 they are clustered around [3, 20]. This discrepancy may be attributed to different lettuce cultivars, different growing conditions (Fig. 1 in soil, Fig. 2 in hydroponics) or, most likely, to different chemical analysis methods.

This paper is an attempt to model quantitatively the data of Fig. 1, by modifying a previously developed dynamic lettuce model, NICOLET (Seginer, 2003aGo), which is valid for the initial stage of growth, where composition is time-invariant (Fig. 2). The modification extends the validity of the model to a later stage of vegetative growth where ontogenetic composition-changes take place. We assume here, following Caloin and Yu (1984)Go, that viewing the plant as comprising distinct metabolic and support components is essentially correct, and that the trend shown in Fig. 1 is real.


   THE BEITEM DATA
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
During the years 1999 and 2000, experiments were conducted in the Provincial Research Centre at Beitem, Belgium (latitude 50°54'N, longitude 03°07'W, 30 m a.s.l.), in an attempt to quantify the effect of various environmental treatments (supplementing nitrogen, light, CO2 and heat) on the nitrate content of butterhead lettuce. The crops were soil-grown in glasshouse compartments, following common agricultural practices, except for the environmental treatments. The second-year experiments produced more consistent data than the first year and were, therefore, selected for the present modelling attempt. A detailed description of these experiments is available upon request (Bleyaert and Breugelmans, 2001Go; unpubl. res.).

Data from experiments where N supply was judged from the results to be ample and for which chemical analyses were available, were used for the current analysis. Since the interest here is not in the effects of the various treatments (which only rarely resulted in significant effects), the data were pooled to produce four distinct data sets, as shown in Table 1. The environmental data were averaged over the two experimental compartments (or four sub-compartments, when appropriate), and the crop data were averaged over the two higher N treatments, over four blocks (replications), as well as over the compartments (climatological treatments). Hence, each point in Fig. 1 and all the averaged fresh mass, dry mass, nitrate content and reduced-N content data (except the initial harvests), are means of at least eight individual determinations. The crops were sampled for analysis five or six times, at approximately uniform intervals of light integral, along the growing period, starting with the day of transplanting and ending at the final harvest.

Table 1 shows that the effects of season and cultivar are partially confounded. Cultivar ‘Flandria’ is associated with spring and summer and cultivar ‘Troubadour’ is associated with autumn and winter (no confounding for seasons within each cultivar). The ‘Flandria’ points in Fig. 1 appear to lie a bit lower than those of ‘Troubadour’. A few apparent irregularities can be found in Table 1 (bold font): (a) the transplants of the spring crop were unusually large; (b) the ‘light efficiency’ of the summer crop was considerably lower than that of the other crops; (c) the reduced-N of the autumn crop is unusually low.


   THE CURRENT NICOLET MODEL
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Description
The NICOLET model is described and justified in detail elsewhere (Seginer, 2003aGo). It consists of three (virtual) compartments, labelled ‘structure’, ‘vacuole’ and ‘excess-carbon’. The soluble organic and mineral constituents of the vacuole complement each other in contributing to the (constant) osmotic potential, while the carbon and organic-N compounds of the other two compartments are osmotically neutral. The structural N : C ratio and the water : structural-carbon ratio are assumed to be constant. The excess-carbon compartment is assumed to contain only ‘dry’ carbon compounds, and is usually activated only when the crop is N-stressed. Since the available Beitem data are for abundant supply of nitrogen, the excess-carbon compartment and the nutrition uptake formulation are omitted from the summary description provided here, as well as from the computations. A few simplifying notational changes have also been introduced. The resulting two-compartment model is a special case of the (original) three-compartment model. It may be obtained from the latter by setting {xi} = 0 in eqns [27] to [29] of Seginer (2003a)Go. The simplified model is shown in Fig. 3.



View larger version (13K):
[in this window]
[in a new window]
 
FIG. 3. A schematic description of the simplified NICOLET model. The two compartments and the carbon and nitrogen fluxes are shown. MC and MN are masses of carbon and nitrogen per unit ground area, where the lower-case subscripts, s for ‘structure’ and v for ‘vacuole’, indicate the compartments. V is the volume of water per unit ground area. FC and FN are fluxes of carbon and nitrogen. The fluxes with double compartment subscripts are between the compartments, and the other fluxes, namely FCp, FCg, FCm and FNu, are the photosynthesis, growth respiration, maintenance respiration and nitrogen uptake fluxes.

 
The fluxes are controlled by: (a) the environmental conditions – photosynthesis by light and CO2, and structural growth by temperature; (b) the ground cover of the crop; and (c) attenuation functions, which reflect the inhibition of certain processes. Photosynthesis is modelled as a Michaelis–Menten process, and respiration and growth are formulated as exponential functions of temperature (Q10 = 2). The ground cover, to which all fluxes are proportional, approaches one asymptotically, and the attenuation functions control the flows in times of stress.

If the environmental conditions are constant, the current NICOLET model predicts that the composition of the crop equilibrates with these conditions, namely remains constant with time. This will be changed below, to allow for ontogenetic effects.

System equations
In the following equations, the symbols denoting masses and fluxes are defined in Fig. 3. The environmental conditions I, T and CC are light flux, ambient temperature and CO2 concentration in the air. The parentheses {} enclose the arguments of functions, and t is time. All other symbols are constants, defined in the notation list and explained in Seginer (2003a)Go.

(i) Carbon balances (used as state equations)

(1)

(2)

(ii) Nitrogen balances

(3)

(4)

(iii) Compositional relationships

(5)

(6)

(7)
Equation (5) formulates the complementing role of nitrate, MNv, and soluble organic compounds, MCv, in maintaining the osmotic potential; eqn (6) formulates the proportionality between nitrogen, FNvs, and carbon, FCvs, in the structure; and eqn (7) formulates the proportionality between water, V, and structural carbon, MCs.

(iv) Normalized content of soluble carbon-compounds

(8)

(v) Carbon fluxes

(9)

(10)

(11)

(12)
Equations (9) to (11) are products of a function of the environment, p, e or g, a function of the size of the crop, f, and an attenuation function, h, which are given in more detail, as follows:

(13)

(14)

(15)

(16)
where f is the fraction of light intercepted by the canopy, and all fluxes are proportional to f.

(vi) Attenuation functions

(17)

(18)
The attenuation functions protect the vacuole from carbon over-spilling and over-draining.


   EXTENSION TO INCLUDE METABOLIC AND SUPPORT COMPONENTS
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Variable r and {lambda}
The originally uniform plant, with structural N : C ratio r and water : structural-C ratio {lambda} (eqns 6 and 7), is now divided into ‘metabolic’ and ‘support’ components. The structural N : C ratios of these components are denoted by r(m) and r(s), where m and s denote ‘metabolic’ and ‘support’, respectively, and r(m) is larger than r(s). The water : structural-C ratios are similarly denoted by {lambda}(m) and {lambda}(s), where {lambda}(m) is smaller than {lambda}(s).

By definition

(19)
The metabolic component is made proportional to the fraction of intercepted light, given by eqn (16). When the crop is young, namely when aMCs << 1,

(20)
and since at that time essentially all leaves are sunlit, the metabolic component may be defined by

(21)
With this, the total reduced-nitrogen content becomes

(22)
Note that is proportional to f{MCs}, which in view of eqns (9) to (12) makes the rate of photosynthesis, and all other metabolic activities, directly proportional to the reduced-nitrogen content of the metabolic compartment. This is in agreement with Ågren's (1985)Go view of ‘nitrogen productivity’.

Dividing eqn (22) by MCs, the N : C ratio of the structure as a whole is obtained:

(23)
Differentiating eqn (22) with respect to time and replacing state derivatives with the appropriate fluxes, the flux of nitrogen into the structure can be related to the flux of carbon as

(24)
which shows that initially the nitrogen flux is a fraction r(m) of the carbon flux, diminishing towards r(s) as the canopy closes (as f -> 1).

The total water content is, similar to eqn (22)

(25)
and the overall V : MCs ratio is, similar to eqn (23)

(26)
This extension to the original model adds two parameters [r(m), r(s), {lambda}(m) and {lambda}(s) instead of just r and {lambda}], but requires no additional state variables. The original, constant-r and constant-{lambda} formulation, is a special case of the new model and can be obtained by setting r(m) = r(s) = r and {lambda}(m) = {lambda}(s) = {lambda}.

Uptake of nitrogen
Differentiating the osmotica balance, eqn (5), with respect to time, and utilizing eqns (1), (2), (3), (6), (7), (23) and (26), results (Appendix A) in the flux form

(27)
where

(28)
Hence, the uptake of nitrogen (when N-supply is abundant) is

(29)


   CONVERSION BETWEEN MEASUREMENTS AND MODEL STATES
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
The comparison of model predictions with experimental data requires that the two can be mapped back and forth into each other. The experimental results of Bleyaert and Breugelmans can be expressed in terms of (a) fresh mass per unit ground area, WF, (b) dry mass per unit ground area, WD, (c) molar nitrate content on dry-mass basis, Cnit-N, and (d) molar reduced-N content on dry-mass basis, Cred-N.

From these, the water volume per unit ground area is determined via

(30)
where {rho} is the density of water. Combining this with eqn (26), an implicit expression for one of the state variables, MCs (eqn 2), is obtained

(31)

From the definition of Cnit-N,

(32)
and substituting eqns (31) and (32) into eqn (5), an expression for the second state variable, MCv (eqn (1)), is obtained

(33)
Equations (31) and (33) express the state variables in terms of the measured quantities WF, WD and Cnit-N as well as some of the model parameters.

The inverse conversion requires, in addition, an assumption regarding the conversion of dry matter, such as

(34)
where the conversion factors, {eta}C and {eta}N, are determined from the molecular masses of typical compounds.


   ESTIMATING THE PARAMETERS FOR THE ‘METABOLIC’ AND ‘SUPPORT’ COMPONENTS
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Analysis
Figure 1 shows a negatively sloping linear relationship between water content and reduced-N content. The empirically fitted line may be expressed as

(35)
where A (negative) and B are the slope and intercept.

The ontogenetic model is in general agreement with this observation, as can be shown by eliminating f{MCs}/aMCs between eqns (23) and (26), and multiplying throughout by {rho}MCs/WD, to obtain

(36)
The slope is strictly constant and negative and the intercept is positive and nearly constant, depending on the constancy of MCs/WD over time.

Equating eqns (35) and (36), results in

(37)
and

(38)
which places two constraints on the four parameters r(m), r(s), {lambda}(m) and {lambda}(s). Hence, given just the line of Fig. 1, two of the four parameters remain free for fitting. Further information can be extracted from the two ends of the data cluster of Fig. 1 if it is assumed that the lower right data points represent young plants with effectively no support tissue and the upper left data points represent mature plants consisting mainly of support tissue. For these data, utilizing eqns (23) and (26) and assuming the validity of eqn (20) for the young plants (subscript y) and MCs -> {infty} for the mature plants (subscript m), the result is

(39)

(40)

Since {rho}V/MNs is the slope of lines radiating from the origin in Fig. 1, denoting this slope for the young plants by y and that for the mature plants by m, eqns (37) to (40), can be used to evaluate the four ‘ontogenetic’ parameters as follows:

(41)

(42)

(43)

(44)

Evaluation for the Beitem data
The values of A, B, y and m can be read directly from Fig. 1, resulting approximately in




The value of m is an underestimate, since even at harvest the metabolic component is not negligible. Estimating that its appropriate value is about m = 70 kg[water] mol–1[N], using {rho} = 1000 kg[water] m–3, and substituting into eqns (41) to (44), the mean results over the two cultivars and over all seasons become




The common factor, WD/MCs, which depends to some extent on model parameters and environmental conditions, increases, as the simulated crop grows, from about 0·04 to about 0·05 kg[DM] mol–1[C], since the increasing water content implies more soluble compounds per unit of structural carbon. Multiplying the metabolic parameters (superscript m) by 0·04 and the support parameters (s) by 0·05 kg[DM] mol–1[C], yields the values of the four parameters as shown in Table 2, columns 1 and 2.


View this table:
[in this window]
[in a new window]
 
TABLE 2. Summary of parameter values used in the various simulations

 


View larger version (20K):
[in this window]
[in a new window]
 
FIG. 5. Comparison between measured (circles) and simulated trajectories (lines) of water content versus reduced-N content, using variable compositional ratios. Results for the four data sets are superimposed, trajectories starting at lower right. The measured points are the same as in Fig. 1.

 


View larger version (15K):
[in this window]
[in a new window]
 
FIG. 7. Comparison between measured (circles) and simulated trajectories (four overlapping lines) of water content versus reduced-N content, using fixed compositional ratios. Results for the four data sets are superimposed, data for young plants at lower right and for mature plants at upper left. The measured points are the same as in Fig. 1.

 


   SIMULATIONS
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Variable compositional ratios
The model was fitted to the available data, by changing some of the original NICOLET parameters, selected according to previous experience with the model (Broadley et al., 2003Go; Linker et al., 2004Go). It turns out that the two cultivars require somewhat different values for two of the parameters, {varepsilon} and {nu}, as shown in columns 1 and 2 of Table 2. Apparently, ‘Troubadour’ has a somewhat higher photosynthetic efficiency, {varepsilon}, and a somewhat lower partitioning to growth, {nu}, than ‘Flandria’; however, the differences are rather small and may result from the partial confounding between cultivar and season. The values of the compositional ratios r and {lambda} are the same for the two cultivars, as derived in the previous section from the common regression line of Fig. 1. All other parameters were the same for all simulations.

The results of the simulations are presented in Fig. 4, showing fair to good agreements between the measured and simulated values. The first two variables, dry and fresh mass, are cumulative and are predicted rather well. The last three variables are ratios, with noisier signals and looser fits. The logarithmic scale of the fresh mass and the linear scale of the dry mass together show how the initially exponential growth changes gradually into linear. The model was able to mimic the rather different shapes of the growth curves, resulting from different weather, cultivar and spacing. The prediction of nitrate content, the main focus of the original NICOLET model, is rather sensitive to short (daily) and medium term (weekly) environmental changes. In view of this, the fit to the various trends, including the nearly constant level of the winter experiment, should be considered as rather good.

The focus of this study is on the changes with time of the water and reduced-N content. The model predicts an increase with time of the former (decrease of DMC) and decrease of the latter. This is generally the trend of both data and simulation, but the details of the fit are less clear, partly because several data points are missing and partly because the scatter around the line in Fig. 1 is rather large. Figure 5 summarizes the simulations of water and reduced-N contents. Although the individual trajectories, starting at the lower right and ending at the upper left, do not closely agree with the data, the general behaviour in Fig. 5 is as expected.

Comparison with fixed compositional ratios
The question now is whether the new model produces better predictions than the original NICOLET model, where the water and reduced-N ratios were fixed (single r and {lambda}; original model). Selecting mean values for these ratios and adjusting slightly {varepsilon} and {nu} (Table 2), the results of Fig. 6 were obtained. The figure shows that the trajectories of dry mass, fresh mass and nitrate content, are generally similar to those of Fig. 4. The dry matter and reduced-N contents, however, are rather different. The difference between the models is particularly striking when viewed in terms of the co-ordinates of Figs 1 and 5, as shown in Fig. 7.

Comparison of Figs 5 and 7 shows clearly that the model with fixed compositional ratios cannot capture the changes occurring in water and reduced-N contents as the crop grows.


   DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Reduced-N and water contents
The new model, with ontogenetically decreasing reduced-N content and increasing water content, is a refinement of a previous model, NICOLET, which had been designed to predict growth and nitrate content of lettuce. The division of the ‘ontogenetic’ plant model, into ‘metabolic’ and ‘support’ tissue with high and low nitrogen contents, respectively, corresponds well with the morphology and nitrogen distribution of ‘ordinary’ crops. However, the extension of the concept to head lettuce and to water content can at the moment, while detailed data (of individual organs over the growing season) are not available, only be supported by the overall trends, such as shown in Fig. 1.

Figures 1 and 2 show the wide range of possible combinations of water and reduced-N content in lettuce. The negative correlation of Fig. 1 is due to ontogeny, while the positive correlation of Fig. 2 is due to N-stress. It may be of theoretical interest to explore what might be possible and impossible trajectories in the water versus reduced-N plane. The data of Fig. 1 are likely to mark the upper bound of the feasible region, while the data of Fig. 2 are likely to mark the lower bound.

Regarding the rate of ‘maturation’, in the sense of rate of N ‘dilution’ (moving from right to left in Fig. 1), this rate presumably depends on the rate of growth (faster growth leading to faster ‘dilution’). To some extent, however, it may also depend on plant spacing, isolated plants having higher reduced-N content than dense canopies of similar age (Seginer, 2003bGo).

Perhaps the most noteworthy point regarding the ontogenetic model is the similarity between eqns (35) and (36). The first is an empirical linear relationship between water and reduced-N contents, and the second is the model-formulated relationship, which is approximately linear. The correspondence between the two expressions makes it possible to effectively uncouple the estimation of the ‘ontogenetic’ parameters r(m), r(s), {lambda}(m) and {lambda}(s) (utilizing Fig. 1), from the estimation of the other parameters.

Uptake of nitrogen
The new model has the potential advantage of estimating nitrogen uptake more accurately than before. Partitioning the total uptake into structural reduced-nitrogen and vacuolar nitrate, NICOLET treats the former strictly as a function of the size of the plant (as measured by MCs) and the latter as a function of both size and environment. Regarding reduced-N, the new version predicts uptake rates for mature plants which are lower than those predicted by the original model, since in the new model the mature plants recycle much of their structural nitrogen (r(s) < r(n)). On the other hand, the increasing water content ({lambda}(s) > {lambda}(n)), results in an increased demand for nitrate to maintain the required osmotic potential. Judging by the potentially large nitrate content of lettuce, changes in the environment may have a significant effect on the uptake of nitrogen.

The accuracy of the uptake calculations depends on the accuracy of the predicted size and nitrogen content of the crop. At the moment, as shown by Fig. 8 (which may be regarded as a cumulative version of Fig. 1), there is some disagreement between model and data with respect to reduced-N accumulation in the crop. According to the model, all the points of Fig. 8 should lie on a single curve, obtained by plotting eqn (25) (times {rho}) against eqn (22). The data, however, show that two of the autumn points are considerably off the trend of the other three seasons (evident also in the reduced-N frames of Fig. 4). This may be blamed on the model (insufficient detail) and/or on the data (insufficient accuracy). An optimistic view would be to emphasize the good agreement of the data of three of the seasons, but definitely more data are required to establish a reliable set of the relevant parameters [r(m), r(s), {lambda}(m), {lambda}(s) and a].



View larger version (15K):
[in this window]
[in a new window]
 
FIG. 8. Cumulative mass of water against cumulative reduced-N, both on unit ground area basis. Squares, spring; triangles, summer; diamonds, autumn; circles, winter. Spring and summer data for cultivar Flandria and autumn and winter data for Troubadour. A cumulative version of Fig. 1.

 
Good estimation of the nitrate content is even more difficult, since it strongly depends on the environmental conditions. Considering the potentially high nitrate content of lettuce, the accuracy of its estimation should also be improved before the model can be used reliably for N-uptake calculations. For that, a good estimate of the ontogenetic change in water content is required.

Further points
The model was able to fit data of different seasons with the same set of (cultivar specific) parameters, even when the differences were as significant as between autumn and winter for the same cultivar (Troubadour; Fig. 4). The fit (Fig. 4) was generally good, despite significant differences in initial size of transplants (Flandria), in light efficiency (Flandria) and in final reduced-N content (Troubadour), as pointed out in conjunction with Table 1.

Comparison of Figs 4 and 6 shows that the new model does not improve significantly the prediction of growth and nitrate content of lettuce. It does, however, improve the prediction of water and reduced-N contents (Figs 5 and 7). Whether this improvement is useful to horticultural practice depends on the effect of these two attributes on lettuce quality and on the ability to make better nitrogen uptake predictions.

Most of the other tests of the NICOLET model have been against data from hydroponics (Broadley et al., 2003Go; Linker et al., 2004Go). The present data are of soil-grown lettuce which, however, was assumed to have a plentiful supply of nitrogen. The predictability for the soil-grown lettuce seems as good as for the hydroponics lettuce.

In conclusion, the new model is an improvement, as a comparison of Figs 5 and 7 demonstrates. It has the potential for improving the estimates of nitrogen uptake, thus leading to more accurate calculation of fertilizer needs and leaching potentials.


   APPENDIX A: DERIVATION OF EQN (27)
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
The time derivative of eqn (5) is

(A1)
From (26) and (16), respectively

(A2)

(A3)
Using eqns (16), (26), (A2) and (A3) to substitute into eqn (A1), yields

(A4)
Substituting from eqns (1), (2), (3) and (12) into eqn (A4)

(A5)
or


{mch155eq1}

(50)
or

(A7)
where

(A8)


   APPENDIX B: NOTATION
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
Main symbols

Symbol Description SI units

A slope of data line in Fig. 1kg[water] mol–1[N]
a light extinction coefficient m2[ground] mol–1[C]
B intercept of data line in Fig. 1 kg[water] kg–1[DM]
bk attenuation border of process k
CC CO2 concentration in the air mol[C] m–3[air]
c exponent in respiration equation K–1
e maintenance respiration rate of a closed-canopy crop mol[C] m–2[ground] s–1
FCk carbon flux due to process k mol[C] m–2[ground] s–1
FJij flux of constituent J from compartment i to compartment j mol[J] m–2[ground] s–1
FNu uptake rate of nitrogen mol[N] m–2[ground] s–1
f fraction of light intercepted by canopy
G a collection of parameters m3[water] Pa mol–1[C]
g uninhibited growth of closed canopy mol[C] m–2[ground] s–1
hk attenuation function of process k(hk = 1 when no attenuation)
I light flux mol[PAP] m–2[ground] s–1
k maintenance respiration rate at Ta = T* mol[C] m–2[ground] s–1
m slope of mature-plants line in Fig. 1 kg[water] mol–1[N]
MJj molar mass of constituent J in compartment j mol[J] m–2[ground]
p uninhibited gross photosynthesis of a closed canopy mol[C] m–2[ground] s–1
r N : C ratio in structure mol[N] mol–1[C]
sk attenuation slope of process k
T temperature K
T* reference temperature K
t time s
V volume of water m3[water] m–2[ground]
WD mass of dry matter of crop kg[DM] m–2[ground]
WF mass of fresh matter of crop kg[FM] m–2[ground]
y slope of young-plants line in Fig. 1 kg[water] mol–1[N]
ßJ osmotic potential associated with one unit of vacuolar J m3[water]Pa mol–1[J]
{Gamma} normalized carbon concentration in vacuole
{varepsilon} photosynthetic efficiency mol[C] mol–1[PAP]
{eta}J dry-mass equivalent of model components kg[DM] mol–1[J]
{theta} growth respiration as fraction of growth
{lambda} water associated with unit of structure m3[water] mol–1[C]
{nu} ratio g : e
{Pi} osmotic potential in vacuole Pa
{sigma} leaf conductance to CO2 m/s
{rho} density of water kg[water] m–3[water]

Subscripts
Compartments

i, j general compartment indices
s structure
v vacuole

Constituents

C carbon
J general constituent index
N nitrogen

Processes

g growth
k general process index
m maintenance
p photosynthesis
u uptake

Superscripts

(m) metabolically active component of plant
(s) support component of plant

Acronyms

DAP days after (trans)planting
DM dry matter
DMC dry matter content
FMfresh matter
PAP photosynthetically active photons

Notes

– indicates dimensionless quantities
{} used exclusively to contain the arguments of functions.


   ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 
This research has been funded by EU project FAIR6-CT98-4362 (NICOLET). We would like to acknowledge the constructive comments offered by Prof. Gerrit van Straten, Dr Fokke Buwalda and Dr Leo Marcelis, Wageningen University, The Netherlands.


   LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 THE BEITEM DATA
 THE CURRENT NICOLET MODEL
 EXTENSION TO INCLUDE METABOLIC...
 CONVERSION BETWEEN MEASUREMENTS...
 ESTIMATING THE PARAMETERS FOR...
 SIMULATIONS
 DISCUSSION
 APPENDIX A: DERIVATION OF...
 APPENDIX B: NOTATION
 ACKNOWLEDGEMENTS
 LITERATURE CITED
 

    Ågren GI. 1985. Theory for growth of plants derived from the nitrogen productivity concept. Physiologia Plantarum 64: 17–28.[CrossRef]

    Alt C, Kage H, Stützel H. 2000. Modelling nitrogen content and distribution in cauliflower (Brassica oleracea L. botrytis). Annals of Botany 86: 963–973.[Abstract/Free Full Text]

    Bleyaert P, Breugelmans M. 2001. Verification of basic model assumptions. In: Lemeur R, Hofman G, Calus A. FAIR6-CT98-4362 Technical report year 2. Ghent University, 65–69.

    Broadley MR, Seginer I, Burns A, Escobar-Gutiérrez AJ, Burns IG, White PJ. 2003. The nitrogen and nitrate economy of butterhead lettuce (Lactuca sativa var. capitata L.). Journal of Experimental Botany 54: 2081–2090.[Abstract/Free Full Text]

    Burns IG. 1992. Influence of plant nutrient concentration on growth rate: use of a nutrient interruption technique to determine critical concentrations of N, P and K in young plants. Plant and Soil 142: 221–233.[CrossRef]

    Caloin M, Yu O. 1984. Analysis of the time course of change in nitrogen content in Dactylis glomerata L. using a model of plant growth. Annals of Botany 54: 69–76.[Abstract/Free Full Text]

    Charles-Edwards DA, Stützel H, Ferraris R, Beech DF. 1987. An analysis of spatial variation in the nitrogen content of leaves from different horizons within a canopy. Annals of Botany 60: 421–426.[Abstract/Free Full Text]

    Chen JL, Reynolds JF, Harley PC, Tenhunen JD. 1993. Coordination theory of leaf distribution in a canopy. Oecologica 93: 63–69.

    Colnenne C, Meynard JM, Reau R, Justes E, Merrien A. 1998. Determination of a critical nitrogen dilution curve for winter oilseed rape. Annals of Botany 81: 311–317.[Abstract/Free Full Text]

    de Pinheiro Henriques AR, Marcelis LFM. 2000. Regulation of growth at steady-state nitrogen nutrition in lettuce (Lactuca sativa L.): interactive effects of nitrogen and irradiance. Annals of Botany 86: 1073–1080.[Abstract/Free Full Text]

    European Commission. 1999. Amending Regulation (EC) no. 194/97 setting maximum levels for nitrate in lettuce and spinach. EC Regulation no. 864/99.

    Greenwood DJ, Draycott A. 1989. Quantitative relationships for growth and N content of different vegetable crops grown with and without ample fertilizer-N on the same soil. Fertilizer Research 18: 175–188.

    Greenwood DJ, Lemaire G, Gosse G, Cruz P, Draycott A, Neeteson JJ. 1990. Decline in percentage N of C3 and C4 crops with increasing plant mass. Annals of Botany, 66: 425–436.[Abstract/Free Full Text]

    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, Cambridge 128: 377–396.[CrossRef]

    Grindlay DJC, Sylvester-Bradley R, Scott RK. 1995. The relationship between canopy green area and nitrogen in the shoot. INRA, Les Colloques No. 82, 53–60.

    Hansen H. 1976.. The content of nitrate and protein in lettuce (Lactuca sativa var. capitata (butterhead lettuce)) grown under different conditions. Tiddskrift Planteavl 80: 370–380 [In Danish with summary, figures and tables in English].

    Heinen M, de Jaeger A, Niers H. 1991. Uptake of nutrients by lettuce on NFT with controlled composition of the nutrient solution. Netherlands Journal of Agricultural Science 39: 197–212.

    Hirose T, Werger MJA. 1987. Nitrogen use efficiency in instantaneous and daily photosynthesis of leaves in the canopy of a Solidago altissima stand. Physiologia Plantarum 70: 215–222.[CrossRef]

    Ingestad T, Ågren GI. 1988. Nutrient uptake and allocation at steady-state nutrition. Physiologia Plantarum 72: 450–459.[CrossRef]

    Ingestad T, Ågren GI. 1992. Theories and methods on plant nutrition and growth. Physiologia Plantarum 84: 177–184.[CrossRef]

    Justes E, Mary B, Meynard JM, Machet JM, Thelier-Huche L. 1994. Determination of a critical nitrogen dilution curve for winter wheat crop. Annals of Botany 74: 397–407.[Abstract/Free Full Text]

    Linker R, Seginer I, Buwalda F. 2004. Calibration of a model for lettuce grown in a nitrate-limiting environment. Mathematical and Computer Modelling (in press).

    Oscarson P, Ingemarsson B, Larsson CM. 1989. Growth and nitrate uptake properties of plants grown at different relative rates of nitrogen supply. I. Growth of Pisum and Lemna in relation to nitrogen. Plant, Cell and Environment 12: 779–785.[CrossRef]

    Plénet D, Lamaire G. 2000. Relationships between dynamics of nitrogen uptake and dry matter accumulation in maize crops. Determination of critical N concentration. Plant and Soil 216: 65–82.[Web of Science]

    Seginer I. 2003a. A dynamic model for nitrogen-stressed lettuce. Annals of Botany 91: 623–635.[Abstract/Free Full Text]

    Seginer I. 2003b. Plant spacing effect on the nitrogen content of a crop. European Journal of Agronomy (in press).

    Seligman NG. 1993. Nitrogen redistribution in crop plants: regulation and significance. In: International Crop Science, Vol. I. Buxton DR et al., eds. Madison, WI: Crop Science Society of America, 813–829.

    Sheehy JE, Dionora MJA, Mitchell PL, Peng S, Cassman KG, Lemaire G, Williams RL. 1998. Critical nitrogen concentrations: implications for high-yielding rice (Oryza sativa L.) cultivars in the tropics. Field Crops Research 59: 31–41.[CrossRef]

    Thornley JHM. 1997. Modelling allocation with transport/conversion processes. Silva Fennica 31: 341–355.

    Ulrich A. 1952. Physiological bases for assessing the nutritional requirements of plants. Annual Review of Plant Physiology 3: 207–228.[Web of Science]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
Appl. Environ. Microbiol.Home page
M. T. Brandl and R. Amundson
Leaf Age as a Risk Factor in Contamination of Lettuce with Escherichia coli O157:H7 and Salmonella enterica
Appl. Envir. Microbiol., April 15, 2008; 74(8): 2298 - 2306.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
94/3/393    most recent
mch155v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (3)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by SEGINER, I.
Right arrow Articles by BREUGELMANS, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by SEGINER, I.
Right arrow Articles by BREUGELMANS, M.
Agricola
Right arrow Articles by SEGINER, I.
Right arrow Articles by BREUGELMANS, M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?