Data-driven forecasting of tropical cloud forest phenology: A dynamical systems and machine learning hybrid approach
Authors: Pedro Garcia, Saul Flores
Organizations: Departamento de Física, Facultad de Ingeniería, Universidad Católica Andrés Bello, Caracas, Venezuela · Laboratorio de Ecología de Suelos, Centro de Ecología, Instituto Venezolano de Investigaciones Científicas, Caracas, Venezuela
Forecasting how species and ecosystems respond to climate change is a critical challenge in ecology, particularly in hyper-diverse tropical biomes. While traditional phenological models often rely on preconceived functional forms, this study introduces a data-driven hybrid approach that integrates dynamical systems theory with machine learning to forecast flowering and fruiting cycles. Using a robust 35-year dataset from a Venezuelan tropical cloud forest, we reconstruct the system's state-space using delay vectors inspired by Takens' theorem, incorporating climate forcing variables such as temperature and precipitation. The model was trained using Gaussian Process Regression, providing both point predictors and explicit uncertainty estimates. Our results clearly suggest that the synergy between dynamical systems and machine learning effectively captures the non-linear complexities of ecological time series. Specifically, the model achieved higher predictive accuracy for flowering and fruiting. This hybrid framework offers a scalable tool for anticipating phenological mismatches and informing conservation strategies in the face of global environmental change.
Figures & tables
No.
Family
Species
GBIF link
Growth form
Reproductive density
(4m -2 )
1
Apocynaceae
Aspidosperma fendleri Woodson
3618540
Tree
0.25 ± 0.18
2
Burseraceae
Protium tovarense Pittier
3992944
Tree
0.29 ± 0.21
3
Phyllanthaceae
Hieronyma fendleri Briq.
3079076
Tree
0.20 ± 0.16
4
Phyllanthaceae
Richeria grandis Vahl
3079144
Tree
0.23 ± 0.27
5
Marcgraviaceae
Ruyschia tremadena (Ernst) Lundell
3915996
Climber
0.10 ± 0.01
Table 1: Focal study species.
Figure 1: State space of the system forest-clima.
Figure 2: Forecasting mean phenology.
Figure 3: GPR flowering and fruiting prediction for each specie vs. actual data with confidence interval.
Species
Flowering
Fruiting
A.fendleri
0.444
0.581
B.reticulata
0.503
0.495
C.tremadena
0.444
0.478
E.amazonicum
0.455
0.637
E.americanus
0.705
0.546
G.latifolia
0.449
0.482
Table 2: Root-mean-square error of GPR by species.
Figure 4: K-Means Clustering of species based on RMSE.
Time series forecasting predicts future values from ordered observations using statistical, neural, and probabilistic models to support decisions under uncertainty.
Accurate prediction of grape phenology is essential for timely vineyard management decisions, such as scheduling irrigation and fertilization, to maximize crop yield and quality. While traditional biophysical models calibrated on historical field data can be used for season-long predictions, they lack the precision required for fine-grained vineyard management. Deep learning methods are a compelling alternative but their performance is hindered by sparse phenology datasets, particularly at the cultivar level. We propose a hybrid modeling approach that combines multi-task learning with a recurrent neural network to parameterize a differentiable biophysical model. By using multi-task learning to predict the parameters of the biophysical model, our approach enables shared learning across cultivars while preserving biological structure, thereby improving the robustness and accuracy of predictions. Empirical evaluation using real-world and synthetic datasets demonstrates that our method significantly outperforms both conventional biophysical models and baseline deep learning approaches in predicting phenological stages, as well as other crop state variables such as cold-hardiness and wheat yield.
Accurate prediction of crop states (e.g., phenology stages and cold hardiness) is essential for timely farm management decisions such as irrigation, fertilization, and canopy management to optimize crop yield and quality. While traditional biophysical models can be used for season-long predictions, they lack the precision required for site-specific management. Deep learning methods are a compelling alternative, but can produce biologically unrealistic predictions and require large-scale data. We propose a hybrid modeling approach that uses a neural network to parameterize a differentiable biophysical model and leverages multi-task learning for efficient data sharing across crop cultivars in data limited settings. By predicting the parameters of the biophysical model, our approach improves the prediction accuracy while preserving biological realism. Empirical evaluation using real-world and synthetic datasets demonstrates that our method improves prediction accuracy by 60% for phenology and 40% for cold hardiness compared to deployed biophysical models.
William Solow, Paola Pesantez-Cabrera, Markus Keller +3
Oregon State University · Washington State University
Plant phenology modelling aims to predict the timing of seasonal phases, such as leaf-out or flowering, from meteorological time series. Reliable predictions are crucial for anticipating ecosystem responses to climate change. While phenology modelling has traditionally relied on mechanistic approaches, deep learning methods have recently been proposed as flexible, data-driven alternatives with often superior performance. However, mechanistic models tend to outperform deep networks when data distribution shifts are induced by climate change. Domain Adaptation (DA) techniques could help address this limitation. Yet, unlike standard DA settings, climate change induces a temporal continuum of domains and involves both a covariate and label shift, with warmer records and earlier start of spring. To tackle this challenge, we introduce Mid-feature Rank-adversarial Domain Adaptation (MIRANDA). Whereas conventional adversarial methods enforce domain invariance on final latent representations, an approach that does not explicitly address label shift, we apply adversarial regularization to intermediate features. Moreover, instead of a binary domain-classification objective, we employ a rank-based objective that enforces year-invariance in the learned meteorological representations. On a country-scale dataset spanning 70 years and comprising 67,800 phenological observations of 5 tree species, we demonstrate that, unlike conventional DA approaches, MIRANDA improves robustness to climatic distribution shifts and narrows the performance gap with mechanistic models.
Yuchang Jiang, Jan Dirk Wegner, Vivien Sainte Fare Garnot