Study on Applicability of WOFOST Model in Xuchang Tobacco Production

doi:10.11923/j.issn.2095-4050.cjas2025-0081

Abstract

Abstract:

To address the issue of relying on experience and lacking quantitative simulation tools in the production management of flue-cured tobacco in central Henan, in order to clarify the applicability of WOFOST model in Xuchang flue-cured tobacco production area, the parameter sensitivity analysis was carried out by OAT (one-at-a-time) method, and the localization calibration and independent verification of the model were completed by ' trial and error method '. The study was based on the field observation and meteorological data of Jian'an District and Xiangcheng County of Xuchang, Henan Province from 2021 to 2022, and the model's simulation accuracy for leaf dry weight (WLV), above-ground biomass (AGP), and leaf area index (LAI) was evaluated using the coefficient of determination (R²), consistency index (d), and normalized root mean square error (NRMSE). The main conclusions are as follows. (1) In the calibration of leaf dry weight (WLV), the simulation values from the WOFOST model showed a significant linear relationship with the observed values, indicating high accuracy in simulating leaf dry weight. The calibration accuracy for above-ground biomass (AGP) was also good. However, the calibration of the leaf area index (LAI) was less accurate, with NRMSE ranging from 20% to 30%. (2) The model was validated using data from Xuchang in 2021, Xuchang and Xiangcheng in 2022, and the simulation results were good for the most important biomass, leaf dry weight. Overall, the WOFOST model showed high accuracy in validation. The R2 and consistency index d of 90 % (8 / 9) validation data were higher than 0.8, nearly 80 % (7 / 9) were higher than 0.9, and NRMSE was between 10 % and 20 %. The correlation and consistency between simulated and measured values were good. (3) After calibration and verification, the WOFOST model could effectively simulate the growth process of tobacco in Xuchang, Henan, providing a foundation for quantitative and digital management of tobacco production. Subsequently, remote sensing data can be coupled to carry out research on crop growth monitoring and yield forecasting.

Key words: tobacco, WOFOST model, quantification, parameter localization, calibration, verification

CLC Number:

LI Wenfeng, JIN Shuyuan, PU Tuanwei. Study on Applicability of WOFOST Model in Xuchang Tobacco Production[J]. Journal of Agriculture, 2026, 16(5): 101-108.

Figures/Tables 5

References 23

[1]	郭建茂, 王琦, 施俊怡, 等. 遥感信息与作物模型结合在冬小麦区域模拟中的应用[J]. 大气科学学报, 2014, 37(2):237-242.
[2]	郑昌玲, 张蕾, 侯英雨, 等. 基于WOFOST模型的冬小麦产量动态预报方法[J]. 干旱地区农业研究, 2022, 40(6):242-250,267.
[3]	郑昌玲, 侯英雨, 吴门新, 等. 基于WOFOST模型的华北冬小麦动态长势评估指标构建[J]. 麦类作物学报, 2020, 40(6):746-753.
[4]	朱津辉, 郭建茂, 毛留喜. 基于模型的河北省保定市冬小麦最佳灌溉方案研究[J]. 气象, 2014, 40(11):1398-1407.
[5]	张阳, 王连喜, 李琪, 等. 基于WOFOST模型的吉林省中西部春玉米灌溉模拟[J]. 中国农业气象, 2018, 39(6):411-420.
[6]	张微玮. 基于WOFOST模型的春玉米生长模型精细化及其应用[D]. 沈阳: 沈阳农业大学, 2018.
[7]	王秀珍. 玉溪烤烟气象灾害风险及作物模型适用性研究[D]. 南京: 南京信息工程大学, 2021.
[8]	杨霏云, 郑秋红, 李文科, 等. 基于WOFOST模型的辽宁省春玉米干旱灾损风险评估[J]. 干旱地区农业研究, 2020, 38(6):218-225.
[9]	郑昌玲, 魏瑞江, 张蕾. 基于WOFOST模型的冬小麦长势评估和产量预报——以2021年度冬小麦为例[J]. 气象科学, 2025, 45(2):284-292.
[10]	赵晓松, 关德新, 吴家兵, 等. 长白山阔叶红松林的零平面位移和粗糙度[J]. 生态学杂志, 2004(5):84-88.
[11]	郭彩云. 基于WOFOST模型的烟叶灌溉方案及生长预测与评估研究[D]. 南京: 南京信息工程大学, 2022.
[12]	ANGSTROM A. Solar and terrestrial radiation. Report to the international commission for solar AquaCrop[J]. Agronomy journal, 2009, 101(3):477-487. doi: 10.2134/agronj2008.0179s URL
[13]	和清华, 谢云. 我国太阳总辐射气候学计算方法研究[J]. 自然资源学报, 2010, 25(2):308-319. doi: 10.11849/zrzyxb.2010.02.015
[14]	何亮, 侯英雨, 赵刚, 等. 基于全局敏感性分析和贝叶斯方法的WOFOST作物模型参数优化[J]. 农业工程学报, 2016, 32(2):169-179
[15]	何亮, 赵刚, 靳宁, 等. 不同气候区和不同产量水平下APSIM-Wheat模型的参数全局敏感性分析[J]. 农业工程学报, 2015, 31(14):148-157.
[16]	DANIEL C. One- at- a-time plans[J]. Journal of the American statisctical association, 1973, 68(342):353-360.
[17]	邬定荣, 欧阳竹, 赵小敏, 等. 作物生长模型WOFOST在华北平原的适用性研究[J]. 植物生态学报, 2003, 27(5):594-602.
[18]	马玉平, 王石立, 张黎. 针对华北小麦越冬的WOFOST模型改进[J]. 中国农业气象, 2005, 26(3):145-149.
[19]	MORRIS M D. Factorial sampling plants for preliminary computational experiments[J]. Technometrics, 1991, 33(2):161-173. doi: 10.1080/00401706.1991.10484804 URL
[20]	MAKOWSKI D, WALLACH D, TREMBLAY M. Using a Bayesian approach to parameter estimation; comparison of the GLUE and MCMC methods[J]. Agronomie, 2002, 22(2):191-203. doi: 10.1051/agro:2002007 URL
[21]	胡雪琼, 徐梦莹, 买苗, 等. WOFOST模型对于云南烤烟的适用性研究[J]. 南京信息工程大学学报(自然科学版), 2015, 7(5):451-457.
[22]	RADER M, KIRSHEN P, RONCOLI C. Agricultural risk decision support system for resource poor farmers in Burkina Faso, West Africa[J] Journal of water resources planning and management, 2009, 135(5):323-333. doi: 10.1061/(ASCE)0733-9496(2009)135:5(323) URL
[23]	姜志伟, 陈仲新, 周清波, 等. CERES-Wheat作物模型参数全局敏感性分析[J]. 农业工程学报, 2011, 27(1):236-242.

参数名称	参数含义	叶重(WLV)	茎重(WST)	地上部总干物重(TAGP)	叶面积指数(LAI)	数值	单位
SLATB0.0	以发育阶段（DVS）为函数的比叶面积	0.56093	0.55689	0.53642	0.7806	0.0017	hm²/kg
SLATB0.5	0.10876	0.03288	0.02777	0.19794	0.0012	hm²/kg hm²/kg
SLATB1.0	1.21008	0.12013	0.08157	1.31822	0.0008	hm²/kg
SLATB2.0	0.21223	0.01581	0.01053	0.26948	0.0003	hm²/kg
TSUM1	出苗到开花所需积温	0.53479	0.11222	0.11929	0.51937	2200	℃·d
TSUM2	开花到成熟所需积温	1.63854	0.4614	0.08856	2.32127	1800	℃·d
RMR	根维持呼吸消耗系数	0.5832	0.09776	0.07437	0.36482	0.01	-
RML	叶片维持呼吸消耗系数	0.00179	6.23E-04	4.23E-04	0.001	0.01	-
Q10	温度变化10℃呼吸作用变化的速率	0.82604	0.0981	0.07354	0.79228	2	-
AMAXTB0.0	以发育阶段(DVS)为函数的 CO₂最大同化速率	0.39479	0.40579	0.37682	0.33764	5.0	kg/（hm²·h）
AMAXTB0.5		0.53479	0.11222	0.11929	0.51937	8.0	kg/（hm²·h）
AMAXTB1.0		0.38859	0.18289	0.12662	0.37778	60.0	kg/（hm²·h）
AMAXTB2.0		1.06445	0.09435	0.06271	0.97152	8.0	kg/（hm²·h）
CVL	叶片的干物质转化效率	0.47957	0.29203	0.23432	0.45277	0.72	kg/kg
CVR	根的干物质转化效率	0.25528	0.16158	0.12921	0.25143	0.72	kg/kg
CVS	茎的干物质转化效率	0.09663	0.06294	0.06113	0.09393	0.45	kg/kg
RMS	茎维持呼吸消耗系数	1.41066	0.20428	0.1385	1.36864	0.045	-
SPAN	叶片衰老系数	2.22831	0.15882	0.11179	2.19934	50.0	kg/ hm²

参数名称	参数含义	叶重(WLV)	茎重(WST)	地上部总干物重(TAGP)	叶面积指数(LAI)	数值	单位
SLATB0.0	以发育阶段（DVS）为函数的比叶面积	0.56093	0.55689	0.53642	0.7806	0.0017	hm²/kg
SLATB0.5	0.10876	0.03288	0.02777	0.19794	0.0012	hm²/kg hm²/kg
SLATB1.0	1.21008	0.12013	0.08157	1.31822	0.0008	hm²/kg
SLATB2.0	0.21223	0.01581	0.01053	0.26948	0.0003	hm²/kg
TSUM1	出苗到开花所需积温	0.53479	0.11222	0.11929	0.51937	2200	℃·d
TSUM2	开花到成熟所需积温	1.63854	0.4614	0.08856	2.32127	1800	℃·d
RMR	根维持呼吸消耗系数	0.5832	0.09776	0.07437	0.36482	0.01	-
RML	叶片维持呼吸消耗系数	0.00179	6.23E-04	4.23E-04	0.001	0.01	-
Q10	温度变化10℃呼吸作用变化的速率	0.82604	0.0981	0.07354	0.79228	2	-
AMAXTB0.0	以发育阶段(DVS)为函数的 CO₂最大同化速率	0.39479	0.40579	0.37682	0.33764	5.0	kg/（hm²·h）
AMAXTB0.5		0.53479	0.11222	0.11929	0.51937	8.0	kg/（hm²·h）
AMAXTB1.0		0.38859	0.18289	0.12662	0.37778	60.0	kg/（hm²·h）
AMAXTB2.0		1.06445	0.09435	0.06271	0.97152	8.0	kg/（hm²·h）
CVL	叶片的干物质转化效率	0.47957	0.29203	0.23432	0.45277	0.72	kg/kg
CVR	根的干物质转化效率	0.25528	0.16158	0.12921	0.25143	0.72	kg/kg
CVS	茎的干物质转化效率	0.09663	0.06294	0.06113	0.09393	0.45	kg/kg
RMS	茎维持呼吸消耗系数	1.41066	0.20428	0.1385	1.36864	0.045	-
SPAN	叶片衰老系数	2.22831	0.15882	0.11179	2.19934	50.0	kg/ hm²

统计指标	2021年襄城叶干重	2021年襄城地上部生物量	2021年襄城叶面积指数
RMSE	363.37	606.59	0.58
NRMSE/%	16.21	15.06	22.20
R²(1:1)	0.89	0.93	0.95
d	0.86	0.92	0.81

统计指标	2021年襄城叶干重	2021年襄城地上部生物量	2021年襄城叶面积指数
RMSE	363.37	606.59	0.58
NRMSE/%	16.21	15.06	22.20
R²(1:1)	0.89	0.93	0.95
d	0.86	0.92	0.81

模拟对象		RMSE	NRMSE/%	R²(1:1)	d
叶干重	2021年许昌	254.67	10.23	0.93	0.92
	2022年许昌	276.02	11.92	0.95	0.94
	2022年襄城	271.46	12.75	0.93	0.92
地上部生物量	2021年许昌	488.34	11.14	0.96	0.95
	2022年许昌	349.31	9.60	0.98	0.98
	2022年襄城	905.81	24.63	0.83	0.80
叶面积指数	2021年许昌	0.34	12.29	0.86	0.83
	2022年许昌	0.36	15.17	0.95	0.90
	2022年襄城	0.35	16.98	0.90	0.90