Agronomy Journal Journal of Natural Resources and Life Sciences Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online 12 March 2007
Published in Agron J 99:511-520 (2007)
DOI: 10.2134/agronj2006.0112
© 2007 American Society of Agronomy
677 S. Segoe Rd., Madison, WI 53711 USA
This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Weisz, R.
Right arrow Articles by Farrer, D. C.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Weisz, R.
Right arrow Articles by Farrer, D. C.
Agricola
Right arrow Articles by Weisz, R.
Right arrow Articles by Farrer, D. C.
Related Collections
Right arrow Other Crop Management
Right arrow Best Management Practices
Right arrow Wheat

Wheat

In-Season Tissue Testing to Optimize Soft Red Winter Wheat Nitrogen Fertilizer Rates

Influence of Wheat Biomass

Randy Weisza,*, Ravi P. Sripadaa, Ronnie W. Heinigera, Jeffrey G. Whiteb and Dianne C. Farrera

a Dep. of Crop Science, North Carolina State Univ., Raleigh, NC 27695
b Dep. of Soil Science, North Carolina State Univ., Raleigh, NC 27695

* Corresponding author (Randy_Weisz{at}ncsu.edu)

Received for publication April 12, 2006.

    ABSTRACT
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
In the southeastern USA, soft red winter wheat (Triticum aestivum L.) N fertilizer recommendations are based on growth stage (GS) 30 tissue testing and models that assume that the relationship between tissue N concentration (Ncon) and optimum N fertilizer rates (MaxN30) is stable across fields differing in GS-30 biomass. However, previous research has indicated this may not be the case. Consequently, it was critical to re-evaluate these models. Using a split-split plot design, six experiments were conducted in North Carolina between 2002 and 2004. Main plots were planting date–seeding rate combinations that produced wheat with different GS-30 biomass. Subplots and sub-subplots were five N rates applied at GS-25 and GS-30, respectively. Wheat yield was responsive to fertilizer N at all site-years. The overall relationship between MaxN30 and Ncon was weak (r2 = 0.43). The relationship between MaxN30 and N uptake (Ncon x biomass) was weaker (r2 = 0.27). However, when the data were divided into different biomass classes, the overall model improved (R2 = 0.75). For biomass < 340 kg ha–1, the Ncon at which no additional N fertilizer was required (Ncritical) was 70.0 g N kg–1. As biomass increased, Ncritical decreased to 33.2 g N kg–1. Intermediate classes had slopes of MaxN30 versus Ncon and Ncritical values that were similar to those previously reported. This study indicates that to use tissue testing to determine N fertilizer recommendations across a range of GS-30 biomass conditions requires information regarding dry matter biomass.

Abbreviations: GS, Zadoks winter wheat growth stage • LPOptN30, economically optimum growth stage-30 nitrogen fertilizer rates determined using linear-plateau functions • MaxN30, the growth stage-30 nitrogen fertilizer rate that resulted in the maximum yield • N25, nitrogen rate applied at growth stage 25 • N30, nitrogen rate applied at growth stage-30 • Ncon, whole plant tissue nitrogen concentration • Ncritical, whole plant tissue nitrogen concentration at which no additional nitrogen fertilizer was required to optimize yield • NUE, nitrogen-use efficiency • OptN30, economically optimum growth stage-30 nitrogen fertilizer rate • PDSR, planting date-seeding rate combination • QOptN30, economically optimum growth stage-30 nitrogen fertilizer rates determined using quadratic functions • Wbiomass, the wheat biomass at growth stage-30


    INTRODUCTION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
THE N-use efficiency (NUE) of winter wheat (Triticum aestivum L.) is generally low with, on average, only about 33% of applied N fertilizer utilized by the crop (Raun and Johnson, 1999). In-season tissue testing has been used to optimize wheat N fertilizer rates and increase NUE (Adams and Chapman, 1984; Alley et al., 1994; Weisz and Heiniger, 2004a). Tissue testing has been an especially valuable tool in warm humid regions of the world where soil mineral N content is too variable to be a good indicator of N availability to crops (Fox and Piekielek, 1978). For example, in the humid southeastern USA, soil mineral N measured in the spring when winter wheat is usually fertilized can range from 27 to 96% of the value measured in the fall before planting (Scharf and Alley, 1994).

In many winter wheat production areas, at least half of the total N fertilizer is applied pre-plant. An exception to this is the southeastern USA where most N fertilizer is applied at growth stage (GS) (Zadoks et al., 1974) 25 or 30. This is because warm temperatures and high precipitation typical of the southeastern USA result in denitrification and leaching of pre-plant N fertilizer (Counce et al., 1984; Scharf et al., 1993). Under these conditions, an N application may be recommended at GS-25 if tiller densities are low (Scharf and Alley, 1993; Weisz et al., 2001), but the majority of N fertilizer is usually applied at GS-30 based on a tissue test (Scharf et al., 1993; Flowers et al., 2004). In on-farm tests in Virginia, in-season N fertilizer rate optimization based on GS-25 tiller density and GS-30 tissue testing increased producer income by an average of $73 ha–1 (Scharf and Alley, 1993) and increased apparent fertilizer use efficiency (Scharf et al., 1993) compared with applying a single fixed N rate. In North Carolina, a similar in-season N fertilizer recommendation system (Weisz and Heiniger, 2004a) was compared with single fixed N-rate applications and was found to increase wheat yields up to 2.27 Mg ha–1 and to reduce N inputs by up to 48.6% (Flowers et al., 2004). The North Carolina and Virginia recommendation systems are based on a relationship between optimum GS-30 N fertilizer rates and GS-30 whole-plant tissue N concentration (Ncon) developed by Scharf et al. (1993).

Given the success in-season wheat N fertilizer rate recommendation systems have had, Flowers et al. (2004) tested the system developed by Weisz and Heiniger (2004a) on a within-field, site-specific basis. In these fields there was a large variance in Ncon. Using the field average Ncon and applying the recommended optimum N rate across the whole field resulted in over-application of N in five fields and under-application in two. Correcting for this by tissue sampling, determining optimum N rates, and applying GS-30 N on a site-specific basis improved fertilizer use efficiency but did not result in higher yields compared with the whole-field optimization. The results demonstrated an improvement in using tissue-based N fertilizer optimization compared with typical growers' practices, but the failure to improve yield on a site-specific basis suggested that the recommended fertilizer rates might not be accurate under conditions where crop growth varied (Flowers et al., 2004).

The question of how reliable tissue-based systems are for N fertilizer rate recommendations is not new. Early research in relating Ncon to N fertilizer requirements focused on finding "critical values" (Ncritical) or "sufficiency levels" for Ncon (i.e., the value of Ncon that is high enough that no further N fertilizer is required). Adams and Chapman (1984) reported an Ncritical of 27.5 g kg–1 in Arkansas. In North Dakota, Engel and Zubriski (1982) found Ncritical ranging from 52 to 56 g kg–1. Vaughan et al. (1990) reported Ncritical for dryland hard winter wheat grown in the west-central Great Plains (USA) of 32 g kg–1. Beringer and Hess (1979) concluded that Ncritical was too variable to be useful in Germany. Focusing specifically on soft red winter wheat in the eastern USA, Ncritical values of 38 to 41, 40, and 35 g N kg–1 were reported by Donohue and Brann (1984), Baethgen and Alley (1989), and Roth et al. (1989), respectively. The consistency of these results led to the conclusion that Ncon could be used to determine N fertilizer rate recommendations in Pennsylvania (Roth et al., 1989) and Virginia (Baethgen and Alley, 1989; Scharf et al., 1993). In the subsequent development of a relationship between optimum GS-30 N fertilizer rates and Ncon, Baethgen and Alley (1989) and Scharf et al. (1993) limited their sampling to experimental plots or to farm fields that were uniform and had high yield potential. Even so, the relationships they found between optimum N rate and Ncon were weak (r2 = 0.59 and 0.51).

We hypothesized that the relationship between soft red winter wheat Ncon and GS-30 Ncritical or optimum N fertilizer rates may not be as stable across the wide range of conditions typically found in wheat fields. If this is true, then the model developed by Scharf et al. (1993) may not be well adapted to low-density wheat with low yield potential, and this might explain the finding of Flowers et al. (2004). We wanted to determine if the relationship between optimum GS-30 N fertilizer rate and Ncon might be more sensitive to the density of the wheat crop (e.g., its biomass) than previously thought. If so, this might explain some of the variation in the model developed by Scharf et al. (1993). This is important because this recommendation system is being used by southeastern USA producers to determine N fertilizer rates across a wide range of field conditions, not only in high-yield potential fields. In addition, such information would be valuable for site-specific N management.


    MATERIALS AND METHODS
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Research Locations and Experimental Design
Experiments were conducted at four sites throughout North Carolina during 2002, 2003, and 2004. Specific locations were the Piedmont Research Station near Salisbury, NC, in 2003 (P2003); the Cunningham Research Station near Kinston, NC, in 2002, 2003, and 2004 (C2002, C2003, and C2004); the Lower Coastal Plain Tobacco Research Station near Kinston, NC, in 2002 (L2002); and the Tidewater Research Station near Plymouth, NC, in 2004 (T2004). Taxonomic classification of the soils at these sites is shown in Table 1.


View this table:
[in this window]
[in a new window]

 
Table 1. Site-year code (Site ID), location, soil series, and soil taxomomic classification for each experimental location.

 
At each site-year, an RCBD in a split-split plot arrangement with five replications was used. Main plots consisted of three planting date–seeding rate (PDSR) combinations designed to produce wheat stands with different amounts of GS-30 biomass (Table 2). The PDSR combinations were designed so that PDSR-1 would produce the thickest wheat stand with the highest biomass, whereas PDSR-2 and PDSR-3 would produce progressively thinner stands with lower biomass. Subplots consisted of five N rates applied at GS-25 (N25). These N25 rates were 0, 22, 45, 67, and 101 kg N ha–1 at P2003 where high N carry-over was anticipated due to very low grain yields in the previous corn (Zea mays L.) crop and 0, 34, 67, 101, and 134 kg N ha–1 at all other locations. Sub-subplots consisted of five N rates (0, 34, 67, 101, and 134 kg N ha–1) applied at GS-30 (N30).


View this table:
[in this window]
[in a new window]

 
Table 2. Site-year code (Site ID), year, cultivar, planting date–seeding rate combination (PDSR), planting date, seeding rate, sub-subplot size, and wheat row spacing used in a split-split plot design at each location.

 
Agronomics
Planting dates, seeding rates, row spacing, sub-subplot sizes, and cultivar used at each site-year are described in Table 2. At C2002, C2003, C2004, and L2002, PDSR-1 consisted of a high seeding rate (480 seeds m–2) and timely planting date (Weisz and Heiniger, 2004b), and PDSR-2 and PDSR-3 were late plantings at high (480 seeds m–2) and low (185 seeds m–2) seeding rates, respectively. Weather conditions prevented timely planting at P2003 and T2004, so all three PDSR treatments were planted on the same date and consisted of three different seeding rates (Table 2). Site-years C2002, C2003, and C2004 received a pre-plant N application of 34 kg N ha–1 as N–P–K: 10–13.2–24.9% (N-P2O5–K2O: 10–30–30%) (N source unknown). Pre-plant N was not applied at other site-years. Growth stage 25 and 30 N treatments were applied as aqueous urea-ammonium nitrate (30% N) except at P2003, where prilled ammonium nitrate (34% N) was applied at GS-25.

Lime and fertilizer rates other than N followed standard recommendations for North Carolina based on annual soil tests (Hardy et al., 2002; Crozier et al., 2004). All study sites followed corn and were conventionally tilled. Pre- and postemergence herbicides were applied as recommended (York, 2004), and weed management was excellent at all site-years except L2002, where weed cover between GS-25 and GS-30 was rated at approximately 10% in PDSR-1, 14% in PDSR-2, and 33% in PDSR-3.

Data Collection
The number of tillers with a minimum of three leaves in a 1-m section of row was determined at two random locations in each subplot before N25 application. This resulted in 10 samples per main plot, and main plot tiller density was estimated as the average of these samples. Plant samples for biomass and Ncon were taken at GS-30 before N30 applications. These samples were taken randomly within each sub-subplot from four 0.46-m sections of row by cutting whole plants just above the soil level. Subplot biomass was estimated as the mean of these 20 samples (i.e., four samples within each sub-subplot multiplied by five sub-subplots per subplot). At site-year L2002, where weed populations were high, weed tissue was removed from the samples before biomass and Ncon measurements were made. Plant samples were dried at 60°C for 48 h for dry matter determination. Tissue N concentration of the dried samples was determined by Waters Agriculture Laboratories (Camilla, GA) using a CHN analyzer (McGeehan and Naylor, 1988). At harvest, the center 2 m of each sub-subplot was cut with a Massey Ferguson MF-8 or Gleaner K2 plot combine (AGCO Corp., Duluth, GA), and yields were measured with a HarvestMaster grain gauge (Juniper Systems, Logan, UT). Yields were adjusted to a moisture content of 135 g kg–1.

Statistical Analysis
The PDSR combinations were designed to produce different GS-30 biomass levels at each site-year and depended on the soils and actual dates when wheat could be planted at each site-year. Consequently, PDSR combinations were not uniform across site-years, and ANOVA for GS-25 tiller density, GS-30 biomass, and grain yield were done individually by site-year. At GS-25, ANOVA for a RCBD was used to determine the effects of PDSR on tiller density. At GS-30, ANOVA for a split-plot design was used to determine the effects of PDSR and N25 on wheat biomass. After harvest, ANOVA for a split-split plot design was used to determine the effects of PDSR, N25, and N30 on grain yield at each site-year. For these analyses, PROC GLM in SAS Version 6 (SAS Institute, Cary, NC) was used.

Determination of Maximum Growth Stage-30 Nitrogen Rates
For each PDSR-N25 combination at each site-year, the N30 fertilizer rate that resulted in the maximum yield (MaxN30) was calculated. Initially, the yield responses to N30 were modeled using linear-plateau, quadratic-plateau, the exponential Mitscherlich, and second-order polynomial functions. In the majority of cases the quadratic-plateau, Mitscherlich, and second-order polynomial functions failed to fit the data (i.e., they were unable to predict a statistically significant MaxN30 or they predicted a value that was higher than what is agronomically reasonable for soft red winter wheat). The yield data were, however, well modeled with the linear-plateau function (for example see Fig. 1 ), and consequently MaxN30 values were estimated with the following rules:

  1. The yield response to N30 was determined by ANOVA in PROC GLM. If the response to N30 was not significant, MaxN30 was assumed to be zero.
  2. For PDSR-N25 combinations for which the response was significant, Fishers Protected LSD ({alpha} = 0.05) was used for means separation. If the highest mean yield was for N30 = 0, then MaxN30 was determined to be zero.
  3. If data sets had significant and positive yield responses to N30, they were analyzed as follows:
    a. A linear-plateau function was fit to the data using PROC NLIN (SAS Institute, Cary, NC). If the 95% confidence interval for the estimated slope of the linear portion of the model, and for the breakpoint (the N30 value beyond which the yield response was flat) did not contain zero, MaxN30 was determined to be equal to the breakpoint (Fig. 1). In four cases, the yield response to N30 was initially positive, reached a plateau, and then decreased at higher N rates. When this occurred, the N30 treatment levels that had significantly decreased yield compared with the N30 rate with the highest yield as determined by the Fisher's protected LSD ({alpha} = 0.05) mean separation were deleted.
    b. If the linear-plateau function did not result in a meaningful fit to the data (i.e., if the 95% confidence interval for the linear slope or for the breakpoint included zero), then treatment means were compared in PROC GLM using a series of linear contrasts to determine MaxN30.


Figure 1
View larger version (15K):
[in this window]
[in a new window]

 
Fig. 1. Response of wheat yield (mg ha–1, solid squares) and partial profit ($ ha–1, open circles) to N fertilizer applied at growth stage (GS) 30. Yield was modeled with a linear-plateau function. The N rate resulting in the highest yield (MaxN30) was estimated as the breakpoint. Partial profit was computed assuming a N/grain price ratio of 10, modeled with a quadratic function, and economic optimal N rate (OptN30) estimated as the inflection point. These example data were from the Lower Coastal Plains Research Station 2002, planting date–seeding rate combination = 3, and GS-25 N application = 67 kg ha–1.

 
Maximum Nitrogen Rates, Tissue Nitrogen Concentration, Biomass, and Critical Nitrogen Values
Once MaxN30 values were determined, linear regression was used to evaluate the overall relationship between these rates and Ncon. Forward and backward stepwise regression in PROC GLM was used to explore potential relationships of MaxN30 (dependent variable) with Ncon and biomass (independent variables). Independent variables in this analysis included Ncon, biomass, Ncon x biomass (usually defined as total N uptake [NUP]), Ncon2, and biomass2. Analysis using biomass as a class variable and Ncon as a continuous variable in PROC GLM was used to further elucidate the relationships between these agronomic variables. Values of Ncritical were estimated as the x axis intercept of significant regressions of MaxN30 versus Ncon.

Determination of Optimum Economic Growth Stage-30 Nitrogen Rates
The economically optimum N30 rate (OptN30) is not necessarily equal to MaxN30. Baethgen and Alley (1989) showed that because wheat yield response to N is sufficiently large, the impact of changing grain and/or N prices on OptN30 is small. In fact, they showed that for N/grain price ratios between 2 and 6, OptN30 was essentially equal to MaxN30. However, since 1989, the price of N has increased significantly, and price ratios between 6 and 10 are more representative of current conditions. Consequently, we evaluated the impact a price ratio of 10 would have on OptN30 calculations.

When the yield response to N30 is modeled with a linear-plateau function, MaxN30 (estimated at the breakpoint) is equal to OptN30 if the linear slope is greater than or equal to the N/grain price ratio. If the slope is lower than the N/grain price ratio, then the linear-plateau model predicts that OptN30 will be equal to 0. This approach was used for data sets that conformed to rule 3a. For data sets that followed rules 1 or 2, OptN30 was determined to be zero. Data sets that followed rule 3b (n = 13) were eliminated from this analysis.

A second and potentially more robust way to estimate OptN30 is to convert the yield data directly to partial profit and model the partial profit response to N30. For each PDSR-N25 combination at each site-year, partial profit was estimated by subtracting the cost of N30 from the income received from the grain yield assuming a N/grain price ratio of 10. These data were well fit with a second-order polynomial function (see Fig. 1). If the inflection point was below zero, OptN30 was set to zero; otherwise, OptN30 was determined to equal the N30 value at the inflection point.


    RESULTS
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Tiller Density, Biomass, and Grain Yield
Planting date–seeding rate combinations successfully produced different tiller densities at all site-years (Table 3), with mean GS-25 tiller density being consistently highest for PDSR-1 and lowest for PDSR-3. Biomass levels at GS-30 could have been affected by PDSR and N25. There was a significant PDSR x N25 interaction only at C2003. The main effects of PDSR and N25 were significant at all site-years except C2004, where only N25 had an effect (Table 3). At all site-years except C2004, GS-30 biomass was higher with earlier planting and/or higher seeding rates (Fig. 2A ). At all site-years, GS-30 biomass increased as N25 increased to 67 kg ha–1, with no further increase at higher N rates (Fig. 2B). Despite the PDSR x N25 interaction at C2003, these trends were still apparent at that site-year (data not shown). These data indicate that the PDSR treatments were successful in establishing a wide range of GS-30 biomass.


View this table:
[in this window]
[in a new window]

 
Table 3. Growth stage (GS)-25 wheat tiller density, GS-30 biomass, and grain yield ANOVA results for main effects and their interactions at six site-years (C2002, C2003, C2004, L2002, T2004, and P2003) for the treatment factors planting date-seeding rate combination (PDSR), N applied at GS-25 (N25), and N applied at GS-30 (N30). For the grain yield ANOVA, site-year means and model CV and R2 values are also shown.

 

Figure 2
View larger version (28K):
[in this window]
[in a new window]

 
Fig. 2. Mean growth stage (GS) 30 biomass for (A) wheat planted with planting date–seeding rate combinations (PDSR) 1, 2, and 3. (B) Plots receiving five different rates of N at GS-25. These data are from five site-years (Cunningham Research Station 2002, and 2004 indicated as C2002 and C2004, respectively; Lower Coastal Plains Research Station 2002, L2002; Piedmont Research Station 2003, P2003; and Tidewater Research Station 2004, T2004) where the PDSR by GS-25 N rate interaction was not statistically significant. Means within site-years with different letters were significantly different at the 0.05 probability level.

 
Mean grain yield at each site-year ranged from 3.0 to 6.1 Mg ha–1. These extremes were recorded at P2003 and T2004, respectively (Table 3). Site-year P2003 was planted very late (Table 2), had an unusually cold winter, and had a moderate infestation of leaf and glume blotch [Stagonospora nodorum (Berk.)]. These factors could have contributed to the low mean yield at that site-year. Grain fill at T2004 was longer than usual due to cool temperatures and high rainfall. That combined with disease-free stands may have been responsible for the high mean yield at T2004. Planting date–seeding rate combinations were significant factors explaining grain yield at five of the six site-years as indicated by a significant main effect or interaction term involving PDSR (Table 3). The main effects and interaction of N25 and N30 were significant at all site-years, demonstrating that wheat yield was responsive to N fertilizer.

Growth Stage-30 Nitrogen Rates for Maximum Yield
At each site-year, grain yield response to N30 for each PDSR–N25 combination was analyzed using ANOVA (data not shown). For 23 PDSR–N25 combinations, N30 was not significant, or the highest yield was attained at N30 = 0 kg ha–1. In these cases, MaxN30 was defined as 0 kg N ha–1. For the remaining PDSR–N25 combinations, response to N30 was modeled using a linear-plateau function (e.g., Fig. 1). Of these, 52 resulted in meaningful fits to the data, and MaxN30 was estimated as the resultant model breakpoint. For 13 PDSR–N25 combinations, MaxN30 was defined using a series of linear contrasts.

The overall relationship between MaxN30 and Ncon (MaxN30 = 156.5 – 2.73Ncon) was weak (r2 = 0.43) (Fig. 3A ). The linear relationships between MaxN30 and Ncon previously reported by Scharf et al. (1993) (MaxN30 = 235 – 4.8 Ncon; n = 136; r2 = 0.51) and by Baethgen and Alley (1989) (MaxN30 = 314 – 7.90 Ncon; n = 16; r2 = 0.59) are shown in Fig. 3A. The relationship reported by Baethgen and Alley (1989) had the steepest slope and lowest Ncritical (x axis intercept) of 39.7 g N kg–1. Scharf et al. (1993) used substantially more data in their analysis, which resulted in a somewhat lower slope and Ncritical close to 49.0 g N kg–1. The slope we found for the overall relationship between MaxN30 and Ncon was substantially lower than those previously reported and resulted in an Ncritical level of 57.3 g N kg–1.


Figure 3
View larger version (27K):
[in this window]
[in a new window]

 
Fig. 3. The growth stage (GS)-30 N fertilizer rate that resulted in the maximum yield (MaxN30) versus GS-30 whole plant tissue N concentration (Ncon) for the entire data set (A, open circles) and for five GS-30 biomass classes (B–F). The entire data set is shown in each graph, but for B through F, individual biomass classes are highlighted using larger solid symbols, and the linear regression model for that biomass class only (solid line and r2) is indicated. The linear models previously reported by Scharf et al. (1993) (dotted line) and by Baethgen and Alley (1989) (dashed line) are shown for comparison in each graph.

 
Our primary objective was to determine if the relationship between MaxN30 and Ncon varied depending on GS-30 biomass. Regression analysis showed that there was no statistically significant relationship between MaxN30 and GS-30 biomass or biomass2. There was a significant relationship between MaxN30 and NUP (the product of Ncon and GS-30 biomass), but the relationship was weak (r2 = 0.27). Allowing forward and backward stepwise regression of MaxN30 versus GS-30 biomass, Ncon, NUP, biomass2, and Ncon2 to find the combination of independent variables with the best model fit resulted in the following relationships, respectively:

Formula 1[1]

Formula 2[2]
where Wbiomass (kg ha–1) was the GS-30 wheat biomass. Both models had an R2 = 0.69. Consequently, biomass was an important parameter in determining MaxN30. To visualize this biomass effect on the relationship between MaxN30 and Ncon, we divided the complete data set into five biomass classes. The relationships between MaxN30 and Ncon for each of these five classes are illustrated graphically in Fig. 3BF.

Biomass Less than 340 kg ha–1
Maximum N rates for wheat with biomass values lower than 340 kg ha–1 fell along the upper right edge of the total data cloud (Fig. 3B). The linear regression model for this biomass class was only weakly significant (r2 = 0.49; p = 0.08). Although the data seemed to be more quadratic than linear, the quadratic term in a second-order polynomial fit to these data was not statistically significant. Maximum N30 rates for this biomass class were generally higher than those reported by Scharf et al. (1993) (dotted line in Fig. 3B) and had a much higher Ncritical (x axis intercept) of 70.0 g N kg–1 (Table 4). To our knowledge, this is higher than any Ncritical value for soft red winter wheat reported in the literature.


View this table:
[in this window]
[in a new window]

 
Table 4. ANOVA results for modeling maximum growth stage (GS) 30 N fertilizer rates (MaxN30), economically optimum GS-30 N rates based on quadratic models of partial profit (QOptN30), and economically optimum GS-30 N rates based on linear-plateau models of yield (LPOptN30) as a linear function of tissue N concentration (Ncon) and biomass, with biomass as a class variable and Ncon as a continuous variable. The biomass classes, the parameter estimates for each biomass class, the SEs of those estimates, and the resultant N sufficiency level are shown.

 
Biomass from 340 to 1000 kg ha–1
As GS-30 biomass increased to the range from 340 to 1000 kg ha–1, the data subset shifted to the left (Fig. 3C). The linear regression of MaxN30 versus Ncon (r2 = 0.75; p < 0.001) indicated a lower Ncritical (x axis intercept) of 55.1 g N kg–1 (Table 4) and had a slope that was similar to that reported by Scharf et al. (1993) (dotted line in Fig. 3C).

Biomass from 1000 to 1500 kg ha–1
The next highest biomass class (1000–1500 kg ha–1) fell in the center of the full data cloud (Fig. 3D). The linear regression model of this data subset had a Ncritical (x axis intercept) of 49.2 g N kg–1 (Table 4). This was nearly equal to the Ncritical found by Scharf et al. (1993) (dotted line in Fig. 3D). Based on the overlap of the SEs of the slopes of this regression line and either of those from the two previous biomass classes (Table 4), the slopes of these three regressions (Fig. 3BD) were not statistically different. Apparently, as the GS-30 biomass increased from below 340 up to 1500 kg ha–1 the major effect was to shift the relationship between MaxN30 and Ncon progressively to the left resulting in decreasing Ncritical values.

Biomass from 1500 to 2000 kg ha–1
This biomass class fell in the left half of the full data cloud (Fig. 3E) and had a steeper slope compared with the lower biomass classes based on the non-overlapping SEs of the slope estimates (Table 4). This biomass class had a linear regression model that was similar to that previously reported by Baethgen and Alley (1989) (dashed line in Fig. 3E). The value of Ncritical for this biomass class (41.4 g N kg–1) continued the trend of lower values with increasing GS-30 biomass.

Biomass Greater than 2000 kg ha–1
The highest biomass class (>2000 kg ha–1 (Fig. 3F) was on the left-most edge of the data cloud. The linear regression model for this biomass class had a slope similar to that reported by Baethgen and Alley (1989) (dashed line in Fig. 3F) and the lowest Ncritical (33.2 g N kg–1) of any biomass class.

Combined Model of Maximum Growth Stage-30 Nitrogen Fertilizer Rates
The stepwise regression models in Eq. [1] and [2] imply that as biomass increases, the slope of the linear relationship between MaxN30 and Ncon continually decreases (becomes more negative). Figures 3BF demonstrate that this was not the case. The slope was relatively stable across the three lower biomass classes (Fig. 3BD) and decreased in the two higher biomass classes (Fig. 3E and 3F and Table 4). To more accurately represent these data, MaxN30 values from the entire data set were modeled using Ncon as a continuous variable and biomass as a class variable (Table 4). Tissue N concentration and the interaction of Ncon and biomass class were significant (p = 0.0001), with an overall model R2 = 0.75. The full model fit for each biomass class is represented by the regression lines shown in Fig. 3BF.

Economically Optimum Growth Stage-30 Nitrogen Rates
Because the relationship between MaxN30 and Ncon was not statistically significant (at the 0.05 probability level) for wheat with biomass <340 kg ha–1, this biomass class was not included in the economic analysis. Following the results of the statistical analysis of MaxN30, the economically optimum N30 rates that were determined using quadratic functions fitting partial profit (QOptN30) and from linear-plateau functions fitting yield (LPOptN30) were modeled using Ncon as a continuous variable and biomass as a class variable (Table 4). For QOptN30 and LPOptN30, tissue N concentration and the interaction of Ncon and biomass class were significant (p = 0.0001), and the overall model R2 values were 0.79 and 0.77, respectively (Table 4). The relationship between MaxN30, QOptN30, and LPOptN30 for each biomass class is illustrated in Fig. 4 . In general, estimated economically optimal N rates at any Ncon within a given biomass class, whether estimated using the linear-plateau models of yield or the quadratic models of partial profit, were nearly identical (Fig. 4). In fact, based on the non-overlap of the SEs of the slopes and intercepts (Table 4), the regression models of MaxN30, QOptN30, and LPOptN30 versus Ncon for a given biomass class were not statistically different.


Figure 4
View larger version (13K):
[in this window]
[in a new window]

 
Fig. 4. Linear regression of optimum growth stage (GS) 30 N fertilizer rates computed using three different methods (LPOptN30, economically optimum N fertilizer rates determined using linear-plateau functions; MaxN30, N fertilizer rates that resulted in the maximum yield; and QOptN30, economically optimum N fertilizer rates determined using quadratic functions) versus GS-30 whole plant tissue N concentration (Ncon). The data are divided into five GS-30 biomass classes.

 

    DISCUSSION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Our primary objective was to determine if the relationship between MaxN30 and Ncon was stable across wheat stands with different GS-30 biomass. To test this hypothesis, we needed to establish plots that differed in GS-30 biomass. Specifically, we wanted to simulate the range in biomass frequently found in wheat fields in the spring across the southeastern USA. This would include fields with high tiller density and biomass at GS-30 typical of timely planted wheat (e.g., following corn or early maturing soybean [Glycine max (L.) Merr.]) and fields with lower tiller density and biomass typical of late-planted wheat (e.g., following late-harvested cotton [Gossypium hirsutum L.]). By using different planting-date and seeding-rate combinations, we were able to create these conditions in our experimental plots. This was a substantial departure from previous research, which generally attempted to evaluate the relationship between MaxN30 and Ncon only in uniform, well established wheat.

Using all the data, the relationship we found between MaxN30 and Ncon was the weakest reported in the literature for soft red winter wheat. However, when the entire data set was subdivided into GS-30 biomass classes, the model relating MaxN30 and Ncon resulted in an R2 = 0.75, which is the highest value we are aware of in the published literature. Wheat in the lowest biomass class had the highest Ncritical and required the highest N30 to optimize yield for any given Ncon. Conversely, wheat in the highest biomass class had an extremely low Ncritical and the lowest N30 requirement to maximize yield for any given Ncon. The GS-30 biomass classes representing 340 to 1000 kg ha–1 (Fig. 3C) and 1500 to 2000 kg ha–1 (Fig. 3E) had a regression models that were similar to those reported by Scharf et al. (1993) and Baethgen and Alley (1989), respectively. The amount of biomass at GS-30 had a direct impact on the prediction of Ncritical and MaxN30. This might explain some of the differences among the relationships previously reported in the literature and might explain some of the unaccounted-for variability in those models. It also indicates that tissue testing, when combined with an estimate of GS-30 biomass, should lead to more accurate GS-30 N rate recommendations.

Baethgen and Alley (1989) suggested that MaxN30 would be best modeled using NUP instead of Ncon. However, their conclusion was based on a data set of 16 observations, and Roth et al. (1989) found that the relationship between MaxN30 and NUP differed across years. Flowers et al. (2003) also concluded that predicted N fertilizer rates based on NUP were highly unreliable. The relationship reported here between MaxN30 and NUP was statistically significant but very weak (r2 = 0.27). For example, consider two wheat fields: the first with biomass of 700 kg ha–1 and Ncon of 55 g kg–1 (lower right corner of Fig. 3C) and the second with biomass of 2500 kg ha–1 and Ncon of 15 g kg–1 (upper left corner of Fig. 3F). Both these hypothetical fields have the same NUP, but MaxN30 would be 0.0 and 122 kg ha–1 for the first and second field, respectively. These results demonstrated that the relationship between MaxN30 and biomass was more complex than can be explained by NUP, which is the simple product of Ncon and biomass. Equation 1 can be rearranged such that:

Formula 3[3]
or

Formula 4[4]
where

Formula 5[5]
One way to interpret Eq. [4] is that MaxN30 can be described as a family of linear models with each using Ncon as the independent variable but where the slope of any given line is a function of the biomass (Eq. [5]). This is consistent with the trend shown in Fig. 3 and demonstrates that there is a more complex relationship between MaxN30, Wbiomass, and Ncon than can be explained by NUP alone.

Models predicting MaxN30 and economically optimal N30 rates (assuming a N/grain price ratio of 10) were not statistically different (Fig. 4). This is consistent with the findings of Baethgen and Alley (1989), who examined lower price ratios than we did. Apparently, the wheat yield response to GS-30 N fertilizer is large compared with current economic considerations.

Figures 3 and 4 indicate that a field with a thick, well developed wheat stand needs a different calibration curve for determining optimum GS-30 N fertilizer rates than one with a thin stand or low GS-30 biomass. This is a critically important consideration in using GS-30 tissue N concentration to optimize N fertilizer rates in farmers' fields. These results also suggest that if a site-specific N fertilizer recommendation system is to be based on GS-30 tissue N concentration, then an estimate of site-specific biomass would be needed to develop an optimal GS-30 N prescription.


    NOTES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
This research was sponsored in part by Initiative for Future Agricultural and Food Systems Grant no. 00-52103-9644 from the USDA Coop. State Research, Education, and Extension Service.


    REFERENCES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 




This article has been cited by other articles:


Home page
Agron. J.Home page
R. P. Sripada, D. C. Farrer, R. Weisz, R. W. Heiniger, and J. G. White
Aerial Color Infrared Photography to Optimize In-Season Nitrogen Fertilizer Recommendations in Winter Wheat
Agron. J., October 15, 2007; 99(6): 1424 - 1435.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Weisz, R.
Right arrow Articles by Farrer, D. C.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Weisz, R.
Right arrow Articles by Farrer, D. C.
Agricola
Right arrow Articles by Weisz, R.
Right arrow Articles by Farrer, D. C.
Related Collections
Right arrow Other Crop Management
Right arrow Best Management Practices
Right arrow Wheat


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Crop Science Vadose Zone Journal
Journal of Natural Resources
and Life Sciences Education
Soil Science Society of America Journal
Journal of Plant Registrations Journal of
Environmental Quality
The Plant Genome