Probabilistic assimilation of optical satellite data with physiologically based growth functions improves crop trait time series reconstruction

Abstract

A sound understanding of plant growth is critical to maintaining future crop productivity under ongoing climate change. Remotely sensed time series of crop functional traits from optical satellite imagery are an invaluable tool for deriving appropriate management practices that facilitate risk mitigation and increase the resilience of agroecosystems. However, the availability of imagery is limited by atmospheric disturbances that cause large temporal gaps and noise in the trait time series. Therefore, time series reconstruction methods are required for accurate crop growth modelling. Physiological priors, such as the fact that plant growth is mainly controlled by a few environmental covariates, among which air temperature plays a prominent role, represent a promising approach to improve the representation of crop growth. Here, a novel approach is proposed that combines Sentinel-2 Green Leaf Area Index (GLAI) observations with three dose response curve approaches describing the a priori physiological relationship between growth and temperature in winter wheat. A probabilistic ensemble Kalman filtering data assimilation scheme allows the combination of high temporal resolution air temperature data and satellite imagery, which also allows quantification of uncertainties. The proposed approach requires a smaller number of satellite observations compared to conventional remote sensing time series algorithms, making it suitable for agricultural areas with high cloud cover, and is considerably less complex than a mechanistic crop growth model. Validation was carried out using in-situ data collected on winter wheat plots in Switzerland in two consecutive years. The validation results suggest that the proposed assimilation of Sentinel-2 GLAI and temperature-response-based growth rates allows the reconstruction of physiologically meaningful GLAI time series. In particular, the systematic underestimation of high in-situ GLAI values (> 5 m^2 m^-2) often prevalent in purely remote sensing driven GLAI time series reconstruction was reduced. Thus, the proposed approach is advantageous compared to state-of-the-art remote sensing approach based on wide-spread logistic functions by means of physiological plausibility, fitting requirements and representation of high in-situ GLAI values. This has great potential to increase the reliability of remotely sensed crop productivity assessment.

Plant phenotyping relevance

衛星光学データと生理モデルを統合して作物GLAI時系列を再構築する手法を提案し、圃場データで検証しており、植物形質推定が研究の中心である。

Code and data availability

The authors explicitly state that code and data to reproduce the entire workflow (DRC fitting, Sentinel-2 GLAI assimilation, and validation) are publicly available on GitHub under GNU GPL v3.0. This is a paper-specific, public, actionable asset. Other URLs in the text are cited references or generic libraries (e.g., NL