Abstract
The performance of crop models in simulating various aspects of the cropping system is sensitive to parameter calibration. Parameter estimation is challenging, especially for time-dependent parameters such as cultivar parameters with 2–3 years of lifespan. Manual calibration of the parameters is time-consuming, requires expertise, and is prone to error. This research develops a new automated framework to estimate time-dependent parameters for crop models using a parallel Bayesian optimization algorithm. This approach integrates the power of optimization and machine learning with prior agronomic knowledge. To test the proposed time-dependent parameter estimation method, we simulated historical yield increase (from 1985 to 2018) in 25 environments in the US Corn Belt with APSIM. Then we compared yield simulation results and nine parameter estimates from our proposed parallel Bayesian framework, with Bayesian optimization and manual calibration. Results indicated that parameters calibrated using the proposed framework achieved an 11.6% reduction in the prediction error over Bayesian optimization and a 52.1% reduction over manual calibration. We also trained nine machine learning models for yield prediction and found that none of them was able to outperform the proposed method in terms of root mean square error and R2. The most significant contribution of the new automated framework for time-dependent parameter estimation is its capability to find close-to-optimal parameters for the crop model. The proposed approach also produced explainable insight into cultivar traits’ trends over 34 years (1985–2018).
Similar content being viewed by others
Introduction
Crop models build upon the physiological understanding of plant response to the soil, weather, and management conditions. This understanding is approximated with a series of non-linear equations while the effects of different soil properties and cultivars are captured by different parameter values (input parameters to the model)1,2. Presently, there are many crop models including APSIM3, DSSAT4, EPIC5, FASSET6, SWAP and WOFOST7, SOYGRO8, and crop modeling groups9. Different models share similar modeling principles but differ in equations and parameter values used. Crop models are widely used in testing crop management strategies10, forecasting crop yield11, evaluating cultivars in breeding programs12, and exploring climate change impacts on productivity13. For example, Archontoulis et al.14 and Togliatti et al.11 used the APSIM crop model to predict crop yields in the US Corn Belt. Dumont et al.15 tested the crop yield performance of the crop model under two different weather scenarios (stochastically generated climatic data and the mean climate data).
Mechanistic crop models aim to mimic natural processes by simulating crop performance for each set of input values. When the model is used outside of a tested environment, a parameter calibration (i.e., change of cultivar or soil input values) is typically needed to improve the simulation output and be closer to observations16. The more parameters that need to be calibrated, the more complex the calibration process is. Manual calibration by trial and error17 is commonly used to estimate parameter values. Crop modelers have developed protocols (step by step approaches) in which calibration starts with phenological parameters, then continues with biomass production and partitioning related parameters14,18,19. Calibration processes need to run a crop model multiple times using an effective and iterative procedure to evaluate each parameter’s effect on crop model outputs20.
Parameter calibration algorithms used to improve crop models’ performance have been classified into two categories: local and global methods21. Local methods, called simple derivative-based methods, evaluate crop model outputs’ response to variation one parameter at a time while keeping the other parameters fixed22. The implementation of these methods is relatively easy with a low computational cost. However, they suffer several drawbacks, including the inability to specify impacts of interactions between parameters on outputs and high dependency on the parameter’s base value21. On the other hand, global methods have been developed to overcome these drawbacks by evaluating the impacts of the entire multi-dimensional parameter space on crop model outputs. These methods provide a more comprehensive view of the sensitivity analysis of parameters on outputs. Several global sensitivity analysis methods have been proposed to calibrate the parameters in crop models, including Fourier amplitude sensitivity test (FAST)23, extended FAST24, multi-normal approximations25, simulated annealing26, shuffled complex evolution method27, least squares28, regression-based model29, interaction-based model30, non-parametric smoothing31, Markov chain Monte Carlo parameter estimation20,32,33,34,35, generalized likelihood uncertainty estimation32, multi-model ensembles36, maximal conditional posterior distribution37, hybrid metropolis Hastings Gibbs algorithm38, differential evolution adaptative metropolis algorithm39,40, Bayesian model33,35,41,42,43,44, Bayesian optimization (BO)45, and Bayesian-based multilevel factorial analysis46. FAST methods have fast convergence with an acceptable computing cost and less accuracy for non-linear and complex models47,48. The shuffled complex evolution, simulated annealing, and differential evolution try to find optimal or close-to-optimal values for parameters and do not provide any information about the parameter estimates’ confidence39. Non-parametric smoothing can be applied to quantify the main effects of parameters on the response variable when there are low order interactions between parameters and response variable. Non-parametric smoothing and approximation models were not reliable for non-linear and complex models10,49. A Bayesian model needs a definition of a prior distribution where, in most cases, there is no reliable prior knowledge. Monte Carlo-based methods are computationally intensive and fail to provide adequate sample density from solution space when the number of parameters increases50.
In capturing historical trends in crop models like maize yield increase from 1985 to 2018 in the US Corn Belt, cultivar parameters need to be time-dependent to reflect changes in crop cultivars due to plant breeding or environment51,52,53,54,55,56,57,58,59,60. Using a static set of cultivar parameter values will not capture the trend of observed outcome or simulate how nature produces output from input over time in the US Corn Belt and other regions61,62. Development of time-dependent parameters is challenging because of the lack of specific data as well as the absence of a robust calibration framework for time-dependent parameters. The more time-dependent parameters that need to be calibrated, the more complex the calibration process is. To our knowledge, most of the previous model calibration exercises using optimization algorithms were not dealing with the time-dependent parameter values16,17,20,39,49. Integrating prior agronomic knowledge with an optimization model to develop a constrained model can reduce the prediction error to zero, but that would be over-fitting on the observed data set. Hence, the calibration structure must be designed in such a way to improve crop model performance and prevent the over-fitting problem. Moreover, this structure must follow some prior agronomic knowledge such that parameter trends have to follow a specific pattern over time but spatially may differ from location to location.
The goal of this paper is to develop a new time-dependent parameter calibration method for crop modeling towards improving the simulation accuracy. We designed a new automated and efficient framework and a parallel Bayesian optimization (PBO) algorithm to estimate time-dependent parameters from limited data that could integrate the prior agronomic knowledge constraints. This research describes the new algorithm and its prediction accuracy and further compares prediction accuracy against the prediction using a Bayesian algorithm and a manual calibration approach. We considered historical yield data (1984 to 2019; 34 years) from 25 locations across the US Corn Belt. Finally, we interpret time-dependent parameter values towards explaining historical yield increase in this region.
Problem definition
Crop models simulate how nature produces outcome performance from appropriate inputs using the physiological understanding of plant processes1,2; however, not all input parameters are available to inform crop models to mimic nature accurately. In capturing historical trends (e.g., yield increases) by crop models, which is the product of concurrent changes in cultivar, management, technological gains, equipment, and environment, the main challenge is to derive time-dependent parameters’ values resulting from changes in plant traits with years of plant breeding or environment58,59,60. Weather and management data are typically available. In terms of cultivar, there are good records on grain yields but not on the underlying traits, which are necessary inputs to crop models. The challenge is to approximate the crop model’s parameters with limited data available. That is why parameter calibration is an essential job in using crop models. The structure of parameter calibration in crop growth and development simulation is illustrated in Fig. 1. Although the manual and automated parameter calibration simulate outcomes that may be very consistent with observed outcomes, they could not integrate prior agronomic knowledge or time dependency of parameters within their structures. Unconstrained calibration approaches at least produce the marginal solution for the future because they ignore the time-dependency between parameters and prior agronomic knowledge during their processes. From an agronomic point of view, a constrained calibration can satisfy the reasonable trend for parameters.
Method
The proposed framework for time-dependent parameter calibration is shown in Fig. 2. The training data set from year 1 to year Y was divided into two different time-windows. Time-window p with \(n>1\) years denotes main time-window and time-window \({\hat{p}}\) with less than n years denotes auxiliary time-window. The main time-windows are defined so that two successive time-windows have \(n-1\) years overlap. For example, the first main time-window contains years 1 to n, where the second one is from year 2 to year \(n+1\) and so on. The final main time-window covers year \(Y-n+1\) to year Y. Some main time-windows have auxiliary time-windows. To estimate the parameters based on n time-windows for each year, the auxiliary time-window is defined for the first and the last training years. After defining the time-windows, we used the PBO as the optimization model in the structure of parameter calibration shown in Fig. 1 to optimize each window’s parameters to minimize the crop model’s prediction error. Then, to calculate the best combination of parameters at year y (\({\bar{P}}_y\)), we take a weighted average of parameters over windows that year y is covered by them as follows:
Taking an average of n time-windows for each year tries to reduce variance in parameter estimation because it avoids over-fitting. To obtain a consistent trend for each parameter, we fit the autoregressive model63 on the average value of the parameters during the simulation period because at most the two successive years have \(n-1\) common time-windows. Also, applying the Augmented Dickey-Fuller test64 reveals that the trend of the parameter over simulation periods has autocorrelation. The parameters’ trend cannot reject the null hypothesis that a unit root is present in a time series sample. The notation AR(t) indicates an autoregressive model of order t. Order t means future value depends on how many years before year y need to be considered.
where \(\beta _0\), \(\beta _i\), and \(\varepsilon _{y}\) are constant, the model’s parameters, and white noise at year y, respectively. For each parameter, we can use the partial autocorrelation function65 to determine the optimal order in the autoregressive model. Then we use this autoregressive model to predict parameter values for test year \(Y+1\) and run APSIM with them to simulate crop yield performance for test year \(Y+1\). In this structure, the autoregressive model is used, while this structure can be extended to other types of time series models. Moreover, the number of years in the main time-window (n) shows the trade-off between bias and variance in predicting parameters trend over simulation periods. Increasing n is to increase the bias and reduce the variance because of under-fitting on the training set and vice versa. Parameter n is the key hyper-parameter that can affect the performance of the algorithm. Optimizing n can minimize the summation of bias and variance. We tried different values for n in the case study, and \(n=5\) resulted in the best performance.
Parallel Bayesian optimization
This section proposes an extended version of BO66 to improve the time-dependent parameter calibration structure in the proposed framework (see Fig. 2). We briefly review the classic BO method and the key statistical and optimization methods on which it relies in the Appendix 1. We use BO because it is one of the most successful approaches in optimizing calibration problems67. BO suffers from some limitations, and thus we modify this algorithm. First, BO has a high sensitivity to the type of acquisition function as an objective function and the nature of the gaussian procedure (GP), including kernel types and kernel hyperparameters68. Second, BO loses its efficiency in exploration and exploitation as the input dimensions increase. Third, the acquisition function is often difficult to optimize, and its performance depends on an optimizer to search the surface.
We develop a new parallel optimization scheme for the optimization part of the time-dependent parameter calibration framework, which addresses the BO method’s drawbacks. In this approach, instead of allocating all searching budget (the number of iteration) to one BO model, we divide the searching budget into several BO algorithms to search a solution space in parallel. To address BO’s sensitivity to the type of acquisition functions, each parallel BO has a specific type of acquisition function. Besides, each BO applies a particular kernel type and hyper-parameters for constructing the GP distribution. When BO models determine the next sample point by maximizing its equations function, they share their findings to help each other to acquire more knowledge about the solution space. The structural differences between BO and PBO are illustrated in Fig. 3. The PBO algorithm is defined in Algorithm 1.
We provide an additional explanation about the algorithm as follows.
-
Hyperparameters Parameter N should reflect the number of parallel BO models such that any increase in N will greatly expand the searching power to search through more enumeration. The maximum number of iterations T indicates the stop criteria of the algorithm. The maximum number of iterations can be replaced with a threshold so that we terminate the algorithm if the difference between the current value of the best new points and the incumbent solution is less than a predefined threshold.
-
Step 0 In this step, we set one type of acquisition function, kernel, and kernel parameters for constructing the GP of each BO model.
-
Step 1 GP updates the posterior probability, and constructs acquisition function in this step.
-
Step 2 We use the limited-memory quasi-Newton algorithm for bound-constrained optimization \(\text {L-BFGS-B}\)69,70 in Python to optimize the acquisition function with a set of random starting points. We run \(\text {L-BFGS-B}\) several times with different arbitrary starting points to improve the efficiency and overcome the difficulty of optimizing the acquisition function. We select the next sample point as the best optimal solution between all found solutions of \(\text {L-BFGS-B}\) algorithm with different starting points.
-
Step 3 The APSIM is run to evaluate a new parameter combination from each BO model in this step. Then, the BO models share their findings to expand their knowledge about posterior distribution such that the BO model applies GP in the next iteration to update its posterior distribution and construct the acquisition function.
-
Step 4 The incumbent solution is updated as the best combinations of parameters evaluated by APSIM in this step.
Case study
A case study with 25 locations from 1984 to 2019 (34 years) has been designed to illustrate how good the proposed time-dependent parameter estimation method is. We used the APSIM maize model71 (version 7.9), and in particular, the APSIM simulation files outlined by Archontoulis et al.14 that comprise simulation of shallow water tables and inhibition of root growth due to excess water stress18, and waterlogging functions72. We used APSIM because it is open-source software, advanced simulator of the cropping systems, and the model is commonly used in the US Corn Belt as well as worldwide73. We selected 5 counties (Logan (IL), Greene (IN), Keokuk (IA), Boone (IA), and Obrien (IA)) across the US Corn Belt and within each county 5 different soils, generating a range of variability in the weather and soil inputs (see Fig. 4). The APSIM simulation files were generated using pSIMS74, which is used to run APSIM across the US Corn Belt at a resolution of 5 arcminutes19,75. Each APSIM simulation file contained a soil profile, a weather file, the historical management, including N-fertilization, planting dates, plant population per year and state, and the simulated yield by year as well as a number of other soil-plant processes. Historical statewide plant density and planting dates per year were derived from USDA NASS quick stats76. County-level N fertilization rate was derived from a combined analysis of Cao et al.77 and USDA NASS76. The same observed yield for five locations at each county is used because the county-level yield data was acquired from the National Agricultural Statistics Service of the United States76. Since the management and the environment (weather \(\times \) soil) is already captured in the model, the only unknown is the cultivar parameters and how these change from year to year or from decade to decade. Climate information derived from a reanalysis weather database based on NASA Power79 and Iowa Environmental Mesonet80, both daily historical weather sources.
For each county, we used a K-medoids method81 as a clustering technique to select the five locations with the most distinct soil information. Figure 4 (right) shows the variability in soil organic carbon (in %) and plant available water (the difference between field capacity and wilting point in cm; depth 0–150 cm) across the 25 soils. We used adapted cultivars within a county in terms of relative maturities, given the 1.5-fold gradient in temperature from southeast to northwest (see Fig. 4).
Experiment settings
In this study, we considered the prediction of corn yield increase from 1985 to 2018 to evaluate the performance of the developed framework in estimating nine time-dependent parameters that can capture the observed historical yield increase: 8160 kg/ha to 12,679 kg/ha. We selected the following parameters summarized in Table 1 to optimize for two reasons: (1) their known sensitivity on yield and (2) evidence that these parameters have changed with plant breeding from 1980 to today51,52,53,54,55,56.
This case study was conducted for two purposes. First, to test the proposed optimization method in the proposed framework, we calibrate the time-dependent APSIM parameters to predict historical corn yield increase with manual calibration, the proposed framework using BO and PBO algorithms. In the manual calibration, the following six parameters were changed every 5-years in the model51: potential kernel number ear, grain growth rate, thermal time from silking to start grain filling, the radiation use efficiency, the transpiration use efficiency, and finally the grain nitrogen critical concentration. The direction and magnitude of changes were informed from literature studies51,52,53,54,55,56. Simultaneously, the results were used to benchmark new predictions using the crop model with tuned parameters from the proposed framework using BO and PBO algorithms. For the APSIM-BO and APSIM-PBO calibration, we follow a similar approach, but we allowed three additional parameters to change over time for a total of 9 parameters (Table 1). The 1985 to 2018 period was divided into 30 time-windows such that each main time-window contains five years. The proposed framework was implemented in Python 3.
Second, we benchmarked the yield prediction performance of the proposed parameter calibration method with machine learning models to illustrate the impact of the proposed time-dependent parameter estimation framework in enhancing the crop model’s performance. Nine machine learning models applied for this benchmark are summarized as follows:
-
Linear regression I model was implemented in Python using the Sklearn package82.
-
Linear regression II model was implemented in Python using the Sklearn package82.
-
Lasso model was trained using the Sklearn package82 in Python.
-
Ridge regression was implemented using the Sklearn package82 in Python.
-
Elastic net was trained using the Sklearn package82 in Python. Its hyper-parameter “alpha” was tuned using five-fold cross-validation (without data leakage).
-
Bayesian ridge was implemented in Python using the Sklearn package82. Hyper-parameters lambda (precision of the weights) and alpha (precision of the noise) were using five-fold cross-validation (without data leakage).
-
Gradient boosting was trained using the Sklearn package82 in Python, which was an efficient and scalable implementation of gradient boosting framework. Three hyper-parameters were tuned using five-fold cross-validation (without data leakage): “max_depth” (depth of tree), “n_estimators” (number of trees), and “max_features” (number of features).
-
Random forest was trained using the Sklearn package82 in Python, which was an ensemble of decision trees and trains with the bagging method. Three hyper-parameters were tuned using five-fold cross-validation: “n_estimators” (number of trees), “max_depth” (depth of tree), and “max_features” (number of features).
-
Deep neural network was implemented in Python using the Sklearn package82. Hidden layer sizes and the number of hidden layers as its hype-parameter were tuned using five-fold cross-validation.
For this test, we considered the last five years. To predict corn yield at test year \(Y+1\) with APSIM crop model, the historical data set, including all the environment (weather and soil), management, and crop information, and crop yield data from 1985 to year Y are applied to tune the parameters using the proposed framework. After fitting the autoregressive model to historical tuned parameters and estimating the parameter values for test year \(Y+1\), we ran APSIM with the estimated parameters at year \(Y+1\) to predict the corn yield.
To train the machine learning models, we create data set with soil profiles from SSURGO78, climate information from NASA Power79 and Iowa Environmental Mesonet80, and historical management information and yield performance data76 for all 293 counties of the states of Illinois, Indiana, and Iowa from 1985 to 2018. Weather data included four weekly average surface weather parameters (converted from daily data): precipitation (mm), solar radiation (MJ/m\(^2\)), maximum temperature (C\(^\circ \)), and minimum temperature (C\(^\circ \)) from weeks 13 (late March) to 52 (late December). Soil variables included dry bulk density (g cm\(^{-3}\)), clay percentage (%), soil pH, drained upper limit (mm mm\(^{-1}\)), saturated hydraulic conductivity (mm/day), drained lower limit (mm mm\(^{-1}\)), organic matter (%), sand percentage (%), and saturation (mm mm\(^{-1}\)) at nine different depths of soil: 0–5, 5–10, 10–15, 15–30, 30–45, 45–60, 60–80, 80–100, and 100–120 cm. Management variables included acres planted at the county level, weekly cumulative percentage of planted and harvested acreages, and nitrogen applications at the county level. Therefore, the total number of variables is 324 (including 160 for weather, 81 for soil, 83 for management).
The time variable included the year of observation to capture the changes in crop parameter values over the years. The created data set was divided into training and test data sets. The training data set contains the historical information from 1985 to year Y for all 293 counties, and the test data set is comprised of information of year \(Y+1\) for five considered counties. We used 293 counties so that the machine learning models have enough data set to train and figure out the relationship between input and response variables. We tuned the machine learning models’ hyper-parameters and trained them using the training data set, and then their performances were evaluated using the test data set. The first linear regression model (Linear regression I) fits the linear line on county-yield over time such that it represents a linear trend of county-yield vs time. However, other machine learning models intend to predict county-yield using weather, soil, management and time information. A feature selection step was performed to select a set of high-quality features before fitting the data set to machine learning models. We applied forward and backward stepwise selection by considering a 10-fold cross-validation schema to identify features.
Three criteria were used for comparing and evaluating the prediction performance of models: root mean square error (RMSE), which computes the difference between predicted and observed yield, relative RMSE (RRMSE), which measures the normalized difference between predicted and observed yield, and coefficient of determination (\(R^2\)), which calculates the proportion of the variance in the yield that is explained by independent variables.
The proposed framework was implemented in Python 3. We performed the experiments on a workstation with a 3.4 GHz CPU and 16 GB memory. The training was done by location such that calibration of these nine parameters at each time-window for one location took 8 CPU hours. The workstation’s CPU allows us to run 6 combinations of locations by time-windows at the same time. Hence, the total computation time for estimating nine parameters in all combinations of location by time-window (30 \(\times \) 25) was approximately 41 CPU days.
Results
Comparison with other parameter calibration methods
To compare with other parameter calibration methods, the median predicted yield of five locations in each county is used instead of the mean because it reduces the outliers’ effect on county-level yield. The result shows that the proposed framework with the PBO algorithm outperformed the BO and the manual approaches in terms of yield simulation from 1985 to 2018 (Table 2). Compared to the manual optimizer, the PBO optimizer decreased the RMSE by 47% at Logan and up to 65% at Obrien. The PBO optimizer has also achieved lower error than the BO optimizer, at least 4% at Greene and up to 16% at Boone. This performance can be attributed to using the proposed optimization algorithm in parameter tuning.
Figure 5 shows the time-series simulations. The median predicted yield of five locations at each county was considered as the predicted yield of counties because of county-level observed yield. Each sub-figure corresponds to one county in this figure, and the error bar illustrates our prediction interval at each year. According to the result, the consistency of the simulated yield performance across the corn belt’s spatial and temporal dimensions indicates the proposed framework’s effectiveness using PBO optimizer in the calibration of time-dependent parameters in APSIM. Another observation is that the proposed framework can determine the time-dependent parameters’ values to capture both the overall increasing trend of yield over simulation years and yield fluctuations from one year to another year.
Prediction accuracy comparison with other machine learning models
Prediction performance of five test years using the crop model with estimated cultivar parameters and machine learning models are reported in Table 3. These results indicate that the proposed framework with the PBO method outperformed other approaches for four out of five test years. It seems that close-to-optimal parameters’ values from the proposed framework help APSIM to simulate the yield more accurately and significantly reduce the crop yield prediction error.
The goodness of fit of the proposed framework with the different optimizer and random forest in crop yield prediction problem suggest that the prediction of APSIM with parameters from the proposed time-dependent parameter is close to the observation (the better distribution of the residuals derived from the distance between predicted and observed yield corresponding to line \(x=y\)), which decreased overall prediction bias.
Time-dependent parameter estimates
The linear trend of tuned parameters for five counties over time is illustrated in Fig. 6 with five dash-dot colored lines such that each color corresponds to a specific county. For each county, the average trend its locations is shown with the dash-dot line. The solid black lines and dashed black lines indicate average and standard deviation of parameters’ trends over five counties. The confidence interval indicates the range of variability across 25 environments. We did not find any correlation between optimized parameters and soil or weather variables to explain this range of variability. Probably, it is the result of a complex interaction. The slope of the linear trend of parameters is reported in Table 4.
The PBO estimated parameter values indicated that the RUE parameter increases linearly from 1985 to 2018 with an average slope of 0.001 g dm/mj. The transpiration efficiency parameter remained the same. Both grain component parameters (grain growth rate and potential kernel number per ear) increased over time with an average slope of 0.03 mg/grain/d and 1.41 kernel/ear, respectively. The leaf appearance rate decreased over time with an average slope of 0.05 deg day, indicating lower values for this parameter and thus faster leaf production. Among three thermal time parameters explored in this study, the emergence to a juvenile parameter that refers to the vegetative phase duration was decreased by 0.73 \(^{\circ }\)Cd per year, the flowering to the maturity parameter that refers to the reproductive phase duration was increased by 2.32 \(^{\circ }\)Cd per year, while the flowering to start seed parameter decreased by 0.33 \(^{\circ }\)Cd per year. Finally, the parameter related to grain protein (n_conc_crit_grain) was very slightly increased. Also, in terms of the estimated parameters’ variation range over the 25 environments, the largest variation was observed for the tt_emerg_to_endjuv parameter.
Discussion
This study developed and demonstrated a new time-dependent parameter estimation framework for the APSIM crop model that captures well historical yield increase in the US Corn Belt. To the best of our knowledge, this is the first study that develops a new automated framework to estimate time-dependent cultivar related parameters. Compared to other calibration methods (manual17 and BO39,45 methods) in the literature, our models can capture the trend of the crop model’s output by determining the temporal dependency of parameters. While the proposed framework was designed to enhance the crop model’s performance in simulating nature by calibrating its model parameters, we demonstrated its effectiveness by comparing its prediction accuracy of crop with that of machine learning models, which were specifically trained for prediction accuracy and not for parameter calibration. Results from our case study showed that the prediction accuracy of the calibrated crop model with the proposed framework was favorably comparable with that of machine learning models as well as what has been reported in the literature of crop yield prediction75,83,84,85,86.
Table 2 and Fig. 5 indicate that APSIM with estimated parameters of the proposed approach with the PBO optimizer was able to capture observed yield trend and yield fluctuations from one year to another year during the simulation period. The improved prediction performance was also more generalizable to the test data set, as demonstrated in Table 3. The proposed framework’s revealed trends of plant traits that have changed from 1985 to 2015 over the year due to crop improvement. Furthermore, results from this study can stimulate similar cultivar parameter calibration efforts in other crop models used for yield prediction and climate change assessments in the US Corn Belt and beyond. Our study offers an alternative way to the traditional field experimentation, era studies, to understand the reasons for the historical yield increase in the US Corn Belt, although ultimate validation of key parameters will still require field experiments.
Although we used the autoregressive model to capture the time-dependent trend of parameters over time in the proposed framework, we can use a linear regression model instead of an autoregressive model to capture the linear parameters’ trend from 1985 to 2018. More details about the autoregressive model’s performance on estimating weighted average parameters values in the proposed framework are provided in Appendix 2.
The estimated trend in RUE parameter increasing from 1985 to 2018 is consistent with era studies, in which cultivars from different eras were compared51,52. In these studies, the authors estimated an RUE increase of \(0.0057 \pm 0.002\) g dm/mj per year , which is close to our estimate. The PBO approach indicated that the transpiration efficiency coefficient did not change over time, however, as it is shown in Fig. 6f there was quite some variation among the 25 environments, with some environment to show an increase. Literature studies52,87,88 have indicated an increase in water use efficiency, a proxy for the transpiration use efficiency parameter. Perhaps the disagreement between our estimate and literature estimates comes from the soil that used (e.g. soil with or without watertable) and the overall simulation of the water balance.
Regarding the grain components, the PBO indicated that both the grain growth rate (a proxy for the harvest index) and the potential kernel number have increased over time, which in general follows literature experimental findings52,55,56,89,90,91,92. For example, the kernels per ear parameter in our study increased from 700 to 750, while the range of increase that has been reported in the literature using different era hybrids and different environments ranges55. The estimated increase in grain filling period by about 100 \(^{\circ }\)Cd (\(\sim 5\) days) over the years agrees with experimental reports53,90, satellite data analysis93,94, and NASA phenological data analysis95. We did not find relevant literature to benchmark our estimates for the other three phenological parameters (leaf appearance rate, vegetative stage duration, and thermal from flowering to start grain filling).
In contrast to previous studies54,56,92,96 that suggest a decline a grain protein (and therefore grain N concentration) with years of hybrid release, the estimated critical grain N concentration parameter by the PBO approach did not show a consistent decreasing trend over time (e.g., Figs. 5, 18, and 21 in Appendix 2). In 14 of the 25 cases (5 soils × 5 counties), the critical grain N concentration increased over time (linear slope is more than \(1 \times 10^{-5}\) g N/g biomass per year), in 6 out of the 25 cases decreased over time (linear slope is less than \(-1 \times 10^{-5}\) g N/g biomass per year) while in 5 out of the 25 cases remained unchanged over time (linear slope is between \(-1 \times 10^{-5}\) g N/g biomass per year and \(1 \times 10^{-5}\) g N/g biomass per year) (e.g., Figs. 8 and 17 in Appendix 2). As the only one out of nine calibrated parameters that showed an inconsistent trend with the literature (as well as manual APSIM calibration), this may be an artifact of “equifinality”, i.e., different combinations of parameter values can yield the same result20,97,98 in the PBO approach, which suggest the need for a) more data on grain N concentration to better estimate this parameter or b) development of alternative algorithms. There may also be inconsistency in its definition as the critical grain N concentration parameter in the crop model sets a target, which influences the final grain N concentration estimation, but it is not equivalent to the final grain N concentration that is measured in field experiments. The grain N concentration is a very complex trait. Many other processes influence its final value, including soil N mineralization, N leaching, denitrification, crop N uptake, and fertilizer input in addition to the critical grain N concentration parameter.
Despite its encouraging performance, the proposed approach has its limitations. While yield variation can be influenced by a large number of parameters in APSIM, we chose the nine parameters we considered most likely to vary in the Midwestern germplasm as free parameters. The possibility remains that there are alternative combinations of traits that underlie the yield changes between 1985 and 2018, which potentially explains why fitted values like critical grain N concentration are not consistent with experimental observations. Choosing alternative combinations of free parameters may suggest alternative explanations for yield changes.
A recent study99 reported that county-level yield data are insufficient to test crop model performance in extreme weather years (e.g., the year 2012). Our simulations did not provide evidence for this. This is because our model setup included shallow water tables that may have positive and negative impacts on productivity14, especially in extreme weather years such as 2012. We used the autoregressive model to capture the parameters’ trends; the modeling framework in Fig. 2 was based on this time-series model. To accommodate other trends, one can replace the autoregressive model with different time series models. Moreover, the PBO algorithm is computationally complex and sensitive to a number of parameters. In our case study, only nine parameters needed to be calibrated; when this number increases to many dozens, the computational time for the proposed framework could significantly increase.
Beyond parameter calibration for APSIM, the new approach can be used for a wide variety of studies that involve parameter calibration, especially those whose parameters demonstrate specific patterns over time. The modeling framework illustrated in Fig. 2 can be readily modified to accommodate other models. Future research should further explore the environmental and management parameters to make more biologically and agronomically insightful discoveries. Moreover, this framework can be used to optimize both short-term and long-term farm management strategies.
Conclusions
We developed and demonstrated a new framework to estimate time-dependent parameters’ values for crop models to improve performance with regards to capturing historical yield trends. We integrated the power of optimization, machine learning, and prior agronomic knowledge to develop a new automated parameter calibration algorithm. The developed PBO algorithm search for the best combination of parameters in high dimensional space to enhance the exploration and exploitation of the optimization step in the proposed framework. A computational case study using 25 environments in the US Corn Belt, each with 34 years of data, confirmed that the proposed framework outperformed the state-of-the-art calibration model (BO and manual calibration) and enhanced crop model performance in corn yield prediction accuracy so that it can outperform accurate machine learning models in the literature. Finally, this work produced insights about plant traits’ changes over the last 34 years, offering an alternative way to deeper understand historical yield increases in the US Corn Belt and inform future experimentation.
References
Heslot, N., Akdemir, D., Sorrells, M. E. & Jannink, J.-L. Integrating environmental covariates and crop modeling into the genomic selection framework to predict genotype by environment interactions. Theoretical and Applied Genetics 127, 463–480 (2014).
Schauberger, B. et al. Consistent negative response of US crops to high temperatures in observations and crop models. Nature Communications 8, 13931 (2017).
Holzworth, D. P. et al. APSIM-evolution towards a new generation of agricultural systems simulation. Environmental Modelling & Software 62, 327–350 (2014).
Basso, B., Liu, L. & Ritchie, J. T. A comprehensive review of the CERES-wheat,-maize and-rice models’ performances. In Advances in Agronomy, vol. 136, 27–132 (Elsevier, 2016).
Wang, X.-C., Li, J., Tahir, M. N. & Fang, X.-Y. Validation of the EPIC model and its utilization to research the sustainable recovery of soil desiccation after alfalfa (medicago sativa l.) by grain crop rotation system in the semi-humid region of the loess plateau. Agriculture, Ecosystems & Environment 161, 152–160 (2012).
Bassu, S., Asseng, S., Motzo, R. & Giunta, F. Optimising sowing date of durum wheat in a variable mediterranean environment. Field Crops Research 111, 109–118 (2009).
Eitzinger, J., Trnka, M., Hösch, J., Žalud, Z. & Dubrovskỳ, M. Comparison of CERES, WOFOST and SWAP models in simulating soil water content during growing season under different soil conditions. Ecological Modelling 171, 223–246 (2004).
Monsi, M. & Saeki, T. On the factor light in plant communities and its importance for matter production. Annals of Botany 95, 549 (2005).
Wallach, D. et al. How well do crop modeling groups predict wheat phenology, given calibration data from the target population?. BioRxiv 708578, (2020).
Zhao, G., Bryan, B. A. & Song, X. Sensitivity and uncertainty analysis of the APSIM-wheat model: Interactions between cultivar, environmental, and management parameters. Ecol. Model. 279, 1–11 (2014).
Togliatti, K., Archontoulis, S. V., Ranae, D., Puntel, L. & VanLoocke, A. How does inclusion of weather forecasting impact in-season crop model predictions? Field Crops Res. 214, 261–272, https://doi.org/10.1016/j.fcr.2017.09.008 (2017).
Lamsal, A. et al. Efficient crop model parameter estimation and site characterization using large breeding trial data sets. Agric. Syst. 157, 170–184 (2017).
Anwar, M. R. et al. Climate change impacts on phenology and yields of five broadacre crops at four climatologically distinct locations in Australia. Agric. Syst. 132, 133–144 (2015).
Archontoulis, S. V. et al. Predicting crop yields and soil-plant nitrogen dynamics in the US corn belt. Crop Sci. 60, 721–738 (2020).
Dumont, B. et al. A comparison of within-season yield prediction algorithms based on crop model behaviour analysis. Agric. For. Meteorol. 204, 10–21 (2015).
Seidel, S. J., Palosuo, T., Thorburn, P. & Wallach, D. Towards improved calibration of crop models-where are we now and where should we go?. Eur. J. Agron. 94, 25–35 (2018).
Ahmed, M. et al. Calibration and validation of APSIM-wheat and CERES-wheat for spring wheat under rainfed conditions: Models evaluation and application. Comput. Electron. Agric. 123, 384–401 (2016).
Ebrahimi-Mollabashi, E. et al. Enhancing APSIM to simulate excessive moisture effects on root growth. Field Crops Res. 236, 58–67 (2019).
Baum, M., Archontoulis, S. & Licht, M. Planting date, hybrid maturity, and weather effects on maize yield and crop stage. Agron. J. 111, 303–313 (2019).
Archontoulis, S. V., Miguez, F. E. & Moore, K. J. 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 (2014).
Saltelli, A. et al. Global Sensitivity Analysis: The Primer (Wiley, New York, 2008).
Turanyi, T. & Rabitz, H. Local methods. In Sensitivity Analysis 81–99 (Wiley, Chister, 2000).
Cukier, R., Schaibly, J. & Shuler, K. E. Study of the sensitivity of coupled reaction systems to uncertainties in rate coefficients. iii. analysis of the approximations. J. Chem. Phys. 63, 1140–1149 (1975).
Saltelli, A., Tarantola, S. & Chan, K.-S. A quantitative model-independent method for global sensitivity analysis of model output. Technometrics 41, 39–56 (1999).
Kuczera, G. & Mroczkowski, M. Assessment of hydrologic parameter uncertainty and the worth of multiresponse data. Water Resour. Res. 34, 1481–1489 (1998).
Sumner, N. R., Fleming, P. M. & Bates, B. C. Calibration of a modified SFB model for twenty-five Australian catchments using simulated annealing. J. Hydrol. 197, 166–188 (1997).
Duan, Q., Sorooshian, S. & Gupta, V. Effective and efficient global optimization for conceptual rainfall-runoff models. Water Resour. Res. 28, 1015–1031 (1992).
Ramirez-Villegas, J., Koehler, A.-K. & Challinor, A. J. Assessing uncertainty and complexity in regional-scale crop model simulations. Eur. J. Agron. 88, 84–95 (2017).
Helton, J. C., Davis, F. & Johnson, J. D. A comparison of uncertainty and sensitivity analysis results obtained with random and latin hypercube sampling. Reliabil. Eng. Syst. Saf. 89, 305–330 (2005).
Ansarifar, J. & Wang, L. New algorithms for detecting multi-effect and multi-way epistatic interactions. Bioinformatics 35, 5078–5085 (2019).
Young, P. C. Stochastic, dynamic modelling and signal processing: Time variable and state dependent parameter estimation. In Nonlinear and Nonstationary Signal Processing 74–114 (Cambridge University Press, Cambridge, 2000).
Sexton, J., Everingham, Y. & Inman-Bamber, G. A theoretical and real world evaluation of two Bayesian techniques for the calibration of variety parameters in a sugarcane crop model. Environ. Model. Softw. 83, 126–142 (2016).
Dumont, B. et al. Parameter identification of the STICS crop model, using an accelerated formal MCMC approach. Environ. Model. Softw. 52, 121–135 (2014).
Iorgulescu, I., Beven, K. & Musy, A. Data-based modelling of runoff and chemical tracer concentrations in the Haute–Mentue research catchment (Switzerland). Hydrol. Process. 19, 2557–2573 (2005).
Makowski, D., Wallach, D. & Tremblay, M. Using a Bayesian approach to parameter estimation; comparison of the GLUE and MCMC methods. Agronomie 22, 191–203 (2002).
Georgakakos, K. P., Seo, D.-J., Gupta, H., Schaake, J. & Butts, M. B. Towards the characterization of streamflow simulation uncertainty through multimodel ensembles. J. Hydrol. 298, 222–241 (2004).
Mara, T. A., Fajraoui, N., Younes, A. & Delay, F. Inversion and uncertainty of highly parameterized models in a Bayesian framework by sampling the maximal conditional posterior distribution of parameters. Adv. Water Resour. 76, 1–10 (2015).
Alderman, P. D., Boote, K. J., Jones, J. W. & Bhatia, V. S. Adapting the CSM-CROPGRO model for pigeonpea using sequential parameter estimation. Field Crops Res. 181, 1–15 (2015).
Laloy, E., Fasbender, D. & Bielders, C. Parameter optimization and uncertainty analysis for plot-scale continuous modeling of runoff using a formal Bayesian approach. J. Hydrol. 380, 82–93 (2010).
Zadeh, F. K., Nossent, J., Woldegiorgis, B. T., Bauwens, W. & van Griensven, A. Impact of measurement error and limited data frequency on parameter estimation and uncertainty quantification. Environ. Model. Softw. 118, 35–47 (2019).
Kuczera, G., Kavetski, D., Franks, S. & Thyer, M. Towards a Bayesian total error analysis of conceptual rainfall-runoff models: Characterising model error using storm-dependent parameters. J. Hydrol. 331, 161–177 (2006).
Beven, K. & Binley, A. The future of distributed models: Model calibration and uncertainty prediction. Hydrol. Process. 6, 279–298 (1992).
Wallach, D., Keussayan, N., Brun, F., Lacroix, B. & Bergez, J.-E. Assessing the uncertainty when using a model to compare irrigation strategies. Agron. J. 104, 1274–1283 (2012).
Ceglar, A., Črepinšek, Z., Kajfež-Bogataj, L. & Pogačar, T. The simulation of phenological development in dynamic crop model: The Bayesian comparison of different methods. Agric. For. Meteorol. 151, 101–115 (2011).
Berkenkamp, F., Krause, A. & Schoellig, A. P. Bayesian optimization with safety constraints: safe and automatic parameter tuning in robotics. arXiv:1602.04450 (2016).
Liu, Y., Li, Y., Huang, G., Zhang, J. & Fan, Y. A Bayesian-based multilevel factorial analysis method for analyzing parameter uncertainty of hydrological model. J. Hydrol. 553, 750–762 (2017).
Song, X. et al. Time-dependent sensitivity of a process-based ecological model. Ecol. Model. 265, 114–123 (2013).
Wang, F., Mladenoff, D. J., Forrester, J. A., Keough, C. & Parton, W. J. Global sensitivity analysis of a modified century model for simulating impacts of harvesting fine woody biomass for bioenergy. Ecol. Model. 259, 16–23 (2013).
Gallagher, M. & Doherty, J. Parameter estimation and uncertainty analysis for a watershed model. Environ. Model. Softw. 22, 1000–1020 (2007).
Nissanka, S. P. et al. Calibration of the phenology sub-model of APSIM-oryza: Going beyond goodness of fit. Environ. Model. Softw. 70, 128–137 (2015).
Reyes Ponce, A. The physiological basis of genetic improvement in maize (zea mays l) yield in the US corn belt (The University of Queensland, Thesis for Doctor of Philosophy, 2016).
Curin, F., Severini, A. D., González, F. G. & Otegui, M. E. Water and radiation use efficiencies in maize: Breeding effects on single-cross argentine hybrids released between 1980 and 2012. Field Crops Res. 246, 107683 (2020).
Duvick, D. N. The contribution of breeding to yield advances in maize (Zea mays L.). Adv. Agron. 86, 83–145 (2005).
Duvick, D. N. & Cassman, K. G. Post-green revolution trends in yield potential of temperate maize in the north-central United States. Crop Sci. 39, 1622–1630 (1999).
DeBruin, J. L., Schussler, J. R., Mo, H. & Cooper, M. Grain yield and nitrogen accumulation in maize hybrids released during 1934 to 2013 in the US midwest. Crop Sci. 57, 1431–1446 (2017).
Haegele, J. W., Cook, K. A., Nichols, D. M. & Below, F. E. Changes in nitrogen use traits associated with genetic improvement for grain yield of maize hybrids released in different decades. Crop Sci. 53, 1256–1268 (2013).
Ansarifar, J., Akhavizadegan, F. & Wang, L. Performance prediction of crosses in plant breeding through genotype by environment interactions. Sci. Rep. 10, 1–11 (2020).
Moeinizade, S., Hu, G., Wang, L. & Schnable, P. S. Optimizing selection and mating in genomic selection with a look-ahead approach: An operations research framework. G3 9, 2123–2133 (2019).
Moeinizade, S., Wellner, M., Hu, G. & Wang, L. Complementarity-based selection strategy for genomic selection. Crop Sci. 60, 149–156 (2020).
Moeinizade, S., Kusmec, A., Hu, G., Wang, L. & Schnable, P. S. Multi-trait genomic selection methods for crop improvement. Genetics 215, 931–945 (2020).
Guo, J., Zhao, J., Wu, D., Mu, J. & Xu, Y. Attribution of maize yield increase in China to climate change and technological advancement between 1980 and 2010. J. Meteorol. Res. 28, 1168–1181 (2014).
Merlos, F. A. et al. Potential for crop production increase in Argentina through closure of existing yield gaps. Field Crops Res. 184, 145–154 (2015).
Akaike, H. Fitting autoregressive models for prediction. Ann. Inst. Stat. Math. 21, 243–247 (1969).
Dickey, D. A. & Fuller, W. A. Distribution of the estimators for autoregressive time series with a unit root. J. Am. Stat. Associ. 74, 427–431 (1979).
Ramsey, F. L. et al. Characterization of the partial autocorrelation function. Ann. Stat. 2, 1296–1301 (1974).
Mockus, J. Bayesian Approach to Global Optimization: Theory and Applications Vol. 37 (Springer Science & Business Media, New York, 2012).
Jones, D. R. A taxonomy of global optimization methods based on response surfaces. J. Glob. Optim. 21, 345–383 (2001).
Snoek, J., Larochelle, H. & Adams, R. P. Practical Bayesian optimization of machine learning algorithms. In Advances in Neural Information Processing Systems 2951–2959 (University of Toronto, Toronto, 2012).
Zhu, C., Byrd, R. H., Lu, P. & Nocedal, J. Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization. ACM Trans. Math. Softw. 23, 550–560 (1997).
Byrd, R. H., Lu, P., Nocedal, J. & Zhu, C. A limited memory algorithm for bound constrained optimization. SIAM J. Sci. Comput. 16, 1190–1208 (1995).
Keating, B. A. et al. An overview of APSIM, a model designed for farming systems simulation. Eur. J. Agron. 18, 267–288 (2003).
Pasley, H. R., Huber, I., Castellano, M. J. & Archontoulis, S. V. Modeling flood-induced stress in soybeans. Front. Plant Sci. 11, 62 (2020).
Keating, B. A. & Thorburn, P. J. Modelling crops and cropping systems-evolving purpose, practice and prospects. Eur. J. Agron. 100, 163–176 (2018).
Elliott, J. et al. The parallel system for integrating impact models and sectors (pSIMS). Environ. Model. Softw. 62, 509–516 (2014).
Shahhosseini, M., Hu, G., Archontoulis, S. V. & Huber, I. Coupling machine learning and crop modeling improves crop yield prediction in the US corn belt. arXiv:2008.04060 (2020).
Agricultural Statistics Service, N. https://quickstats.nass.usda.gov.
Cao, P. et al. Historical nitrogen fertilizer use in agricultural ecosystems of the contiguous united states during 1850–2015: application rate, timing, and fertilizer types. Earth Syst. Sci. Data Discuss. 10, 969 (2018).
Soil Survey Geographic Database, G. https://websoilsurvey.nrcs.usda.gov.
Prediction of Worldwide Energy Resource, N. https://power.larc.nasa.gov.
Environmental Mesonet, I. https://mesonet.agron.iastate.edu.
Park, H.-S. & Jun, C.-H. A simple and fast algorithm for K-medoids clustering. Expert Syst. Appl. 36, 3336–3341 (2009).
Pedregosa, F. et al. Scikit-learn: Machine learning in python. J. Mach. Learn. Res. 12, 2825–2830 (2011).
You, J., Li, X., Low, M., Lobell, D. & Ermon, S. Deep Gaussian process for crop yield prediction based on remote sensing data. In Thirty-First AAAI Conference on Artificial Intelligence (2017).
Marko, O., Brdar, S., Panic, M., Lugonja, P. & Crnojevic, V. Soybean varieties portfolio optimisation based on yield prediction. Comput. Electron. Agric. 127, 467–474 (2016).
Shahhosseini, M., Martinez-Feria, R. A., Hu, G. & Archontoulis, S. V. Maize yield and nitrate loss prediction with machine learning algorithms. Environ. Res. Lett. 14, 124026 (2019).
Romero, J. R. et al. Using classification algorithms for predicting durum wheat yield in the province of Buenos Aires. Comput. Electron. Agric. 96, 173–179 (2013).
Reyes, A. et al. Soil water capture trends over 50 years of single-cross maize (Zea mays L.) breeding in the US corn-belt. J. Exp. Bot. 66, 7339–7346 (2015).
Wu, Y. et al. Evolution of roots to improve water and nitrogen use efficiency in maize elite inbred lines released during different decades in China. Agric. Water Manag. 216, 44–59 (2019).
Lorenz, A., Gustafson, T., Coors, J. & De Leon, N. Breeding maize for a bioeconomy: A literature survey examining harvest index and stover yield and their relationship to grain yield. Crop Sci. 50, 1–12 (2010).
Chen, K. et al. Genetic improvement in density and nitrogen stress tolerance traits over 38 years of commercial maize hybrid release. Field Crops Res. 196, 438–451 (2016).
Campos, H. et al. Changes in drought tolerance in maize associated with fifty years of breeding for yield in the US corn belt. Maydica 51, 369 (2006).
Mueller, S. M., Messina, C. D. & Vyn, T. J. Simultaneous gains in grain yield and nitrogen efficiency over 70 years of maize genetic improvement. Sci. Rep. 9, 1–8 (2019).
Zhu, P. et al. The important but weakening maize yield benefit of grain filling prolongation in the US midwest. Glob. Change Biol. 24, 4718–4730 (2018).
Butler, E. E., Mueller, N. D. & Huybers, P. Peculiarly pleasant weather for US maize. Proc. Natl. Acad. Sci. USA 115, 11935–11940 (2018).
Sacks, W. J. & Kucharik, C. J. Crop management and phenology trends in the US corn belt: Impacts on yields, evapotranspiration and energy balance. Agric. For. Meteorol. 151, 882–894 (2011).
Scott, M. P., Edwards, J. W., Bell, C., Schussler, J. & Smith, J. Grain composition and amino acid content in maize cultivars representing 80 years of commercial maize varieties. Maydica 51, 417 (2006).
Beven, K. & Freer, J. Equifinality, data assimilation, and uncertainty estimation in mechanistic modelling of complex environmental systems using the glue methodology. J. Hydrol. 249, 11–29 (2001).
Jones, J. W. et al. Estimating DSSAT cropping system cultivar-specific parameters using Bayesian techniques. In: Methods of Introducing System Models Into Agricultural Research 365–394 (2011).
Guarin, J. R., Asseng, S., Martre, P. & Bliznyuk, N. Testing a crop model with extreme low yields from historical district records. Field Crops Res. 249, 107269 (2020).
Funding
This work was partially supported by the National Science Foundation under the LEAP HI and GOALI programs (Grant number 1830478) and under the EAGER program (Grant number 1842097) and by the Plant Sciences Institute at Iowa State University.
Author information
Authors and Affiliations
Contributions
F.A., J.A., L.W, and S.V. designed the research questions. F.A., J.A., I.H., and S.V. prepared and cleaned the database. F.A., and J.A. performed the experiment, statistical analysis, and analyzed the data set. F.A., and J.A. designed and implemented a new algorithm. F.A. created the figures. F.A., J.A, L.W., and S.V. interpreted experiment results. F.A. and J.A. wrote the manuscript. L.W. and S.V. thoroughly reviewed the manuscript and provided feedback. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Akhavizadegan, F., Ansarifar, J., Wang, L. et al. A time-dependent parameter estimation framework for crop modeling. Sci Rep 11, 11437 (2021). https://doi.org/10.1038/s41598-021-90835-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-90835-x
This article is cited by
-
How to build a crop model. A review
Agronomy for Sustainable Development (2023)
-
An interaction regression model for crop yield prediction
Scientific Reports (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.