Original Research ARTICLE
Modeling Long-Term Corn Yield Response to Nitrogen Rate and Crop Rotation
- 1Department of Agronomy, Iowa State University, Ames, IA, USA
- 2Commonwealth Scientific and Industrial Research Organisation Agriculture, St Lucia, QLD, Australia
Improved prediction of optimal N fertilizer rates for corn (Zea mays L.) can reduce N losses and increase profits. We tested the ability of the Agricultural Production Systems sIMulator (APSIM) to simulate corn and soybean (Glycine max L.) yields, the economic optimum N rate (EONR) using a 16-year field-experiment dataset from central Iowa, USA that included two crop sequences (continuous corn and soybean-corn) and five N fertilizer rates (0, 67, 134, 201, and 268 kg N ha-1) applied to corn. Our objectives were to: (a) quantify model prediction accuracy before and after calibration, and report calibration steps; (b) compare crop model-based techniques in estimating optimal N rate for corn; and (c) utilize the calibrated model to explain factors causing year to year variability in yield and optimal N. Results indicated that the model simulated well long-term crop yields response to N (relative root mean square error, RRMSE of 19.6% before and 12.3% after calibration), which provided strong evidence that important soil and crop processes were accounted for in the model. The prediction of EONR was more complex and had greater uncertainty than the prediction of crop yield (RRMSE of 44.5% before and 36.6% after calibration). For long-term site mean EONR predictions, both calibrated and uncalibrated versions can be used as the 16-year mean differences in EONR’s were within the historical N rate error range (40–50 kg N ha-1). However, for accurate year-by-year simulation of EONR the calibrated version should be used. Model analysis revealed that higher EONR values in years with above normal spring precipitation were caused by an exponential increase in N loss (denitrification and leaching) with precipitation. We concluded that long-term experimental data were valuable in testing and refining APSIM predictions. The model can be used as a tool to assist N management guidelines in the US Midwest and we identified five avenues on how the model can add value toward agronomic, economic, and environmental sustainability.
The economic optimum nitrogen (N) rate (EONR) is the fertilizer rate at which crop yield increase is not large enough to pay for additional N application, and therefore more N would only result in unnecessary costs (Sawyer et al., 2006). Optimal N input needs to be considered when making N recommendations since it has the potential to improve N use efficiency, crop yield, and profitability as well as to reduce environmental impacts (Wang et al., 2003; Lawlor et al., 2008; Kyveryga et al., 2009; Basso et al., 2016). Nitrogen losses by leaching are proportional to the N rate applied and tend to increase rapidly at rates greater than optimal for crop use (Haghiri et al., 1978; Cooper and Cooke, 1984; Andraski et al., 2000; Randall et al., 2000).
There is tremendous uncertainty and risk associated with prediction of the EONR in corn–based systems, both at the field and sub-field scale (Paz et al., 1999; Scharf et al., 2005; Tremblay et al., 2012). Farmers may attempt to protect corn yield potential with high fertilizer N inputs, which leads to decreased profitability (Lambert et al., 2006) and increased likelihood of environmental contamination (Andraski et al., 2000; Jaynes et al., 2001; Robertson and Groffman, 2007).
A number of approaches have been developed to predict optimal N application rates. These include yield goal-based N recommendations and N budgets (Stanford, 1973, 1982; Stanford and Legg, 1984), pre-plant and pre-sidedress soil nitrate test (PPNT and PSNT, Bundy and Andraski, 1995; Shapiro et al., 2008), Illinois soil nitrogen test (ISNT, Mulvaney et al., 2001), crop canopy sensing (NDVI, Schmidt et al., 2009 and chlorophyll meter, Blackmer and Schepers, 1995; Varvel et al., 1997), and economic maximum return to N (MRTN, Sawyer et al., 2006). Some of these tools are static in that they give the same recommendation regardless of yearly weather or crop/fertilizer prices, or evaluate N status after grain crop harvest. Soil tests or hand-held crop meters are often time consuming, expensive, and/or require periodic and intense sampling (Blackmer et al., 1997; Ma and Dwyer, 1999; Grove and Schwab, 2006; van Es et al., 2007; Lemaire et al., 2008; Franzen et al., 2016). Most current and widely adopted N management practices also assume field-uniformity, recommending N applications that ignore variation in landscape factors such as topography, soil texture, and organic matter (Cassman et al., 2002; Mamo et al., 2003; Scharf et al., 2005), as well as interactions with plant population and hybrid (Ciampitti and Vyn, 2012). Use of precision agriculture technologies (real-time remote sensing, unmanned aerial images, soil mapping, etc.) combined with variable N application have the potential to increase N use efficiency by matching the N requirements within field zones (Dobermann and Cassman, 2002; Ferguson et al., 2002; Mamo et al., 2003; Mulla, 2013). However, the selection of a site-specific optimum N rate is difficult to predict based on the large temporal and spatial variability of the N supply and demand (van Es et al., 2007; Setiyono et al., 2011). Unfortunately, the above approaches have not fully resolved needed improvements from N management and gains in N use efficiency (Raun and Johnson, 1999; Fageria and Baligar, 2005) since N losses from corn-based systems are still high with negative environment impacts (Jaynes et al., 2001; Mitsch et al., 2001).
The challenge in managing N and estimating the optimum N fertilization rate comes from the complex interactions that exist in the dynamic soil-plant-atmosphere system and uncertainty in weather (Havlin et al., 2005; Tremblay and Belec, 2006; Brady and Weil, 2008). Soil N mineralization from SOC and crop N uptake, and N losses are three important components defining the optimum N rate, however, these processes are dynamic and difficult to predict (Cassman et al., 2002). Therefore N management tools that simultaneously consider dynamics in soil organic carbon mineralization, crop growth, weather conditions, and agronomic practices may greatly improve site- and year-specific EONR estimates (Basso et al., 2012, 2016; Dumont et al., 2016). Dynamic cropping system simulation models such as Agricultural Production Systems sIMulator (APSIM; Holzworth et al., 2014), DSSAT (Jones et al., 2003), RZWQM (Ahuja et al., 2000), CropSyst (Stockle et al., 2003), SALUS (Basso et al., 2006), and others have been used to investigate soil-crop-weather dynamics, however, model use has been limited to address long-term optimum N rates (Ma et al., 2007; Basso et al., 2010). The scientific literature is also rich with examples of model applications to improve our understanding of N dynamics and to answer questions that cannot be addressed with field research due to time and cost constraints (Batchelor et al., 2002; Schnebelen et al., 2004; Fountas et al., 2006; Malone et al., 2010; Basso et al., 2012, 2016; Anapalli et al., 2014). However, use of models in practical applications to assist real-life challenges such as N rate guidance is limited because models typically require: (a) a large number of input parameters, which are usually not available (Wallach, 2006; Basso et al., 2012); (b) particular skills to develop model specific input parameters and cultivar coefficients from internet databases; and (c) intensive training for use.
Over the last few years web-applications have been developed to simplify the use of models (e.g., Yield Prophet, Carberry et al., 2009). Furthermore, digital soil and weather databases such as web soil survey1 (Soil Survey Staff, 2006) and daymet (Daymet, 1980–2008; Thornton et al., 2012) provide free access to high-resolution input parameters. As a result, the potential of using simulation models to assist with real-life practical problems and especially to predict the risk associated with selecting specific N fertilizer rates has received strong industrial interest (Thorp et al., 2007; Gowda et al., 2008; Nangia et al., 2008). The next challenge to applying models across different scales (within fields, regions, and cropping systems) is to determine prediction accuracy; e.g., how well cropping system models can predict crop yield, N dynamics, and EONR. And if they can predict corn response to N rate, how can this information be used to develop better N rate guidelines.
In this study we used a 16-year field research dataset from a site in central Iowa, USA that included five N rates and two crop sequences to test the ability of the APSIM model (Holzworth et al., 2014) to predict crop yields and optimal N rate for corn. Our specific objectives were to: (a) quantify model prediction accuracy before and after calibration, and report calibration steps; (b) compare crop model-based techniques in estimating optimal N rate for corn; and (c) utilize the calibrated model to explain factors causing year to year variability in yield and optimal N. The APSIM model was selected for use in this study because of its flexibility and easy use in specifying crop rotations via the user interface, capability in simulating long-term dynamics in both soil and crop processes, advanced flexibility in simulating the effect of shallow water table dynamics that are important in this geographic region (Helmers et al., 2012) and previously determined good performance in this geographic region (Malone et al., 2007; Hammer et al., 2009; Lobell et al., 2013; Archontoulis et al., 2014a,b, 2016; Basche et al., 2016; Dietzel et al., 2016; Martinez-Feria et al., 2016).
Materials and Methods
Site, Weather, and Experimental Datasets
The field-experiment was conducted at the Agricultural Engineering and Agronomy Research Farm near Ames, Iowa, USA (42° 0′37.50″N, 93°47′22.98″W) on a Clarion loam soil (fine-loamy, mixed, superactive, mesic Typic Hapludoll). The experiment was initiated in 1999 and continuous to the present. For this study we used data from 1999 to 2014 (16-years). The climate at the site is humid continental (warm, rainy summers) with annual precipitation of 900 mm and a mean temperature of 9°C (Supplementary Figure S1). 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; Supplementary Figure S1). Years 2008, 2010, and 2014 were the wettest and years 2000, 2011, 2012, and 2013 the driest. Mean annual air temperatures were 16 and 23°C for spring and summer, respectively. Year 2012 was the warmest and year 2008 the coolest (Supplementary Figure S1).
The long-term experiment was designed to study the effect of five N fertilizer rates (0, 67, 134, 201, and 268 kg N ha-1; hereafter N0, N67, N134, N201, and N268, respectively) on corn yield in continuous corn (CC) and soybean-corn rotation (SC). The experimental design was a randomized complete block design with four replications. Nitrogen fertilizer was applied near planting (± 10–15 days). Specific information on the fertilizer type and application dates are provided in Supplementary Table S1. Within the SC rotation, corn and soybean phases were present each year in the rotation: thus a simulation starting with corn in year one and another simulation starting with soybeans on year one were set up. Hereafter SC when the rotation starts with corn in year one (odd numbered years) and a validation set (SC_val) when the rotation starts with soybean in year 1 (even numbered years). Each treatment had four replications. Nitrogen fertilizer was only applied to corn. Supplementary Table S1 provides management information by year and rotation. Measurements included corn and soybean grain yields each year (expressed at 15.5 and 13% moisture content, respectively). Soil organic carbon measurements were available at 0–15 cm in 1999, 2009, and 2014, and at 0–30 cm in 2009 for CC (Brown et al., 2014; Poffenbarger et al., unpublished).
The APSIM Modeling Platform
The APSIM (Keating et al., 2003; Holzworth et al., 2014) is an open-source advanced simulator of agricultural systems that combines several process-based models in a modular design. APSIM is a field-scale model that operates mainly on a daily time step. Details about APSIM and its performance across a range of studies can be found at http://www.apsim.info.
APSIM Configuration and Calibration
Two rounds of APSIM model evaluations were performed; a blind phase (uncalibrated model) where management and cultivar information were used, and a calibrated phase (calibrated model) where crop yield and SOC data were provided into the model. Similar protocols have been used in the AgMIP project (Agricultural Model Inter-Comparison and Improvement Project; Rosenzweig et al., 2013).
Blind-Phase Model Parameters and Set-Up
For the blind phase, we first incorporated available management information into APSIM (Supplementary Table S1). When required management information was unavailable, we used typical values from the literature relevant to the research site (Abendroth et al., 2011; Pedersen and Licht, 2014). The following input parameters were held constant across the 16-years: planting depth of 5 cm for both crops, plant populations of 8 and 38 plants m-2 for corn and soybean, respectively, and November 10th and April 10th dates for fall and spring tillage operations; and corn hybrid (106-day) and soybean variety (2.5 maturity group) values derived from previous studies in the region (Archontoulis et al., 2014a,b; Supplementary Table S2). Daily weather data were obtained from the Iowa Environmental Mesonet (2014). Soil profile information was taken from Web Soil Survey (Soil Survey Staff, 2006) and soil-root related parameters were developed following the methodology described in Archontoulis et al. (2014a). The maximum rooting depths for corn and soybean were set to 1.5 and 1.2 m, respectively.
We set up APSIM by connecting the following models: corn and soybean crop models (Keating et al., 2003), Soil N (soil N and C cycling model with default soil temperature model; Probert et al., 1998), SoilWat (a tipping bucket soil water model; Probert et al., 1998); SURFACEOM (residue model; Probert et al., 1998; Thorburn et al., 2001, 2005), and the following management rules: planting, harvesting, fertilizer, tillage, and rotations (Keating et al., 2003). In addition we implemented within the MANAGER module an N deposition rule that simulates atmospheric N deposition as a function of daily precipitation (N deposition in kg N ha-1 d-1 = 0.01∗ precipitation in mm; Holland et al., 2005). On average this added about 7 kg N ha-1 year-1 into the system. Initial model conditions such as root mass, surface residue mass, soil water, soil nitrate, and SOC pool partitioning were obtained by starting the model 6 years prior to the start of the experiment (Supplementary Table S3). Experience using APSIM in this geographic region for simulating corn-soybean production systems has indicated that the fast microbial SOC pool (BIOM) of APSIM requires at least 4 years to stabilize (Basche et al., 2016; Dietzel et al., 2016; Martinez-Feria et al., 2016). Having this pool stabilized is important to remove confounding effects of microbial SOC buildup or decline which affects N dynamics. The APSIM version 7.6 was used on a daily time step. The simulation process was consecutive to account for carry-over effects from year to year, such as soil inorganic nitrogen, soil moisture, root and residue carbon and nitrogen inputs from previous crops.
Model Calibration and Testing
In the calibration and testing phase, we used end-of-season grain yields and SOC data to improve predictions. The long-term (end-of-season) data are powerful in detecting weakness in the model (i.e., years with low prediction accuracy), but do not provide guidance on which of the model’s processes or parameters needed to be improved. Therefore, to inform the calibration process, additional information was used: knowledge gained from other APSIM calibration studies in Iowa (Archontoulis et al., 2014a, 2016; Basche et al., 2016; Dietzel et al., 2016), sensitivity techniques, and model behavior analysis coupled with expert judgment (Supplementary Figure S2). The odd-numbered years (for CC and SC) were used for calibration and even numbered years for validation (SC_val dataset).
During calibration the following changes in APSIM were made: first, we replaced the default APSIM soil temperature model that uses EPIC model equations (Williams et al., 1984) with a more mechanistic soil temperature model (Campbell, 1985) available in APSIM (soiltemp2). The reason was twofold: (a) soiltemp2 has been found to perform better in Iowa (Archontoulis et al., 2014a; Basche et al., 2016; Dietzel et al., 2016); and (b) soiltemp2 better represents reality than the default model as it accounts for soil temperature changes due to tillage, residue cover, and management practices. Second, we replaced SoilWat with the SWIM soil water model (Huth et al., 2012) available in APSIM. This model allowed simulation of fluctuating shallow groundwater tables, which in this region varies from about 80 to 200 cm (Groundwater, USGS, Iowa Water Science Center). Third, we improved the simulation of soybean residue C:N ratio at harvest because the simulated C:N ratio was low when compared to published data (Johnson et al., 2007) and caused an over-prediction of corn yields in the SC rotation with no N applied. We improved soybean C:N ratio by decreasing the critical N concentration of different plant tissues at physiological maturity by about 20% (Supplementary Table S2). Additionally, we decreased the potential N fixation rate (Supplementary Table S2) to better match seasonal N fixation estimates to those observed in the literature for this region (Salvagiotti et al., 2008). No changes were made in the corn crop model, although various options were explored via sensitivity analysis. Given all these changes we re-initialized conditions at the start of the simulation on year 1999 (Supplementary Table S3).
Estimation of the Annual Economic Optimum Nitrogen Rate
The relationship between observed or simulated yield and N rate was fit using the quadratic
In these equations, y represents corn yield (either observed or simulated), x is the fertilizer N rate, a is the intercept, b is the linear coefficient, c is the quadratic coefficient, and x0 is the N rate at the join point. The PROC NLIN procedure in SAS (Version 9.4, SAS, 2013). Equations were deemed significant at p < 0.05 and the equations with the smallest sums of squares and largest R2 were selected.
Corn EONR and the yield at the EONR (YEONR) were calculated from the N response equations by setting the first derivative of the fitted response curve equal to the historical price ratio of 5.6:1 N:corn grain price (US$ kg-1 N:US$kg-1 grain) ratio (Cerrato and Blackmer, 1990; Bullock and Bullock, 1994). The impact of the N:corn grain price ratio on EONR has been well documented in the literature (Cerrato and Blackmer, 1990; Sawyer et al., 2006). In this study, we used a fixed ratio across years similarly to other modeling studies (Basso et al., 2012). Using this approach, we calculated EONR and YEONR values for: (a) the observed data (EONR-Obs, YEONR-Obs); (b) the simulated data from the uncalibrated model (EONR-APSIM-Unc, YEONR-APSIM-Unc); and (c) the simulated data from the calibrated model (EONR-APSIM-Cal, YEONR-Cal).
Additionally, a different technique to calculate an optimal N rate was used (Basso et al., 2016). The calibrated APSIM model was ran for every 5 kg N ha-1 increments from 0 to 350 kg N ha-1 to simulate corn yields. Then the N rate at which the economic return on N was maximized [hereafter RTN (return to N approach)-APSIM] was estimated by difference: simulated yield times corn price minus fertilizer rate times N cost between two levels of N rate. A value of zero (or near zero) corresponds to the optimum N rate. Same prices for corn grain and N fertilizer was used as with the EONR technique.
The RTN-APSIM technique differs from the EONR-APSIM-Cal in the following way: EONR-APSIM-Cal estimates the economic optimum N rate through regression equations (Eq. 1–3) fitted to five simulated corn yields at 0, 68, 134, 200, and 268 kg N ha-1. The RTN-APSIM uses the ability of APSIM to run on any desired N rate increment to predict corn yield (every 5 kg N, from 0 to 350 kg N ha-1) and therefore the economic optimum N rate can be identified without use of regression equations. The RTN approach follows a similar methodology to that is currently used for corn N rate recommendations in the USA Midwest (known as the MRTN approach; Sawyer et al., 2006). The difference between RTN-APSIM and MRTN is that the regression equations for MRTN are within a database with extensive N rate response trials and associated regression equations, while RTN-APSIM generates a synthetic database, which depends on the accuracy of the model to predict yields and N response.
Estimation of Site Mean Economic Optimum Nitrogen Rate
Two methods were used to estimate the site mean EONR and YEONR: (a) we first averaged individual annual estimates of EONR and YEONR for each rotation, and then calculated the associated standard deviation (SD; across years mean); and (b) we averaged corn yields across years for each rotation and then we estimated EONR and YEONR using regression equation fitting and EONR calculation, which is an approach used for N recommendations (pooled mean, Sawyer et al., 2006). The same methods were used for RTN-APSIM, except no regression equation fitting was required.
Statistical Evaluation of Model Performance
To evaluate APSIM model goodness of fit, we used graphical and statistical methods. For the statistical evaluation, we computed the root mean square error (RMSE),
where Ō is the mean observed value, Si is the model estimated value, Oi is the observed value, and n is the number of data pairs. The RMSE summarizes the average difference between observed and predicted values, while RRMSE provides the relative difference. In both cases, the lower the value of the index the better the model performance. 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).
Factors Affecting Optimal Nitrogen Rate Inter-Annual Variability
Regression analysis was performed to identify statistical significant relationships between simulated EONR and explanatory factors. We considered three explanatory factors: yield at optimum N rate, time of N application rate relative to corn planting date, and precipitation sums over different time periods. We used R2 to evaluate predictability of optimum N rate based on the factors mentioned above.
Observed Corn and Soybean Yield Response to N Fertilizer, and Crop Rotation
Observed corn yield varied across years, N rates, and crop sequences (Figures 1 and 2). Yearly variability CC, corn yield averages across years ranged from 4.2 (N0) to 11.6 (N268) Mg ha-1 with a maximum yield response to N (difference between N0 and N268 treatment) of 7.6 Mg ha-1 (Figure 3). In SC, corn yield averages were greater for all N treatments compared to CC, and varied from 7.8 (N0) to 12.9 (N268) Mg ha-1 with a maximum yield response to N of 5.1 Mg ha-1. At N0, for individual years the largest yield difference between CC and SC was 3.6 Mg ha-1. Greater yearly variability in corn yield was observed in CC (coefficient of variation, CV = 17.8%) than in SC (CV = 12.9%). The CV decreased with increasing N rate in CC (from 24.8 to 15.4%), but was consistent across N rate in SC. Across rotations, high corn yields under non-limited N condition were obtained in wet years (precipitation above 1100 mm, e.g., 2008 and 2010; Figures 1–3) and low corn yields in dry years (precipitation below 600 mm precipitation, e.g., 2000, 2012, and 2013). Observed soybean yields varied from 2.1 to 4.8 Mg ha-1 across years and N rates (Figure 4; Supplementary Figure S2). The yearly variability in soybean yield had a CV of 19.5%. Soybean yields were not affected by N rates applied to corn (Figure 3; Supplementary Figure S2).
FIGURE 1. Corn yield response to N fertilizer for the continuous corn (CC) cropping system. The blue points with standard errors (n = 4) indicate the observations. The gray and red points are Agricultural Production Systems sIMulator (APSIM) simulations before and after calibration, respectively. Continuous lines are regression fits from Eqs. 1–3. When lines are not shown it means that Eqs. 1–3 did not converge. Relative root mean square error for both calibrated (RRMSE) and uncalibrated model (RRMSE_un) are shown for each year.
FIGURE 2. Corn yield response to N fertilizer for the soybean-corn (SC) cropping system. The blue points with standard errors (n = 4) indicate the observations. The gray, red, and green points indicate uncalibrated, calibrated and validated simulations from the Agricultural Production Systems sIMulator (APSIM) model. Continuous lines are regression fits from Eqs. 1–3. When lines are not shown it means that Eqs. 1–3 did not converge. Relative root mean square error for both calibrated (RRMSE) and uncalibrated model (RRMSE_un) are shown for each year.
FIGURE 3. Sixteen year mean crop yield response to N fertilizer rate (A, D, and G panels), and observed versus simulated crop yields across years and N rate (B, C, E, F, H, and I). Points are observations or simulations, continuous lines are regression fits from Eqs. 1–3, and broken lines show 1:1 relationship.
FIGURE 4. Economic optimum N rate (EONR) and corn yield at the EONR (YEONR) for every year in CC and SC. The EONR and YEONR estimates from observations using Eqs. 1-3 are shown as bars. Different color symbols show Agricultural Production Systems sIMulator (APSIM) model simulations: red points calibrated model, gray points uncalibrated model, and green points return to N approach (RTN) from the calibrated model.
Model Accuracy before and after Calibration
Simulation of Corn Yields
Overall, across years, N rates and crop sequences, APSIM explained from 50–69% (before calibration) to 67–88% (after calibration) of the observed variability in corn yield (Figure 3). The model agreement improved during calibration from moderate (RRMSE = 19.6%, uncalibrated) to good (RRMSE = 12.3%, calibrated) for corn yield prediction (Figure 3). In CC, the uncalibrated model simulated corn yield response to N well in 7 years (RRMSE < 15%), moderately well in 6 years (RRMSE 15–30%), and poorly in 3 years (RRMSE > 30%); while after calibration the model simulated yields well in 14 years and moderately well in 3 years (Figure 1). In SC, the uncalibrated model simulated corn yield response to N well in 10 years, moderately well in 3 years, and poorly in 2 years; while after calibration the model simulated yields well in 11 years and poorly in 4 years (Figure 2). In general the calibrated model captured the trends in the observed variability in corn yields across years (Supplementary Figure S4A) as well as the annual yield response to N rates (Figures 1 and 2).
Simulation of Soybean Yields
Given that the simulation setup was sequential and soybean was part of the CS rotation, the ability of APSIM in simulating soybean yields was also tested. The model simulated no response to N rate applied to the previous corn crop, which agrees with the observed data (Figure 4; for individual years see Supplementary Figure S3). The agreement in simulated soybean yields was moderate before and after calibration (calibrated RRMSE = 19%; Figure 3).
Simulation of Optimum N Rate and Methods Comparison
Site Mean Optimum N Rate
The calibration process improved the prediction of the site mean EONR in the SC but not in CC (Table 1; Figure 3). The simulated EONR (both calibrated and uncalibrated versions) was overestimated in CC and underestimated in SC (Table 1). The absolute difference in site mean EONR between simulated and observed values was smaller in SC; -39 and 18 kg N ha-1 for CC and SC, respectively, before calibration and -41 and 10 kg N ha-1 for CC and SC, respectively, after calibration (Table 1). In addition, the simulated EONR SD was high with the APSIM-Unc, largely due to mis-estimation of some years as non-N responsive.
TABLE 1. Mean economic optimum N rate (EONR, kg N ha-1) across 16-years for: observed, Obs; un-calibrated Agricultural Production Systems sIMulator (APSIM) model, Unc; calibrated model, Cal; and the return to N approach from the calibrated model, RTN.
Annual Optimum N Rate
The calculated EONR-Obs (from observations) was highly variable from year to year and ranged from 123 to 268 kg N ha-1 in CC and from 42 to 241 kg N ha-1 in SC (Figure 4). The inter-annual variability in EONR-Obs was greater in SC than in CC (CV of 32 vs. 22%, respectively).
The calculated EONR from the APSIM model followed some of the observed annual trends (Figure 4), with the prediction error to be larger in SC than CC (Supplementary Figure S7). In CC, the RMSE ranged from 63 kg N ha-1 before calibration to 56 kg N ha-1 after calibration. In SC, the RMSE ranged from 83 kg N ha-1 before calibration to 68 kg N ha-1 after calibration. Interestingly, the two methods of calculating EONR from modeled yields (via regression Eqs. 1–3 or via the RTN approach) had similar RMSE and RRMSE values across years, but the annual predictions of optimum N using the RTN approach were less variable across years (Figure 4). These results show that for year-to-year simulation of EONR, the calibrated version should be used either via Eqs. 1–3 with regression analysis or the RTN approach. Overall the calibration process reduced the RRMSE in annual EONR predictions by 14.2% in CC and 10.3% in SC (Supplementary Figure S8).
The calculated yearly YEONR-Obs (from observation) was less variable compared to the EONR variability (CV of 17 and 12% for CC and SC, respectively, Figure 4). The simulated YEONR followed the observed annual trends well (Figure 4; RMSE of 1.88 Mg ha-1 before calibration and 1.41 Mg ha-1 after calibration). The model simulated YEONR was more accurate than EONR. In relative terms, the error in YEONR prediction was about four times lower than the error in EONR prediction (Supplementary Figures S6 and S7). However, there was no correlation between these errors (Supplementary Figure S6).
Use of the RTN approach to compute the optimum N rate, and compared to the simulated calibrated values (Table 1), produced a closer EONR in CC to the observed EONR (-8 kg N ha-1), but a greater difference in SC (18 kg N ha-1). Unlike the APSIM-Cal and APSIM-Unc simulations, the RTN-APSIM did not over-estimate EONR in CC, but underestimated in SC (Table 1).
Factors Causing Yearly Variability in Optimal Nitrogen Rate
The YEONR-Obs (Supplementary Figure S5), precipitation (Figure 5), and the time of N application (Supplementary Figure S8) were explored as possible factors to explain inter-annual variability in EONR. There was a significant positive relationship between spring precipitation and EONR-Obs but the relationship had low predictive power (p < 0.05; R2 = 0.27–0.45; Figure 5). Spring precipitation, defined here as precipitation accumulated from April 1 to June 31, was selected from among many other precipitation intervals explored in this study as the best predictor of inter-annual EONR variability (Supplementary Figure S8). The YEONR, time of N rate application, the July precipitation (15 days window around corn silking), and combinations of those factors (including spring precipitation) via multi-factor regression modeling did not result in any significant correlation.
FIGURE 5. Cumulative spring precipitation from April 1st to June 30th (every year) versus economic optimum N rate (observed and simulated EONR, circles and squares, respectively; top panels), simulated spring soil N supply (from soil organic carbon mineralization; middle panels), and simulated spring N loss (denitrification and leaching; lower panels) for CC and SC crop sequences.
The calibrated APSIM version showed a similar relationship between EONR and spring precipitation as with EONR-Obs (Figure 5) and therefore the model was used to provide insights into factors causing this relationship. Soil net N mineralization (simulated N supply), and the sum of denitrification and N leaching below 1 m depth (simulated N loss) were used as explanatory variables. The model indicated that the relationship between EONR and spring precipitation was primarily caused by an exponential increase in simulated N loss and to some extent by a reduction in simulated N supply with increasing spring precipitation (Figure 5). The model also showed that the rate of the simulated N supply reduction with increased precipitation was similar between rotations. Furthermore model analysis showed that the level of simulated N supply was 50% higher in SC than CC, which explains the lower EONR values typically found in SC systems (Table 1; Figure 3).
Calibration Strategy and Steps
Evaluating a model against long-term data is critical when the model is to be used for N management. This is because processes such as N mineralization, require several years to be sufficiently evaluated (Jenkinson et al., 1994; Leigh et al., 1994; Körschens, 2006) and can differentially affect N response among years. Our study is among a few in the literature that tests a process-based model in the long-term (Ma et al., 2007). The long-term data were powerful in detecting weakness in the model, but did not provide guidance on which of the model’s processes or parameters needed to be improved (Kersebaum et al., 2015). Therefore, during calibration we aimed to improve the overall representation of the system based on previous knowledge of the site (for example, C:N ratio of soybean and corn residue, phenology, etc.) rather than just optimizing cultivar parameters by year to better fit the observed data within the study range. This strategy is robust and allows the calibrated model to be used outside the study period (future years) with confidence at this site.
During calibration we implemented the alternate soil water (SWIM) and temperature (soiltemp) models available in the framework, and changed parameters influencing soybean residue C:N ratio (Table 2; Figure 3). Among changes made in the model, the activation of fluctuating water table via the SWIM soil water model was found to be the most important (e.g., see improvements in yield prediction from 2012 drought in Figures 2 and 3). Yet, few models have this capability, despite the great importance of water table depth on water balance, N dynamics, and crop growth (Kalita and Kanwar, 1992; Varvel et al., 1997; Hefting et al., 2004; Kahlown et al., 2005; Nosetto et al., 2009; Portela et al., 2009). Simulations of the groundwater table depth (Supplementary Figure S2I) were reasonable judging measurements in nearby sites (Hatfield et al., 1999; Helmers et al., 2012; Archontoulis et al., 2016).
TABLE 2. Statistical evaluation of the uncalibrated and calibrated APSIM model performance in simulating corn and soybean yields by N rate across 16-years.
The simulated soybean residue C:N ratio was initially low (∼20, Supplementary Figure S2) compared to literature values (25–40; Al-Kaisi et al., 2005; Bichel, 2013; Li et al., 2013). The low C:N ratio occurred mainly because APSIM supplies enough N through fixation to ensure non-N limiting soybean growth, and thus no response of soybean yield to prior-year N application to corn (potential residual inorganic-N (Figure 4; Supplementary Figure S2). This effect resulted in simulated luxurious N uptake in plant tissues and therefore low C:N ratio of the soybean residue. After calibration, the soybean residue C:N ratio increased to reasonable values (around 30; Supplementary Figure S2D), the simulation of the annual N fixation decreased to realistic estimates (around 180 kg N ha-1 year-1, Salvagiotti et al., 2008; Christianson et al., 2012), while the model maintained good performance in terms of N fixation and yield response to prior-year corn N fertilization (Figure 4; Supplementary Figures S2G and S3). We believe these changes improve the representation of N fixation (Chen et al., 2016) and soybean residue in the model. However, more experimental work is needed to verify these changes and improve the simulation of soybean rotation effects in APSIM.
APSIM Performance in Simulating Yields before and after Calibration
This study quantified accuracy of both calibrated and uncalibrated versions in order to show the degree of improvement possible with the use of long-term data. The ability of a model to predict crop yield in any environment depends on the given inputs (soil, weather, and management), how well the model structure represents reality, and the model parameters. To evaluate and test APSIM model performance, we used a system approach that explicitly considered available experimental data (corn and soybean yield, and SOC measurements) and also using literature information to evaluate additional model processes such as N fixation, root/shoot ratio, N concentrations, phenology, and others (Supplementary Figures S3 and S4). Interestingly, the model simulated long-term SOC change equally before and after calibration (RRMSE < 5%). The reason is likely that for modeling SOC, the cumulative carbon input change over time is most important (e.g., Luo et al., 2011). Annual over- and under-prediction of yields and corresponding carbon input are compensated over time if the long-term site mean yield prediction is similar before and after calibration.
The fact that APSIM simulated well yields and crop yield response to N (Figures 1–3), provided strong evidence that important soil and crop processes were being accounted for in the model. Otherwise, the model would consistently produce large under- or over-estimated yields values, resulting in different patterns across time compared to the observations (Supplementary Figure S4A). It is important to note the good simulation of grain yield with no N fertilizer input across the 16-years (Supplementary Figure S4A), which provides evidence that soil N supply and N uptake were well simulated by the model. Furthermore, the model simulated greater net soil N mineralization in the SC rotation than in CC, which is in line with literature reports (Bundy et al., 1993; Schoessow et al., 2010). Difference in net soil N mineralization (caused by residue amount and C:N ratio) was the main cause of EONR difference between CC and SC as the simulated N loss was found to be about the same in both rotations (Figure 5).
As expected, APSIM performance in simulating crop yields improved after calibration: RRMSE decreased from 15–30% to below 15%; Figure 3; Table 2). These performance evaluation results for crop yields are comparable to those reported in the literature for other models (Ahmed et al., 2007; Thorp et al., 2008; Liu et al., 2011; Yang et al., 2014). For a fair judgment of model performance, we should also mention the following assumptions, and those unknowns that may have an impact on model results: (a) we used a fixed cultivar focused on representing well the phenology across the 16-years (Supplementary Figure S2); however, different cultivars were used in the experiment (Supplementary Table S1) and some probably had different physiological characteristics; (b) there were unknowns in plant population at harvest, tillage date, and depth; (c) there was likely to be abiotic stresses in some years, hail storm damage (2013, Figure 2), and lodging issues (2002 and 2004 in CC, Figure 1) that were not considered within the model.
Corn yield predictability with calibrated APSIM (in particular for CC) increased at high N rates in particular for CC (see RMSE, Table 2). This occurred mainly because at high N rates the model has to account for only water limitations to crop growth, while at low N rates both water and N (and their interactions) become limiting factors to crop growth. This is also in accordance with published results from other crop models (Timsina and Humphreys, 2006; Liu et al., 2011; Yang et al., 2013; Li et al., 2015a,b).
Modeling Optimal Nitrogen Rate
Simulating EONR was more sensitive and complex, and had more associated uncertainty, than simulating yields (Supplementary Figures S6 and S7). This occurred because identification of the optimum N rate and associated yield in the yield-N response relationship – were quite dependent on the small incremental change (slope) in yield as N rate approached the maximum response. Over- or under-estimation of simulated yields around the optimum N rate resulted in deviations in model-derived EONR values (Eqs. 1–3). For example, in year 2002 the RRMSE for CC yield predictions by the calibrated model across N rates was 8.6% (Figure 1). This variation resulted in a 50% RRMSE in EONR prediction and in a 9.3% RRMSE in YEONR prediction. The difficulty in accurately predicting EONR from five simulated yield points is no different than uncertainties associated with the selection of regression equations to describe yield response to N with observed yields, and thus can affect APSIM estimation of the agronomic and economic optimum N rate (Waugh et al., 1973; Cerrato and Blackmer, 1990; Kyveryga, 2005; Scharf et al., 2005; Hernandez and Mulla, 2008).
However, the RTN-APSIM technique that did not use regression fitted equations had similar RMSE and RRMSE values as with EONR-APSIM-Cal, and was in close agreement with the across years mean (pooled) EONR-Obs (Table 1). A main difference between the two techniques was that the RTN-APSIM was less variable from year to year, especially for CC (Figure 4). The lower inter-annual variability in optimum N estimates from the RTN-APSIM method could be attributed to the N rate increments used in the calculations (5 kg vs. 67 kg N increments) and computation method differences, which can affect the identification of the optimum N point in yield response to N rate (Bachmaier, 2012). We concluded that the differences between simulated and observed annual optimum N rate values (Supplementary Figure S7) are due to over- and under-estimates of corn yields, especially surrounding the N rate inflection point, and to a smaller extent, due to the sensitivity of Eqs. 1–3 used to estimate EONR.
The 16-year mean differences in EONRs, especially for RTN-APSIM (Table 1), which could be called estimated errors, are acceptable within historical and current N rate ranges (46–56 kg N ha-1; ± 23–28 kg N ha-1) suggested for corn (Voss and Shrader, 1979; Sawyer et al., 2006) that includes uncertainty in estimation of optimal N (note that the range also depends on the fertilizer: corn price ratio). This means that the APSIM model can be used as a tool to assist optimum N rate recommendations in this USA region.
An important question is how the model can be used to add value within existing N rate guidelines. The main problem with EONR estimation is that the determination is made after crop harvest when yields are known. However, rate guidance is needed before N application in the fall or spring before and after planting. Thus farmers and crop advisers use guidelines based on extensive numbers of N rate research trials (MRTN; Sawyer et al., 2006). This makes the estimation of the site mean EONR very important in this study, given also that a large portion of Midwestern farmers apply N before crop planting. The APSIM model can assist N rate decisions via the following pathways. First, if the objective of long-term experiments is to derive site mean EONR recommendations, then the model can assist in this task (Figure 3). Given that the calibration processes improved more the yearly EONR prediction (14.2 and 10.3% reduction in the RRMSE for CC and SC, respectively, see Supplementary Figure S7) than site mean EONR prediction (no improvement for CC and 5% reduction in the RRMSE for SC rotation using the pooled mean), this study provides an encouraging result for model usability if only minimum site-information is available.
Second, APSIM has the potential to predict in real-time soil nitrate dynamics within the soil profile and this information could be used to adjust early-to-mid season N application rates (Archontoulis et al., 2016). This approach is currently being used by commercial companies. Third, since APSIM can predict grain yields early in the season using a range of possible weather conditions (actual, historical, future; Archontoulis et al., 2016), it could also predict needed N rates based on yield predictions as it is currently being applied in Australia as decision-support tool (Yield Prophet; Carberry et al., 2009). Nitrogen rate accuracy from yield prediction would be highly dependent on model yield predictability, and needs to be confirmed with an additional study. The value added by models and the accuracy in predicting needed N rates will always be greater when models are supported by local experimental data to periodically check performance and allow updates in the model algorithms or parameters to deal with new genetics and changes in soil and weather over time (Ahuja et al., 2014).
Causes of Optimal Nitrogen Rate Variability
In addition to predictability, deeper understanding of the factors causing EONR inter-annual variability is important for optimizing agronomic, economic, and environmental outcomes. Among three factors explored with data available for this study (time of N application, optimum yield, and precipitation periods), the cumulative May to June (spring) precipitation explained yearly EONR variability (Figure 5). However, the relatively low predictability of this relationship might be due to other confounding factors such as time of N application and planting date, which were not constant over the 16-year period in this study. The relationship between spring precipitation and EONR found in this study agreed with other studies conducted in rainfed environments (Vanotti and Bundy, 1994; Piekielek et al., 1995; Kachanoski et al., 1996; Bundy, 2000; Lory and Scharf, 2003; Sawyer et al., 2006; Scharf et al., 2006), but not with studies conducted in irrigated regions where yield level (optimum yield) was the main driver for the inter-annual variability in EONR (Dobermann et al., 2003; Gehl et al., 2005). Interestingly, the July precipitation which reflects the ± 15 day period around corn silking (see APSIM diagnostics Supplementary Figure S2) was not correlated with EONR yearly variability (R2< 0.25, p < 0.05; Supplementary Figure S9) despite the great importance of this period for kernel number determination and corn grain yields (Edmeades et al., 2000; Andrade et al., 2002; Calviño et al., 2003). This would be attributed to high soil moisture capacity of the soil and shallow groundwater tables in this region that can compensate for period of water stress and also due to the fact that corn takes up about 70% of its total N uptake by silking (Ciampitti and Vyn, 2012; Woli et al., 2016).
Variability in EONR and its relationship with spring precipitation in the US Midwest has typically been associated in previous research with an increase in N loss with high spring precipitation but previous studies lacked comprehensive measurements (Meisinger, 1984; Eghball and Varvel, 1997; Kay et al., 2006; Lawlor et al., 2008). The APSIM model analysis explicitly quantified the shape and magnitude of N loss per mm of precipitation and indicated that the shape of the relationship is similar in CC and SC systems (Figure 5).
The ability of APSIM and other mechanistic process-based models to simulate and explain the effect of precipitation on simulated N loss and supply, and thus, the impact on N response (Figure 5) becomes even more relevant with future climate change scenarios. For the US Midwest several studies have predicted higher frequency of both drought and flood events (Schoof et al., 2010; Kunkel et al., 2013; Dai et al., 2015). In this context, long-term simulations with different weather allow to capture the ranges of yield responses to N rates and choose the optimal rate that more frequently provides the best outcomes in terms of higher yield and lower nitrate leaching (Basso et al., 2016). Furthermore, the APSIM model can generate predictions across different weather scenarios that could be used to inform potential need for changes in future N management decisions.
Model analysis of a 16-year field-experiment dataset that included crop yields and SOC values with five N fertilizer rates and two crop sequences revealed the following main findings:
(1) The fact that APSIM simulated well crop yields and crop yield response to N rate, provided strong evidence that important soil and crop processes were being accounted for in the model;
(2) Model calibration (implementation of SWIM soil water model with activation of soil water table, use of soil temperature 2 model, and improvements in soybean residue C:N ratio) reduced the simulation error (RRMSE) in crop yield prediction by 9% and the annual EONR prediction by 12%. We also found that SOC prediction was insensitive to calibration when long-term mean crop yield was simulated well;
(3) The optimum N rate was higher for CC than SC and according to the model analysis this is associated with higher SOC net mineralization in the SC rotation.
(4) Simulation of EONR was more sensitive and complex than simulating crop yield. Results suggest that for long-term site mean EONR predictions both versions (calibrated and uncalibrated) can be used, while for accurate year-by-year simulation of EONR the calibrated version should be used. Use of the RTN-APSIM approach (small N rate increment with no regression fit) for optimal rate estimation had similar performance compared to EONR-APSIM-Cal approach (five N rate-points and regression fit). A main difference in optimal N rate estimation between the two techniques was that the RTN-APSIM output was less variable from year to year.
(5) Five potential applications were identified where the model could assist N management: (a) estimation of long-term mean EONR; (b) simulation of N dynamics (soil N available and crop N demand); (c) prediction of optimal N using a range of possible weather conditions; (d) simulation of climate change impact on optimal N need;
(6) The APSIM model can be used to explore and explain factors causing inter-annual variability in EONR. For example, the model showed that in rainfed corn-based systems in Iowa, the higher the spring precipitation (April to June) the higher the EONR because simulated N loss via denitrification and leaching increased exponentially while simulated N supply via mineralization tended to decrease.
Finally, for rainfed corn-based systems in the USA Midwest, a combination of process-based modeling, coupled with existing N rate recommendation methods and field data, may be the best approach to fine tune optimal N rate guidance for corn and to develop future management-based strategies under climate change scenarios for maximizing agronomic, economic, and environmental outcomes.
JS and DB collected crop data. HP and MC collected soil data. LP, JS, and SA designed the modeling study. LP and SA performed the model analysis. LP synthesized results and constructed tables, figures for the paper. LP and SA created the first draft and all co-authors (JS, DB, RD, HP, MC, PT, and KM) contributed to the final version of the manuscript.
This work was part of the Agriculture and Food Research Initiative Hatch project No. 1004346 and was also partially supported by the Plant Science Institute, and the Brown Graduate Fellowship program of Iowa State University.
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.
We thank the Agronomy Research Farm managers and personnel for their assistance with the research over the years. We also thank the Iowa State University College of Agriculture and Life Sciences for partial financial support of the long-term field research trial.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01630/full#supplementary-material
Ahmed, I., Rudra, R., McKague, K., Gharabaghi, B., and Ogilvie, J. (2007). Evaluation of the root zone water quality model (RZWQM) for southern ontario: part I. Sensitivity analysis, calibration, and validation. Water Qual. Res. J. Canada 42, 202–218.
Ahuja, L. R., Ma, L., Lascano, R. J., Saseendran, S. A., Fang, Q. X., Nielsen, D. C., et al. (2014). “Syntheses of the current model applications for managing water and needs for experimental data and model improvements to enhance these application,” in Practical Applications of Agricultural System Models to Optimize the Use of Limited Water (Madison, WI: American Society of Agronomy), 399–438.
Al-Kaisi, M. M., Yin, X. H., and Licht, M. A. (2005). Soil carbon and nitrogen changes as affected by tillage system and crop biomass in a corn–soybean rotation. Appl. Soil Ecol. 30, 174–191. doi: 10.1016/j.apsoil.2005.02.014
Anapalli, S., Ahuja, L., Ma, L., and Trout, T. (2014). Modeling the effects of irrigation frequencies, initial water and nitrogen on corn yield responses for best management. Adv. Agric. Sys. Model. 1, 26–52.
Andrade, F. H., Echarte, L., Rizalli, R., Della Maggiora, A., and Casanovas, M. (2002). Kernel number prediction in maize under nitrogen or water stress. Crop Sci. 42, 1173–1179. doi: 10.2135/cropsci2002.1173
Archontoulis, S. V., 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 Proceedings of the International Crop Modelling Symposium, 15-17, Berlin, 188–189.
Archontoulis, S. V., Miguez, F. E., and Moore, K. J. (2014a). Evaluating APSIM maize, soil water, soil nitrogen, manure, and soil temperatures modules in the Midwestern United States. Agron. J. 106, 1025–1040. doi: 10.2134/agronj2013.0421
Archontoulis, S. V., Miguez, F. E., and Moore, K. J. (2014b). A methodology and an optimization tool to calibrate phenology of short-day species included in the APSIM PLANT model: application to soybean. Environ. Model. Softw. 62, 465–477. doi: 10.1016/j.envsoft.2014.04.009
Basche, A. D., Archontoulis, S. V., Kaspar, T. C., Jaynes, D. B., Parkin, T. B., and Miguez, F. E. (2016). Simulating long-term impacts of cover crops and climate change on crop production and environmental outcomes in the Midwestern United States. Agric. Ecosyst. Environ. 218, 95–106. doi: 10.1016/j.agee.2015.11.011
Basso, B., Cammarano, D., Troccoli, A., Chen, D., and Ritchie, J. T. (2010). Long-term wheat response to nitrogen in a rainfed Mediterranean environment: field data and simulation analysis. European J. Agron. 33, 132–138. doi: 10.1016/j.eja.2010.04.004
Basso, B., Dumont, B., Cammarano, D., Pezzuolo, A., Marinello, F., and Sartori, L. (2016). Environmental and economic benefits of variable rate nitrogen fertilization in a nitrate vulnerable zone. Sci. Total Environ. 54, 227–235. doi: 10.1016/j.scitotenv.2015.12.104
Basso, B., Ritchie, J. T., Grace, P. R., and Sartori, L. (2006). Simulation of tillage systems impact on soil biophysical properties using the SALUS model. Italian J. Agron. 4, 677–688. doi: 10.4081/ija.2006.677
Basso, B., Sartori, L., Cammarano, D., Fiorentino, C., Grace, P. R., Fountas, S., et al. (2012). Environmental and economic evaluation of N fertilizer rates in a maize 15 crop in Italy: a spatial and temporal analysis using crop models. Biosyst. Eng. 113, 103–111. doi: 10.1016/j.biosystemseng.2012.06.012
Batchelor, W. D., Basso, B., and Paz, J. O. (2002). Examples of strategies to analyze spatial and temporal yield variability using crop models. Eur. J. Agron. 18, 141–158. doi: 10.1016/S1161-0301(02)00101-6
Brown, K. H., Bach, E. M., Drijber, R. A., Hofmockel, K. S., Jeske, E. S., Sawyer, J. E., et al. (2014). A long-term nitrogen fertilizer gradient has little effect on soil organic matter in a high-intensity maize production system. Glob. Change Biol. 20, 1339–1350. doi: 10.1111/gcb.12519
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. (2000). “Nitrogen recommendations and optimum nitrogen rates: how do they compare?,” in Proceedings of the North Central Ext.-Industry Soil Fertil. Conference, Vol. 16, (St. Louis, MO: Potash and Phosphate Inst), 5–13.
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
Chen, C., Lawes, R., Fletcher, A., Oliver, Y., Robertson, M., Bell, M., et al. (2016). How well can APSIM simulate nitrogen uptake and nitrogen fixation of legume crops? Field Crop Res. 187, 35–48. doi: 10.1016/j.fcr.2015.12.007
Christianson, L., Castellano, M., and Helmers, M. (2012). “Nitrogen and phosphorus balance in Iowa cropping systems: sustaining Iowa’s soil resource,” in Proceedings of the Integrated Crop Management Conference, (Ames, IA: Iowa State University Extension).
Ciampitti, I. A., and Vyn, T. J. (2012). Physiological perspectives of changes over time in maize yield dependency on nitrogen uptake and associated nitrogen efficiencies: a review. Field Crops Res. 133, 48–67. doi: 10.1016/j.fcr.2012.03.008
Dai, S., Shulski, M. D., Hubbard, K. G., and Takle, E. S. (2015). A spatiotemporal analysis of Midwest US temperature and precipitation trends during the growing season from 1980–2013. Int. J. Climatol. 36, 517–525. doi: 10.1002/joc.4354
Dietzel, R., Liebman, M., Ewing, R., Helmers, M., Horton, R., Jarchow, M., et al. (2016). How efficiently do corn- and soybean-based cropping systems use water? A systems modeling analysis. Glob. Change Biol. 22, 666–681.
Dobermann, A., and Cassman, K. G. (2002). Plant nutrient management for enhanced productivity in intensive grain production systems of the United States and Asia. Plant Soil 247, 153–175. doi: 10.1023/A:1021197525875
Dobermann, A., Witt, C., Abdulrachman, S., Gines, H. C., Nagarajan, R., Son, T. T., et al. (2003). Estimating indigenous nutrient supplies for site-specific nutrient management in irrigated rice. Agron. J. 95, 924–935. doi: 10.2134/agronj2003.0924
Dumont, B., Basso, B., Bodson, B., Destain, J. P., and Destain, M. F. (2016). Assessing and modeling economic and environmental impact of wheat nitrogen management in Belgium. Environ. Model. Softw. 79, 184–196. doi: 10.1016/j.envsoft.2016.02.015
Edmeades, G. O., Bänziger, M., and Ribaut, J. M. (2000). “Maize improvement for drought-limited environments,” in Physiological Bases for Maize Improvement, eds M. E. Otegui and G. A. Slafer (New York: Food Products Press), 75–111.
Eghball, B., and Varvel, G. E. (1997). Fractal analysis of temporal yield variability of crop sequences: implications for site-specific management. Agron. J. 89, 851–855. doi: 10.2134/agronj1997.00021962008900060001x
Ferguson, R. B., Hergert, G. W., Schepers, J. S., Gotway, C. A., Cahoon, J. E., and Peterson, T. A. (2002). Site-specific nitrogen management of irrigated maize: yield and soil residual nitrate effects. Soil Sci. Soc. Am. J. 66, 544–553. doi: 10.2136/sssaj2002.5440
Fountas, S., Wulfsohn, D., Blackmore, B. S., Jacobsen, H. L., and Pedersen, S. M. (2006). A model of decision-making and information flows for information-intensive agriculture. Agric. Sys. 87, 192–210. doi: 10.1016/j.agsy.2004.12.003
Franzen, D. W., Sharmab, L. K., Bu, H., and Dentond, A. (2016). Evidence for the ability of active-optical sensors to detect sulfur deficiency in corn. Agron. J. 108, 2158–2162. doi: 10.2134/agronj2016.05.0287
Gowda, P. H., Chavez, J. L., Colaizzi, P. D., Evett, S. R., Howell, T. A., and Tolk, J. A. (2008). ET mapping for agricultural water management: present status and challenges. Irrigation Sci. 26, 223–237. doi: 10.1007/s00271-007-0088-6
Haghiri, R., Miller, R. H., and Logan, T. J. (1978). Crop response and quality of leachate as affected by land application of beef cattle waste. J. Environ. Qual. 7, 406–412. doi: 10.2134/jeq1978.00472425000700030023x
Hammer, G. L., Dong, Z., and McLean, G. (2009). Can changes in canopy and/or root system architecture explain historical maize yield trends in the U.S. Corn Belt? Crop Sci. 49, 299–312. doi: 10.2135/cropsci2008.03.0152
Hatfield, J. L., Jaynes, D. B., Burkart, M. R., Cambardella, C. A., Moorman, T. B., Prueger, J. H., et al. (1999). Water quality in walnut creek watershed: setting and farming practices. J. Environ. Qual. 28, 11–24. doi: 10.2134/jeq1999.28111x
Hefting, M. M., Clement, J. C., Dowrick, D., Cosandry, A. C., Bernal, S., Cimpian, C., et al. (2004). Water table elevation controls on soil nitrogen cycling in riparian wetlands along a European climatic gradient. Biogeochemistry 67, 113–134. doi: 10.1023/B:BIOG.0000015320.69868.33
Helmers, M., Christianson, R., Brenneman, G., Lockett, D., and Pederson, C. (2012). Water table, drainage, and yield response to drainage water management in southeast Iowa. J. Soil Water Conserv. 67, 495–501. doi: 10.2489/jswc.67.6.495
Holland, E. A., Braswell, B. H., Sulzman, J., and Lamarque, J. (2005). Nitrogen deposition on to the United States and Western Europe: synthesis of observations and models. Ecol. Appl. 15, 38–57. doi: 10.1890/03-5162
Holzworth, D. P., Huth, N. I., deVoil, P. G., Zurcher, E. J., Herrmann, N. I., McLean, G., 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
Iowa Environmental Mesonet (2014). Iowa Ag Climate Network. Ames, IA: Iowa State University. Available at: http://mesonet.agron.iastate.edu/agclimate/ (accessed 10 May, 2014).
Jenkinson, D. S., Bradbury, N. J., and Coleman, K. (1994). “How the rothamsted classical experiments have been used to develop and test models for the turnover of carbon and nitrogen in soil,” in Long-term Experiments in Agricultural and Ecological Sciences, eds R. A. Leigh and A. E. Johnston (Wallingford, CT: CAB International), 117–137.
Jones, J. W., Hoogenboom, G., Porter, C. H., Boote, K. J., Batchelor, W. D., Hunt, L. A., et al. (2003). The DSSAT cropping system model. Euro. J. Agron. 18, 235–265. doi: 10.1016/S1161-0301(02)00107-7
Kay, B. D., Mahloubi, A. A., Beauchamp, E. J., and Dharmakeerthi, R. S. (2006). Integrating soil and weather data to describe variability in plan available nitrogen. Soil Sci. Soc. Am. J. 70, 1210–1221. doi: 10.2136/sssaj2005.0039
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., Boote, K. J., Jorgenson, J. S., Nendel, C., Bindi, M., Frühauf, C., et al. (2015). Analysis and classification of data sets for calibration and validation of agro-ecosystem models. Environ. Modell. Softw. 72, 402–417. doi: 10.1016/j.envsoft.2015.05.009
Kunkel, K. E., Stevens, L. E., Stevens, S. E., Sun, L., Janssen, E., Wuebbles, D., et al. (2013). “Regional climate trends and scenarios for the U.S. national climate assessment: part 3,” in Climate of the Midwest U.S. NOAA Technical Report NESDIS 142-3 (Washington, D.C: National Oceanic and Atmospheric Administration), 103.
Kyveryga, P. M., Blackmer, A. M., and Zhang, J. (2009). Characterizing and classifying variability in corn yield response to nitrogen fertilization on subfield and field scales. Agron. J. 101, 269–277. doi: 10.2134/agronj2008.0168
Lambert, D. M., Lowenberg-DeBoer, J., and Malzer, G. L. (2006). Economic analysis of spatial–temporal patterns in corn and soybean response to nitrogen and phosphorus. Agron. J. 98, 43–54. doi: 10.2134/agronj2005.0005
Lawlor, P. A., Helmers, M. J., Baker, J. L., Melvin, S. W., and Lemke, D. W. (2008). Nitrogen application rate effect on nitrate-nitrogen concentration and loss in subsurface drainage for a corn-soybean rotation. Trans. ASABE 51, 83–94. doi: 10.13031/2013.24229
Leigh, R. A., Prew, R. D., and Johnston, A. E. (1994). “The management of long-term agricultural field experiments: procedures and policies evolved from the Rothamsted classical experiments,” in Long-Term Experiments in Agricultural and Ecological Sciences, eds R. A. Leigh and A. E. Johnston (Wallingford: CAB International), 253–268.
Li, L. J., Han, X. Z., You, M. Y., Yuan, Y. R., Ding, X. L., and Qiao, Y. F. (2013). Carbon and nitrogen mineralization patterns of two contrasting crop residues in a Mollisol: effects of residue type and placement in soils. Eur. J. Soil. Biol. 54, 1–6. doi: 10.1016/j.ejsobi.2012.11.002
Li, T., Hasegawa, T., Yin, X., Zhu, Y., Boote, K., Adam, M., et al. (2015a). Uncertainties in predicting rice yield by current crop models under a wide range of climatic conditions. Glob. Change Biol. 21, 1328–1341. doi: 10.1111/gcb.12758
Li, Z. T., Yang, J. Y., Drury, C. F., and Hoogenboom, G. (2015b). Evaluation of the DSSAT-CSM for simulating yield and soil organic C and N of a long-term maize and wheat rotation experiment in the Loess Plateau of Northwestern China. Agric. Syst. 135, 90–104. doi: 10.1016/j.agsy.2014.12.006
Liu, H. L., Yang, J. Y., Drury, C. F., Reynolds, W. D., Tan, C. S., Bai, Y. L., et al. (2011). Using the DSSAT-CERES-Maize model to simulate crop yield and nitrogen cycling in fields under long-term continuous maize production. Nutr. Cycl. Agroecosys. 89, 313–328. doi: 10.1007/s10705-010-9396-y
Liu, S., Yang, J. Y., Zhang, X. Y., Drury, C. F., Reynolds, W. D., and Hoogenboom, 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
Lobell, D. B., Hammer, G. L., McLean, G., Messina, C., and Roberts, M. J. W. S. (2013). Understanding the critical role of extreme heat for maize production in the United States. Nat. Clim. Change 3, 497–501. doi: 10.1038/nclimate1832
Luo, Z., Wang, E., Sun, O. J., Smith, C. J., and Probert, M. E. (2011). Modeling long-term soil carbon dynamics and sequestration potential in semi-arid agro-ecosystems. Agric. Forest Meteorol. 151, 1529–1544. doi: 10.1016/j.agrformet.2011.06.011
Ma, B. L., and Dwyer, L. M. (1999). Within plot variability in available soil mineral nitrogen in relation to leaf greenness and yield. Commun. Soil Sci. Plant Anal. 30, 1919–1928. doi: 10.1080/00103629909370342
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
Malone, R. W., Jaynes, D. B., Ma, L., Nolan, B. T., Meek, D., and Karlen, D. (2010). Soil test N recommendation augmented with PEST-optimized RZWQM simulations. J. Environ. Qual. 39, 1711–1723. doi: 10.2134/jeq2009.0425
Malone, R. W., Ma, L., Heilman, P., Karlen, D. L., Kanwar, R. S., and Hatfield, J. L. (2007). Simulated N management effects on corn yield and tile-drainage nitrate loss. Geoderma 140, 272–283. doi: 10.1016/j.geoderma.2007.04.011
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.0958
Martinez-Feriaa, R. A., Dietzela, R., Liebmana, M., Helmersb, M. J., and Archontoulisa, S. V. (2016). Rye cover crop effects on maize: a system-level analysis. Field Crop Res. 196, 145–159. doi: 10.1016/j.fcr.2016.06.016
Meisinger, J. J. (1984). “Evaluating plant-available nitrogen in soil–crop systems,” in Nitrogen in Crop Production, American Society of Agronomy, Chap. 26, ed. R. A. Hauck (Madison, WI: Crop Science Society of America, Soil Science Society of America), 391–417.
Mitsch, W. J., Gilliam, J. W., Groffman, P. M., Hey, D. L., Randall, G. W., et al. (2001). Reducing nitrogen loading to the Gulf of Mexico from the Mississippi river basin: strategies to counter a persistent large-scale ecological problem. Bioscience 51, 373–388. doi: 10.1641/0006-3568(2001)051[0373:RNLTTG]2.0.CO;2
Mulvaney, R. L., Khan, S. A., Hoeft, R. G., and Brown, H. M. (2001). A soil organic nitrogen fraction that reduces the need for nitrogen fertilization. Soil Sci. Soc. Am. J. 65, 1164–1172. doi: 10.2136/sssaj2001.6541164x
Nangia, V., Gowda, R. H., Mulla, D. J., and Sands, G. R. (2008). Water quality modeling of fertilizer management impacts on nitrate losses in tile drains at the field scale. J. Environ. Qual. 37, 296–307. doi: 10.2134/jeq2007.0224
Nosetto, M. D., Jobbágy, E. G., Jackson, R. B., and Sznaider, G. (2009). Reciprocal influence between crops and shallow ground water in sandy landscapes of the Inland Pampas. Field Crop Res. 113, 138–148. doi: 10.1016/j.fcr.2009.04.016
Paz, J. O., Batchelor, W. D., Colvin, T. S., Logsdon, S. D., Kaspar, T. C., Karlen, D. L., et al. (1999). Model based techniques to determine variable rate nitrogen for corn. Agric. Syst. 60, 69–75. doi: 10.1016/S0308-521X(99)00035-9
Piekielek, W. P., Fox, R. H., Toth, J. D., and Macneal, K. E. (1995). Use of a chlorophyll meter at the early dent stage of corn to evaluate nitrogen sufficiency. Agron. J. 87, 403–408. doi: 10.2134/agronj1995.00021962008700030003x
Portela, S. I., Andriulo, A. E., Jobbágy, E. G., and Sasal, M. C. (2009). Water and nitrate exchange between cultivated ecosystems and groundwater in the Rolling Pampas. Agric. Ecosyst. Environ. 134, 277–286. doi: 10.1016/j.agee.2009.08.001
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. doi: 10.1016/S0308-521X(97)00028-0
Randall, G. W., Iragavarpu, T. K., and Schmitt, M. A. (2000). Nutrient losses in subsurface drainage water from dairy manure and urea applied for corn. J. Environ. Qual. 29, 1244–1252. doi: 10.2134/jeq2000.291350x
Rosenzweig, C., Jones, J. W., Hatfield, J. L., Ruane, A. C., Boote, K. J., Thorburn, P., et al. (2013). The agricultural model intercomparison and improvement project (AgMIP): protocols and pilot studies. Agric. For. Meteorol. 170, 166–182. doi: 10.1016/j.agrformet.2012.09.011
Salvagiotti, F., Specht, J. E., Cassman, K. G., Walters, D. T., Weiss, A., and Dobermann, A. (2008). Growth and nitrogen fixation in high yielding soybean: impact of nitrogen fertilization. Agron. J. 101, 958–970. doi: 10.2134/agronj2008.0173x
Sawyer, J., Nafziger, E., Randall, G., Bundy, L., Rehm, G., and Joern, B. (2006). Concepts and Rationale for Regional Nitrogen Guidelines for Corn, PM2015. Ames, IA: Iowa State University Extension, 27.
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., 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
Schmidt, J. P., Dellinger, A. E., and Beegle, D. B. (2009). Nitrogen recommendations for corn: an on-the-go sensor compared with current recommendation methods. Agron. J. 101, 916–924. doi: 10.2134/agronj2008.0231x
Schnebelen, N., Nicoullaud, B., Bourennane, H., Couturier, A., Verbeque, B., and Revalier, C. (2004). The STICS model to predict nitrate leaching following agricultural practices. Agronomie 24, 423–435. doi: 10.1051/agro:2004039
Schoessow, K. A., Kenneth, C. K., and Bundy, L. G. (2010). Soybean residue management and tillage effects on corn yields and response to applied nitrogen. Agron. J. 102, 1186–1193. doi: 10.2134/agronj2010.0033
Schoof, J. T., Pryor, S. C., and Suprenant, J. (2010). Development of daily precipitation projections for the United States based on probabilistic downscaling. J. Geophys. Res. 115, 1–13. doi: 10.1029/2009JD013030
Setiyono, T. D., Yangb, H. D., Waltersa, T., Dobermannc, A., Fergusona, R. B., Robertsd, D. F., et al. (2011). Maize-N: a decision tool for nitrogen management in maize. Agric. J. 103, 1276–1283. doi: 10.2134/agronj2011.0053
Thorburn, P. J., Meier, E. A., and Probert, M. E. (2005). Modelling nitrogen dynamics in sugarcane systems: recent advances and applications. Field Crops 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 Crops Res. 70, 223–232. doi: 10.1016/S0378-4290(01)00141-1
Thornton, P. E., Thornton, M. M., Mayer, B. W., Wilhelmi, N., Wei, Y., and Cook, R. B. (2012). Daymet: Daily Surface Weather on a 1 km Grid for North America, 1980–2012. Oak Ridge, TN: Oak Ridge National Laboratory Distributed Active Archive Center. doi: 10.3334/ORNLDAAC/Daymet_V2
Thorp, K. R., DeJonge, K. C., Kaleita, A. L., Batchelor, W. D., and Paz, J. O. (2008). Methodology for the use of DSSAT models for precision agriculture decision support. Comput. Electron. Agric. 64, 276–285. doi: 10.1016/j.compag.2008.05.022
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
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., and Sogbedji, J. M. (2007). “Nitrogen management for maize in humid regions: case for a dynamic approach,” in Managing Crop Nutrition for Weather, ed. T. Bruulsema (Peachtree Corners, GA: International Plant Nutrition Institute).
Varvel, G. E., Schepers, J. S., and Francis, D. D. (1997). Ability for in-season correction of nitrogen deficiency in corn using chlorophyll meters. Soil Sci. Soc. Am. J. 61, 1233–1239. doi: 10.2136/sssaj1997.03615995006100040032x
Wang, D., Prato, T., Qiu, Z., Kitchen, N. R., and Sudduth, K. A. (2003). Economic and environmental evaluation of variable rate nitrogen and lime application for claypan soil fields. Precis. Agric. 4, 35–52. doi: 10.1023/A:1021858921307
Waugh, D. L., Cate, R. B., and Nelson, L. A. (1973). Discontinuous Models for Rapid Correlation, Interpretation, and Utilization of Soil Analysis and Fertilizer Response Data. Technical Bulletin 7. Raleigh, NC: North Carolina State University.
Williams, J. R., Jones, C. A., and Dyke, P. T. (1984). “The EPIC model and its application,” in Proceedings of the International Symposium On Minimum Data Sets for Agrotechnology Transfer, March 21–26, 1983 (Hyderabad: ICRISAT Center).
Woli, K. P., Boyer, M. J., Elmore, R. W., Sawyer, J. E., Abendroth, L. J., and Barker, D. W. (2016). Corn era hybrid response to nitrogen fertilization. Agron. J. 108, 473–486. doi: 10.2134/agronj2015.0314er
Yang, J. M., Yang, J. Y., Dou, S., Yang, X. M., and Hoogenboom, G. (2013). Simulating the effect of long-term fertilization on maize yield and soil C/N dynamics in northeastern China using DSSAT and CENTURY-based soil model. Nutr. Cycl. Agroecosyst. 95, 287–303. doi: 10.1007/s10705-013-9563-z
Keywords: maize, economic optimum N rate, soybean, soil organic carbon, modeling, APSIM
Citation: Puntel LA, Sawyer JE, Barker DW, Dietzel R, Poffenbarger H, Castellano MJ, Moore KJ, Thorburn P and Archontoulis SV (2016) Modeling Long-Term Corn Yield Response to Nitrogen Rate and Crop Rotation. Front. Plant Sci. 7:1630. doi: 10.3389/fpls.2016.01630
Received: 16 July 2016; Accepted: 17 October 2016;
Published: 11 November 2016.
Edited by:Urs Feller, University of Bern, Switzerland
Reviewed by:Dimitris L. Bouranis, Agricultural University of Athens, Greece
Ignacio Antonio Ciampitti, Kansas State University, USA
Copyright © 2016 Puntel, Sawyer, Barker, Dietzel, Poffenbarger, Castellano, Moore, Thorburn 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) or licensor 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.