江蘇省農(nóng)村牲畜棚三帶喙庫蚊發(fā)生消長的中長期影響初探*

楊維芳 褚宏亮 徐 燕 張育富 劉 慧 劉大鵬 吳治明

(江蘇省疾病預(yù)防控制中心，南京 210009)

三帶喙庫蚊Culextritaeniorhynchus屬于雙翅目(Diptera)蚊科(Culicidae)，是我國大多數(shù)地區(qū)流行性乙型腦炎(簡稱“乙腦”)的主要傳播媒介，也是東南亞和遠(yuǎn)東地區(qū)包括越南、印度、柬埔寨、日本、韓國、朝鮮等國乙腦的主要媒介。該蚊幼蟲主要孳生于稻田、沼澤、池塘等城鄉(xiāng)清潔或稍污染的大中型水體中(李春曉，2007；褚宏亮等，2020)；成蚊野棲，主要吸食豬、牛等動(dòng)物血液，兼吸人血(錢萬紅等，2008；Lytraetal.，2014)。三帶喙庫蚊的發(fā)生消長與乙腦發(fā)病密切相關(guān)，其密度高峰出現(xiàn)1個(gè)月后為乙腦發(fā)病高峰，而該蚊消失后1個(gè)月，不再有乙腦病例發(fā)生(褚宏亮，2017；張育富等，2021)。我國蚊蟲監(jiān)測(cè)的5類生境中，牲畜棚是三帶喙庫蚊活動(dòng)棲息的絕對(duì)優(yōu)勢(shì)生境(劉小波等，2014；胡雅劼等，2017；楊維芳等，2010，2021；郭玉紅等，2019；趙寧等，2019)，該生境三帶喙庫蚊的發(fā)生消長對(duì)乙腦控制有重要意義(張洪江等，2016)。為了解江蘇省農(nóng)村牲畜棚三帶喙庫蚊發(fā)生消長及其前期密度對(duì)后期密度的影響，以有效控制三帶喙庫蚊的叮刺騷擾和傳病危害，特開展了相關(guān)研究。

1 材料與方法

1.1 數(shù)據(jù)來源

數(shù)據(jù)來自江蘇省2009—2021年省級(jí)病媒生物監(jiān)測(cè)網(wǎng)絡(luò)中省級(jí)監(jiān)測(cè)點(diǎn)報(bào)送的牲畜棚生境中三帶喙庫蚊監(jiān)測(cè)數(shù)據(jù)。江蘇省下轄13個(gè)設(shè)區(qū)市、21個(gè)縣級(jí)市和19個(gè)縣共計(jì)53個(gè)監(jiān)測(cè)點(diǎn)全部納入省級(jí)病媒生物監(jiān)測(cè)網(wǎng)絡(luò)，其中蘇南15個(gè)監(jiān)測(cè)點(diǎn)、蘇中13個(gè)監(jiān)測(cè)點(diǎn)、蘇北25個(gè)監(jiān)測(cè)點(diǎn)。

1.2 監(jiān)測(cè)方法

選用武漢市吉星環(huán)?？萍加邢挢?zé)任公司生產(chǎn)的功夫小帥誘蚊燈(波長2537?，功率8W)，采用誘蚊燈法開展成蚊監(jiān)測(cè)。在蚊蟲活動(dòng)季節(jié)，設(shè)區(qū)市監(jiān)測(cè)點(diǎn)選4個(gè)農(nóng)村牲畜棚生境，市(縣)監(jiān)測(cè)點(diǎn)選2個(gè)農(nóng)村牲畜棚生境，日落前1 h至次日晨連續(xù)監(jiān)測(cè)整夜(12 h)，每月監(jiān)測(cè)2次，分類鑒定、計(jì)數(shù)三帶喙庫蚊。

1.3 統(tǒng)計(jì)學(xué)分析

1.3.1密度計(jì)算 設(shè)區(qū)市、市(縣)監(jiān)測(cè)點(diǎn)，蘇南、蘇中、蘇北，以及全省農(nóng)村牲畜棚三帶喙庫蚊密度，為行政區(qū)域(省、蘇南、蘇中、蘇北、設(shè)區(qū)市、市、縣)匯總密度，計(jì)算公式如下，

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

1.3.2消長分析 年際消長與季節(jié)消長，在excel 2016中進(jìn)行。

1.3.3貝葉斯相關(guān)Bayesian correlation分析 在統(tǒng)計(jì)分析軟件JASP 0.11.1.0中，采用貝葉斯相關(guān)矩陣法Bayesian correlation matrix，分析全省、蘇南、蘇中、蘇北農(nóng)村牲畜棚三帶喙庫蚊年密度的相關(guān)性。JASP是一款免費(fèi)、開源的統(tǒng)計(jì)軟件，使用 R 語言工具包進(jìn)行數(shù)據(jù)處理,但使用者不需要安裝 R，就能輕松地進(jìn)行描述統(tǒng)計(jì)、t檢驗(yàn)、方差分析、相關(guān)與回歸分析、因子分析等。

1.3.4貝葉斯配對(duì)樣本t檢驗(yàn)(Bayesian paired samplest-test) 在統(tǒng)計(jì)分析軟件JASP 0.11.1.0中，采用貝葉斯配對(duì)樣本t檢驗(yàn)，分析蘇南、蘇中和蘇北農(nóng)村牲畜棚三帶喙庫蚊年密度是否存在差異。

1.3.5多元線性逐步回歸分析 在統(tǒng)計(jì)分析軟件JASP 0.11.1.0中，根據(jù)2009—2021年的全省年密度和年度月密度，13個(gè)設(shè)區(qū)市本級(jí)年密度，蘇南、蘇中和蘇北的年密度，設(shè)定因變量Yn、自變量Xn，Yn和Xn均為年密度，其中n為年份；采用多元線性逐步回歸法，分析2009—2021年期間，前期密度對(duì)后期密度的影響趨勢(shì)。

2 結(jié)果

2.1 牲畜棚三帶喙庫蚊發(fā)生消長

2.1.1年際消長 2009—2021年全省農(nóng)村牲畜棚三帶喙庫蚊密度的年際消長總體呈下降趨勢(shì)；年密度在12.11～243.12只/(燈·夜)波動(dòng)，其中2009、2010、2015年3個(gè)年份的年密度≥60只/(燈·夜)，2016和2021年的年密度<24只/(燈·夜)；2010、2015、2017年3個(gè)密度波峰，隨時(shí)間波峰密度顯著下降；2011—2013、2016年處于明顯的密度波谷時(shí)期，且位于年際消長趨勢(shì)線下方(圖1)。

圖1 2009—2021年江蘇省農(nóng)村牲畜棚三帶喙庫蚊密度的年際消長Fig.1 Annual density dynamic of Cx. tritaeniorhynchus in rural livestock sheds of Jiangsu province 2009-2021

2.1.2蘇南、蘇中、蘇北的年際消長 2009—2021年，蘇南、蘇中、蘇北農(nóng)村牲畜棚的年際消長總體均呈下降趨勢(shì)；蘇南的年密度整體高于蘇北、蘇中(圖2)。蘇南、蘇中、蘇北和全省年密度的貝葉斯相關(guān)矩陣法分析顯示，蘇南的年密度占全省年密度的權(quán)重大，兩者的相關(guān)程度最密切(肯德爾等級(jí)相關(guān)系數(shù)Kendall′s tau-b =0.821，貝葉斯因子BF10=307.988)，蘇中、蘇北與年密度的相關(guān)程度不密切(Kendall′s tau-b分別為0.026、0.205；BF10<10，分別為0.351、0.540)。蘇南、蘇中、蘇北年密度的貝葉斯配對(duì)樣本t檢驗(yàn)顯示，蘇南與蘇中年密度差異顯著(BF10=26.261)，蘇南與蘇北、蘇中與蘇北年密度差異不顯著(BF10<10)。

圖2 2009—2021年江蘇省蘇南、蘇中和蘇北農(nóng)村牲畜棚三帶喙庫蚊密度的年際消長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季節(jié)消長 2009—2021年，全省農(nóng)村牲畜棚三帶喙庫蚊密度的季節(jié)消長總體呈單峰型(圖3)，3月份偶見三帶喙庫蚊活動(dòng)，6月份開始活動(dòng)明顯增強(qiáng)，7月份為其活動(dòng)高峰，高峰期平均密度為313.45只/(燈·夜)；隨后密度逐步顯著下降，但8、9月份的密度仍處于較高水平，其中8月份密度為176.94只/(燈·夜)[>60只/(燈·夜)]；10月份及后期密度均較低。密度季節(jié)消長還顯示，平均密度越高的月份，其年際間該月的密度波動(dòng)幅度越大，如全省各年7月份的密度為36.99～1 470.85只/(燈·夜)。

圖3 2009—2021年江蘇省農(nóng)村牲畜棚三帶喙庫蚊密度的季節(jié)消長Fig.3 Seasonal density dynamic of Cx. tritaeniorhynchus in rural livestock sheds of Jiangsu province 2009-2021

2.2 前期密度對(duì)后期密度的影響

2009—2021年，江蘇省全省農(nóng)村牲畜棚的年密度和年度月密度，13個(gè)設(shè)區(qū)市本級(jí)年密度，蘇南、蘇中和蘇北的年密度，通過多元線性逐步回歸分析顯示(表1)，13個(gè)年份中，2010—2013、2015—2016、2019—2021年9個(gè)年份的密度與該年份前5年內(nèi)的1～2個(gè)年份的密度呈顯著線性相關(guān)；2014、2017—2018年共3個(gè)年份的密度與該年份前5年內(nèi)及前5～10年的多個(gè)年份的密度呈顯著線性相關(guān)。

表1 2009—2021年江蘇省農(nóng)村牲畜棚三帶喙庫蚊前期密度對(duì)后期密度的影響Tab.1 Impact of early density of Cx. tritaeniorhynchus on later density in rural livestock sheds of Jiangsu province 2009-2021

3 討論

3.1 農(nóng)村三帶喙庫蚊的防制成效

江蘇省農(nóng)村牲畜棚2009—2021年共計(jì)13年的三帶喙庫蚊密度整體呈下降趨勢(shì)，這與江蘇省農(nóng)村地區(qū)近20年來衛(wèi)生城鎮(zhèn)(村)建設(shè)，以及“金山銀山就是綠水青山”的生態(tài)環(huán)境保護(hù)等措施密不可分。2005 — 2018 年江蘇省共報(bào)告乙腦病例 495 例，年均發(fā)病率為 0.045/10 萬，乙腦發(fā)病率總體呈下降趨勢(shì)(χ2= - 16.11，P<0.01)(張育富等，2021)，顯示江蘇省農(nóng)村三帶喙庫蚊的防制取得重大成效。

2009—2021年間，隨著GB/T 27771—2011《病媒生物密度控制水平 蚊蟲》發(fā)布，衛(wèi)生城鎮(zhèn)(村)創(chuàng)建工作對(duì)蚊蟲防制提出了更高、更具體的要求，江蘇省在防蚊滅蚊措施強(qiáng)化、蚊蟲生境科學(xué)管理方面，做了大量工作，蚊蟲密度出現(xiàn)顯著下降，形成了2011—2013年期間的第1個(gè)農(nóng)村牲畜棚三帶喙庫蚊的年密度波谷；在此期間，江蘇省愛衛(wèi)會(huì)也命名了一大批的“江蘇省衛(wèi)生鎮(zhèn)(鄉(xiāng)、街道)”，如“蘇愛衛(wèi)[2011]第3號(hào)”中命名了22個(gè)“江蘇省衛(wèi)生鎮(zhèn)(鄉(xiāng)、街道)”、380個(gè)“江蘇省衛(wèi)生村”，蘇愛衛(wèi)〔2013〕11號(hào)中重新確認(rèn)命名了1個(gè)“江蘇省衛(wèi)生縣城”、19個(gè)“江蘇省衛(wèi)生鄉(xiāng)鎮(zhèn)”和390個(gè)“江蘇省衛(wèi)生村”。2016年是江蘇省農(nóng)村牲畜棚三帶喙庫蚊13年中的第2個(gè)年密度波谷，可能跟2016年江蘇省畜禽養(yǎng)殖污染治理專項(xiàng)工作的大力實(shí)施，規(guī)?；B(yǎng)殖的關(guān)閉、搬遷，以及農(nóng)戶畜禽養(yǎng)殖大規(guī)模減少有一定關(guān)系。

3.2 三帶喙庫蚊密度存在顯著的區(qū)域差異

蘇南農(nóng)村牲畜棚三帶喙庫蚊年密度與蘇中差異顯著，蘇南與蘇北、蘇中與蘇北差異不顯著，同時(shí)蘇南整體高過蘇北、蘇中。這種情況可能與蘇南地區(qū)河湖水網(wǎng)密布、雨水充沛、農(nóng)村水田比例較高等，適于三帶喙庫蚊幼蟲孳生，同時(shí)蘇南夏季高溫持續(xù)時(shí)間普遍多于蘇北、蘇中，更適合三帶喙庫蚊生存等因素有關(guān)。牲畜棚周邊的地面遙感信息情況，如水面情況，特別是稻田情況，對(duì)區(qū)域蚊蟲密度的影響已被證實(shí)(王曉東等，2008；何隆華等，2014；楊維芳等，2015)。

3.3 三帶喙庫蚊前期密度對(duì)后期密度的影響

江蘇省農(nóng)村牲畜棚三帶喙庫蚊連續(xù)13年的年密度顯示，多數(shù)年份的年密度僅對(duì)后期年密度有短期(<5年)影響，少數(shù)年份的年密度對(duì)后期年密度有中期(5～10年)影響，未見有長期(≥10年)影響的情況。由此可見，三帶喙庫蚊早期密度多對(duì)后期密度起中、短期影響，很難有長期影響，這也說明，某個(gè)年份的低密度，并不能保證后續(xù)年份的持續(xù)低密度，只有從控制三帶喙庫蚊孳生地、活動(dòng)棲息地等治本清源的角度出發(fā)，才能真正控制三帶喙庫蚊叮刺騷擾、傳病危害。

3.4 農(nóng)村三帶喙庫蚊的防制應(yīng)常抓不懈

牲畜棚中三帶喙庫蚊密度總體下降，顯示該蚊種叮刺血和傳播疾病的平均風(fēng)險(xiǎn)總體降低，但部分牲畜棚三帶喙庫蚊高峰期7月份的密度可高達(dá)>1 000只/(燈·夜)，8月份密度≥60只/(燈·夜)的牲畜棚也非常常見，同時(shí)江蘇省在三帶喙庫蚊攜帶乙腦病原的監(jiān)測(cè)中，也發(fā)現(xiàn)多地乙腦病原體陽性比例高的情況(Chuetal.，2017)，因此，在三帶喙庫蚊密度高峰季節(jié)、部分三帶喙庫蚊密度高或乙腦病原體陽性比例高的牲畜棚及地區(qū)，仍要高度重視農(nóng)村，特別是牲畜棚三帶喙庫蚊的防制工作。