Towards Optimising the Derivation of Phenological Phases of Different Crop Types over Germany Using Satellite Image Time Series

Abstract

Operational crop monitoring applications, including crop type mapping, condition monitoring, and yield estimation, would benefit from the ability to robustly detect and map crop phenology measures related to the crop calendar and management activities like emergence, stem elongation, and harvest timing. However, this has proven to be challenging due to two main issues: first, the lack of optimised approaches for accurate crop phenology retrievals, and second, the cloud cover during the crop growth period, which hampers the use of optical data. Hence, in the current study, we outline a novel calibration procedure that optimises the settings to produce high-quality NDVI time series as well as the thresholds for retrieving the start of the season (SOS) and end of the season (EOS) of different crops, making them more comparable and related to ground crop phenological measures. As a first step, we introduce a new method, termed UE-WS, to reconstruct high-quality NDVI time series data by integrating a robust upper envelope detection technique with the Whittaker smoothing filter. The experimental results demonstrate that the new method can achieve satisfactory performance in reducing noise in the original NDVI time series and producing high-quality NDVI profiles. As a second step, a threshold optimisation approach was carried out for each phenophase of three crops (winter wheat, corn, and sugarbeet) using an optimisation framework, primarily leveraging the state-of-the-art hyperparameter optimization method (Optuna) by first narrowing down the search space for the threshold parameter and then applying a grid search to pinpoint the optimal value within this refined range. This process focused on minimising the error between the satellite-derived and observed days of the year (DOY) based on data from the German Meteorological Service (DWD) covering two years (2019–2020) and three federal states in Germany. The results of the calculation of the median of the temporal difference between the DOY observations of DWD phenology held out from a separate year (2021) and those derived from satellite data reveal that it typically ranged within ±10 days for almost all phenological phases. The validation results of the detection of dates of phenological phases against separate field-based phenological observations resulted in an RMSE of less than 10 days and an R-squared value of approximately 0.9 or greater. The findings demonstrate how optimising the thresholds required for deriving crop-specific phenophases using high-quality NDVI time series data could produce timely and spatially explicit phenological information at the field and crop levels.

Plant phenotyping relevance

衛星時系列から作物のフェノロジー（SOS・EOS等）を抽出するNDVI再構成法と閾値最適化法を開発し、独立した圃場観測で検証しており、植物状態の取得手法が研究の中心である。

Code and data availability

No paper-specific public assets qualify. The Data Availability Statement states that phenological in-situ data are confidential/restricted and that all other data will be made available only upon request, with no public repository URL or deposited code/models provided anywhere in the supplied blocks.

No evidence-backed public reproduction asset is currently recorded.