[논문 리뷰] Natural human mobility patterns and spatial spread of infectious diseases

Vitaly Belik, T. Geisel외 1명|DSpace@MIT (Massachusetts Institute of Technology)|2011. 03. 31.

COVID-19 epidemiological studies인용 수 33

한 줄 요약

이 논문은 실생활의 이동 패턴을 반영하는 양방향 이동을 명시적으로 포함한 공간적 전염병 확산을 위한 새로운 메타집단 모델을 제안한다. 기존의 반응-확산 모델이나 감염력 모델과 달리, 이동 빈도가 증가함에 따라 유행 전파 속도가 포화 상태에 도달함을 예측하고, 전 세계 유행의 새로운 임계 조건을 도출한다. 이 임계 조건은 총 이동 유량이 아닌 환승률에 따라 결정되며, 고전적 유행 모델링의 가정에 도전한다.

ABSTRACT

We investigate a model for spatial epidemics explicitly taking into account bi-directional movements between base and destination locations on individual mobility networks. We provide a systematic analysis of generic dynamical features of the model on regular and complex metapopulation network topologies and show that significant dynamical differences exist to ordinary reaction-diffusion and effective force of infection models. On a lattice we calculate an expression for the velocity of the propagating epidemic front and find that in contrast to the diffusive systems, our model predicts a saturation of the velocity with increasing traveling rate. Furthermore, we show that a fully stochastic system exhibits a novel threshold for attack ratio of an outbreak absent in diffusion and force of infection models. These insights not only capture natural features of human mobility relevant for the geographical epidemic spread, they may serve as a starting point for modeling important dynamical processes in human and animal epidemiology, population ecology, biology and evolution.

연구 동기 및 목표

• 기존의 반응-확산 모델과 감염력 모델이 실생활 인간 이동 패턴을 충분히 반영하지 못하는 데서 비롯된 한계를 해결하기 위해.
• 거주지와 목적지 사이의 공간적으로 제약된 이중 이동이 유행 역학에 어떤 영향을 미치는지 조사하기 위해.
• 개별 이동 네트워크가 전파 속도와 유행 임계 조건과 같은 주요 유행 특성에 어떻게 영향을 미치는지 파악하기 위해.
• 정규 격자와 복잡한 메타집단 네트워크에서 제안된 모델을 기존 표준 모델과 비교하기 위해.
• 이중 이동에 기반한 유니버설 역학적 특성들이 네트워크 구조에 관계없이 공통적으로 나타나는지 규명하기 위해.

제안 방법

• 모델은 각 집단에 대해 병원체-감염자-면역(SIR) 프레임워크를 사용하며, 기초(거주지)와 목적지 사이의 이중 이동을 명시적으로 고려한다.
• 개별 이동 네트워크를 명시적으로 모델링하여, 각 개인이 주로 거주지와 직장 사이를 이동하는 구조를 반영한다.
• 정규 격자와 복잡한 메타집단 네트워크(예: 스케일프리 및 에르되시-레니 네트워크)에서 시스템을 분석하며, 위치 간 유량을 정의한다.
• 주요 매개변수로는 위치 간 이동률 $\omega$, 환승률 $\omega^{-}$, 회복률 $\beta$ 가 있으며, $\omega = 2\omega^{+}\omega^{-}/(2\omega^{+}+\omega^{-})$ 로 정의된다.
• 전역 공격 비율을 분석하고 침입 임계 조건을 규명하기 위해 스토케스틱 시뮬레이션 프레임워크를 사용한다.
• 이론적 분석을 통해 네트워크 구조와 함께 감염자 수 $\lambda \approx N\omega / (\beta + \omega^{-})$ 를 기반으로 한 유행 전파 임계 조건을 유도한다.

실험 결과

연구 질문

• RQ1 거주지와 목적지 사이의 이중 이동이 유행 전파의 전파 속도에 어떤 영향을 미치는가?
• RQ2 공간적으로 제약된 이동 네트워크 존재가 전 세계 유행 전파에 대한 새로운 임계 조건을 초래하는가?
• RQ3 이 모델의 행동은 기존의 반응-확산 모델과 감염력 모델과 비교해 유행 역학 측면에서 어떻게 다를까?
• RQ4 환승률 $\omega^{-}$ 는 유행이 전 세계로 퍼지기 위해 필수적인 역할을 하는가?
• RQ5 관찰된 임계 조건 행동은 라티스와 복잡한 네트워크와 같은 다양한 네트워크 구조에서 공통적으로 나타나는가?

주요 결과

• 이동률이 증가함에 따라 유행 전파 속도가 포화 상태에 도달함을 예측하며, 이는 확산 모델에서 관찰되는 무한 증가와 대조된다.
• 전 세계 유행에 대한 새로운 임계 조건이 존재하며, 이는 총 이동 유량 $\omega$ 가 아닌 환승률 $\omega^{-}$ 에 따라 결정되며, 반응-확산 및 감염력 모델에서는 관찰되지 않는다.
• 낮은 환승률($\omega^{-}$)에서는 공격 비율이 1에 가까워지며, 높은 환승률에서는 급격히 감소하고 사라지며, 이는 전 세계 확산 실패를 의미한다.
• 일차원 격자에서 유행 전파의 임계 조건은 $\frac{\omega^{-}}{\beta} \lesssim 2N(\mathcal{R}_{0}-1)\frac{\omega^{+}}{\beta} - 1$ 으로 유도되며, 환승률과 지역 인구 수에 의존함을 보여준다.
• 다양한 인구 규모에서 공격 비율이 $\rho \sim \rho(\omega^{-}/N)$ 에 따라 유니버설 곡선으로 수렴함을 확인하여, 네트워크 구조와 무관한 유니버설 메커니즘이 있음을 시사한다.
• 한 지점당 무한히 많은 개체가 존재하는 극한에서 임계 효과는 사라지며, 이는 임계 효과가 유한한 인구 규모에 의존함을 확인한다.