Published in Agron. J. 96:494-501 (2004).
© American Society of Agronomy
677 S. Segoe Rd., Madison, WI 53711 USA
NITROGEN MANAGEMENT
Midseason Nitrogen Fertility Management for Corn Based on Weather and Yield Prediction
Nathan E. Derby*,a,
Francis X. M. Caseya,
Raymond E. Knightonb and
Dean D. Steelec
a Dep. of Soil Sci., North Dakota State Univ., Fargo, ND 58105-5638
b USDA-CSREES-NRE, Washington, DC 20250-2210
c Dep. of Agric. and Biosyst. Eng., North Dakota State Univ., Fargo, ND 58105-5626
* Corresponding author (nathan.derby{at}ndsu.nodak.edu).
Received for publication March 25, 2003.
 |
ABSTRACT
|
|---|
Nitrogen fertilizer applications for irrigated corn (Zea mays L.) based on yield goals established before planting may result in under- or overapplication of N because of weather-induced variations in yield potential from year to year. This study was conducted to develop and evaluate a regression model to predict corn grain yield at a midpoint in the growing season based on the current year's cumulative thermal factors and N fertility levels. The relationship among early-season growing conditions, N fertility, and corn grain yield under continuous corn production and sprinkler irrigation was investigated from 19901995 in southeastern North Dakota. Fertility levels on small plots ranged from 0 to 224 kg ha1 applied N. There was a wide range of cumulative growing degree days (GDD) and evapotranspiration (ET) during the study period. Least-squares regression was used to develop equations to predict, on 10 July, grain yield based on cumulative ET or GDD from 1 May to 10 July and soil N plus applied N. Yields predicted on 10 July corresponded well to individual observed yields (r2 = 0.80), and predicted optimum yield (YOPT) was highly correlated to observed YOPT (r2 = 0.885). The model could provide a season-specific yield potential used to modify N application during the growing season, resulting in fertilizer savings in the extreme years when cool early-season weather limited yield potential or fertilizer increases to take advantage of optimum growing conditions and increase yields. Including meteorological measurements can improve fertilizer management decisions by providing midseason adjustments to fertilizer recommendations.
Abbreviations: ET, evapotranspiration ETR, reference evapotranspiration GDD, growing degree days SA-N, soil nitrogen plus applied nitrogen YOPT, optimum yield
 |
INTRODUCTION
|
|---|
PROPER N FERTILIZER MANAGEMENT is crucial for both economic and environmental reasons. Inadequate N fertility results in low yields and lower economic returns while overfertilization can adversely affect surface and ground water resources and also reduce net economic returns. In addition, overfertilization may also result in higher emissions of nitrous oxidean important greenhouse gas (Kauppi and Sedjo, 2001). Current North Dakota recommendations for fertilizing corn are based on the yield goal concept. This is done by applying 21.4 kg N ha1 for each megagram of corn per hectare expected (yield goal), reduced by the amount of soil test nitrate N (Franzen and Cihacek, 1996). Timely soil test data and reasonable yield goals are required to apply the appropriate amount of N fertilizer using this concept. While the yield goal concept is the accepted approach in North Dakota as well as Minnesota and South Dakota, it may not be suitable for other environments where high soil organic matter and high mineralization reduce the need for N fertilizer to obtain YOPT. Others have stressed the importance of accurately accounting for the amount of soil N, particularly late-spring nitrate N, when making fertilizer recommendations (Binford et al., 1992; Blackmer et al., 1989; Vanotti and Bundy, 1994) and have determined relationships between soil nitrate N and yield. Similarly, Bundy and Malone (1988) found no increase in yield to applied N when soil nitrate levels were high, emphasizing the need for soil test data for accurate fertilizer recommendations.
Previous work indicates that yield goal selection should also be based, to some degree, on soil factors (not simply the maximum desired yield) to produce a more site-specific recommendation (Vanotti and Bundy, 1994; Bundy and Andraski, 1995). Blackmer et al. (1997) recommend, on a site-specific basis, the use of a late-spring soil test to determine the amount of N available before rapid crop N uptake. They also recommended that yield goals no longer be used because much of the N taken up by the crop is supplied by the soil. Additionally, they address weather's role in modifying N recommendation where heavy spring rainfall may leach plant available N from the root zone or warm conditions may increase mineralization and available N.
The relationship of corn maturity to air temperature and cumulative GDD in the Northern Great Plains is well established. Runge (1968) found a high correlation between corn yield and ET (rainfall and maximum daily temperature). Wienhold et al. (1995) noted no difference in yield between two fertility levels when growing season temperatures were below normal and observed the lowest yields in a year with below-average cumulative GDD. Several studies have indicated that a given corn variety requires a minimum number of heat units throughout the growing season to reach maturity (Andrew et al., 1956; Roth and Yocum, 1997; Sutton and Stucker, 1974). This relationship of yield and air temperature (in particular, cumulative GDD) has been used in models to predict corn grain yield (Swan et al., 1990; Bollero et al., 1996; Bauder and Randall, 1982). Taking the prediction one step further, Duchon (1986) used CERES-Maize to predict corn yield at midseason by using historical weather records to fill in the remainder of the year until crop maturity to predict yield. This method assumed N was not limiting.
Knowledge of the effects of N fertility and yearly weather fluctuations in the form of ET or heat unit accumulation needs to be combined to develop better fertilizer management practices. The objective of this study was to develop and evaluate a regression model to predict final corn grain yield at a midpoint in the growing season based on the current year's cumulative thermal factors and N fertility levels. This predicted yield could then be used as a new projected yield goal for that year. Adjustments in the midseason N application recommendations could then be based on this new projected yield goal. The result would be reduced N applications when climatic conditions are less than favorable for maximum yield. Conversely, increased N application may be warranted in years with ideal growing conditions in anticipation of higher yields.
 |
METHODS
|
|---|
Site Description
The site was a 64-ha sprinkler-irrigated field located in southeastern North Dakota (46°2'60'' N, 98°6'36'' W; elevation 397 m). The site was dominated by Hecla loamy fine sand soil (sandy, mixed, frigid Oxyaquic Hapludoll) in the area of the field where the plots described here were located (Fig. 1)
. Eighteen small plots were located in the southeast quadrant of the field, and 18 were located in the northeast quadrant. Each plot was 9.1 m wide (12 corn rows) and 18.3 m long. Six plots were arranged in two-plot-wide by three-plot long blocks. The plots in the southeast quadrant were used for model development while the northeast plots were used for validation. The organic matter content of the top 30 cm is 1.8%. The site was not irrigated before 1989 and was planted to corn (cv. Pioneer 3737) in 19891995. Each year, irrigation was scheduled according to different methods in each quasi-quadrant (see Fig. 1), each based on measured or predicted soil moisture status and crop water use (Steele et al., 2000).

View larger version (31K):
[in this window]
[in a new window]
|
Fig. 1. Map of plot locations and soil types. One block of plots was expanded to show randomized N treatments. A = 0, B = 45, C = 90, D = 135, E = 180, and F = 225 kg ha1 fertilizer N.
|
|
Fertilizer Application
From 19901995, the plots received N fertilizer treatments ranging from 0 (Treatment A) to 225 (Treatment F) kg N ha1 in 45 kg N ha1 increments in single or split applications, depending upon the rate. The B and C treatments (45 and 90 kg N ha1) were applied in mid-June at the six-leaf (V6) stage. All other treatments (D, E, and F) also received 90 kg N ha1 at that time. In mid-July at the 15-leaf stage (V15), an additional 45 kg N ha1 was applied to the D treatments and 90 kg N ha1 to the E and F treatments. In early August at silking (R1), the final application of 45 kg N ha 1 was made to the F plots. Applications at V6 were done by injecting urea ammonium nitrate solution adjacent to the rows. Urea was hand-applied in the rows at V15 and R1 and incorporated by irrigation. Starter fertilizer was applied at planting at the rate of 13 kg N ha1 in 19911993, 48 kg N ha1 in 1994, and 11 kg N ha1 in 1995. This amount was considered in the total N applied.
Soil and Plant Sampling
Soil nitrate N was measured in the center of each plot to a depth of 61 cm each year in early June. Two between-row soil cores from each plot were composited and subsampled. Samples were placed in plastic zip-top bags and frozen until 2 M KCl extracts and nitrate N analysis (EPA Method 353.2, USEPA, 1983) could be done. The sum of soil test nitrate N and applied fertilizer N will now be referred to as soil N plus applied N (SA-N).
Representative areas 4.6 m long and two rows wide were hand-harvested in the center four rows of each plot for final grain yield after maturity (R6). Ears were separated from the plants, dried, shelled, and weighed for yield determination (corrected to 15.5% moisture on a dry weight basis). Optimum grain yields (YOPT) were calculated from the yield vs. SA-N relationship for each year by differentiating the quadratic function of the relationship. A quadratic model for yield vs. SA-N was used since it is the most widely used for describing plant nutrient responses (Blackmer et al., 1989).
Growing Degree Days and Evapotranspiration Calculations
Maximum and minimum air temperature, precipitation, and solar radiation were measured at an automated weather station located 2.7 km from the site. Irrigation and precipitation were also measured on-site. The normal (19601990) maximum and minimum temperatures for 1 May through 30 September were 25.2 and 10.4°C, respectively, and normal 1 May through 30 September rainfall was 316 mm.
The method for GDD calculation was the mean daily air temperature minus 10°C, with bases of 10 and 30°C for the minimum and maximum temperatures, respectively (Cross and Zuber, 1972; Gilmore and Rogers, 1958). Daily ET was calculated by five methods: (i) reference ET (denoted here as ETR) by the method of Jensen and Haise (1963), which uses mean daily air temperature and daily solar radiation to calculate ETR; (ii) estimates of corn ET by the water balance or checkbook (CBK ET) approach of Lundstrom and Stegman (1988), which calculates average corn water use based on daily maximum air temperature and weeks after emergence; (iii) the ET algorithm of Stegman and Coe (1984), which is a modification of the JensenHaise ETR using a fourth-order polynomial ET crop curve based on days past emergence, percentage available soil water, and adjustments for a wet soil surface after rain or irrigation (denoted as Stegman DPE ET); (iv) an ET crop curve developed by Sajid (see Steele et al., 1996) that is another modification of the JensenHaise ETR that uses a fifth-order polynomial crop curve with a days-past-planting base developed from 11 yr of nonweighing lysimeter data (denoted as Sajid DPP ET); and (v) Sajid GDD ET, which is the same as the Sajid DPP ET except that GDD since planting is used as the time base for the ET crop curve. The ET calculations, with the exception of JensenHaise, were specific to southeastern ND.
Growing degree day and ET values accumulated beginning at either 1 May, planting, or emergence until 10 July and until crop maturity or 30 Sept, depending on the method, are listed in Table 1. The cutoff date chosen for midseason GDD and ET accumulation was 10 July because it was the approximate date of the second (V15) fertilizer split application and the latest time in the season to reasonably apply additional N. This is also during a period of very rapid N uptake (Ritchie et al., 1986). A revised yield estimate must be available at this time to make an adjustment, if needed, in the remaining balance of N fertilizer to be applied.
View this table:
[in this window]
[in a new window]
|
Table 1. Cumulative growing degree days (GDD) and evapotranspiration (ET) values through 10 July and end of season for different methods of calculation, observed optimum yield (YOPT), and water application amounts for each year.
|
|
Statistics
Statistical analysis was performed using the SAS JMP software (SAS Inst., 1994) Fit Model and Fit Y by X procedures. Visual inspection of plots of residuals (the difference between observed values and predicted values from the linear regression of yield vs. SA-N) was done to determine if relationships between yield and other variables were linear or quadratic. If the residuals for a linear fit of the data are plotted and the points fall evenly and randomly around the line of zero difference, then the linear method of fitting is valid. If the residual points form a quadratic-type pattern relative to the zero difference line, then the relationship is quadratic. This inspection indicated that SA-N was related to yield in a quadratic fashion, 10 July accumulations of GDD and ET were linearly related to yield (YOPT), and the interaction term of ET x SA-N was linearly related to yield. Hence, the model for prediction of yield included the SA-N term, SA-N squared term, the GDD or ET term, and the ET or GDD x SA-N interaction term in the general form of Eq. [1]
 | [1] |
where ypred is predicted yield, C0 is the y-axis intercept, and C1, C2, C3, and C4 are the regression estimates.
 |
RESULTS AND DISCUSSION
|
|---|
Corn Variety and Irrigation-Scheduling Methods
The corn variety planted each year, Pioneer 3737, is a 98-d relative maturity hybrid that requires 1304°C GDD to reach maturity, which is appropriate for this latitude (Sutton and Stucker, 1974). Therefore, in normal years, the crop not reaching maturity will not limit the yield. In extremely cool years, however, the crop may not mature fully, resulting in reduced yield.
All irrigation-scheduling methods provided adequate but not excess moisture to replenish the soil water deficit as determined by direct measurement or model prediction. Although yield was slightly affected by irrigation scheduling method over the 19901995 study period (at the 0.05 level, F = 4.13), the effect of year (weather) was found to be much more significant (at the 0.01 level, F = 425.8) (Steele et al., 2000). Irrigation and precipitation amounts are listed in Table 1.
Model Development and Evaluation
Fitting the polynomial yield model (Eq. [1]) using individual plot yields and SA-N data for all years (n = 108) and 10 July accumulations of GDD or ET calculated by each of the methods showed that using 1 May through 10 Jul ETR in the model resulted in the highest r2 (Table 2). This model will be discussed further while the models that used GDD and the other methods of ET calculation will not be discussed due to the lack of significance or relatively data-intensive calculations. However, the model parameter estimates and summary statistics for all ET and GDD models are included in Table 2 for those who may be interested in the data.
View this table:
[in this window]
[in a new window]
|
Table 2. Standard regression model estimates and summary statistics for models to predict yield from soil N plus applied N (SA-N) and growing degree days (GDD) or evapotranspiration (ET), via different methods, accumulated to 10 July. The model is of the form: ypred = C0 + C1(SA-N) + C2(SA-N)2 + C3(GDD10July or ET10July) + C4(SA-N x GDD10July or ET10July). n = 108.
|
|
Optimum yields determined from least-squares regression of the observed yield vs. SA-N (the maximum value from the quadratic yield vs. SA-N curve) were positively correlated to the maximum yields predicted midseason in each year of the study, with r2 = 0.89 and p = 0.0051 (Table 3). Predicted maximum yields (ETR method) were not significantly different than observed YOPT's (
= 0.01 level) over all years (t = 0.085, t0.01 = 3.36). Predicted vs. observed optimum SA-N was not significantly correlated (r2 = 0.367, p = 0.279).
View this table:
[in this window]
[in a new window]
|
Table 3. Comparison of optimum yields (YOPT ) from observed N response curve and predicted YOPT from reference evapotranspiration (ETR) model. Number of observation each year was 18.
|
|
Midseason yield prediction over the range of SA-N for each year was compared with the observed yields in Fig. 2
. Predicted yield using the ETR 10 July accumulation values showed increased yield with increased SA-N and yearly variations due to early-season temperatures. Paired t-test statistics indicated that the predicted yields were not significantly different than the observed yields for 1993 and 1994 (t < t0.05). Although statistically the predicted yields were different than the observed yields for years other than 1993 and 1994, the predicted yields were still visually very close to observed yields for all but 1992 and 1995.

View larger version (40K):
[in this window]
[in a new window]
|
Fig. 2. Yield vs. soil plus applied N (SA-N) regressions for observed yield (n = 18) and yield predicted midseason with the reference evapotranspiration (ETR) model for the range of SA-N values. Predicted yields not significantly different than observed yields at the 0.05 level.
|
|
Example of Model Implementation
Since crop growth in 1993 was adversely affected by the cool weather in the spring and the crop did not fully mature before frost on 21 September, we used that year as an example of how the yield prediction model could be used to modify N application midway through the growing season. Based on the farmer's past experience in this portion of the field, the yield goal was set at 12 Mg ha1. Using the yield goal relationship of 21.4 kg N ha1 for each megagram of corn per hectare resulted in 257 kg N ha1 to be applied under current recommendations (Franzen and Cihacek, 1996). Subtracting the average spring soil test nitrate N from the 18 plots for 1993 (27 kg ha1) left a balance of 230 kg N ha1 required. The first fertilizer application of 90 kg N ha1 was applied at V6 for a remaining balance of 140 kg N ha1 to be applied. On 10 July 1993, the cumulative ETR since 1 May was 287 mm. Using ETR = 287 mm along with 257 kg ha1 for SA-N in the yield prediction equation resulted in a new yield goal of 5.2 Mg ha1. With the same yield goal relationship of 21.4 kg N ha1 for each megagram of corn per hectare, only 112 kg N ha1 should be applied during the entire growing season to reach maximum yield due to cool early-season growing conditions compared with the original recommendation of 257 kg N ha1. This meant a reduction in fertilizer N application of 145 kg N ha1 compared with the original recommendation while still achieving YOPT for that year. In other words, the original recommendation would have called for 140 kg N ha1 as the required remaining application while use of the yield prediction model indicated that no further applications were needed. This is validated by the fact that the yield from the C treatment plots, i.e., those receiving no further N at the V15 and R1 splits, was basically the same as the 10 July adjusted yield goal and the observed YOPT (Table 4). Predicted yield (adjusted yield goal) and potential fertilizer savings for all study years are also shown in Table 4. Negative fertilizer savings in 1994 (Table 4) indicated that the yield potential was not met by the level of N fertility and more fertilizer could have been applied to maximize yield.
View this table:
[in this window]
[in a new window]
|
Table 4. Adjusted yearly yield goals based on midseason prediction using reference evapotranspiration (ETR) model, associated fertility levels needed to reach yield goals, and potential fertilizer savings. Fertility levels were based on yield goals where 21.4 kg N ha1 was applied for each 1 Mg corn ha1. The "C" rate plots received no additional fertilizer at V15 of R1 splits.
|
|
Model Validation
The model was applied to yield and SA-N data collected from plots in the northeast quadrant of the field for validation. Table 5 shows the observed YOPT's and quadratic regression statistics for the northeast plots and how they compared to the predicted YOPT. The r2 for predicted vs. observed YOPT's was 0.853, and the relationship was significant at the 0.01 level, which indicated that the model was able to predict yield in a different part of the field. Paired t-test analysis also indicated that the predicted and observed YOPT's were not significantly different (t = 2.84, t0.01 = 3.36).
View this table:
[in this window]
[in a new window]
|
Table 5. Comparison of model-predicted optimum yield (YOPT) to YOPT from the quadratic yield vs. soil N plus applied N (SA-N) response curve from the northeast plots. Number of observations for each year is 18.
|
|
The model was also applied to data collected from an independent producer's field (denoted as Test Area) located approximately 3 km from the study site. The average yield and applied N data from the Test Area field was collected each year and reported in a data set known as the Oakes Irrigation Test Area Agricultural Practices Inventory (API). The API serves as an information base for research and monitoring activities in the Oakes Irrigation Test Area (Esser and Weigel, 1998).
Predicted YOPT's for the Test Area field were calculated for years 19901999 based on the fertilizer applied and ETR values for those years. The predicted YOPT values are compared with the yearly average yields from the Test Area field (Table 6). The linear regression of observed vs. predicted yield for the Test Area field showed a significant correlation at the 0.01 level (r2 = 0.67, p = 0.0038), and paired t-test analysis indicated that the predicted and observed yields were not significantly different (t = 1.96, t0.01 = 2.82). These analyses show that the model is capable of predicting yield at different locations and in different years for different varieties.
View this table:
[in this window]
[in a new window]
|
Table 6. Measured average yield from an independent producer's field (Test Area) for 19901999 compared with model-predicted yield based on applied N and reference evapotranspiration (ETR) data for those years. Linear regression statistics of observed vs. predicted yield are r2 = 0.67, RMSE = 0.934, and p = 0.0038.
|
|
 |
SUMMARY
|
|---|
Corn grain yield was found to be highly correlated to cumulative ET and was adversely affected by cool spring growing conditions in southeastern North Dakota. A midseason yield prediction model, which used N fertility and a measure of early-season growing conditions (i.e., ETR), was developed to determine the attainable yield goal before final fertilizer application. It was shown that N application can be adjusted based on an adjusted yield goal, which reduces N applications in extreme years when growing conditions limit yield potential. It can also result in an increase in N application to take advantage of optimum growing conditions and increase yields. Fine-tuning the N application would not necessarily increase yields but would reduce leaching losses of excess nitrate N in extremely cool years and increase profits by reducing input costs. It is important to note that use of this method of fertilizer recommendation adjustment requires a method of fertilizer application that can be done when the crop is 2 m or more tall, such as fertigation through a sprinkler irrigator. Future studies need to be done to further address the implementation of this method.
 |
REFERENCES
|
|---|
- Andrew, R.H., F.P. Ferwerda, and A.M. Strommen. 1956. Maturation and yield of corn as influenced by climate and production technique. Agron. J. 48:231236.[Abstract/Free Full Text]
- Bauder, J.W., and G.W. Randall. 1982. Regression models for predicting corn yields from climatic data and management practices. Soil Sci. Soc. Am. J. 46:158161.[Abstract/Free Full Text]
- Binford, G.D., A.M. Blackmer, and M.E. Cerrato. 1992. Relationship between corn yields and soil nitrate in late spring. Agron. J. 84:5359.
- Blackmer, A.M., D. Pottker, M.E. Cerrato, and J. Web. 1989. Correlations between soil nitrate concentrations in late spring and corn yields in Iowa. J. Prod. Agric. 2:103109.
- Blackmer, A.M., R.D. Voss, and A.P. Mallarino. 1997. Nitrogen fertilizer recommendations for corn in Iowa. Iowa State Univ. Ext. Publ. Pm-1714. Iowa State Univ. Ext., Ames.
- Bollero, G.A., D.G. Bullock, and S.E. Hollinger. 1996. Soil temperature and planting date effects on corn yield, leaf area, and plant development. Agron. J. 88:385390.[Abstract/Free Full Text]
- Bundy, L.G., and T.W. Andraski. 1995. Soil yield potential effects on performance of soil nitrate tests. J. Prod. Agric. 8:561568.
- Bundy, L.G., and E.S. Malone. 1988. Effect of residual profile nitrate on corn response to applied nitrogen. Soil Sci. Soc. Am. J. 52:13771383.[Abstract/Free Full Text]
- Cross, H.Z., and M.S. Zuber. 1972. Prediction of flowering dates in maize based on different methods of estimating thermal units. Agron. J. 64:351355.[Abstract/Free Full Text]
- Duchon, C.E. 1986. Corn yield prediction using climatology. J. Clim. Appl. Meteorol. 25:581590.
- Esser, D.R., and J.K. Weigel. 1998. Agricultural Practices Inventory for the Oakes Test Area study program. p. 351365. In J. Schaack, A.W. Freitag, and S.S. Anderson (ed.) Proc. 1997 Water Manage. Conf., Fargo, ND. 1619 July 1997. U.S. Committee on Irrigation and Drainage, Denver.
- Franzen, D.W., and L.J. Cihacek. 1996. North Dakota fertilizer recommendation tables and equations based on soil test levels and yield goals. SF-882. North Dakota State Univ. Ext. Serv., Fargo.
- Gilmore, E.C., Jr., and J.S. Rogers. 1958. Heat units as a method of measuring maturity in corn. Agron. J. 50:611615.[Abstract/Free Full Text]
- Jensen, M.E., and H.R. Haise. 1963. Estimating evapotranspiration from solar radiation. J. Irrig. Drain. Div. Am. Soc. Civ. Eng. 89:1541.
- Kauppi, P., and R. Sedjo. 2001. Technological and economic potential of options to enhance, maintain, and manage biological carbon reservoirs and geo-engineering. p. 301343. In Metz et al. (ed.) Climate change 2001: Mitigation. Cambridge Univ. Press, Cambridge, UK.
- Lundstrom, D.R., and E.C. Stegman. 1988. Irrigation scheduling by the checkbook method. Bull. AE-792 (rev.). North Dakota State Univ. Ext. Serv., Fargo.
- Ritchie, S.W., J.J. Hanway, and G.O. Benson. 1986. How a corn plant develops. Spec. Rep. 48 (revised January 1986). Iowa State Univ. of Sci. and Technol. Coop. Ext. Serv., Ames.
- Roth, G.W., and J.O. Yocum. 1997. Use of hybrid growing degree day ratings for corn in the northeastern USA. J. Prod. Agric. 10:283288.
- Runge, E.C. 1968. Effects of rainfall and temperature interactions during the growing season on corn yield. Agron. J. 60:503507.[Abstract/Free Full Text]
- SAS Institute. 1994. JMP user's guide. Version 3.1.5. SAS Inst., Cary, NC.
- Steele, D.D., A.H. Sajid, and L.D. Prunty. 1996. New corn evapotranspiration crop curves for southeastern North Dakota. Trans. ASAE 39(3):931936.
- Steele, D.D., E.C. Stegman, and R.E. Knighton. 2000. Irrigation management for corn in the northern Great Plains, USA. Irrig. Sci. 19:107114.
- Stegman, E.C., and D.A. Coe. 1984. Water balance irrigation scheduling based on JensenHaise equation: Software for Apple II, II+ and IIE computers. Res. Rep. 100. North Dakota Agric. Exp. Stn., North Dakota State Univ., Fargo.
- Sutton, L.M., and R.E. Stucker. 1974. Growing degree days to black layer compared to Minnesota relative maturity rating of corn hybrids. Crop Sci. 14:408412.[Abstract/Free Full Text]
- Swan, J.B., J.A. Staricka, M.J. Shaffer, W.H. Paulson, and A.E. Peterson. 1990. Corn yield response to water stress, heat units, and management: Model development and calibration. Soil Sci. Soc. Am. J. 54:209216.[Abstract/Free Full Text]
- USEPA. 1983. Methods for chemical analysis of water and wastes. EPA-6000/479020 (PB84128677). U.S. Environ. Monitoring and Support Lab., Cincinnati, OH.
- Vanotti, M.B., and L.G. Bundy. 1994. Corn nitrogen recommendations based on yield response data. J. Prod. Agric. 7:249256.
- Wienhold, B.J., T.P. Trooien, and G.A. Reichman. 1995. Yield and nitrogen use efficiency of irrigated corn in the northern Great Plains. Agron. J. 87:842846.[Abstract/Free Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
K. Girma, S. L. Holtz, D. B. Arnall, L. M. Fultz, T. L. Hanks, K. D. Lawles, C. J. Mack, K. W. Owen, S. D. Reed, J. Santillano, et al.
Weather, Fertilizer, Previous Year Yield, and Fertilizer Levels Affect Ensuing Year Fertilizer Response of Wheat
Agron. J.,
November 6, 2007;
99(6):
1607 - 1614.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
N. E. Derby, F. X. M. Casey, and D. W. Franzen
Comparison of Nitrogen Management Zone Delineation Methods for Corn Grain Yield
Agron. J.,
February 6, 2007;
99(2):
405 - 414.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
B. L. Ma, K. D. Subedi, and T. Q. Zhang
Pre-Sidedress Nitrate Test and Other Crop-Based Indicators for Fresh Market and Processing Sweet Corn
Agron. J.,
January 1, 2007;
99(1):
174 - 183.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
N. E. Derby, D. D. Steele, J. Terpstra, R. E. Knighton, and F. X. M. Casey
Interactions of Nitrogen, Weather, Soil, and Irrigation on Corn Yield
Agron. J.,
August 17, 2005;
97(5):
1342 - 1351.
[Abstract]
[Full Text]
[PDF]
|
 |
|