孟慶禮，黃利亞

（1. 長白山自然保護(hù)管理中心，吉林 安圖 133613；2. 長白山科學(xué)研究院，吉林 延邊 133613）

口蹄疫是由小RNA 病毒科（Picornaviridae）口瘡病毒屬的口蹄疫病毒（Foot-and-mouth disease virus，F(xiàn)MDV）引起的急性傳染病，可感染包括牛、鹿、豬、山羊和羚羊在內(nèi)的多種偶蹄類動(dòng)物，對(duì)野生動(dòng)物種群的生態(tài)健康造成了一定威脅。由于口蹄疫是一種高度接觸性傳染病，在物種間傳播速度快、傳播途徑多、并且具有較強(qiáng)的變異性，導(dǎo)致其一度在全球范圍廣泛流行。與東北虎豹國家公園東南部相毗鄰的朝鮮，曾在全國范圍內(nèi)暴發(fā)口蹄疫疫情，因此東北虎豹國家公園對(duì)口蹄疫的風(fēng)險(xiǎn)評(píng)估與防控極為必要。由于國家公園內(nèi)自然環(huán)境和氣候條件多樣，景觀多樣性較高，因此如何實(shí)現(xiàn)口蹄疫的精準(zhǔn)防控，已經(jīng)成為現(xiàn)階段亟待解決的問題之一。

1 材料與方法

1.1 FMDV 分布數(shù)據(jù)的收集和處理FMDV 病例的地理位點(diǎn)數(shù)據(jù)，在NCBI、CNKI 以及OIE（世界動(dòng)物衛(wèi)生組織）中分別檢索“Foot-and-mouth disease virus”、“口蹄疫”，追蹤相關(guān)文獻(xiàn)及口蹄疫疾病報(bào)告。排除無法確定地理坐標(biāo)和重復(fù)數(shù)據(jù)。共收集到97個(gè)分布點(diǎn)。

1.2 環(huán)境數(shù)據(jù)收集和處理環(huán)境數(shù)據(jù)來源于“野生動(dòng)物疫源疫病研究平臺(tái)”數(shù)據(jù)庫，從中獲取高分辨率氣候數(shù)據(jù)，空間分辨率為30arc-seconds 1 km。所使用的67個(gè)環(huán)境變量分別為全年12個(gè)月的平均最低溫度（Tmin1～12）、全年12個(gè)月的平均最高溫度（Tmax1～12）、 全年12個(gè)月 的 平 均 溫 度（Temp1～12）、全年12個(gè)月的降水量（Prec1～12）和19個(gè)生物氣候變量（Bio1～19）。

1.3 環(huán)境數(shù)據(jù)篩選使用ArcGis 中Spatial Analyst工具提取97個(gè)分布點(diǎn)對(duì)應(yīng)的環(huán)境屬性值，并采用SPSS 22 軟件對(duì)FMDV 的環(huán)境數(shù)據(jù)進(jìn)行主成分分析（Principal Component Analysis，PCA）。根據(jù)主成分個(gè)數(shù)的選取原則，主成分的特征值是主要影響主成分大小的指標(biāo)，如果特征值小于1，則表明該主成分的綜合解釋力度不足。本研究中選取特征值大于1 的主成分，并選擇該主成分中載荷值大于0.95 的變量。通過上述準(zhǔn)則對(duì)環(huán)境變量進(jìn)行篩選，將篩選后得到的氣象變量使用SPSS 22 軟件進(jìn)行共線性診斷，若依然存在多重共線性，則將篩選后得到的氣象變量輸入MaxEnt 模型中進(jìn)行分析，依次剔除貢獻(xiàn)率小于1、5、10、20 的環(huán)境變量，將篩選后得到的環(huán)境變量再次進(jìn)行共線性診斷，直到VIF小于10，即可停止。本研究中投影坐標(biāo)系為GCS_WGS_1984。

1.4 FMDV 風(fēng)險(xiǎn)分析方法使用MaxEnt 3.4.1 版 本（http://www.cs.princeton.edu/～schapire/maxent/），將采集的數(shù)據(jù)按測(cè)試集30%和訓(xùn)練集70%隨機(jī)分組，選擇刀切法來評(píng)估環(huán)境因子的貢獻(xiàn)率，重復(fù)次數(shù)為3次；并采用ROC 曲線和AUC 值來評(píng)價(jià)MaxEnt 模型的精確度，運(yùn)行得到FMDV 空間分布預(yù)測(cè)圖。

2 結(jié)果與分析

2.1 環(huán)境變量確定及統(tǒng)計(jì)分析通過PCA（主成分分析法）共篩選得到9個(gè)相關(guān)性較大的變量進(jìn)行MaxEnt 模型運(yùn)算（見表1）。環(huán)境因子分為2 類，分別是溫度相關(guān)和降水量相關(guān)的環(huán)境變量。

表1 相關(guān)性≥0.95 的環(huán)境變量

2.2 AUC 檢測(cè)值及環(huán)境變量貢獻(xiàn)率結(jié)果分析

（1）MaxEnt 模型檢驗(yàn)預(yù)測(cè) 利用MaxEnt 最大熵模型多次建模，利用方差膨脹因子VIF 對(duì)環(huán)境因子進(jìn)行檢驗(yàn)，保證VIF＜10。經(jīng)過多次運(yùn)算，最終保留Temp6、Tmax5 以及Bio5 進(jìn)行預(yù)測(cè)。MaxEnt 模型重復(fù)運(yùn)算3次，最終得到的模型AUCtest 的平均值為0.919＞0.9，具有較高可信度（見圖1）。

圖1 MaxEnt 模擬預(yù)測(cè)口蹄疫的精度分析曲線

（2）各變量貢獻(xiàn)率及單一因子分析 如表2所示，Temp6、Tmax5、Bio5 這3個(gè)環(huán)境變量對(duì)FMDV具有較大的貢獻(xiàn)率，貢獻(xiàn)率均在10%以上。其中，Temp6、Tmax5 均為溫度變量，貢獻(xiàn)率累積達(dá)到85.3%。

表2 MaxEnt 模型中環(huán)境變量相對(duì)貢獻(xiàn)率分析

由響應(yīng)曲線以及頻數(shù)分布圖可知，在6月平均溫度為14℃左右時(shí)，口蹄疫發(fā)病率最高（圖2-a）；在5月最高溫度為28℃左右時(shí)，口蹄疫發(fā)病率最高（圖2-b）；在最熱月的最高溫度為32℃左右時(shí)，口蹄疫發(fā)病率最高（圖2-c）。

圖2 -c Bio5 環(huán)境變量與口蹄疫關(guān)系反應(yīng)曲線圖

圖2 -b Tmax5 環(huán)境變量與口蹄疫關(guān)系反應(yīng)曲線圖

圖2 -a Temp6 環(huán)境變量與口蹄疫關(guān)系反應(yīng)曲線圖

2.4 東北虎豹公園FMDV風(fēng)險(xiǎn)分布分析由圖3可以看出，虎豹公園中部及南部部分區(qū)域處于高風(fēng)險(xiǎn)水平（＞0.7，深黃色至紅色區(qū)域）；西北部分地區(qū)處于中風(fēng)險(xiǎn)水平（0.5～0.7，淡黃色至黃色區(qū)域）；南部地區(qū)（琿春市南部）處于低風(fēng)險(xiǎn)水平（＜0.4，淺藍(lán)色至藍(lán)色區(qū)域）；西部邊界部分地區(qū)、北部部分地區(qū)以及東北部分地區(qū)處于低風(fēng)險(xiǎn)水平（＜0.4，淺藍(lán)色至藍(lán)色區(qū)域）。

圖3 虎豹公園口蹄疫疫情風(fēng)險(xiǎn)水平分布圖

3 討 論

本研究所用方法是采取PCA 法對(duì)環(huán)境因子進(jìn)行篩選，在景觀尺度上利用MaxEnt 等手段，將口蹄疫風(fēng)險(xiǎn)分布區(qū)環(huán)境選擇機(jī)制與環(huán)境因素耦合起來，綜合評(píng)估流行風(fēng)險(xiǎn)。

由主成分分析可知與口蹄疫分布相關(guān)的第一主成分為溫度主成分。溫度對(duì)FMDV 的發(fā)生均成正相關(guān)性。Temp6、Tmax5、Bio5 這3個(gè)環(huán)境變量對(duì)FMDV 影響最大。3個(gè)環(huán)境變量均與溫度有關(guān)，其中，Tmax7、Tmax9 累計(jì)貢獻(xiàn)率達(dá)到78.4%。由模型最終結(jié)果可得，口蹄疫分布主要與溫度有關(guān)，且以溫暖季節(jié)發(fā)病較多，這主要與野生動(dòng)物在溫暖季節(jié)行為活動(dòng)比較頻繁、患病動(dòng)物接觸密切以及口蹄疫具有較高的傳染性有關(guān)。東北虎豹公園應(yīng)全面掌握口蹄疫流行態(tài)勢(shì)、分布規(guī)律，并根據(jù)風(fēng)險(xiǎn)評(píng)估制定針對(duì)性的口蹄疫防治計(jì)劃，建立完善的綜合防治措施與政策，并在國家公園內(nèi)進(jìn)行定期監(jiān)測(cè)、分區(qū)防治以及生態(tài)健康綜合治理，并重視對(duì)口蹄疫跨境傳播的防治與監(jiān)測(cè)，降低疫病的輸入性風(fēng)險(xiǎn)。同時(shí)，國家公園應(yīng)當(dāng)加強(qiáng)對(duì)管理人員的技能培訓(xùn)，提高野生動(dòng)物源疫病的防治水平，最終實(shí)現(xiàn)對(duì)口蹄疫有效控制的目標(biāo)。