论文部分内容阅读
关键词： 枣;DNDC模型;干旱胁迫;敏感性分析;参数调整
中图分类号： S165 .27 文献标识码： A 文章编号： 1000-4440（2021）03-0812-05
Adjustment of jujube yield parameters based on DNDC model under drought stress
WANG De-juan， SUI Li-chun， LIU Di， WANG Jin， YUAN Huan-huan， XU Jia-li
（School of Geological Engineering and Geomatics， Chang’an University， Xi’an 710054，China）
Key words： jujube;DNDC model;drought stress;sensitivity analysis;parameter adjustment
红枣作为南疆林果业发展的优势树种，在新疆绿洲生产、生态建设中起着重要作用[1]。新疆地区应用的估产方法主要有气象估计，经验估计等[2-3]。段丁丁等[4]利用DSSAT内部的马铃薯模型对马铃薯进行了区域化模拟估产研究;吴立峰[5]等通过Morris方法和EFAST方法对不同灌溉水平下CROPGRO棉花模型进行了敏感性分析和不确定性分析;纪甫江[6]等通过将CASA模型和WOFOST模型进行耦合从而提高了大豆单产估算的精度。而目前对于枣的产量估算研究较少，因此结合新疆地区的生态条件，利用更加精准的作物生长模型将成为在枣类估产中的有效方法。
经过多年的发展和应用，作物模型主要分为三大派系，分别是以ORYZA2000模型为代表的荷兰作物模型体系，以CERES模型为代表的美国作物模型体系，以CCSODS模型为代表的中国作物模型体系[7-9]。但在众多作物生长模型中用于果树产量模拟的模型较少，其中DNDC模型是由美国新罕布什尔大学的李长生教授团队研发的生物地球化学模型，该模型是农业生态系统中一系列控制碳和氮迁移转化的生物化学及地球化学反应机制的计算机模拟表达[10]。DNDC模型以日为步长，通过气候、土壤、植被和管理4个基本生态因子对作物生长进行机理性模拟，并通过DNDC模型特色的硝化反硝化作用模拟作物在生长过程中的碳、氮元素的变化，在进行产量预报的同时进一步实现碳、氮元素在地球生态系统中循环的模拟[11-13]。DNDC模型包含简单易操作的新作物类构建方法，并且对水分处理具有单独的页面，因此对枣树的产量模拟可操作性更高。
作物生长模型作为机理性模型需要大量的作物和实际的环境参数作为输入参数，但只有部分参数对模型的模拟结果起决定作用，因此参数的准确性将会直接影响模型模拟结果的真实性[14-15]。同时试验地区的特殊环境因素也成为参数敏感性分析的关键。新疆地区特殊的地理位置及天气条件，使得水分成为作物生长过程中不容忽略的问题。参数敏感性分析主要分为全局参数敏感性分析和局部参数敏感性分析[16]。局部敏感性分析是指在一定范围内调整某一个参数的取值从而检验该参数对输出结果的影响，但它却无法分析参数间的交互作用对模拟结果的影响[17]。全局敏感性分析则检验多个参数对模型输出结果的总影响，同时分析参数之间的交互作用对模拟结果的影响[18]。因此在使用多参数输入的非线性作物模型模拟时，越来越多的学者选择全局敏感性分析的方法对模型进行敏感性分析。目前常见的全局敏感性分析方法包括Sobol法、傅里叶振幅敏感性分析法（FAST）、可扩展傅里叶振幅敏感性分析法（EFAST）、LH-OAT法、普适似然不确定性估计法（GLUE）等[19]。目前基于可扩展傅里叶振幅敏感性分析法已应用于DASST[4]、EPIC[20]、WOFOST[6]等作物生长模型的全局敏感性分析中，但对于DNDC模型的敏感性分析研究较少。
基于上述情况，本研究以红枣为例，采用可扩展傅里叶振幅敏感性分析法，针对干旱胁迫和正常条件下的DNDC模型几类输入参数进行定量的全局敏感性分析。研究枣的作物参数、土壤因素在干旱胁迫和正常条件下对DNDC模型模拟枣产量的影响，旨在为DNDC枣类模型参数本地化、区域化以及未来与遥感数据的耦合提供选择敏感性分析变量的基础。
1 材料与方法
1.1 DNDC作物生长模型与试验数据
1.1.1 DNDC作物生长模型 DNDC模型是以模拟作物反硝化-分解反应为主的集合多种作物的地球生物化学过程的模型[21]。模型共包括两大部分，第一部分主要是通过气象因素、土壤因素、植被信息以及田间管理数据等对作物所处的土壤环境的动态变化进行模拟;第二部分通过输入的土壤环境因素预测模型输入参数对微生物活动影响，从而计算土壤系统中碳、氮元素的变化（主要为气体的排放）[10，22-24]。模型通过植物中碳元素的变化计算有机物的变化，模拟产量则通过最终果實内部的碳含量进行展示，其中模型默认作物产量的40%是果实内部的含碳量。本研究采用的模型为DNDC 95版本。
与其他模型不同，DNDC模型提供了简单易懂的作物类别开发界面，使用者可在数据库中建立新作物类别。红枣模型即为本研究在果树模型的基础上结合试验数据构造的新作物类别模型。
1.1.2 试验数据 本研究的试验区位于新疆昆玉第十四师224团地块。试验区共涵盖5块实验田，共包含采样点10个。
DNDC模型的输入数据包括气象参数、土壤参数、作物参数和田间管理数据。选取的作物参数包括：最大果实产量、全株生物量中果实比例、全株生物量中叶片比例、全株生物量中茎比例、全株生物量中根比例、果实中碳氮含量比例、茎中碳氮含量比例、叶中碳氮含量比例、根中碳氮含量比例、固氮指数、需水量、适宜温度、生长积温、管束结构指数、最大作物产量、来自土壤的氮量、来自大气的氮量、需氮量、植株最大高度、黏土含量、容质量、饱和导水率、田间持水率、萎蔫点、孔隙度、有机碳量、起始硝酸根含量、起始氨气含量、酸碱度、顶部均匀土层密度和下层土壤中有机碳（SOC）沉降速率。 研究所需的气象输入数据来自中国气象数据网的中国地面气候资料日值数据集（V3.0），共包含2015-2019年新疆和田地区逐日最高和最低气温，日降水量，日照时数。此格式为DNDC模型中特有的最简格式，可避免计算太阳辐射量带来的误差。
土壤数据主要来源于当地农业局实地测量数据和《新疆土种志》中数据，土壤主要为砂性土，透水性强，土壤肥力适中，有机质含量在1.0%～1.5%，pH值为8.05。测量数据主要为土壤容质量、饱和导水率、田间持水率、有机碳含量等。模型输入的田间管理数据主要来自于当地种植标准化操作手册和实地调查数据。实测产量数据来源于当地农业科学研究所多年的统计数据。
1.2 可扩展傅里叶振幅敏感性分析法
可扩展傅里叶振幅敏感性分析法 [25]是由Saltelli等结合Sobol法和傅里叶振幅敏感度测试法的优点所提出的基于模型方差分析的全局敏感性分析方法。EFAST方法的基本原理是模型的输入参数的变化对模型的输出结果的方差存在影响，通过计算每一个参数（Xi）的变化对输出结果（Y）的方差的贡献率得到敏感性[26]。其中输出结果Y=f（X）=f（x1，x2，x3…xn）。
对于一个多参数相互作用的模型，ST.i（参数xi的总敏感度）反映了Si（参数的直接贡献率）以及Sij～S1，2…i…k（参数间的互相作用）对模型输出总方差的贡献率，可表示为各阶敏感度之和：
ST.i=Si Sij Sijm … S1，2，…i…，k（1）
EFAST将敏感度划分为2个指标，一个为一阶敏感度即单参数的影响;另一个为全局敏感度即多个参数之间的相互作用的影响[27]。
DNDC模型包含多参数相互作用的过程，因此全局敏感性分析方法更适合DNDC模型的参数敏感性分析。EFAST方法通过对模型方差的分解，从而定量地得到每个参数的一阶敏感度及总敏感度。
1.3 试验设计
本研究的目的是为模型本地化选择敏感参数，试验设计方案如下：
在干旱胁迫和正常条件下对DNDC模型作物参数进行全局敏感性分析，并对比两者对模型最终模拟结果的影响，进一步调整模型参数。其中根据试验区域的实际情况，该地区的年降水总量平均为100 mm，因此大部分水分来源于农田灌溉，设置干旱胁迫灌溉量为正常灌溉量的75%和50%，在这种情况下分析红枣的19个作物参数和12个土壤参数，从中筛选最敏感参数，从而为模型本地化调整提供依据。
方案一：在正常水分条件和干旱胁迫条件下对DNDC模型的作物参数、土壤参数、田间管理参数进行全局敏感性分析。通过EFAST方法筛选出在正常条件和干旱胁迫条件下对模型有关键影响的敏感参数。
方案二：在方案一筛选出敏感参数的基础上，调整模型的参数，进一步优化输入数据，获取模拟结果。将调整前后的模拟结果与实测结果进行误差分析，通过相对误差这一指标进一步评价模型的精度。
研究采用的全局敏感性分析主要借助敏感性分析软件Simlab对模型的各个参数进行敏感性分析。该软件主要是利用蒙特卡洛方法对红枣的作物参数以及土壤参数进行随机采样，并利用全局敏感性分析法分解方差从而对输入参数进行敏感性分析。首先利用Simlab选择采样方法（本研究选取的为EFAST法，该方法认为分析结果有效时采样数需大于参数个数的65倍）[20]进行采样，生成输入参数 ，采样次数为4 000次。利用Python编写程序将参数编写为DNDC模型规定的输入格式，随后将参数输入到模型中得到模拟结果，将各参数和输出结果利用EFAST方法进行全局敏感性分析。本研究将全局敏感度大于0.1、一阶敏感度大于0.05的参数确定为敏感参数。
2 结果与分析
2.1 正常条件下和干旱胁迫下作物参数全局敏感性分析对比
本试验共针对DNDC模型中19个作物参数分别进行正常条件下和干旱胁迫条件下的敏感性分析。由图1可知，在干旱胁迫下和正常条件下的全局敏感参数筛选中，筛选的参数一致性较高，但参数的敏感度却不相同。其中正常条件下和75%水分条件下筛选出的敏感参数的敏感度由大到小为：全株生物量中果实比例
中图分类号： S165 .27 文献标识码： A 文章编号： 1000-4440（2021）03-0812-05
Adjustment of jujube yield parameters based on DNDC model under drought stress
WANG De-juan， SUI Li-chun， LIU Di， WANG Jin， YUAN Huan-huan， XU Jia-li
（School of Geological Engineering and Geomatics， Chang’an University， Xi’an 710054，China）
Key words： jujube;DNDC model;drought stress;sensitivity analysis;parameter adjustment
红枣作为南疆林果业发展的优势树种，在新疆绿洲生产、生态建设中起着重要作用[1]。新疆地区应用的估产方法主要有气象估计，经验估计等[2-3]。段丁丁等[4]利用DSSAT内部的马铃薯模型对马铃薯进行了区域化模拟估产研究;吴立峰[5]等通过Morris方法和EFAST方法对不同灌溉水平下CROPGRO棉花模型进行了敏感性分析和不确定性分析;纪甫江[6]等通过将CASA模型和WOFOST模型进行耦合从而提高了大豆单产估算的精度。而目前对于枣的产量估算研究较少，因此结合新疆地区的生态条件，利用更加精准的作物生长模型将成为在枣类估产中的有效方法。
经过多年的发展和应用，作物模型主要分为三大派系，分别是以ORYZA2000模型为代表的荷兰作物模型体系，以CERES模型为代表的美国作物模型体系，以CCSODS模型为代表的中国作物模型体系[7-9]。但在众多作物生长模型中用于果树产量模拟的模型较少，其中DNDC模型是由美国新罕布什尔大学的李长生教授团队研发的生物地球化学模型，该模型是农业生态系统中一系列控制碳和氮迁移转化的生物化学及地球化学反应机制的计算机模拟表达[10]。DNDC模型以日为步长，通过气候、土壤、植被和管理4个基本生态因子对作物生长进行机理性模拟，并通过DNDC模型特色的硝化反硝化作用模拟作物在生长过程中的碳、氮元素的变化，在进行产量预报的同时进一步实现碳、氮元素在地球生态系统中循环的模拟[11-13]。DNDC模型包含简单易操作的新作物类构建方法，并且对水分处理具有单独的页面，因此对枣树的产量模拟可操作性更高。
作物生长模型作为机理性模型需要大量的作物和实际的环境参数作为输入参数，但只有部分参数对模型的模拟结果起决定作用，因此参数的准确性将会直接影响模型模拟结果的真实性[14-15]。同时试验地区的特殊环境因素也成为参数敏感性分析的关键。新疆地区特殊的地理位置及天气条件，使得水分成为作物生长过程中不容忽略的问题。参数敏感性分析主要分为全局参数敏感性分析和局部参数敏感性分析[16]。局部敏感性分析是指在一定范围内调整某一个参数的取值从而检验该参数对输出结果的影响，但它却无法分析参数间的交互作用对模拟结果的影响[17]。全局敏感性分析则检验多个参数对模型输出结果的总影响，同时分析参数之间的交互作用对模拟结果的影响[18]。因此在使用多参数输入的非线性作物模型模拟时，越来越多的学者选择全局敏感性分析的方法对模型进行敏感性分析。目前常见的全局敏感性分析方法包括Sobol法、傅里叶振幅敏感性分析法（FAST）、可扩展傅里叶振幅敏感性分析法（EFAST）、LH-OAT法、普适似然不确定性估计法（GLUE）等[19]。目前基于可扩展傅里叶振幅敏感性分析法已应用于DASST[4]、EPIC[20]、WOFOST[6]等作物生长模型的全局敏感性分析中，但对于DNDC模型的敏感性分析研究较少。
基于上述情况，本研究以红枣为例，采用可扩展傅里叶振幅敏感性分析法，针对干旱胁迫和正常条件下的DNDC模型几类输入参数进行定量的全局敏感性分析。研究枣的作物参数、土壤因素在干旱胁迫和正常条件下对DNDC模型模拟枣产量的影响，旨在为DNDC枣类模型参数本地化、区域化以及未来与遥感数据的耦合提供选择敏感性分析变量的基础。
1 材料与方法
1.1 DNDC作物生长模型与试验数据
1.1.1 DNDC作物生长模型 DNDC模型是以模拟作物反硝化-分解反应为主的集合多种作物的地球生物化学过程的模型[21]。模型共包括两大部分，第一部分主要是通过气象因素、土壤因素、植被信息以及田间管理数据等对作物所处的土壤环境的动态变化进行模拟;第二部分通过输入的土壤环境因素预测模型输入参数对微生物活动影响，从而计算土壤系统中碳、氮元素的变化（主要为气体的排放）[10，22-24]。模型通过植物中碳元素的变化计算有机物的变化，模拟产量则通过最终果實内部的碳含量进行展示，其中模型默认作物产量的40%是果实内部的含碳量。本研究采用的模型为DNDC 95版本。
与其他模型不同，DNDC模型提供了简单易懂的作物类别开发界面，使用者可在数据库中建立新作物类别。红枣模型即为本研究在果树模型的基础上结合试验数据构造的新作物类别模型。
1.1.2 试验数据 本研究的试验区位于新疆昆玉第十四师224团地块。试验区共涵盖5块实验田，共包含采样点10个。
DNDC模型的输入数据包括气象参数、土壤参数、作物参数和田间管理数据。选取的作物参数包括：最大果实产量、全株生物量中果实比例、全株生物量中叶片比例、全株生物量中茎比例、全株生物量中根比例、果实中碳氮含量比例、茎中碳氮含量比例、叶中碳氮含量比例、根中碳氮含量比例、固氮指数、需水量、适宜温度、生长积温、管束结构指数、最大作物产量、来自土壤的氮量、来自大气的氮量、需氮量、植株最大高度、黏土含量、容质量、饱和导水率、田间持水率、萎蔫点、孔隙度、有机碳量、起始硝酸根含量、起始氨气含量、酸碱度、顶部均匀土层密度和下层土壤中有机碳（SOC）沉降速率。 研究所需的气象输入数据来自中国气象数据网的中国地面气候资料日值数据集（V3.0），共包含2015-2019年新疆和田地区逐日最高和最低气温，日降水量，日照时数。此格式为DNDC模型中特有的最简格式，可避免计算太阳辐射量带来的误差。
土壤数据主要来源于当地农业局实地测量数据和《新疆土种志》中数据，土壤主要为砂性土，透水性强，土壤肥力适中，有机质含量在1.0%～1.5%，pH值为8.05。测量数据主要为土壤容质量、饱和导水率、田间持水率、有机碳含量等。模型输入的田间管理数据主要来自于当地种植标准化操作手册和实地调查数据。实测产量数据来源于当地农业科学研究所多年的统计数据。
1.2 可扩展傅里叶振幅敏感性分析法
可扩展傅里叶振幅敏感性分析法 [25]是由Saltelli等结合Sobol法和傅里叶振幅敏感度测试法的优点所提出的基于模型方差分析的全局敏感性分析方法。EFAST方法的基本原理是模型的输入参数的变化对模型的输出结果的方差存在影响，通过计算每一个参数（Xi）的变化对输出结果（Y）的方差的贡献率得到敏感性[26]。其中输出结果Y=f（X）=f（x1，x2，x3…xn）。
对于一个多参数相互作用的模型，ST.i（参数xi的总敏感度）反映了Si（参数的直接贡献率）以及Sij～S1，2…i…k（参数间的互相作用）对模型输出总方差的贡献率，可表示为各阶敏感度之和：
ST.i=Si Sij Sijm … S1，2，…i…，k（1）
EFAST将敏感度划分为2个指标，一个为一阶敏感度即单参数的影响;另一个为全局敏感度即多个参数之间的相互作用的影响[27]。
DNDC模型包含多参数相互作用的过程，因此全局敏感性分析方法更适合DNDC模型的参数敏感性分析。EFAST方法通过对模型方差的分解，从而定量地得到每个参数的一阶敏感度及总敏感度。
1.3 试验设计
本研究的目的是为模型本地化选择敏感参数，试验设计方案如下：
在干旱胁迫和正常条件下对DNDC模型作物参数进行全局敏感性分析，并对比两者对模型最终模拟结果的影响，进一步调整模型参数。其中根据试验区域的实际情况，该地区的年降水总量平均为100 mm，因此大部分水分来源于农田灌溉，设置干旱胁迫灌溉量为正常灌溉量的75%和50%，在这种情况下分析红枣的19个作物参数和12个土壤参数，从中筛选最敏感参数，从而为模型本地化调整提供依据。
方案一：在正常水分条件和干旱胁迫条件下对DNDC模型的作物参数、土壤参数、田间管理参数进行全局敏感性分析。通过EFAST方法筛选出在正常条件和干旱胁迫条件下对模型有关键影响的敏感参数。
方案二：在方案一筛选出敏感参数的基础上，调整模型的参数，进一步优化输入数据，获取模拟结果。将调整前后的模拟结果与实测结果进行误差分析，通过相对误差这一指标进一步评价模型的精度。
研究采用的全局敏感性分析主要借助敏感性分析软件Simlab对模型的各个参数进行敏感性分析。该软件主要是利用蒙特卡洛方法对红枣的作物参数以及土壤参数进行随机采样，并利用全局敏感性分析法分解方差从而对输入参数进行敏感性分析。首先利用Simlab选择采样方法（本研究选取的为EFAST法，该方法认为分析结果有效时采样数需大于参数个数的65倍）[20]进行采样，生成输入参数 ，采样次数为4 000次。利用Python编写程序将参数编写为DNDC模型规定的输入格式，随后将参数输入到模型中得到模拟结果，将各参数和输出结果利用EFAST方法进行全局敏感性分析。本研究将全局敏感度大于0.1、一阶敏感度大于0.05的参数确定为敏感参数。
2 结果与分析
2.1 正常条件下和干旱胁迫下作物参数全局敏感性分析对比
本试验共针对DNDC模型中19个作物参数分别进行正常条件下和干旱胁迫条件下的敏感性分析。由图1可知，在干旱胁迫下和正常条件下的全局敏感参数筛选中，筛选的参数一致性较高，但参数的敏感度却不相同。其中正常条件下和75%水分条件下筛选出的敏感参数的敏感度由大到小为：全株生物量中果实比例