Prediction of Corn Zea mays L. Phenology Based on Cardinal Temperature Estimation, Spline Interpolation, and Numerical Analysis

Izzar Suly Nashrudin, Nur Afiiqotul Ula, Imro Atul Khoir

Abstract


Accurate crop phenology prediction is essential for modern agricultural management, irrigation scheduling, and climate change adaptation. This study develops a numerical-analysis-based framework to predict maize (Zea mays L.) growth stages using daily meteorological data. The proposed workflow integrates: (i) the non-linear Wang–Engel formulation to compute daily thermal units, (ii) cubic spline interpolation for data reconstruction under a missing-data validation scenario, (iii) Simpson’s 3/8 rule for numerical integration of cumulative thermal units, (iv) the central difference method to analyze the accumulation-rate dynamics, and (v) Taylor series expansion for local approximation of the Wang–Engel function around the optimum temperature. Daily meteorological data were obtained from the Open-Meteo Historical API for Jakarta, Indonesia in 2025, comprising 348 observation days. Numerical integration yields a cumulative thermal unit of 112.37 over the first 120 days. Derivative analysis identifies the maximum accumulation rate of 0.9784 per day at day 44. Using the adopted thermal thresholds, the model predicts the V3 stage at day 127 and the V6 stage at day 342. Furthermore, the second-order Taylor approximation attains a maximum error of approximately 8.4 × 10⁻³ within a ±7°C range around the optimum, while the third-order approximation reduces the error to the order of 10⁻⁴ over the tested range. This numerically robust framework can be extended for future integration with machine learning approaches.

Keywords


crop phenology; numerical analysis; growing degree days; spline interpolation; simpson integration; indonesian meteorology

Full Text:

PDF

References


  1. [1] F. Liu, C. R. Jensen, A. Shahanzari, M. N. Andersen, and S.-E. Jacobsen,
    "Aba regulated stomatal control and photosynthetic water use efficiency of potato (Solanum tuberosum L.) during progressive soil drying,"
    Plant Science, vol. 168, no. 3, pp. 831–836, 2005.
    doi: https://doi.org/10.1016/j.plantsci.2004.10.016.

    [2] B. Jamsheed, J. A. Wani, and A. Rehman,
    "Estimation of yield, phenology and agro-meteorological indices of quality protein maize (Zea mays L.) under different growing environments,"
    Journal of Krishi Science and Sustainable Systems, vol. 5, no. 3, pp. 1–12, 2023.

    [3] Y. Liu, Y. Qin, Q. Ge, J. Dai, and Q. Chen,
    "Responses and sensitivities of maize phenology to climate change from 1981 to 2009 in Henan Province, China,"
    Journal of Geographical Sciences, vol. 27, no. 9, pp. 1072–1084, 2017.
    doi: 10.1007/s11442-017-1422-4.

    [4] G. Zhou and Q. Wang,
    "A new nonlinear method for calculating growing degree days,"
    Scientific Reports, vol. 8, p. 10149, 2018.
    doi: 10.1038/s41598-018-28392-z.

    [5] I. M. Saiyed, P. R. Bullock, H. D. Sapirstein, G. J. Finlay, and C. K. Jarvis,
    "Thermal time models for estimating wheat phenological development and weather-based relationships to wheat quality,"
    Canadian Journal of Plant Science, vol. 89, no. 3, pp. 429–439, 2009.
    doi: 10.4141/CJPS07114.

    [6] C. F. Shaykewich,
    "An appraisal of cereal crop phenology modelling,"
    Canadian Journal of Plant Science, vol. 75, no. 2, pp. 329–341, 1995.
    doi: 10.4141/cjps95-057.

    [7] A. A. Sadiq, G. E. Njiddari, and P. A. Musa,
    "An overview of growing degree days and heat units application in crop production,"
    International Journal of Agriculture and Earth Science, vol. 11, no. 8, pp. 83–102, 2025.
    doi: 10.56201/ijaes.vol.11.no8.2025.pg83.102.

    [8] L. Roth, H.-P. Piepho, and A. Hund,
    "Phenomics data processing: Extracting dose–response curve parameters from high-resolution temperature courses and repeated field-based wheat height measurements,"
    in silico Plants, vol. 4, no. 1, diac007, May 2022.
    doi: 10.1093/insilicoplants/diac007.

    [9] Open-Meteo,
    "Open-Meteo Weather API Database," Registry of Open Data on AWS, 2023.
    Accessed: 15 Dec 2025.

    [10] E. Wang and T. Engel,
    "Simulation of phenological development of wheat crops,"
    Agricultural Systems, vol. 58, no. 1, pp. 1–24, 1998.
    doi: https://doi.org/10.1016/S0308-521X(98)00028-6.

    [11] N. A. Streck, I. Lago, L. F. Gabriel, and F. K. Samboranha,
    "Simulating maize phenology as a function of air temperature with a linear and a nonlinear model,"
    Pesquisa Agropecuária Brasileira, vol. 43, no. 4, pp. 449–455, 2008.
    doi: 10.1590/S0100-204X2008000400002.

    [12] M. Vrugt and B. Merkus,
    "Quantifying the relationship between thermal indices and maize phenology,"
    International Journal of Agriculture and Food Science, vol. 5, pp. 134–137, 2023.
    doi: 10.33545/2664844X.2023.v5.i2b.181.

    [13] K. Djaman, S. Irmak, W. Rathje, D. Martin, and D. Eisenhauer,
    "Maize evapotranspiration, yield production functions, biomass, grain yield, harvest index, and yield response factors under full and limited irrigation,"
    Transactions of the ASABE, vol. 56, pp. 373–393, 2013.
    doi: 10.13031/2013.42676.

    [14] P. Paredes et al.,
    "Estimating the lengths of crop growth stages to define the crop coefficient curves using growing degree days (GDD): Application of the revised FAO56 guidelines,"
    Agricultural Water Management, vol. 319, p. 109758, 2025.
    doi: https://doi.org/10.1016/j.agwat.2025.109758.

    [15] L. J. Abendroth, R. W. Elmore, M. J. Boyer, and S. K. Marlay,
    Corn Growth and Development. Ames, Iowa: Iowa State University Extension and Outreach, 2011.

    [16] M. Jamhuri, M. I. Irawan, I. Mukhlash, M. Iqbal, and N. N. T. Puspaningsih,
    "Neural networks optimization via Gauss–Newton based QR factorization on SARS-CoV-2 variant classification,"
    Systems and Soft Computing, vol. 7, p. 200195, 2025.

    [17] M. Jamhuri, S. P. Sari, S. Amiroch, J. Juhari, and V. A. Fitria,
    "Inexact generalized Gauss–Newton–CG for binary cross-entropy minimization,"
    Jurnal Riset Mahasiswa Matematika, vol. 5, no. 2, pp. 102–122, 2025.

    [18] C. R. Chikwendu, H. K. Oduwole, and S. I. Okoro,
    "An application of spline and piecewise interpolation to heat transfer (cubic case),"
    Mathematical Theory and Modeling, vol. 5, pp. 28–39, 2015.

    [19] M. Mulyono,
    "Evaluasi dari metode: Trapesium, Simpson 1/3, Simpson 3/8 dan Newton-Cotes orde 4–10 untuk menghitung integral tertentu secara numerik,"
    AKSIOMA: Jurnal Matematika dan Pendidikan Matematika, vol. 13, pp. 466–479, 2022.
    doi: 10.26877/aks.v13i3.12908.

    [20] MathForCollege.com,
    "Simpson’s 3/8 Rule for Integration," Online textbook chapter, MathForCollege.com, 2025.
    Accessed: 10 Dec 2025.

    [21] M. Jamhuri, M. I. Irawan, A. Kusumastuti, K. C. Mondal, and J. Juhari,
    "On the approximation capabilities of deep neural networks for multivariate time series modeling,"
    CAUCHY: Jurnal Matematika Murni dan Aplikasi, vol. 10, no. 2, pp. 1401–1416, 2025.




DOI: https://doi.org/10.18860/jrmm.v5i3.39988

Refbacks

  • There are currently no refbacks.