曹玲,张庆瑞
(河南林业职业学院,河南洛阳 471002)
中华松干蚧(Matsucoccus sinensisChen),又名松针蚧,属半翅目(Hemiptera)松干蚧科(Matsucoccidae)松干蚧属(Matsucoccus),在河南省主要分布于三门峡、洛阳、南阳、信阳等地,主要危害油松(Pinus tabulaeformis)和马尾松(P.massoniana),不仅破坏自然景观,还严重影响树木的生长,制约了农村经济发展和林业生态建设,属河南省林业有害生物补充检疫对象,该虫在豫西地区的危害面积和危害程度呈逐年提高趋势。为确定中华松干蚧的发生对油松生长量的影响,笔者在2014—2018年,对豫西地区中华松干蚧的发生量和油松生长量进行了系统调查。
河南省三门峡市卢氏县东湾林场海拔1 035~1 585 m,属大陆性季风气候,年平均气温12.6℃,1月最冷,7月最热,无霜期平均184 d,年均日照时数2 118.0 h,年均降雨量647.8 mm,集中在7—9月。
在中华松干蚧最早传入地卢氏县官坡镇东湾林场草房林区及徐家湾乡设代表性标准地4个,标准地面积为0.667 hm2,每个标准地内设标准株不少于30株,标准地立地条件相同,标准株树龄相同。
从2014年开始,每年3月对标准地内样株进行每木检尺、标号,调查样株虫口密度,记录各项调查数据。并于次年同期作生长量调查,记录各标准株对应虫口密度的年度胸径生长量,并计算对应材积生长量(材积计算依据1998年7月出版的《河南省主要树种立木材积表》中的山区树种一元立木材积公式),至2018年度共记录4次,取得4组数据。将所有4组数据按虫口密度分级(每百束针叶50头为1级)进行统计。
根据调查样株胸径,计算立木材积,根据防治成本和防治效果,按照经济允许危害水平和经济阈值,分析确定若虫危害期防治指标。
山区树种一元立木材积计算公式为:
对于油松,式中:a=0.000 093 059,b=-0.315 3,c=0.983 17,e=1.992 682,f=0.450 60,g=0.621 517,h=0.005 174 8,i=0.671 166 1。
对表1和表2进行方差分析显示,各虫口密度对应生长量间(组间)有显著差异,其中表1经反正弦变换分析,F=39.223 4>F0.05(8,24)=2.36,进一步q检验得D=6.17,分析显示,当虫口密度达到250~300头/百束针叶时,胸径生长量与对照间有显著差异。
表1 油松胸径调查分析
2.2.1建立回归方程式。对表2中的系列虫口密度区间中值与年度材积平均增长量间进行回归分析,二者间有关系式(抛物线,二次回归):y=0.007 513 5+0.000 016 354x-0.000 000 088 679x2(r=-0.856 1)。
表2 油松材积分析
2.2.2确定防治指标。令Vz′=V0-Vz,为每公顷的损失量,是虫口密度Z的函数;Si为油松价值(元/m3);Ni为防治效果;Mi为油松每公顷的立木数;Ci为防治费用(元/hm2)。
根据Ci=Vz×Si×Ni×Mi和虫口密度(Z)与油松材积生长量(V)的关系模型(Vz=0.007 513 5+0.000 016 354Z-0.000 000 088 679Z2)得:
Ci=Vz′×Si×Ni×Mi=(V0-Vz)×Si×Ni×Mi
当Z=0时,V0=0.007 513 5,则对应虫口密度损失材积公式为:
现油松平均价值(Si)按500元/m3计,每公顷立木数(Mi)按1 000株/hm2计,防治效果(Ni)为80%,防治费用(Ci)为300元/hm2,则 有:300=(-0.000 016 354Z+0.000 000 088 679Z2)×500×0.8×1 000,即有:0.037 689Z2-6.950 45Z-300=0。解得油松中华松干蚧防治指标Z=220.5(头/百束针叶)(2龄若虫)。
根据中华松干蚧的发生发展规律、生活习性和虫口密度影响油松胸径与材积生长量的相关关系,确定2龄若虫(主要危害虫期)理论防治指标为220.5头/百束针叶,对生产实践具有指导意义。