引言

菌核病是油菜生产中三大病害之首，从苗期到近熟期都可以发生，以中、后期发病最为普遍，叶、茎、荚都可被害，其中茎部被害损失最重(李运良等，2009)。菌核病可引起油菜植株早枯，角果减少，籽粒皱瘪，严重影响油菜产量和品质，油菜感病后可减产10%~70%，含油量降低1%~5%(陈欣欣等, 2019; 齐永霞等, 2006)。油菜是安徽省种植面积最大的油料作物，常年种植面积为90万~95万hm²(刘瑞娜等, 2016)，安徽也是油菜菌核病的重发区，该病害已成为制约当地油菜生产的主要因子(许大凤等, 2014)。开展油菜菌核病预报预测和监测评估，对做好油菜菌核病防治工作，优化油菜种植布局都具有重要意义。

油菜菌核病是一种典型的“气象型”病害，其发生程度与气象条件密切相关，徐森富等(2012)研究表明“浙油50”开花期是菌核病主要流行期，其平均气温和相对湿度对油菜菌核病发生具有显著影响；易红娟等(2008)指出江苏省南通市通州区影响油菜菌核病流行程度的主要气象因子是油菜开花至终花期的雨湿条件；张立良等(2000)研究认为，安徽和县地区花期降雨量直接影响病害的流行，5月中旬的降雨量影响病害的发病程度；张莉等(2015)分析表明3月和4月的雨日和雨量直接决定了湖北荆州市当年油菜菌核病的发病程度；长江中下游地区一般4月降水偏多，当年菌核病偏重发生(陆均天, 2002; 叶殿秀, 2003)。区域间菌核病发病关键期存在一定差异性，但降水日数、降水量、相对湿度等与降水相关的气象要素是影响各地区油菜菌核病发生程度的主要气象因子。目前关于油菜菌核病发生气象等级的研究多集中在预测预报方面(朱金良等, 2012; 罗泽青等, 2016; 郑露等, 2017; 居为民等, 2000; Ficke et al, 2018; Naresh, 2019)，而针对监测评估的研究相对较少。

在利用气象条件开展油菜菌核病发生程度的研究中，多以降水日数、降水量、相对湿度等气象要素为主，综合考虑同一时段降水日数和降水强度对菌核病影响的研究尚不多见，且由于温度与油菜菌核病发生程度相关性不显著常被忽略(易红娟等, 2008; 常彭阳, 1999;朱金良等, 2012)，常彭阳, 1999但刘勇等(2019)认为低温、多雨天气会直接延长油菜的开花时间，为油菜菌核病的侵染、发生、蔓延和严重危害带来非常有利的气象条件。陈士华等(2005)也指出菌核病发病过程中温度影响子囊孢子释放侵入及菌丝再侵染。所以油菜菌核病发生流行不仅需要充沛的降水条件，适宜的温度也能促进病害的发生。在前人研究基础上，综合考虑降水日数、降水等级和温度对菌核病的影响，形成综合气象条件指数，并基于综合气象条件指数建立了油菜菌核病气象等级评价模型。以期为安徽省油菜菌核病的监测评估提供技术支撑，进而为安徽省油菜的种植结构调整和菌核病的科学防控提供参考依据。

1 资料与方法 1.1 数据资料

研究以茎病株率作为油菜菌核病发生程度指标。1995-2019年安徽省池州市和桐城市油菜菌核病茎病株率资料分别来自当地植物保护站，对应年份的气象资料来源于安徽省气象信息中心，主要包括逐日平均气温(单位：℃)和降水量(单位：mm)。1995-2015年的资料作为历史样本用于评估模型建立，2016-2019年的资料作为独立样本用于模型模拟验证。

1.2 油菜菌核病气象等级划分

根据安徽省植保部门关于油菜菌核病发生程度划分标准，结合安徽省油菜菌核病发生特点，将菌核病发生气象等级划分为5个等级，即茎病株率 < 10%为1级，表示气象条件不适宜，对应油菜菌核病轻发生；10%≤茎病株率 < 20%为2级，表示气象条件基本适宜，对应中等偏轻发生；20%≤茎病株率 < 30%为3级，表示气象条件较适宜，对应中等发生；30%≤茎病株率 < 40%为4级，表示气象条件适宜，对应中等偏重发生；茎病株率≥40%为5级，表示气象条件非常适宜，对应大发生。

1.3 研究方法 1.3.1 研究思路

油菜菌核病的发生流行与气象条件密切相关，油菜开花期为菌核病病原侵染的敏感期，花期降水多利于子囊孢子的释放、侵染和菌丝再侵染(张立良等, 2000)。角果发育期为茎秆发病的敏感期，冯兰萍等(1999)通过田间调查发现油菜菌核病从油菜终花期开始表现症状，终花后7~10 d发病率明显上升，茎秆发病在终花后20~25 d达到发病高峰；杜晓宇等(2009)指出油菜茎病株发生一般出现在终花后7~15 d；苏跃和吴沿友(2016)研究表明油菜茎部菌核病发病率从盛花期开始到成熟期均呈上升趋势。所以初步确定油菜开花至成熟前为菌核病茎发病关键期。相关研究表明雨日、雨量是影响油菜菌核病发生的重要因素(刘振忠等, 1999; 陈士华等, 2005; 张莉等, 2015)，气温可影响子囊孢子的侵入和菌丝再侵染。

根据油菜菌核病发生流行与气象条件的关系，以降水日数作为影响菌核病的主要气象因子，通过相关分析确定菌核病发病关键期，同时考虑降水等级和温度对菌核病的影响，引入雨量系数和温度系数形成综合气象条件指数。利用回归分析，建立基于综合气象条件指数的油菜菌核病气象等级评估模型。

1.3.2 雨量系数计算方法

由于油菜菌核病发病过程中不同降水量所起的作用不同，参照天气预报中24 h降水量级标准，将菌核病发生关键期日降水量(R)分为相应6个不同等级，并对不同等级的降水赋予不同灾害系数(f_R)，计算方法如下：

$ f_{R}= \begin{cases}r_{1} & 0.1 \leqslant R<10 \\ r_{2} & 10 \leqslant R<25 \\ r_{3} & 25 \leqslant R<50 \\ r_{4} & 50 \leqslant R<100 \\ r_{5} & 100 \leqslant R<250 \\ r_{6} & R \geqslant 250\end{cases} $

(1)

式中: R为油菜菌核病发生关键期日降水量，单位：mm；r₁、r₂、r₃、r₄、r₅、r₆为不同降水等级对应的雨量系数。

1.3.3 温度系数计算方法

适温高湿的环境条件是油菜菌核病发生的重要因素，适宜的温度利于油菜菌核病子囊孢子侵染和菌丝再侵染，温度过高、过低对菌丝生长均有抑制作用。根据油菜菌核病发生流行对温度的需求特点，采用下列函数计算温度系数。

$ f_{T}= \begin{cases}0 & T<T_{1} \\ \frac{T-T_{1}}{T_{\mathrm{ol}}-T_{1}} & T_{1} \leqslant T<T_{\mathrm{ol}} \\ 1 & T_{\mathrm{ol}} \leqslant T<T_{\mathrm{ou}} \\ \frac{T_{\mathrm{u}}-T}{T_{\mathrm{u}}-T_{\mathrm{ou}}} & T_{\mathrm{ou}} \leqslant T<T_{\mathrm{u}} \\ 0 & T \geqslant T_{\mathrm{u}}\end{cases} $

(2)

式中: f_T为温度系数；T为菌核病发生关键期日平均气温；T_l为菌核病发育下限温度，低于这一温度时，温度系数为0；T_u为菌核病发育上限温度，超过这一温度时，温度系数为0；T_ol、T_ou分别为菌核病发育最适温度的下限和上限。

1.3.4 综合气象条件指数

综合气象条件指数是反映降水和温度对病害的综合影响程度，为关键期逐日雨量系数与温度系数乘积的累计值，计算方法如下：

$ R_{\mathrm{C}}=\sum\limits_{i=1}^{n} f_{R_{i}} f_{T_{i}} $

(3)

式中: R_C为影响油菜菌核病发生的综合气象条件指数; i为菌核病发生关键生育期降水日数，i=1，2，…，n; f_Ri为关键期雨量系数; f_Ti为关键期温度系数。综合气象条件指数越大，即降水和温度相互作用越好，气象条件越适宜菌核病的发生流行; 反之，则不适宜菌核病的发生流行。

1.4 模型检验方法

油菜菌核病气象等级评估模型的回代检验和模拟验证的准确率参照王志伟等(2010)提出的计算方法，即模型模拟等级与实际等级相差0~0.5级为完全符合，准确率计为100%；相差0.5~1.0级为比较符合，计为80%；相差1.0~1.5级为基本符合，计为70%；相差1.5级以上为不符合，计为0；历史符合率为各年符合程度准确率的平均值。

1.5 数据处理

本研究采用SPSS 19.0统计软件对数据进行相关分析和回归分析，利用Visual Basic语言对雨量系数进行循环处理。

2 结果与分析 2.1 油菜菌核病发生关键期

根据安徽省农业气象作物发育期观测资料，池州市和桐城市油菜一般3月上旬始花，5月中旬成熟，但不同年份发育期略有差异。为确定油菜菌核病发生关键期，以旬为单位分别统计3月上旬至5月中旬逐5旬、6旬、7旬、8旬的降水日数，再与菌核病发生气象等级进行相关分析。结果表明，池州市和桐城市油菜菌核病发生气象等级均与3月下旬至5月上旬的降水日数相关系数最大，分别为0.618 1和0.679 1(表 1)，均通过显著性水平检验(P < 0.01)。从常年发育期看，3月下旬至5月上旬池州和市桐城市油菜处于盛花期至角果成熟期，与油菜菌核病茎秆发病敏感时段较为一致，所以确定池州市和桐城市油菜菌核病发生关键期为3月下旬至5月上旬。

2.2 油菜菌核病气象等级评估模型 2.2.1 雨量系数

通过对1995-2015年池州市和桐城市3月下旬至5月上旬日降水等级分布进行统计分析，由图 1可以看出，池州市和桐城市菌核病发生关键期的降水以小雨和中雨为主，小雨与中雨日数之和分别占两市总雨日的91.8%和85.6%，暴雨及以上等级的雨日较少。考虑不同降水等级样本数量，本研究将池州市和桐城市油菜菌核病发生关键期日降水量划分为3个等级，即小雨、中雨、大雨及以上。

将不同降水等级在油菜菌核病发生过程中所起到的作用定义为雨量系数f_R，关键期雨量系数之和定义为实效雨日。利用Visual Basic语言作为计算工具，将小雨、中雨和大雨及以上对应的雨量系数设定在一定范围内，以0.1作为步长进行循环，计算出不同雨量系数对应的实效雨日，再与菌核病发生气象等级进行相关分析，其中相关系数最大值对应的实效雨日即为不同降水等级的雨量系数。通过上述方法，得到池州市和桐城市油菜菌核病发病关键期雨量系数如下：

$ {\rm{池洲}}:{f_R} = \left\{ {\begin{array}{*{20}{l}} {0.7}&{0.1 \le R < 10}\\ {1.0}&{10 \le R < 25}\\ {2.6}&{R \ge 25} \end{array}} \right. $

(4)

$ 桐城:{f_R} = \left\{ {\begin{array}{*{20}{l}} {0.6}&{0.1 \le R < 10}\\ {1.1}&{10 \le R < 25}\\ {1.9}&{R \ge 25} \end{array}} \right. $

(5)

由表 2可以看出，池州市和桐城市实效雨日与菌核病气象等级的相关性分别为0.674 1和0.755 0，均高于降水量、降水日数与菌核病气象等级的相关性，说明雨量系数能更好地反映不同降水等级对油菜菌核病的影响。

2.2.2 温度系数

倪守延等(1995)研究表明，油菜菌核病菌丝的生长温度范围为0~35℃，最适温度为23~25℃，所以式(2)中T_l、T_ol、T_ou、T_u分别取值为0、23、25、35℃，油菜菌核病温度系数计算式(2)可改为公式(6)。当菌核病发生关键期降水日日平均气温低于0℃或高于35℃时，温度系数f_T值为0；当日平均气温在23~25℃时，f_T值为1；当日平均气温在0~23℃时，随着温度的升高，f_T逐渐增加；当日平均气温在25~35℃时，随着温度的升高，f_T逐渐减小。

$ f_{T}= \begin{cases}0 & T<0 \\ \frac{T-0}{23-0} & 0 \leqslant T<23 \\ 1 & 23 \leqslant T<25 \\ \frac{35-T}{35-25} & 25 \leqslant T<35 \\ 0 & T \geqslant 35\end{cases} $

(6)

2.2.3 综合气象条件指数

根据池州市和桐城市1995-2015年3月下旬至5月上旬逐日降水和气温资料，利用式(3)~式(6) 计算出不同年份菌核病发生综合气象条件指数。结果显示，池州市和桐城市综合气象条件指数与菌核病气象等级相关系数分别为0.810 5和0.784 0，均高于降水日数、实效雨日与菌核病发生气象等级的相关性，所以本研究提出的雨量系数和温度系数的观点不仅符合油菜菌核病发生流行生理学特性，且综合气象条件指数能更好地反映降水日数、降水等级和温度对菌核病发生的综合影响，对提高菌核病气象等级评估精准度具有重要意义。

2.2.4 气象等级评估模型

运用SPSS 19.0统计软件中的曲线回归方法，以综合气象条件指数为自变量，以油菜菌核病发生气象等级为因变量，建立池州市最优回归模型：Y₁=3.017 6-0.784 9R_C1+0.081 11R_C1²-0.001 894R_C1³，相关系数r=0.815 9，通过显著性水平检验(P < 0.01)，拟合曲线见图 2a；桐城市最优回归模型：Y₂=-3.159 2+1.199 2R_C2-0.087 4R_C2²+0.002 637R_C2³，相关系数r=0.787 6，通过显著性水平检验(P < 0.01)，拟合曲线见图 2b。式中，Y₁和Y₂分别为池州市和桐城市油菜菌核病发生气象等级，R_C1和R_C2分别为池州市和桐城市综合气象条件指数。

2.3 模型检验 2.3.1 模型回代检验

利用1995-2015年池州市和桐城市油菜菌核病茎病株率资料和气象观测资料，对菌核病气象等级评估模型进行回代检验。由表 3可见，回代结果为池州市模型回代等级与实际等级平均误差为0.604，回代平均准确率为86.2%，21个样本中回代准确率为100%的有9 a，准确率为80%的7 a，准确率为70%的5 a。桐城市模型回代等级与实际等级平均误差为0.685，回代平均准确率为83.8%，21个样本中回代准确率为100%的有5 a，准确率为80%的14 a，准确率为70%的2 a。从回代检验结果看，本文建立的池州市和桐城市的油菜菌核病气象条件等级评估模型模拟效果较好，总体上能反映气象条件对菌核病发生的影响。

2.3.2 模型模拟验证

2016-2019年利用池州市和桐城市的油菜菌核病茎病株率和对应的逐日气象观测资料，对油菜菌核病气象等级评估模型进行模拟验证。由表 4可以看出，2016-2019年池州市和桐城市模型模拟结果与实际发生情况平均误差分别为0.582和0.635，平均模拟准确率均为85.0%。其中池州市模拟结果与实际情况相比，2年偏多，2年偏少，误差幅度在0.345~0.740，桐城市模拟结果为3年偏多，1年略偏少，误差幅度在0.121~0.890。

3 结论与讨论

本文利用1995-2015年池州市、桐城市菌核病茎病株率资料和同期气象观测数据，以降水日数作为主要影响因子，通过相关分析，明确池州市和桐城市油菜菌核病茎秆发病关键期均为3月下旬至5月上旬，该时段池州市和桐城市油菜主要处于盛花至角果成熟期，与冯兰萍等(1999)、蒋耀培等(2000)关于油菜菌核病茎秆发病敏感期观点基本一致。气象条件对菌核病的影响是一个连续的过程，油菜开花期是菌核病侵染的关键时期(朱金良等, 2012; 陈庭华, 1990)，降水多利于子囊孢子的释放、侵染和菌丝再侵染，致使叶病株率高(张立良等, 2000)，染病的叶片又为后期茎秆发病提供了大量的菌源。油菜终花至成熟前是菌核病茎秆显症的重要时期(冯兰萍等, 1999; 杜晓宇等, 2009; 苏跃和吴沿友, 2016)，该时段的降水条件是影响茎秆发病程度的主要因素(蒋耀培等, 2000)。所以研究确定的关键期符合油菜菌核病茎秆发病的生理学特性。

基于3月下旬至5月上旬的降水和气温要素，建立了沿江地区油菜菌核病气象等级评估模型，通过回代检验和模拟验证，池州和桐城地区的回代准确率和模拟准确率均在80%以上，基本满足了利用气象条件开展菌核病发生程度的监测评估需求。模型中降水和温度作为气象部门常规观测和预报要素，在未来中长期预报水平达到一定程度精细化条件下，建立的评估模型也可用于菌核病发生气象等级的预测预报。降水的连续性在“气象型”病害的预测评估模型中起到了较好的作用(苏荣瑞等, 2015; 张旭晖等, 2009)，但本研究建立的菌核病气象等级评估模型未引入降水连续性。因为统计分析发现，引入降水连续性不能提高模型拟合的准确率，其原因可能是油菜菌核病发生关键期降水达到一定日数，降水分散与集中对菌核病的影响无明显差别。

油菜菌核病发生和流行除与气象条件密切相关外，还受菌源、品种和栽培条件等影响。Dunker and von Tiedemann(2004)田间试验证实油菜菌核病发生程度与种植密度和播期相关。不同连作类型田的油菜菌核病发生特点也存在差异，苏永新等(2013)研究发现棉花-油菜旱地连作田花期菌核萌发数量较水稻-油菜水田连作油菜田多62.5%，而杜晓宇等(2009)发现油菜菌核病发生偏重年份，旱地和水田发病程度趋于一致，但在中等及以下发生年份，水稻-油菜连作田发生程度普遍重于旱地连作发生程度。不同类型和品种以芥菜型抗性较好，甘蓝型次之，白菜型最感病(李明桃, 2012)。此外田间菌核密度、子囊盘密度与茎秆病株率呈极显著正相关关系(秦虎强等, 2018)。建立的基于气象条件的菌核病气象等级评估模型，仅能从区域尺度上评估病害发生程度平均状况。所以开展精细化田块尺度的油菜菌核病评估，需要结合田间气象条件监测、品种抗性、菌源量、栽培管理等多因素建立综合评估模型，以提高评估精细化和精准化程度。