基于EFAST方法的AquaCrop作物模型参数全局敏感性分析

doi:10.3864/j.issn.0578-1752.2017.01.006

中国农业科学 ›› 2017, Vol. 50 ›› Issue (1): 64-76.doi: 10.3864/j.issn.0578-1752.2017.01.006

• 耕作栽培·生理生化·农业信息技术 • 上一篇下一篇

基于EFAST方法的AquaCrop作物模型参数全局敏感性分析

邢会敏^1,2,3,4，相诗尧¹，徐新刚^2,3，陈宜金¹，冯海宽^2,3，杨贵军^2,3，陈召霞^2,3

¹中国矿业大学（北京）地球科学与测绘工程学院，北京 100083；²北京农业信息技术研究中心遥感技术部，北京 100097；³国家农业信息化工程技术研究中心遥感技术部，北京 100097；⁴商丘师范学院环境与规划学院，河南商丘476000

收稿日期:2016-06-01 出版日期:2017-01-01 发布日期:2017-01-01
通讯作者: 徐新刚，Tel：010-51503215；E-mail：xxgpaper@126.com
作者简介:邢会敏，Tel：010-51503215；E-mail：hmxing1980a@163.com
基金资助:
国家自然科学基金（41571416）、北京市农林科学院创新能力建设专项（KJCX20150409）、北京市自然科学基金（4152019）

Global Sensitivity Analysis of AquaCrop Crop Model Parameters Based on EFAST Method

XING HuiMin^{1, 2, 3, 4}, XIANG ShiYao¹, XU XinGang^{2, 3}, FENG HaiKuan^{2, 3}, YANG GuiJun^{2, 3}, CHEN ZhaoXia^{2, 3}

¹College of Geoscience and Surveying Engineering, China University of Mining & Technology, Beijing 100083; ²Remote Sensing Mintech, Beijing Research Center for Information Technology in Agriculture, Beijing 100097; ³Remote Sensing Mintech, National Engineering Research Center for Information Technology in Agriculture, Beijing 100097; ⁴Department of Environment and Planning, Shangqiu Normal University, Shangqiu 476000, Henan

Received:2016-06-01 Online:2017-01-01 Published:2017-01-01

摘要/Abstract

摘要： 【目的】敏感性分析是作物模型本地化过程中的重要环节，对作物模型的校正与应用有重要的意义。【方法】本研究以国家精准农业示范研究基地2012—2013、2013—2014和2014—2015年冬小麦试验为研究对象，采用全局敏感性分析方法扩展傅里叶幅度检验法（Extended Fourier Amplitude Sensitivity Test，EFAST）对AquaCrop模型42个作物参数进行敏感性分析，以评估模型在北京地区的敏感参数。【结果】（1）对干生物量敏感作物参数是：水分和温度胁迫参数（生物量生产的最小生长度（stbio）, 引起冠层早衰的土壤水分消耗上限（psen））、生物量和产量参数（归一化水分生产力（wp））、蒸散参数（作物冠层形成后到衰老之前的作物系数（kcb））、作物冠层和物候发展参数（冠层生长系数（cgc），从播种到出苗时长（eme），最大冠层覆盖度（mcc），冠层衰老系数（cdc），从播种到成熟的时长（mat），产量形成过程中收获指数的建立长度（hilen））。其中stbio，kcb，wp和cgc 4个作物参数敏感性指数最大；（2）对冠层覆盖度最敏感的参数是：作物冠层和物候发展参数（cgc，mcc，每公顷株数（den），出苗率达到90%时的土壤覆盖度（ccs），mat和cdc）、根区发展参数（最大有效根深（rtx））、水分和温度胁迫参数（psen）、蒸散参数（kcb）；（3）对产量最敏感的参数是作物冠层和物候发展参数（从播种到开花时长（flo），mat，cdc，hilen和从播种到开始衰老时长（sen））、水分和温度胁迫参数（psen）、生物量和产量参数（参考收获指数（hi）和wp）、蒸散参数（kcb）。【结论】利用EFAST方法对AquaCrop模型中的作物参数进行一阶和全局敏感分析，最大干物量的敏感性分析结果以及干生物量随时间变化的敏感性分析结果显示，敏感性参数的选择上差异不大，但排序上存在较大的差异，最大干生物量的敏感性分析不能分析作物参数对干生物量在整个生育期的影响，结果不全面；冠层覆盖度随时间变化的一阶和全局敏感性分析结果显示，在敏感参数的选择和排序上均有较好的一致性，全局敏感性分析中作物参数的敏感性指数更高，对冠层覆盖度的影响表现得更明显。本研究结果用于AquaCrop模型本地化，可提高该模型在北京地区的模拟效率和模拟精度。

关键词: 冬小麦, AquaCrop模型, 敏感性分析, EFAST方法, 干生物量

Abstract: 【Objective】Sensitivity analysis is an important link in crop model localization, and it plays an important role in AquaCrop model calibration and application.【Method】In this study, in order to identify the sensitivity parameters, the 2012-2013, 2013-2014 and 2014-2015 winter wheat experiments were conducted in National Precision Agriculture Demonstration Research Base in Beijing, China, the Extended Fourier Amplitude Sensitivity Test (EFAST) method was used to carry out sensitivity analysis of 42 crop parameters of AquaCrop model.【Result】The sensitivity parameters were: (1) For dry biomass: water and temperature stress (minimum growing degrees required for full biomass production (stbio), upper threshold of soil water depletion factor for canopy senescence (psen)), biomass and yield production (water productivity normalized (wp)), transpiration (crop coefficient when canopy is complete but prior to senescence (kcb)), canopy and phaenological development (GGD-increase in canopy cover (cgc), GDD-from sowing to emergence (eme), maximum canopy cover in fraction soil cover (mcc), GGD-decrease in canopy cover (cdc), total length of crop cycle in growing degree-days (mat), building-up of harvest index during yield formation (hilen)). stbio, kcb, wp and cgc were the four most sensitive parameters; (2) For canopy cover: canopy and phaenological development (cgc, mcc, number of plants per hectare (den), soil surface covered by an individual seedling at 90% emergence (ccs), mat and cdc), root development (maximum effective rooting depth (rtx)), water and temperature stress (psen), transpiration (kcb); (3) For yield: canopy and phaenological development (GDD-from sowing to flowering (flo), mat, cdc, hilen and GDD-from sowing to start senescence (sen)), water and temperature stress (psen), biomass and yield production (reference harvest index (hi) and wp), transpiration (kcb).【Conclusion】The results of first order and total order sensitivity analysis for AquaCrop model of winter wheat maximum dry biomass and dry biomass time-varying showed that there was a little difference in the choice of sensitivity parameters, but many differences in the ranking. The sensitivity analysis of maximum dry biomass was not comprehensive, which could not analyze the effect of crop parameters on dry biomass during the whole growth period. The results of the first order and total order sensitivity analysis for AquaCrop model of winter wheat canopy cover time-varying showed that there was a good consistency in the selection and ranking of sensitive parameters. The values of total order sensitivity indices of crop parameters were higher than first order, and the influences on canopy cover were more obvious. This study provides guidelines for AquaCrop model calibration and application in Beijing, China, as well providing guidance to simplify the AquaCrop model and improve its precision, especially when many parameters are used.

Key words: winter wheat, AquaCrop model, sensitivity analysis, EFAST method, biomass

邢会敏，相诗尧，徐新刚，陈宜金，冯海宽，杨贵军，陈召霞. 基于EFAST方法的AquaCrop作物模型参数全局敏感性分析[J]. 中国农业科学, 2017, 50(1): 64-76.

XING HuiMin, XIANG ShiYao, XU XinGang, FENG HaiKuan, YANG GuiJun, CHEN ZhaoXia. Global Sensitivity Analysis of AquaCrop Crop Model Parameters Based on EFAST Method[J]. Scientia Agricultura Sinica, 2017, 50(1): 64-76.

0
/ / 推荐

导出引用管理器 EndNote|Reference Manager|ProCite|BibTeX|RefWorks

链接本文: https://www.chinaagrisci.com/CN/10.3864/j.issn.0578-1752.2017.01.006

https://www.chinaagrisci.com/CN/Y2017/V50/I1/64

参考文献

[1] ITTERSUM M K V, DONATELLI M. Modelling cropping systems—highlights of the symposium and preface to the special issues. European Journal of Agronomy, 2003, 18(3): 187-197.

[2] BUTTERBACH-BAHL K, KESIK M, MIEHLE P, PAPEN H, LI C. Quantifying the regional source strength of N-trace gases across agricultural and forest ecosystems with process based models. Plant & Soil, 2004, 260(1/2): 311-329.

[3] KLEPPER O. Multivariate aspects of model uncertainty analysis: Tools for sensitivity analysis and calibration. Ecological Modelling, 1997, 101(1): 1-13.

[4] CONFALONIERI R, BELLOCCHI G, TARANTOLA S, ACUTIS M, DONATELLI M, GENOVESE G. Sensitivity analysis of the rice model WARM in Europe: Exploring the effects of different locations, climates and methods of analysis on model sensitivity to crop parameters. Environmental Modelling & Software, 2010, 25(4): 479-488.

[5] SALTELLI A, RATION M, ANDRES T, CAMPOLONGO F, CARIBONI J, GATELLI D, SAISANA M, TARANTOLA S. Global Sensitivity Analysis: The Primer. Chichester: John Wiley & Sons Inc, 2008.

[6] GRIENSVEN A V, MEIXNER T, GRUNWALD S, BISHOP T, DILUZIO M, SRINIVASAN R. A global sensitivity analysis tool for the parameters of multi-variable catchment models. Journal of Hydrology, 2006, 324(1/4): 10-23.

[7] SALTELLI A, ANNONI P. How to avoid a perfunctory sensitivity analysis. Environmental Modelling & Software, 2010, 25(12): 1508-1517.

[8] SALTELLI A, TARANTOLA S, CAMPOLONGO F. Sensitivity analysis as an ingredient of modeling. Statistical Science, 2000, 15(4): 377-395.

[9] MORRIS M D. Factorial sampling plans for preliminary computational experiments. Technometrics, 1991, 33(2): 161-174.

[10] SOBOL’ I M. Sensitivity estimates for nonlinear mathematical models. Mathematical Modelling and Computational Experiments, 1990, 2(1): 112-118.

[11] CUKIER R I, LEVINE H B, SHULER K E. Nonlinear sensitivity analysis of multiparameter model systems. Journal of Computational Physics, 1978, 26(1): 1-42.

[12] SALTELLI A. Sensitivity analysis: Could better methods be used? Journal of Geophysical Research Atmospheres, 1999, 104(D3): 3789-3793.

[13] SALTELLI A, TARANTOLA S, CHAN P S. A quantitative model-independent method for global sensitivity analysis of model output. Technometrics, 1999, 41(1): 39-56.

[14] KIOUTSIOUKIS I, TARANTOLA S, SALTELLI A, GATELLI D. Uncertainty and global sensitivity analysis of road transport emission estimates. Atmospheric Environment, 2004, 38(38): 6609-6620.

[15] SALTELLI A, TARANTOLA S, CAMPOLONGO F, RATTO M. Sensitivity analysis in practice: A guide to assessing scientific models. Journal of the American Statistical Association, 2006, 101(473): 398-399.

[16] 姜志伟, 陈仲新, 周清波, 任建强. CERES-Wheat作物模型参数全局敏感性分析. 农业工程学报, 2011, 27(1): 236-242.

JIANG Z W, CHEN Z X, ZHOU Q B, REN J Q. Global sensitivity analysis of CERES-Wheat model parameters. Transactions of the Chinese Society of Agricultural Engineering, 2011, 27(1): 236-242. (in Chinese)

[17] DEJONGE K C, ASCOUGHⅡ J C, AHMADI M, ANDALES A A, ARABI M. Global sensitivity and uncertainty analysis of a dynamic agroecosystem model under different irrigation treatments. Ecological Modelling, 2012, 231(4): 113-125.

[18] 宋明丹, 冯浩, 李正鹏, 高建恩. 基于Morris和EFAST的CERES-Wheat模型敏感性分析. 农业机械学报, 2014, 45(10): 124-131.

SONG M D, FENG H, LI Z P, GAO J E. Global sensitivity analyses of DSSAT-CERES-wheat model using Morris and EFAST methods. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(10): 124-131. (in Chinese)

[19] 吴立峰, 张富仓, 范军亮, 周罕觅, 邢英英, 强生才. 不同灌水水平下CROPGRO棉花模型敏感性和不确定性分析. 农业工程学报, 2015(15): 55-64.

WU L F, ZHANG F C, FAN J L, ZHOU H M, XING Y Y, QIANG S C. Sensitivity and uncertainty analysis for CROPGRO-cotton model at different irrigation levels. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(15): 55-64. (in Chinese)

[20] 何亮, 赵刚, 靳宁, 庄伟, 于强. 不同气候区和不同产量水平下APSIM-Wheat模型的参数全局敏感性分析. 农业工程学报, 2015, 31(14): 148-157.

HE L, ZHAO G, JIN N, ZHUANG W, YU Q. Global sensitivity analysis of APSIM-Wheat parameters in different climate zones and yield levels. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(14): 148-157. (in Chinese)

[21] WANG J, LI X, LU L, Feng F. Parameter sensitivity analysis of crop growth models based on the extended Fourier Amplitude Sensitivity Test method. Environmental Modelling & Software, 2013, 48(5): 171-182.

[22] 何亮, 侯英雨, 赵刚, 邬定荣, 于强. 基于全局敏感性分析和贝叶斯方法的WOFOST作物模型参数优化. 农业工程学报, 2016, 32(2): 169-179.

HE L, HOU Y Y, ZHAO G, WU D R, YU Q. Parameters optimization of WOFOST model by integration of global sensitivity analysis and Bayesian calibration method. Transactions of the Chinese Society of Agricultural Engineering, 2016, 32(2): 169-179. (in Chinese)

[23] 吴锦, 余福水, 陈仲新, 陈晋. 基于EPIC模型的冬小麦生长模拟参数全局敏感性分析. 农业工程学报, 2009, 25(7): 136-142.

WU J, YU F S, CHEN Z X, CHEN J. Global sensitivity analysis of growth simulation parameters of winter wheat based on EPIC model. Transactions of the Chinese Society of Agricultural Engineering, 2009, 25(7): 136-142. (in Chinese)

[24] VANUYTRECHT E, RAES D, WILLEMS P. Global sensitivity analysis of yield output from the water productivity model. Environmental Modelling & Software, 2014, 51(1): 323-332.

[25] 金秀良. 基于AqoaCrop模型与多源遥感数据的北方冬小麦水分利用效率估算[D]. 扬州: 扬州大学, 2015.

JIN X L. Estimation of water use efficiency of winter wheat based on AquaCrop model and multi-source remote sensing data in Northern[D]. Yangzhou: Yangzhou University, 2015. (in Chinese)

[26] 黄健熙, 武思杰, 刘兴权, 马冠南, 马鸿元, 吴文斌, 邹金秋. 基于遥感信息与作物模型集合卡尔曼滤波同化的区域冬小麦产量预测. 农业工程学报, 2012, 28(4): 142-148.

HUANG J X, WU S J, LIU X Q, MA G N, MA H Y, WU W B, ZOU J Q. Regional winter wheat yield forecasting based on assimilation of remote sensing data and crop growth model with Ensemble Kalman method. Transactions of the Chinese Society of Agricultural Engineering, 2012, 28(4): 142-148. (in Chinese)

[27] 黄健熙, 马鸿元, 田丽燕, 王鹏新, 刘峻明. 基于时间序列LAI 和ET 同化的冬小麦遥感估产方法比较. 农业工程学报, 2015, 31(4): 197-203.

HUANG J X, MA H Y, TIAN L Y, WANG P X, LIU J M. Comparison of remote sensing yield estimation methods for winter wheat based on assimilating time-sequence LAI and ET. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(4): 197-203. (in Chinese)

[28] 姜志伟, 陈仲新, 任建强, 周清波. 粒子滤波同化方法在CERES- Wheat作物模型估产中的应用. 农业工程学报, 2012, 28(14): 138-146.

JIANG Z W, CHEN Z X, REN J Q, ZHOU Q B. Estimation of crop yield using CERES-Wheat model based on particle filter data assimilation method. Transactions of the Chinese Society of Agricultural Engineering, 2012, 28(14): 138-146. (in Chinese)

[29] 解毅, 王鹏新, 刘峻明, 李俐. 基于四维和变分和集合卡尔曼滤波同化方法的冬小麦单产估测. 农业工程学报, 2015, 31(1): 187-195.

XIE Y, WANG P X, LIU J M, LI L. Winter wheat yield estimation based on assimilation method combined with 4DVAR and EnKF. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(1): 187-195. (in Chinese)

[30] STEDUTO P, HSIAO T C, RAES D, FERERES E. AquaCrop—The FAO crop model to simulate yield response to water: I. Concepts and underlying principles. Agronomy Journal, 2009, 101(3): 426-437.

[31] RAES D, STEDUTO P, HSIAO T C, FERERES E. AquaCrop —The FAO crop model to simulate yield response to water: II. Main algorithms and software description. Agronomy Journal, 2009, 101(3): 438-447.

[32] HSIAO T C, HENG L, STEDUTO P, ROJAS L B, RASE D, FERERES E. AquaCrop—The FAO crop model to simulate yield response to water: III. Parameterization and testing for maize. Agronomy Journal, 2009, 101(3):448-459.

[33] DOORENBOS J, KASSAM A H. Yield response to water. Irrigation And Drainage Paper no. 33. FAO, Rome, 1979.

[34] FAO. AquaCrop –the FAO crop model to simulated yield response to water: Reference manual annex I – AquaCrop crop parameters. www. Fao.org/nr/water/Aquacrop.html, 2012.

[35] 李振海. 基于遥感数据和气象预报数据的DSSAT模型冬小麦产量和品质预报[D]. 杭州: 浙江大学, 2016.

LI Z H. Predicting winter wheat yield and quality by integrating of remote-sensing data and the weather forecast data into the DSSAT model[D]. Hangzhou: Zhejiang University, 2016. (in Chinese)

基于EFAST方法的AquaCrop作物模型参数全局敏感性分析

Global Sensitivity Analysis of AquaCrop Crop Model Parameters Based on EFAST Method

RichHTML

PDF

赞

可视化

被引次数

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

Metrics

本文评价

推荐阅读 0

[1]	徐久凯, 袁亮, 温延臣, 张水勤, 李燕婷, 李海燕, 赵秉强. 畜禽有机肥氮在冬小麦季对化肥氮的相对替代当量[J]. 中国农业科学, 2023, 56(2): 300-313.
[2]	王洋洋,刘万代,贺利,任德超,段剑钊,胡新,郭天财,王永华,冯伟. 基于多元统计分析的小麦低温冻害评价及水分效应差异研究[J]. 中国农业科学, 2022, 55(7): 1301-1318.
[3]	伊英杰,韩坤,赵斌,刘国利,林佃旭,陈国强,任昊,张吉旺,任佰朝,刘鹏. 长期不同施肥措施冬小麦-夏玉米轮作体系周年氨挥发损失的差异[J]. 中国农业科学, 2022, 55(23): 4600-4613.
[4]	刘丰,蒋佳丽,周琴,蔡剑,王笑,黄梅,仲迎鑫,戴廷波,曹卫星,姜东. 美国软麦籽粒品质变化趋势及对我国弱筋小麦标准达标度分析[J]. 中国农业科学, 2022, 55(19): 3723-3737.
[5]	韩守威,司纪升,余维宝,孔令安,张宾,王法宏,张海林,赵鑫,李华伟,孟鈺. 山东省冬小麦产量差与氮肥利用效率差形成机理解析[J]. 中国农业科学, 2022, 55(16): 3110-3122.
[6]	孟雨,温鹏飞,丁志强,田文仲,张学品,贺利,段剑钊,刘万代,冯伟. 基于热红外图像的小麦品种抗旱性鉴定与评价[J]. 中国农业科学, 2022, 55(13): 2538-2551.
[7]	高志源,许吉利,刘硕,田汇,王朝辉. 大田群体冬小麦氮收获指数变异特征研究[J]. 中国农业科学, 2021, 54(3): 583-595.
[8]	毛安然,赵护兵,杨慧敏,王涛,陈秀文,梁文娟. 不同覆盖时期和覆盖方式对旱地冬小麦经济和环境效应的影响[J]. 中国农业科学, 2021, 54(3): 608-618.
[9]	向晓玲,陈松鹤,杨洪坤,杨永恒,樊高琼. 秸秆覆盖与施磷对丘陵旱地小麦产量和磷素吸收利用效应的影响[J]. 中国农业科学, 2021, 54(24): 5194-5205.
[10]	高兴祥,张悦丽,安传信,李美,李健,房锋,张双应. 山东省冬小麦田杂草群落调查及其变化原因分析[J]. 中国农业科学, 2021, 54(24): 5230-5239.
[11]	宗毓铮,张函青,李萍,张东升,林文,薛建福,高志强,郝兴宇. 大气CO₂与温度升高对北方冬小麦旗叶光合特性、碳氮代谢及产量的影响[J]. 中国农业科学, 2021, 54(23): 4984-4995.
[12]	王金凤,王壮壮,谷丰序,牟海萌,王宇,段剑钊,冯伟,王永华,郭天财. 氮密调控对两个冬小麦品种碳氮代谢及产量的影响[J]. 中国农业科学, 2021, 54(19): 4070-4083.
[13]	李敏, 苏慧, 李阳阳, 李金鹏, 李金才, 朱玉磊, 宋有洪. 黄淮海麦区小麦耐热性分析及其鉴定指标的筛选[J]. 中国农业科学, 2021, 54(16): 3381-3392.
[14]	费帅鹏,禹小龙,兰铭,李雷,夏先春,何中虎,肖永贵. 基于高光谱遥感和集成学习方法的冬小麦产量估测研究[J]. 中国农业科学, 2021, 54(16): 3417-3427.
[15]	马红梅,曹寒冰,谢英荷,李廷亮,刘凯,张奇茹,姜丽伟,曹静,邵靖琳,武文玥,栗文琪. 晋南黄土旱塬小麦养分投入与化肥减施经济环境效应评价[J]. 中国农业科学, 2021, 54(13): 2804-2817.