土地覆盖变化对恙虫病传播风险的影响：基于泰国南府17年数据的生态流行病学研究

　　本文通过分析泰国南府2003-2019年间土地覆盖变化与恙虫病（Scrub typhus）发病率的关联，揭示了植被类型转换（如灌丛向阔叶林转化）对媒介螨（chigger mites）生态的显著影响。研究发现灌丛地（shrubland）和镶嵌林地（mosaic land）是高风险区域，而再造林过程呈现倒U型关系，为景观扰动促进病原体（Orientia tsutsugamushi）传播提供了重要证据。

　　

引言

恙虫病是由恙虫病东方体（Orientia tsutsugamushi）引起、通过恙螨属（Leptotrombidium）媒介传播的人兽共患病，临床表现为发热、焦痂及多器官损伤，未经治疗死亡率可达6%。全球年发病约100万例，主要流行于"恙虫病三角区"（日本、东南亚及澳大利亚）。泰国北部是重要 endemic area，2003-2017年报告103,345例，部分地区血清阳性率高达77%。现有研究表明土地覆盖变化可能通过改变媒介栖息地微环境（温度、湿度）及宿主（啮齿类）分布影响疾病传播，但缺乏村级尺度的长期生态流行病学研究。

研究方法

研究区域为泰国南府（19.15°N, 100.83°E），面积12,100 km²，包含905个村庄。采用2500米缓冲区分析村庄周边环境，整合以下数据：

1. 1.

   恙虫病病例数据：来自泰国国家疾病监测系统（2003-2019年）

2. 2.

   土地覆盖数据：欧洲空间局气候变化倡议（ESA CCI）300米分辨率遥感数据，重分类为9类（阔叶林、针叶林、灌丛地等）

3. 3.

   地形与人口数据：海拔、坡度及村庄男性人口（因性别差异显著）

   采用广义可加模型（GAM）分析土地覆盖组成与变化对年发病数的非线性影响，并通过空间+方法（spatial+ method）控制空间混杂效应。

结果

疾病负担时空特征

2003-2019年累计报告8,045例，发病呈上升趋势（2003年119例→2019年499例），2012年（906例）和2016年（721例）出现峰值。空间分析显示病例集中分布于西北部山地与丘陵区，南河河谷平原发病率较低。

土地覆盖动态

南府以镶嵌林地（median=4,203 km2）、阔叶林（3,249 km2）和雨养农田（2,394 km2）为主。2003-2019年间灌丛地减少25%，转为镶嵌林地、雨养农田及阔叶林。针叶林面积自2010年持续增长，城镇区域2015年后趋于稳定。

土地覆盖与发病关联

GAM模型（伪R2=0.54）显示：

• •

  灌丛地（p<0.001）与镶嵌林地（p<0.001）面积与病例数呈正相关

• •

  阔叶林呈现倒U型关系（p<0.001），中等覆盖率时风险最高

• •

  针叶林具有保护效应（p<0.001）

• •

  男性人口数量与发病显著正相关（p<0.001）

土地覆盖变化的影响

580个村庄发生显著土地覆盖变化，其累计发病率高于稳定村庄（中位数2例 vs 1例，p<0.05）。过渡模型（GAM2）显示：

• •

  灌丛地向阔叶林转化呈现倒U型关系（p<0.01），中等转化强度时风险最高

• •

  稳定阔叶林（p=0.08）与灌丛地向草地转化（p=0.28）无显著影响

• •

  空间自相关通过经纬度交互项有效控制

讨论

本研究首次在村级尺度揭示土地覆盖变化对恙虫病的非线性影响：

1. 1.

   景观扰动假说：灌丛地与镶嵌林地作为过渡性生态系统，可能提供适宜媒介螨（如Leptotrombidium deliense）的微环境及宿主（Rattus tanezumi）资源

2. 2.

   森林覆盖悖论：阔叶林中等覆盖率时风险最高，可能与生物多样性阈值及生态网络复杂度有关

3. 3.

   再造林风险：从灌丛向森林的转化过程（而非稳定森林）是关键风险因素，提示生态系统过渡期增加人类暴露风险

   局限性包括遥感数据分辨率限制（无法区分商业种植园）、气候变量缺失及被动监测导致的漏报。建议未来结合媒介-宿主-病原体系统生态学调查，深化景观配置（如边缘效应）对传播动力学的影响机制研究。