江苏省农村牲畜棚三带喙库蚊发生消长的中长期影响初探*

杨维芳 褚宏亮 徐 燕 张育富 刘 慧 刘大鹏 吴治明

(江苏省疾病预防控制中心，南京 210009)

三带喙库蚊Culextritaeniorhynchus属于双翅目(Diptera)蚊科(Culicidae)，是我国大多数地区流行性乙型脑炎(简称“乙脑”)的主要传播媒介，也是东南亚和远东地区包括越南、印度、柬埔寨、日本、韩国、朝鲜等国乙脑的主要媒介。该蚊幼虫主要孳生于稻田、沼泽、池塘等城乡清洁或稍污染的大中型水体中(李春晓，2007；褚宏亮等，2020)；成蚊野栖，主要吸食猪、牛等动物血液，兼吸人血(钱万红等，2008；Lytraetal.，2014)。三带喙库蚊的发生消长与乙脑发病密切相关，其密度高峰出现1个月后为乙脑发病高峰，而该蚊消失后1个月，不再有乙脑病例发生(褚宏亮，2017；张育富等，2021)。我国蚊虫监测的5类生境中，牲畜棚是三带喙库蚊活动栖息的绝对优势生境(刘小波等，2014；胡雅劼等，2017；杨维芳等，2010，2021；郭玉红等，2019；赵宁等，2019)，该生境三带喙库蚊的发生消长对乙脑控制有重要意义(张洪江等，2016)。为了解江苏省农村牲畜棚三带喙库蚊发生消长及其前期密度对后期密度的影响，以有效控制三带喙库蚊的叮刺骚扰和传病危害，特开展了相关研究。

1 材料与方法

1.1 数据来源

数据来自江苏省2009—2021年省级病媒生物监测网络中省级监测点报送的牲畜棚生境中三带喙库蚊监测数据。江苏省下辖13个设区市、21个县级市和19个县共计53个监测点全部纳入省级病媒生物监测网络，其中苏南15个监测点、苏中13个监测点、苏北25个监测点。

1.2 监测方法

选用武汉市吉星环保科技有限责任公司生产的功夫小帅诱蚊灯(波长2537Å，功率8W)，采用诱蚊灯法开展成蚊监测。在蚊虫活动季节，设区市监测点选4个农村牲畜棚生境，市(县)监测点选2个农村牲畜棚生境，日落前1 h至次日晨连续监测整夜(12 h)，每月监测2次，分类鉴定、计数三带喙库蚊。

1.3 统计学分析

1.3.1密度计算 设区市、市(县)监测点，苏南、苏中、苏北，以及全省农村牲畜棚三带喙库蚊密度，为行政区域(省、苏南、苏中、苏北、设区市、市、县)汇总密度，计算公式如下，

密度[只/(灯·夜)]=

1.3.2消长分析 年际消长与季节消长，在excel 2016中进行。

1.3.3贝叶斯相关Bayesian correlation分析 在统计分析软件JASP 0.11.1.0中，采用贝叶斯相关矩阵法Bayesian correlation matrix，分析全省、苏南、苏中、苏北农村牲畜棚三带喙库蚊年密度的相关性。JASP是一款免费、开源的统计软件，使用 R 语言工具包进行数据处理,但使用者不需要安装 R，就能轻松地进行描述统计、t检验、方差分析、相关与回归分析、因子分析等。

1.3.4贝叶斯配对样本t检验(Bayesian paired samplest-test) 在统计分析软件JASP 0.11.1.0中，采用贝叶斯配对样本t检验，分析苏南、苏中和苏北农村牲畜棚三带喙库蚊年密度是否存在差异。

1.3.5多元线性逐步回归分析 在统计分析软件JASP 0.11.1.0中，根据2009—2021年的全省年密度和年度月密度，13个设区市本级年密度，苏南、苏中和苏北的年密度，设定因变量Yn、自变量Xn，Yn和Xn均为年密度，其中n为年份；采用多元线性逐步回归法，分析2009—2021年期间，前期密度对后期密度的影响趋势。

2 结果

2.1 牲畜棚三带喙库蚊发生消长

2.1.1年际消长 2009—2021年全省农村牲畜棚三带喙库蚊密度的年际消长总体呈下降趋势；年密度在12.11～243.12只/(灯·夜)波动，其中2009、2010、2015年3个年份的年密度≥60只/(灯·夜)，2016和2021年的年密度<24只/(灯·夜)；2010、2015、2017年3个密度波峰，随时间波峰密度显著下降；2011—2013、2016年处于明显的密度波谷时期，且位于年际消长趋势线下方(图1)。

图1 2009—2021年江苏省农村牲畜棚三带喙库蚊密度的年际消长Fig.1 Annual density dynamic of Cx. tritaeniorhynchus in rural livestock sheds of Jiangsu province 2009-2021

2.1.2苏南、苏中、苏北的年际消长 2009—2021年，苏南、苏中、苏北农村牲畜棚的年际消长总体均呈下降趋势；苏南的年密度整体高于苏北、苏中(图2)。苏南、苏中、苏北和全省年密度的贝叶斯相关矩阵法分析显示，苏南的年密度占全省年密度的权重大，两者的相关程度最密切(肯德尔等级相关系数Kendall′s tau-b =0.821，贝叶斯因子BF10=307.988)，苏中、苏北与年密度的相关程度不密切(Kendall′s tau-b分别为0.026、0.205；BF10<10，分别为0.351、0.540)。苏南、苏中、苏北年密度的贝叶斯配对样本t检验显示，苏南与苏中年密度差异显著(BF10=26.261)，苏南与苏北、苏中与苏北年密度差异不显著(BF10<10)。

图2 2009—2021年江苏省苏南、苏中和苏北农村牲畜棚三带喙库蚊密度的年际消长Fig.2 Annual density dynamic of Cx. tritaeniorhynchus in rural livestock sheds of Southern,Central and Northern Jiangsu from 2009 to 2021

2.1.3季节消长 2009—2021年，全省农村牲畜棚三带喙库蚊密度的季节消长总体呈单峰型(图3)，3月份偶见三带喙库蚊活动，6月份开始活动明显增强，7月份为其活动高峰，高峰期平均密度为313.45只/(灯·夜)；随后密度逐步显著下降，但8、9月份的密度仍处于较高水平，其中8月份密度为176.94只/(灯·夜)[>60只/(灯·夜)]；10月份及后期密度均较低。密度季节消长还显示，平均密度越高的月份，其年际间该月的密度波动幅度越大，如全省各年7月份的密度为36.99～1 470.85只/(灯·夜)。

图3 2009—2021年江苏省农村牲畜棚三带喙库蚊密度的季节消长Fig.3 Seasonal density dynamic of Cx. tritaeniorhynchus in rural livestock sheds of Jiangsu province 2009-2021

2.2 前期密度对后期密度的影响

2009—2021年，江苏省全省农村牲畜棚的年密度和年度月密度，13个设区市本级年密度，苏南、苏中和苏北的年密度，通过多元线性逐步回归分析显示(表1)，13个年份中，2010—2013、2015—2016、2019—2021年9个年份的密度与该年份前5年内的1～2个年份的密度呈显著线性相关；2014、2017—2018年共3个年份的密度与该年份前5年内及前5～10年的多个年份的密度呈显著线性相关。

表1 2009—2021年江苏省农村牲畜棚三带喙库蚊前期密度对后期密度的影响Tab.1 Impact of early density of Cx. tritaeniorhynchus on later density in rural livestock sheds of Jiangsu province 2009-2021

3 讨论

3.1 农村三带喙库蚊的防制成效

江苏省农村牲畜棚2009—2021年共计13年的三带喙库蚊密度整体呈下降趋势，这与江苏省农村地区近20年来卫生城镇(村)建设，以及“金山银山就是绿水青山”的生态环境保护等措施密不可分。2005 — 2018 年江苏省共报告乙脑病例 495 例，年均发病率为 0.045/10 万，乙脑发病率总体呈下降趋势(χ2= - 16.11，P<0.01)(张育富等，2021)，显示江苏省农村三带喙库蚊的防制取得重大成效。

2009—2021年间，随着GB/T 27771—2011《病媒生物密度控制水平 蚊虫》发布，卫生城镇(村)创建工作对蚊虫防制提出了更高、更具体的要求，江苏省在防蚊灭蚊措施强化、蚊虫生境科学管理方面，做了大量工作，蚊虫密度出现显著下降，形成了2011—2013年期间的第1个农村牲畜棚三带喙库蚊的年密度波谷；在此期间，江苏省爱卫会也命名了一大批的“江苏省卫生镇(乡、街道)”，如“苏爱卫[2011]第3号”中命名了22个“江苏省卫生镇(乡、街道)”、380个“江苏省卫生村”，苏爱卫〔2013〕11号中重新确认命名了1个“江苏省卫生县城”、19个“江苏省卫生乡镇”和390个“江苏省卫生村”。2016年是江苏省农村牲畜棚三带喙库蚊13年中的第2个年密度波谷，可能跟2016年江苏省畜禽养殖污染治理专项工作的大力实施，规模化养殖的关闭、搬迁，以及农户畜禽养殖大规模减少有一定关系。

3.2 三带喙库蚊密度存在显著的区域差异

苏南农村牲畜棚三带喙库蚊年密度与苏中差异显著，苏南与苏北、苏中与苏北差异不显著，同时苏南整体高过苏北、苏中。这种情况可能与苏南地区河湖水网密布、雨水充沛、农村水田比例较高等，适于三带喙库蚊幼虫孳生，同时苏南夏季高温持续时间普遍多于苏北、苏中，更适合三带喙库蚊生存等因素有关。牲畜棚周边的地面遥感信息情况，如水面情况，特别是稻田情况，对区域蚊虫密度的影响已被证实(王晓东等，2008；何隆华等，2014；杨维芳等，2015)。

3.3 三带喙库蚊前期密度对后期密度的影响

江苏省农村牲畜棚三带喙库蚊连续13年的年密度显示，多数年份的年密度仅对后期年密度有短期(<5年)影响，少数年份的年密度对后期年密度有中期(5～10年)影响，未见有长期(≥10年)影响的情况。由此可见，三带喙库蚊早期密度多对后期密度起中、短期影响，很难有长期影响，这也说明，某个年份的低密度，并不能保证后续年份的持续低密度，只有从控制三带喙库蚊孳生地、活动栖息地等治本清源的角度出发，才能真正控制三带喙库蚊叮刺骚扰、传病危害。

3.4 农村三带喙库蚊的防制应常抓不懈

牲畜棚中三带喙库蚊密度总体下降，显示该蚊种叮刺血和传播疾病的平均风险总体降低，但部分牲畜棚三带喙库蚊高峰期7月份的密度可高达>1 000只/(灯·夜)，8月份密度≥60只/(灯·夜)的牲畜棚也非常常见，同时江苏省在三带喙库蚊携带乙脑病原的监测中，也发现多地乙脑病原体阳性比例高的情况(Chuetal.，2017)，因此，在三带喙库蚊密度高峰季节、部分三带喙库蚊密度高或乙脑病原体阳性比例高的牲畜棚及地区，仍要高度重视农村，特别是牲畜棚三带喙库蚊的防制工作。