高级检索+

基于数据同化及参数优化的WOFOST模型性能改进

Improvements in the performance of the WOFOST model based on data assimilation and parameter optimization

  • 摘要: 针对作物模型对干旱过程模拟能力不足问题,基于在锦州农业气象试验站开展的多年春玉米分期播种试验观测数据,利用集合卡尔曼滤波(EnKF)同化方法,采用叶面积指数(leaf area index, LAI)数据同化和土壤湿度(soil moisture, SM)数据同化、光合参数优化和水分胁迫函数(reduction factor for water stress, RFWS)改进开展多种组合方案模拟。各方案模拟结果表明:在水分正常年份(2011—2013年),原模型的地上生物量和产量模拟平均绝对相对误差(MARE)分别为42.4%和38.3%,而LAI数据同化与光合参数优化的组合方案最优,其MARE分别降至14.2%和11.5%;在干旱年份(2014—2015年),原模型的地上生物量和产量MARE分别为32.7%和54.1%,SM数据同化、光合参数优化与RFWS改进的组合方案最优,其MARE分别为18.0%和34.1%。利用水分正常年的最优方案对2011—2015及2023年进行滚动模拟,确定拔节期后约60d为产量预测趋稳的“最佳节点”,并确定该阶段累计降水量149 mm为干旱阈值。基于此提出差异化模拟策略,利用2015、2018、2019和2024年开展滚动预报验证,表明该方案在不同气象条件下均取得较好效果。该研究为发展作物模型气象年型自适应动态模拟新框架提供参考。

     

    Abstract: This study aimed to enhance the performance of the WOFOST crop model in simulating spring maize growth and yield under varying water conditions, with particular emphasis on improving the representation of drought-induced water stress. The research focused on addressing the model's limitations in capturing soil moisture effects and developing a dynamic, adaptive framework for yield forecasting based on data assimilation and parameter optimization. Multi-year, multi-sowing-date field experiments were conducted at the Jinzhou Agrometeorological Experimental Station in Northeast China, covering three distinct moisture scenarios: normal water years (2011-2013), typical drought years (2014, 2015, and 2018), and mild drought years (2019, 2023, and 2024). This comprehensive dataset provides a robust foundation for systematically evaluating the model's performance across different water regimes. An Ensemble Kalman Filter (EnKF) data assimilation approach was implemented using Leaf Area Index (LAI) and Soil Moisture (SM) observations obtained from field experiments. Key photosynthetic parameters, including the maximum leaf CO2 assimilation rate (AMAXTB), were optimized through iterative calibration against observed biomass and yield data. Additionally, the soil water stress response function (RFWS) was substantially modified to better reflect the nonlinear effects of soil moisture deficits on photosynthetic processes and biomass accumulation. The modification addressed the oversensitivity of the original function to slight moisture reductions below critical thresholds, which had previously led to excessive yield underestimation under drought conditions. Multiple simulation schemes combining different assimilation and optimization strategies were systematically designed and evaluated. These included: S1 (original model), S2 (LAI assimilation only), S3 (LAI assimilation with parameter optimization), S4 (SM assimilation only), S5 (SM assimilation with parameter optimization), and S9 (SM assimilation with both parameter optimization and RFWS modification). Each scheme was rigorously evaluated by comparing simulated above-ground biomass and final yield with comprehensive field observations under both normal and drought conditions. Model performance assessment revealed substantial improvements through the integrated approach. In normal-precipitation years (2011-2013), LAI assimilation combined with parameter optimization (Scheme S3) achieved the highest accuracy, reducing mean absolute relative errors (MARE) of yield and above-ground biomass from 38.3% and 42.4% in the original scheme to 11.5% and 14.2%, respectively. This improvement demonstrates the particular effectiveness of LAI assimilation under adequate moisture conditions where leaf development primarily drives yield formation. In drought years (2014-2015), SM assimilation combined with parameter optimization and modified RFWS function (Scheme S9) performed best, lowering yield and biomass MARE from 54.1% and 32.7% to 34.1% and 18.0%, respectively. The 2014 case study was particularly insightful, where despite severe meteorological drought, actual yield exceeded expectations due to carry-over soil moisture from the previous wet year (2013). Scheme S9 captured this anomaly more effectively than other schemes, though some yield underestimation persisted, indicating remaining challenges in extreme drought simulation. Rolling forecasts implemented throughout the growing season indicated that approximately 60 days after the jointing stage represented the optimal prediction node, corresponding to the critical period of water demand for maize and the phase when yield components become largely determined. A differentiated assimilation strategy was developed based on a cumulative precipitation threshold of 149 mm during this period: Scheme S3 (LAI assimilation with photosynthetic parameter optimization) was applied when accumulated precipitation exceeded 149 mm, while Scheme S9 (SM assimilation with photosynthetic parameter optimization and modified water stress function) was implemented when accumulated precipitation fell below or equal to this threshold. Validation with 2015, 2018, 2019, and 2024 sowing dates showed that the differences between simulated and observed yields were -649 kg/hm2 in 2018 and 58 kg/hm2 in 2019, demonstrating the model's capability in capturing yield trends across varying water conditions. The combined framework of data assimilation, parameter optimization, and water stress function modification substantially enhanced WOFOST’s ability to simulate spring maize growth and yield across diverse water conditions. This framework provides a practical tool for yield prediction, drought impact assessment, and decision support under varying water conditions. Future work could focus on refining drought response mechanisms, incorporating multi-factor parameterization including temperature and radiation effects, and extending the framework to other crops and ecological regions to improve model stability and applicability.

     

/

返回文章
返回