日光温室黄瓜霜霉病初侵染阶段关键预测因子的筛选及验证

纪涛 刘慧英 许建平 柳瑞 刘冉 李明

摘要：温室黄瓜霜霉病的暴发依赖于环境信息和栽培管理措施等多种因素的相互作用，生产上测报依赖各种经验模型，但输入因子种类较多，需要进行简化处理。笔者基于田间调查试验，采用主成分分析的方法，从14组日光温室黄瓜霜霉病初侵染阶段的预测因子中，筛选出了反映湿度综合信息、温度综合信息和温室管理措施的3个主成分，累计贡献率达到80.76%，并结合前人的研究成果构建经验模型。该模型对黄瓜霜霉病初侵染阶段的预测效果较好（R²=0.94），可为日光温室黄瓜霜霉病的提早防治提供决策参考。

关键词：黄瓜;霜霉病菌;预警系统;主成分分析;经验模型

由古巴假霜霉菌[Pseudoperonospora cubensis（Berk.&Curt.;） Rostov.]引起的黄瓜霜霉病，是目前葫芦科作物危害性最严重的叶部病害之一。该病害具有潜育期短、发展迅速、再侵染频繁的特点。在美国、中国、日本、以色列和一些欧洲国家造成严重的经济损失。在我国，设施作物栽培在全国范围内大面积推广，日光温室相对封闭的结构特点，极易满足黄瓜霜霉病对“中温高湿”的环境要求，一旦发病未及时处理，后期防治非常困难。传统的防治主要依赖杀菌剂的定期使用，通常7～14d就要施药1次，而且见到症状后还要增加施药次数或用量，由于温室黄瓜采收频繁，容易造成产品农药残留[5]，因此必须做好对初侵染的预测和防治工作。

近年来国内外众多学者开发了以环境因素为驱动的病害預测模型来指导生产，对病害发生时间和风险进行预测，通过生态防治，在减少施药的基础上控制病害的发展。如何自福等人运用距始病期的时间、前10d累积相对湿度与气温的比值、前10d累积雨日与气温的比值、前10d累积雨量与相对湿度的比值、前10d累积雨日与相对湿度的比值等多个输入因子构建模型，对广州地区黄瓜霜霉病的发展进行预测。徐宁通过接种日至发病日时间、标准累积温度、标准累积湿度温、标准累积湿度、标准累积温湿度积之倒数等多个变量分别构建了塑料大棚黄瓜白粉病、霜霉病预测模型。但预测过程中常会面临模型输入变量过多造成某些重要信息被覆盖的现象，通过主成分分析对预测因子适当的筛选和化简，能够使模型运行更加方便、迅速。主成分分析就是将多个变量通过线性变换，剔除具有相关性的变量，从而筛选出少数重要综合变量的多元统计分析方法。其原理是通过数学降维处理，设法将多个原变量重新组合成新的相互无关的几个综合变量，这些新变量尽可能多地保留原来较多变量所反映的信息。该方法在植物病害管理方面具有很好的应用。

寄主、病原、环境和栽培措施共同构成了温室黄瓜病害的四面体，病害的发生、发展受到这四类因素的共同作用。鉴于寄主、病原一般在生产上比较稳定，本研究仅考虑了易变、易控的环境和栽培管理措施作为预测因子筛选的范围，利用SPSS.19进行主成分分析，对影响温室黄瓜霜霉病初侵染阶段的多种因素进行筛选，为温室黄瓜霜霉病预警系统的构建提供参考。

1 材料与方法

1.1 温室黄瓜种植

试验时间为2017年9-11月，试验地点位于北京市昌平区小汤山精准农业试验基地（40.18°N，116.47°E）。试验在（11号）日光温室中进行，温室规格为30m×7m，由聚乙烯薄膜覆盖，东西延长坐北朝南。于3月种植霜霉病中抗品种‘京研迷你2号，并建立生产档案，如实记录灌溉、施肥、夜间通风等栽培管理措施。

1.2 温室环境信息采集

在调查对象温室中设置自动气象站，每30min对温室内的温度、相对湿度、露点温度等数据实时观测并储存。室外也配有相应的自动气象站，每h对温度、相对湿度、雨量等天气要素进行记录。数据通过设备上的无线传输模块和手动下载相结合来收集。

1.3 霜霉病调查

每次调查在清晨露水未干时进行，定植后每天进行全棚普查，直到典型霜霉病早期症状（叶片正面为边缘模糊的黄褐色病斑，叶片背面为多角形水渍状病斑）出现，记录首次发病的日期（本试验初显症日期为10月1-2日）。此后采取对角线5点取样，定点定株调查，每点5～10株，周期改为3～4d调查1次，统计发病率并根据GB/T 17980.26-2000对病害严重度进行分级和记录。并在调查期间记录浇水情况、夜间通风等生产管理措施。

1.4 数据记录

选取发病前15d温室内外的环境信息和生产管理信息，X旧平均气温/℃（24h温度的平均值）、X₂日平均相对湿度/%（24h RH的平均值）、X₃气温-露点差/℃、X₄相对湿度90%的时间/h、X₅相对湿度≥80%的时间/h、X₆温度介于15～20℃的时间/h、X₇温度介于20～25℃的时间/h、X₈夜间温度/℃（20：00-8：00）、X₉夜间相对湿度/%（20：00-8：00）、X₁₀距上一次浇水日的时间/d、X₁₁每次浇水时长/h、X₁₂室外平均温度/℃、X₁₃室外降雨状况（降雨记作“1”、未降雨记作“0”）、X₁₄距上一次夜间通风的时间/d，作为日光温室黄瓜霜霉病初侵染阶段的预测因子进行主成分分析。参数计算采用Excel 2007进行。

1.5 数据分析

进行主成分分析的主要步骤如下：收集并选取需进行筛选的变量和数据;数据标准化处理;相关性分析;确定主成分个数和各个变量的载荷量。

1.6 经验模型构建

经验模型是一种传统的农业预测手段，具有输入简单、运行迅速、应用便捷的特点。笔者通过定性的方法从己有文献中归纳前人的研究成果和生产经验，并结合预测因子筛选的结果提出日光温室黄瓜霜霉病初侵染预测规则。

1.7 模型验证

选用2006-2007年北京地区日光温室黄瓜霜霉病实际发生情况作为模型的验证数据（4栋温室，每栋温室选择5个调查点，共20个调查点）。对模拟值和观测进行线性回归分析，计算线性回归斜率k、回归截距b以及决定系数R²。另外，为了进一步对模型的准确度、精度、偏差率等指标作出评价，分别引入拟合指数（WAI），置信指数（CI）来描述模型预测的准确度和精度;采用平均偏离差（MBE）、平均离差（MAE）来评估模型的平均误差幅度和偏离方向计算均方根误差（RUSE），RUSE的值越小，说明模拟值与观测值之间的偏差越小，模型的精度越高。计算方法如下：

式中：Oi—觀测值;Si—模拟值;Qi—观测值的离均差;Si—模拟值的离均差;n—样本总量。

2 结果与分析

2.1 数据的标准化处理

在进行分析前，为了消除各监测数据原始量纲和数量级不相同的影响，要对原始变量数据进行标准化的处理，方法是对同一变量减去其均值再除以其标准差，得到这14个原始变量标准化后的标准变量Z，然后采用SPSS 19.0中的“降维”过程进行主成分分析。

2.2 相关分析

通过相关分析（表1）可以发现，一些变量之间相关系数较大，例如日均温度X₁与RH≥80%的时间X₅、RH≥90%的时间X₆、浇水时长X₁₁等预测因子存在较强的相关性，其相关系数分别为-0.81、-0.69、0.83;平均相对湿度X₂与气温-露点差X₃的相关系数也达到了-1.0，说明这些变量存在信息上的重叠，较适合主成分分析方法。

2.3 主成分分析

主成分是以最少的个数反映尽量多的信息为原则，主成分个数提取原则有3个：①只取λ>1的特征值对应的主成分;②累计百分比达到80%～85%以上的λ值对应的主成分;③根据特征根变化的突变点决定主成分的数量。如表2所示，第1至第3个主成分的贡献率分别为44.45%、23.52%和12.79%，此时累计贡献率己达到80.76%，介于80%～85%，说明前3个主成分能够反映原始变量提供约80.76%的信息。因此，选取前3个特征值，并计算相应的主成分载荷（表3）。

主成分载荷矩阵体现了各变量与主成分之间的紧密度，某一主成分与某一变量的载荷系数的绝对值越大，则该主成分与变量之间的联系越紧密，因此，从主成分贡献率及载荷情况可以看出影响黄瓜霜霉病初侵染发展的预测因子类型及其所占比重。

第1主成分占总载荷量的44.45%，由表3可以看出，RH≥80%所持续的时间（0.960）、RH≥90%所持续的时间（0.959）、日平均相对湿度（0.929）、气温一露点差（-0.929）、夜间相对湿度（0.854）具有较大的载荷量，且依次减小;可以把第1主成分看成是由RH≥80%的时间、日平均相对湿度、气温一露点差、夜间相对湿度所反映的温室内湿度信息综合指标。

第2主成分占总荷载量的23.52%，其中夜间温度（0.867）、20～25℃的时间（0.785）所占的载荷量较大;因此第2主成分主要反映的是由夜间温度、20～25℃的时间所构成的温度信息综合指标。

第3主成分占总载荷量的12.79%，其中浇水时长和距上一次浇水的时间具有较大的载荷量，分别为-0.725和0.707，这表明第3主成分与温室的灌溉措施关系密切。

由结果可以分析出，从以上14组预测因子中，筛选出了反映温室湿度综合信息、温度信息和温室管理措施的3项综合指标（即3个主成分）。其中对日光温室黄瓜霜霉病初侵染影响最大的因子分别是RH≥80%所持续的时间、夜间温度以及每次浇水时长。

2.4 经验模型构建

根据前人的经验：赵胜荣等发现当RH≥80%的累计时间达到140h，霜霉病的发病率约为20%，病情指数的发展也与RH≥80%累积时间呈显著正相关;乔晓军发现夜间相对湿度大于83%，发病率约为28%;且霜霉病的侵染一般发生在夜间，侵染最适温度为15～20℃。因此，笔者结合前人的研究成果、生产经验，并综合预测因子筛选的结果，归纳总结出以下适宜温室黄瓜霜霉病初侵染阶段的预警规则：①自定植1周后开始RH≥80%的累积时间达到200h左右;②夜间相对湿度≥83%，夜间温度介于15～20℃;③若夜间或傍晚灌水，增加发病概率;满足以上条件，病害将会在3d内发生。

2.5 模型验证

采用2006-2007年北京地区日光温室黄瓜霜霉病发生情况作为模型的验证数据。如表4所示，20个调查点中满足生产要求（即命中）的预测频数为16次，占总体样本的80%;预测发病日期在实际发病日期之后（即延迟漏报）的预测频数为4次，占总体样本的20%，模型预测的准确率较高。

将预测值与实际值进行回归分析，进一步验证模型的预测质量。如表5所示，线性回归的决定系数R²=0.94，置信指数C=0.72，拟合指数W=0.75其绝对值均趋近于1，表明模型拟合效果较好;MBE=-0.7<0说明该模型能够提供早于实际发病日的预测;预测值与实际值相比平均离差MAE在1d左右，基于1：1散点图，均方根误差RMSE=1.34，说明模型总体预测质量较好，但仍然需要对延迟、漏报的现象进行优化。

3 讨论与结论

黄瓜霜霉病的发生、发展与温、湿度等环境因素关系密切，高湿和叶片上的游离水是抱囊萌发、形成初级侵染结构的必要条件;而温度则决定了病害发展的速率和严重程度。叶片自由水存在的条件下，抱子囊5～32℃均可萌发，适温在15～20℃;当接种物浓度较高时，最适温度15～20℃下保湿2h，便可完成侵染过程;如果低于巧℃或高于28℃，则不利于病害的发生发展;且温、湿度对抱子囊的存活也有显著影响，在相同的温度条件下，与45%和100%相比，75%左右的中等湿度抱子囊的存活率最低。高湿是发病的最重要条件，或叶片湿润时间6h以上可诱导抱子囊形成，且相对湿度越大，产抱越多;如果相对湿度低于60%时，则不利于病害的发生。相对湿度在93%以上时，发病严重。若阴天或雨天，将迫使种植者减少通风时间以维持棚内温度，这将产生较高的相对湿度和较长的叶片湿润时间，导致黄瓜霜霉病发生严重;若浇水次数多，或采用大水漫灌、喷灌、夜间灌水等灌溉方式，也造成棚内相对湿度增大，利于霜霉病发生、发展。采取膜下滴灌或膜下沟灌则能显著降低黄瓜霜霉病的危害程度。

通过对日光温室黄瓜霜霉病初侵染阶段预测因子的分析，可以看出利用主成分分析的方法能充分的展现出原始数据所包含的信息，通过主成分贡献率来考量各个变量重要性，具有基于数据的客观性和全面性，能够全面地反映各变量对病害发生发展的影响程度。但是病害的发生与否，除了环境条件和管理措施外，还与病原物、易感寄主等密切相关。因此，接下来的研究中，还需采取适当的方法，在预测模型中逐步加入病原、寄主生育时期、施肥频率、施肥量、温室揭盖帘情况、累积太阳辐射以及累积有效温度等因素综合考虑，以提高预测的准确度。

笔者通过主成分分析，筛选出了反映温室湿度综合信息、温度信息和温室管理措施的3项综合指标，为预测模型的建立及预警系统的开发奠定了基础。根据主成分分析结果并结合前人的生产经验提出黄瓜霜霉病预警规则：自定植后1周开始，RH≥80%的累积时间达到200h;夜间平均RH≥83%，夜间平均温度介于15～20℃;若当天灌溉时间较晚，则霜霉病可能发生。笔者使用了2006-2007年北京地区实际发病情况作为模型的验证数据，预测效果较好，但模型的适用范围还有待进一步验证。在今后的研究工作中，应尽量使用不同地区，不同年份的实际发病数据进行验证，从而提高模型的普适性。