Skip to main content

Methodological approach for predicting and mapping the phenological adaptation of tropical maize (Zea mays L.) using multi-environment trials



The phenological development of the maize crop from emergence through flowering to maturity, usually expressed as a rate (i.e. 1/duration), is largely controlled by temperature in the tropics. Maize plant phenological responses vary between varieties and quantifying these responses can help in predicting the timing and duration of critical periods for crop growth that affect the quality and quantity of seed. We used routine multi-environment trials data of diverse tropical maize varieties to: (1) fit 82 temperature dependent phenology models and select the best model for an individual variety, (2) develop a spatial framework that uses the phenology model to predict at landscape level the length of the vegetative and reproductive phases of diverse varieties of maize in different agro-ecologies. Multi-environment trial data of 22 maize varieties from 16 trials in Kenya, Ethiopia, and Sudan was analyzed and the Levenberg–Marquardt algorithm combined with statistical criteria was applied to determine the best temperature-dependent model.


The Briere model, which is not often used in plant phenology, provided the best fit, with observed and predicted days to flowering showing good agreement. Linking the model with temperature and scaling out through mapping gave the duration from emergence to maturity of different maize varieties in areas where maize could potentially be grown.


The methodology and framework used in the study provides an opportunity to develop tools that enhance farmers’ ability to predict stages of maize development for efficient crop management decisions and assessment of climate change impacts. This methodology could contribute to increase maize production if used to identify varieties with desired maturity for a specific agro-ecology in in the targeted regions.


The method of using temperature and time as factors to describe the development of plants and insects was proposed and has been continuously studied for over 300 years [1]. The concept began with the French naturalist [1] who suggested that temperature variations were probably one of the causes of the changes in plant and insect phenology. They proposed summing up the mean daily air temperature for an identical number of months in a location, which led to the birth of the degree-day concept as the values of the developmental rate of plants obtained from year to year were roughly constant [1]. The degree-day equation was later modified by several authors [2,3,4]. As the concept of degree-day evolved, a new hypothesis arose stating that the rate of chemical reactions is either doubled or tripled for each 10 °C rise in temperature [5]. This idea was expressed in form of a constant called Q10 [5]. Towards the end of the 18th century, the principle started to influence the formulation of plant and insect development rate functions [5].

Mathematical, statistical and machine learning approaches are commonly used to enhance the understanding of plant architecture, growth, development, and interactions with the environment. A single-parameter compartmental model was used to describe the transport of fluoride in living plants [6]; and an ensemble of statistical metrics was applied to quantify large sets of plant transcription factor binding sites [7]. In Qiongyan et al. [8] a neural network was developed to detect the spikes of wheat plants, and an algorithm with deep learning method was applied for counting the leaf of rosette plants [9]. Increasingly, the formulation of mathematical expressions to represent the development of plants and insects continued to progress centered on the postulation that “within certain ranges, as temperature (Te) decreases, the rate of development slows and, if the temperature drops low enough, development will cease at the organism’s lower developmental threshold; often called the base temperature (Tb). As temperature increases, development rate increases until temperature reaches the optimum temperature (To), above which the development rate decreases and eventually ceases at a value called maximum temperature (Tm)”. Inspired by this hypothesis, various models have been developed, each with strengths and weaknesses [10]. Largely, these models derived the relationship between the development of plants and insects and temperature either empirically or through process-based methodology. Empirical functions are formulated and parameterized using the same measured phenomenon as the phenomenon to be derived [11]. In contrast process-based functions are mathematical expressions formulated and parameterized using biological knowledge such as the enzyme kinetics to predict the time from planting to anthesis [11].

Climate has a considerable impact on the distribution and abundance of plants and other organisms, and the mathematical depiction of the climatic effect on crop phenology has been of significant interest among scientists [12, 13]. Temperature is the most important climate variable determining plant phenology (via rate of development) and plant distribution (via base and maximum temperature limits to survival) [14, 15]. Understanding the phenology of maize can therefore help to define crop adaption to a region or site (i.e. its ability to mature and set seed within a growing season). In maize, it can also help hybrid seed production by determining appropriate planting dates of lines to ensure flowering synchrony [15, 16]. Therefore, the accurate prediction of phenological development is fundamental to define crop adaption and yield potential [15, 17].

Despite the abundance of literature [15, 18, 19] on the use of temperature dependent models for the prediction of maize developmental phases, few studies have been conducted on tropical maize in sub Saharan Africa (SSA). Furthermore, in order to get the best possible prediction of phenology for individual varieties across space and time, a range of well-known and well-tested temperature-dependent models for plant and insect development should be evaluated and made available through an open-source modeling framework. If this framework can be based on existing and annually or routinely conducted breeding or variety-testing trials, then the process can become institutionalized. Models can then be used to predict the adaptation of different maize varieties spatially and temporally across major maize growing agro-ecologies, information that is important in efforts to improve maize yields. The objectives of this study were therefore to: (1) to evaluate 82 different temperature dependent models to predict the phenological development of 22 diverse tropical maize varieties using existing multi-environment trial data; and (2) develop a spatial framework that uses the phenology model to predict the period from sowing to flowering or maturity across sites and agro-ecologies, i.e. to map the adaptation of different varieties.

Materials and methods

Maize varieties and test locations

Twenty-two (22) early and intermediate maturity open-pollinated maize varieties were planted in nine different locations in Kenya, Ethiopia, and Sudan in 2004 (a total of 16 trials) as described in Table 1 and Fig. 1 (adapted from [16]). In each of these trials the date when 50% of the plants in a plot: (1) emerged; (2) had a tassel (male); (3) had a silk (female) and (4) were physiologically mature (black layer formation) was recorded. Selection of experimental field data for model development and calibration was conducted using altitudinal change (400–1600 m) to represent a wide range of temperature values (11–20 °C minimum; 24–37 °C maximum). Six sites covering the maximum range of temperature were selected for model development and calibration, namely Alupe, Embu, Kibos and Bungoma, Pawe, Wad Madani. The remaining three sites, namely Kiboswa, Nyahera and Vihiga, were used for independent model evaluation. In model development about 1/3 of data are typically used for independent validation.

Table 1 Description of the locations where data used to develop, calibrate and validate the models for the selected maize varieties was collected as described in [16]
Fig. 1
figure 1

Locations of early and intermediate open-pollinated maize varieties trials conducted in Kenya, Ethiopia and Sudan

Temperature datasets and shape file of suitable soil and weather conditions

Two types of temperature data were used, daily temperatures for model development and monthly gridded temperature for mapping. Daily minimum and maximum temperature were obtained from weather stations in the vicinity of field trials [16]. Maize phenology data were collected during two growing seasons; March to August 2004 characterized by higher rainfall, and September 2004 to January 2005 characterised by lower rainfall. The data were used to develop the phenology model. Gridded average monthly minimum and maximum temperature for Kenya were obtained from WorldClim ( and the Climate Change, Agriculture and Food Security (CCAFS) ( databases. The data are organized in layers (grids) with a spatial resolution of 30 s that is equal to 0.9 × 0.9 km.

The shape file of suitable soil and weather conditions for maize production used in this study was obtained from HarvestChoice [20]. It is a spatially disaggregated production statistics database of derived from the Spatial Production Allocation Model (SPAM). The resolution is 5 arc-minute grid cells [20].


  1. 1.

    Temperature and time were considered as the primary factors affecting the rate of maize development. Photoperiod or daylength can also affect development (maize is a quantitative short-day plant), but breeders select against this trait and most tropical maize varieties are effectively insensitive to photoperiod. Furthermore, Kenya straddles the equator and the variation in photoperiod is very small [21, 22].

  2. 2.

    It was assumed that the trials were grown with adequate nutrients and water and did not suffer adverse stress based on the yields obtained [16]. In general phenological development is only adversely affected by very severe stress [16].

  3. 3.

    Maize plant developmental cycle was divided into two phases: the vegetative phase (VP), starting from emergence to tasseling/silking (flowering), and the reproductive phase (RP), which is from flowering to physiological maturity (PM). It was assumed that the duration of the VP and RP were the same in any given variety based on previous experience. The overall developmental rate per phase for the whole crop was estimated by accumulating the daily development rate values:

    $$D = \sum\limits_{0}^{t} {r(T(t),P)\Delta t}$$

    where development D is a function of temperature T which in turn is a function of time t, r is the development rate, Δt is the time increment and P the vector of parameter values [23]. The time step for the model was fixed to a day.

  1. 4.

    Phenology was predicted spatially only in areas where maize is known to be grown, based on suitability maps and known maize areas.


According to [24] development is the process of cellular differentiation manifested by different life-phases of an organism. Development time is the duration between life-phases and growth is the increase in biomass resulting from the development. Selecting a mathematical expression for representing the phase of development of a crop is challenging [25]. This study applied 82 temperature-dependent nonlinear equations (Additional file 1: Table S1) which have been used in agricultural production, either in the context of insect phenology modeling or crop development [25, 26].

Estimation of the model parameters

An important step in conducting this study was to evaluate how well the selected models (out of the 82 tested) fit the observed data. We used the Levenberg–Marquardt (LM) algorithm [27] because it combines both the ‘steepest descent’ and the Gauss–Newton method to evaluate the models. This algorithm iteratively locates the minimum of a function that is expressed as the sum of squares of the nonlinear model [27]. After initialization of the models, the LM algorithm was launched and a goodness of fit procedure was used to find the parameters of the model. The model parameters were estimated by fitting equations to the recorded data from multi-environment trials. A program written in R computing language [28] provided an interactive processing, in which initial values of the model parameter are entered; and once the LM algorithm is launched, the values of the parameters are optimized based on the matching of model output to observed data following the ‘steepest descent’ and the Gauss–Newton methods. Overall, the process is similar to the functioning of LEAF-E; a tool developed to analyze grass leaf growth using function fitting [29].

Goodness of fit and selection criteria of the model

No single method exists to best assess the goodness of fit of an individual model to specific data. Therefore visual assessment and the coefficient of determination R2_Adj [25, 30] were used. The R2_Adj is estimated as follows:

$$R^{2} \_Adj = 1 - \left( {1 - R_{z}^{2} } \right)\left( {\frac{n - 1}{n - k - 1}} \right)$$

where, n is the number of observations, k is the number of parameters of the zth sub-model, and \(R_{z}^{2}\) is the correlation coefficient obtained by the following expression:

$$R_{z}^{2} = 1 - \sum\limits_{j = 1}^{n} {\left( {y_{j} - \hat{y}_{jz} } \right)^{2} - } \sum\limits_{j = 1}^{n} {\left( {y_{j} - \bar{y}} \right)^{2} }$$

where, \(\bar{y}\) is the observed median, y the observed mean and \(\hat{y}_{ji}\) is the jth predicted value from the ith function.

The best-fitted model is selected by examining the residuals and comparing Akaike’s Information Criterion (AIC) and the Model Selection Criterion (MSC) [25]. The mathematical expressions of the two criteria are:

$$AIC = n\ln - \left( {\sum\limits_{i = 1}^{n} {w_{i} \left( {Y_{{obs_{i} }} - Y_{{est_{i} }} } \right)^{2} } } \right) + 2p$$
$$MSC = \ln \left( {\frac{{\sum\nolimits_{i = 1}^{n} {w_{i} \left( {Y_{{obs_{i} }} - \bar{Y}_{est} } \right)^{2} } }}{{\sum\nolimits_{i = 1}^{n} {w_{i} \left( {Y_{{obs_{i} }} - \bar{Y}_{{est_{i} }} } \right)^{2} } }}} \right) - \frac{2p}{n}$$

where, n is the number of observations, \(Y_{{obs_{i} }}\) and \(Y_{{est_{i} }}\) are the observed and estimated values for the ith observation, p is the number of parameters, and wi is the weight required for each observation. In the case of nested models, F test was used to check whether the addition of parameters has a statistically significant contribution to the model before selection [25].

Best-fit model: Briere et al. [31]

Although all the 82 models were fitted with the available multi-environment trial datasets of different maize varieties, herein we report only two models, which provided the best fit to the selected maize varieties based on criteria outlined above. The full set of models is available in a in the Additional file 1: Table S1.

The Briere model was developed based on four principles: (1) the option to estimate the lower and upper temperature thresholds; (2) the inclusion of an asymmetry in relation to the optimum value of temperature where the developmental rate is highest; (3) the presence of an inflection point; and (4) the introduction of a decay in development rate at temperatures above the optimum temperature [31]. With these considerations, the lower and the upper thermal limits were unequivocally integrated in the equation to represent two important parameters of the model, Tb and Tm, respectively [31]. To get decay at high temperatures, a square root was included to allow a high slope when the values of temperatures approach Tm. By combining the products of different powers of temperature, an inflection point occurs yielding the Briere_1 model with the following mathematical expression.

$$r(T) = a*T\left( {T - T_{o} } \right)\left( {\sqrt {T_{max} - T} } \right)$$

A second model, Briere_2, was derived from Briere_1 model by replacing the square root with a general power equal to d = 1/µ, where µ is the new parameter and a is an empirical constant [31]. The Briere_2 equation is as follows:

$$r(T) = a*T\left( {T - T_{o} } \right)\left( {\sqrt {T_{max} - T} } \right)^{d}$$

Evaluation of the model

The evaluation of the model in this study was defined as the level to which the selected model correctly predicted the number of days a maize variety takes from emergence to flowering and from flowering to physiological maturity at different, independent locations. To carry out this procedure, independent data recorded at three typical sites in Kenya (Kiboswa, Nyahera and Vihiga) were used. Flowering and maturity dates were predicted for a standard and common sowing date for each variety from the models developed using data from the other six sites and observed and predicted values compared.

Spatial predictions of the duration of maize phenology

A Kenya boundaries shape file was divided into grids or cells. To predict the number of days each phase of maize development takes to be completed in an individual cell, the temperature-dependent mathematical expression obtained during the modeling step was run at each individual cell of Kenya using a typical sowing date. The gridded temperature datasets were loaded into the computer, simultaneously extracted from the database and then organized in matrix format using longitude as column and latitude as a row. A point object picks the temperature-dependent mathematical expression of the maize vegetative and reproductive phases, and these are consecutively applied in each geographical coordinate of the grid. After replacing the values of temperature in the phenology model in each grid, a new matrix with the values of the development rate in the respective geographical coordinates was computed. These values are then inversed (1/development rate) to estimate the number of days used to complete each phase of maize development. The results are converted into ASCII files and transferred to an open source software Q-GIS [32] for visualization.


The multi-environment datasets were used to fit a total of 82 temperature driven models (Additional file 1: Table S1). However, only the model of Briere et al. [31] was selected to calculate the duration of vegetative and reproductive phases of 22 different maize varieties in Kenya. In addition to describing the data accurately (R2_Adj and R2 > 0.82 for all tested varieties), the Briere model was able to provide biologically meaningful estimates of Tb and Tm for development of 8–10 °C, and 38–40 °C, respectively.

Parameter estimates from the Briere_1 and Briere_2 models of eight-selected maize varieties are presented in Table 2. The estimated Tb ranged from 8.1 to 10.1 °C while Tm ranged from 34.3 to 40.4 °C. The adjusted R2 was high for all varieties (0.83–0.91). The root mean square error (RMSE) for each variety’s development rate ranged from 0.0014 to 0.0022, which further indicated the ability of the model to predict the duration of individual maize variety’s phases of phenology. Although the selected model did not explicitly estimate the optimum temperature at which the developmental rate is highest, graphical visualization (Fig. 2) showed that these values were between 23 and 32 °C, consistent with the literature, depending on the variety.

Table 2 Estimated parameters of Briere—1 and Briere—2 models [23] for eight maize varieties covering the vegetative and reproductive phases of the crop development
Fig. 2
figure 2

Temperature dependent models of development rate (1/development time) for maize varieties VE 201, VE 203, VE 206, VE 208, VE 210, VE 212, VE218, and VE 220 respectively

The number of days taken from emergence to flowering, and from flowering to physiological maturity, for the 22 maize varieties at three independent data sets (locations) used for model evaluation is presented in Table 3. The range in observed flowering durations was between 48 and 80 across sites and varieties. Comparison of the duration obtained from model outputs with field trial results showed good agreement (r = 0 .892**, Adj R2 = 0.786 for Kiboswa; r = 0.712**, Adj R2 = 0.482 for Nyahera; and r = 0.920**, Adj R2 = 0.838 for Vihiga) suggesting that the model predicted phenology well and could be used with confidence to predict phenology spatially. The number of days from sowing to emergence is approximately 7–10 days and also varies with temperature. Adding this to the duration of the vegetative and reproductive phases provides gives the length for the growing period of the 22 varieties of maize in any location.

Table 3 Observed and predicted number of days taken by each maize variety during the vegetative and reproductive phases of the crop development at three locations used for model evaluation

As expected, maize development is predicted to take longer in the cooler mid-altitudes of Western Kenya and to be shorter in the warmer lowlands of the coastal areas. Out of the 22 varieties of maize analyzed here, the longest growing period from sowing to flowering was for VE 212, where most predicted values were > 105 days. The earliest variety was VE 220 at < 90 days in the majority of regions in Kenya. These maps show that all selected varieties of maize can be grown in most parts of Kenya, with a range of adaptation among varieties available (Figs. 3, 4). The phenology models for each variety were mapped spatially and then filtered or masked to represent areas where soil and weather conditions are suitable for growing maize (Fig. 5 for selected varieties).

Fig. 3
figure 3

Maps of Kenya showing the number of days 10 maize varieties (VE 212, VE 213, VE 214, VE 215, VE 216, 217, VE 218, VE 219, VE 220 and VE 221) would take during each phase (vegetative or reproductive) of development in different agro-ecologies. The small white dots are areas with missing values of temperature

Fig. 4
figure 4

Maps of Kenya showing the number of days 10 maize varieties (VE 201, VE 203, VE 206, VE 208, VE 210, VE 212, VE 218, and VE 220) would take during each phase (vegetative or reproductive) of development in different agro-ecologies. This figure was obtained by masking the distribution map of suitable soil and weather conditions for growing maize with the phenology maps presented in Fig. 3 and in this figure

Fig. 5
figure 5

Maps of Kenya showing the number of days 12 maize varieties (VE 201, VE 202, VE 203, VE 204, VE 205, VE 206, VE 207, VE 208, VE 209, VE 210 and VE 211) would take during each phase (vegetative or reproductive) of development in different agro-ecologies. The small white dots are areas with missing values of temperature


This study had two major objectives, namely to develop a model framework to rigorously predict maize phenology from routine yield trials and to map adaptation spatially. Maize breeders conduct many yield trials each year in a wide range of locations representing target environments as well as ‘hot-spots’ for biotic or abiotic stresses [16]. In most, if not all, trials phenology, a key trait, is recorded. These standard data sets are available for predictive modeling and, as we have shown, are resources that can be successfully used to add value to the breeding process.

The maize plant requires a certain amount of heat units or thermal energy to transition between different phases of development [13, 33]. Because of yearly and within season variation in weather patterns, measuring the heat units accumulated over time provides a physiological time scale that is biologically more accurate than calendar days for predicting stages of development [13, 33]. Phenology models, characterized by a rate of development, are often applied to predict the timing of events in plant development. However, it is difficult to directly measure the rate of plant development [4]. Usually a relationship is established between the development rate (development rate is calculated as the inverse of development time), and the development time as the latter (e.g. flowering date) can be easily and routinely measured from field studies. The most critical step in the current study was the choice of a model among numerous thermal mathematical models that exist to describe the developmental response of plants and insects to temperature. No standard method exists to choose between competing models and the modeling framework developed here provides a solution. The decision on choice of a model should not be based on statistics alone but a combination of statistics and scientist’s experience and skills to identify an appropriate model that provides the best biologically meaningful and statistically significant parameters. The model of Briere et al. [23] described well the development rate of the selected maize varieties. This finding was interesting because, to the best of our knowledge, the Briere model has not been used in the context of plant development. The model has always been applied to describe the temperature dependent development rate of insect life stages [31]. The model provided estimates of number of days from emergence to the end of the vegetative phase that closely matched the recorded data at the three locations. This suggested that this model is suitable for modeling growth of tropical maize. The model has several features: it requires a small number of parameters making it easy to fit and the formulation allows the global minimum of the loss function to be rapidly reached. The parameters Tb and Tm have biological definitions that provide smooth curves with adequate approximation of biological processes. The curve captured maize response to temperature within a large range of temperature values. The selected model includes numerous features such as the dependency at low temperatures, the positive linear dependence at intermediate temperatures, the parabola response across the optimum temperatures and the negative linear dependency at high temperatures [10, 34]. Many of these features are not found among the eight temperature dependent models used for predicting the phenology of maize as described in [11]. This study showed that the accuracy of the models applied was mainly associated with the temperature response across a selected range rather than Tb and Tm. Herein we argue that such association may not always be adequate as the physiological process of maize crop development is illustrated by gradual variation of organs and tissues stemming in gradients [19]. During maize development, the temperature responses are absolute and not relative; therefore shifting temperature beyond a set limit can stimulate an immediate effect, which may not be evident if the value of the temperature was unchanged or remained below or above the threshold. In other words the threshold value of temperature is an important variable for accurately predicting maize phenology and should not be omitted but carefully estimated through the smoothness of the model we have proposed [19].

Although numerous efforts have been made by researchers to suggest mathematical expressions to describe the relationship between the timing of events during maize development and environmental temperature [15, 18], few attempts were made on the mapping of phenological development in SSA. Simulation of the phenological development of wheat and maize at global scale was conducted by [15]. The authors applied the concept of heat units, which assumes the accumulation of daily temperature above Tb. Because the model was primarily developed for the temperate regions, the equation included an incremental modification of the heat units due to the effect of photoperiod. Using crop calendars and large-scale pattern of phenological characteristics of varieties, a spatial prediction of the length of growing season for each crop was produced. The results revealed an over and under estimation by 0.5–1.5 months of the duration of cropping period of maize [15]. Such outcome may be explained by the fact that a linear model (degree-day), which lacks the ability to capture the nonlinearity in the processes that govern plant development, was used. In addition, the model was projected globally; ignoring that in tropical regions the effects of the sensitivity to day length (photoperiod) is negligible on maize phenology. We suggest that the calibration of model should be conducted for a defined area of interest, as a single model is unlikely to work everywhere.

Taking into account that maize crop is grown across large agro-ecologies in Kenya, and elsewhere in Africa, it was important in this study to adopt a framework that could provide some level of knowledge at agro-ecological and national level. The mapping framework used readily available open-access temperature (and other variables if desired) data and could be used for other crops and research questions. Therefore, to scale up our results, a bottom-up approach in which a number of sites were selected for the analysis of the model followed by mapping was adopted. During the mapping step the temperature applied contained wider ranges than those used at individual trial locations, which helped to estimate the length of growing period for each variety of maize in the whole of Kenya. Such outcome is of high importance to farmers as it provides location specific and accurate duration of events during the growing of maize, which could guide in management planning. The predictive mapping framework could be used at a range of scales depending on the user and their needs. For example, extension workers may want to know what variety performs best at a particular location. Input dealers, on the other hand, may want to know what range of varieties to stock for their market area.

Implications of the study and conclusions

Phenology models are important analytical tools for predicting, evaluating, and understanding the length of crop growing duration under diverse environmental conditions. Linking the model with climate drivers such as temperature, and scaling out through predictive mapping, permitted us to estimate the phase duration of specific maize varieties throughout areas where maize is and can be grown in Kenya. Furthermore, these predictive models can be developed for other species and also used with future climate scenarios, to assess climate change impacts [35] and future breeding targets.

In East Africa, information on maize seed packages is usually very limited; with varieties being categorized in imprecise altitude classes (low, medium, high) and/or maturity groups (early, medium, late), neither of which are predictive. This study provided predictive, and hence location specific length of growing periods for 22 maize varieties. This information can help farmers in two ways. First, information on when varieties will mature can help farmers select the best variety for their location and growing season. Second, phonological predictions can be used to predict the timing of key stages of development that can be used to provide advice on the timing key events such weeding and application of fertilizer that will lead to improved yields. In future we intend to compact the maps and include additional features for each variety of maize (name, type, color, potential yield, maturity class, ecology, resistance to diseases/pests, tolerance to abiotic stresses, year and company of release) to be provided by breeders and specific country variety release institutions to develop a tool and smart phone applications that we hope will contribute to increased maize production in sub-Saharan Africa.


  1. Egerton FN. A history of the ecological sciences, part 21: Réaumur and his history of insects. Ecol Soc Am. 2006;87:212–24.

    Article  Google Scholar 

  2. Wang JY. A critique of the heat unit approach to plant response studies. Ecology. 1960;41:785–90.

    Article  Google Scholar 

  3. Allen JC. A modified sine wave method for calculating degree days. Environ Entomol. 1976;5:388–96.

    Article  Google Scholar 

  4. DeGaetano AT, Knapp WW. Standardization of weekly growing degree day accumulations based on differences in temperature observation time and method. Agric For Meteorol. 1993;66:1–19.

    Article  Google Scholar 

  5. Montagnes DJ, Kimmance SA, Atkinson D. Using Q10: can growth rates increase linearly with temperature? Aquat Microb Ecol. 2003;32:307–13.

    Article  Google Scholar 

  6. Converse AK, Ahlers EO, Bryan TW, Hetue JD, Lake KA, Ellison PA, et al. Mathematical modeling of positron emission tomography (PET) data to assess radiofluoride transport in living plants following petiolar administration. Plant Methods. 2015;11:18.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Hosseini P, Ovcharenko I, Matthews BF. Using an ensemble of statistical metrics to quantify large sets of plant transcription factor binding sites. Plant Methods. 2013;9:12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Qiongyan L, Cai J, Berger B, Okamoto M, Miklavcic SJ. Detecting spikes of wheat plants using neural networks with Laws texture energy. Plant Methods. 2017;13:83.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Ubbens J, Cieslak M, Prusinkiewicz P, Stavness I. The use of plant models in deep learning: an application to leaf counting in rosette plants. Plant Methods. 2018;14:6.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Yan W, Hunt LA. An equation for modelling the temperature response of plants using only the cardinal temperatures. Ann Bot. 1999;84:607–14.

    Article  Google Scholar 

  11. Kumudini S, Andrade FH, Boote KJ, Brown GA, Dzotsi KA, Edmeades GO, et al. Predicting maize phenology: intercomparison of functions for developmental response to temperature. Agron J. 2014;106:2087–97.

    Article  Google Scholar 

  12. Willis CG, Ruhfel BR, Primack RB, Miller-Rushing AJ, Losos JB, Davis CC. Favorable climate change response explains non-native species’ success in Thoreau’s woods. PLoS ONE. 2010;5:e8878.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Mohammed AR, Tarpley L. Instrumentation enabling study of plant physiological response to elevated night temperature. Plant Methods. 2009;5:7.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Damos P, Savopoulou-Soultani M. Temperature-driven models for insect development and vital thermal requirements. Psyche J Entomol. 2012;2012:123405.

    Article  Google Scholar 

  15. Bussel LGJ, Stehfest E, Siebert S, Müller C, Ewert F. Simulation of the phenological development of wheat and maize at the global scale. Glob Ecol Biogeogr. 2015;24:1018–29.

    Article  Google Scholar 

  16. Makumbi D, Kanampiu F, Mugo S, Karaya H, et al. Agronomic performance and genotype × times × environment interaction of herbicide-resistant maize varieties in Eastern Africa. Crop Sci. 2015;55:540–55.

    Article  Google Scholar 

  17. Schröder W, Schmidt G, Schönrock S. Modelling and mapping of plant phenological stages as bio-meteorological indicators for climate change. Environ Sci Eur. 2014;26:5.

    Article  Google Scholar 

  18. Streck NA, Lago I, Gabriel LF, Samboranha FK. Simulating maize phenology as a function of air temperature with a linear and a nonlinear model. Pesqui Agropecuária Bras. 2008;43:449–55.

    Article  Google Scholar 

  19. Sánchez B, Rasmussen A, Porter JR. Temperatures and the growth and development of maize and rice: a review. Glob Change Biol. 2014;20:408–17.

    Article  Google Scholar 

  20. Maize Production (mt, 2005) | HarvestChoice

  21. Bonhomme R, Derieux M, Edmeades GO. Flowering of diverse maize cultivars in relation to temperature and photoperiod in multilocation field trials. Crop Sci. 1994;34:156–64.

    Article  Google Scholar 

  22. Gouesnard B, Rebourg C, Welcker C, Charcosset A. Analysis of photoperiod sensitivity within a collection of tropical maize populations. Genet Resour Crop Evol. 2002;49:471–81.

    Article  Google Scholar 

  23. Worner SP. Performance of phenological models under variable temperature regimes: consequences of the Kaufmann or rate summation effect. Environ Entomol. 1992;21:689–99.

    Article  Google Scholar 

  24. Wagner TL, Olson RL, Willers JL. Modeling arthropod development time. J Agric Entomol USA. 1991;8(4):251–70.

    Google Scholar 

  25. Archontoulis SV, Miguez FE. Nonlinear regression models and applications in agricultural research. Agron J. 2015;107:786–98.

    Article  Google Scholar 

  26. Tonnang EZH, Juarez H, Carhuapoma P, Gonzales JC, Mendoza D, Sporleder M, et al. ILCYM-insect life cycle modeling. A software package for developing temperature-based insect phenology models with applications for local, regional and global analysis of insect population and mapping. Int Potato Cent Lima Peru. 2013;193.

  27. Lampton M. Damping–undamping strategies for the Levenberg–Marquardt nonlinear least-squares method. Comput Phys. 1997;11:110–5.

    Article  Google Scholar 

  28. Ihaka R, Gentleman R. R: a language for data analysis and graphics. J Comput Graph Stat. 1996;5:299–314.

    Google Scholar 

  29. Voorend W, Lootens P, Nelissen H, Roldán-Ruiz I, Inzé D, Muylle H. LEAF-E: a tool to analyze grass leaf growth using function fitting. Plant Methods. 2014;10:37.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Wallach D, Makowski D, Jones JW. Working with dynamic crop models: evaluating, analyzing, parameterizing, and applications. New York: Elsevier; 2006.

    Google Scholar 

  31. Briere J-F, Pracros P, Le Roux A-Y, Pierre J-S. A novel rate model of temperature-dependent development for arthropods. Environ Entomol. 1999;28:22–9.

    Article  Google Scholar 

  32. Steiniger S, Hunter AJ. Free and open source GIS software for building a spatial data infrastructure. Geospat Free Open Sour Softw 21st Century. 2012;247–61.

  33. Yoo SY, Kim Y, Kim SY, Lee JS, Ahn JH. Control of flowering time and cold response by a NAC-domain protein in Arabidopsis. PLoS ONE. 2007;2:e642.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Shi P-J, Reddy GV, Chen L, Ge F. Comparison of thermal performance equations in describing temperature-dependent developmental rates of insects:(I) empirical models. Ann Entomol Soc Am. 2016;109:211–5.

    Article  Google Scholar 

  35. Xu H, Twine TE, Girvetz E. Climate change and maize yield in Iowa. PLoS ONE. 2016;11:e0156083.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

DM—Conceived and designed the experiments; HEZT—analyzed the data; PC—contributed to materials/analysis tools; HEZT, DM, PC—wrote the paper. All authors read and approved the final manuscript.


This work was implemented by CIMMYT as part of TAMASA (Taking Maize Agronomy to Scale) project in collaboration with IITA and IPNI, made possible by the generous support from Bill and Melinda Gates Foundation (contract OPP1113374). Any opinions, finding, conclusion, or recommendations expressed in this publication are those of the authors and do not necessarily reflect the view of the donor.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Henri E. Z. Tonnang.

Additional file

Additional file 1: Table S1.

Summary of the mathematical expressions of the 82 models tested in the study.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tonnang, H.E.Z., Makumbi, D. & Craufurd, P. Methodological approach for predicting and mapping the phenological adaptation of tropical maize (Zea mays L.) using multi-environment trials. Plant Methods 14, 108 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: