南充市水稻稻瘟病区划和发生流行规律研究

彭昌家,白体坤,丁攀,冯礼斌,郭建全,杨宇衡,尹怀中,龙维国,文旭,肖立,崔德敏,苟建华,李鸿韬,何海燕,郑艳

中国农学通报. 2016, 32(6): 182-192

PDF(2789 KB)
PDF(2789 KB)
中国农学通报 ›› 2016, Vol. 32 ›› Issue (6) : 182-192. DOI: 10.11924/j.issn.1000-6850.casb15090025
植物保护 农药

南充市水稻稻瘟病区划和发生流行规律研究

  • 彭昌家1,白体坤1,丁攀1,冯礼斌1,郭建全2,杨宇衡3,尹怀中4,龙维国5,文旭6,肖立7,崔德敏8,苟建华9,李鸿韬9,何海燕10,郑艳5
作者信息 +

Regionalization and Epidemic Regularity of Rice Blast in Nanchong City

  • Peng Changjia1, Bai Tikun1, Ding Pan1, Feng Libin1, Guo Jianquan2, Yang Yuheng3, Yin Huaizhong4,Long Weiguo5, Wen Xu6, Xiao Li7, Cui Demin8, Gou Jianhua9, Li Hongtao9, He Haiyan10, Zheng Yan5
Author information +
History +

摘要

为提高南充市水稻稻瘟病监测预警和综合防控水平,切实减轻稻瘟病危害,采用GPS和GIS定位、系统监测和普查、品种抗稻瘟性鉴定、稻瘟病菌生理小种监测和气象资料分析等方法,开展了南充市水稻稻瘟病发生流行区划、精准勘界和发生流行规律研究。首次明确了南充市稻瘟病的流行区划和精准勘界,探明了水稻稻瘟病菌越冬后的初始侵染源(即始见期)地点、发生特点、流行趋势和特点,明确了南充市水稻稻瘟病在田间存在4个流行高峰期,其中第1个流行高峰期发生面积对当年病害流行起着决定作用,5月底的累计发生面积与年发生面积成正相关,相关系数为0.8172;建立了孕穗末期病田率和蜡熟期病株率同年发生面积预测数学模型。对其成因探讨发现,导致南充市近年水稻稻瘟病重发及流行的内因是水稻品种抗稻瘟能力的降低或丧失、稻瘟病菌致病力强的生理小种增加与更强致病力的小种出现;温度适宜和阴雨寡照偏多是造成稻瘟病在南充市流行的外因。2001—2015年,全市水稻稻瘟病短期预报准确率达100%,中、长期预报准确率分别达到98%和95%以上,比1997年以前提高5~15个百分点,促进稻瘟病防效提高15%~30%。

Abstract

In order to improve the level of monitoring, early warning and comprehensive control of rice blast in Nanchong City, and reduce the harm of rice blast, the methods of GPS and GIS location, systematical monitoring and survey, rice blast resistance variety identification, physiologic variety monitoring and weather data analysis were adopted to study the regionalization, accurate boundary settlement and occurrence and epidemic regularity of the disease. In this study, regionalization and boundary settlement were analyzed for the first time, the primary infection location, occurrence characteristics, trends and characteristics of the disease were verified. There were four epidemic peak periods and the damage areas of the first peak period played a decisive role, the accumulated occurrence area at the end of May was positively correlated with the year occurrence area, and the correlation coefficient was 0.8172. A prediction model of occurrence area based on the disease field rate at the end of boot stage and the disease plant rate at dough stage was established. Through investigation of the causes of blast disease in Nanchong in recent years, the authors found that the internal causes were the decrease or lose of rice blast resistance, the increase of high pathogenicity physiological race and the appearing of stronger pathogenicity physiological race, the external causes were suitable temperature and more overcast and rain and less sunlight. During 2010 to 2015, the short-term forecast accuracy for rice blast in Nanchong was up to 100%, medium-term and long-term forecast accuracy were also up to 98% and 95%, respectively, which increased by 5-15 percentage points compared with that before 1997, thus increased the control effect of rice blast in Nanchong by15%-30%.

关键词

水稻;稻瘟病;区划;发生特点;流行规律;预警水平

Key words

rice; rice blast; regionalization; occurrence characteristics; epidemic regularity; early warning level

引用本文

导出引用
彭昌家,白体坤,丁攀,冯礼斌,郭建全,杨宇衡,尹怀中,龙维国,文旭,肖立,崔德敏,苟建华,李鸿韬,何海燕,郑艳. 南充市水稻稻瘟病区划和发生流行规律研究. 中国农学通报. 2016, 32(6): 182-192 https://doi.org/10.11924/j.issn.1000-6850.casb15090025
Peng Changjia,Bai Tikun,Ding Pan,Feng Libin,Guo Jianquan,Yang Yuheng,Yin Huaizhong,Long Weiguo,Wen Xu,Xiao Li,Cui Demin,Gou Jianhua,Li Hongtao,He Haiyan and Zheng Yan. Regionalization and Epidemic Regularity of Rice Blast in Nanchong City. Chinese Agricultural Science Bulletin. 2016, 32(6): 182-192 https://doi.org/10.11924/j.issn.1000-6850.casb15090025

0 引言

土壤水是连接大气水、地表水和地下水的重要枢纽,在四水转化过程中发挥着至关重要的作用,精确地刻画和描述土壤水运动状况有益于更深刻地理解四水转化内涵。土壤水作为载体,还承担着碳氮磷钾等营养元素的运移转化。此外,与土壤含水量密切相关的诸如盐碱化、盐渍化、植被退化、农业干旱、生态系统退化、生态环境恶化等地质和生态环境问题很大程度限制了农业、生态系统和社会经济的健康发展[1,2,3,4],受到越来越多的关注。因此,准确地获取土壤水分的时空分布信息,精确刻画土壤水分运动过程和规律,针对性提出管理措施对新时代生态文明建设背景下山水林田湖草生态系统的健康尤为重要。本研究将介绍不同尺度土壤水分监测的主要研究手段,并就此研究方向做出进一步展望。

1 田间、场地等中小尺度

1.1 取土烘干法

取土烘干法是测定土壤含水率最直接的方法,应用最广,也最为准确和直观,一般用于对其他监测方法的标定。该方法利用烘箱等仪器将取得的土壤烘干至重量不再变化,将烘干损失的水量除以烘干前的重量获得土壤的质量含水率,进一步测量该样品的土壤容重就可以获得体积含水率。其工作量大,需要消耗大量的时间和精力,获得的数据具有延时性,难以长时间连续地观测,且取样的钻孔会对土壤的结构进行破坏,影响作物的根系吸水和土壤水分运动,进一步影响后续的观测结果的精度。

1.2 介电常数法

介电常数表征的是某种介质传递电磁波或者脉冲的能力,土和水的介电常数相差较大,土壤水分含量的细微变化将可以引起土-水介质介电常数的剧烈变化,因此可以通过测定土-水介质的介电常数得到土壤含水率。目前应用较广的有频域反射法系统(FDR)和时域反射法系统(TDR)。该方法布设和操作简单,且可以在线获得连续精确的土壤水分监测数据,保持土体相对的完整性,但测量结果容易受到土壤孔隙率、温度和容重的影响。研究者也相继提出了修正方法使FDR系统可以精确地测量不同温度下的土壤含水量,例如高磊等[5]引入温度修正系数来提高FDR系统在5~60℃条件下的观测精度;李元寿等[6]将FDR系统利用烘干法标定以后应用于测量高寒草地的土壤含水量。此外,研究者基于此技术引入了分频技术,实现了农田土壤水盐的连续监测,获得了精确的基础水盐监测数据[7]

1.3 中子仪法

中子仪通过测量快中子与土壤水分中氢原子碰撞而转化为慢中子的数量来感知土壤水分状况,通过标定可以将计数器的计数转换为土壤水分含量[8]。该方法操作简单,只需将探头安放在预先设置好的观测点便可以获得连续的监测数据,上下移动探头便可以获得不同深度的土壤含水率,且可以与计算机进行联网,实现数据的实时获得,受外界扰动的影响较小,可以长时间地观测。
值得注意的是,建立中子数与土壤含水率的回归方程进行仪器的标定是该方法的一大难题,陈明松和邵明安[9]对中子仪分层标定,实现了坡面土壤水分含水率的监测;田昌玉等[10]针对不同土壤层次、土壤质地、土壤有机质含量和土壤剖面含水量变化等条件提出了标定曲线的绘制方法,为该方法的推广提供了数据基础;刘巧玲等[11]开展中子水分仪测试土壤水分的标定试验,系统分析了中子仪计数时间、中子管安装稳定时间以及降雨、灌水等提高土壤含水率的方法对标定结果的影响。该方法的弊端在于费用较高,且存在辐射,可能会对试验人员的健康造成危险,需要在安装和取样过程中注意安全。

1.4 示踪剂法

当前用来指示表征土壤水分运移转化的示踪剂主要有氢氧等稳定同位素示踪剂、温度示踪剂、化学示踪剂(如溴示踪剂)。同一元素、不同同位素原子(或分子)之间由于其质量或自旋等核性质不同而造成的物理和化学性质相异,使不同时间和空间范围内土壤水氢氧稳定同位素的丰度有很大的差异[12]。气候条件、植被类型、下垫面条件、包气带岩性和土壤水补给来源特征等因子也极大地影响着土壤水氢氧稳定同位素的分布特征[13]。同位素分馏效应的不同使得其可以用来指示水循环的路径、揭示各种水体的成因、赋存条件和水循环机理等[14,15]。该方法应用简单,且不存在放射性危险,目前已广泛应用于研究植物蒸散发的分割、土壤入渗规律、植物需水等水循环过程的溯源和分析:例如马斌等[16]和刘添文[17]利用氢氧稳定同位素分别研究了华北平原石家庄包气带土壤水的入渗规律及年补给量和孝感市土壤水的入渗补给;张翔等[13]借助氢氧稳定同位素对鄱阳湖湿地中土壤水的组成进行了溯源;Jean等[18]利用氧同位素对植物蒸腾和土壤蒸发进行了分割;吴友杰和杜太生[19]利用氧同位素研究分析了土壤水、蒸发水汽、蒸腾水汽和大气背景混合水汽氧同位素组成分布特征,并对玉米农田蒸散发进行估算和区分。
温度作为示踪剂的理论基础是地质体中受各种环境因素影响的对流传热问题[20]。土壤水分运动过程中会传递能量,同时温度场的差异也会反过来影响土壤水分的运动。温度示踪的研究也从早期的定性分析发展到基于渗流场与温度场耦合模型的定量分析,其应用范围也越来越广,地质环境条件和边界条件也更趋复杂[21,22,23]。霍思远等[24]运用温度示踪法探寻了土壤水入渗补给速率,并发现浅层土壤含水率波动性是影响温度示踪精度的主要因子。潘维艳等[25]利用温度示踪法研究了干旱-半干旱地区农田渠道的入渗规律。该方法易于操作,且可以实现原位观测,即时地对数据进行采集和分析,但如何削减表层温度场的剧烈变化给观测带来的误差仍需要进一步地研究。
化学示踪剂投放简单,经济实用,可以定时、定位、定量观测,观测示踪剂峰值的变化,就可以计算水分通量,且无需测定水文地质参数,已经广泛应用于水循环研究过程中。例如谭畅等[26]借助KBr示踪剂,研究了干旱地区的潜水蒸发规律,为潜水蒸发的研究提供了新的思路。也有研究表明,应用该示踪剂时需要排除降雨对试验结果的影响,因此其更适用于干旱-半干旱地区[27]

1.5 地球物理探测技术反演法

地球物理探测技术是一种间接观测方法,其通过建立物理参数与土壤含水率的经验关系或物理模型来实现对土壤含水率的反演,常用的土壤物理参数有电导率和电磁波反射率等。目前常用的技术方法主要有高密度电阻成像法、探地雷达、电阻层析成像技术、核磁共振技术[28]
地球物理探测技术测量较快,而且可应用于大尺度监测,但投入代价高,图像的解译难度较高,而且解析出的土壤含水率的精度很大程度依赖于建立的土壤含水率与物探观测参数的经验关系和物理模型,给监测的结果带来了很大的误差和不确定性,其应用前景仍有待进一步地研究。

2 区域流域等大尺度

区域尺度的土壤水分监测影响因素众多,且具有很大的时空变异性,传统监测很难获得较大区域的土壤水分信息,遥感技术的引入给大尺度土壤水分监测提供了新的思路,研究者借助可见光-近红外遥感、热红外遥感、微波遥感等技术方法获得空间大尺度遥感数据,建立遥感数据与土壤水分信息之间的统计或经验模型,便可以实现土壤水分信息的间接观测。遥感直接监测的数据指标主要包括反射率、地表温度、土壤介电特性等。

2.1 反射率法

研究者发现特定土壤类型对可见光-近红外波段的反射率随土壤水分增加而减小[29],基于此现象进一步提出了很多经验模型,如式(1)所示[30]
R=aebP
(1)
式中R是光谱反射率,P是土壤含水率,ab是与土壤质地和类型相关的待给定的参数。基于该模型,刘培君等引入了“光学植被盖度”参数,并进一步借助土壤水分光谱信息,构建了土壤水分遥感AVHRR模型,并在阜北农场进行了验证,与实际情况较为符合,给大面积土壤水分监测提供了新的方法[31]。然而,该模型缺乏坚实的物理基础,需要进行大量的前置工作来获取不同土壤类型准确的参数值,耗时较长。
有植被覆盖的地区,植被生长状况与土壤含水率值呈现密切的正相关,研究者通过对红和近红外波段的反射信息的线性或非线性组合,进一步提出了植被指数的概念,并将其广泛应用于评估农田水分的干旱状况,目前应用最广的是归一化植被指数NDVI (normalized difference vegetation index),计算公式如式(2)所示。
NDVI=(Rnr-Rr)/Rnr+Rr
(2)
式中RnrRr分别是近红外和红波段的反射率。NDVI反映的是特定年份的土壤干旱状况的绝对值,其与区域的气候状况、人类活动和管理水平密切相关,指导管理措施制定时,还需要根据地域气候特征进行调整。因此,研究者进一步引入了NDVI多年平均值和极值作为参照物,进一步提出了距平植被指数AVI (anomoly vegetation index)和植被条件指数VCI (vegetation condition index)等参数[32],借助上述参数就可以量化特定年份的干旱状况在较长时间序列中所处的水平和特定区域干旱状况的空间分布,例如吕潇然利用VCI评估了云南省2004—2013年的干旱状况,获得了干旱易发的时段和地区[33],Alamdarloo等利用NDVIVCI评估了极度干旱区、干旱区、半干旱区和半湿润区4个气候条件下干旱对植被生态脆弱性的影响[34]。相比于NDVIAVIVCI作为干旱评价指标更能客观地表征某一地区的干旱状况,更方便地指导生产实践和制定干旱应对措施[35,36,37]。上述指数也存在一定的局限性,只能应用到植被旺盛的地区,应用到植被较为稀疏的地区或裸土时,精度会大幅下降。

2.2 反射率-温度协同反演法

土壤的水热传输密切相关,土壤含水率的增加会引起土壤比热容增加,改变土壤表层温度的变化幅度,因此温度的空间分布也一定程度上反映了土壤含水率的空间分布。然而,土壤温度和含水率的影响因素众多,难以确定普适性的关系来刻画土壤温度与含水率的关系。研究者经过系统的研究,分析了热惯量遥感成像的原理,提出了表观热惯量(ATI)的概念,并进一步忽略地理纬度的影响,获得了表观热惯量的简化计算方法[38,39],如式(2)所示。
ATI=(1-P)/Tmax-Tmin
(3)
式中,P是反射率,TmaxTmin分别是昼夜的最高和最低温度。例如吴黎借助ATI反演了河北省3个县的土壤含水率并取得了良好的结果,并进一步研究发现了应用ATI反演土壤水分的NDVI阈值为0.35,NDVI大于0.35时,应用该方法会使得反演精度大幅下降[40],杨树聪也发现了类似结果[41]
另外,研究者经过对地表温度和植被指数分布的散点图进行分析,发现地表温度和不同的植被指数之间存在不同的特征空间,其中应用较广的有Price提出的地表温度和NDVI的三角形特征空间和Moran引入空气温度以后提出的梯形空间[42]。基于此,研究者提出了条件温度植被干旱指数(TVDI),并广泛应用于干旱的预测。例如齐述华等[43]利用TVDI预测了全国的旱情状况,并借助实测的气象站0~10 cm土壤含水率数据对预测结果进行了验证,结果表明利用基于地表温度和植被指数的旱情参数评价指标优于基于植被指数的评价指标;王纯枝等[44]利用TVDI反演了黄淮海平原的土壤湿度,也获得了合理的结果。值得注意的是该方法中地表温度与植被指数的特征空间的构建和模型参数的获取是一大难点,需要大量的精确数据来支撑。

2.3 介电特性法

土壤和水各自介电常数的巨大差异是构建介电常数与土壤含水率数学关系的物理基础,且试验结果表明,对于特定土壤,其介电常数的实部由土壤含水率决定[45],而介电常数难以通过遥感获得。研究者通过对微波遥感器发射的信号和地表反射的微波信号进行处理,提出了后向散射系数的概念,而后向散射系数与介电常数密切相关,通过试验数据分析后向散射系数和土壤含水率的关系,就可以反演得到土壤含水率,其中线性关系最为简单实用。Ulaby等[46]研究发现线性关系不宜应用于过干或者过湿的土壤。Ulaby等[47]同时提出了土壤含水率在田间持水率50%~150%时去除植被干扰项的方法。随后研究者提出了土壤水分反演Dobson模型和Hallikainen简化实部模型,且都通过了实地观测数据的验证,成为目前应用很广的土壤水分遥感反演模型[48]。例如王学等[49]在河套灌区解放闸灌域分别构建了Dobson半经验模型、Hallikainen简化实部经验模型和统计回归模型,与土壤墒情的实测结果对比表明,上述模型表现良好,且可以应用于盐渍土[49]
介电常数的获得一般采用微波遥感,其物理机制明确,且微波可以全天时无间断地发射,穿透力强,不受云层的干扰。但是后向散射系数对诸多影响因素都比较敏感,需要排除诸如粗糙度,植被覆盖等的影响[50]
遥感技术给大尺度土壤水的观测提供了便利的方法,对难以获得数据的地区,提出了干旱系数的替代方法来间接表征土壤含水率状况,其精确值的获得仍是一大难题。此外,现在的观测精度较高的结果深度局限于地表0~20 cm,最深不超过40 cm。土壤表层含水率与深层含水率之间的关系影响因素众多,难以通过构建简单的数学模型来刻画表层含水率与深层含水率之间的关系,已有研究者将微波遥感和土壤水分运动模型相结合,获得了根区土壤较大深度的含水率[51,52]。此外,遥感数据和水文模型的尺度匹配问题仍然有待进一步研究。

3 结论和展望

本研究梳理了场地等中小尺度和区域流域等大尺度常用的土壤水分监测的技术手段和方法,并就其优缺点和适用范围进行了简单的介绍,研究者可根据自己的实际情况选用合适的方法。然而,现有的高精度观测技术手段包括无人机遥感大多局限于中小尺度。卫星遥感和航空遥感等遥感技术方法是大尺度观测的主要手段,随着遥感卫星可获取数据的增加,如何有效整合多源遥感数据,更好地服务于土壤水分反演也是一大挑战。此外,遥感技术大多依赖实地观测高精度数据的验证,且利用遥感信息反演土壤水分所需时间较长,存在一定的滞后性,其观测精度也随土壤深度的增加逐渐降低。因此,如何在保持精度的前提下将中小尺度观测手段、数据产品扩展到大尺度或探寻遥感信息与土壤水分相互作用的物理机制并最终实现大尺度土壤水分高精度低延时性观测仍有待进一步研究。
另一方面,伴随计算机计算能力高速提升,土壤水分运动数值模型的快速发展也给土壤水分观测带来更大的便利,其以较低的成本和耗时就可以获得现行和变化环境下土壤水分的响应规律,将实地观测与数值模拟有机结合,以高标准高精度的观测数据校核数值模型,利用数值模型分析人类活动和气候变化对土壤水分的影响,反过来指导土壤水分观测在关键区、生态脆弱区和重大战略区的应用和实践,发挥各自优势,更有效地制定土壤水分管理方案。
此外,土壤水分变化诱发的生态环境、地质环境问题频发,亟需研究其发生机理、揭示其物理机制和内在联系。同步实现土壤水分观测与生态环境地质问题预警也是未来一大趋势。

参考文献

[1] Zhang Q F. Strategies for developing green super rice[J].Proceedings of the national Academy of Sciences,2007,104(42):16402-16409.
[2] 彭昌家,丁攀,白体坤,等.生物农药和生化复配制剂预防水稻穗颈瘟的效果[J].贵州农业科学,2015,43(8):95-98.
[3] 彭昌家,白体坤,冯礼斌,等.南充市水稻稻瘟病综合防控技术探讨[J].中国稻米,2015,21(3):62-67.
[4] 何秀英,廖耀平,陈钊明,等.水稻稻瘟病抗病育种研究进展与展望[J].广东农业科学,2011,38(1):30-33.
[5] 刘国权,孟昭河,任艳军,等.水稻抗稻瘟病研究进展与对策[J].中国农学通报,2004,20(1):211-214.
[6] 何秀英,刘新琼,王丽,等.稻瘟病新隐性抗病基因pi55(t)的遗传及定位[J].中国科学:生命科学,2012,42(2):125-134.
[7] 李彬,邓元宝,颜学海,等.一个粳稻来源抗稻瘟病基因的鉴定、遗传分析和基因定位[J].作物学报,2014,40(1):54-62.
[8] 王玲,左示敏,张亚芳,等.四川省稻瘟病菌群体遗传结构分析[J].中国水稻科学,2015,29(3):327-334.
[9] 何秀英,王玲,吴伟怀,等.水稻稻瘟病抗性基因的定位、克隆及育种应用研究进展[J].中国农学通报,2014,30(6):1-12.
[10] 孙漱沅,孙国昌.我国稻瘟病研究的现状和展望[J].植保技术与推广,1996,16(3):39-40.
[11] 温小红,谢明杰,姜健,等.水稻稻瘟病防治方法研究进展[J].中国农学通报,2013,29(3):190-195.
[12] 彭昌家,白体坤,冯礼斌,等.南充市水稻稻瘟病综合防控技术研究[J].中国农学通报,2015,31(11):190-199.
[13] 彭昌家.南充市水稻稻瘟病重发原因及综合防治研究[J].中国农学通报,2002,18(1):78-81.
[14] 何会流,毛建辉,卢代华,等.四川水稻品种对稻瘟病和稻曲病抗性评价[J].西南大学学报:自然科学版,2008,30(7):104-109.
[15] 杜晓宇,丁攀,陈晓娟,等.南充市水稻稻瘟病重发成因及治理对策[J].中国农学通报,2009,25(5):218-222.
[16] 殷世才,杨和生,保世森,等.陆良稻瘟病发生规律及预测预报模型应用[J].云南农业大学学报:自然科学版,1995,10(2):172.
[17] 金星,谈孝凤,杨光灿,等.贵州省稻瘟病发生规律研究[J].耕作与栽培,2000,(3):48-49,57.
[18] 宋成艳.寒地水稻稻瘟病发生规律研究[D].哈尔滨:东北农业大学,2007.
[19] 夏敬源.公共植保、绿色植保的发展与展望[J].中国植保导刊,2010,30(1):5-9.
[20] 张跃进.农作物有害生物测报技术手册[M].北京:中国农业出版社,2006.
[21] 中华人民共和国国家质量监督检验检疫总局,中国国家标准化管理委员会.中华人民共和国国家标准,GB 15790-2009-T稻瘟病测报调查规范[EB/OL].(2009-03-27)[2015-01-26].http://www.doc88.com/p-282362196466.html.
[22] 柳武革,王丰,刘振荣,等.利用分子标记技术聚合Pi-1和Pi-2基因改良三系不育系荣丰A的稻瘟病抗性[J].分子植物育种,2012,10(5):575-582.
[23] 彭昌家,丁攀,白体坤,等.南充市主要粮油作物重大病虫害发生流行及灾变规律研究[J].安徽农学通报,2014,20(13):78-83.
[24] 刘彬,叶慧丽,姚琳,等.不同类型水稻鉴别品种对四川稻瘟病菌生理小种的鉴定与评价研究[J].西南农业学报,2007,20(3):400-407.
[25] 张雪梅,冯慧,白玉莲,等.水稻抗稻瘟病基因和杂交稻品种在四川的抗性评价[J].西南农业学报,2012,25(4):1266-1272.
[26] 中国农业科学研究院植物保护研究所,中国植物保护学会.中国农作物病虫害(第三版,上册)[M].北京:中国农业出版社,2015:1-9.
[27] 四川省植农牧厅物保护站,四川省植物保护学会.植保专业队员手册[M].成都:四川科技出版社,1992:17-22.
[28] 农业部全国植物保护总站.植物医生手册[M].北京:化学工业出版社,1994:20-24.
PDF(2789 KB)

文章所在专题

农业气象

Accesses

Citation

Detail

段落导航
相关文章

/