|
|
||||||||
a Div. of Plant Science, Univ. of Missouri, Columbia, MO 65211
b USDA-NRCS, Columbia, MO 65203
c USDA-ARS, Cropping Syst. and Water Quality Res. Unit, Columbia, MO 65211
* Corresponding author (hongn{at}missouri.edu)
Received for publication May 2, 2006.
| ABSTRACT |
|---|
|
|
|---|
EONR (total N applied EONR). This relationship was best modeled by a plateau-linear function, with a low RSN plateau at N rates well below the EONR. A linear increase in RSN began anywhere from 65 kg N ha1 below the EONR to 20 kg N ha1 above the EONR at the three sites with good data resolution near the EONR. Applying N rates in excess of the EONR produced elevated RSN values in all six experiments. Our results suggest that applying the EONR will produce environmental benefits in an economically sound manner, and that continued attempts to develop methods for accurately predicting EONR are justified.
Abbreviations: CP, claypan soil region DL, deep loess soil region ECa, apparent bulk soil electrical conductivity EONR, economically optimal N rate MD, Mississippi delta soil region PNR, the producers' N rates RSN, post-harvest residual soil NO3N
EONR, the difference between total N applied and the EONR 00, 2000 01, 2001
| INTRODUCTION |
|---|
|
|
|---|
Nitrate losses from annual row crops such as corn are greater than from perennial forages (Randall et al., 1997). Most N fertilizer in the USA, and especially in the Mississippi River basin, is applied to corn. Corn is also grown more widely on tile-drained land than other N-receiving crops, creating a rapid pathway for nitrate transport to surface waters. This makes corn the crop that loses the greatest amount of nitrate to water resources in the midwestern USA (Randall et al., 1997; Burkart and James, 1999).
Most percolation in the midwestern USA occurs during the fall/winter/spring recharge months, and this is when potential for nitrate movement out of the root zone is greatest. Timing and magnitude of nitrate movement probably vary significantly across the region due to variation in both average fall/winter/spring percolation and in periods with frozen soil. The percolated nitrate can transfer to surface waters via subsurface tile drainage or baseflow or both (Steinheimer et al., 1998). Nitrate remaining in soil after harvest (referred to as post-harvest residual soil nitrate [RSN]) is probably the main source of nitrate found in percolating water. Fertilizer N (always as ammonia) applied in the fall for the following year's corn crop can also be vulnerable, but must first convert to nitrate (a temperature-dependent process) before it can move.
Reducing nitrate movement from agricultural fields to water resources thus requires an understanding of the conditions that lead to high levels of post-harvest RSN. As fertilizer N and mineralizable soil N are the largest N pools in the Mississippi River basin (Burkart and James, 1999), the answer seemingly lies in how these two N pools are managed.
Evidence from N fertilizer response trials suggests that there is a great deal of variability in the amount of N supplied to a corn crop by the soil. High yields with no N fertilizer applied are not uncommon (Bundy and Andraski, 1995), and the amount of additional yield that can be produced due to N fertilizer is highly variable from field to field (Lory and Scharf, 2003).
This variability in soil N supply (mineralizable organic N and residual mineral N) is rarely accounted for in current N fertilizer management practices. The dominant practice for agricultural producers in the Midwest is to apply a single rate of N fertilizer over whole fields and often whole farms. Extensive research documents that crop N needs vary widely between fields (e.g., Schmitt and Randall, 1994; Scharf et al., 2005) and within fields (e.g., Blackmer and White, 1998; Scharf et al., 2005). Uniform N application rates across fields and farms with varying N needs lead to frequent mismatches between fertilizer N applied and crop N needs. When fertilizer N plus mineralized soil N exceed crop needs, this may lead to the accumulation of RSN (Roth and Fox, 1990; Mitsch et al., 2001), which is susceptible to transfer to water resources.
During the 20th century, annual precipitation increased by 10 to 20% in the Midwest, mostly due to an increase in the number of days with heavy and very heavy precipitation events (NAST, 2001). The increase in precipitation is projected to continue across much of the region during the 21st century (NAST, 2001), so under the current N management practices the potential for nitrate losses to water resources is likely to increase.
Two important scientific syntheses have suggested that applying spatially appropriate amounts of fertilizer N can help to reduce N movement from cropping systems to water resources (Mitsch et al., 2001; Power et al., 2001). An important question in this context is, "What is the appropriate amount?" From the economic standpoint, the EONR is by definition the appropriate amount, but additional evidence is needed as to whether the EONR is the appropriate amount of N from an environmental standpoint. A few investigators have provided evidence for environmental benefits of EONR at the small-plot scale (Andraski et al., 2000; Bélanger et al., 2003). Additional research is needed to determine whether these environmental benefits can be obtained by applying EONR across variable landscapes. Our objectives were (i) to determine if RSN can be reduced if N fertilizer is applied at EONR as it varies across fields, as compared to uniform producers' N rates (PNR) in the same fields, and (ii) to compare RSN levels for N fertilizer rates below, at, and above the EONR.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Six different producers' fields (three in 2000 and three in 2001) were chosen where the row direction appeared to cross the greatest variability in soil type and landscape (Fig. 1, Table 1). Universal Transverse Mercator coordinates and experimental layouts are given in Fig. 6 and 8 of Scharf et al. (2005). All fields had been cropped to soybean [Glycine max (L.) Merr.] the year before the study year. Corn was planted by cooperating producers using their equipment. Planting date, hybrid, planting population, and tillage practices were selected by cooperating producers but were representative of practices used for corn production in these soil regions. The fields in the MD soil region were irrigated using center-pivot irrigation systems.
|
|
|
|
At each sample site within a transect, we used a hydraulically driven probe (Giddings Machine Co., Windsor, CO) to take three soil cores to 0.9-m depth. Each of the three soil cores was sectioned into three depth increments: 0 to 30, 30 to 60, and 60 to 90 cm. We mixed each section of the three soil cores, and then air-dried at 25°C for 72 h, and crushed to pass a 2-mm sieve. Soil nitrate N concentrations were determined using a colorimetric autoanalyzer to perform cadmium reduction followed by the sulfanilamide reaction (Keeney and Nelson, 1982). Soil nitrate N concentrations were then converted to NO3N mass assuming a soil bulk density of 1.5 g cm3.
A total of 120 soil samples were collected (six fields x four transects per field x five plots per transect). Three RSN observations were removed from our dataset because two were missampled (global positioning system data indicated the sample location was not in the correct plot) and one had unreasonably high RSN (256 kg RSN ha1 was associated with 28 kg excess fertilizer N ha1 [total N applied EONR]).
Data Analysis
Corn grain was harvested using a combine instrumented with an AL2000 grain yield monitor (AgLeader Technology, Ames, IA). A detailed description of yield monitor data collection and post-cleaning processes is given in Scharf et al. (2005). Each field was divided into cells 20-m long (in the direction of the corn rows) and 40-m wide containing all treatments. In all experimental cells where RSN was measured, N fertilizer increased corn yield. In each cell, this yield increase was modeled as a function of N fertilizer rate using a quadratic-plateau function. These functions were then used to calculate EONR values using a corn price of $0.08 kg1 and a N fertilizer price of $0.55 kg1, which were typical prices during 2000 and 2001 (USDA, 2005a, 2005b). The EONR is the N rate at which profit is optimized. Yield at this N rate is slightly below maximum yield, but the cost of the additional fertilizer to achieve maximum yield exceeds the value of the additional yield produced. For quadratic-plateau yield response functions, EONR = [(N price/corn price)-b]/2c, where b and c are the linear and quadratic coefficients of the response function, respectively. If the soil nitrate sample transect was across the center of a cell, the EONR of that cell was used. If the transect was along the border of two or three cells, we averaged the EONR of these cells and assigned the mean to this transect. The determination, distributions, and discussion of both within-field and field-wide EONR in these six fields are detailed in Scharf et al. (2005).
The RSN content at the EONR and the PNR was estimated for each field from the relationship between the RSN content and
EONR. The term
EONR is defined as:
![]() | [1] |
We fitted linear, exponential, plateau-linear, and quadratic-plateau functions to describe the RSN content response to
EONR. We chose these functions based on the literature, our knowledge of the relationship between the RSN content and the
EONR, and the shapes of the scatterplots of the RSN content vs. the
EONR. We assessed the goodness-of-fit of these functions based on the magnitude, randomness, and normality of the model-fit residuals. An ideal model would have the smallest residuals that exhibit a random pattern and are normally distributed. We fitted the data in SAS PROC REG (SAS Institute, 2002) for the linear function and SAS PROC NLIN (SAS Institute, 2002) for the exponential, plateau-linear, and quadratic-plateau functions. In this study, we reported all results based on the RSN content at the 0- to 0.9-m depth and mean comparisons were conducted using ANOVA or a paired t test.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
The
EONR in sampled plots ranged from 228 to 231 kg N ha1 averaging 45 ± 11 kg N ha1. The
EONR at 84 out of the 117 sampling sites was greater than zero indicating that about 72% of sampled plots received an overapplication of fertilizer N (Fig. 4).
|
EONR) applied to the fertilized plots in the six fields was 95 ± 7 kg N ha1, so about half of the excess N was recovered as RSN. The fate of the other half of the excess N was probably divided among immobilization in soil organic matter, luxury consumption of N by corn (e.g., Binford et al., 1992), NH3 volatilization from the canopy (e.g., Sharpe and Harper, 1995), loss via denitrification, and in-season leaching to below 0.9-m depth. Uncertainties are associated with our RSN values due to differences between actual soil bulk density and the assumed value of 1.5 g cm3. Average moist bulk density values from soil surveys for the soil map units represented in this study ranged from 1.35 to 1.6 g cm3, suggesting that errors associated with this assumption would be 10% or less.
Residual Soil Nitrate Nitrogen Response to Difference between Total Nitrogen Applied and Economically Optimal Nitrogen Rate
The plateau-linear function provided the overall best description of the relationship between
EONR and RSN (Fig. 5), and was used to estimate the RSN content at both EONR and PNR for each field. The plateau-linear function had the lowest model-fit residuals in four of the six fields. For the MD01 (MD-2001) and CP00 (CP-2000) fields, the plateau-linear function had equivalent and 26% greater residuals, respectively, as compared to the functions with the lowest residuals. The model fitting procedure was problematic with the DL00 (DL-2000), MD00 (MD-2000), and MD01 fields. These fields had large gaps (averaging 120 kg N ha1) in the data near
EONR = 0, while the average gap in the other three fields with good data resolution near
EONR = 0 was 23 kg N ha1. The unconstrained plateau-linear models for the DL00, MD00, and MD01 fields were nearly identical to simple linear models. We felt that this was an artifact due to the data gap at these locations, given the clear plateau-linear behavior at all three locations with better data resolution. We chose to constrain the models of the DL00, MD00, and MD01 fields to have the joint point of the plateau-linear function at 65 kg N ha1 because this was the lowest value observed among the three fields with good data resolution near
EONR = 0. Constraining the joint point in this way reduced R2 by very little (0.010, 0.003, and 0.004 for the three fields) relative to the unconstrained models.
|
EONR for each field, but the strength of such relationships varied among fields as reflected in R2 and p values of the regression analyses (Fig. 5). The response of the RSN to the
EONR exhibited two contrast features: when
EONR < 0, the RSN content was lower with less variability; when
EONR > 0, the RSN content was higher with greater variability. Similar contrast features were observed by others (e.g., Andraski et al., 2000, Fig. 3; Bélanger et al., 2003, Fig. 2).
When averaged over all fields, the RSN content with
EONR < 0 ranged from 2 to 40 kg N ha1 with a standard deviation of 11 kg N ha1 and mean of 17 ± 2 kg N ha1. If the data from the check plots were excluded, the average RSN was 21 ± 4 kg N ha1. When
EONR > 0, average RSN increased to 69 ± 5 kg N ha1, indicating a direct link between excessive N inputs and NO3N accumulation in soil. Variability in RSN was considerable when
EONR > 0, ranging from 4 to 237 kg N ha1 with a standard deviation of 48 kg N ha1. This greater variability is partly due to our limited sampling in that we used only three cores per sample. Because of the likely small-scale spatial variability in RSN (Ruffo et al., 2005), using a small number of subsamples might have resulted in substantial uncertainty. This is especially true when nitrate levels in soil are high. The high variability in RSN with
EONR > 0 creates difficulty in selecting a function to predict the RSN content based solely on the amount of excess fertilizer N. This fact is reflected by low R2 for models of RSN as a function of
EONR at some locations (e.g., the fields labeled MD00 and MD01 in Fig. 5).
Residual Soil Nitrate Nitrogen Response to Difference between Total Nitrogen Applied and Economically Optimal Nitrogen Rate Classes
Figure 5 illustrates the distribution of the RSN content over
EONR classes, and includes
EONR from all transects and fields. The RSN increased with increasing
EONR. When
EONR < 0 (with N applied), average RSN was 21 ± 4 kg N ha1, and was not significantly different from that in the check plots (zero N applied). The plots of this class received an average of 174 ± 6 kg total fertilizer N ha1, which was 28 ± 7 kg N ha1 less than the mean EONR for these plots. This suggests that if fertilizer N is applied at a rate below the EONR, there might be little fertilizer-derived NO3N in soil after harvest even if the N rate applied is high.
The RSN content at the EONR (i.e., the class labeled " = 0" in Fig. 6), estimated for each field using the plateau-linear functions shown in Fig. 5, varied from 12 to 57 with a mean of 33 ± 7 kg N ha1. This mean is well less than the 108 kg N ha1 found in Wisconsin by Andraski et al. (2000, Fig. 3) for corn after corn. Part of this difference may be due to rotation. Bundy (2004) reported that RSN for corn after soybean was 37 kg N ha1 at the apparent EONR (142 kg N ha1), while RSN for corn after corn was 104 kg N ha1 at the apparent EONR (189 kg N ha1). All of our fields had soybean as the previous crop, and agree well with the results of Bundy (2004) for RSN at the EONR for corn following soybean. Average RSN at the EONR was not significantly greater than that in the fertilized plots with
EONR < 0 (p = 0.12), but was significantly greater than that in the check plots (p = 0.003).
|
EONR < 50, average RSN increased to 39 ± 8 kg N ha1, and was significantly greater than when
EONR < 0 (with N applied) (p = 0.09), but was not significantly greater than when
EONR = 0 (p = 0.64). When 50 <
EONR < 100, average RSN increased to 49 ± 6 kg N ha1, and was not significantly greater than when
EONR = 0 (p = 0.23). When
EONR > 100, average RSN increased to 91 ± 7 kg N ha1, which was significantly greater than all groups with
EONR < 100 kg N ha1 (p < 0.006). It should be noted that about 21% of the PNR plots in the six fields received more than 100 kg excess fertilizer N ha1. Consequently, such excessive RSN would mainly be found in a small proportion of the total area (i.e., "hotspots").
Comparisons of Residual Soil Nitrate Nitrogen at Economically Optimal Nitrogen Rate and Producers' Nitrogen Rate
The mean RSN content at EONR averaged over all fields was less by 12 or 15 kg N ha1 (Table 3) than that at PNR (paired t test, p = 0.04 if the plateau-linear functions were used to estimate RSN at PNR, or p = 0.08 if the RSN observations were used). However, in fields where average EONR was greater than or approximately equal to the PNR (i.e., the fields CP00, CP01 [CP-2001], and DL01 [DL-2001] in Table 3), there was little difference in the RSN content associated with the two strategies. In fields where the PNR was at least 43 kg N ha1 above the mean EONR (i.e., the fields DL00, MD00, and MD01), applying fertilizer N at the EONR could reduce the RSN content by 22 kg N ha1 (paired t test, p = 0.003) if the plateau-linear functions were used or 26 kg N ha1 (p = 0.09) if the RSN observations were used. This suggests that identifying whole fields or large portions of a field where N rates can be reduced below the current PNR without economic loss will be a beneficial first step toward reducing RSN and N loss. Additional environmental benefits may be possible by adopting variable-rate N applications.
Timing Effect of Nitrogen Application on Residual Soil Nitrate Nitrogen
Time of N fertilization influenced RSN in only two fields, but with opposite outcomes. For the CP01 field, average RSN with fertilizer N applied at V1 was significantly less by 41 kg N ha1 than when fertilizer N was applied at V7 (paired t test, p = 0.09). The opposite trend was found in the DL01 field where average RSN with fertilizer N applied at V1 was significantly greater by 34 kg N ha1 than when fertilizer N was applied at V7 (paired t test, p = 0.10). The opposite timing effects between the two fields might be associated with the presence of the rainfall events right before and after N applications. In the CP01 field, 3.70 cm of rain accumulated in the 5 d after N application at V1, but there was only 2.20 cm of rain in the 5 d after N application at V7 (Table 2 and Fig. 2). In the DL01 field, there was no rain at V1 until 9 d after N application. At V7, 4.45 cm of rain occurred 1 d before N application and 3.45 cm of rain accumulated in the 5 d after N application. We speculate that significant rainfall events right before and/or after N application might have increased nitrate loss via denitrification and other pathways like leaching via preferential flow, thus causing the opposite timing effects. In the CP fields, high-clay, low hydraulic conductivity subsoils result in poorly or somewhat poorly drained soils, which are vulnerable to saturation and N loss via denitrification. Jokela and Randall (1997, Table 5) reported a similar inconsistent timing effect (between N applications at planting and V8) on RSN. Weather changes and denitrification were considered as two important factors contributing to their inconsistent timing effects.
While timing effects occurred in two fields, they were opposite in direction and no consistent effect was seen. Thus, when averaged over the six fields, there were no significant differences in the RSN content between fertilizer N applied at V1 and V7 (paired t test, p = 0.67). Thus, time of application had little effect on RSN in these six fields. Nitrogen rate was a much more important determinant of RSN.
| SUMMARY AND CONCLUSIONS |
|---|
|
|
|---|
EONR. When
EONR < 0, average RSN was 21 kg N ha1, and was not significantly different than if no fertilizer was applied or if fertilizer N was applied at the EONR. When 0 <
EONR < 50 kg N ha1, average RSN increased to 39 kg N ha1, but was not significantly greater than when
EONR = 0. When 50 <
EONR < 100, average RSN increased to 49 kg N ha1, and was not significantly greater than when
EONR = 0. When
EONR > 100 kg N ha1, average RSN increased to 91 kg N ha1, which was significantly greater than all groups with
EONR < 100 kg N ha1. Applying fertilizer N at the EONR or less can achieve environmental benefits by reducing RSN. Economically optimal N rates at sampling sites varied widely both among and within these six fields, suggesting the need to accurately diagnose EONR at both whole- and sub-field scales. Increasing global N use increases the need to minimize environmental impacts of N fertilizer, and in North America increasing natural gas and fertilizer costs increase the economic need to avoid overapplication of N. Both needs can be addressed by applying the EONR. Further improvements in techniques for diagnosing EONR at the sub-field scale are justified in these production environments. These techniques might be based on crop reflectance sensors, aerial imagery, soil tests, and/or soil/landscape attributes.
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
F. Li, Y. Miao, F. Zhang, Z. Cui, R. Li, X. Chen, H. Zhang, J. Schroder, W. R. Raun, and L. Jia In-Season Optical Sensing Improves Nitrogen-Use Efficiency for Winter Wheat Soil Sci. Soc. Am. J., July 14, 2009; 73(5): 1566 - 1574. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. P. Schmidt, A. E. Dellinger, and D. B. Beegle Nitrogen Recommendations for Corn: An On-The-Go Sensor Compared with Current Recommendation Methods Agron. J., July 7, 2009; 101(4): 916 - 924. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. C. Scharf and J. A. Lory Calibrating Reflectance Measurements to Predict Optimal Sidedress Nitrogen Rate for Corn Agron. J., May 8, 2009; 101(3): 615 - 625. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. J. Boquet, B. S. Tubana, H. J. Mascagni Jr., M. Holman, and S. Hague Cotton Yield Responses to Fertilizer Nitrogen Rates in a Cotton-Corn Rotation Agron. J., March 4, 2009; 101(2): 400 - 407. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. M. Kyveryga, A. M. Blackmer, and J. Zhang Characterizing and Classifying Variability in Corn Yield Response to Nitrogen Fertilization on Subfield and Field Scales Agron. J., February 4, 2009; 101(2): 269 - 277. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Hernandez and D. J. Mulla Estimating Uncertainty of Economically Optimum Fertilizer Rates Agron. J., September 1, 2008; 100(5): 1221 - 1229. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Vadose Zone Journal | ||||
| Soil Science Society of America Journal | Journal of Plant Registrations | The Plant Genome | |||