分析2005-2014年盐城市沿海滩涂恙虫病的时空分布特征，并对气象、鼠类种群分布及生物学特点等影响因素进行相关性研究。
收集整理2005-2014年盐城市沿海滩涂地区恙虫病疫情、人口学以及气象学资料，并于2011年4月至2013年12月监测该地区鼠类种群分布及其密度、季节消长情况。采用Excel 2003、SPSS 16.0、MapInfo 11.0、Satscan 9.0和Stata/SE 10.0等软件，利用描述性方法、时空重排扫描统计量法、自相关及互相关等方法分析该地区恙虫病时空分布特征，及与鼠密度消长、气象因素变化的相关性，根据恙虫病发病数据建立Poisson回归模型。
（1） 2005-2014年盐城市沿海滩涂地区恙虫病的发病呈间歇性逐年升高，2014年发病率为最高（5.81/10万）。呈现秋季高峰，11月发病率最高（12.02/10万）；滨海、大丰、响水县／市年均发病率最高，分别为3.30/10万、3.21/10万和2.79/10万，沿海地区的12个乡镇为高发，发病率为4.41/10万至10.03/10万。（2）恙虫病可能存在3个发病聚集区域，分别为2012年10-11月大丰、东台和射阳县／市的25个乡镇、2014年10-11月响水县的5个乡镇以及2006年11月滨海县的6个乡镇。（3）优势鼠种为黑线姬鼠（89.19%）。鼠密度呈冬季（1月）与夏季（5-8月）2个高峰，分别为5.05%（1月）和8.05%（8月）。发病月分布为每年的10-11月有一次发病高峰，每次发病高峰前1~ 2个月有一个鼠密度升高趋势，每次发病高峰前均有一个气温、降水量高峰。鼠密度消长、气温、降水量与恙虫病发病有统计学相关性。根据ZIP回归模型分析发病例数与降水、平均最低气温、滞后3个月平均最低气温呈正相关，与日照时长和相对湿度呈负相关。
盐城市沿海滩涂地区存在恙虫病疫区，且有时空聚集性，并存在媒介生物和气象等流行因素，应采取有效防控措施。
版权归中华医学会所有。
未经授权，不得转载、摘编本刊文章，不得使用本刊的版式设计。
除非特别声明，本刊刊出的所有文章不代表中华医学会和本刊编委会的观点。
恙虫病是一种自然疫源性疾病，其主要贮存宿主为鼠类，传播媒介为恙螨。目前研究发现，我国的恙虫病按照流行季节、流行范围、宿主动物、媒介恙螨、病原体毒力和临床表现分为夏季型、秋季型和冬季型三种季节类型[1]。并已发现，其流行特征受气象参数等环境因素[2]和宿主、传播媒介及社会经济因素等影响。盐城市沿海滩涂除有繁茂的自然植被外，还拥有140多种脊椎动物及300余种昆虫资源[3,4,5,6]，覆盖东台、大丰、射阳、滨海和响水县／市的53个乡镇200万人。近年来上述地区恙虫病常年发生，偶有暴发[7]。为此对2005-2014年盐城市沿海滩涂地区恙虫病的时空分布特征及其鼠类宿主种群分布、密度消长和气象影响因素进行分析。
从中国疾病预防控制信息系统以乡镇为单位收集2005-2014年东台、大丰、射阳、滨海和响水县／市恙虫病的发病数据，人口数据源于盐城市统计部门。从盐城市乡镇电子地图（1∶10 000）中提取上述县／市相关乡镇空间数据。收集整理2011年4月至2013年12月该地区鼠类种群分布与密度消长监测资料。收集2011年1月至2013年12月盐城市气象资料，主要气象指标包括月平均气温（某月逐日平均气温的平均值）、月平均最高气温（某月逐日最高气温的平均值）、月平均最低气温（某月逐日最低气温的平均值）、相对湿度、降水量、降水天数、日照时间等。
统计盐城市沿海滩涂地区县和乡镇两级的恙虫病疫情与人口资料，将每个乡镇构建唯一编号，建立空间分析的属性数据库。根据盐城市沿海乡镇空间数，建立与属性数据库中乡镇相匹配的编码，构建空间数据库。按照盐城市沿海滩涂分布情况，随机选择滨海县滨海港范围内的港区、镇区、二罾（振东镇），盐城市国家级丹顶鹤自然保护区，大丰市大丰港区、大丰国家级麋鹿自然保护区，东台市金东台农场笆斗高滩围垦区、三仓垦区、国华风电场和港城作为病媒生物种群分布与密度消长监测点，按照《全国病媒生物监测方案（试行）》，于2011年4月至2013年12月监测沿海滩涂地区鼠类种群分布、密度、季节消长。
采用SPSS 16.0软件描述性统计方法从时间和空间方面分析2005-2014年盐城市沿海滩涂地区恙虫病流行特征；利用MapInfo 11.0软件建立空间地理发病率数据，并采用Satscan 9.0软件分析恙虫病时空聚集性；利用Stata/SE 10.0软件互相关方法分析恙虫病时空分布与鼠密度消长、气象指标变化的相关性；利用SPSS 16.0软件描述性分析、零检验和自相关分析等方法分析2011-2013年恙虫病月发病报告数据特征，并建立气象因素对恙虫病发病影响的回归模型。
2005-2014年盐城市沿海滩涂地区恙虫病总发病率为21.74/10万，10年间发病呈间歇性逐步升高，最高发病为2014年，发病率为5.81/10万（图1）。恙虫病呈现秋季高峰，且病例主要集中在10月和11月，10年累计发病分别为183例（42.07%）和241例（55.40%）。
滨海、大丰、响水县／市年均发病率高，分别为3.30/10万、3.21/10万、2.79/10万，东台市年均发病率为1.46/10万，射阳县最低为0.91/10万，差异有统计学意义（χ2＝100.94，P<0.01）。沿海地区的12个镇发病率较高（4.41/10万至10.03/10万）。
将恙虫病发病数据和空间数据库分别导入SaTScan 9.4软件。分析显示恙虫病可能存在3个发病聚集区域：①2012年10-11月东台区域2个乡镇、大丰区域12个乡镇、射阳区域11个乡镇恙虫病发病存在聚集性，发病66例，似然比（LLR）＝103.44，RR＝12.95，P<0.001。②2014年10-11月响水区域的老舍、南河、双港、七套、大有镇恙虫病发病存在聚集性，发病42例，LLR＝85.66，RR＝20.87，P<0.001。③2006年11月滨海区域的淮海农场及振东、临海、八滩、滨海港、五汛镇恙虫病发病存在聚集性，发病33例，LLR＝55.35 ，RR＝ 14.01，P<0.001（图2）。
注：★代表聚集区域的乡镇
2011-2013年沿海滩涂地区的鼠种分布有褐家鼠、小家鼠、黑线姬鼠。其中黑线姬鼠为优势鼠种，构成比为89.19%（429/481）；褐家鼠、小家鼠和鼩鼱的构成比分别为4.57%（22/481）、2.70%（13/481）、3.53%（17/481）。鼠密度呈现冬季与夏季2个高峰（图3）。冬季高峰月份为1月，鼠密度为5.05%；夏季高峰自5月持续至8月，鼠密度最高的8月为8.05%。2011-2013年沿海地区恙虫病发病的时间分布，以每年的10-11月有1个高峰，每次的发病高峰前1～ 2个月均有1个鼠密度升高趋势。
利用互相关分析鼠密度的消长与恙虫病发病关系，结果显示鼠密度消长与恙虫病发病（滞后3个月）互相关关系（相关系数为0.439）有统计学意义，而滞后的3个月相当于媒介螨虫生活周期与恙虫病最短潜伏期之和的时长（恙螨的生活周期约3个月，恙虫病的最短潜伏期为4 d），见图4。
2011 － 2013年月均气温、平均最高气温、平均最低气温（tmin）的高峰月份均在7月，分别为28.13 ℃、31.73 ℃、25.13 ℃；平均相对湿度和平均日照时间的最高月份均在8月，相对湿度（hum）和日照时间（sun）分别为81.00%、169.37 h；平均降水量和降水天数最高月份均在7月，分别为311.67 mm和13.00 d。从气象因素与恙虫病月发病率分布图所示（图5），每次恙虫病发病高峰前均有一个气温、降水量（rain）高峰，而sun、hum则无上述一致性。
恙虫病月发病病例与气象因素的相关分析显示，气象因素中月均气温、tmin、平均最高气温3个变量满足正态分布，月均气温与月均最低气温的Pearson相关系数为0.997，月均气温与月均最高气温的Pearson相关系数为0.996，月均最低气温与月均最高气温的Pearson相关系数为0.988，变量间高度相关（P<0.001）。由于最低气温为恙螨卵生长发育的限制因素，故以tmin和rain作为气象影响因素，对2011-2013年恙虫病发病报告病例与气象因素进行互相关分析。结果显示，rain与恙虫病病例（滞后4个月）相互关系有统计学意义，而滞后的4个月相当于媒介螨虫生活周期与恙虫病最长潜伏期之和的时长（恙螨的生活周期约3个月，恙虫病最长潜伏期21 d）；tmin与恙虫病病例（滞后3个月）相互关系有统计学意义（表1）。
2011-2013年盐城市沿海滩涂地区恙虫病月报告病例与气象因素互相关关系
2011-2013年盐城市沿海滩涂地区恙虫病月报告病例与气象因素互相关关系
| 滞后月数 | rain | sun | hum | tmin |
|---|---|---|---|---|
| 0 | 0.012 | 0.023 | -0.116 | 0.057 |
| 1 | -0.101 | 0.191 | -0.095 | 0.009 |
| 2 | -0.197 | -0.135 | 0.184 | 0.105 |
| 3 | 0.063 | 0.205 | 0.161 | 0.374a |
| 4 | 0.763a | -0.097 | -0.062 | 0.220 |
注：aP<0.05
2011-2013年盐城市沿海滩涂地区报告恙虫病病例204例，均数为5.67例，方差为198.743，病例出现最多的1个月为54例，方差远大于均数，过度离散零检验统计量为29 137.84（P< 0.001），提示数据存在过度离散。2011-2013年恙虫病报告病例的累计月数中有25个月无病例报告，在全部观测值中零计数占69.44%，说明该数据中过度离散可能与零计数过多有关。报告病例的自相关分析显示月发病数间无相关性。
根据恙虫病发病报告病例的数据特征，选择建立具有Poisson分布类型的模型，通过非嵌套模型选择的Vuong检验方法，得出V＝1.98>1.96，可以认为选择ZIP模型建立2011-2013年盐城市沿海滩涂地区气象因素对恙虫病发病病例影响的回归模型更合适（表2）。
ZIP回归模型
ZIP回归模型
| 变量 | 系数 |  | Z值 | P值 | 95%CI | ||
|---|---|---|---|---|---|---|---|
| 方程一 | |||||||
| 常数项 | 0.605 | 0.398 | 1.520 | 0.128 | -0.175 ～1.385 | ||
| 方程二 | |||||||
| rain | 0.036 | 0.009 | 3.940 | 0.000 | 0.018～0.053 | ||
| Lag(rain, 4) | -0.002 | 0.001 | -1.150 | 0.248 | -0.004～0.001 | ||
| sun | -0.030 | 0.008 | -3.920 | 0.000 | -0.045～-0.015 | ||
| hum | -0.456 | 0.094 | -4.870 | 0.000 | -0.639 ~ -0.272 | ||
| tmin | 0.260 | 0.076 | 3.400 | 0.001 | 0.110 ～0.409 | ||
| Lag(tmin, 3) | 0.133 | 0.032 | 4.170 | 0.000 | 0.071 ～0.196 | ||
| 常数项 | 33.766 | 7.072 | 4.770 | 0.000 | 19.906 ～47.627 | ||
注：Lag （rain，4）为滞后4个月降水；Lag(tmin, 3)为滞后3个月平均最低温度
最终ZIP回归模型：
由此可知，恙虫病发病例数与降水量、平均最低气温、滞后3个月平均最低气温呈正相关，与日照时间和相对湿度呈负相关。
2005-2014年盐城市沿海滩涂地区恙虫病发病呈间歇性逐步升高，其中2014年发病率最高，并呈现秋季高峰，最高发病为11月（12.02/10万）。与文献报道一致[8]。采用时空聚集性分析表明[9,10]，该地区恙虫病可能存在3个发病聚集区域，即2012年10-11月在东台区域的2个乡镇、大丰区域的12个乡镇、射阳区域的11个乡镇，2014年10-11月在响水区域的5个乡镇，2006年11月在滨海区域的淮海农场及5个乡镇，并与高发乡镇疫情的分布具有一致性。近年来这些区域农业种植、工业规模化生产发展较快，土地利用率高，人员流动较大，提示存在传播恙虫病的流行因素。
盐城市沿海滩涂地区主要有褐家鼠、小家鼠、黑线姬鼠分布，其中黑线姬鼠为优势鼠种（89.19%），与秋冬型恙虫病的宿主动物分布吻合[11,12]。2011－ 2013年该地区鼠密度呈现冬季与夏季2个高峰。其中冬季高峰在1月，鼠密度为5.05%；夏季高峰持续时间为5-8月，8月鼠密度为最高（8.05%），且黑线姬鼠为绝对优势种群，结果与以往文献报道不同[13]。其差异可能是沿海滩涂地区鼠类生存环境与生物学特点有别于其他地区，尚有待进一步研究。
对盐城市沿海滩涂地区鼠类密度消长、气象参数与恙虫病发病的相互关系分析发现，2011-2013年沿海滩涂地区恙虫病每年10-11月有一次发病高峰，每次发病高峰前1~ 2个月有一个鼠密度升高。而每次恙虫病发病高峰前均有一个气温、降水量高峰，日照时间与相对湿度则无上述一致性。互相关分析鼠密度的消长、气象因素与恙虫病发病关系表明，鼠密度消长与恙虫病发病（滞后3个月）、降水量与恙虫病发病（滞后4个月）、平均最低气温与恙虫病发病（滞后3个月）互相关关系有统计学意义。ZIP回归模型可知，恙虫病发病例数与降水量、平均最低气温、滞后3个月平均最低气温呈正相关，与日照时间和相对湿度呈负相关。此外有研究显示，恙虫病除与媒介宿主密切相关及受气象环境因素影响外，还与人群社会活动方式等社会非环境因素有关[14]。
综上所述，盐城市沿海滩涂地区恙虫病呈现秋季高峰，在一些乡镇存在聚集性，其发生发展与鼠类种群分布、密度消长及气象因素可能有关。因此，在沿海滩涂开发过程中，要采取有针对性的预防控制措施，有效控制虫媒传染病的危害。
无
