Next Article in Journal
Cytochrome c Deficiency Differentially Affects the In Vivo Mitochondrial Electron Partitioning and Primary Metabolism Depending on the Photoperiod
Next Article in Special Issue
Modeling of Cowpea (Vigna unguiculata) Yield and Control Insecticide Exposure in a Semi-Arid Region
Previous Article in Journal
Both Allene Oxide Synthases Genes Are Involved in the Biosynthesis of Herbivore-Induced Jasmonic Acid and Herbivore Resistance in Rice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of Phenology Algorithms of the Spring Barley Model in APSIM 7.9 and APSIM Next Generation: A Case Study for High Latitudes

1
Department of Agricultural Research for Northern Sweden, Swedish University of Agricultural Science, 90183 Umeå, Sweden
2
Department of Crop Production and Ecology, Swedish University of Agricultural Science, 705007 Uppsala, Sweden
3
Natural Resources Institute Finland (Luke), FI-00790 Helsinki, Finland
4
CSIRO Agriculture and Food, Canberra 2601, Australia
5
The New Zealand Institute for Plant and Food Research Limited, Private Bag 4704, Christchurch 8140, New Zealand
*
Author to whom correspondence should be addressed.
Plants 2021, 10(3), 443; https://doi.org/10.3390/plants10030443
Submission received: 27 January 2021 / Revised: 18 February 2021 / Accepted: 22 February 2021 / Published: 26 February 2021
(This article belongs to the Special Issue Agroecosystem Modeling)

Abstract

:
Phenology algorithms in crop growth models have inevitable systematic errors and uncertainties. In this study, the phenology simulation algorithms in APSIM classical (APSIM 7.9) and APSIM next generation (APSIM-NG) were compared for spring barley models at high latitudes. Phenological data of twelve spring barley varieties were used for the 2014–2018 cropping seasons from northern Sweden and Finland. A factorial-based calibration approach provided within APSIM-NG was performed to calibrate both models. The models have different mechanisms to simulate days to anthesis. The calibration was performed separately for days to anthesis and physiological maturity, and evaluations for the calibrations were done with independent datasets. The calibration performance for both growth stages of APSIM-NG was better compared to APSIM 7.9. However, in the evaluation, APSIM-NG showed an inclination to overestimate days to physiological maturity. The differences between the models are possibly due to slower thermal time accumulation mechanism, with higher cardinal temperatures in APSIM-NG. For a robust phenology prediction at high latitudes with APSIM-NG, more research on the conception of thermal time computation and implementation is suggested.

1. Introduction

Process-based crop models simulate dynamic and complex interactions between environment, genotype, and management factors. Various algorithms and parameters in crop models simulate different plant and soil processes on several interactions and linkages. The processes related to phenology, dry matter accumulation and partitioning, and soil hydrology and chemistry are simulated with the various algorithms and parameters. Model algorithms and parameters are still simplifications of real systems [1], which makes crop models contain unavoidable systematic errors.
Additionally, at the user level, the quality of input data [2] and the choice of model parameterization approach [3] create further uncertainties. Since crop models are being extensively applied to a wide range of agricultural research questions and hypothesis testing, such as assessments of climate change effects [4], decision making and planning [5], farmer advisory [6], crop–livestock systems [7,8], agronomic management [9], physiological mechanisms and traits [10,11], linking phenotype to genotype [12,13], and plant breeding [14,15], it is critical to find ways to reduce the uncertainties.
Most of the commonly and widely used crop models (e.g., DSSAT, APSIM, CROPSYST, EPIC, STICS, WOFOST, DAISY, ORYZA, GLAM, and INFOCROP) were constructed at least two decades ago using crop data obtained under controlled growing conditions and limited datasets [16]. Since the development of the models, improving their prediction accuracy has been the main emphasis [17,18] to better adapt them for new environments and future research questions. Suggestions for identifying and improving the prediction accuracy of various processes have been an important part of several recent crop model ensembles and inter-comparison studies [19,20,21,22,23].
More specifically, studies have highlighted that robust phenological submodels are crucial tools to improve the accuracy of crop models [24,25], given the importance of the timing of anthesis and crop duration in crop yield determination [26,27]. Such improvements have great importance and can increase confidence and reliability when simulating the response of plants, for example, to climate change [28] and temperature [24], to better plan agronomic activities to maximize crop yield while reducing the risks. Interestingly five different phenological models/modules (i.e., CERES-Rice, ORYZA2000, RCM, Beta Model, SIMRIW) were compared and reported to simulate varied plant responses to climate change and variability in different regions [29]. The variations were more evident when the temperature in the regions was above or below the optimum temperature, as defined for optimum growth in the models. In crop models, thermal time accumulation is sometimes implemented with piecewise non-linear response functions with temperature-limit values called cardinal temperatures. With different cardinal temperatures and response functions used to compute thermal time for the targeted developmental stages, ORYZA2000 and CERES-Rice simulated varied phenology [24].
Various crop growth processes implemented in crop models, such as dry matter accumulation and partitioning, leaf number, and leaf area dynamics, are strongly linked with phenological development [25,30]. Rötter et al. [31] found incongruities in the predictions of crop phenology and yield. The study compared and applied nine crop models at seven different locations in Northern and Central Europe and reported that models such as WOFOST and STICS were more accurate in simulating phenology; however, they were not necessarily the best predictors of the grain yields. This can be associated with the fact that a small error in phenology can have a large influence on the yield, especially if errors occur during sensitive growth stages, such as fertilization. Thus, to increase the overall accuracy of crop models, it is critical to better synchronize phenological simulations with other processes in the model.
APSIM is currently one of the most widely used crop and farming systems models globally [32]. With the increasing challenges and demand for agricultural modelling, the APSIM initiative and its predecessor The Agricultural Production Systems Research Unit (APSRU) have been constantly improving and building APSIM model components since its release in 1990, in order to improve its prediction accuracy and increase its applicability to a wider range of farming systems [33]. An example of this effort is the implementation of a new phenology simulating mechanism in the recently released APSIM next-generation (APSIM-NG) models compared to the one used in the classical versions of APSIM (APSIM 7.x).
Spring barley is ranked fourth in world cereal production and is particularly important in boreal regions with short growing seasons, such as at high latitudes [34]. Although there have been several studies on barley using APSIM and other crop models (e.g., [21,31]), achieving a robust calibration is still a challenging task when less detailed experimental datasets are available [3]. The availability of good datasets generated at high latitudes are particularly poor in terms of data on critical growth stages, which limits the robustness of calibration and makes models prone to having poor prediction accuracy. For the same reason, the strengths and weaknesses of the APSIM barley model still need to be critically assessed at high latitudes for further application as a research tool. Therefore, the data on days to anthesis and physiological maturity on recently bred barley varieties were used in the study to evaluate the phenology simulation mechanisms of APSIM 7.9 and APSIM-NG. Focusing on one aspect of the simulating mechanism (in this case, phenology) allows for a more systematic assessment of the contributing factors to the uncertainty and the targets for improvement [31]. The relevance of assessing the models using current varieties, which were developed considering the short and risky cropping seasons experienced at high latitudes, particularly for phenology-simulating mechanisms, is significant and applicable in terms of their potential use in the region [35,36].
This study aimed to evaluate the phenological algorithms of the APSIM classical version 7.9 (hereafter APSIM 7.9) and the APSIM-NG (version APSIM2019.11.27.4417) models of barley using experimental data for recently developed barley varieties from high latitudes (north and central Sweden and central Finland). The specific objectives of the study were to 1) calibrate and evaluate phenology algorithms of APSIM 7.9 and APSIM-NG barley models using a smaller and larger dataset on days to anthesis and physiological maturity, and 2) identify potential sources of error in the algorithms and suggest means for future improvement.

2. Results

2.1. Crop Phenology Overview

All the varieties, on average, needed more days to reach physiological maturity in the colder cropping season in 2017 than in the warm and dry cropping season of 2018 (Table 1). On average, across the locations, varieties reached physiological maturity 33 days earlier in 2018 than 2017. At Röbäcksdalen in 2017, the variation in days to anthesis among the varieties was less than the variation in days to physiological maturity. Compared with 2017, the variation and duration of days to anthesis and physiological maturity were lower in 2018. The soil moisture and phenology data at Röbäcksdalen suggest that barley development was rapid as the warm and dry conditions intensified 40 days after sowing in 2018 (Figure 1 and Table 1. The years 2014 and 2018 were similar with respect to days to physiological maturity, whereas the years 2015 and 2016 were intermediate in terms of maturity in 2017 and 2018. On average, barley matured latest at Ås and earliest at Offer and Ruukki.

2.2. Calibration and Evaluation of APSIM 7.9

For Calibration 1_AN, RMSE was between 0 and 1.3 d, depending on the variety (Table 2). These RMSEs were achieved with photop_sens between 0 and 1, vern_sens, 0 and 0.5; tt_end_of_juvenile, 200 or 300; tt_floral_initiation, 300 or 320. For Calibration 2_PM, when the object variable was days to physiological maturity, the RMSE range among varieties was between 2.2 and 8.6 d, a wider range than for Calibration 1_AN. The wider range of RMSE was because of the different combinations of the parameters, a more variable photop_sens (0–6); the same range of vern_sens, the same range of tt_end_of_juvenile but with a few entries of 250, depending on the variety, same range of tt_floral_initiation, and tt_start_grain_fill between 500 and 625. It should be noted that different combinations of the parameters resulted in the same RMSEs, particularly for Calibration 1_AN, which is a common phenomenon, termed equifinality, when a large number of combinations are tested [37].
For Calibration 3_PM, RMSE for different varieties was between 7.5 and 11 days (Table 2). Although the dataset for this calibration was larger than for Calibration 2_PM, the range of RMSE was broader. Compared to the previous two calibrations, here, photop_sens was between 0 and 3; vern_sens, 0, tt_end_of_juvenile, 250–450 with more higher values of 450; tt_floral_initiation, 240–320, and tt_start_grain_fill, 300 or 450.
Based on the combination of parameters that resulted in the best possible RMSEs in the above calibrations, it can be inferred that the combinations were different because of the target of the growth stage and the size of the dataset used for calibration.
Parameters from Calibration 2_PM in the Evaluation 2_PM overestimated the days to maturity particularly for environments of early maturity and underestimated for late maturing environments (Figure 1A). However, in the Evaluation 3_PM for the parameters calibrated under Calibration 3_PM, the performance of the model was better, with lower RMSE terms and better r2 (Figure 2B and Table 3).

2.3. Calibration and Evaluation of APSIM-NG

For Calibration 1_AN, APSIM-NG showed very good calibration performance (RMSE zero days) for all varieties (Table 4). This RMSE was achieved with the following parameter values: PpSensitivity, zero, VrnSensitivity, 0.5–1.5; BasePhyll, 55 or 60; MinimumLeafNumber, 6; EarlyReproductivePpSensitivity, zero.
Similarly, for Calibration 2_PM, RMSE was zero days, indicating a better performance for APSIM-NG than APSIM 7.9 for the same calibration using the same dataset. This RMSE was achieved with the following parameter values: PpSensitivity equal to zero, VrnSensitivity, 0.5–1.5; BasePhyll, 50–60; MinimumLeafNumber, 5–6.5; GrainFill, 450–625; and EarlyReproductivePpSensitivity, zero. Similar to APSIM 7.9 different combinations of the parameters resulted in the same RMSEs under Calibration 1_AN. Therefore, the parameter combinations in Table 4 were selected with the same process as for APSIM 7.9.
For Calibration 3_PM, RMSE was between 4.6 and 7.5 days, depending on the variety (Table 4). Compared to the previous two calibrations, here the PpSensitivity was 0; there was a slightly wider range of VrnSensitivity (0–3) and BasePhyllo (50–65); the MinimumLeafNumber was 5–6.5, and the values for GrainFill were lower and narrower, 400–575; and EarlyReproductivePpSensitivity was 0.
As with APSIM 7.9, the best combination of the selected parameters for APSIM-NG was different depending on the target of the growth stage and the size of the dataset used for calibration.
Evaluation 2_PM showed that days to physiological maturity was consistently overestimated (Figure 2A). Twenty evaluation simulations did not reach physiological maturity and stopped at the end date of 30 November. The days after sowing of those incomplete simulations, together with the observed data, are presented in Figure 2A. Based on the completed simulations, the intercept was closer to zero compared to the same evaluation of APSIM 7.9 (Table 3). The consistent overestimation for the completed simulations resulted in a higher RMSE than for APSIM 7.9. However, with similar RMSEsys, and higher RMSEnos, the overall correlation (r2) was higher.
When the APSIM-NG calibration was performed with a larger dataset than in the evaluation (Evaluation 3_PM) RMSE was still higher; however, RMSEsys and RMSEnos were lower compared with the same evaluation and dataset for APSIM 7.9. Intercept and slope showed the inclination of overestimation; however, this was less obvious than for Evaluation 2_PM. The overall relationship (r2) of simulated values with the observed ones was lower than for APSIM 7.9 (Figure 3B and Table 3). Similar to APSIM 7.9, RMSE, RMSEsys, and RMSEnos under Evaluation 3_PM were lower than for Evaluation 2_PM.

3. Discussion

3.1. Performance of Phenology Simulation Algorithms in the Calibration and Evaluation of APSIM 7.9 and APSIM-NG

Although the calibration of the parameters for simulation of days to physiological maturity for APSIM-NG was better when based on the RMSE, in both evaluations, slopes < 1 suggested that APSIM-NG with these calibrations was inclined to overestimate the days to physiological maturity, particularly for later maturing varieties. This was affirmed by the consistent overestimation of days to physiological maturity with higher RMSE, and a more variable range of r2 for APSIM-NG (Table 3). In contrast, APSIM 7.9 overestimated the duration of the varieties when they matured earlier and underestimated them when they matured later. This indicated that the thermal time accumulation for different growth stages in APSIM-NG, determined by the cardinal temperatures applied, might be more sensitive to colder days in the late crop-maturing phase.
Similarly, for the calibration of anthesis, a lower RMSE for APSIM-NG than for APSIM 7.9 indicated that APSIM-NG’s phenology algorithm and the parameters for simulating days to anthesis were more adjustable to capture the responses of varieties to warm and dry (2018) and wet and cold (2017) cropping seasons. However, neither of the models was evaluated for days to anthesis due to lack of data.

3.2. Differences between APSIM 7.9 and APSIM-NG Algorithms for Simulating Phenology

For broader reliability of the model prediction, a robust calibration is indispensable for building confidence when predicting the targeted variables in different environments [38,39,40,41]. In this study, three calibrations with the targets of different growth stages, using different data sizes, provided a more robust scrutinization of phenology-regulating parameters. The parameters in APSIM-NG simulated days to physiological maturity less accurately in the evaluations compared with APSIM 7.9. As mentioned above, the reasons for this could be linked with cardinal temperatures (which were not calibrated) and their usage in the thermal time computations. In addition, the implementation of vernalization and photoperiod effects in the algorithms could be equally important.
The implementation of Kirby’s Framework [42,43] in APSIM-NG to capture the effects of vernalisation and photoperiod in sync with the Haun stage might provide better synchronicity and continuity in processes such as leaf emergence, floral primordia initiation, and anthesis, which could be the reason for the better calibration of APSIM-NG.
Thermal time computation in both APSIM-NG and APSIM7.9 models is achieved with the same approach (mean of maximum and minimum temperature during the day, and from which the base temperature is subtracted), but with different cardinal temperatures. Cardinal temperatures in both models, like other crop models, are crop-specific parameters. The new phenology simulation mechanism in APSIM-NG utilizes higher cardinal temperatures compared to APSIM 7.9. These parameters were not calibrated in this study, for two reasons: (1) the complete dataset on both growth stages was only for two seasons and one location, which we think were not enough to manipulate crop-specific parameters; (2) we aimed to test the phenology algorithms and parameters in their original form.
Apart from the usage of different cardinal temperatures, another difference that can lead to different responses in these models is employing the crown temperature for the computation of thermal time. In APSIM 7.9, air temperature is first converted into crown temperature, and then the thermal time computation is made. However, in APSIM-NG, instead of using crown temperature for the computations, higher base temperature (Tb), optimum temperature (Topt), and maximum temperature (Tmax) are used as proxies for crown temperature. In APSIM-NG, the response function of thermal time computation has two slopes, with three air cardinal temperatures: Tb = 4 °C, Topt = 26 °C, (a positive slope from Tb to Topt) and maximum Tm = 37 °C, with a negative slope from Topt to Tm that results in 0, 22 and 0 °Cd thermal time for the respective cardinal temperatures. In APSIM 7.9, the air cardinal temperatures are, Tb = 0 °C, Topt = 26 °C and Tm = 34 °C, with the corresponding thermal times as 0 °Cd, 26 °Cd, and 0 °Cd, with the same response function as in APSIM-NG. This suggests that APSIM-NG has a slower thermal time accumulation mechanism compared with APSIM 7.9, which resulted in the overestimation of later maturing varieties when the temperature was relatively lower in the later crop maturing phase. Due to the slower thermal time accumulation, twenty simulations did not complete physiological maturity. For APSIM 7.9, the accumulation was faster, hence it underestimated maturity for some of the late-maturing varieties.
Similar to other models, such as DSSAT, STICS, GLAM, WOFOST, DAISY [44,45,46,47,48], the linear response functions in APSIM-NG and APSIM 7.9 that guide the developmental rate based on the cardinal temperatures assume a rapid increase in the developmental rate above the optimum temperature compared to the rate between base and optimum temperatures. Such an incorporation can lead to poor predictions and systematic errors [49,50]. In this study, the average temperature was below the optimum temperature during the studied cropping seasons, and the prediction of days to maturity in the evaluations was still not good, particularly for APSIM-NG. This indicated that the rate of thermal time accumulation with the default optimum temperature or the response function was not the primary reason in the models for errors in simulating phenology.
The average temperature at the studied location was below 4 °C (Tb in APSIM-NG) for many days during the crop maturation phase. Earlier studies have reported that when the cardinal temperature at the study locations falls outside of the range of temperature used for algorithm development, prediction discrepancies can be observed [24,51]. Thus, it is possible that the higher Tb in APSIM-NG and, most likely, the computation of thermal time without considering the crown temperature, could be the primary reason for having a worse performance than APSIM 7.9 in some aspects in the evaluations.
The two different phenology simulation mechanisms in both versions of APSIM essentially make them two completely different models. Models with relatively simpler phenological simulating mechanisms, accounting for the effects of temperature and daylength (e.g., WOFOST, and FASSET), to more complex mechanisms, accounting for the effects of temperature, photoperiod, and vernalisation (e.g., DAISY and HERMES), have been reported to simulate variable phenology [20,29]. Due to such responses, much work has been conducted recently to improve phenology submodules by altering the cardinal temperatures, particularly at higher temperatures [52,53]. The current study is in line with other findings [24,51], which suggest that caution is needed in the use of cardinal temperatures, particularly when a model is to be used in a new environment where the temperatures may fall outside the range for which the phenological model was developed. This consideration is critical at high latitudes, where the night temperature during seedling and maturity stages can be low, as observed in the current study, where APSIM-NG, with a higher base temperature, took more days to simulate the days to physiological maturity than APSIM 7.9. However, the broader implications of different cardinal temperatures in the models can be assessed when they are calibrated and evaluated with comprehensive datasets. Therefore, further studies are recommended in this direction.

3.3. Equifinality and Selection of the Parameters and Their Influence On Phenology

The combinations of the best parameters for the varieties varied under each calibration. The different combinations are not surprising, since the targeted set of parameters that was calibrated under Calibration 1_AN was different than under Calibration 2_PM and Calibration 3_PM. Besides the different parameters, the dataset for each calibration was different. For Calibration 1_AN, data were used only on days to anthesis from two environments. For Calibration 2_PM and Calibration 3_PM, the data were on days to physiological maturity for two and twenty environments, respectively.
When several parameters regulate a mechanism, different combinations of the parameter values can result in the same solution, termed Equifinality [37]. In this study, the phenomenon was also observed. For example, a thermal type parameter in APSIM 7.9, tt_end_of Juvenile, has a potential range between 300 and 450 for Alvari (Table S3) in the quest to find the best solution. The parameter tt_end_of Juvenile and others that control phenological development do not stand alone in terms of simulating days to anthesis. Their role is coupled with the other parameters, i.e., tt_floral_initiation, photop_sens, and vern_sens. Therefore, specific combinations of the different parameters can result in the same output, and hence the same RMSE. However, our hypothesis for the equifinality in this study is that several parameters that were used to calibrate days to anthesis using the smaller dataset resulted in the same solution. The data on days to anthesis for the varieties were available only for two seasons, which, in mathematical terms, has two solutions for one growth stage for each variety, and with more than 14,560 (in the case of APSIM 7.9) parameter combinations, it was very likely that several combinations would achieve those two solutions. If the data on anthesis and earlier growth stages were available for more years and locations, then the number of best parameter combinations might be lower.
The selection of the best set of parameters based on high frequency corroborates the fact that their relevance for regulating phenology is higher, and they hence have a greater chance to more accurately simulate days to anthesis under several environments (years and locations) than others. Although the parameters regulating the growth stages before anthesis were calibrated on arbitrary assumptions due to lack of data for those stages, their sequential dependency on simulating days to anthesis makes them useful for the cases when the simulation target is to achieve accurate days to anthesis.
For calibration of days to anthesis under Calibration 1_AN, the parameter PpSesnitivity in APSIM-NG was more stable, with only one value across the varieties. This indicated that the parameter was more stable in the study locations. However, PpSesnitivity’s synonym parameter photop_sens in APSIM 7.9 was relatively more variable. Contrastingly, VrnSensitivity in APSIM-NG was more variable compared to vern_sens in APSIM 7.9. The parameter MinimumLeafNumber, regulating the leaf number in APSIM-NG, had only one value for all varieties: six leaves. Since VrnSensitivity changes with Huan stage (or leaf number) and given the dependence of leaf number on BasePhylo, a stable leaf number for the varieties was the result of the variable BasePhylo in combination with variable VrnSensitivity.
For the calibration of days to physiological maturity of the varieties under Calibration 2_PM and Calibration 3_PM, parameter combinations and variability were different within each calibration and with Calibration 1_AN for both models (Table 2 and Table 4). For Calibration 2_PM with the smaller dataset, on days to physiological maturity, PpSensitivity was the same, and VrnSensitivity had the same variability but had different values for the varieties compared with Calibration 1_AN. BasePhylo and MinimumLeafNumber became more variable. However, for most of the varieties, the algorithm selected one of the lower values for these parameters, which was used for the calibration. For GrainFill, these were mid-range values. With the larger dataset on days to physiological maturity under Calibration 3_PM, the combinations of the parameter changed again compared to the other two calibrations. Notable variations were observed with VrnSensitivity, BasePhylo, and GrainFill. This suggested that the large variations in the parameter VrnSensitivity, BasePhylo, MinimumLeafNumber, and GrainFill were more influential in regulating the duration of growth stages than PpSensitivity and EarlyReproductivePpSensitivity, which were unchanged.
Similarly, the combinations and variations in the parameters of APSIM 7.9 with different datasets suggested that all calibrated parameters were influential in regulating the duration of the growth stages. Contrastingly, a larger variability of photop_sens suggested that it was more influential in APSIM 7.9 than its synonymous parameter PpSensitivity in APSIM-NG. It was evident in both models that combinations of and variations in the best parameters changed under the calibrations with different sized datasets.
The responses of the parameters in the evaluations suggested that the error terms (RMSE, RMSEsys, and RMSEnos) were lower when they were calibrated with larger datasets compared to their calibration with smaller datasets. However, by assessing both models’ performance based on the error terms and r2 across the evaluations, it was difficult to point out which model was better.

3.4. Further Development of APSIM Barley Models for Northern Regions

In most crop models, thermal time accumulation is adjusted by simple functions to incorporate the effects of vernalisation and photoperiod. Such adjustments are usually linear, which makes the models rigid (biased) to predict accurate phenology at the locations and conditions where the linearity does not hold [24,51]. Palosuo et al. [20] reported that models that simulate crop phenology as a function of temperature, photoperiod, and vernalisation (e.g., DAISY and DSSAT) simulated wheat phenology better than the models that simulated crop phenology as a function of temperature and photoperiod only (e.g., WOFOST and FASSET) for northwestern, central and southeastern Europe. While comparing the same models as in [20], Rötter et al. [31] observed that phenology prediction discrepancies were linked with the temperature and photoperiod factors in barley models for northern and central Europe. To ease such linear rigidity, APSIM 7.9 models were incorporated with broken linear functions to accommodate the responses to temperature and photoperiod [54], which could be the reason it performed better with a large dataset in the evaluation in this study. Vernalisation accumulation is dependent on the temperature ranges during the cropping season. The duration to satisfying the vernalisation requirements could be different if the base temperature is different in the model, which was the case in the two models, and thus a different response was observed in this study simulating days to anthesis using the same dataset. Spring barley varieties, which are typically assumed to have no vernalisation requirement, were used in this study. By calibrating the vernalisation parameters, we did not automatically assume that there was no vernalisation effect, like in some some studies [55], and unlike others, such as Rötter et al. [31]. This is because although breeding programs have selected for the no vernalisation requirement, it is not completely eliminated. The magnitude of vernalisation requirements would be different to winter barley varieties using winter barley APSIM models. In both versions of APSIM, a three-hourly interpolation of maximum and minimum temperature is computed from the daily maximum and minimum temperature based on the method described [56]. The mean of such interpolation could be the same, irrespective of the different temperature amplitudes [57]. The three-hour estimate has an advantage, particularly at high latitudes where the effect of maximum temperature is greater than the minimum temperature during the spring cropping season due to long days. To further improve the prediction accuracy of phenology at higher latitudes in APSIM, and particularly APSIM-NG, addressing the temperature amplitude can be a target.
As previously mentioned, the excessive rain and low temperature during the crop-maturing phase in 2017 caused barley varieties to lodge or not to mature, whereas, in the hot and dry 2018 cropping season, varieties developed rapidly and produced the poorest yield since the 1950s [58]. By capturing such uneven climates at high latitudes, using non-linear functions could further increase the prediction capacity of the phenology algorithms [57]. Unlike crown temperature usage in APSIM 7.9, APSIM-NG does not directly use this for thermal time computation; instead, higher cardinal temperatures are considered, which could be another research area to reconsider, besides base temperature, to increase the robustness of APSIM-NG at high latitudes.

4. Materials and Methods

4.1. Location and Agronomic Management

Crop growth and phenological data for twelve spring barley varieties were collected or acquired from on-going official variety trials for four locations in northern Sweden, and one in northern Finland, for five cropping seasons (2014–2018) (Table 1). The varieties are modern and bred for short growing seasons at high latitudes. The varieties include both two-row and six-row barley. The locations for the study in northern Sweden were: Röbäcksdalen (63.80 °N and 20.18 °E), Öjebyn (lat. 65.34 °N and 21.39 °E), Offer (63.11 °N and 17.74 °E) and Ås (63.24 °N and 14.56 °E), and, in Finland, the location was Ruukki (lat. 64.36 °N and 25.22 °E) (Figure 4). The sowing dates and agronomic practices for five seasons at the study locations are provided in Table 1.

4.2. Phenological Data Collection

Frequent phenological developmental monitoring was conducted at Röbäcksdalen to obtain the precise days to 50% anthesis, and Zadoks stage 65 [59], (hereafter, days to anthesis) of twelve barley varieties for two cropping seasons (2017 and 2018) (Table 1). Yellow to brownish stamens in the panicles were observed as indicators of anthesis. Ten random panicles were chosen and dissected for the stamens from the barley plots in two replications. For the days to maturity (Zadoks stage 90), plant samples were collected from a 0.5 m2 quadrat from the middle rows of the plots. At each sampling date, at least 1 m of space was left while taking the samples. Plant samples were threshed and clean grains were obtained for recording fresh weight, and subsequently dried at 65 °C until there was no weight change. Using the date of sampling and dry matter content in the grains, a linear or 2nd-order regression was used to estimate the date of 70% dry matter in the grains to mark the loss of grain color of the peduncle and, hence, end of grain filling [60], which was recorded as physiological maturity for the crop. This measurement is in accordance with APSIM’s identification of the stages, i.e., end of grain filling and physiological maturity for barley.
The data for the varieties were complete for both phenological stages for the 2017 and 2018 cropping seasons at Röbäcksdalen. For other seasons and locations, the data were only for physiological maturity. Since not all twelve varieties were sown in every season and location, as they were in 2017 and 2018 at Röbäcksdalen, the data on the available varieties that were used for model calibration and evaluation are presented in Table 1. Missing values indicate that the variety was not grown or the measurement was not performed on time to compute the maturity date.

4.3. Soil Characteristics

Soil moisture (volumetric water content, %) was measured at Röbäcksdalen using Diviner2000 (Sentek technologies, Australia) from 11 to 138 days after sowing (DAS) in 2017, and 21 to 90 DAS in 2018, with an interval of 3–15 days, depending on the canopy development (Figure 1). The measurements were less frequent during early growth stages and more frequent during and after the closing of the canopy. In 2017 and 2018, two and one access tubes, respectively, were installed in each barley plot. There was no chronic or long water stress period during any of the cropping seasons, neither in the measurements nor in simulated data. Field capacity, saturation, wilting point and the pH of the soil were measured at Röbacksdalen, Offer, Ås, and Ruukki and used in the APSIM simulations (Table S1). The Swedish soil data are taken from existing studies [61,62]. The soils at all study sites are classified as Podzols (FAO). The Röbäcksdalen soil has a clayey silt loam topsoil [63], the Offer soil has a silty clay loam topsoil, the Ås soil has gravelly loam topsoil, the Öjebyn soil has a sandy topsoil [64] and the Ruukki soil has a sandy mull topsoil [65]. Figure 1 shows the simulation of the models with the input data for only one location (Röbäcksdalen) and two years. The results provided a base to assume that both models can simulate the patterns of soil water dynamics at the study locations.

4.4. Climate

Daily weather data for Swedish locations were obtained from the online database of the Swedish Meteorological and Hydrological Institute (www.smhi.se, access date: 28 November 2019 for Röbäcksdalen; 29 November 2019 for Öjebyn; 22 October 2019 for Offer; 10 Spetember 2019 for Ås) and the Swedish University of Agricultural Sciences (www.ffe.slu.se/lm/LMHome, access date: 28 November 2019 for Röbäcksdalen; 29 November 2019 for Öjebyn; 22 October 2019 for Offer; 10 Spetember 2019 for Ås). For Finland, the data were obtained from the Finnish Meteorological Institute (https://en.ilmatieteenlaitos.fi/, access date: 20 November 2019). The preference was to obtain measured station data. However, for missing days, the data were obtained from the gridded system from the same sources. The weather factors affecting phenological development are presented in Figure 5.

4.5. Description of Phenology Modules of APSIM 7.9 and APSIM-NG

The term “APSIM Classic” is used to refer to APSIM7.x versions; in this study, it specifically refers to APSIM 7.9. APSIM Classic simulates barley phenology using a thermal time approach derived from [66,67], with the temperature response characterized by a base, optimum, and maximum, with linear relationships between the critical temperatures. The thermal time target between emergence and floral initiation is adjusted by day length [68,69]. The thermal time from emergence to floral initiation, divided by the plastochron, determines the total leaf number. Vernalisation affects phenology from emergence to floral initiation by regulating the duration of thermal time with a vernalization factor. The vernalization factor is computed from daily maximum and minimum temperature, which, based on the temperature requirements for vernalization, can prolong the vernalization sensitive phase until the requirements are met. Water and nitrogen stresses affect leaf appearance rate, which, depending on intensity, can delay phenology during the vegetative stages. The duration of the developmental stages from anthesis to end of grain filling is solely simulated through thermal time targets. The mechanism was implemented in most of the APSIM 7.x cereal models, including the spring barley model, by adapting the APSIM-Wheat module. A detailed description of the APSIM-Wheat module with the parameters and factors affecting phenology is available online (www.apsim.info, access date: 25 February 2020). Differences between APSIM-Wheat and APSIM-spring barley models are reported by [70].
As the plant modelling framework was imported from APSIM Classic [71] to develop the APSIM-NG models [72], many processes of the growth and development of the APSIM-Barley NG model are the same. The major difference arises in APSIM-NG’s use of Kirby’s Framework [42,43] to capture the effects of vernalisation and photoperiod. According to the framework, the simulation of the timing of the anthesis is dependent on the timing of the flag leaf appearance. It assumes that flag leaf appearance is dependent on the Final Leaf Number, which sets a target, and leaf appearance rate (phyllochron), which sets the rate to proceed toward the target to accumulate thermal time. The accumulation of thermal time and a cultivar specific phyllochron, which changes with the Haun stage, regulate the Leaf appearance rate [73]. The Haun system refers to the leaf development stage by expressing leaf length for each emerging leaf during the crop developmental cycle [74]. The requirements of low temperature for vernalisation and day length for photoperiod decrease with successive Haun stages as the temperature and day length increase and the crop proceeds towards anthesis. The APSIM-NG approach assumes more synchronicity and continuity among the processes, such as leaf emergence, floral primordia initiation, and anthesis, as they are simulated in parallel and are dependent on each other. However, APSIM 7.x simulates phenology with a target-based thermal time approach, independent of leaf appearance, and other processes related to anthesis, and hence there is less synchronicity and continuity among the processes. A detailed description of modelling the final leaf number, different developmental stages, and the effect of photoperiod and vernalisation in APSIM-NG, is available online (https://apsimnextgeneration.netlify.app/, access date: 25 February 2020).
Brown et al. [71] recently outlined the modernisation approaches for the advancements of interface and execution in APSIM-NG with APSIM “Classic”. Brown et al. [72] described the development of crop models in APSIM-next generation from the software, interface and the usage of larger datasets as a suggestion for the advancement of the modern crop models using the wheat module as a case study. The generic processes and functions in the APSIM-NG barley module are the same as the wheat module; however, they are not yet published [72].

4.6. Model Ccalibration: Using APSIM’s In-Built Factorial-Based Approach

The models were investigated in two steps: (1) calibration: variety-specific parameters regulating the phenology were calibrated using one dataset, and (2) evaluation: the model performance with the calibrated parameters was statistically evaluated using another dataset. For calibration, a factorial based approach applying an in-built system to create factorial combinations of pre-determined parameter value ranges in the Windows-based modelling software of APSIM-NG was performed. For APSIM 7.9, the same approach was implemented, using R programming with the apsimr package. Recently, this calibration approach was applied to calibrate the APSIM-NG maize model [32]. After executing the calibration set-ups, the best set of parameter values was extracted based on minimizing root mean square error (RMSE).
The calibration for each variety was performed with a large number of combinations of phenology-related parameter values (see Table 5). Descriptions of the parameters regulating specific growth stages are provided in Table S2. The selected phenology parameters and their ranges for creating the combinations were chosen based on a manual sensitivity analysis changing one parameter at a time, values of the default cultivars in APSIM, and expert knowledge of the crop phenology in the region. The values of parameters that were not calibrated were kept as in the base cultivars in the APSIM systems. The parameter tt_emergence, which regulates days to emergence after sowing in APSIM 7.9, was kept constant, and the value was determined based on the observations recorded at Röbäckdalen during 2017 and 2018. The total number of combinations of the parameters was 14,560 for APSIM 7.9 and 21,600 for APSIM-NG. The numbers for the two models were different because the calibration values of some parameters were higher in APSIM-NG. The simulations with all the parameter combinations were performed for all studied years and locations, as shown in Table 6. RMSEs of the estimated days to anthesis and physiological maturity were calculated using the R programming language. All the codes are provided in Supplementary Material (Codes in R programming).
Three independent calibrations were performed; calibration of days to anthesis (Calibration 1_AN) and two calibrations for days to physiological maturity (Calibration 2_PM and Calibration 3_PM). The dataset for these calibrations is presented in Table 6. Calibration 1_AN was for parameters that regulate days to anthesis, i.e., tt_end_of_juvenile, tt_floral_initiation, vern_sens, and photop_sens for APSIM 7.9 and BasePhyllo, MinimumLeafNumber, VrnSensitivity, PpSensitivity, and EarlyReproductivePpSensitivity for APSIM-NG (Table 5). Calibration 2_PM and Calibration 3_PM covered all the parameters that regulate days to anthesis plus the parameters that regulate days from anthesis to physiological maturity (tt_start_grain_fill for APSIM 7.9 and GrainFill for APSIM-NG). These parameters regulate thermal time accumulation and aim, directly or indirectly, to simulate different phenological stages. The models respond differently to calibrations using different datasets [75]; in this study, both versions of APSIM were tested with two different data sizes. Therefore, the two calibrations differed in terms of the data available for calibrations: Calibration 2_PM with fewer data and Calibration 3_PM with more data. A total of 20 combinations of locations and seasons generated 291,200 and 432,000 independent simulations for APSIM 7.9 and APSIM-NG in Calibration 3_PM.
The calibration and evaluation datasets were categorized based on (1) the available data on two growth stages, and (2) temperature variation at the locations. Data on anthesis and physiological maturity at Röbacksdalen for 2017–2018 were used in Calibration 1_AN and Calibration 2_PM, and the remaining dataset in Evaluation 2_PM was for the parameters which were in calibrated in Calibration 2_PM. Offer was the warmest location in terms of maximum temperature and Ruukki was the warmest in terms of minimum temperature. Because of this, the years 2015–2016 for Offer and Ruukki were selected in Evaluation 3_PM for the parameters calibrated in Calibration 3_PM using other locations and years of data. In Calibration 1_AN, to obtain a better calibration for anthesis, we used both years of available data, instead of splitting one year for calibration and another for evaluation. Although this approach did not allow us to evaluate the parameters calibrated under Calibration 1_AN, we chose this approach because the calibration of the parameters using two years of data is more robust compared with that using a single year of data.
The parameters which were calibrated under Calibration 1_AN could not be statistically evaluated, since all two seasons of data were used for the calibration of anthesis. Because of that, the statistical evaluation, using an independent dataset, was performed only for the calibrated parameters under Calibration 2_PM and Calibration 3_PM, which were denoted as Evaluation 2_PM and Evaluation 3_PM, respectively. For Evaluation 2_PM, data were used from 22 environments and for Evaluation 3_PM, data were from four different environments. The evaluation datasets and calibration datasets are summarized in Table 6.

4.7. Description of the Selection of Best Parameter Combinations

Both models were run with all the factorial combinations of the parameters for two years (2017 and 2018) for the Röbäcksdalen site. The days to anthesis were extracted for all the individual combinations for the two years. Further, the RMSE was computed using the observed data on 12 varieties for the corresponding site and years. In cases when several parameter combinations resulted in RMSE, the single set of parameter combinations was chosen for each variety based on the frequency of individual parameter values among the equally good parameter sets. For a detailed example of the procedure, please see Supplementary Materials.

4.8. Statistical Determinants to Assess the Model Calibration and Evaluation

The best sets of parameter combinations, i.e., those simulating a phenology closest to the observed data, were identified by minimizing root mean square error (RMSE) between estimated and observed dates of anthesis and physiological maturity. Model performances with the calibrated parameters were statistically evaluated with independent datasets. RMSE is a commonly used indicator to evaluate models but it does not highlight the source or type of error, which can be important in refining the models. Therefore, the systematic (RMSEsys) and non-systematic (RMSEnos) error components of RMSE were also computed [76]. RMSE indicates the mean deviation of the predicted values of the variables under evaluation with respect to the observed ones. RMSEsys indicates the systematic error associated with the model’s predictive ability. A higher magnitude of RMSEsys than zero indicates bias in the model performance. RMSEnos indicates unexplained error. Besides the RMSE error terms, the coefficient of determination (r2) was also computed for the overall relationship of simulation and measured datapoints.
As RMSE and r2 are the commonly used determinants, only computations of RMSEsys and RMSEnos are provided below
R M S E s y s = 1 n i = 1 n ( o i s ^ i ) 2
R M S E n o s = 1 n i = 1 n ( s i s ^ i ) 2
where o i is observed values of the tested variable, s i is simulated values of the tested variable, n is the number of entries for o i and s i , s ^ i = a + b o i , with a and b as the intercept and the slope of the regression model, respectively.

5. Conclusions

The better performance of the phenology algorithm in APSIM-NG in the calibration indicated its potential suitability at high latitudes compared to the algorithm in APSIM 7.9. However, in the evaluation APSIM-NG showed a greater tendency to overestimate days to physiological maturity. The combination of the best parameters simulating different growth stages close to the observed data changed with the calibrations and size of the datasets used for the calibration. The differences in both algorithms are potentially associated with the methods used to compute thermal time, using different cardinal temperatures. To further improve the performance and suitability of APSIM-NG at high latitudes, we suggest reviewing the cardinal temperatures in the current algorithm and reviewing the elimination of crown temperature in thermal time computation. Hence, more research is suggested on the conception of thermal time computations and their implementation in APSIM-NG.

Supplementary Materials

The following are available online at https://www.mdpi.com/2223-7747/10/3/443/s1. Codes in R programming, Table S1: Soil water characteristics at Röbäcksdalen, Öjebyn, Offer, Ås and Ruukki used to run APSIM 7.9 and APSIM-NG barley models. The bolded figures were assumed to complete the entries for the soil profile up to 1 m, Table S2: Description of the parameters used for the factorial combinations for calibration of APSIM 7.9 and APSIM-NG barley phenology models, Table S3: Combinations of the parameters resulted in the same lowest RMSE for Alvari and Anneli. The highlighted parameter combinations for the varieties were selected based on the frequency of individual parameters (Figure S2), which are presented in Table 5, Figure S1: Frequency of RMSE with the combinations of all calibrated parameters and varieties for Calibration 1_AN and RMSEs of Alvari and Anneli with the respective legend, RMSE_Alvari and RMSE_Anneli, Figure S2: Frequency of selection of individual parameter values in APSIM 7.9 resulted in the lowest RMSE for the varieties.

Author Contributions

Conceptualization, U.K., and D.P.; Methodology, U.K., D.P., J.M., T.P., and H.B.; Software, U.K., J.M., D.P., A.P., and H.B.; Validation, U.K., D.P., A.-M.G., T.P., A.P., M.A., H.B.; Formal Analysis, U.K., D.P., A.-M.G., and H.B.; Investigation, U.K., and D.P.; Resources, U.K., and D.P.; Data curation, U.K., A.-M.G., and T.P.; Writing—original draft preparation, U.K.; Writing—review and editing, U.K., D.P., J.M., A.P., T.P., G.B., and M.A.; Visualization, U.K., D.P., and A.-M.G.; Supervision, U.K., and D.P.; Project administration, U.K., and D.P.; Funding acquisition, U.K., and D.P. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by Swedish University of Agricultural Sciences and Regional jordbruksforskning för norra Sverige (RJN).

Acknowledgments

The study was made possible by the Swedish Infrastructure for Ecosystem Sciences (SITES) at Röbäcksdalen and official variety trials in northern Sweden and Finland.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Makowski, D.; Naud, C.; Jeuffroy, M.-H.; Barbottin, A.; Monod, H. Global sensitivity analysis for calculating the contribution of genetic parameters to the variance of crop model prediction. Reliab. Eng. Syst. Saf. 2006, 91, 1142–1147. [Google Scholar] [CrossRef]
  2. Confalonieri, R.; Bregaglio, S.; Acutis, M. Quantifying uncertainty in crop model predictions due to the uncertainty in the observations used for calibration. Ecol. Model. 2016, 328, 72–77. [Google Scholar] [CrossRef]
  3. 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. 2018, 94, 25–35. [Google Scholar] [CrossRef]
  4. Asseng, S.; Martre, P.; Maiorano, A.; Rötter, R.P.; O’Leary, G.J.; Fitzgerald, G.J.; Girousse, C.; Motzo, R.; Giunta, F.; Babar, M.A.; et al. Climate change impact and adaptation for wheat protein. Glob. Chang. Biol. 2019, 25, 155–173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Hunt, L.; Ash, A.; MacLeod, N.; McDonald, C.; Scanlan, J.; Bell, L.; Cowley, R.; Watson, I.; McIvor, J. Research Opportunities for Sustainable Productivity Improvement in the Northern Beef Industry: A Scoping Study; Meat & Livestock Australia: Sydney, Australia, 2019; pp. 157–210. [Google Scholar]
  6. Oliver, D.M.; Fish, R.D.; Winter, M.; Hodgson, C.J.; Heathwaite, A.L.; Chadwick, D.R. Valuing local knowledge as a source of expert data: Farmer engagement and the design of decision support systems. Environ. Model. Softw. 2012, 36, 76–85. [Google Scholar] [CrossRef]
  7. Parsons, D.; Nicholson, C.; Blake, R.W.; Ketterings, Q.; Ramírez-Aviles, L.; Fox, D.G.; Tedeschi, L.O.; Cherney, J.H. Development and evaluation of an integrated simulation model for assessing smallholder crop–livestock production in Yucatán, Mexico. Agric. Syst. 2011, 104, 1–12. [Google Scholar] [CrossRef] [Green Version]
  8. Phelan, D.C.; Harrison, M.T.; McLean, G.; Cox, H.; Pembleton, K.G.; Dean, G.J.; Parsons, D.; Richter, M.E.A.; Pengilley, G.; Hinton, S.J.; et al. Advancing a farmer decision support tool for agronomic decisions on rainfed and irrigated wheat cropping in Tasmania. Agric. Syst. 2018, 167, 113–124. [Google Scholar] [CrossRef]
  9. Masuya, Y.; Shimono, H. Mining a yield-trial database to identify high-yielding cultivars by simulation modeling: A case study for rice. J. Agric. Meteorol. 2017, 73, 51–58. [Google Scholar] [CrossRef] [Green Version]
  10. Kumar, U.; Laza, M.R.; Soulié, J.-C.; Pasco, R.; Mendez, K.V.S.; Dingkuhn, M. Compensatory phenotypic plasticity in irrigated rice: Sequential formation of yield components and simulation with SAMARA model. Field Crops Res. 2016, 193, 164–177. [Google Scholar] [CrossRef]
  11. Kumar, U.; Laza, M.R.; Soulié, J.C.; Pasco, R.; Mendez, K.V.S.; Dingkuhn, M. Analysis and simulation of phenotypic plasticity for traits contributing to yield potential in twelve rice genotypes. Field Crops Res. 2017, 202, 94–107. [Google Scholar] [CrossRef]
  12. Dingkuhn, M.; Pasco, R.; Pasuquin, J.M.; Damo, J.; Soulié, J.C.; Raboin, L.M.; Dusserre, J.; Sow, A.; Manneh, B.; Shrestha, S.; et al. Crop-model assisted phenomics and genome-wide association study for climate adaptation of indica rice. 1. Phenology. J. Exp. Bot. 2017, 68, 4369–4388. [Google Scholar] [CrossRef] [Green Version]
  13. Dingkuhn, M.; Pasco, R.; Pasuquin, J.M.; Damo, J.; Soulié, J.C.; Raboin, L.M.; Dusserre, J.; Sow, A.; Manneh, B.; Shrestha, S.; et al. Crop-model assisted phenomics and genome-wide association study for climate adaptation of indica rice. 2. Thermal stress and spikelet sterility. J. Exp. Bot. 2017, 68, 4389–4406. [Google Scholar] [CrossRef] [PubMed]
  14. Hammer, G.L.; Van Oosterom, E.; McLean, G.; Chapman, S.C.; Broad, I.; Harland, P.; Muchow, R.C. Adapting APSIM to model the physiology and genetics of complex adaptive traits in field crops. J. Exp. Bot. 2010, 61, 2185–2202. [Google Scholar] [CrossRef] [Green Version]
  15. Messina, C.D.; Podlich, D.; Dong, Z.; Samples, M.; Cooper, M. Yield-trait performance landscapes: From theory to application in breeding maize for drought tolerance. J. Exp. Bot. 2011, 62, 855–868. [Google Scholar] [CrossRef] [PubMed]
  16. Angulo, C.; Rötter, R.; Lock, R.; Enders, A.; Fronzek, S.; Ewert, F. Implication of crop model calibration strategies for assessing regional impacts of climate change in Europe. Agric. For. Meteorol. 2013, 170, 32–46. [Google Scholar] [CrossRef]
  17. Lobell, D.B.; Sibley, A.; Ivan Ortiz-Monasterio, J. Extreme heat effects on wheat senescence in India. Nat. Clim. Chang. 2012, 2, 186–189. [Google Scholar] [CrossRef]
  18. Rötter, R.P.; Carter, T.R.; Olesen, J.E.; Porter, J.R. Crop–climate models need an overhaul. Nat. Clim. Chang. 2011, 1, 175–177. [Google Scholar] [CrossRef]
  19. Martre, P.; Wallach, D.; Asseng, S.; Ewert, F.; Jones, J.W.; Rötter, R.P.; Boote, K.J.; Ruane, A.C.; Thorburn, P.J.; Cammarano, D.; et al. Multimodel ensembles of wheat growth: Many models are better than one. Glob. Chang. Biol. 2015, 21, 911–925. [Google Scholar] [CrossRef]
  20. Palosuo, T.; Kersebaum, K.C.; Angulo, C.; Hlavinka, P.; Moriondo, M.; Olesen, J.E.; Patil, R.H.; Ruget, F.; Rumbaur, C.; Takáč, J.; et al. Simulation of winter wheat yield and its variability in different climates of Europe: A comparison of eight crop growth models. Eur. J. Agron. 2011, 35, 103–114. [Google Scholar] [CrossRef] [Green Version]
  21. Tao, F.; Rötter, R.P.; Palosuo, T.; Gregorio Hernández Díaz-Ambrona, C.; Mínguez, M.I.; Semenov, M.A.; Kersebaum, K.C.; Nendel, C.; Specka, X.; Hoffmann, H.; et al. Contribution of crop model structure, parameters and climate projections to uncertainty in climate change impact assessments. Glob. Chang. Biol. 2018, 24, 1291–1307. [Google Scholar] [CrossRef]
  22. Rosenzweig, C.; Jones, J.W.; Hatfield, J.L.; Ruane, A.C.; Boote, K.J.; Thorburn, P.; Antle, J.M.; Nelson, G.C.; Porter, C.; Janssen, S.; et al. The Agricultural Model Intercomparison and Improvement Project (AgMIP): Protocols and pilot studies. Agric. For. Meteorol. 2013, 170, 166–182. [Google Scholar] [CrossRef] [Green Version]
  23. Ewert, F.; Rötter, R.P.; Bindi, M.; Webber, H.; Trnka, M.; Kersebaum, K.C.; Olesen, J.E.; van Ittersum, M.K.; Janssen, S.; Rivington, M.; et al. Crop modelling for integrated assessment of risk to food production from climate change. Environ. Model. Softw. 2015, 72, 287–303. [Google Scholar] [CrossRef]
  24. Van Oort, P.A.J.; Zhang, T.; de Vries, M.E.; Heinemann, A.B.; Meinke, H. Correlation between temperature and phenology prediction error in rice (Oryza sativa L.). Agric. For. Meteorol. 2011, 151, 1545–1555. [Google Scholar] [CrossRef]
  25. White, M.A.; Thornton, P.E.; Running, S.W. A continental phenology model for monitoring vegetation responses to interannual climatic variability. Glob. Biogeochem. Cycles 1997, 11, 217–234. [Google Scholar] [CrossRef]
  26. Peake, A.S.; Das, B.T.; Bell, K.L.; Gardner, M.; Poole, N. Effect of variable crop duration on grain yield of irrigated spring-wheat when flowering is synchronised. Field Crops Res. 2018, 228, 183–194. [Google Scholar] [CrossRef]
  27. Stapper, M.; Fischer, R.A. Genotype, sowing date and plant spacing influence on high-yielding irrigated wheat in southern New South Wales. III. Potential yields and optimum flowering dates. Aust. J. Agric. Res. 1990, 41, 1043. [Google Scholar] [CrossRef]
  28. Baldocchi, D.; Falge, E.; Wilson, K. A spectral analysis of biosphere–atmosphere trace gas flux densities and meteorological variables across hour to multi-year time scales. Agric. For. Meteorol. 2001, 107, 1–27. [Google Scholar] [CrossRef]
  29. Zhang, S.; Tao, F. Modeling the response of rice phenology to climate change and variability in different climatic zones: Comparisons of five models. Eur. J. Agron. 2013, 45, 165–176. [Google Scholar] [CrossRef]
  30. Cao, W.; Moss, D.N. Modelling phasic development in wheat: A conceptual integration of physiological components. J. Agric. Sci. 1997, 129, 163–172. [Google Scholar] [CrossRef]
  31. Rötter, R.P.; Palosuo, T.; Kersebaum, K.C.; Angulo, C.; Bindi, M.; Ewert, F.; Ferrise, R.; Hlavinka, P.; Moriondo, M.; Nendel, C.; et al. Simulation of spring barley yield in different climatic zones of Northern and Central Europe: A comparison of nine crop models. Field Crops Res. 2012, 133, 23–36. [Google Scholar] [CrossRef]
  32. Morel, J.; Parsons, D.; Halling, M.A.; Kumar, U.; Peake, A.; Bergkvist, G.; Brown, H.; Hetta, M. Challenges for Simulating Growth and Phenology of Silage Maize in a Nordic Climate with APSIM. Agronomy 2020, 10, 645. [Google Scholar] [CrossRef]
  33. Holzworth, D.; Huth, N.I.; Fainges, J.; Brown, H.; Zurcher, E.; Cichota, R.; Verrall, S.; Herrmann, N.I.; Zheng, B.; Snow, V. APSIM Next Generation: Overcoming challenges in modernising a farming systems model. Environ. Model. Soft. 2018, 103, 43–51. [Google Scholar] [CrossRef]
  34. Saldivar, S.O.S. Cereals: Types and Composition. In Encyclopedia of Food and Health, 1st ed.; Benjamin, C., Finglas, P.M., Toldra, F., Eds.; Acedemic Press: Cambridge, MA, USA; Elsevier Ltd.: Oxford, UK, 2015; Volume 1, p. 703. [Google Scholar]
  35. Olesen, J.E.; Trnka, M.; Kersebaum, K.C.; Skjelvåg, A.O.; Seguin, B.; Peltonen-Sainio, P.; Rossi, F.; Kozyra, J.; Micale, F. Impacts and adaptation of European crop production systems to climate change. Eur. J. Agron. 2011, 34, 96–112. [Google Scholar] [CrossRef]
  36. Ruosteenoja, K.; Räisänen, J.; Venäläinen, A.; Kämäräinen, M. Projections for the duration and degree days of the thermal growing season in Europe derived from CMIP5 model output. Int. J. Climatol. 2015, 36, 3039–3055. [Google Scholar] [CrossRef]
  37. Beven, K.; Freer, J. Equifinality, data assimilation, and uncertainty estimation in mechanistic modelling of complex environmental systems using the GLUE methodology. J. Hydrol. 2001, 249, 11–29. [Google Scholar] [CrossRef]
  38. Jagtap, S.S.; Jones, J.W. Adaptation and evaluation of the CROPGRO-soybean model to predict regional yield and production. Agric. Ecosyst. Environ. 2002, 93, 73–85. [Google Scholar] [CrossRef]
  39. Challinor, A.J.; Ewert, F.; Arnold, S.; Simelton, E.; Fraser, E. Crops and climate change: Progress, trends, and challenges in simulating impacts and informing adaptation. J. Exp. Bot. 2009, 60, 2775–2789. [Google Scholar] [CrossRef]
  40. Ahmed, M.; Ahmad, S.; Raza, M.A.; Kumar, U.; Ansar, M.; Shah, G.A.; Parsons, D.; Hoogenboom, G.; Palosuo, T.; Seidel, S. Models Calibration and Evaluation. In Systems Modeling; Ahmed, M., Ed.; Springer: Singapore, 2020; pp. 149–176. [Google Scholar]
  41. Wolf, J. Comparison of two potato simulation models under climate change. I. Model calibration and sensitivity analyses. Clim. Res. 1996, 7, 253–270. [Google Scholar] [CrossRef] [Green Version]
  42. Kirby, E.J.M. Analysis of leaf, stem and ear growth in wheat from terminal spikelet stage to anthesis. Field Crops Res. 1988, 18, 127–140. [Google Scholar] [CrossRef]
  43. Kirby, E.J.M. Co-ordination of leaf emergence and leaf and spikelet primordium initiation in wheat. Field Crops Res. 1990, 25, 253–264. [Google Scholar] [CrossRef]
  44. Boote, K.J.; Prasad, V.; Allen, L.H.; Singh, P.; Jones, J.W. Modeling sensitivity of grain yield to elevated temperature in the DSSAT crop models for peanut, soybean, dry bean, chickpea, sorghum, and millet. Eur. J. Agron. 2018, 100, 99–109. [Google Scholar] [CrossRef]
  45. Brisson, N.; Ruget, F.; Gate, P.; Lorgeou, J.; Nicoullaud, B.; Tayot, X.; Plenet, D.; Jeuffroy, M.-H.; Bouthier, A.; Ripoche, D.; et al. STICS: A generic model for simulating crops and their water and nitrogen balances. II. Model validation for wheat and maize. Agronomie 2002, 22, 69–92. [Google Scholar] [CrossRef]
  46. Challinor, A.J.; Wheeler, T.R.; Craufurd, P.Q.; Slingo, J.M.; Grimes, D.I.F. Design and optimisation of a large-area process-based model for annual crops. Agric. For. Meteorol. 2004, 124, 99–120. [Google Scholar] [CrossRef]
  47. Supit, I.; Hooijer, A.A.; van Diepen, C.A. System Description of the WOFOST 6.0 Crop Simulation Model Implemented in CGMS; CGMS Publication EUR 15956 EN; Office for Official Publications of the E.U.: Luxembourg, 1994; p. 150. [Google Scholar]
  48. Abrahamsen, P.; Hansen, S. Daisy: An open soil–crop–atmosphere system model. Environ. Model. Softw. 2000, 15, 313–330. [Google Scholar] [CrossRef]
  49. Sadras, V.O.; Monzon, J.P. Modelled wheat phenology captures rising temperature trends: Shortened time to flowering and maturity in Australia and Argentina. Field Crops Res. 2006, 99, 136–146. [Google Scholar] [CrossRef]
  50. Challinor, A.J.; Wheeler, T.R. Crop yield reduction in the tropics under climate change: Processes and uncertainties. Agric. For. Meteorol. 2008, 148, 343–356. [Google Scholar] [CrossRef] [Green Version]
  51. Zhang, T.; Zhu, J.; Yang, X. Non-stationary thermal time accumulation reduces the predictability of climate change effects on agriculture. Agric. For. Meteorol. 2008, 148, 1412–1418. [Google Scholar] [CrossRef]
  52. McMaster, G.S.; White, J.W.; Hunt, L.A.; Jamieson, P.D.; Dhillon, S.S.; Ortiz-Monasterio, J.I. Simulating the influence of vernalization, photoperiod and optimum temperature on wheat developmental rates. Ann. Bot. 2008, 102, 561–569. [Google Scholar] [CrossRef]
  53. White, J.W.; Kimball, B.A.; Wall, G.W.; Ottman, M.J.; Hunt, L.A. Responses of time of anthesis and maturity to sowing dates and infrared warming in spring wheat. Field Crops Res. 2011, 124, 213–222. [Google Scholar] [CrossRef]
  54. Hammer, G.L.; Carberry, P.S.; Muchow, R.C. Modelling genotypic and environmental control of leaf area dynamics in grain sorghum. I. Whole plant level. Field Crops Res. 1993, 33, 293–310. [Google Scholar] [CrossRef]
  55. Tao, F.; Rötter, R.P.; Palosuo, T.; Díaz-Ambrona, C.G.H.; Mínguez, M.I.; Semenov, M.A.; Kersebaum, K.C.; Nendel, C.; Cammarano, D.; Hoffmann, H.; et al. Designing future barley ideotypes using a crop model ensemble. Eur. J. Agron. 2017, 82, 144–162. [Google Scholar] [CrossRef]
  56. Jones, C.A.; Kiniry, J.R. CERES-Maize: A simulation model of Maize Growth and Development; Texas A&M University Press: College Station, TX, USA, 1986; p. 194. [Google Scholar]
  57. Porter, J.R.; Semenov, M.A. Crop responses to climatic variation. Philos. Trans. R. Soc. B Biol. Sci. 2005, 360, 2021–2035. [Google Scholar] [CrossRef] [PubMed]
  58. Swedish Board of Agriculture. Production of Cereal Crops, Dried Pulses and Oilseed Crops in 2018 Preliminary Results for the Whole Country; Report Jo 19 SM 1801; Swedish Board of Agriculture: Örebro, Sweden, 2018.
  59. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stages of cereals. Weed Res. 1974, 14, 415–421. [Google Scholar] [CrossRef]
  60. Juskiw, P.E.; Jame, Y.-W.; Kryzanowski, L. Phenological development of spring barley in a short-season growing area. Agron. J. 2001, 93, 370–379. [Google Scholar] [CrossRef] [Green Version]
  61. Morel, J.; Kumar, U.; Parsons, D. Laboratory Measured Characteristics of 6 Swedish Soil Profiles; Mendeley Data, V1; Mendeley Data: New York, NY, USA, 2019. [Google Scholar] [CrossRef]
  62. Andersso, S.; Wiklert, P. Studier av markprofiler i svenska åkerjordar. Del II.; Norrbottens, Västerbottens-, Västernorrlands- och Jämtlands län; Report 104; Institutionen För Markvetenskap, Sveriges Lantbruksuniversitet: Uppsala, Sweden, 1977; p. 72. [Google Scholar]
  63. Ericson, L. Soil physical properties, organic carbon and trends in barley yield in four different crop rotations. In Proceedings of the 1st Circumpolar Agricultural Conference, Whitehorse, YT, Canada, 26 September–2 October 1992; Scott Smith, C.A., Ed.; pp. 189–193. [Google Scholar]
  64. Palmborg, C. Kol och Kväve i Mark och Grödor i försök med Monokulturer och Växtföljder; En delstudie inom projektet Climate CAFE; Institutionen För Markvetenskap, Sveriges Lantbruksuniversitet: Umeå, Sweden, 2019; p. 12. [Google Scholar]
  65. Saarela, I.; Järvi, A.; Hakkola, H.; Rinne, K. Phosphorus status of diverse soils in Finland as influenced by long-term P fertilisation I. Native and previously applied P at 24 experimental sites. Agric. Food Sci. Finl. 2003, 12, 117–132. [Google Scholar] [CrossRef]
  66. Muchow, R.C.; Carberry, P.S. Phenology and leaf-area development in a tropical grain sorghum. Field Crops Res. 1990, 23, 221–237. [Google Scholar] [CrossRef]
  67. Hammer, G.L.; Muchow, R.C. Assessing climatic risk to sorghum production in water-limited subtropical environments II. Effects of planting date, soil water at planting, and cultivar phenology. Field Crops Res. 1994, 36, 235–246. [Google Scholar] [CrossRef]
  68. Hammer, G.L.; Vanderlip, R.L.; Gibson, G.; Wade, L.J.; Henzell, R.G.; Younger, D.R.; And, J.W.; Dale, A.B. Genotype-by-Environment Interaction in Grain Sorghum. II. Effects of Temperature and Photoperiod on Ontogeny. Crop Sci. 1989, 29, 376–384. [Google Scholar] [CrossRef]
  69. Ravi Kumar, S.; Hammer, G.L.; Broad, I.; Harland, P.; McLean, G. Modelling environmental effects on phenology and canopy development of diverse sorghum genotypes. Field Crops Res. 2009, 111, 157–165. [Google Scholar] [CrossRef]
  70. Manschadi, A.M.; Hochman, Z.; McLean, G.; DeVoil, P.; Holzworth, D.; Meinke, H. APSIM-Barley model e adaptation of a wheat model to simulate barley growth and development. In Proceedings of the 13th Australian Agronomy Conference, Perth, Australia, 10–14 September 2006; Turner, N.C., Acuna, T., Eds.; 2006. [Google Scholar]
  71. Brown, H.E.; Huth, N.I.; Holzworth, D. Crop model improvement in APSIM: Using wheat as a case study. Eur. J. Agron. 2018, 100, 141–150. [Google Scholar] [CrossRef]
  72. Brown, H.E.; Huth, N.I.; Holzworth, D.; Teixeira, E.I.; Zyskowski, R.F.; Hargreaves, J.N.G.; Moot, D.J. Plant Modelling Framework: Software for building and running crop models on the APSIM platform. Environ. Model. Soft. 2014, 62, 385–398. [Google Scholar] [CrossRef] [Green Version]
  73. Jamieson, P.D.; Brooking, I.R.; Porter, J.R.; Wilson, D.R. Prediction of leaf appearance in wheat: A question of temperature. Field Crops Res. 1995, 41, 35–44. [Google Scholar] [CrossRef]
  74. Haun, J.R. Visual Quantification of Wheat Development. Agron. J. 1973, 65, 116–119. [Google Scholar] [CrossRef]
  75. He, D.; Wang, E.; Wang, J.; Robertson, M.J. Data requirement for effective calibration of process-based crop models. Agric. For. Meteorol. 2017, 234–235, 136–148. [Google Scholar] [CrossRef]
  76. Willmott, C.J. On the validation of models. Phys. Geogr. 1981, 2, 184–194. [Google Scholar] [CrossRef]
Figure 1. Soil moisture dynamics of 1-m soil profile with ten layers of equal depths during the 2017–2018 cropping seasons. The upper part of the figure represents the whole profile (0 to 100 cm) and the lower represents the topsoil layers (0–30 cm). Volumetric water content (%) (VWC) is the volume of water per cubic unit of soil, shown as available water (mm). The points represent measured data and continuous lines represent simulated data. Thick lines are for APSIM 7.9, and thin lines are for APSIM-NG. Wilting point is defined as 6.39% of VWC based on the data in Table S1. 2017: Red lines and points, 2018: green lines and points.
Figure 1. Soil moisture dynamics of 1-m soil profile with ten layers of equal depths during the 2017–2018 cropping seasons. The upper part of the figure represents the whole profile (0 to 100 cm) and the lower represents the topsoil layers (0–30 cm). Volumetric water content (%) (VWC) is the volume of water per cubic unit of soil, shown as available water (mm). The points represent measured data and continuous lines represent simulated data. Thick lines are for APSIM 7.9, and thin lines are for APSIM-NG. Wilting point is defined as 6.39% of VWC based on the data in Table S1. 2017: Red lines and points, 2018: green lines and points.
Plants 10 00443 g001
Figure 2. Evaluation of APSIM 7.9 barley model. (A) Evaluation 2_PM represents the evaluation with the data from 22 environments; (B) Evaluation 3_PM represents evaluation with the data from 4 environments. Each point represents simulation and observed data for one variety. Observed DAS and Simulation DAS are observed and simulated days after sowing outputs, respectively.
Figure 2. Evaluation of APSIM 7.9 barley model. (A) Evaluation 2_PM represents the evaluation with the data from 22 environments; (B) Evaluation 3_PM represents evaluation with the data from 4 environments. Each point represents simulation and observed data for one variety. Observed DAS and Simulation DAS are observed and simulated days after sowing outputs, respectively.
Plants 10 00443 g002
Figure 3. Evaluation of APSIM-NG barley model. (A) Evaluation 2_PM represents the evaluation with the data from 22 environments; (B) Evaluation 3_PM represents evaluation with the data from 4 environments. Each point represents simulation and observed data for one variety. Simulaitions that did not reach maturity were not included in the regression line.
Figure 3. Evaluation of APSIM-NG barley model. (A) Evaluation 2_PM represents the evaluation with the data from 22 environments; (B) Evaluation 3_PM represents evaluation with the data from 4 environments. Each point represents simulation and observed data for one variety. Simulaitions that did not reach maturity were not included in the regression line.
Plants 10 00443 g003
Figure 4. Study locations in northern Sweden and Finland.
Figure 4. Study locations in northern Sweden and Finland.
Plants 10 00443 g004
Figure 5. Weather data from 2014 to 2018 for the studied locations in northern Sweden and Finland. (A) Black points represent global radiation and the red line represents the radiation for clear sky. (B) The black circles represent maximum temperature and olive green circles represent minimum temperature. (C) The bars represent the daily average precipitation for each month, and error bars are standard errors. Area graphs are cumulated precipitation from January to December with the magnitudes on the right-hand y-axis.
Figure 5. Weather data from 2014 to 2018 for the studied locations in northern Sweden and Finland. (A) Black points represent global radiation and the red line represents the radiation for clear sky. (B) The black circles represent maximum temperature and olive green circles represent minimum temperature. (C) The bars represent the daily average precipitation for each month, and error bars are standard errors. Area graphs are cumulated precipitation from January to December with the magnitudes on the right-hand y-axis.
Plants 10 00443 g005aPlants 10 00443 g005bPlants 10 00443 g005c
Table 1. Phenology data for twelve spring barley varieties for 2014–2018 cropping seasons with the sowing dates at the studied locations. AN: days to anthesis after sowing; PM: days to physiological maturity after sowing; 6R and 2R refer to 6- and 2-row barley varieties.
Table 1. Phenology data for twelve spring barley varieties for 2014–2018 cropping seasons with the sowing dates at the studied locations. AN: days to anthesis after sowing; PM: days to physiological maturity after sowing; 6R and 2R refer to 6- and 2-row barley varieties.
VarietiesRöbäcksdalenÖjebynOfferÅsRuukki
201420152016201720182014201520162017201820142015201620172018201420152017201820142015201620172018
PMPMPMANPMANPMPMPMPMPMPMPMPMPMPMPMPMPMPMPMPMPMPMPMPM
Alvari (6R) 591165384 11085 9975 131 84
Anneli (2R) 100591165487 10410887 10710377 13683
Aukusti (6R)8410587541095184741009910581711008998748192125 781008710478
GN10063 (6R) 98521225184 10010879 10374 141
Judit (6R)76108935410851837398971057870988993748083126
Kaarle (6R) 117100591305487 112126112 109116109 105146 80
Kannas (2R)88111105571245484841101071088679108 10475979814283
Rödhette (6R) 591325489 11692 10977 14482
Severi (6R)811169159118518776106102109807210695102748897130
Vertti (6R) 9891521175181 10391107 1018494 88124 78
Vilde (6R)83103905712054857510497108827210388101758696126 8010189104
Vilgott (2R)96121104591185588871141121129081112991107710310913989
Mean84110965711953857810610410984741059610275899613484801018810479
SD7863722661035551161998331101
Sowing date28
May
9
Jun
2
Jun
3
Jun
23
May
23
May
18
Jun
7
Jun
10
Jun
12
Jun
28
May
27
May
30
May
1
Jun
19
May
30
May
28
May
25
May
20
May
20
May
22
May
9
May
30
May
16
May
Agronomic practices: Fertiliser application, recommended; seed rate, 500 seeds m−2; Plot size, 1.5 × 9 m; Number of rows/plots, 12; row spacing, 12.5 cm.
Table 2. Calibrated parameters for APSIM 7.9. Calibration 1_ AN, Calibration 2_ PM, and Calibration 3_PM represent calibration for days to anthesis, days to physiological maturity, and days to physiological maturity using data from 2, 2, and 22 environments, respectively. oCd and RMSE are degree days and root mean square error, respectively.
Table 2. Calibrated parameters for APSIM 7.9. Calibration 1_ AN, Calibration 2_ PM, and Calibration 3_PM represent calibration for days to anthesis, days to physiological maturity, and days to physiological maturity using data from 2, 2, and 22 environments, respectively. oCd and RMSE are degree days and root mean square error, respectively.
Variety
ParameterAlvariAnneli AukustiGN10063JuditKaarleKannasRodhetteSeveriVerttiVildeVilgot
Calibration 1_ AN
photop_sens110001011000
vern_sens0010.5100000.500.5
tt_end_of_juvenile (oCd)300300200200200300300300300200300300
tt_floral_initiation (oCd)320320300320300320320320300320320300
tt_start_grain_fill (oCd)
RMSE (d)1.30.7010.40.700.720.800
Calibration 2_PM
photop_sens610313160600
vern_sens000.5000000.500.50.5
tt_end_of_juvenile (oCd)250300200200200350300300250250300250
tt_floral_initiation (oCd)300320300320300300320320300300320300
tt_start_grain_fill (oCd)525525575625575550575550575525500575
RMSE (d)4.532.27.52.58.68.28.546.763.5
Calibration 3_PM
photop_sens113313130113
vern_sens000000000000
tt_end_of_juvenile (oCd)450400250450400450400450450450450450
tt_floral_initiation (oCd)260320300240280240320240240240240260
tt_start_grain_fill (oCd)400400450400350450400450400350350400
RMSE (d)8.37.87.510.48.71110.1119.98.78.57.6
Table 3. Model evaluation determinants. Evaluation 2_PM represents statistical evaluation of the parameter calibrated under Calibration 2_ PM; Evaluation 3_PM represents statistical evaluation of the parameter calibrated under Calibration 3_ PM. Intercept, Slope, RMSE, RMSEsys, RMSEnos, and r2 represent regression intercept and slope, root mean square error, systematic component of RMSE, non-systematic component of RMSE, and coefficient of determination, respectively.
Table 3. Model evaluation determinants. Evaluation 2_PM represents statistical evaluation of the parameter calibrated under Calibration 2_ PM; Evaluation 3_PM represents statistical evaluation of the parameter calibrated under Calibration 3_ PM. Intercept, Slope, RMSE, RMSEsys, RMSEnos, and r2 represent regression intercept and slope, root mean square error, systematic component of RMSE, non-systematic component of RMSE, and coefficient of determination, respectively.
InterceptAPSIM 7.9APSIM-NG
Evaluation 2_PM *Evaluation 3_PMEvaluation 2_PM ¤Evaluation 3_PM
−32.7−2.22.318
Slope1.20.990.870.81
RMSE11.44.7313.65.8
RMSEsys10.173.110.192
RMSEnos20.647.1522.635.81
r20.740.810.780.63
* Based on complete 159 pairs of simulated and observed datapoints (Figure 2A). ¤ Based on complete 139 pairs of simulated and observed datapoints. 20 pairs that did not reach physiological maturity were excluded (Figure 3A).
Table 4. Calibrated parameters for APSIM Next Generation (APSIM-NG). Calibration 1_ AN, Calibration 2_ PM, and Calibration 3_ PM represent calibration for days to anthesis, days to physiological maturity, and days to physiological maturity using data from 2, 2, and 22 environments, respectively.
Table 4. Calibrated parameters for APSIM Next Generation (APSIM-NG). Calibration 1_ AN, Calibration 2_ PM, and Calibration 3_ PM represent calibration for days to anthesis, days to physiological maturity, and days to physiological maturity using data from 2, 2, and 22 environments, respectively.
Variety
ParameterAlvariAnneli AukustiGN10063JuditKaarleKannasRodhetteSeveriVerttiVildeVilgot
Calibration 1_AN
PpSensitivity000000000000
VrnSensitivity1.51.50.50.511.51.51.50.50.51.51.5
BasePhyllo (oCd)606060605560606060556060
MinimumLeafNumber666666666666
GrainFill (oCd)
EarlyReproductivePpSensitivity000000000000
RMSE (d)000000000000
Calibration 2_PM
PpSensitivity000000000000
VrnSensitivity1.51.51.51110.50.51.51.50.51.5
BasePhyllo (oCd)505550506050505050505060
MinimumLeafNumber556.5555556655
GrainFill (oCd)525525450550450600600625500500575500
EarlyReproductivePpSensitivity000000000000
RMSE (d)000000000000
Calibration 3_PM
PpSensitivity000000000000
VrnSensitivity0.50000.530.50.530.50.50
BasePhyllo (oCd)656050656555605055606050
MinimumLeafNumber566.56.555555556.5
GrainFill (oCd)450450450400400450500575400450450525
EarlyReproductivePpSensitivity000000000000
RMSE (d)4.64.15.56.27.06.67.07.57.26.86.74.9
Table 5. Parameters and their levels applied in factorial calibration of APSIM 7.9 and APSIM-NG barley phenology models. The total number of combinations for APSIM 7.9 and APSIM-NG were 14,560 and 21,600, respectively.
Table 5. Parameters and their levels applied in factorial calibration of APSIM 7.9 and APSIM-NG barley phenology models. The total number of combinations for APSIM 7.9 and APSIM-NG were 14,560 and 21,600, respectively.
Factorial-Based Calibration
APSIM 7.9APSIM-NG
ParameterDefault ValueLevelParameterDefault ValueLevel
tt_emergence (°Cd)140BasePhyllo (oCd) *5050, 55, 60, 65, 70, 75
tt_end_of_juvenile (°Cd) *400150, 200, 250, 300, 350, 400, 450MinimumLeafNumber *95, 5.5, 6, 6.5, 7
tt_floral_initiation (°Cd) *230180, 200, 220, 240, 260, 280, 300, 320 GrainFill (oCd)¤540350, 400, 450, 500, 525, 550, 575, 600, 625, 650, 675, 700
tt_start_grain_fill (°Cd)¤545350, 400, 450, 500, 525, 550, 575, 600, 625, 650, 675, 700, 750VrnSesnsitivity*00, 0.5, 1, 1.5, 3
vern_sens *1.50, 0.5, 1, 1.5, 3PpSesnsitivity *30, 1, 3, 6
photop_sens *30, 1, 3, 6EarlyReproductivePpSesnsitivity *00, 1, 2
* Regulating days to anthesis. ¤ Regulating days after anthesis to physiological maturity.
Table 6. Calibration and evaluation datasets for APSIM 7.9 and APSIM-NG spring barley models. Calibration 1_AN represents calibration of days to anthesis and Calibration 2_PM, and Calibration 3_PM represent calibration of days to physiological maturity. Evaluation 2_PM and Evaluation 3_PM represent the independent evaluation of Calibration 2_PM, and Calibration 3_PM, respectively.
Table 6. Calibration and evaluation datasets for APSIM 7.9 and APSIM-NG spring barley models. Calibration 1_AN represents calibration of days to anthesis and Calibration 2_PM, and Calibration 3_PM represent calibration of days to physiological maturity. Evaluation 2_PM and Evaluation 3_PM represent the independent evaluation of Calibration 2_PM, and Calibration 3_PM, respectively.
Growth StageLocationCalibration DatasetEvaluation Dataset
Calibration 1_ANCalibration 2_PMCalibration 3_PM Evaluation 2_PMEvaluation 3_PM
SeasonSeasonSeasonLocationSeasonSeason
Days to anthesisRöbäcksdalen2017–2018
Days to physiological maturityRöbäcksdalen 2017–2018 Röbäcksdalen2014–2016
Days to physiological maturityRöbäcksdalen 2014–2018
Ås 2014, 2016–2018Ås2014, 2016–2018
Öjebyn 2014–2018Öjebyn2014–2018
Offer 2014, 2017–2018Offer2014–20182015–2016
Ruukki 2014, 2017–2018Ruukki2014–20182015–2016
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kumar, U.; Morel, J.; Bergkvist, G.; Palosuo, T.; Gustavsson, A.-M.; Peake, A.; Brown, H.; Ahmed, M.; Parsons, D. Comparative Analysis of Phenology Algorithms of the Spring Barley Model in APSIM 7.9 and APSIM Next Generation: A Case Study for High Latitudes. Plants 2021, 10, 443. https://doi.org/10.3390/plants10030443

AMA Style

Kumar U, Morel J, Bergkvist G, Palosuo T, Gustavsson A-M, Peake A, Brown H, Ahmed M, Parsons D. Comparative Analysis of Phenology Algorithms of the Spring Barley Model in APSIM 7.9 and APSIM Next Generation: A Case Study for High Latitudes. Plants. 2021; 10(3):443. https://doi.org/10.3390/plants10030443

Chicago/Turabian Style

Kumar, Uttam, Julien Morel, Göran Bergkvist, Taru Palosuo, Anne-Maj Gustavsson, Allan Peake, Hamish Brown, Mukhtar Ahmed, and David Parsons. 2021. "Comparative Analysis of Phenology Algorithms of the Spring Barley Model in APSIM 7.9 and APSIM Next Generation: A Case Study for High Latitudes" Plants 10, no. 3: 443. https://doi.org/10.3390/plants10030443

APA Style

Kumar, U., Morel, J., Bergkvist, G., Palosuo, T., Gustavsson, A. -M., Peake, A., Brown, H., Ahmed, M., & Parsons, D. (2021). Comparative Analysis of Phenology Algorithms of the Spring Barley Model in APSIM 7.9 and APSIM Next Generation: A Case Study for High Latitudes. Plants, 10(3), 443. https://doi.org/10.3390/plants10030443

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop