Laila A. Puntel1*
John E. Sawyer1
Daniel W. Barker1
Peter J. Thorburn2
Michael J. Castellano1
Kenneth J. Moore1Andrew VanLoocke1
Emily A. Heaton2
Sotirios V. Archontoulis1*- 1Department of Agronomy, Iowa State University, Ames, IA, United States
- 2CSIRO Agriculture, St Lucia, QLD, Australia
Historically crop models have been used to evaluate crop yield responses to nitrogen (N) rates after harvest when it is too late for the farmers to make in-season adjustments. We hypothesize that the use of a crop model as an in-season forecast tool will improve current N decision-making. To explore this, we used the Agricultural Production Systems sIMulator (APSIM) calibrated with long-term experimental data for central Iowa, USA (16-years in continuous corn and 15-years in soybean-corn rotation) combined with actual weather data up to a specific crop stage and historical weather data thereafter. The objectives were to: (1) evaluate the accuracy and uncertainty of corn yield and economic optimum N rate (EONR) predictions at four forecast times (planting time, 6th and 12th leaf, and silking phenological stages); (2) determine whether the use of analogous historical weather years based on precipitation and temperature patterns as opposed to using a 35-year dataset could improve the accuracy of the forecast; and (3) quantify the value added by the crop model in predicting annual EONR and yields using the site-mean EONR and the yield at the EONR to benchmark predicted values. Results indicated that the mean corn yield predictions at planting time (R2 = 0.77) using 35-years of historical weather was close to the observed and predicted yield at maturity (R2 = 0.81). Across all forecasting times, the EONR predictions were more accurate in corn-corn than soybean-corn rotation (relative root mean square error, RRMSE, of 25 vs. 45%, respectively). At planting time, the APSIM model predicted the direction of optimum N rates (above, below or at average site-mean EONR) in 62% of the cases examined (n = 31) with an average error range of ±38 kg N ha−1 (22% of the average N rate). Across all forecast times, prediction error of EONR was about three times higher than yield predictions. The use of the 35-year weather record was better than using selected historical weather years to forecast (RRMSE was on average 3% lower). Overall, the proposed approach of using the crop model as a forecasting tool could improve year-to-year predictability of corn yields and optimum N rates. Further improvements in modeling and set-up protocols are needed toward more accurate forecast, especially for extreme weather years with the most significant economic and environmental cost.
Introduction
Over and under fertilization of nitrogen (N) in corn production affects the farmer's profitability and the environment (Shanahan et al., 2008). Predicting the economic optimum N rate (EONR) before crop planting is an ongoing research effort. The challenge persists because of multiple dynamic factors influencing the EONR. In brief, genotypic inputs (cultivars), environment (soil × weather, especially rainfall and its distribution), and management choices (tillage, N application time, etc.) affect soil and crop processes in various ways. The result of all these dynamic processes and their interactions (soil supply vs. crop demand) determine the yield at any N fertilization level (Figure 1). There are several examples in the literature where a single component of the system was studied in detail without acknowledging other system's components and their inherent feedbacks (effect of tillage, effect of residue removal, Kwaw-Mensah and Al-Kaisi, 2006; Coulter and Nafziger, 2008). Several tools and methodologies have been developed over time to assist farmers with N rate decisions (e.g., yield goal approach, Stanford, 1973; soil nitrate test, Bundy and Andraski, 1995; Shapiro et al., 2008), while other new tools such as sensor technologies and simulation models are currently being developed and tested (Scharf, 2015; Banger et al., 2017).
Figure 1. Overview of the main factors influencing the economic optimum nitrogen fertilizer (EONR) rate and their interactions. Soil organic matter (SOM).
Most of today's N-rate decision tools focus on a single component of the soil-plant system to predict EONR rather than utilizing all of the N dynamics and interactions that occur among processes in the system (Figure 1; Arbuckle and Lasley, 2013). For example, the single time soil testing approach around corn at 6th leaf stage provides an estimate of soil N supply (Magdoff et al., 1984; Blackmer et al., 1989; Binford et al., 1992). However, the prediction accuracy and usage of this approach is low because the soil nitrate varies greatly in time and space, especially in rainfed production regions with unpredictable rain events (Jemison and Lytle, 1996; Ma et al., 2007; Arbuckle and Lasley, 2013). Other examples are crop sensors and remote sensing technologies that provide an estimate of crop N status during the season (Yuan et al., 2016). Although promising as a N diagnostic tool this approach has not yet been widely adopted by farmers (Mamo et al., 2003; Scharf et al., 2005; Hawkins et al., 2007) because they typically delay N application until the plant N status can be reliably determined. That requirement increases the risk of not being able to apply N or adds extra application costs (Van Es et al., 2006; Tremblay et al., 2012; Franzen et al., 2016).
Another N tool, the yield goal (Stanford, 1973; Stanford and Legg, 1984), requires inputs such as grain yield, N concentration, nitrogen use efficiency, and N credits coming from manure, legume crops, and soil organic matter to determine corn N rates. These inputs are difficult to estimate because they are derived properties (outputs) of many interactive processes occurring simultaneously within the soil-plant-atmosphere system (Figure 1). Thus, farmers use prior knowledge or guesswork to provide these inputs. According to Lory and Scharf (2003) and Shanahan (2011), the yield goal approach usually results in over fertilization as a form of “insurance” against uncertain soil N supply. In contrast, the Maximum Return To Nitrogen approach (MRTN; Sawyer et al., 2006) requires simple inputs such as location, rotation, crop and fertilizer prices. By using an extensive experimental network of derived N-response datasets across the USA Midwest, it provides farmers with a N-rate recommendation per geographical areas (Sawyer et al., 2006).
In contrast to the aforementioned tools, process-based cropping system models that account for different soil-crop processes, and their interactions with management, cultivar and environmental conditions (Figure 1), have demonstrated capabilities explaining causes of EONR variability and perform scenario analyses (Thorp et al., 2007; Gowda et al., 2008; Nangia et al., 2008; Puntel et al., 2016). However, most of the model applications in N research have been applied ex-post (after harvest), which is of limited use to farmers (Kersebaum et al., 2005; Nendel et al., 2013).
Successful crop yield forecasting approaches using process-based models and historical weather in the USA (Morell et al., 2016; Togliatti et al., 2017) and Australia (Carberry et al., 2009), have indicated the potential to complement the explanatory power of cropping system models with the forecasting component that is needed for N decisions. Theoretically, by running a process-based model for different N rates using actual and historical weather, data the model can provide end-of-season yields (the yield-N response curve) at any time during the growing season. Thus, the EONR could be estimated as early as planting time, providing a new approach to make N rate decisions and supporting information on crop yields, N supply, and N demand (Figure 1). To our knowledge, the validity of this approach and the uncertainty around the year-by-year (annual) EONR prediction at planting and during the growing season have not been previously investigated.
A critical aspect in crop yield forecasting (and ultimately N-forecasting) is the use of weather information to drive model simulations. Traditionally, historical weather data is used to fill the unknown weather for the remainder of the growing season and calculate yield probabilities (Hammer et al., 1996; Quiring and Legates, 2008). Thus, forecasting yield and optimal N-rate at planting time is a challenge given the uncertainty in weather. As weather information becomes available during the season, the uncertainty around crop yield prediction decreases, but it follows different patterns from year to year and cropping systems (Archontoulis et al., 2016; Togliatti et al., 2017). Furthermore, use of analogous historical weather years (i.e., weather with similar precipitation or temperature patterns as the year to be forecasted) instead of using 35-year record has been used as an approach to improve yield forecast (i.e., Hammer et al., 2001; Hansen et al., 2004). The time of the forecast may also affect yield and EONR predictions and uncertainty patterns.
In this study, we tested the hypothesis that the use of a calibrated cropping system model coupled with an assembly of actual and historical weather datasets can predict EONR as early as planting time with similar accuracy as the prediction at harvest with known weather (ex-post). We build upon Puntel et al. (2016) in which the APSIM model was calibrated using 16 years of corn yield response to N data from two crop rotation systems at a site in central Iowa, USA. Our objectives were to: (1) evaluate the accuracy and uncertainty of corn yield and EONR predictions at four forecast times (planting time, 6th and 12th leaf, and silking crop physiological stages) compared with the observed and simulated values at harvest; (2) investigate whether the use of selective historical weather records (e.g., years with similar precipitation) will increase accuracy of yield and EONR predictions as opposed to the 35-year historical record; and (3) quantify value added by the crop model in predicting annual EONR and yields using the site-mean EONR and the yield at the EONR (average of yearly values) to benchmark predicted values.
Methods
Experimental Data and Site Description
We used 16 years (1999–2014) of corn yield response to N fertilizer rate data from a field experiment conducted in central Iowa, USA (details in Puntel et al., 2016). The experiment was designed to study the effect of five N fertilizer rates (0, 67, 134, 201, and 268 kg N ha−1) on corn yield in continuous corn (CC) and corn following soybean (SC) cropping systems. Application of N was either pre-plant or side-dress; see details in Puntel et al. (2016). Corn grain yield was reported at 15.5% moisture content. The corresponding EONR values per year and rotation were calculated using measured yields (see section Estimation of the Annual Economic Optimum Nitrogen Rate). The weather at the experimental site is humid continental with warm rainy summers with an average annual precipitation of 900 mm and annual temperature of 9°C. Over the 16-year experimental period, crops experienced warm and wet conditions (3 years), cool and wet conditions (3 years), warm and dry conditions (5 years), and cool and dry conditions (5 years; Figure S1a). The soil at the site is a deep fertile loamy (Clarion soil series) with topsoil organic matter of 3.4% and profile plant available water of 250 mm.
The APSIM Modeling Platform, Set Up, and Calibration
The APSIM model is an open-source advanced simulator of agricultural systems that combines several process-based models in a modular design (Holzworth et al., 2014). In this study, we used the recently calibrated version of the model for this site with no additional changes (see Puntel et al., 2016 for detailed calibration information). Model performance is also provided in this study (see yield predictions at maturity, results section). The simulation process was continuous, starting in 1999 and ending in 2014 without annual re-initialization of inputs to capture carryover effects on soil N, water, residue, and soil organic matter dynamics from 1 year to the next. The following APSIM models were used: corn and soybean crop models (Keating et al., 2003), Soil N (soil N and C cycling model; Probert et al., 1998), SWIM (soil water model using the Richard equation and fluctuated water tables; Huth et al., 2012); SURFACEOM (residue model; Probert et al., 1998; Thorburn et al., 2001; Thorp et al., 2005), soil temperature (Campbell, 1985), and the following management rules: planting, harvesting, fertilizer, tillage, and rotations (Keating et al., 2003).
Forecasting Yields and EONR
To forecast yields and EONR in each study year (1999–2014) we combined actual weather data up to a specific crop stage and historical weather data thereafter (Figure S2). The following forecast times were considered in this study: planting, 6th leaf (V6), 12th leaf (V12), and silking (R1; Abendroth et al., 2011). We selected these times because farmers in rainfed production regions can make delayed N applications during vegetative and early reproductive stages (R1) without a significant negative effect on yield (Scharf et al., 2002). In this environment, there is not much evidence of yield response to N fertilizer application after corn silking, thus we did not explore yield and N forecasts during grain filling period.
Except for planting time, which only includes observed weather data until the planting date, the rest of the forecast times used observed in-season weather data until the forecast time (Figure S2). Historical daily weather data from 1980 to 2015 (35 years, maximum available for this site) was obtained from the Iowa Environmental Mesonet (2014) to fill the unknown weather. The combination of 35 years of simulations, four forecast times, 16 study years, five N rates, and two crop rotations resulted in more than 22,000 simulations. We calculated the yearly mean and standard deviation of yield and EONR from the 35 simulated values per year and forecasting time similar to Togliatti et al. (2017). The simulated site-mean EONR was calculated by averaging individual annual estimates of EONR and yield at EONR (YEONR) for each rotation, and then calculating the associated standard deviation (SD; across years mean).
Once the simulations were completed, we grouped the data to represent the following five weather scenarios and quantify the impact on yield and EONR forecasts:
• Use simulated yields from calendar years with similar annual precipitation and temperature patterns as the year of interest (e.g., 2010 and 1984). These groups of years can be found in Figure S1a.
• Use simulated yields from calendar years with similar summer (June to August) precipitation and temperature patterns as the year of interest. These groups of years can be found in Figure S1b.
• Use simulated yields from the five calendar years prior to the study year (e.g., for 2010 we used 2005–2009 years).
• Use simulated yields from the 10 calendar years prior to the study year.
• Use simulated yields from the 20 calendar years prior to the study year.
Scenarios I and II were investigated because previous findings showed that the used of analogous weather years instead of the full weather record could increase accuracy of yield predictions and thus, derived EONR predictions (Hammer et al., 2001; Hansen et al., 2004). Scenarios III to V were investigated to quantify the impact of using limited amount of weather data because long term (35-years) weather data are not available in every location (Hansen et al., 2004; Grassini et al., 2015).
Data Analysis
Estimation of the Annual Economic Optimum Nitrogen Rate
The relationship between observed or simulated yield and the five N rates was fit using the quadratic and quadratic-plus-plateau model following the methodology described in Puntel et al. (2016). Models were deemed significant at p < 0.05 and the equations with the smallest sums of squares and largest R2 were selected. The EONR and YEONR was calculated from the N response equations by setting the first derivative of the fitted response curve equal to a common price ratio of 5.6:1 N: corn grain price (US$ kg−1 N: US$ kg−1 grain) ratio during the study years (Cerrato and Blackmer, 1990; Bullock and Bullock, 1994). Variations in N: grain price ratios during the period of this study (1999–2014) were not taken into consideration because of its potential confounding effect on the simulated EONR and YEORN evaluation. Different price ratios will need to be considered for historical evaluation (Amatya et al., 2008; Sawyer, 2015).
Both the observed EONR and YEONR for each year are presented in Puntel et al. (2016).
Statistical Evaluation of Model Performance
To evaluate the APSIM model simulations goodness of fit, we used graphical and statistical methods. For the statistical evaluation, we computed the root mean square error (RMSE), and the relative root mean square error (RRMSE; see equations in Archontoulis and Miguez, 2015) between observed and predicted values. The lower the value of RMSE and RRMSE the better the model performed. In this study, we considered RRMSE ≤15% as “good” agreement; 15–30% as “moderate” agreement; and ≥30% as “poor” agreement (Liu et al., 2013; Yang et al., 2014).
We quantified the accuracy of yield and EONR predictions by comparing the closeness agreement between the simulated and observed means. The accuracy of EONR prediction was considered good when the error was < ± 30 kg N ha−1. This threshold represents the historical and current suggested N fertilization rate range for corn (Voss and Shrader, 1979; Sawyer et al., 2006).
The uncertainty around the simulated yields and EONR was calculated as the standard deviation of the 35 estimates (Togliatti et al., 2017). The standard deviation characterized the range of values within which the mean prediction is asserted to lie. We used the standard deviation among forecasting times to evaluate the impact of known weather on the uncertainty around the mean prediction for EONR and yield. Finally, we evaluated the proximity of the observed EONR and yield values to the acceptable range within the standard deviation.
To calculate the site-mean EONR and YEONR, we averaged annually observed values for each rotation similar to Puntel et al. (2016). These site-mean values were used to benchmark model predictions, i.e., above, below or at average. We then counted the number of years in which the model correctly predicted the direction being above, below or at average of EONR and YEONR values at planting time. In each year, we calculated the absolute differences between the observed and the predicted EONR and we counted the years with an average difference of ± 30 kg N ha−1.
We also examined APSIM predictions of EONR in extreme N rate years defined as years where the observed EONR is at least 30 kg N ha−1 above or below the site-mean EONR. We specified three categories: (1) 30 kg N ha−1 greater than the site-mean EONR (high N need), (2) within ± 30 kg N ha−1 of the site-mean EONR, and (3) 30 kg N ha−1 lower than the site-mean EONR (low N need). In addition, we explored the performance of the model by measuring the agreement (correlation, R2) between the observed and predicted EONR and YEONR at each forecasting time and by weather conditions based on Figure S2.
Statistical Evaluation of Weather Scenarios
The impact of the weather scenarios on yield and EONR predictions was evaluated by calculating the RRMSE for each scenario and then subtracting the RRMSE value of each scenario from the RRMSE of the standard approach using the 35-years of weather data. A positive difference means that a scenario (I–V) performed worse than the 35-year standard approach and a negative difference means that selection of weather years was better than the 35-year standard approach.
Results
Yield Predictions at Different Forecasting Times Using 35 Years of Historical Weather
Yield prediction at different forecast times did not significantly change during the growing season compared with yield prediction at maturity (Figure 2). On average, the RRMSE decreased by 1.2% from planting to maturity. Across five N rates, two crop rotations, and 16 years APSIM explained 77% of the observed variability of corn yield (Figure 2) when predictions were made early in the season (planting and V6 stage), about 79% of the variability during mid-season (V12 and R1 stages), and 81% of the variability at the end of the season (R6 stage). The simulated mean yield at the four forecasting times (with unknown weather after the corresponding growth stage) was similar to the final yield prediction at R6 simulated with known weather (Figure 2).
Figure 2. Simulated vs. observed corn yield (Top panel), economic optimum N rate (Bottom panel, EONR), and yield at EONR (Central panel, YEONR) at different corn stages using 35-yr weather data. In the top panels data presented by N-rate while in the bottom panels data presented by cropping system: CC, continuous corn and CS, soybean-corn. The V6 and V12 are the 6th and 12th leaf stage, respectively; R1 and R6 are silking and physiological maturity stages, respectively (Abendroth et al., 2011).
At early stages (planting to V12), APSIM was able to explain ~40% of the variability in observed YEONR and 57% of the variability at maturity (Figure 2). The APSIM model was better at explaining the observed variability in the YEONR during warm weather years than cold weather years across the different forecasting times (R2 ~ 0.8 vs. 0.5, respectively; Table S3). Absolute differences between observed and predicted YEONR were 1.55 and 0.63 Mg ha−1 for cold and warm weather years, respectively. Overall model performance improved (R2 = 0.6; data not shown) when looking at this relationship while excluding the 2008, 2010, 2012, and 2014 extreme weather condition years in the 16-years of experiments (Figure S2).
Between the two cropping systems, yield predictions were slightly better for SC than CC rotations across all four forecasting times (RRMSE of 12.9 vs. 14.2%, respectively; Table S1). The standard deviation of the mean yield, a measure of the uncertainty around yield predictions, did not decrease from planting to R1 (Figures S3, S4). At maturity, there is no uncertainty (standard deviation = 0) because the actual weather was known (Figures S3, S4). In agreement with yield predictions, the uncertainty around the predicted mean YEONR did not decrease from planting to R1 in a significant number of the years (Figures 3, 4). However, the mean deviation (precision) of yield predictions across all years was 27% lower for R1 than at planting (Figures 3, 4).
Figure 3. Predicted and observed yield at the economic optimum N rate (YEONR) for continuous corn. The connected red dots indicate APSIM model predictions at planting time, V6 (6th leaf), V12 (12th leaf), and R1 (silking) stages using actual and historical weather data (vertical bars indicate standard deviation, n = 35). At maturity (R6), the weather is known so there is no uncertainty in yield prediction. Blue squares represent the observed YEONR. The green dashed line represents the mean YEONR for the study period at this site.
Figure 4. Predicted and observed yield at the economic optimum N rate (YEONR) for soybean-corn. The connected red dots indicate APSIM model predictions at planting time, V6 (6th leaf), V12 (12th leaf), and R1 (silking) stages using actual and historical weather data (vertical bars indicate standard deviation, n = 35). At maturity (R6), the weather is known so there is no uncertainty in yield prediction. Blue squares represent the observed YEONR. The green dashed line represents the mean YEONR for the study period at this site.
EONR Prediction at Different Forecasting Times Using 35 Year of Historical Weather
In contrast to yield predictions, the EONR predictions were less accurate (Figure 2), and at different forecasting times showed different patterns across rotations and years (Figures 5, 6). In some cases, the EONR prediction was more accurate at planting or V6 than at maturity (e.g., 2004, Figure 5). In other cases, the EONR prediction was more accurate at maturity than early in the season, but the differences in terms of actual value were small (e.g., Figure 5, 2007, the difference between simulated EONR at planting and maturity was 18 kg N ha−1 for CC). Overall, the mean deviation (precision) of EONR predictions across all years was 25% lower for R1 than at planting (Figures 5, 6). The forecasted site-mean EONR and the standard deviation was similar to the observed values, but low for SC and high for CC (Table 1).
Figure 5. Predicted and observed economic optimum N rate (EONR) for corn-corn. The connected red dots indicate APSIM model predictions at planting time, V6 (6th leaf), V12 (12th leaf), and R1 (silking) stages using actual and historical weather data (vertical bars indicate standard deviation, n = 35). At maturity (R6), the weather is known so there is no uncertainty in yield prediction. Blue squares represent the observed EONR. The green dashed line represents the mean EONR for the study period at this site.
Figure 6. Predicted and observed economic optimum N rate (EONR) for soybean-corn. The connected red dots indicate APSIM model predictions at planting time, V6 (6th leaf), V12 (12th leaf), and R1 (silking) stages using actual and historical weather data (vertical bars indicate standard deviation, n = 35). At maturity (R6), the weather is known so there is no uncertainty in yield prediction. Blue squares represent the observed EONR. The green dashed line represents the mean EONR for the study period at this site.
Table 1. Site-mean and the associated standard deviation economic optimum N rate (EONR, units: kg N ha−1) for each crop rotation, forecast time, and the measured site-mean across years.
At early forecasting times (planting, V6, and V12), the absolute average differences between observed and predicted EONR for both rotations were lower for warm than for cold weather years (37 vs. 50 kg N ha−1, respectively; Table S3). The model better explained observed variability in the EONR in cold/dry than warm/wet seasons across forecasting times and crop rotations (R2 = 0.32).
In terms of uncertainty, the standard deviation of the mean EONR prediction did not show a consistent pattern of decrease during the growing season (Figures 5, 6). Overall, the EONR predictions were more accurate in CC (RRMSE = 25%) than SC (RRMSE = 45%) across forecast times (Figure 2).
The APSIM model predicted EONR with an average error of ±38 kg N ha−1 in 62% of the study cases (n = 31). Prediction error was below the threshold value of ±30 kg N ha−1 in about 50% of the cases.
Assessing the Impact of Different Weather Scenarios
Using 35 years of weather data as input to the model resulted in marginally lower RRMSE (on average 1.4%) values and therefore better yield predictions across forecast times compared to the use of 5, 10, or 20 years of weather data or weather data that included only years with similar weather patterns (Figure 7A). Similar results were found for the EONR prediction; that is, use of 35-year weather data performed better than all other scenarios (Figure 7B). Across all forecasting times, RRMSE was 3.5% lower when using 35-years of weather data. Among the scenarios compared against the use of 35-years, the use of 20 years of weather data had the lowest RRMSE for yield and EONR at all forecasting times (Figure 7A). The exception was the yield predictions at V12 and silking stage where the RRMSE was the lowest for scenarios with 5, 10, or 20 years of weather data (Figure 7A).
Figure 7. Differences in relative root mean square error (RRMSE) between weather scenarios (I-V) and standard approach (use of 35-years of weather data) for predicted corn yield (A) and economic optimum N rate (B) at four forecasted times: planting, V6 (6th leaf), V12 (12th leaf), and R1 (silking) stage. The positive values indicate that using the 35-years of weather history is better than use of other weather scenarios (I–V). Continuous corn (CC) and soybean-corn (SC). Vertical bars represents the standard deviation.
Independent of the weather scenarios, RRMSE values for EONR predictions were four times higher compared to yield predictions. The selection of a specific weather scenario tended to have a greater impact on predictions early in the season than toward the end-of-season (Figure 7). For example, at flowering the RRMSE was on average across weather scenarios 0.8 and 2.3% lower than RRMSE at planting for yield and EORN, respectively.
Comparison between APSIM and the 16 Year Site-Mean EONR and YEONR
We use the site-mean EONR and YEONR (see green horizontal line in Figures 3–6) to benchmark the direction and magnitude of error in annual APSIM predictions. Across 16-years in CC and 15-years in SC (n = 31 cases), the observed EONR values were in 11, 14, and 6 cases above, below and at site- mean EONR values respectively. APSIM predictions of EONR at planting time were in 11, 16, and 4 cases, above, below and at the site-mean EONR values respectively (Figures 5, 6). In 19 of 31 (62%) cases, APSIM correctly predicted the direction of annual EONR being above, below, or at average. For those cases, the average error was 38 kg N ha−1 which represents 20 and 25% error base on the average N rate for CC (observed mean = 184 kg N ha−1) and SC (149 kg N ha−1), respectively.
The observed YEONR values were in 14, 11, and 6 cases above, below and at site-mean YEONR values respectively, across rotations and years. APSIM predictions of YEONR at planting time were in 8, 14, and 9 cases, above, below and at the site-mean YEONR values respectively (Figures 3, 4). In 21 of 31 (67%) cases, APSIM correctly predicted at planting time the direction of annual YEONR being above, below, or at average.
In low or high N need years the simulated EONR deviated more from observations than in average N need years (Figures 5, 6; Table 2). In 2 out of 16 years (extreme wet years; 2008 and 2010) the APSIM model greatly underestimated the EONR when > 30 kg N ha−1. In low fertilization need years, APSIM over predicted the N rate on average by 36 and 69 kg N ha−1 in CC and SC rotation, respectively (Table 2). Overall the distribution of the differences between APSIM predicted EONR and observed EONR was skewed to the right (~ + 30 kg N ha−1) for CC and skewed to the left for SC, underestimating the observed EONR (Figure 8).
Table 2. Comparison of economic optimum nitrogen rate (EONR) estimated by APSIM to the site mean EONR and observed EONR for continuous corn and soybean-corn over three categorical EONR ranges.
Figure 8. Distribution of differences between simulated economic optimum N rate (EONR) at planting time, V6 (6th leaf), V12 (12th leaf), R1 (silking), and R6 (maturity) growth stages, and the site-mean EONR minus the yearly observed EONR for continuous corn (CC) and soybean corn (SC) rotation.
Discussion
We demonstrated an alternative way to utilize the power of mechanistic cropping systems models to assist N-rate decisions in real time as opposed to current ex-post use of models that are of little interest to farmers (Quiring and Legates, 2008; Thompson et al., 2015). This is the first study that concurrently forecast corn yields and their N requirements in the USA Corn Belt.
Yield Forecasts
This study showed that end-of season corn yields and YEONR can be predicted within an acceptable error (RRMSE < 17%) as early as planting time in most of the study years (Figures 2–4 and Figures S3, S4). The 16-years of data used to test our yield forecasts accounted for different weather years (wet, dry, warm, cold, Figure S1a) and management practices, which further increased our confidence that our process-based forecasting approach is robust (He et al., 2017). This means that our science-based yield forecasting approach has the potential to inform corn producers in the USA Midwest that are using the yield-goal approach to predict N-rates by “guessing” end-of-season yields (Arbuckle and Lasley, 2013; Raun et al., 2017; Morris et al., 2018). Further testing of our approach across multiple environments is needed before application. This would minimize expected errors that are often encountered when using a conventional yield approach based on yield averages across years (Raun et al., 2017).
In most cases, the simulated yield at maturity was within the standard deviation of yield prediction at planting (Figures S3, S4). The uncertainty in yield predictions did not substantially decrease from planting to silking even though that the weather uncertainty decreased (Figure S3). We believe that if we had run the forecast at 150°C-days (or 15 calendar days) after silking we would have seen a decrease in the uncertainty as two of the major determinants of corn yield (kernel number and potential kernel weight) is set (Andrade et al., 1999; Borrás and Vitantonio-Mazzini, 2018). Previous studies have shown a decrease in the uncertainty around corn yield predictions at or about 150°C-days after silking (Thornton et al., 1997; Hansen et al., 2004; Togliatti et al., 2017). By R1 stage, about 50% of the weather is still unknown and the weather during grain fill period has substantial effects on yield. This weather uncertainty introduces yield variability in model predictions as different weather variables affects various plant and soil processes (and thus the final product that is yield) in different ways that are also phenologically time dependent (Semenov et al., 1993). For a thorough sensitivity analysis of weather effects on APSIM simulated grain yield for central Iowa, USA we refer to Togliatti et al. (2017).
In two of the years, yield predictions at maturity deviated substantially from the simulated yields at planting (2004 and 2008 in both rotations; Figures 3, 4). In these years, yield predictions changed during the season as a response to extremely low temperature conditions during grain fill compared to the 35-year average temperature (Figures S2A, S5). The model responded to this weather event by increasing the length of grain fill resulting in higher yield predictions at maturity than at early stages. However, in our study the response of the model was not enough to match the observed yields; perhaps use of a newer maize version (Soufizadeh et al., 2018) of the APSIM model may capture these dynamics better. Nevertheless, this example shows that use of crop model offers both predictability and the reasons behind yield predictions (Banger et al., 2017).
EONR Forecast
The APSIM model predictions at planting time were directionally correct in 62% of the study cases. This means that our N forecasting approach is promising and has future potential to directly and/or indirectly aid N rate decisions by providing a year-to-year opportunity to adjust N rate recommendations early in the season (Figures 5, 6). In terms of prediction accuracy, our modeling approach forecasted annual EONR values at all forecast times with an error range of ±38 kg N ha−1 in about 62–69% of the simulated cases (Figures 5, 6, 8). These results emphasize the difficulty in predicting annual EONR in corn production. Given a potential threshold error range of ±30 kg N ha−1, this means that further improvements are needed in modeling algorithms as well as development of more precise soil and crop management inputs to the model.
To our knowledge, there is no other tool that can provide both yield and EONR forecasts as early as at planting time. The soil nitrate test and remote sensing approaches are in-season tools. Pre-season tools such as the yield-goal approach rely on yield guesses (Van Es et al., 2006; Sela et al., 2017); while the MRTN tool incorporates N response variation and accounts for price fluctuations, but does not directly adjust for year-to-year variability (Sawyer et al., 2006).
As illustrated in Figure 1, there are many dynamic and interactive factors involved in EONR prediction (Scharf, 2001; Mamo et al., 2003; Scharf et al., 2006; Dhital and Raun, 2016). Every year, management, environment, and genotype interactions determine the shape of the yield response to N fertilizer relationship that is used to estimate EONR. Our modeling approach accounted for the majority of the factors by simulating residue decomposition, soil organic matter dynamics, changes in soil temperature, moisture, and nitrate levels during fallow periods and accounted for management such as planting date, row spacing and date of N application. All these factors affect soil N supply and crop N demand. Low EONR predictions from the model were usually associated with high estimates of soil residual nitrate levels at planting and delays in planting. But, there are many interactions that simultaneously occur among soil-plant processes within the system that makes it difficult to simplify and generalize the reasons of yield and EONR variability (see Figure S5 and Table S2). For example, in a previous work using APSIM model, we found that spring precipitation was related to N losses and EONR (the higher the spring precipitation the higher the EONR; Puntel et al., 2016). In contrast, a later study showed that the amount of precipitation had little influence on simulated grain yield in central Iowa because of the effect of groundwater tables (Togliatti et al., 2017).
The inter-annual prediction accuracy of the EONR forecast was lower than that of the yield forecast (Figures 3–fig66, Figures S3, S4). One reason is that EONR values are not a direct output of the model but a result from a regression analysis of model outputs (yields) which incorporates further uncertainty (Puntel et al., 2016). Across years, the simulated site-mean EONR value by the model was similar to the observed site-mean (Table 1). We believe that by running the model sequentially, annual over-and under-predictions are canceled out and the model is able to predict the site-mean values more accurately.
Future Improvements Toward More Accurate N Rate Forecasts
This study identified three areas for future research: (1) simulation set up, continuously vs. annual reset; (2) YEONR and EONR predictions in extreme years, and (3) weather data to drive simulations.
Our simulation protocol is characterized by a high degree of difficulty as we ran the model sequentially from 1999 to 2014 to avoid annual re-initialization of soil input parameters and to capture the carry-over effects on N dynamics (Constantin et al., 2011; Basso and Ritchie, 2015). That approach was followed because of the lack in data to update the model every year and because there was evidence from other studies that models can simulate yields and organic matter trends at different N rates in the long term (Ma et al., 2007; Puntel et al., 2016). The drawback of this approach is that simulation errors can accumulate over time and affect the next year's simulation, especially if the model fails to predict crop yield and N dynamics in one of the years (Salo et al., 2016).
We believe that yield and EONR forecasts at planting time can be further improved if additional information on soil nitrate, water, and surface residue was available to eliminate model uncertainties in simulating carry-over effects during fallow periods (Hansen et al., 2006; Carberry et al., 2009; Ines et al., 2013; Yin et al., 2017). Previous modeling work in Iowa has shown that the available soil water and N status at planting time largely affected predicted mean yields and the range of yield level probabilities (Archontoulis et al., 2016). Therefore, when the model is primarily used for accurate year-by-year yield and EONR forecasting purposes, an annual reset approach may be more appropriate than sequential (Ines et al., 2013; Iqbal et al., 2017). This information can be derived from emerging technologies such as remote sensing, and other common tools such as soil N testing (Basso et al., 2001; Jin et al., 2017; Reimer et al., 2017). Testing and possible improvements in the APSIM surface organic matter model for corn and soybean residues in future studies may improve the simulation of the carry over effects and predictions of EONR, especially in the SC rotation.
There is a need for further research to improve the annual EONR predictability in extreme weather years with precipitation above or below historical average precipitations. Differences between predicted and observed EONR tended to be higher for cold weather years (temperature below the historical average), and their agreement was largely affected by the inaccurate prediction on extremely wet (i.e., 2008) and dry (i.e., 2000) weather years (Figure S2). In this study, we noticed that EONR forecast was more accurate in years that had EONR values near the long-term site average and the prediction accuracy decreased in years with extreme high or low EONR measured values (±20% from long-term site average, Table 2).
For the farmer, both extreme high or low EONR predictions result in economic loss, with the greatest being yield loss when under-fertilizing in high N rate requirement years. Furthermore, when high N is predicted in low N need years the risk of N loss to the environment is greatest (Raun et al., 2017). In particular, the model failed to forecast accurate EONR values in extreme high precipitation years such as 2008 and 2010. We believe the main reason is the incorrect simulation of the carry-over effects from previous year (i.e., soil inorganic nitrogen, soil moisture, root and carbon and nitrogen inputs from previous crops) or possible underestimation of N losses from the system (He et al., 2017; Yin et al., 2017). Prediction of EONR is very complex and less accurate than yield (Figure 2). Perhaps an alternative way of using the model as a forecast tool may be to predict the key-components of the EONR such as yields and soil supply and plant N uptake dynamics, in which the model performs well (Archontoulis and Licht, 2017). In that way, we avoid accumulation of errors that leads to a low prediction accuracy of the EONR.
For weather uncertainty affecting yield and derived EONR forecasts (Hansen et al., 2006), our results indicated that use of the entire historical record (35 years) in the simulation process is better than selecting years based on similar weather events or using only a smaller set number of previous years (Figure 7). This was more evident in EONR predictions and less in yield predictions. Although the reasons are not sufficiently clear, based on other studies, we believe that simulated crop growth depends more on the distribution of weather within the season than the season average which was used to categorize the years used for the weather scenarios (Hansen and Indeje, 2004; Figure S1b).
We also found that the use of at least 10 and 20 years of weather data is associated with small error (2.5% less accuracy than the 35-year; Figure 7) which is an important finding given that not all sites have accurate weather records for 35 years (Hansen et al., 2004; Grassini et al., 2015). Our findings agree with Van Wart et al. (2013), who showed that 6–10 years of weather data is needed for sites with similar rain patterns as the one used in this study. However, we also recognize there are other ways of selecting years to be used in forecasting studies and this is a topic for further investigation (Hansen et al., 2004).
Conclusions
This study provided evidence that use of a calibrated cropping systems model can aid yield and N rate forecasts. At planting time, model predictions were directionally correct in predicting whether the optimum N-rate for corn would be above, below or at site-mean value. The associated prediction error was within an error range of ±30 kg N ha−1 in ~60% of the years. Predictability of corn yields was more accurate than optimum N-rates at planting time. In most years, in-season yield and optimum N-rate forecasts were not better than the predictions made at the beginning of the season. Use of 35-year historical weather data was found to be more accurate than using analogous weather years based on similarity of current weather conditions. Process-based modeling forecast of corn yields and optimum N-rates is already promising and has potential for realizing further improvements in achieving even more accurate early-season predictions.
Author Contributions
LP, JS, and SA: conceived the research idea and designed the modeling study; LP and SA: performed the model analysis and drafted the manuscript; All coauthors (JS, MC, PT, KM, EH, AV) reviewed the manuscript and contributed to the final version.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work was part of the Agriculture and Food Research Initiative Hatch project No. 1004346 and partially funded by the Iowa State University College of Agriculture and Life Sciences and Plant Science Institute. We also thank the Agronomy Research Farm manager and personnel for their assistance with the research and data collection over the years.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00436/full#supplementary-material
References
Abendroth, L. J., Elmore, R. W., Boyer, M. J., and Marlay, S. K. (2011). Corn growth and development. Publ. PMR 1009. Iowa State University Extension and Outreach, Ames, IA.
Amatya, P., Yu, M., and Ewell, F. (2008). Economic analysis of optimal nitrogen application in corn production. Texas J. Agric. Nat. Res. 21, 101–108.
Andrade, F. H., Vega, C., Uhart, S., Cirilo, A., Cantarero, M., and Valentinuz, O. (1999). Kernel number determination in maize. Crop Sci. 39, 453–459. doi: 10.2135/cropsci1999.0011183X0039000200026x
Arbuckle, J., and Gordon Lasley, P. (2013). Iowa Farm and Rural Life Poll: 2013 Summary Report. Extension Community and Economic Development Publications. 26. Available online at: https://lib.dr.iastate.edu/extension_communities_pubs/26
Archontoulis, S., and Licht, M. (2017). Forecasting and Assessment of Cropping Systems in Northwest Iowa. Farm Progress Reports. Available online at: https://lib.dr.iastate.edu/farmprogressreports/vol2016/iss1/118
Archontoulis, S., Licht, M., Dietzel, R., Catellano, M., VanLoocke, A., Moore, K., et al. (2016). “In-Season forecasting of crop yields, soil water-nitrogen, and weather using APSIM and WRF models in Iowa, USA,” in International Crop Modelling Symposium, 15-17 March. Berlin, 188–189.
Archontoulis, S. V., and Miguez, F. E. (2015). Nonlinear regression models and applications in agricultural research. Agron. J. 107, 786–798. doi: 10.2134/agronj2012.0506
Banger, K., Yuan, M., Wang, J., Nafziger, E. D., and Pittelkow, C. M. (2017). A Vision for Incorporating environmental effects into nitrogen management decision support tools for US maize production. Front. Plant Sci. 8:1270. doi: 10.3389/fpls.2017.01270
Basso, B., and Ritchie, J. T. (2015). “Simulating crop growth and biogeochemical fluxes in response to land management using the SALUS model,” in The Ecology of Agricultural Landscapes: Long-Term Research on The Path to Sustainability, eds S. K. Hamilton, J. E. Doll, and G. P. Robertson (New York, NY: Oxford University Press), 252–274.
Basso, B., Ritchie, J. T., Pierce, F. J., Braga, R. P., and Jones, J. W. (2001). Spatial validation of crop models for precision agriculture. Agric. Syst. 68, 97–112. doi: 10.1016/S0308-521X(00)00063-9
Binford, G. D., Blackmer, A. M., and Meese, B. G. (1992). Optimal concentrations of nitrate in cornstalks at maturity. Agron. J. 84:881–887. doi: 10.2134/agronj1992.00021962008400050022x
Blackmer, A. M., Pottker, D., Cerrato, M. E., and Webb, J. (1989). Correlations between soil nitrate concentrations in late spring and corn fields in Iowa. J. Prod. Agric. 2, 103–109. doi: 10.2134/jpa1989.0103
Borrás, L., and Vitantonio-Mazzini, L. N. (2018). Maize reproductive development and kernel set under limited plant growth environments. J. Exp. Bot. erx452. doi: 10.1093/jxb/erx452
Bullock, D. G., and Bullock, D. S. (1994). Quadratic and quadratic plus plateau models for predicting optimal nitrogen rate of corn: a comparison. Agron. J. 86, 191–195. doi: 10.2134/agronj1994.00021962008600010033x
Bundy, L. G., and Andraski, T. W. (1995). Soil potential yield effects on performance of soil nitrate tests. J. Prod. Agric. 8, 561–568.
Campbell, G. S. (1985). Soil Physics With BASIC. Transport Models for Soil-Plant Systems. New York, NY: Elsevier Science Publishers.
Carberry, P. S., Hochman, Z., Hunt, J. R., Dalgliesh, N. P., McCown, R. L., Whish, J. P. M., et al. (2009). Re-inventing model-based decision support with Australian dryland farmers. 3. Relevance of APSIM to commercial crops. Crop Past. Sci. 60, 1044–1056. doi: 10.1071/CP09052
Cerrato, M. E., and Blackmer, A. M. (1990). Comparison of models for describing corn yield response to nitrogen fertilizer. Agron. J. 82, 38–143. doi: 10.2134/agronj1990.00021962008200010030x
Constantin, J., Beaudoin, N., Laurent, F., Cohan, J. P., Duyme, F., and Mary, B. (2011). Cumulative effects of catch crops on nitrogen uptake, leaching and net mineralization. Plant Soil 341, 137–154. doi: 10.1007/s11104-010-0630-9
Coulter, J. A., and Nafziger, E. D. (2008). Continuous corn response to residue management and nitrogen fertilization. Agron. J. 100, 1774–1780. doi: 10.2134/agronj2008.0170
Dhital, S., and Raun, W. R. (2016). Variability in optimum nitrogen rates for maize. Agron. J. 108, 2165–2173. doi: 10.2134/agronj2016.03.0139
Franzen, D., Kitchen, N., Holland, K., Schepers, J., and Raun, W. (2016). Algorithms for in-season nutrient management in cereals. Agron. J. 108, 1775–1781. doi: 10.2134/agronj2016.01.0041
Gowda, P. H., Mulla, D. J., and Jaynes, D. B. (2008). Simulated long-term nitrogen losses for a midwestern agricultural watershed in the United States. Agric. Water Manage. 95, 616–624. doi: 10.1016/j.agwat.2008.01.004
Grassini, P., van Bussel, L. G., Van Wart, J., Wolf, J., Claessens, L., Yang, H., et al. (2015). How good is good enough? Data requirements for reliable crop yield simulations and yield-gap analysis. Field Crops Res. 177, 49–63. doi: 10.1016/j.fcr.2015.03.004
Hammer, G. L., Hansen, J. W., Phillips, J., Mjelde, J. W., Hill, H. S. J., Love, A., et al. (2001). Advances in application of climate prediction in agriculture. Agric. Syst. 70, 515–553. doi: 10.1016/S0308-521X(01)00058-0
Hammer, G. L., Holzworth, D. P., and Stone, R. (1996). The value of skill in seasonal climate forecasting to wheat crop management in a region with high climatic variability. Crop Past. Sci. 47, 717–737. doi: 10.1071/AR9960717
Hansen, J. W., Challinor, A., Ines, A., Wheeler, T., and Moron, V. (2006). Translating climate forecasts into agricultural terms: advances and challenges. Climate Res. 33, 27–41. doi: 10.3354/cr033027
Hansen, J. W., and Indeje, M. (2004). Linking dynamic seasonal climate forecast with crop simulation for maizeyield prediction in semi-arid Kenya. Agricult. Forest Meteorol. 125, 143–157. doi: 10.1016/j.agrformet.2004.02.006
Hansen, J. W., Potgieter, A., and Tippett, M. K. (2004). Using a general circulation model to forecast regional wheat yields in northeast Australia. Agric. For. Meteorol. 127, 77–92. doi: 10.1016/j.agrformet.2004.07.005
Hawkins, J. A., Sawyer, J. E., Barker, D. W., and Lundvall, J. P. (2007). Using relative chlorophyll meter values to determine nitrogen application rates for corn. Agron. J. 99, 1034–1040. doi: 10.2134/agronj2006.0309
He, D., Wang, E., Wang, J., and Robertson, M. J. (2017). Data requirement for effective calibration of process-based crop models. Agric. For. Meteorol. 234, 136–148. doi: 10.1016/j.agrformet.2016.12.015
Holzworth, D. P., Huth, N. I., Zurcher, E. J., Herrmann, N. I., McLean, G., Chenu, K., et al. (2014). APSIM–evolution towards a new generation of agricultural systems simulation. Environ. Model. Softw. 62, 327–350. doi: 10.1016/j.envsoft.2014.07.009
Huth, N. I., Bristow, K. L., and Verburg, K. (2012). SWIM3: model use, calibration, and validation. Trans. Am. Soc. Agric. Eng. 55, 1303–1313. doi: 10.13031/2013.42243
Ines, A. V., Das, N. N., Hansen, J. W., and Njoku, E. G. (2013). Assimilation of remotely sensed soil moisture and vegetation with a crop simulation model for maize yield prediction. Remote Sens. Environ. 138, 149–164. doi: 10.1016/j.rse.2013.07.018
Iqbal, J., Necpalova, M., Archontoulis, S. V., Anex, R. P., Bourguignon, M., Herzmann, D., et al. (2017). Extreme weather-year sequences have non-additive effects on environmental nitrogen losses. Glob. Change Biol. 24, e303–e317. doi: 10.1111/gcb.13866
Jemison, J. M., and Lytle, D. E. (1996). Field evaluation of two nitrogen testing methods in Maine. J. Prod. Agric. 9, 108–113. doi: 10.2134/jpa1996.0108
Jin, Z., Prasad, R., Shriver, J., and Zhuang, Q. (2017). Crop model-and satellite imagery-based recommendation tool for variable rate N fertilizer application for the US Corn system. Precision Agric. 18, 779–800. doi: 10.1007/s11119-016-9488-z
Keating, B. A., Carberry, P. S., Hammer, G. L., Probert, M. E., Robertson, M. J., Holzworth, D., et al. (2003). An overview of APSIM, a model designed for farming systems simulation. Eur. J. Agron. 18, 267–288. doi: 10.1016/S1161-0301(02)00108-9
Kersebaum, K. C., Lorenz, K., Reuter, H., Schwarz, J., Wegehenkel, M., and Wendroth, O. (2005). Operational use of agro-meteorological data and GIS to derive site-specific nitrogen fertilizer recommendations based on the simulation of soil and crop growth processes. Phys. Chem. Earth 30, 59–67. doi: 10.1016/j.pce.2004.08.021
Kwaw-Mensah, D., and Al-Kaisi, M. (2006). Tillage and nitrogen source and rate eff ects on corn response in corn-soybean rotation. Agron. J. 98, 507–513. doi: 10.2134/agronj2005.0177
Liu, S., Yang, J. Y., Zhang, X. Y., Drury, C. F., and Reynolds, H. G. (2013). Modelling crop yield, soil water content and soil temperature for asoybean– maize rotation under conventional and conservation tillage systems in Northeast China. Agric. Water Manage. 123, 32–44. doi: 10.1016/j.agwat.2013.03.001
Lory, J. A., and Scharf, P. C. (2003). Yield goal versus delta yield for predicting fertilizer nitrogen need in corn. Agron. J. 95, 994–999. doi: 10.2134/agronj2003.0994
Ma, L., Malone, R. W., Heilman, P., Karlen, D. L., Kanswar, R. S., Cambardella, C. A., et al. (2007). RZWQM Simulation of long term crop production, water and nitrogen balances in Northeast Iowa. Geoderma 140, 247–259. doi: 10.1016/j.geoderma.2007.04.009
Magdoff, F. R., Ross, D., and Amadon, J. (1984). A soil test for nitrogen availability to corn. Soil Sci. Soc. Am. J. 48, 1301–1304. doi: 10.2136/sssaj1984.03615995004800060020x
Mamo, M., Malzer, G. L., Mulla, D. J., Huggins, D. R., and Strock, J. (2003). Spatial and temporal variation in economically optimum nitrogen rate for corn. Agron. J. 95, 958–964. doi: 10.2134/agronj2003.09580
Morell, F. J., Yang, H. S., Cassman, K. G., Wart, J. V., Elmore, R. W., Licht, M., et al. (2016). Can crop simulation models be used to predict local to regional maize yields and total production in the U.S. Corn Belt? Field Crops Res. 192, 1–12. doi: 10.1016/j.fcr.2016.04.004
Morris, T. F., Murrell, T. S., Beegle, D. B., Camberato, J. J., Ferguson, R. B., Grove, J., et al. (2018). Strengths and Limitations of nitrogen rate recommendations for corn and opportunities for improvement. Agron. J. 110, 1–37. doi: 10.2134/agronj2017.02.0112
Nangia, V., Turral, H., and Molden, D. (2008). Increasing water productivity with improved N fertilizer management. Irrig. Drainage Syst. 22, 193–207. doi: 10.1007/s10795-008-9051-9
Nendel, C., Wieland, R., Mirschel, W., Specka, X., Guddat, C., and Kersebaum, K. C. (2013). Simulating regional winter wheat yields using input data of different spatial resolution. Field Crop. Res. 145, 67–77. doi: 10.1016/j.fcr.2013.02.014
Probert, M. E., Dimes, J. P., Keating, B. A., Dalal, R. C., and Strong, W. M. (1998). APSIM's water and nitrogen modules and simulation of the dynamics of water and nitrogen in fallow systems. Agric. Syst. 56, 1–28.
Puntel, L. A., Sawyer, J. E., Barker, D. W., Dietzel, R., Poffenbarger, H., and Castellano, M. J. (2016). Modeling long-term corn yield response to nitrogen rate and crop rotation. Front Plant Sci. 7:630. doi: 10.3389/fpls.2016.01630
Quiring, S. M., and Legates, D. R. (2008). Application of CERES-Maize for within-season prediction of rainfed corn yields in Delaware, USA. Agric. For. Meteorol. 148, 964–975. doi: 10.1016/j.agrformet.2008.01.009
Raun, W., Figueiredo, B., Dhillon, J., Fornah, A., Bushong, J., Zhang, H., et al. (2017). Can yield goals be predicted? Agron. J. 109, 2389–2395. doi: 10.2134/agronj2017.05.0279
Reimer, A. J. E., Doll, B., Basso, S. T., Marquart-Pyatt, G. P., Robertson, D. Zhao, et al. (2017). Moving toward sustainable farming systems: Insights from private and public sector dialogues on nitrogen management. J. Soil Water Conserv. 72, 5A−9A. doi: 10.2489/jswc.72.1.5A
Salo, T. J., Palosuo, T., Kersebaum, K. C., Nendel, C., Angulo, C., Ewert, F., et al. (2016). Comparing the performance of 11 crop simulation models in predicting yield response to nitrogen fertilization. J. Agric. Sci. 154, 1218–1240. doi: 10.1017/S0021859615001124
Sawyer, J. E. (2015). 45th Annual North Central Extension-Industry Soil Fertility Conference. Integrated Crop Management News. Available online at: https://lib.dr.iastate.edu/cropnews/1488
Sawyer, J., Nafziger, E., Randall, G., Bundy, L., Rehm, G., and Joern, B. (2006). Concepts and Rationale for Regional Nitrogen Guidelines for Corn. Iowa State Univ. Extension Publ. PM2015. Available online at: http://www.extension.iastate.edu/Publications/PM2015.pdf
Scharf, P. C. (2001). Soil and plant tests to predict optimum nitrogen rates for corn. J. Plant Nutr. 24, 805–826. doi: 10.1081/PLN-100103775
Scharf, P. C. (2015). “Managing nitrogen,” in Managing Nitrogen in Crop Production (Madison, WI, USA: American Society of Agronomy, Inc., Crop Science Society of America, Inc., and Soil Science Society of America, Inc.), 25–76.
Scharf, P. C., Kitchen, N. R., Sudduth, K. A., Davis, J. G., Hubbard, V. C., and Lory, J. A. (2005). Field-scale variability in optimal nitrogen fertilizer rate for corn. Agron. J. 97, 452–461. doi: 10.2134/agronj2005.0452
Scharf, P. C., Kitchen, N. R., Sudduth, K. A., and Davis, J. G. (2006). Spatially variable corn yield is a weak predictor of optimal nitrogen rate. Soil Sci. Soc. Am. J. 70, 2154–2160. doi: 10.2136/sssaj2005.0244
Scharf, P. C., Scmidt, J. P., Kitchen, N. R., Sudduth, K. A., Young, S. Y., Lory, J. A., et al. (2002). Remote sensing for nitrogen management. J. Soil Water Cons. 57, 518–524.
Sela, S., van Es, H. M., Moebius-Clune, B. N., Marjerison, R., Moebius-Clune, D., Schindelbeck, R., et al. (2017). Dynamic model improves agronomic and environmental outcomes for maize nitrogen management over static approach. J. Environ. Qual. 46, 311–319. doi: 10.2134/jeq2016.05.0182
Semenov, M. A., Porter, J. R., and Dele' cole, R. (1993). Climatic change and the growth and development of wheat in the UK and France. Eur. J. Agron. 2, 293–304 doi: 10.1016/S1161-0301(14)80177-9
Shanahan, J. F., Kitchen, N. R., Raun, W. R., and Schepers, J. S. (2008). Responsive in-season nitrogen management for cereals. Comput. Electron. Agric. 61, 51–62. doi: 10.1016/j.compag.2007.06.006
Shapiro, C. A., Ferguson, R. B., Hergert, G. W., Wortman, C. S., and Walters, D. T. (2008). Fertilizer Suggestions for Corn. Lincoln, NE: University of Nebraska.
Soufizadeh, S., Munaro, E., McLean, G., Massignam, A., van Oosterom, E. J., Chapman, S. C., et al. (2018). Modelling the nitrogen dynamics of maize crops – Enhancing the APSIM maize model. Eur. J. Agric. doi: 10.1016/j.eja.2017.12.007
Stanford, G. (1973). Rationale for optimum nitrogen fertilization in corn production. J. Environ. Qual. 2, 159–166. doi: 10.2134/jeq1973.00472425000200020001x
Stanford, G., and Legg, J. O. (1984). “Nitrogen and yield potential,” in Nitrogen in Crop Production ed R. D. Hauck (Madison, WI: ASA, CSSA, and SSSA), 263–272.
Thompson, L. J., Ferguson, R. B., Kitchen, N., Frazen, D. W., Mamo, M., Yang, H., et al. (2015). Model and sensor-based recommendation approaches for in-season nitrogen management in corn. Agron. J. 107, 2020–2030. doi: 10.2134/agronj15.0116
Thornton, P. K., Bowen, W. T., Ravelo, A. C., Wilkens, P. W., Farmer, G., Brock, J., et al. (1997). Estimating millet production for famine early warning: an application of crop simulation modeling using satellite and ground based data in Burkina Faso. Agric. For. Meteorol. 83, 95–112.
Thorburn, P. J., Meier, E. A., and Probert, M. E. (2005). Modelling nitrogen dynamics in sugarcane systems: recent advances and applications. Field Crop Res. 92, 337–352. doi: 10.1016/j.fcr.2005.01.016
Thorburn, P. J., Probert, M. E., and Robertson, F. A. (2001). Modelling decomposition of sugarcane surface residues with APSIM-Residue. Field Crop Res. 70, 223–232. doi: 10.1016/S0378-4290(01)00141-1
Thorp, K. R., Malone, R. W., and Jaynes, D. B. (2007). Simulating long-term effects of nitrogen fertilizer application rates on corn yield and nitrogen dynamics. Trans. ASABE 50, 1287–1303. doi: 10.13031/2013.23640
Togliatti, K., Archontoulis, V. S., Dietzel, R., Puntel, A., and Laila VanLoocke, A. (2017). How does inclusion of weather forecasting impact in-season crop model predictions? Field Crop Res. 214, 261–272. doi: 10.1016/j.fcr.2017.09.008
Tremblay, N., Bouroubi, Y., Bélec, C., Mullen, R., Kitchen, N., and Thomason, W. (2012). Corn response to nitrogen is influenced by soil texture and weather. Agron J. 104, 1658–1671. doi: 10.2134/agronj2012.0184
Van Es, H. M., Kay, B. D., Melkonian, J. J., Sogbedji, J. M., and Bruulsma, T. W. (2006). “Nitrogen management for maize in humid regions: Case for a dynamic modeling approach,” in Managing Crop Nitrogen for Weather: Proceedings of the Symposium “Integrating Weather Variability into Nitrogen Recommendations,” (Indianapolis, IN).
Van Wart, J., Kersebaum, K. C., Peng, S., Milner, M., and Cassman, K. G. (2013). Estimating crop yield potential at regional to national scales. Field Crops Res. 143, 34–43. doi: 10.1016/j.fcr.2012.11.018
Voss, R. D., and Shrader, W. D. (1979). Crop rotations: Effect on Yields and Response to Nitrogen. Cooperative Extension Service, Iowa State University.
Yang, J. M., Yang, J. Y., Liu, S., and Hoogenboom, G. (2014). An evaluation of the statistical methods for testing the performance of crop models with observed data. Agric. Syst. 127, 81–89. doi: 10.1016/j.agsy.2014.01.008
Yin, X., Kersebaum, K. C., Kollas, C., Manevski, K., Baby, S., Beaudoin, N., et al. (2017). Performance of process-based models for simulation of grain N in crop rotations across Europe. Agric. Sys. 154, 63–77. doi: 10.1016/j.agsy.2017.03.005
Keywords: corn, economic optimum N rate, forecast, modeling, APSIM, in-season nitrogen management, nutrient recommendation
Citation: Puntel LA, Sawyer JE, Barker DW, Thorburn PJ, Castellano MJ, Moore KJ, VanLoocke A, Heaton EA and Archontoulis SV (2018) A Systems Modeling Approach to Forecast Corn Economic Optimum Nitrogen Rate. Front. Plant Sci. 9:436. doi: 10.3389/fpls.2018.00436
Received: 28 September 2017; Accepted: 21 March 2018;
Published: 13 April 2018.
Edited by:
Urs Feller, Universität Bern, SwitzerlandReviewed by:
Jose M. Garcia-Mina, Universidad de Navarra, SpainIgnacio Antonio Ciampitti, Kansas State University, United States
Tony John Vyn, Purdue University, United States
Copyright © 2018 Puntel, Sawyer, Barker, Thorburn, Castellano, Moore, VanLoocke, Heaton and Archontoulis. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Laila A. Puntel, lapuntel@iastate.edu
Sotirios V. Archontoulis, sarchont@iastate.edu