Scientia Geographica Sinica  2013 , 33 (1): 123-128 https://doi.org/10.13249/j.cnki.sgs.2013.01.123

Orginal Article

湘江中下游肾综合征出血热传播风险预测和环境危险因素分析

肖洪1, 林晓玲1, 高立冬2, 代翔宇1, 贺新光1, 陈碧云2, 张锡兴3, 赵暕4, 田怀玉1

1. 湖南师范大学资源与环境科学学院, 湖南 长沙 410081
2. 湖南省疾病预防控制中心, 湖南 长沙 410005
3. 长沙市疾病预防控制中心, 湖南 长沙 410001
4. 北京大学医学部, 北京 100191

Environmental Factors Contributing to the Spread of Hemorrhagic Fever with Renal Syndrome and Potential Risk Areas Prediction in Midstream and Downstream of the Xiangjiang River

XIAO Hong1, LIN Xiao-ling1, GAO Li-dong2, DAI Xiang-yu1, HE Xin-guang1, CHEN Bi-yun2, ZHANG Xi-xing3, ZHAO Jian4, TIAN Huai-yu1

1. College of Resources and Environment Science, Hunan Normal University, Changsha, Hunan 410081, China
2. Hunan Provincial Center for Disease Control and Prevention, Changsha, Hunan 410002, China
3. Changsha Municipal Center for Disease Control and Prevention, Changsha, Hunan 410001, China
4. Peking University Health Science Center, Beijing 100191, China

中图分类号:  R188

文献标识码:  A

文章编号:  1000-0690(2013)01-0123-06

收稿日期: 2012-05-5

修回日期:  2012-07-22

网络出版日期:  2013-01-20

版权声明:  2013 《地理科学》编辑部 本文是开放获取期刊文献,在以下情况下可以自由使用:学术研究、学术交流、科研教学等,但不允许用于商业目的.

基金资助:  湖南省重点学科建设项目(2008001)、湖南省自然科学基金(11JJ3119)、湖南省高校创新平台开放基金(11K037)、湖南师范大学地图学与地理信息系统校级重点学科(2011001)资助世界野生动物基金会

作者简介:

作者简介:肖 洪(1964-),女,湖南湘潭县人,博士,副教授,主要研究方向为空间流行病学。E-mail:xiaohong.hnnu@gmail.com

展开

摘要

利用2005~2010年湘江中下游地区HFRS病例数据和相关环境数据,结合地理信息系统、遥感技术和最大熵值法生态位模型探索湘江中下游地区HFRS传播风险和主要危险因素。研究结果表明,望城县和长沙县北部,衡东县、株洲县和湘潭县交界处以及衡东县南部等区域为潜在风险区。5和7月的NDVI值在 0.3~0.4之间的区域HFRS传播风险高,城镇和建筑用地是主要风险用地类型。生态位模型结合地理景观、气象条件及人类活动数据对湘江中下游地区HFRS传播进行分析,较其他模型获得的结果准确率更高。

关键词: 肾综合征出血热 ; 生态位模型 ; 最大熵值法 ; 湘江中下游

Abstract

HFRS (human hemorrhagic fever with renal syndrome) occurs at very high frequency in the midstream and downstream of the Xiangjiang River. It is of great significance for understanding the regulation of the transmission of HFRS in this area due to its major impact on local economy development and public health. The risk factors and potential risk of HFRS in the midstream and downstream of theXiangjiang River were explored by combining the ecologic niche modeling, geographic information systems (GIS) and remote sensing technique in this article. The average area under the Receiver Operating Characteristic curve was 0.785. The 64 cases in 2010 were predicted correctly except for one case in a low risk area. These showed that the model possess good prediction performance. The modeling results indicated that the north of Wangcheng County and Changsha County, the border of Hengdong County , Zhuzhou County and Xiangtan Countiy, and the south of Hengdong County were predicted as the potential risk area. The following factors were found to be closely relative to the highest risk of HFRS transmission: mean annual temperature of 18℃, annual precipitation of about 1 500 mm, at urban and constructive land, and normalized difference vegetation index between 0.3 and 0.4 in May and July and with low elevation. The transmission of HFRS in the midstream and downstream of the Xiangjiang River was analyzed by using the ecologic niche model combining the geographical environment, climatic condition and human activity data. Our results are more accurate than those from other models and these findings in the study. Therefore, they are applicable for targeting control and prevention efforts.

Keywords: hemorrhagic fever with renal syndrome ; ecologic niche modeling ; maximum entropy ; midstream and downstream of the Xiangjiang River

0

PDF (722KB) 元数据 多维度评价 相关文章 收藏文章

本文引用格式 导出 EndNote Ris Bibtex

肖洪, 林晓玲, 高立冬, 代翔宇, 贺新光, 陈碧云, 张锡兴, 赵暕, 田怀玉. 湘江中下游肾综合征出血热传播风险预测和环境危险因素分析[J]. , 2013, 33(1): 123-128 https://doi.org/10.13249/j.cnki.sgs.2013.01.123

XIAO Hong, LIN Xiao-ling, GAO Li-dong, DAI Xiang-yu, HE Xin-guang, CHEN Bi-yun, ZHANG Xi-xing, ZHAO Jian, TIAN Huai-yu. Environmental Factors Contributing to the Spread of Hemorrhagic Fever with Renal Syndrome and Potential Risk Areas Prediction in Midstream and Downstream of the Xiangjiang River[J]. Scientia Geographica Sinica, 2013, 33(1): 123-128 https://doi.org/10.13249/j.cnki.sgs.2013.01.123

肾综合征出血热(hemorrhagic fever with renal syndrome,以下简称HFRS)主要是以携带汉坦病毒(hantavirus,HV)的鼠类为主要传染源的自然疫源性疾病,中国是受汉坦病毒危害最严重的国家,近10 a来报告发病例数占全球病例总数的90%,其中感染流行的主要为HFRS[1]。作为一种典型的以鼠类为主要传染源的自然疫源性疾病,HFRS的分布、传播和流行受到地理景观[2~5]、气候[6~8]和人类活动[9]等自然和社会因素的影响。生态位模型(Ecologic niche models,ENMs)根据疾病传播中媒介、宿主、病原体或人群病例信息,结合相关地理景观和生态要素的变化,探索疾病传播中潜在的地理和生态需求,近年来逐渐被用于分析地理与生态因素的变化和疾病传播流行的关系[10,11],取得了很好的效果[12~14]。生态位模型可以有效地预测区域内疾病的潜在分布区,以便对其进行提前预防,也可以筛查影响疾病在区域内传播的主要风险因素,为疾病的侧重防控提供决策支持,这对自然疫源性疾病的防控工作具有很高的理论与现实意义,可以有效地降低各类损失。由于复杂多变的地形地貌、生态环境和气候特点[15],以及各种人类活动的影响,中国的HFRS疫情存在显著的地区性差异[16],掌握疫情高发区的实际情况十分必要。本研究以湘江中下游区域为例,结合地理信息系统(Geographic Information System,GIS)和遥感(Remote Sensing,RS)构建最大熵值法生态位模型探索近年来湘江中下游HFRS的传播风险,并分析相关的危险因素。

1 资料与方法

1.1 研究区概况

湘江发源于广西东北部,东流经湖南省永州市、衡阳市、株洲市、湘潭市、长沙市,在岳阳市湘阴县入洞庭湖后归于长江,全长817 km,流域面积92 300 km2,为湖南省最大的河流。湖南省处于海拔小于500 m的多水带和过渡带,地势以丘陵为主,气候温和多雨,农作物主要以亚热带双季稻为主,湘江流域为湖南省最发达的区域,湿润的地理环境和密集的人类活动为鼠类的生存活动和HFRS的人间传播提供天然的平台。另外,湖南省是中国HFRS高发病率的典型疫区,其中湘潭市和长沙市均为湖南省HFRS的高发区之一[17~20]。因此,本文选取湘江中下游区域,包括衡阳县,衡东县,株洲县,株洲市区,湘潭县,湘潭市区,望城县,长沙县,长沙市区,湘阴县和汨罗市等8个县级行政单位,3个市辖区为研究区域(图1),对HFRS的传播规律进行研究,具有重要的理论和现实意义。

图1   研究区地理位置

Fig.1   Location of the study area

1.2 研究数据

2005~2010年HFRS病例数据来自湖南省疾病预防控制中心疾病监测信息报告管理系统,所有病例均为实验室诊断数据[21],病例信息包括性别、出生日期、年龄、详细居住地、职业、发病日期等,病例数据未区分汉坦病毒类型。

使用ArcGIS 9.3(Environmental Sciences Research Institute,Redlands,CA)软件建立湘江中下游地理空间数据库(表1),包括地理编码、HFRS病例信息数据、复合地形指数(Compound topographic index,CTI)、人迹指数、主要生境(Major habitat types,MHT)、2005 ~2009年年平均气温和年降水量、逐年土地利用变化和逐月NDVI遥感影像数据,以及湘江中下游数字高程模型(Digital elevation model,DEM),其中NDVI数据取2005~2009年的月平均值,所有栅格数据空间分辨率为1 km。

表1   生态位模型环境变量

Table 1   Environmental variables of ecologic niche models

变量名称数据来源数据类型
土地利用MODIS遥感影象提取与分类分类型
CTI美国地质调查局连续型
MHT世界野生动物基金会分类型
人迹指数国际地球科学信息中心连续型
生态系统世界环保组织分类型
温度中国气象科学数据服务共享网连续型
降水量中国气象科学数据服务共享网连续型
NDVITM遥感影响提取连续型
高程国际科学数据服务平台连续型
坡度高程数据中提取连续型

1.3 最大熵值法生态位模型通过构建生态位模型[10,11,22],根据已知病例分布数据与相关环境变量(如用地类型和气象因子),探索疾病传播对地理景观和气象因子的潜在需求,对区域内的疾病传播风险进行模拟和预测[23]。本研究使用最大熵值法生态位模型探索已知病例分布区环境特征与研究区的非随机关系,根据最大熵模型,选择最优分布,对疫情进行预测[23,24]。假设研究区X是由有限数量的网格单元所组成,π为疫病在研究区内的分布,π(x)为π的分布赋予每个单元x的值,总和为1,在满足已知各种环境变量fjj=1,2,…,n)的条件下(式2和3),根据π的近似期望分布π^的熵(式1)选择熵最大的分布作为疫病的最优分布。

新窗口打开

Hπ^=-xXπ^(x)lnπ^(x)(1)

π^[fj]=π˜[fj](2)

π˜[fj]=1mi=1mfj(xi)(3)

其中, Hπ^为期望分布 π^的熵, xi为研究区X中的第i个单元,m为单元个数,fj为环境变量, π˜[fj]为样本点环境变量在 π˜分布下的先验均值。

由上可知最大熵值法概率分布与吉布斯概率分布(式4)是相同的,均是最大化所有样本点的相似性,降低损失函数(式3)。由于Maxent所用的变量为经验性变量,而不是确定的实际值,可能会导致其对训练数据的过度拟合,为减小过度拟合,需要适当的放宽环境变量的限制(式6)。同时,为了使Maxent能有效选取重要的环境变量,使用最小化对数损失和限制环境变量过大权重λj(式7)的吉布斯分布来表示最大熵值概率分布。

qλx=eλf(x)Zλ(4)

π˜-ln(q)λ(5)

其中,λ是表示所有环境变量权重的向量,f是所有环境变量的向量,Zλ为标准化常数,用于确保qλ的总和为1。

π^[fj]-π˜[fj]βj(6)

π˜-ln(q)λ+jβjλj(7)

其中, βj为常数, π˜-ln(qλ)是损失函数, jβjλj表示对环境变量权重的限制。

2 生态位模型预测结果

根据第一次运算结果,选择月平均NDVI值(包括3、5、7、8和9月,累计贡献率为37.0%),土地利用类型(包括2005、2006、2007和2008年,累计贡献率13.2%),年平均温度(贡献率19.8%),年降水量(贡献率11.7%),人迹指数(贡献率7.4%),坡度(贡献率4.8%)和CTI(贡献率1.1%)等影响HFRS传播的主要生态环境变量,对湘江中下游HFRS的传播风险进行预测,接受者操作特性曲线(Receiver operating characteristic,ROC)平均线下面积(Area under the curve,AUC)值为0.785(SD = 0.032),表明模型预测效果较好。取10个预测模型的平均预测风险等级图为最佳预测图(图2),其中黑线内区域表示潜在风险区,即尚无病例出现但风险较高的区域[25],结果表明,传播高风险区预测与实际病例分布相符合,2010年预留的64例测试病例数据,除1例预测为低风险和13例预测为中风险,其余均位于高风险区,望城县和长沙县的北部,衡东县、株洲县和湘潭县交界处以及衡东县南部等被预测为潜在风险区(图2)。

图2   湘江中下游HFRS发病情况与生态位模型预测结果

Fig. 2   The distribution of human hemorrhagic fever with renal syndrome (HFRS) cases in 2005-2009 and the predicted HFRS risk in 2010 by Maximum Entropy Model in midstream and downstream of the Xiangjiang River

3 环境危险因素分析

研究结果表明,HFRS的传播风险在年平均温度16℃左右开始上升,在18℃左右达到最大,然后一直下降到18.3℃左右保持不变,在年降水量为1 350 mm和1 500 mm左右急剧上升,在1 500 mm左右达到最大后持续下降(图3a和b,图中温度为真实值乘以10后得到的值)。同时HFRS的传播风险随着坡度的上升而下降,在坡度大于10°的区域,HFRS的传播风险基本为0,随着人迹指数的增大而持续增加(图3e和f)。城镇和建筑用地(用地类型4)是HFRS传播的主要用地风险(图3g、h、i和j)。当5月和7月的NDVI 值在0.3~0.4之间时HFRS传播风险最高,且HFRS在NDVI低值区域的传播风险5月大于7月,5月高峰期后在NDVI值为0.7左右有个小高峰(图3c和d)。

降水量可以影响湿度条件和鼠类生存环境,湿润的环境可以增强汉坦病毒的生命力和感染性,适量的降水有利于植被的生长,对鼠类的存活和繁殖有着促进作用,因此适度的降水有助于HFRS的传播,但过多或持续的降水则会破坏鼠类的生存环境,导致其数量的减少,降低鼠类之间的相互接触[8,26],HFRS传播风险在年降水量1 350和1 500 mm的区域急剧上升(图3b)进一步证明鼠类对降水量的敏感。温度通过影响怀孕率、胎仔数、出生率和成活率影响鼠类的种群数量,相关研究表明,作为HFRS主要动物传染源的黑线姬鼠和褐家鼠的怀孕盛期均受到温度影响[7] ,适宜的温度有利于鼠类数量增加。研究表明,年均温度和年降水量是影响湘江中下游HFRS传播的重要气象因子,HFRS高传播风险主要集中在年均温度18℃,降水量1 500 mm左右的区域(图3a、b)。

NDVI是反应地表植被覆盖程度的变量,可综合反应区域的地理环境特征[27, 28],HFRS的传播和NDVI有着密切关系[29,30]。5月和7月的NDVI值对HFRS传播有重要影响,且HFRS在NDVI低值区传播风险5月高于7月(图3c和d),这可能是湘江区域传统的农耕方式和其对鼠类生存环境的影响,7月是双季稻早稻收获的时候,植被覆盖率较高,同时也有利于野外优势鼠种黑线姬鼠的生存和活动。人类活动,包括森林砍伐、道路修建、农业活动、水坝修建、开矿等造成的用地类型的改变,对传染病的暴发和重现有重要影响[31, 32],城市化带来的密集的人口也可以增加感染者和易感染者的接触加快传染病的传播。HFRS的传播风险在地势平坦的区域最高(图3e),这可能与低坡度区域的土地利用情况相关。城镇和建筑用地是影响HFRS传播的主要用地风险,HFRS的传播风险随着人迹指数的增加而增加(图3f、g、h、i、j),说明人类活动对HFRS持续传播的潜在驱动。

图3   湘江中下游环境对肾综合征出血热传播的响应响应曲线

Fig. 3   Response curves for the variables related to presence of HFRS(Red lines are mean values for the 10 Maxent runs and blue bars represent ±1 SD)

4 结果与讨论

本文结合GIS技术、遥感技术和Maxent生态位模型有效地获得了湘江中下游HFRS传播的主要危险因素,同时也成功地预测了HFRS在湘江中下游的潜在风险地区。研究结果表明,湘江中下游区域Maxent生态位模型预测HFRS传播高风险区与实际病例分布相符合,对望城县和长沙县的北部,衡东县、株洲县和湘潭县交界处以及衡东县南部等HFRS潜在风险区的预测(图2),可为区域HFRS疫情的提前防控提供理论依据,同时,NDVI和土地利用对HFRS的传播有着重要影响,可结合NDVI和土地利用对HFRS的传播规律做进一步分析。本文将生态位应用于传染病研究中,有效地解决传统传染病研究在流行强度预测和传播区域预测相结合上的不足[34],但是,本研究单纯考虑地理景观和生态环境对HFRS的传播影响,未结合区域宿主动物的种群分布和带病毒情况,同时传染病的流行与传播的存在时空差异性,如何选择适宜的环境变量也有待进一步深入思考研究。

The authors have declared that no competing interests exist.


参考文献

[1] 张永振,肖东楼,王玉,.

中国肾综合征出血热流行趋势及其防治对策

[J].中华流行病学杂志, 2004, 25(6): 466~469.

[本文引用: 1]     

[2] Chen H X, Qiu F X.

Epidemiologic surveillance on the hemorrhagic fever with renal syndrome in China

[J]. Chin Med J (Engl), 1993, 106(11): 857-863.

[本文引用: 1]     

[3] Chen H X,Qiu F X.

Studies on the environment structure of natural nidi and epidemic areas of hemorrhagic fever with renal syndrome in China

[J]. Chin Med J (Engl), 1994, 107(2):107-112.

[4] Chen H X,Qiu F X,Dong B J,et al.

Epidemiological studies on hemorrhagic fever with renal syndrome in China

[J]. J Infect Dis,1986,154(3): 394-398.

[5] 肖洪,田怀玉,代翔宇,.

地理景观对长沙市肾综合征出血热传播的影响

[J].中华预防医学杂志,2011,46(3): 246~251.

[本文引用: 1]     

[6] Zhang W Y, Fang L Q, Jiang J F, et al.

Predicting the risk of hantavirus infection in Beijing, People's Republic of China

[J]. Am J Trop Med Hyg, 2009,80(4):678-683.

[本文引用: 1]     

[7] 刘静,王洁贞,薛付忠,.

肾综合征出血热发病率与气象因素关系的研究

[J].中国卫生统计,2006, 23(4): 326~329.

[本文引用: 1]     

[8] 肖洪,田怀玉,张锡兴,.

气候变化对长沙市肾综合征出血热发病的影响与预警模型

[J].中华预防医学杂志, 2011, 45(10): 881~885.

[本文引用: 2]     

[9] 阮玉华,柳炜,徐校平,.

流行性出血热发病危险因素病例对照研究

[J].中华流行病学杂志,2000, 21(5): 344~346.

[本文引用: 1]     

[10] Peterson A T.

Ecologic niche modeling and spatial patterns of disease transmission

[J].Emerg Infect Dis, 2006, 12(12): 1822-1826.

[本文引用: 2]     

[11] Peterson A T.

Ecological niche modelling and understanding the geography of disease transmission

[J].Vet Ital, 2007, 43(3): 393-400.

[本文引用: 2]     

[12] Peterson A T, Bauer J T, Mills J N.

Ecologic and geographic distribution of filovirus disease

[J].Emerging Infectious Diseases,2004,10(1): 40-47.

[本文引用: 1]     

[13] Peterson A T,Sánchez-Cordero V, Beard C B,et al.

Ecologic niche modeling and potential reservoirs for Chagas disease, Mexico

[J]. Emerging Infectious Diseases, 2002, 8(7): 662.

[14] Fang L Q,

De Vlas S J, Liang S, et al. Environmental factors contributing to the spread of H5N1 avian influenza in mainland China

[J]. PLoS One, 2008, 3(5): e2268.

[本文引用: 1]     

[15] 罗成旺,陈化新.

肾综合征出血热流行因素的影响研究

[J].中国媒介生物学及控制杂志,2003,14(6): 451~454.

[本文引用: 1]     

[16] 王芹,周航,李德新,.

2009年中国肾综合征出血热监测分析

[J].疾病监测,2010,25(12):938~943.

[本文引用: 1]     

[17] Yan L,Fang L Q,Huang H G,et al.

Landscape elements and Hantaan virus-related hemorrhagic fever with renal syndrome, People's Republic of China

[J]. Emerg Infect Dis,2007,13(9):1301-1306.

[本文引用: 1]     

[18] 刘富强,高立冬,戴德芳,.

2002~2005年湖南省人群肾综合征出血热监测分析研究

[J].疾病监测, 2007, 22(7): 452~455.

[19] 刘富强,高立冬,戴德芳,.

湖南省2007年肾综合征出血热监测结果分析

[J].中国自然医学杂志, 2008, 10(3): 199~202.

[20] 罗成旺,陈化新.

中国1998~2007年肾综合征出血热流行病学特征及疫苗接种策略探讨

[J].中华流行病学杂志, 2008, 29(10):1017~1019.

[本文引用: 1]     

[21] 中华人民共和国卫生部.

扩大国家免疫规划实施方案

[M].2007.

[本文引用: 1]     

[22] Guisan A,Zimmermann N E.

Predictive habitat distribution models in ecology

[J].Ecol Model,2000,135(2-3):147-186.

[本文引用: 1]     

[23] Wei L,Qian Q,Wang Z Q,et al.

Using Geographic Information System-based Ecologic Niche Models to Forecast the Risk of Hantavirus Infection in Shandong Province,China

[J].Am J Trop Med Hyg, 2011,84(3):497-503.

[本文引用: 2]     

[24] Phillips S J,Anderson R P,Schapire R E.

Maximum entropy modeling of species geographic distributions

[J].Ecological Modelling, 2006,190(3-4): 231-259.

[本文引用: 1]     

[25] Peterson A T,Papes M,Soberón J.

Rethinking receiver operating characteristic analysis applications in ecological niche modeling

[J].Ecological Modelling,2008,213(1):63-72.

[本文引用: 1]     

[26] Fang L Q,Wang X J,Liang S, et al.

Spatiotemporal trends and climatic factors of hemorrhagic fever with renal syndrome epidemic in Shandong Province, China

[J]. PLoS Negl Trop Dis,2010, 4(8): e789.

[本文引用: 1]     

[27] 李惠敏,刘洪斌,武伟,.

近10年重庆市归一化植被指数变化分析

[J].地理科学,2010, 30(1): 119~123.

[本文引用: 1]     

[28] 南颖,刘志锋,董叶辉,.

2000~2008年长白山地区植被覆盖变化对气候的响应研究

[J].地理科学, 2010, 30(6): 921~928.

[本文引用: 1]     

[29] 方立群,曹务春,陈化新,.

应用地理信息系统分析中国肾综合征出血热的空间分布

[J].中华流行病学杂志, 2003, 24(4): 265~268.

[本文引用: 1]     

[30] 闫磊,黄华国,张文义,.

肾综合征出血热疫情与NDVI 的时间关系——以内蒙古自治区大杨树镇为例

[J].遥感学报, 2009, 13(5): 880~886.

[本文引用: 1]     

[31] 杨坤,周晓农.

景观流行病学研究现状及其进展

[J].中华流行病学杂志, 2008, 29(2): 198~201.

[本文引用: 1]     

[32] 李雪铭,晋培育.

中国城市人居环境质量特征与时空差异分析

[J].地理科学,2012,32(5): 521~529.

[本文引用: 1]     

[33] 吴玉鸣,柏玲.

广西城市化与环境系统的耦合协调测度与互动分析

[J].地理科学,2011,31(12): 1474~1479.

[34] 肖洪,田怀玉,赵暕,.

传染病模型分析与预测方法研究进展

[J].中华流行病学杂志,2011,32(1): 81~85.

[本文引用: 1]     

/