译文——2019年中国新型冠状病毒国内外传播的初步风险分析

中国春节伊始,没想到我们就面临着一场严峻的疫情考验,作为一名ISLE团队的研究者,献上一篇相关风险分析的译文。希望为抗击疫情贡献一份微薄之力。

本文原始链接


图片来自WorldPop

1 文章相关信息

翻译原文pdf

Title: Preliminary risk analysis of 2019 novel coronavirus spread within and beyond China

题目: 2019年中国新型冠状病毒国内外传播的初步风险分析。

Author:

  • Shengjie Lai-1
  • Isaac I. Bogoch-2,
  • Alexander Watts-3,4
  • Kamran Khan-2,3,4
  • Zhongjie Li-5
  • Andrew Tatem-1

注: 作者姓名-后的数字代表对应的单位

Affiliation:

  • 1 WorldPop, School of Geography and Environmental Science, University of
    Southampton, UK,
  • 2 Department of Medicine, University of Toronto, Toronto, Canada
  • 3 Li Ka Shing Knowledge Institute, St. Michael’s Hospital, Toronto, Canada
  • 4 Bluedot, Toronto, Canada
  • 5 Division of Infectious Diseases, Chinese Centre for Disease Control and Prevention

联系原作者方式: Shengjie.Lai@soton.ac.uk; A.J.Tatem@soton.ac.uk

2 正文翻译

截止至北京时间2020年1月25日,中国报告确诊2019新型冠状病毒(2019-nCoV)感染确诊病例1409人,疑似病例2032人,死亡病例42人,大部分病例来自湖北省武汉市。几乎所有省份都确诊了来自武汉的输入性病例,同时一些省份还报告了部分二次传播的病例。由于中国正处于长达40天的春节长假期间(从2020年1月10日到2月18日),国内迁徙旅游人数激增——这是世界上规模最大的年度人类迁徙,包括数亿人在整个国家范围内的迁徙,这可能加剧了病毒的传播。

本文使用消除了个人隐私信息并聚合过的2013到2015年人口流动数据,该数据来自百度基于位置信息的服务(LBS服务)数据,以及从国际航空运输组织(IATA)获得的2018年国际航空出行数据来探究从武汉到中国其他城市的迁徙者的流动模式,并分析在春运迁徙期间的2019-nCoV病毒在整个中国以及中国之外传播的风险。

本文使用2013到2015的LBS数据发现大规模的迁徙者在春节前倾向于从武汉前往相邻城市或者其他中国大城市(图1到图3以及表1到表3)。输入性病例的数量与在春节前两周内来自武汉的输入风险具有高度相关性(图4)。此外,大部分外地的输入性病例在疾病暴发初期都出现了病征。尽管从2020年1月23日起,湖北省武汉市和其他一些城市已经实施了封城,但是这已经是在从武汉返乡的春运高峰期后期(图1)。如果在接收了大量武汉迁徙者的省份和城市(如北京,上海和广州)发生二次暴发,这些二次暴发很可能会在7天的全国春节假期结束后通过人口流动进一步传播到中国其他与这些城市有高度联系的城市(图5到图7)。此外,基于航空出行的历史数据,中国高风险城市与其他国家的连接性被定义为春节假期前后的三个月(表4和表5)。本文一开始关注的是非洲的一些特定目的地城市,因为这些区域疫情监测和卫生系统都很薄弱(表6到表7以及图8到图9)。但是后续会将类似的风险评估扩展到世界的其他区域。

鉴于当前流行病学以及对该疫情流行病学的有限理解,我们通过历史数据得到的出行模式的发现可能有助于调整公共卫生干预措施。然而,必须注意的是,我们的分析是基于以往非疫情暴发年份的基础上假设出行模式“与往年相同”进行的迁徙,而目前我们正处于前所未有的状况,整个中国的人群迁徙可能会发生重大变化。我们正在密切监测疫情,并会做进一步的分析,用来估计在春节与接下来的几个月期间,病毒在中国与国际间传播的风险。此外,我们也想尝试评估中国城市交通封锁对于疫情控制的有效性,以及度假归来的人群迁徙对于2019-nCoV病毒传播的影响。

3 附图

3.1 第一部分:国内迁徙

图1 湖北省和武汉市的县6个月内每天人群流动的模式。阴影颜色含义:绿色:春节前两周,红色:春节后两周,紫色:元宵节,粉色:清明节,红线:湖北省城市封城日;相对净流量=(流入量-流出量)/人口,基于2013-2014年百度人口流动数据。

表1 中国大陆春节前两周内接待来自武汉人群最多的前30个城市。

排名 城市 人口 (百万人) 省份 流量 (%)
1 荆州 5.7 湖北 7.26
2 襄樊 5.6 湖北 7.11
3 咸宁 2.5 湖北 6.17
4 北京 21.7 北京 6.10
5 黄冈 6.3 湖北 5.94
6 宜昌 4.1 湖北 5.86
7 黄石 2.5 湖北 5.18
8 孝感 4.9 湖北 4.73
9 省直辖 3.5 湖北 4.62
10 十堰 3.4 湖北 4.14
11 上海 24.2 上海 3.52
12 恩施 3.3 湖北 3.51
13 荆门 2.9 湖北 3.06
14 随州 2.2 湖北 2.49
15 广州 14 广东 2.45
16 郑州 9.6 河南 2.22
17 鄂州 1.1 湖北 1.93
18 天津 15.6 天津 1.56
19 嘉兴 4.6 浙江 1.25
20 杭州 9 浙江 1.23
21 长沙 7.6 湖南 1.13
22 西安 8.3 陕西 1.02
23 南京 8.3 江苏 0.97
24 深圳 10.2 广东 0.96
25 重庆 30.9 重庆 0.82
26 福州 7.6 福建 0.58
27 南昌 5.4 江西 0.57
28 成都 14.3 四川 0.56
29 合肥 7.9 安徽 0.52
30 东莞 8.3 广东 0.45
其他城市 1115.8 12.09
总共 1371.5 100.00
人口:2016年人口,国家统计局。

流量:2014年和2015年春节前两周内离开武汉市的人群百分比。数据来自中国互联网公司百度,专注于互联网相关和基于位置的服务,每天百度能接收到近90亿个位置的请求。

表2 中国大陆春节后两周内接待来自武汉人群最多的前30个城市。

排名 城市 人口 (百万人) 省份 流量 (%)
1 宜昌 4.1 湖北 7.48
2 荆州 5.7 湖北 6.65
3 襄樊 5.6 湖北 6.48
4 黄冈 6.3 湖北 5.91
5 北京 21.7 北京 5.56
6 孝感 4.9 湖北 5.16
7 咸宁 2.5 湖北 4.25
8 省直辖 3.5 湖北 4.22
9 上海 24.2 上海 3.97
10 十堰 3.4 湖北 3.89
11 荆门 2.9 湖北 3.51
12 黄石 2.5 湖北 3.46
13 广州 14 广东湖北 3.07
14 恩施 3.3 湖北 3.01
15 随州 2.2 湖北 2.50
16 鄂州 1.1 湖北 2.26
17 郑州 9.6 河南 2.13
18 长沙 7.6 湖南 1.78
19 天津 15.6 天津 1.65
20 深圳 10.2 广东 1.24
21 西安 8.3 陕西 1.24
22 南京 8.3 江苏 1.13
23 杭州 9 浙江 1.12
24 嘉兴 4.6 浙江 1.04
25 南昌 5.4 江西 0.83
26 重庆 30.9 重庆 0.82
27 福州 7.6 福建 0.82
28 合肥 7.9 安徽 0.78
29 苏州 10.6 江苏 0.51
30 东莞 8.3 广东 0.47
Others 1119.5 13.04
Total 1371.5 100.00
人口:2016年人口,国家统计局。

流量:2014年和2015年自正月初一后离开武汉市的人群百分比。数据来自中国互联网公司百度,专注于互联网相关和基于位置的服务,每天百度能接收到近90亿个位置的请求。

表3 中国大陆春节前后接待来自武汉人群的省份排名。

春节前两周内 春节后两周内
排名 省份 人口 (百万人) 流量 (%) 省份 人口 (百万人) 流量 (%)
1 北京 21.5 16.07 北京 21.5 13.5
2 广东 113.5 12.19 广东 113.5 13.3
3 河南 96.1 9.48 上海 24.2 9.6
4 上海 24.2 9.25 河南 96.1 7.9
5 浙江 57.4 8.19 浙江 57.4 7.2
6 江苏 80.5 5.51 江苏 80.5 6.8
7 湖南 69 4.8 湖南 69 6.2
8 陕西 38.6 4.54 陕西 38.6 4.9
9 天津 15.6 4.11 天津 15.6 4
10 山东 100.5 3.66 山东 100.5 3.8
11 四川 83.4 3.13 福建 39.4 3.7
12 江西 46.5 2.75 安徽 63.2 3.2
13 福建 39.4 2.72 江西 46.5 2.9
14 安徽 63.2 2.62 四川 83.4 2.1
15 重庆 31 2.15 重庆 31 2.0
16 河北 75.6 1.94 河北 75.6 1.7
17 云南 48.3 1.22 辽宁 43.6 1.2
18 广西 49.3 1.10 云南 48.3 1
19 辽宁 43.6 1.06 广西 49.3 1
20 海南 9.3 0.58 山西 37.2 0.6
21 陕西 37.2 0.54 海南 9.3 0.4
22 贵州 36 0.47 贵州 36 0.4
23 黑龙江 37.7 0.4 黑龙江 37.7 0.4
24 新疆 24.9 0.4 新疆 24.9 0.3
25 甘肃 26.4 0.32 吉林 27 0.3
26 吉林 27.0 0.31 甘肃 26.4 0.2
27 内蒙古 25.3 0.29 内蒙古 25.3 0.2
28 宁夏 6.9 0.11 宁夏 6.9 0.1
29 青海 6 0.07 青海 6 0.1
30 西藏 3.4 0.03 西藏 3.4 0
除西藏和青海外,所有省份都报告了输入性或本地确诊病例。

人口:2016年人口,国家统计局。

流量:2014年和2015年春节前两周内以及自正月初一后离开武汉市的人群百分比。数据来自中国互联网公司百度,专注于互联网相关和基于位置的服务,每天百度能接收到近90亿个位置的请求。

图2 中国大陆城市在春节期间接待来自武汉的2019-nCoV感染人群的风险。
根据百度公司的人口流动数据,初步将城市级别的输入风险定义为每个城市在春节前后2周内接待来自武汉的人群在从武汉出发的人群总数中所占的百分比。

图3 中国内地省份在农历新年期间接待来自武汉的2019-nCoV感染人群的风险。
根据百度公司的人口流动数据,初步将省级输入风险定义为每个省在春节前后2周内接待来自武汉的人群在从武汉出发的人群总数中所占的百分比。

图4。输入性病例的时间分布及输入性病例数目与来自武汉人群在春节前前两星期内输入的风险之间的关系。(A)-(C)来自武汉的输入病例的时间分布、发病情况、入院情况和各省诊断情况。(D)各省报告的输入病例数与通过迁徙者输入的风险之间的关系。根据百度公司的人口流动数据,初步将省级输入风险定义为每个省在春节前后2周内接待来自武汉的人群在从武汉出发的人群总数中所占的百分比。截至2020年1月24日,病例数据来自中国国家和地方卫生委员会网站。

图5 北京、上海和广东省的每个县6个月内每天人群流动的模式。阴影颜色含义:绿色:春节前两周,红色:春节后两周,紫色:元宵节,粉色:清明节;相对净流量=(流入量-流出量)/人口,基于2013-2014年百度人口流动数据。

图6。中国大陆城市在春节后的未来四周内接待来自18个高风险城市(蓝色圆圈,有2019-nCoV感染或者输入性病例的城市)的旅客的输入风险。城市级别的输入风险初步定义为接待来自高风险城市的人群在从高风险城市出发的人群总数中所占的百分比。基于从百度获取的人口流动数据计算。18个高风险城市包括武汉和其他省份春节前接待大量的来自武汉人群的17个城市:北京、上海、广州、郑州、天津、杭州、嘉兴、长沙、西安、南京、深圳、重庆、南昌、成都、合肥、福州、东莞。

图7 在春节后的四周内,中国大陆接待来自18个具有2019-ncov感染或输入性病例的高危城市(蓝色圆圈)人群的城市的估计连接性。箭头表示城市级别输入风险的联系与方向,初输入风险。城市级别的输入风险初步定义为接待来自高风险城市的人群在从高风险城市出发的人群总数中所占的百分比。基于从百度获取的人口流动数据计算。18个高风险城市包括武汉和其他省份春节前接待大量的来自武汉人群的17个城市:北京、上海、广州、郑州、天津、杭州、嘉兴、长沙、西安、南京、深圳、重庆、南昌、成都、合肥、福州、东莞。

3.2 第二部分:大陆外出行

表4 在三个月内接待了来自中国大陆18个高风险城市(图6)的航空旅客最多的全球前30位城市,时间为从春节前15天到春节后两个半月。

排名 城市 国家或地区 旅客流量(1000 人) 风险 (%)
1 曼谷 泰国 1062.9 7.86
2 香港 中国香港特别行政区 1001.7 7.41
3 台北 中国台湾 857.5 6.34
4 首尔 韩国 757.9 5.61
5 东京 日本 714.3 5.28
6 新加坡 新加坡 568.1 4.2
7 普吉岛 泰国 492.8 3.65
8 大阪 日本 434.1 3.21
9 吉隆坡 马来西亚 382.7 2.83
10 澳门 中国澳门特别行政区 260.4 1.93
11 巴厘岛 印度尼西亚 222.2 1.64
12 悉尼 澳大利亚 207.4 1.53
13 清迈 泰国 156.9 1.16
14 墨尔本 澳大利亚 154.5 1.14
15 洛杉矶 美国 154.5 1.14
16 纽约 美国 145.9 1.08
17 迪拜 阿联酋 144.9 1.07
18 芽庄 越南 143 1.06
19 伦敦 英国 142.1 1.05
20 胡志明市 越南 142 1.05
21 名古屋 日本 140.1 1.04
22 亚庇 马来西亚 130.4 0.96
23 金边 柬埔寨 127.5 0.94
24 甲米府 泰国 125.2 0.93
25 马尼拉 菲律宾 121.9 0.9
26 暹粒 柬埔寨 121.4 0.9
27 巴黎 法国 119.5 0.88
28 雅加达 印度尼西亚 107.6 0.8
29 高雄 中国台湾 107.6 0.8
30 法兰克福 德国 103.3 0.76
其他城市 4158.2 30.77
总计 13514.9 100
相对风险的初步定义为各个城市接待来自高风险城市(18个城市)航空旅客占高风险城市总出行旅客的百分比,基于从国际航空交通协会(IATA)获得的2018年2月到4月航空出行数据。注,2018年的春节是从2月16日开始。

表5 在三个月内接待了来自中国大陆18个高风险城市(图6)的航空旅客最多的全球前30位国家或地区,时间为从春节前15天到春节后两个半月。

排名 国家或地区 旅客流量(1000 人) 风险 (%)
1 泰国 2031.9 15.03
2 日本 1563.3 11.57
3 中国香港特别行政区 1001.7 7.41
4 中国台湾 979.7 7.25
5 韩国 936.6 6.93
6 美国 773.3 5.72
7 马来西亚 634.3 4.69
8 新加坡 568.1 4.20
9 越南 468.4 3.47
10 澳大利亚 455.6 3.37
11 印度尼西亚 412.5 3.05
12 柬埔寨 262.9 1.95
13 中国澳门特别行政区 260.4 1.93
14 菲律宾 250.3 1.85
15 德国 234.9 1.74
16 加拿大 208.5 1.54
17 英国 190.7 1.41
18 阿联酋 162.3 1.20
19 意大利 152.9 1.13
20 俄罗斯 151.3 1.12
21 法国 137.9 1.02
22 新西兰 120.7 0.89
23 印度 106.7 0.79
24 西班牙 105.8 0.78
25 土耳其 66.5 0.49
26 埃及 57.5 0.43
27 斯里兰卡 55.7 0.41
28 马尔代夫 50.7 0.37
29 荷兰 44.9 0.33
30 缅甸 43.3 0.32
其他国家或地区 1025.6 7.59
总计 13514.9 100
相对风险的初步定义为各个国家或地区接待来自高风险城市(18个城市)航空旅客占高风险城市总出行旅客的百分比,基于从国际航空交通协会(IATA)获得的2018年2月到4月航空出行数据。注,2018年的春节是从2月16日开始。

图8 在春节前15天到春节后两个半月,这三个月时间内全球接待中国大陆18个高风险城市(图6)航空出行游客的城市地理分布。航空出行游客流量最大的30个城市被标注在图上。基于国际航空交通协会(IATA)获得的2018年2月到4月的航空出行数据。2018年春节从2月16日开始。

表6 在三个月内接待了来自中国大陆18个高风险城市(图6)的航空旅客最多的非洲前30位城市,时间为从春节前15天到春节后两个半月。

排名 城市 国家或地区 旅客流量 风险 (%)
1 开罗 埃及 56735 20.49
2 约翰内斯堡 南非 20530 7.42
3 毛里求斯 毛里求斯 18297 6.61
4 亚的斯 埃塞俄比亚 17882 6.46
5 卡萨布兰卡 摩洛哥 15787 5.70
6 内罗毕 肯尼亚 12859 4.64
7 恩德培 乌干达 8246 2.98
8 阿克拉 加纳 8211 2.97
9 拉各斯 尼日利亚 8087 2.92
10 卢萨卡 赞比亚 7672 2.77
11 达累斯萨拉姆 坦桑尼亚 6769 2.44
12 阿尔及利亚 阿尔及利亚 6074 2.19
13 罗安达 安哥拉 5994 2.16
14 喀土穆 苏丹 5412 1.95
15 阿布贾 尼日利亚 4193 1.51
16 卢本巴沙 刚果(金沙萨) 3546 1.28
17 阿比让 科特迪瓦 3461 1.25
18 开普敦 南非 3455 1.25
19 科纳克里 几内亚 3455 1.25
20 突尼斯 突尼斯 2912 1.05
21 利伯维尔 加蓬 2786 1.01
22 哈拉雷 津巴布韦 2665 0.96
23 达喀尔 塞内加尔 2659 0.96
24 马普托 莫桑比克 2560 0.92
25 塔那那利佛 马达加斯加 2515 0.91
26 努瓦克肖特 毛里塔尼亚 1955 0.71
27 马拉博赤道 几内亚 1864 0.67
28 马赫岛 塞舌尔 1850 0.67
29 德班 南非 1815 0.66
30 恩多拉 赞比亚 1796 0.65
相对风险的初步定义为每个非洲城市接待来自高风险城市(18个城市)航空旅客占高风险城市总出行旅客的百分比,基于从国际航空交通协会(IATA)获得的2018年2月到4月航空出行数据。注,2018年的春节是从2月16日开始。

表7 在三个月内接待了来自中国大陆18个高风险城市(图6)的航空旅客的非洲国家或地区,时间为从春节前15天到春节后两个半月。

排名 国家或地区 旅客流量 % 排名 国家或地区 旅客流量 %
1 埃及 57516 20.77 27 塞舌尔 1863 0.67
2 南非 26405 9.54 28 博茨瓦纳 1627 0.59
3 埃塞俄比亚 18393 6.64 29 吉布提 1602 0.58
4 毛里求斯 18297 6.61 30 马里 1587 0.57
5 摩洛哥 16974 6.13 31 刚果(布拉柴维尔) 1500 0.54
6 尼日利亚 13734 4.96 32 乍得 1425 0.51
7 肯尼亚 13185 4.76 33 卢旺达 1386 0.5
8 赞比亚 9471 3.42 34 塞拉利昂 1330 0.48
9 坦桑尼亚 8388 3.03 35 纳米比亚 1207 0.44
10 乌干达 8246 2.98 36 马拉维 1139 0.41
11 加纳 8211 2.97 37 贝宁 890 0.32
12 阿尔及利亚 7887 2.85 38 多哥 858 0.31
13 安哥拉 5994 2.16 39 莱索托 853 0.31
14 苏丹 5433 1.96 40 留尼汪 809 0.29
15 刚果(金沙萨) 5248 1.90 41 尼日尔 790 0.29
16 莫桑比克 3928 1.42 42 利比里亚 711 0.26
17 科特迪瓦 3511 1.27 43 南苏丹 711 0.26
18 几内亚 3455 1.25 44 布基纳法索 406 0.15
19 突尼斯 2912 1.05 45 冈比亚 365 0.13
20 加蓬 2786 1.01 46 中非共和国 339 0.12
21 喀麦隆 2734 0.99 47 佛得角 276 0.10
22 津巴布韦 2716 0.98 48 厄立特里亚 246 0.09
23 塞内加尔 2659 0.96 49 布隆迪 232 0.08
24 马达加斯加 2515 0.91 50 科摩罗 178 0.06
25 毛里塔尼亚 1955 0.71 51 索马里 68 0.02
26 赤道几内亚 1864 0.67 52 几内亚比绍 52 0.02
相对风险的初步定义为每个非洲国家或地区接待来自高风险城市(18个城市)航空旅客占高风险城市总出行旅客的百分比,基于从国际航空交通协会(IATA)获得的2018年2月到4月航空出行数据。注,2018年的春节是从2月16日开始。

图9 在春节前15天到春节后两个半月,这三个月时间内非洲接待中国大陆18个高风险城市(图6)航空出行游客的城市地理分布。基于国际航空交通协会(IATA)获得的2018年2月到4月的航空出行数据。2018年春节从2月16日开始。

4 参考文献

  1. Zhu N, et al. (2020) A Novel Coronavirus from Patients with Pneumonia in China, 2019. N Engl J Med. DOI: 10.1056/NEJMoa2001017
  2. Chan J, et al. (2020) A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. Lancet. DOI: 10.1016/S0140-6736(20)30154-9
  3. Huang C, et al. (2020) Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. DOI: 10.1016/S0140-6736(20)30183-5
  4. Baidu Migration. http://qianxi.baidu.com/
  5. Bogoch I, et al. (2020) Pneumonia of Unknown Etiology in Wuhan, China: Potential for International Spread Via Commercial Air Travel. Journal of Travel Medicine, taaa008, https://doi.org/10.1093/jtm/taaa008

5 后记

本译文已授权一览众山小·可持续城市与交通发布,详情链接见抗击疫情特刊⑥ ︱2019年中国新型冠状病毒国内外传播的初步风险分析,此外可以参加下面相关抗击疫情特刊。前两篇笔者均有参与翻译。另外还有一篇我们ISLE团队发布的社区倡议书。

抗击疫情特刊① ︱武汉新型冠状病毒肺炎通过航空国际传播的风险评估

抗击疫情特刊② ︱Science:传染病是这么传染的!传染病者的“有效”距离

抗击疫情特刊③︱bioRxiv:导致近期人病毒性肺炎暴发的新型冠状病毒及其源于蝙蝠的相关研究

抗击疫情特刊④ ︱面对新型肺炎来袭你应该怎么做?听听美国疾控中心的建议

抗击疫情特刊⑤ ︱从2003年SARS到2009年甲型H1N1流感:台湾预防流行病的经验

ISLE团队的社区倡议书。

倡议书 | 面对新型冠状病毒肺炎:社区卫生服务中心(家庭医生/基层党群)行动起来!