Joint Assimilation of Leaf Area Index and Soil Moisture from Sentinel-1 and Sentinel-2 Data into the WOFOST Model for Winter Wheat Yield Estimation

Haizhu Pan, Zhongxin Chen*, Allard de Wit, Jianqiang Ren

*Corresponding author for this work

Research output: Contribution to journalArticleAcademicpeer-review

Abstract

It is well known that timely crop growth monitoring and accurate crop yield estimation at a fine scale is of vital importance for agricultural monitoring and crop management. Crop growth models have been widely used for crop growth process description and yield prediction. In particular, the accurate simulation of important state variables, such as leaf area index (LAI) and root zone soil moisture (SM), is of great importance for yield estimation. Data assimilation is a useful tool that combines a crop model and external observations (often derived from remote sensing data) to improve the simulated crop state variables and consequently model outputs like crop total biomass, water use and grain yield. In spite of its effectiveness, applying data assimilation for monitoring crop growth at the regional scale in China remains challenging, due to the lack of high spatiotemporal resolution satellite data that can match the small field sizes which are typical for agriculture in China. With the accessibility of freely available images acquired by Sentinel satellites, it becomes possible to acquire data at high spatiotemporal resolution (10-30 m, 5-6 days), which offers attractive opportunities to characterize crop growth. In this study, we assimilated remotely sensed LAI and SM into the Word Food Studies (WOFOST) model to estimate winter wheat yield using an ensemble Kalman filter (EnKF) algorithm. The LAI was calculated from Sentinel-2 using a lookup table method, and the SM was calculated from Sentinel-1 and Sentinel-2 based on a change detection approach. Through validation with field data, the inverse error was 10% and 35% for LAI and SM, respectively. The open-loop wheat yield estimation, independent assimilations of LAI and SM, and a joint assimilation of LAI + SM were tested and validated using field measurement observation in the city of Hengshui, China, during the 2016-2017 winter wheat growing season. The results indicated that the accuracy of wheat yield simulated by WOFOST was significantly improved after joint assimilation at the field scale. Compared to the open-loop estimation, the yield root mean square error (RMSE) with field observations was decreased by 69 kg/ha for the LAI assimilation, 39 kg/ha for the SM assimilation and 167 kg/ha for the joint LAI + SM assimilation. Yield coefficients of determination (R2) of 0.41, 0.65, 0.50, and 0.76 and mean relative errors (MRE) of 4.87%, 4.32%, 4.45% and 3.17% were obtained for open-loop, LAI assimilation alone, SM assimilation alone and joint LAI + SM assimilation, respectively. The results suggest that LAI was the first-choice variable for crop data assimilation over SM, and when both LAI and SM satellite data are available, the joint data assimilation has a better performance because LAI and SM have interacting effects. Hence, joint assimilation of LAI and SM from Sentinel-1 and Sentinel-2 at a 20 m resolution into the WOFOST provides a robust method to improve crop yield estimations. However, there is still bias between the key soil moisture in the root zone and the Sentinel-1 C band retrieved SM, especially when the vegetation cover is high. By active and passive microwave data fusion, it may be possible to offer a higher accuracy SM for crop yield prediction.

Original languageEnglish
Article number3161
JournalSensors (Basel, Switzerland)
Volume19
Issue number14
DOIs
Publication statusPublished - 18 Jul 2019

Fingerprint

leaf area index
wheat
soil moisture
Soil moisture
assimilation
food
moisture
winter
Triticum
Soil
Joints
Crops
Food
crops
crop growth
Growth
China
cyhalothrin
Satellites
Monitoring

Keywords

  • data assimilation
  • EnKF
  • LAI
  • Sentinel-1 and Sentinel-2
  • SM
  • winter wheat yield
  • WOFOST

Cite this

@article{d22189556bfe4b34b86ffa81d62bd8c4,
title = "Joint Assimilation of Leaf Area Index and Soil Moisture from Sentinel-1 and Sentinel-2 Data into the WOFOST Model for Winter Wheat Yield Estimation",
abstract = "It is well known that timely crop growth monitoring and accurate crop yield estimation at a fine scale is of vital importance for agricultural monitoring and crop management. Crop growth models have been widely used for crop growth process description and yield prediction. In particular, the accurate simulation of important state variables, such as leaf area index (LAI) and root zone soil moisture (SM), is of great importance for yield estimation. Data assimilation is a useful tool that combines a crop model and external observations (often derived from remote sensing data) to improve the simulated crop state variables and consequently model outputs like crop total biomass, water use and grain yield. In spite of its effectiveness, applying data assimilation for monitoring crop growth at the regional scale in China remains challenging, due to the lack of high spatiotemporal resolution satellite data that can match the small field sizes which are typical for agriculture in China. With the accessibility of freely available images acquired by Sentinel satellites, it becomes possible to acquire data at high spatiotemporal resolution (10-30 m, 5-6 days), which offers attractive opportunities to characterize crop growth. In this study, we assimilated remotely sensed LAI and SM into the Word Food Studies (WOFOST) model to estimate winter wheat yield using an ensemble Kalman filter (EnKF) algorithm. The LAI was calculated from Sentinel-2 using a lookup table method, and the SM was calculated from Sentinel-1 and Sentinel-2 based on a change detection approach. Through validation with field data, the inverse error was 10{\%} and 35{\%} for LAI and SM, respectively. The open-loop wheat yield estimation, independent assimilations of LAI and SM, and a joint assimilation of LAI + SM were tested and validated using field measurement observation in the city of Hengshui, China, during the 2016-2017 winter wheat growing season. The results indicated that the accuracy of wheat yield simulated by WOFOST was significantly improved after joint assimilation at the field scale. Compared to the open-loop estimation, the yield root mean square error (RMSE) with field observations was decreased by 69 kg/ha for the LAI assimilation, 39 kg/ha for the SM assimilation and 167 kg/ha for the joint LAI + SM assimilation. Yield coefficients of determination (R2) of 0.41, 0.65, 0.50, and 0.76 and mean relative errors (MRE) of 4.87{\%}, 4.32{\%}, 4.45{\%} and 3.17{\%} were obtained for open-loop, LAI assimilation alone, SM assimilation alone and joint LAI + SM assimilation, respectively. The results suggest that LAI was the first-choice variable for crop data assimilation over SM, and when both LAI and SM satellite data are available, the joint data assimilation has a better performance because LAI and SM have interacting effects. Hence, joint assimilation of LAI and SM from Sentinel-1 and Sentinel-2 at a 20 m resolution into the WOFOST provides a robust method to improve crop yield estimations. However, there is still bias between the key soil moisture in the root zone and the Sentinel-1 C band retrieved SM, especially when the vegetation cover is high. By active and passive microwave data fusion, it may be possible to offer a higher accuracy SM for crop yield prediction.",
keywords = "data assimilation, EnKF, LAI, Sentinel-1 and Sentinel-2, SM, winter wheat yield, WOFOST",
author = "Haizhu Pan and Zhongxin Chen and {de Wit}, Allard and Jianqiang Ren",
year = "2019",
month = "7",
day = "18",
doi = "10.3390/s19143161",
language = "English",
volume = "19",
journal = "Sensors",
issn = "1424-8220",
publisher = "MDPI",
number = "14",

}

Joint Assimilation of Leaf Area Index and Soil Moisture from Sentinel-1 and Sentinel-2 Data into the WOFOST Model for Winter Wheat Yield Estimation. / Pan, Haizhu; Chen, Zhongxin; de Wit, Allard; Ren, Jianqiang.

In: Sensors (Basel, Switzerland), Vol. 19, No. 14, 3161, 18.07.2019.

Research output: Contribution to journalArticleAcademicpeer-review

TY - JOUR

T1 - Joint Assimilation of Leaf Area Index and Soil Moisture from Sentinel-1 and Sentinel-2 Data into the WOFOST Model for Winter Wheat Yield Estimation

AU - Pan, Haizhu

AU - Chen, Zhongxin

AU - de Wit, Allard

AU - Ren, Jianqiang

PY - 2019/7/18

Y1 - 2019/7/18

N2 - It is well known that timely crop growth monitoring and accurate crop yield estimation at a fine scale is of vital importance for agricultural monitoring and crop management. Crop growth models have been widely used for crop growth process description and yield prediction. In particular, the accurate simulation of important state variables, such as leaf area index (LAI) and root zone soil moisture (SM), is of great importance for yield estimation. Data assimilation is a useful tool that combines a crop model and external observations (often derived from remote sensing data) to improve the simulated crop state variables and consequently model outputs like crop total biomass, water use and grain yield. In spite of its effectiveness, applying data assimilation for monitoring crop growth at the regional scale in China remains challenging, due to the lack of high spatiotemporal resolution satellite data that can match the small field sizes which are typical for agriculture in China. With the accessibility of freely available images acquired by Sentinel satellites, it becomes possible to acquire data at high spatiotemporal resolution (10-30 m, 5-6 days), which offers attractive opportunities to characterize crop growth. In this study, we assimilated remotely sensed LAI and SM into the Word Food Studies (WOFOST) model to estimate winter wheat yield using an ensemble Kalman filter (EnKF) algorithm. The LAI was calculated from Sentinel-2 using a lookup table method, and the SM was calculated from Sentinel-1 and Sentinel-2 based on a change detection approach. Through validation with field data, the inverse error was 10% and 35% for LAI and SM, respectively. The open-loop wheat yield estimation, independent assimilations of LAI and SM, and a joint assimilation of LAI + SM were tested and validated using field measurement observation in the city of Hengshui, China, during the 2016-2017 winter wheat growing season. The results indicated that the accuracy of wheat yield simulated by WOFOST was significantly improved after joint assimilation at the field scale. Compared to the open-loop estimation, the yield root mean square error (RMSE) with field observations was decreased by 69 kg/ha for the LAI assimilation, 39 kg/ha for the SM assimilation and 167 kg/ha for the joint LAI + SM assimilation. Yield coefficients of determination (R2) of 0.41, 0.65, 0.50, and 0.76 and mean relative errors (MRE) of 4.87%, 4.32%, 4.45% and 3.17% were obtained for open-loop, LAI assimilation alone, SM assimilation alone and joint LAI + SM assimilation, respectively. The results suggest that LAI was the first-choice variable for crop data assimilation over SM, and when both LAI and SM satellite data are available, the joint data assimilation has a better performance because LAI and SM have interacting effects. Hence, joint assimilation of LAI and SM from Sentinel-1 and Sentinel-2 at a 20 m resolution into the WOFOST provides a robust method to improve crop yield estimations. However, there is still bias between the key soil moisture in the root zone and the Sentinel-1 C band retrieved SM, especially when the vegetation cover is high. By active and passive microwave data fusion, it may be possible to offer a higher accuracy SM for crop yield prediction.

AB - It is well known that timely crop growth monitoring and accurate crop yield estimation at a fine scale is of vital importance for agricultural monitoring and crop management. Crop growth models have been widely used for crop growth process description and yield prediction. In particular, the accurate simulation of important state variables, such as leaf area index (LAI) and root zone soil moisture (SM), is of great importance for yield estimation. Data assimilation is a useful tool that combines a crop model and external observations (often derived from remote sensing data) to improve the simulated crop state variables and consequently model outputs like crop total biomass, water use and grain yield. In spite of its effectiveness, applying data assimilation for monitoring crop growth at the regional scale in China remains challenging, due to the lack of high spatiotemporal resolution satellite data that can match the small field sizes which are typical for agriculture in China. With the accessibility of freely available images acquired by Sentinel satellites, it becomes possible to acquire data at high spatiotemporal resolution (10-30 m, 5-6 days), which offers attractive opportunities to characterize crop growth. In this study, we assimilated remotely sensed LAI and SM into the Word Food Studies (WOFOST) model to estimate winter wheat yield using an ensemble Kalman filter (EnKF) algorithm. The LAI was calculated from Sentinel-2 using a lookup table method, and the SM was calculated from Sentinel-1 and Sentinel-2 based on a change detection approach. Through validation with field data, the inverse error was 10% and 35% for LAI and SM, respectively. The open-loop wheat yield estimation, independent assimilations of LAI and SM, and a joint assimilation of LAI + SM were tested and validated using field measurement observation in the city of Hengshui, China, during the 2016-2017 winter wheat growing season. The results indicated that the accuracy of wheat yield simulated by WOFOST was significantly improved after joint assimilation at the field scale. Compared to the open-loop estimation, the yield root mean square error (RMSE) with field observations was decreased by 69 kg/ha for the LAI assimilation, 39 kg/ha for the SM assimilation and 167 kg/ha for the joint LAI + SM assimilation. Yield coefficients of determination (R2) of 0.41, 0.65, 0.50, and 0.76 and mean relative errors (MRE) of 4.87%, 4.32%, 4.45% and 3.17% were obtained for open-loop, LAI assimilation alone, SM assimilation alone and joint LAI + SM assimilation, respectively. The results suggest that LAI was the first-choice variable for crop data assimilation over SM, and when both LAI and SM satellite data are available, the joint data assimilation has a better performance because LAI and SM have interacting effects. Hence, joint assimilation of LAI and SM from Sentinel-1 and Sentinel-2 at a 20 m resolution into the WOFOST provides a robust method to improve crop yield estimations. However, there is still bias between the key soil moisture in the root zone and the Sentinel-1 C band retrieved SM, especially when the vegetation cover is high. By active and passive microwave data fusion, it may be possible to offer a higher accuracy SM for crop yield prediction.

KW - data assimilation

KW - EnKF

KW - LAI

KW - Sentinel-1 and Sentinel-2

KW - SM

KW - winter wheat yield

KW - WOFOST

U2 - 10.3390/s19143161

DO - 10.3390/s19143161

M3 - Article

VL - 19

JO - Sensors

JF - Sensors

SN - 1424-8220

IS - 14

M1 - 3161

ER -