Novel travel time aware metapopulation models: A combination with multi-layer waning immunity to assess late-phase epidemic and endemic scenarios

Henrik Zunker, René Schmieding, David Kerkmann, Alain Schengen, Sophie Diexer, Rafael Mikolajczyk, Michael Meyer-Hermann, Martin J. Kühn
doi: https://doi.org/10.1101/2024.03.01.24303602

Abstract

In the realm of infectious disease control, accurate modeling of the transmission dynamics is pivotal. As human mobility and commuting patterns are key components of communicable disease spread, we introduce a novel travel time aware metapopulation model. Our model aims to enhance estimations of disease transmission. By providing more reliable assessments on the efficacy of interventions, curtailing personal rights or human mobility behavior through interventions can be minimized. The proposed model is an advancement over traditional compartmental models, integrating explicit transmission on travel and commute, a factor available in agent-based models but often neglected with metapopulation models.

Our approach employs a multi-edge graph ODE-based (Graph-ODE) model, which represents the intricate interplay between mobility and disease spread. This granular modeling is particularly important when assessing the dynamics in densely connected urban areas or when heterogeneous structures across entire countries have to be assessed. The given approach can be coupled with any kind of ODE-based model.

In addition, we propose a novel multi-layer waning immunity model that integrates waning of different paces for protection against mild and severe courses of the disease. As this is of particular interest for late-phase epidemic or endemic scenarios, we consider the late-phase of SARS-CoV-2 in Germany.

The results of this work show that accounting for resolved mobility significantly influences the pattern of outbreaks. The improved model provides a refined tool for predicting outbreak trajectories and evaluating intervention strategies in relation to mobility by allowing us to assess the transmission that result on traveling. The insights derived from this model can serve as a basis for decisions on the implementation or suspension of interventions, such as mandatory masks on public transportation. Eventually, our model contributes to maintaining mobility as a social good while reducing exuberant disease dynamics potentially driven by travel activities.

Author summary As human contacts and contact networks are key to the development and prediction of infectious disease spread, travel and commuting activities are important components to be considered in mathematical-epidemiological modeling. Two, often contrasting modeling approaches, based on subpopulations and based on individuals can provide insights of different granularity but also come at different levels of complexity.

With this article, we extend a recently introduced Graph-ODE-based model by the explicit introduction of mobility-based infection models in which we allow focused nonpharmaceutical interventions, like face mask mandates in public transport, and in which we can explicitly keep track of secondary cases induced by travel activities, a component mostly not available with equation-based models. In addition, we introduce a novel multi-layer waning immunity model particularly suitable for late-phase epidemic or endemic scenarios.

On a daily level and geographically small scale, the newly proposed model often develops similarly, although our results show that complex mobility networks can lead to substantially different disease dynamics in the entirety of a federal state or country. The proposed model thus enables a better understanding of infectious disease dynamics through mobility. It allows for targeted numerical investigations and thus leads to more appropriate real-world interventions.

Introduction

According to the World Health Organization (WHO), the COVID-19 pandemic resulted in 14.9 million excess deaths in 2020 and 2021 and, although fighting infectious diseases, WHO predictions foresee communicable diseases to still account for 6 % of global deaths by 2048 [1]. The unprecedented development and administration of COVID-19 vaccines has undoubtedly saved many lives and years of live [2]. Past experiences have illustrated that waning immunity [2] and reduced vaccine protection against simple transmission (i.e., any infection) may require continuous pairing of vaccination with nonpharmaceutical interventions (NPIs) such as face masks or physical distancing, in particular in situations of high transmission (risk) [3]. For instance, after several vaccination campaigns, most SARS-CoV-2 related interventions had been discontinued in Germany until summer or early autumn 2022. To further protect vulnerable groups, FFP2 masks became mandatory in all long-distance trains on October 1st, 2022 and federal states could implement additional mask mandates in local or regional public transport [4].

In order to proactively react to infectious disease propagation, mathematical models are of great aid. Models of different types have been largely used to guide policy-makers towards evidence-based decisions, see, e.g. [5] for autumn and winter scenarios in Germany in 2022/2023. Over the last years, to predict SARS-CoV-2 development in Germany, contributions have been made by a variety of different approaches. Among these approaches are agent-based models [6, 7], models based on ordinary differential equations (ODE) [8–10] or delay differential equations [11], advanced ODE-based models [12] using the linear chain trick [13] to waive the implicit assumption of exponentially distributed transition times, and in particular metapopulation or graph-ODE models [14].

Mathematical models based on systems of ordinary differential equations (ODE), often denoted compartmental models, are popular tools used in the context of modeling infectious diseases, see, e.g., [15]. This is due to their well-understood character and established methods for model analysis but also because of the low computational cost. However, transmission dynamics of communicable diseases naturally follow the complex network structures of human mobility. Standard ODE-based metapopulation (see, e.g., [16]) or Graph-ODE models [17] can be seen as a simple compromise between too simple ODE-based and complex agent-based models. They combine the low computational effort of ODE-based modeling on small- or medium-sized regional entities with realistic mobility networks between these regions.

Although we see certain advantages of Graph-ODE models [17] over classic ODE-based metapopulation approaches, the model in [17] neglected travel times and explicit transmission during commuting. In the current paper, we propose an extended travel time-aware Graph-ODE model that allows for explicit modeling of travel time and transmission on travel and commute.

Aside from mobility, disease dynamics can be driven or mitigated by non-protection or protection against a given pathogen. In particular for late-phase epidemic or endemic scenarios, waning immunity is an important factor to be considered. However, a recent review study [18] on models used in different forecast and scenarios hubs showed that only four out of 90 models provided information or modeling on waning immunity. We thus propose a multi-layer waning immunity model based on ordinary differential equations. In this model, we can consider three different subpopulations with corresponding immunity layers and protection factors. Furthermore, protection against mild and severe courses of the disease wane with different paces.

The current paper is organized as follows. In the first part of materials and methods, we introduce a simple ODE-SIR model and extend it by a mobility-model. We then extend a recently proposed hybrid Graph-ODE ansatz to use our novel travel-time aware model locally. Then, we additionally introduce an advanced ODE-SECIRS model containing three different immunity layers and two different paces for waning immunity, against any and severe infection, respectively. Then, we provide parameters for mobility, contact patterns, and transmission dynamics. In the results section, we show how the new model advances the previous one. Eventually, we discuss advantages and limitations and draw a conclusion.

Materials and methods

In this section, we will introduce a novel travel time aware metapopulation model which will be implemented in a multi-edge graph based on [17] to account for spatially heterogeneous disease dynamics. For a county-level resolution of Germany, we provide travel and commute patterns based on [19]. Furthermore, we introduce a multi-layer waning immunity model based on ordinary differential equations (ODEs) and which is of particular interest in late-phase epidemic or endemic scenarios. Finally, we will summarize parameters used in our estimations for the late-phase SARS-CoV-2 demonstrator case.

A novel travel time aware metapopulation model

For compartment-based modeling, we need to define a list of disease or infection states. In the simplest SIR case, we have the states Susceptible for persons who can get infected, Infected or Infectious for persons who can infect others, and Recovered or Removed for individuals that cannot get reinfected. In order to introduce our travel time aware metapopulation model, we start from this simple ODE-SIR model which writes Embedded Image and where ϕ(k)(t) represents the (mean) daily contacts of a person at time t, ρ(k)(t) the (mean) transmission risk, Embedded Image the (mean) time a person stays infected, and N (k) = N (k)(t) = S(k)(t) + I(k)(t) + R(k)(t) the (constant) total population size. Note that in these model, we always use mean or averaged values so that we drop the mean or average description to parameters in the following.

For the sake of a consistent development of the new model, we have already introduced the index k that will later be used to refer to a particular region k but which is without meaning for the presentation in this subsection.

We suppose some initial conditions at t0 to be given by (S(k)(t0), I(k)(t0), R(k)(t0)). In the following paragraphs, we will introduce our mobility-model approach for ODE-based models. In the end, the model (1) can be replaced by any kind of complex equation-based model.

In infectious disease modeling, the number of daily contacts is often estimated from surveys or contact diaries like [20] and we assume that we can estimate the number of contacts in traffic modes Embedded Image.

Classic traffic contacts integration and traffic-related NPIs

A straightforward approach modeling (local) traffic-related contacts in compartmental models is given with a separation of contacts in traffic and nontraffic. We can then write Embedded Image where Embedded Image represents all contacts that happen elsewhere, i.e., in home, school, workplace, or leisure activity settings. If Eq (2) is plugged into Eq (1), we can model a reduction in travel contacts from day t1 > t0 with factor r ∈ [0, 1] by using Embedded Image where we define δ as a transition interval and Embedded Image as a transition contact rate between the change of the contact rate to ensure that Embedded Image with 0 < ϵ ≪ 1 for Eq (1) to be well defined.

This type of contact reduction for contact locations such as home, school, work, and other has been considered in [14, 17, 21]. For noncontinuously differentiable parameters, we actually compute the solution of a new initial value problems (IVP) from t1.

Explicit traffic or mobility modeling

In order to separate disease dynamics on travel and commute from nontraffic (i.e., home or work) dynamics, we propose a new way by introducing a mobility-based infection model in which dynamics follow the same mechanisms but in which the population distribution as well as transmission or mitigation parameters differ; see Fig. 1. In the following and for the sake of simplicity, we will refer to the mobility-based infection model simply as mobility model or traffic model – although the model is a pure infection dynamics model, just parameterized for mobility settings.

Fig 1.
Fig 1. Compartmental transmission model (white) complemented by compartmental transmission mobility model (gray).

Note that the novel model is indeed different and expected to behave different where it relaxes assumptions such as homogeneous mixture. However, for meaningful comparisons and in order to be consistent, contact patterns across both models have to be adjusted. In the following, we will explain how to scale contacts that are obtained from classical diary studies to our new model.

Without loss of generality, let us assume that all individuals have contacts in nontraffic locations (which is trivial, i.e., in home or work) while only a portion Embedded Image also has contacts in traffic locations (which is nontrivial as people might not commute at all or only use individual traffic where there is no contact with others). For demonstration purposes, we consider the classic and the new model on a timescale of one day, where each day has unit length, and neglect the daily return trip for now. In the final implementation, this approach is then just down-scaled to half-a-day.

As we are still working with a mean-value model, we assume that mobility completely happens either at the beginning or the end of the day, respectively. We will provide graphs for both realizations, but provide the following description with mobility taking place at the end of the day. To that end, we introduce the daily travel time Embedded Image, such that the mobility process starts near the end of the day. Note that we will later use a time scale of half a day, with a length of stay in another region and return trip to the home region, once we introduced a second spatial region. For a visualization on the final implementation with round trips and stay times in other nodes, see Fig. 4 (right).

We further assume that the number of contacts in nontraffic locations, Embedded Image, is independent of the additional contacts during travel. In order to compute the number of contacts in traffic modes per individual traveler, we compute Embedded Image

The end-of-day approximation of our novel model is then given after solving three initial value problems. We solve Embedded Image with initial values (S0, I0, R0) := (S(k)(0), I(k)(0), R(k)(0)). The scaling of contacts comes from the fact that Embedded Image many individuals now have all their nontraffic related contacts within Embedded Image many days, Embedded Image. This means that the value Embedded Image, given on daily scale, needs to be rescaled. Let Embedded Image be the solution of (5) at Embedded Image. We define the IVP Embedded Image with initial values Embedded Image and a total population of Embedded Image, i.e., the people that do not travel, in the denominator. Furthermore, we define the IVP in the travel node by Embedded Image with initial values Embedded Image and a total population of Embedded Image in the denominator. The complete solution at t = 1 is then given by the summed solutions of (6) and (7). Geometrically, IVPs (5) and (6) are assigned to the white local model from Fig. 1, while (7) corresponds to the grey traffic or mobility model.

Here, with Embedded Image and Embedded Image, we used an equal and potentially unrealistic distribution of the population compartments to traffic and nontraffic node, respectively. In a more complex model, one would also reduce the relative share of infected commuters, letting Embedded Image susceptible but only Embedded Image, infected individuals commute. This reduction has already been used in [17] and naturally translates to our new model.

In order to give the motivation for our scaling, we assume constant contact patterns over the considered day, i.e., Embedded Image and Embedded Image. For the total number of contacts, we then directly obtain Embedded Image where we recognize Eq (2).

In Fig. 2, we provide two examples of our new model for the parameters Embedded Image and ρ(t) = 0.05 (left) and ρ(t) = 0.1 (right). In the new model, even if the total number of contacts stays the same, the contact structure changes from an average of ϕ to a bimodal contact distribution. We see that, both, the mobility-first and mobility-last approach end up with roughly the same number of infections at the end of the day but with, of course, slightly different numbers of susceptible and recovered individuals. We demonstrate the application of the scaling procedure to real contact data in the section Mobility and contact patterns

Fig 2.
Fig 2. Infection curves for simple (blue) and mobility-based infection model (orange and green) approach.

The orange curve shows the dynamics with mobility first while the green curve shows the mobility last approach. The example’s parameters are Embedded Image and ρ(t) = 0.05 (left) and ρ(t) = 0.1 (right). In the legend of the plots, we abbreviate nontraffic with NT and traffic with T.

Mobility complemented metapopulation models through a multi-edge graph

In order to make use of the novel model from the previous section, assume n different geographic units (here, denote regions) to be given. We now combine the mobility model extension of a simple ODE model with the graph approach proposed in [17]. For the paper to be self-contained, we visualize the previous approach in Fig. 3 (left). There, each region is represented by the node of a graph while the (multi-)edge ℰ ij between nodes 𝒩 i and 𝒩 j represents the (outgoing) mobility and mobility-based exchange is realized instantly twice a day (round trip). In practice, each combination of sociodemographic group gl, l = 1, …, G and infection state zl, l = 1, …, Z can be assigned a number of travelers such that the multi-edge consists of G × Z single edges. Return trips are realized by a mapping on the same edge, as the vector of weights on ℰij represents the number of outgoing travelers. For the sake of a simple visualization in Fig. 3 (left), each region was assigned an ODE-SIS model with only two compartments (S and I) and two (sociodemographic) groups (e.g., age groups).

Fig 3.
Fig 3. Hybrid graph-ODE approach (left) and extension with mobility models (right).

For simplification, we visualized an SIS model with two sociodemographic groups (e.g. age groups). In the old approach (left), no travel time was used while the novel approach models (right) travel time through compartmental mobility models.

Fig 4.
Fig 4. Timing of the mobility activities.

The hybrid Graph-ODE approach from [17] (left) synchronizes population exchanges at two common daily junctures. The new approach (right) defines an individual structure of mobility activities for each commuter group, which is considering the travel time and the length of stay in the destination. Within this daily schedule, all commuters still return to their home location at the same time, exactly at the end of the day.

This approach extends classic ODE-based models but simplifies the mobility pattern by instantaneously transferring commuters to their intended destinations, omitting the duration of transit. Upon reaching their destinations, individuals remain there for the day’s remainder before they immediately return to their origin. In Fig. 4 (left), we provide a schematic representation of the mobility distribution under this mobility scheme. This assumption offers the advantage of facilitating the coordination and incorporation of commuters into and out of various regions by standardizing the timing of exchanges. It affords flexibility in selecting larger time steps for the applied numerical integration scheme. However, it also abstracts away critical variables such as the travel time. The instantaneous nature of exchanges precludes the possibility of transmission events during transit, an aspect that can become significant in disease spread modeling.

Deviating from this approach, our novel approach introduces a refined segmentation of the mobility activity schedule by delineating distinct mobility profiles based on local travel and stay times; see Fig. 4 (right). The length of stay can be chosen differently for each region. An even more refined selection, e.g., an individual selection for each edge (i.e., a subpopulation of a predefined socioecomic group with a given infection state), would make the scaling of the contacts considerably more complex. This granular approach not only enhances the precision of our method but also mirrors the complex, real-world patterns of human movement more closely. However, including these details means that mobility activities are much more spread out throughout the daily schedule. Therefore, we have to choose substantially smaller time steps in the numerical integration scheme in order to take into account outgoing or incoming commuters.

In the novel approach, traveling is realized via mobility or traffic models; see Fig. 3 (right). Again, each region gets assigned a local compartmental model, now complemented by a local compartment mobility model.

A key advantage of augmenting the framework with local compartment mobility models lies in the creation of realistic travel chains. The idea is that commuters pass through other regions during the trip and that the individuals spend time in the mobility models of the region they are passing through. This temporary inclusion enables potential interactions with other commuter groups that are simultaneously represented in the same mobility model.

The trip chain 𝒯 (j,k) between two arbitrary nodes 𝒩j and 𝒩k is defined as an ordered tuple of regions, Embedded Image where j is the index of the starting node, k is the index of the destination node. The transit nodes Embedded Image describing the way of travel between these regions.

The construction of these travel chains is achieved by calculating the centroids of both the origin and destination regions, represented by (multi-)polygons. A line is then drawn between these centroids; see Fig. 5. The sequence in which these regions are encountered along the line is crucial, as it establishes the order of the travel itinerary. In the final step, the overall travel duration is allocated proportionally across the regions in the defined trip chain.

Fig 5.
Fig 5. Determination of the path chains.

Given several regions, we draw a line between the start and destination region (left). Afterwards we consider all regions that are hit by this line and determine the path chain based on the order (right). Using geodata “Verwaltungsgebiete 1:2 500 000, Stand 01.01. (VG2500)” from https://gdz.bkg.bund.de, copyright; GeoBasis-DE / BKG 2021, license dl-de/by-2–0, see https://www.govdata.de/dl-de/by-2-0.

Given that the visualized trip from Region 𝒩i to Region 𝒩k in Fig. 3 (right) is a direct trip, a traveler spends half of the travel time in each of the corresponding two mobility models. For a trip from Region 𝒩i to Region 𝒩j, dissected to a trip chain 𝒯 (i,j) = (𝒩i, 𝒩k, 𝒩j), a third of the trip is spent in the single mobility models. This forms our extended hybrid graph-ODE model. The arrows in Fig. 3 (right) are representing mobility between the different mobility models.

Finally, we need to determine who is commuting at all. We would like to note here that different travel restrictions can be active for each region, but there should also be the possibility of restrictions on individual edges, e.g., for the subpopulations of detected or symptomatic individuals. We determine the number of commuters in compartment z ∈ Ƶ leaving region 𝒩j for region 𝒩k by Embedded Image where Ƶ represents the compound of all compartments of the considered local model. Furthermore, d(j,k) denotes a local reduction factor in daily commuting activity which can be defined for each edge, i.e., any combination of sociodemographic group and infection state. Region-based restrictions, e.g., isolation of particular infection states or general stay home restrictions, can be chosen with Embedded Image. For example, we can use Embedded Image to isolate symptomatic individuals and to restrict their commuting activities. Finally, we have the total number of people moving between nodes 𝒩j and 𝒩k, denoted by C(j,k). In summary, we break the total given number of people moving from one region to another, into subpopulations of different infection states for which different restrictions can be applied. These numbers are time dependent.

A important challenge of the approach is tracing the compartmental affiliation of the traveling individuals. We address this issue by employing and extending the method described in [17], where we perform an approximation step using a single-step integration method. With this step, we solely advance the infection states of a focus group of individuals (i.e., travelers) while considering all other persons within the same model (travelers from other destinations and local population) as contact persons only. This approach is independent of the model selected and offers some kind of subpopulation-tracing in a population-level model. As the local travel time is generally small, we also apply a single step method for traffic nodes where only one subpopulation is present. Fig. 6 demonstrates the movement patterns of individual groups and how interactions between different groups can lead to additional infections during a trip. Please note that visualization is based on individuals while in practice only subpopulation shares are known.

Fig 6.
Fig 6. Demonstration of the mobility across three nodes and four edges.

The travel time between region i and k, Embedded Image, is two times the travel time from region i to j, i.e., Embedded Image. For simplification, there is no exchange between nodes k and j. The images read from top left to bottom right and all subpictures visualize the developments between (and at the end of) the considered time frame, i.e., t0 → t1 for top left. At time t0, for all graph nodes, travelers change from their local infection model to their local traffic infection model. During the interval [t0, t1] (top left), state changes in traffic infection models are approximated by a single step time integration method while any other (generally an adaptive high-precision scheme like Cash Karp 5(4) [22]) method can be used in the local nodes. Subpopulations that come from or go to a different location are only considered as contact populations; see [17]. This is the first part of the functionality implemented on the graph’s (multi) edges. At t1, the exchange between nodes i and j happens. By t1 (top right), the exchange along the edge between nodes i and j has already happened and three infected plus one recovered individual have moved from travel model of node 𝒩j to travel model of node 𝒩i and vice versa. In [t1, t2), the individuals present in travel node 𝒩i and subject to move to k get in contact with the individuals arriving from node 𝒩j. After another time step until t2, the travelers between node 𝒩i and node 𝒩k are exchanged. From t2 to t3 (bottom left), we further advance the newly exchanged travelers to finally arrive in the local models at t3, i.e., after they have completed their predefined travel time within the mobility model. For the local models in nodes i and j, we must update the states of the arrived individuals using the introduced method, as both local residents and travelers are present, and we will integrate new individuals to the region in the next step. By t3 + ε, ε > 0 (bottom right), all individuals have reached their destination regions, i.e. the travel models are empty, and local models can be advanced in parallel.

A multi-layer waning immunity model of SECIRS-type

Waning immunity is a paramount feature in epidemiological modeling if late-phase epidemic or endemic scenarios are considered and immunity of individuals is not everlasting. This, for instance, is the case for SARS-CoV-2 or influenza, where neither infection nor vaccination promises lasting protection against (any) re-infection; see, e.g., [23].

In the following, we present a SECIRS-type model with waning immunity; see Fig. 7. The motivation for the particular model is threefold and not only in waning immunity. First, we use a differentiation in nonsymptomatic and symptomatic transmission to separate their related contributions to the disease dynamics. Second, we use three different layers of immunity or subpopulations with different protection factors against mild and severe courses of the disease. Finally, we introduce the compartments of temporary immunity to realize different paces of immune waning against mild and severe infections, respectively. For instance, for SARS-CoV-2, this was observed in [2]. In the consequence, our model realizes quickly waning protection effects against transmission or any infection while it preserves longer-lasting immunity against severe or critical infections. To account for age-specific transmission and infection parameters, our model will also be stratified for age groups. Our open-source implementation in MEmilio [24] allows for a variable number of age groups as well as for flexible integration of other sociodemographic factors such as income.

Fig 7.
Fig 7. SECIRS-type model with high temporary immunity upon immunization and waning immunity over time.

The three different immunity layers are visualized by different colors. The waning is represented by orange arrows. Vaccinations are visualized through dashed green arrows and reduce susceptibility for transmission for a certain (short) period of time. Recovery happens from any kind of infection state to a temporary immunity state with normal black arrow. We have highlighted the disease states responsible for transmission by red boxes. For the sake of simplification, we do not visualize the age group model here.

Our model is based on the previously developed SECIR-type model in [14], which allowed the integration of vaccinations and different immunity layers. However, the introduction of the temporary immunity states changes the meaning of recovery. Recovery always means to enter a temporary immunity compartment and then a susceptible compartment with different protection layer. Furthermore, the states enable the simulation of waning immunity with different paces to suitably model endemic or late-phase epidemic scenarios. The basic infection process of our model follows an susceptible (S), exposed (E), nonsymptomatic infectious (INS or C for carrier), and, potentially, symptomatic infectious (ISy), infected severe (ISev), infected critical (ICr), and Dead (D) structure.

We refer to the first subpopulation (shown in gray in Fig. 7) as the naïve subpopulation. Individuals of this group have either not seen vaccination or infection, did not have a substantial immune reaction to these or immunity has waned completely, with the last infection or vaccination event long ago. Individuals with partial immunity (shown in purple) have an average protection, while individuals with improved immunity (shown in blue) are best protected against the pathogen. The arrows in Fig. 7 show a permeable model where individuals can change the subpopulation during simulation.

The different states of the model and their meanings, independent of the particular immunity, are shown in Table 1. To further motivate our model development, we have sketched out the idea of different protection layers and waning immunity paces on the individual layer in Fig. 8.

View this table:
Table 1. Model compartments and meanings for disease transmission modeling.
Fig 8.
Fig 8. Motivation on different immunity layers and waning immunity waning paces on the individual layer.

Recent infections or vaccinations protect against transmission or any infection for a short time. Less-recent immune boosters may not adequately protect against transmission but may result in good protection against a severe course of the disease.

In order to cover age group-specific characteristics of the viral dynamics, we stratify the population by different groups i ∈ {1, …, G}, where G is the number of age groups. For all parameters and variables, we add the subscript i to distinguish the different age groups.

In order to introduce the model equations, we define Embedded Image as people from age group i ∈ {1, …, G} which are not in the dead compartment. Furthermore, denote by Embedded Image, the j-th disease state and by Embedded Image the worsened disease state, e.g., Embedded Image: Infectious, No Symptoms, Naive (age group i) and Embedded Image: Infectious, Symptoms, Naive (age group i). We order the 18 disease states ƵI as follows: Exposed Naive Not Infectious, …, Dead Naive, Exposed Partial Immunity Not Infectious, …, Dead Improved Immunity.

Let Embedded Image be the probability of transition from disease compartment Embedded Image to Embedded Image. For the sake of simplicity, we dropped the second superindex (k) on zj,l, l ∈ {i, i + 1} here. Then, Embedded Image is the probability to recover from disease state Embedded Image. Let further Embedded Image be the time in days an individual stays in a compartment Embedded Image. For zj referring to a disease state with partial or improved immunity, we introduce the simplified factor κ as the counterpart of the relative reduction 1 − κ between Embedded Image and the corresponding time spent in the naive compartment, e.g., for nonsymptomatic individuals Embedded Image and

Additionally, we define reduction factors p* for the extended protection against severe courses in the partial and improved immunity layers. The basic functioning of these parameters can be described by Embedded Image where i ∈ {1, …, G} again denotes to the age group and zj,* ∈ {INS,*, ISy,*, ISev,*, ICr,*} refers to the disease states. Here, we added the index M ∈ {PI, II} describing the affiliation to the immunity layer. Summarized, we use reduction factors to reduce the probabilities, which we define for the naive immunity layer, to use them with two other layers. For more details about the reduction factors, we refer to [14]. Furthermore, we define Embedded Image and Embedded Image as the time spent in susceptible compartment Embedded Image and Embedded Image, respectively, if no infection event occurs, i.e., these times define the denominator of the waning rates. The (local) vaccination rates of the different immunity layers M ∈ {N, PI, II} are given by Embedded Image. Eventually, we define by Embedded Image and Embedded Image, the share of (locally) nonsymptomatic and symptomatic infectious individuals.

In this section, we now provide the elementary model processes in a general representation. For the full model equations, see S1 Appendix. SECIRS-type model equations and initialization. As the susceptible compartments contain the transmission processes to the exposed compartment as well as the waning immunity from the temporary immunity compartments and to the susceptible compartment of the inferior protection layer, we provide these equations first. For the three different groups of susceptibles M ∈ {N, PI, II}, i ∈ {1, …, G}, these write Embedded Image

Here, we used the simplified notation of M + 1 = PI if M = N and M + 1 = II if M = PI and for the sake of a shorter presentation, we skipped the dependence on t.

Furthermore, except for the exposed and dead compartments, we can either recover from a particular disease state to a temporary immunity state or move on to an aggravated state of the disease. For any given disease state zj ∈ ƵI excluding the different exposed states, we can generically write Embedded Image where we assume that zj−1,i is the previous milder infection state.

In the full model equations, we also provide confirmed compartments. However, there is no flow from undetected to detected compartments in the model equations. Detection within a node is modeled implicitly via Embedded Image and Embedded Image. The detected compartments are only used and filled via testing on traveling and commuting. For more details, see [21]. However, this mechanism is not used in the underlying study.

An explanation of all parameters for defining the model in Eq (13) - (14) is provided in Table 2. In order to solve the local models, we use an adaptive Runge-Kutta Cash Karp 5(4) integration scheme [22] to ensure a small discretization error.

View this table:
Table 2. Parameters used to define our multi-layer waning immunity model of SECIRS-type; Eq (13) - (14).

Parameterization

We will showcase our novel model with a synthetic outbreak scenario and according to the developments of SARS-CoV-2 in autumn 2022 in Germany. As reported case data is divided into the six age groups: 0-4 years, 5-14 years, 15-34 years, 35-59 years, 60-79 years, and 80 years and older, we stratify our model accordingly.

Mobility and contact patterns

The dataset user for inter-county commuting in Germany is based on a macroscopic transport model and a synthetic population for Germany. The transport model called DEMO [25] (short for DEutschlandMOdell, i.e., German Transport Model) is aimed to forecast transport in Germany when changing external factors like political policies or technological novelties. The model consists of multiple origin/destination (OD) matrices for each transport mode containing trip counts, average travel distances and average travel times for each origin/destination relation. Origins and destinations are described as traffic analysis zones. In this case the study area (Germany) is divided in 6633 traffic analysis zones which translates to matrices of size 6633×6633. The synthetic population is generated from the micro census 2017 and is spatially allocated to the traffic analysis zones. Additionally each synthetic person is binned into a behavioral homogeneous person group considering attributes like age and gender. In order to generate trip chains, e.g., home → work → leisure → shopping → home, daily routines including schedules, average travel distances and times are extracted from the national household survey Mobility in Germany 2017 [26]. Counts of each routine are transformed into discrete probability distributions for each person group; for a more detailed description on the process, see [27]. Further disaggregation iterates through each person where a daily routine is picked out of the according distribution. Since routines start and end at home, the start and end zone for the first respectively the last activity is known. For all activities in between, possible destination zones are selected based on travel distance from the starting zone and the trip count from the macroscopic model. The output is a list of trips between zones for each person. Finally the data is aggregated again to an origin/destination relation level with trip counts, transport mode and person group distribution. The preliminary age-group related OD matrices have been published freely in [28]. Note that in this study age group related traffic has been recomputed by (10).

For the underlying study with an aggregated mean-value model, we only allow one travel per day and thus take workplace-related commuting from the described data set. We assume that all trips are possible in the time line of one day each. For the given data, we provide the ratio of internal and inbound commuters in Fig. 9 (left) for each county in Germany, sorted in ascending order of the internal commuter ratios. In Fig. 9 (right), we present the ratio of inbound commuters to the total local population of the considered county. The figure shows, on the one hand, the influence of inbound commuters on the number of total workers and that the share compared to the total population is relevant, if not even substantial.

Fig 9.
Fig 9. Analysis of workplace commuting mobility in Germany.

Ratio of internal commuters to inbound commuters (left) and incoming commuters against the respective local population (right). The counties are sorted in ascending order.

Local contact patterns

In the previous sections, we have already explained how we model contact patterns in mobility and nonmobility settings, given survey numbers of contact ϕtr and ϕnt. The contact rate is an elementary part of the calculation of the transmission or infection rate of the persons in the susceptible compartments. As in [17], we select contact locations Home, Work, School, and Other and choose our baseline contact patterns based on [29] and [30]. To explicitly model, the contacts in transportation, we additionally subdivide the group of Other. This is done based on the contact surveys in [31], which numbers have been provided with [32] in an accessible format. From this data we infer the relative share of contacts in transportation (with respect to other), to apply these factors element-wise to the contact matrix from [29]. Inter- and extrapolation of the aggregated contact data to the given age groups is done with recent demographic data for Germany; [17]. The final contact baseline is shown in Fig. 10. Local contact reduction through NPIs is modeled as in [17] – selecting a zero minimum contact pattern.

Fig 10.
Fig 10. Age-resolved contact matrices for Germany.

The images show the average number of contacts in each of our five contact location. Values that are zero are displayed in white.

Epidemiological parameters for SARS-CoV-2 Omicron variant BA.1/2

In the last part of materials and methods, we focus on the epidemiological parameters specific to SARS-CoV-2 omicron variant BA.1/2, our demonstration case. The introduced model is a highly detailed model. We have explained the required parameters in Table 2. Depending on the available data, we provide stratified parameter ranges. Parameters have been obtained from an extensive literature research and prior findings with our previous models in [14, 17, 21]. For a summary see Table 3. For most parameters, we provide a range from which we uniformly take parameter samples. By making a large number of Monte Carlo runs, we obtain good estimate of the uncertainty in the output.

View this table:
Table 3. Model parameters for the Omicron Variant BA.1/2.

If parameters are stratified by age, vertical lines separate the corresponding columns. Column Reference provides used references in addition to the findings in [14, 17, 21].

To determine the initial states for the individual compartments of our model, we use the German case numbers published by the Robert Koch-Institute [53]. The idea of transforming the confirmed case data to initial states for each compartment was given in [14]. We follow this approach while implementing the changed immunity layers in this approach. For more details, we refer to S1 Appendix. SECIRS-type model equations and initialization.

Results

In this section, we use the introduced SECIRS-type model with three different immunity layers and two paces for waning immunity inside the novel travel-time aware mobility model.

We will consider two different settings. First, we provide a retrospective analysis of SARS-CoV-2 transmissions in the late phase of the pandemic. Second, we provide a synthetic outbreak scenario with a infection hotspot in Cologne and show how infection numbers spread differently with the newly introduced method.

Late-phase pandemic

We first consider the late-phase SARS-CoV-2 pandemic in Germany. We begin our simulation on August 1st, 2022 and examine the official face mask obligation in public transport, which was maintained until February 2023 and allowed both FFP2 and surgical masks [54]. Beginning October 1, 2022, nationwide enforcement of the use of FFP2 masks on long-distance public transportation [54] was implemented. However, there were differences in local public transportation between the German federal states. In Lower Saxony [55] and Hamburg [56] the FFP2 mask obligations also existed in local public transport but were dropped at the start of October 2022.

The newly introduced mobility scheme offers us a flexible approach to implement actions within the mobility node, without necessitating changes to the underlying local model. Within our model, compliance with the mask mandate can be represented by altering the effective contact rate through a damping of the contacts. As shown in [57], the proper usage of FFP2 masks is crucial. For both individuals correctly wearing FFP2 masks, the risk of infection during a 20-minute encounter with an infected person was estimated to be as small as 0.1 %. Notably, [57] also demonstrates that without any form of mask-wearing, the probability of infection approaches 100 %. Therefore, we implement the effect of FFP2 masks by a reduction factor of 0.1 % to the effective contact.

Throughout much of the pandemic, the German national testing strategy was predominantly focused on testing individuals displaying symptoms [58]. Therefore, we assume that the number of confirmed cases corresponds directly to the daily growth of symptomatic individuals if the detection ratio stays constant. This is also in line with the overall decreasing accuracy in positive-testing asymptomatic cases with regard to the Omicron variant [59]. As official reporting was suspected to have large underdetection factors, we based our comparisons not only on the official reported case numbers but also on the self-reported cases in the DigiHero [60] study. We extrapolated the reported infections from 33730 participants in autumn 2022 to the total population of Germany. In Fig. 11 (left), we provide the results of two different ensemble runs with 500 individual runs each where either no face masks (blue) or FFP2 (green) masks in transportation have been simulated. We visualize shaded areas for the regions between the p25 and p75 percentiles of the 500 runs and solid lines for the median results.

Fig 11.
Fig 11. Daily number of symptomatic infections and transmissions in the mobility model.

The simulations provide results using correctly worn FFP2 masks (green) and no face masks (blue) in the mobility models. The number of daily symptomatic infections (left) is compared to data provided by the DigiHero study (red) where the number of participants is extrapolated to the population in Germany and the official reported numbers (yellow). The colored areas between the dotted lines represent the p25 and p75 percentiles of the particular ensemble runs. The mean values are indicated by the solid lines. In the right plot, we visualize the daily transmission in the mobility models during the course of the simulation on a logarithmic scale.

Until beginning of October, Fig. 11 (left) shows substantially larger numbers of infected for the scenario without face masks. The quick decline in numbers in this scenario is due to the very large number of infected before and the temporal immunity saturation obtained. We see that face masks in transportation alone cannot prevent peaking of case numbers. However, the overall peak is considerably smaller (around 20 %). Face masks in transportation can thus contribute in ensuring functionality of critical infrastructure and preventing hospital bed shortages. Compared to reported data, we still see a nonnegligible difference. However, we also assume important underdetection, as test-positive rates range between 30 and 50 %; see [61]. In all curves, we similarly see a first stagnation in numbers and decline of infections towards end of the considered time frame.

An important advantage of our approach lies in its ability to quantify infections occurring during participation in traffic, achieved through the utilization of the traffic models. This offers a fundamental advantage over conventional structures. The significantly different probabilities of infection due to the implementation of a mask obligation in the transport models are clearly reflected in the daily number of transmissions in the traffic models obtained by our simulation; see Fig. 11(right).

Finally, we want to illustrate the spread of the disease on a geographical basis. In Fig. 12, we provide the (median) daily new numbers of symptomatic infected individuals on county level. This clearly emphasizes once again that the masks, as used only in the mobility models, have a significant influence on the spread of infection. The far higher peak in the case where no face masks are worn is clearly recognizable around day 50. Again, due to a temporal immunity seturation, we find a decline from day 70 to day 90 in the nonrestricted scenario while the face mask scenario numbers decline slower.

Fig 12.
Fig 12. Daily new number of symptomatic infected individudals on county level.

This figure illustrates the change in the daily number of symptomatic infected German counties over the course of 90 days beginning with the 1st August 2022. Each subfigure displays a different day with the two different modes (No masks or FFP2 masks). Using geodata “Verwaltungsgebiete 1:2 500 000, Stand 01.01. (VG2500)” from https://gdz.bkg.bund.de, copyright; GeoBasis-DE / BKG 2021, license dl-de/by-2–0, see https://www.govdata.de/dl-de/by-2-0.

Outbreak scenario

In this section, we want to show the impact of the mobility on the evolution of the disease using a fictional outbreak scenario. Therefore, we assume that in the initial stage to only have a single county which has a high amount of exposed people. In this context, we assume 5 % of Cologne’s population to be exposed to the virus. Although is number is ridiculously large, it was chosen by intention to directly visualize the qualitative propagation of new cases to connected counties.

To emphasize the influence of mobility, we compared our results with the mobility scheme presented in [17] (here denoted Existing method, both times using our multi-layer waning immunity model. The advantage of the novel mobility approach accounting for travel and transmission during commuting, contrasting the immediate transfer of travelers in [17].

In the maps of Fig. 13, we illustrate the progression of the total number of currently infected individuals using both approaches over 90 days (excluding Cologne to restrict the colorbar). In particular in the left-most part of the figure, we see that the virus attains counties which have not or almost not been affected with the existing method. We see that the novel approach is leading to a different propagation pattern also easily attaining indirect or loosely coupled regions. The very large number of exposed in Cologne leads to the rapid increase in the upper right part. The second and higher peak in October is due to the waning against any transmission and our seasonality factor in the transmission processes. If we look at the number of infections in all counties except Cologne, a comparison of the median values of the two methods shows that far more individuals will get infected using the novel approach. This is in accordance with the exponential increase in case number in so far virus-free regions. Similar process had been observed for the beginning of the pandemic when only a small number of regions was affected at all. The corresponding process is also reflected in the plots of the maps, where we already see a much more divided infection pattern on the first day.

Fig 13.
Fig 13. Propagation of infections with exiting [17] and novel model for an outbreak scenario.

Median propagation of infections in German counties (except Cologne) over 90 days with the existing method (top row maps) and the novel method (bottom row maps). In the plots of the right hand side, we provide the evolution of case numbers in Cologne and the rest of Germany. The shaded area visualizes the p25 and p75 percentile results of the 500 ensemble runs. The median value is given as solid line. Using geodata “Verwaltungsgebiete 1:2 500 000, Stand 01.01. (VG2500)” from https://gdz.bkg.bund.de, copyright; GeoBasis-DE / BKG 2021, license dl-de/by-2–0, see https://www.govdata.de/dl-de/by-2-0.

Discussion

The results presented in this study are based on a mean-value model, a multitude of assumptions and different parameters. Although based on an extensive literature review and previous findings in [14, 17, 21], there is still uncertainty about many (virus-specific) parameters. However, mathematical models can help to explore the mechanisms in infectious disease spread and reasonable assumptions are necessary, even for models up to digital twins. We have waived a lot of unrealistic assumptions such as homogeneous mixture in space and in demography. Our model has been resolved for meaningful regions and demography has been taken into account to model different vulnerabilities for different age groups.

A critical factor in infectious disease modeling is mobility. As the ground truth for daily mobility is difficult to grasp, again, assumptions have to be used. We have used the output of a state-of-the-art macroscopic mobility model to be used in our simulations. Different mobility patterns clearly lead to different virus distribution patterns.

A precise quantification of waning immunity against different courses of the disease given particular immunization histories is still an open research question. We have used measured patient values from recent literature and our models could correctly predict stagnation and decline after an initial phase of increased case numbers. The effect of waning immunity to SARS-CoV-2 over years, the interplay with new variants or vaccinations is still an open research question our model cannot answer yet.

In order to get meaningful predictions for future developements, our model needs to initialized with reasonable values of the population’s immunity. As the reporting of vaccinations is not conducted on personal level and as the dark figure of unreported transmission numbers increased with the end of the pandemic, the current status of the population is difficult to obtain. Serological studies may be needed to provide approximative starting values for future predictions. These limitations in data and assumptions should be kept in mind when interpreting the results of this work.

Conclusion

In this paper, we have presented a novel travel-time aware metapopulation model combined with a novel multi-layer waning immunity model based on ordinary differential equations. Both models advance already existing models for infectious disease dynamics by taking into account yet neglected but essential properties for virus propagation.

The travel-time aware metapopulation model advances existing metapopulation or hybrid graph-ODE approaches by considering travel time as well as transmissions during transport. As unconstrainted mobility clearly is a driver for infectious disease dynamics, more realistic mobility modeling increases the reliability of model outcomes. The travel-time aware metapopulation model can directly be combined with any kind of ODE-based transmission model and thus enhances global predictions by correctly resolving heterogeneously different spreading dynamics.

The novel multi-layer waning immunity model copes with the problem of heterogeneous and hybrid immunity in a population in mid- or late-phase pandemics as well as in endemic scenarios. The novel model does not only provide three different populations with different protection factors but also two different paces of waning immunity, i.e., against any transmission and against a severe course of the disease – a property that has been observed for SARS-CoV-2. Given reparametrization, it can also be used for other human-to-human transmittable diseases that show similar developments.

With the combination of both models, we could correctly assess developments in the late phase of the pandemic in Germany in 2022. Although vaccination could already mitigate the number of severe courses of the diseases, large numbers of mild cases can also have a desastrous impact on economy and critical infrastructure. Our model serves as a good basis for future waves SARS-CoV-2 in upcoming autumn or winter seasons and adaptation to other viruses is possible. The suggested mobility model can be of great help when particular interventions in public transportations are envisaged and expected transmissions can be assessed beforehand. As it can be combined with any kind of ODE-based model, it can also be used for other diseases such as Influenza or RSV.

Supporting information

S1 Appendix. SECIRS-type model equations and initialization.

Acknowledgments

The authors HZ, RS, AS, and MJK have received funding by the German Federal Ministry for Digital and Transport under grant agreement FKZ19F2211A, RS has received funding by the German Federal Ministry for Digital and Transport under grant agreement FKZ19F2211B. HZ, DK, MMH, and MJK have received funding from the Initiative and Networking Fund of the Helmholtz Association (grant agreement number KA1-Co-08).

Competing interests

The authors declare to not have any competing interests.

Data availability

The code is publicly available on Github under https://github.com/SciCompMod/memilio. The total number of produced simulation data sums up to 200 GB and cannot be uploaded easily. Particular data can be shared upon request.

Author Contributions

Conceptualization: Henrik Zunker, Martin Kühn, David Kerkmann

Data Curation: Henrik Zunker, Martin Kühn, Alain Schengen, Sophie Diexer

Formal Analysis: Henrik Zunker, Martin Kühn, David Kerkmann

Funding Acquisition: Martin Kühn, Michael Meyer-Hermann, Rafael Mikolajczyk

Investigation: Henrik Zunker, Martin Kühn, David Kerkmann, Alain Schengen

Methodology: Henrik Zunker, Martin Kühn, David Kerkmann

Project Administration: Martin Kühn

Resources: Martin Kühn, Michael Meyer-Hermann, Rafael Mikolajczyk

Software: Henrik Zunker, René Schmieding

Supervision: Martin Kühn, Michael Meyer-Hermann

Validation: All authors

Visualization: Henrik Zunker, Martin Kühn

Writing – Original Draft: Henrik Zunker, Martin Kühn

Writing – Review & Editing: All authors

References

  1. 1.↵
    World Health Organization Team Data, Analytics & Delivery. World health statistics 2023: monitoring health for the SDGs, sustainable development goals. Geneva: World Health Organization; 2023. Available from: https://www.who.int/publications/i/item/9789240074323.
  2. 2.↵
    Nordström P, Ballin M, Nordström A. Risk of infection, hospitalisation, and death up to 9 months after a second dose of COVID-19 vaccine: a retrospective, total population cohort study in Sweden. The Lancet. 2022;399(10327):814–823. doi:10.1016/S0140-6736(22)00089-7.
    OpenUrlCrossRefPubMed
  3. 3.↵
    Marcelin JR, Pettifor A, Janes H, Brown ER, Kublin JG, Stephenson KE. COVID-19 Vaccines and SARS-CoV-2 Transmission in the Era of New Variants: A Review and Perspective. Open Forum Infectious Diseases. 2022;9(5):ofac124. doi:10.1093/ofid/ofac124.
    OpenUrlCrossRef
  4. 4.↵
    Tagesschau. Diese Corona-Regeln gelten ab Oktober; 2022. Available from: https://www.tagesschau.de/inland/innenpolitik/corona-infektionsschutzgesetz-115.html.
  5. 5.↵
    Berndt JO, Conrad T, Hasenauer J, Karch A, Kheifetz Y, Kirsten H, et al. Szenarien für den Verlauf der SARS-CoV-2-Pandemie im Winter 2022/23 - Ergebnisse eines Workshops des Modellierungsnetzes für schwere Infektionskrankheiten (Modellierungsnetz). Zenodo; 2022. Available from: https://zenodo.org/record/7126032.
  6. 6.↵
    Müller SA, Balmer M, Charlton W, Ewert R, Neumann A, Rakow C, et al. Predicting the effects of COVID-19 related interventions in urban settings by combining activity-based modelling, agent-based simulation, and mobile phone data. PLOS ONE. 2021;16(10):e0259037. doi:10.1371/journal.pone.0259037.
    OpenUrlCrossRef
  7. 7.↵
    Tapp L, Kurchyna V, Nogatz F, Berndt JO, Timm IJ. School’s out? Simulating schooling strategies during COVID-19. In: International Workshop on Multi-Agent Systems and Agent-Based Simulation. Springer; 2022. p. 95–106.
  8. 8.↵
    Bauer S, Contreras S, Dehning J, Linden M, Iftekhar E, Mohr SB, et al. Relaxing restrictions at the pace of vaccination increases freedom and guards against further COVID-19 waves. PLOS Computational Biology. 2021;17(9):e1009288. doi:10.1371/journal.pcbi.1009288.
    OpenUrlCrossRef
  9. 9.
    Schüler L, Calabrese JM, Attinger S. Data driven high resolution modeling and spatial analyses of the COVID-19 pandemic in Germany. PLOS ONE. 2021;16(8):e0254660. doi:10.1371/journal.pone.0254660.
    OpenUrlCrossRef
  10. 10.↵
    Rodiah I, Vanella P, Kuhlmann A, Jaeger VK, Harries M, Krause G, et al. Age-specific contribution of contacts to transmission of SARS-CoV-2 in Germany. European Journal of Epidemiology. 2023;38(1):39–58. doi:10.1007/s10654-022-00938-6.
    OpenUrlCrossRef
  11. 11.↵
    Mohring J, Burger M, Feßler R, Fiedler J, Leithäuser N, Schneider J, et al. Starker Effekt von Schnelltests (Strong effect of rapid tests); 2023. Available from: http://arxiv.org/abs/2304.05938.
  12. 12.↵
    Kheifetz Y, Kirsten H, Scholz M. On the Parametrization of Epidemiologic Models—Lessons from Modelling COVID-19 Epidemic. Viruses. 2022;14(7):1468. doi:10.3390/v14071468.
    OpenUrlCrossRef
  13. 13.↵
    MacDonald N. Time lags in biological models. vol. 27 of Lecture Notes in Biomathematics. Berlin Heidelberg New York: Springer Science & Business Media; 1978.
  14. 14.↵
    Koslow W, Kühn MJ, Binder S, Klitz M, Abele D, Basermann A, et al. Appropriate relaxation of non-pharmaceutical interventions minimizes the risk of a resurgence in SARS-CoV-2 infections in spite of the Delta variant. PLOS Computational Biology. 2022;18(5):e1010054. doi:10.1371/journal.pcbi.1010054.
    OpenUrlCrossRef
  15. 15.↵
    Brauer F, Castillo-Chavez C, Feng Z. Mathematical models in epidemiology. vol. 32. Springer; 2019.
  16. 16.↵
    Liu J, Ong GP, Pang VJ. Modelling effectiveness of COVID-19 pandemic control policies using an Area-based SEIR model with consideration of infection during interzonal travel. Transportation Research Part A: Policy and Practice. 2022;161:25–47. doi:10.1016/j.tra.2022.05.003.
    OpenUrlCrossRef
  17. 17.↵
    Kühn MJ, Abele D, Mitra T, Koslow W, Abedi M, Rack K, et al. Assessment of effective mitigation and prediction of the spread of SARS-CoV-2 in Germany using demographic information and spatial resolution. Mathematical Biosciences. 2021; p. 108648. doi:10.1016/j.mbs.2021.108648.
    OpenUrlCrossRef
  18. 18.↵
    Müller SA, Paltra S, Rehmann J, Nagel K, Conrad TOF. Explicit modeling of antibody levels for infectious disease simulations in the context of SARS-CoV-2. iScience. 2023;26(9):107554. doi:10.1016/j.isci.2023.107554.
    OpenUrlCrossRef
  19. 19.↵
    Winkler C, Mocanu T. Methodology and application of a German national passenger transport model for future transport scenarios. In: Proceedings of the 45th European Transport Conference; 2017.
  20. 20.↵
    Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al. Social Contacts and Mixing Patterns Relevant to the Spread of Infectious Diseases. PLOS Medicine. 2008;5(3):e74. doi:10.1371/journal.pmed.0050074.
    OpenUrlCrossRefPubMed
  21. 21.↵
    Kühn MJ, Abele D, Binder S, Rack K, Klitz M, Kleinert J, et al. Regional opening strategies with commuter testing and containment of new SARS-CoV-2 variants in Germany. BMC Infectious Diseases. 2022;22(1):333. doi:10.1186/s12879-022-07302-9.
    OpenUrlCrossRef
  22. 22.↵
    Cash JR, Karp AH. A variable order Runge-Kutta method for initial value problems with rapidly varying right-hand sides. ACM Transactions on Mathematical Software. 1990;16(3):201–222. doi:10.1145/79505.79507.
    OpenUrlCrossRef
  23. 23.↵
    Patel MM, York IA, Monto AS, Thompson MG, Fry AM. Immune-mediated attenuation of influenza illness after infection: opportunities and challenges. The Lancet Microbe. 2021;2(12):e715–e725. doi:10.1016/s2666-5247(21)00180-4.
    OpenUrlCrossRef
  24. 24.↵
    Kühn MJ, Abele D, Kerkmann D, Korf SA, Zunker H, Wendler AC, et al. MEmilio v1.0.0 - A high performance Modular EpideMIcs simuLatIOn software; 2023. Available from: https://elib.dlr.de/201660/.
  25. 25.↵
    Winkler C, Mocanu T. Impact of political measures on passenger and freight transport demand in Germany. Transportation Research Part D: Transport and Environment. 2020;87:102476. doi:10.1016/j.trd.2020.102476.
    OpenUrlCrossRef
  26. 26.↵
    Nobis C, Kuhnimhof T. Mobilität in Deutschland – MiD: Ergebnisbericht; 2018. Available from: https://www.mobilitaet-in-deutschland.de/archive/pdf/MiD2017_Ergebnisbericht.pdf.
  27. 27.↵
    Mocanu T, Joshi J. Generating a microscopic dataset of passenger mobility from a national macroscopic travel demand model and a household travel survey. In: European Transport Conference 2022; 2022.Available from: https://elib.dlr.de/188443/.
  28. 28.↵
    Kühn M, Schengen A, Mocanu T, Joshi J. Vorläufige bundesweite Verkehrsströme; 2023. Available from: https://mobilithek.info/offers/573360269906817024..
  29. 29.↵
    Prem K, Cook AR, Jit M. Projecting social contact matrices in 152 countries using contact surveys and demographic data. PLOS Computational Biology. 2017;13(9):e1005697. doi:10.1371/journal.pcbi.1005697.
    OpenUrlCrossRefPubMed
  30. 30.↵
    Fumanelli L, Ajelli M, Manfredi P, Vespignani A, Merler S. Inferring the Structure of Social Contacts from Demographic Data in the Analysis of Infectious Diseases Spread. PLoS Computational Biology. 2012;8(9):e1002673. doi:10.1371/journal.pcbi.1002673.
    OpenUrlCrossRefPubMed
  31. 31.↵
    Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al. Social Contacts and Mixing Patterns Relevant to the Spread of Infectious Diseases. PLoS Medicine. 2008;5(3):e74. doi:10.1371/journal.pmed.0050074.
    OpenUrlCrossRefPubMed
  32. 32.↵
    Tomori DV, Rübsamen N, Berger T, Scholz S, Walde J, Wittenberg I, et al. Individual social contact data and population mobility data as early markers of SARS-CoV-2 transmission dynamics during the first wave in Germany—an analysis based on the COVIMOD study. BMC Medicine. 2021;19(1). doi:10.1186/s12916-021-02139-6.
    OpenUrlCrossRef
  33. 33.
    Kühn MJ, Abele D, Binder S, Rack K, Klitz M, Kleinert J, et al. Regional opening strategies with commuter testing and containment of new SARS-CoV-2 variants in Germany. BMC Infectious Diseases. 2022;22(1). doi:10.1186/s12879-022-07302-9.
    OpenUrlCrossRef
  34. 34.
    Águila-Mejía JD, Wallmann R, Calvo-Montes J, Rodríguez-Lozano J, Valle-Madrazo T, Aginagalde-Llorente A. Secondary Attack Rate, Transmission and Incubation Periods, and Serial Interval of SARS-CoV-2 Omicron Variant, Spain. Emerging Infectious Diseases. 2022;28(6):1224–1228. doi:10.3201/eid2806.220158.
    OpenUrlCrossRef
  35. 35.
    Manica M, Bellis AD, Guzzetta G, Mancuso P, Vicentini M, Venturelli F, et al. Intrinsic generation time of the SARS-CoV-2 Omicron variant: An observational study of household transmission. The Lancet Regional Health - Europe. 2022;19:100446. doi:10.1016/j.lanepe.2022.100446.
    OpenUrlCrossRef
  36. 36.
    Menni C, Valdes AM, Polidori L, Antonelli M, Penamakuri S, Nogal A, et al. Symptom prevalence, duration, and risk of hospital admission in individuals infected with SARS-CoV-2 during periods of omicron and delta variant dominance: a prospective observational study from the ZOE COVID Study. The Lancet. 2022;399(10335):1618–1624. doi:10.1016/S0140-6736(22)00327-0.
    OpenUrlCrossRef
  37. 37.
    Tobin RJ, Wood JG, Jayasundara D, Sara G, Walker CR, Martin GE, et al. Real-time analysis of hospital length of stay in a mixed SARS-CoV-2 Omicron and Delta epidemic in New South Wales, Australia. BMC Infectious Diseases. 2023;23(1):28. doi:10.1186/s12879-022-07971-6.
    OpenUrlCrossRef
  38. 38.
    Stein C, Nassereldine H, Sorensen RJD, Amlag JO, Bisignano C, Byrne S, et al. Past SARS-CoV-2 infection protection against re-infection: a systematic review and meta-analysis. The Lancet. 2023;401(10379):833–842. doi:10.1016/S0140-6736(22)02465-5.
    OpenUrlCrossRefPubMed
  39. 39.
    Bobrovitz N, Ware H, Ma X, Li Z, Hosseini R, Cao C, et al. Protective effectiveness of previous SARS-CoV-2 infection and hybrid immunity against the omicron variant and severe disease: a systematic review and meta-regression. The Lancet Infectious Diseases. 2023;doi:10.1016/S1473-3099(22)00801-5.
    OpenUrlCrossRefPubMed
  40. 40.
    Wang B, Andraweera P, Elliott S, Mohammed H, Lassi Z, Twigger A, et al. Asymptomatic SARS-CoV-2 Infection by Age: A Global Systematic Review and Meta-analysis. Pediatric Infectious Disease Journal. 2022;42(3):232–239. doi:10.1097/inf.0000000000003791.
    OpenUrlCrossRef
  41. 41.
    Nyberg T, Ferguson NM, Nash SG, Webster HH, Flaxman S, Andrews N, et al. Comparative analysis of the risks of hospitalisation and death associated with SARS-CoV-2 omicron (B.1.1.529) and delta (B.1.617.2) variants in England: a cohort study. The Lancet. 2022;399(10332):1303–1312. doi:10.1016/s0140-6736(22)00462-7.
    OpenUrlCrossRef
  42. 42.
    Ferguson NM, Laydon D, Nedjati-Gilani G, Imai N, Ainslie K, Baguelin M, et al. Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand. Imperial College; 2020. Available from: https://www.gov.uk/government/publications/impact-of-non-pharmaceutical-interventions-npis-to-reduce-covid-19-mortality-and-healthcare-demand-16-march-2020.
  43. 43.
    Verity R, Okell LC, Dorigatti I, Winskill P, Whittaker C, Imai N, et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. The Lancet Infectious Diseases. 2020;20(6):669–677. doi:10.1016/S1473-3099(20)30243-7.
    OpenUrlCrossRefPubMed
  44. 44.
    Wang J, Olsson S, Wehmeyer C, Pérez A, Charron NE, de Fabritiis G, et al. Machine learning of coarse-grained molecular dynamics force fields. ACS Central Science. 2019;5(5):755–767. doi:10.1021/acscentsci.8b00913.
    OpenUrlCrossRefPubMed
  45. 45.
    Karagiannidis C, Mostert C, Hentschker C, Voshaar T, Malzahn J, Schillinger G, et al. Case characteristics, resource use, and outcomes of 10 021 patients with COVID-19 admitted to 920 German hospitals: an observational study. The Lancet Respiratory Medicine. 2020;8(9):853–862. doi:10.1016/S2213-2600(20)30316-7.
    OpenUrlCrossRef
  46. 46.
    Stålcrantz J, Kristoffersen AB, Bøås H, Veneti L, Seppälä E, Aasand N, et al. Milder disease trajectory among COVID-19 patients hospitalised with the SARS-CoV-2 Omicron variant compared with the Delta variant in Norway. Scandinavian Journal of Public Health. 2022;50(6):676–682. doi:10.1177/14034948221108548.
    OpenUrlCrossRef
  47. 47.
    Relan P, Motaze NV, Kothari K, Askie L, Le Polain de Waroux O, Van Kerkhove MD, et al. Severity and outcomes of Omicron variant of SARS-CoV-2 compared to Delta variant and severity of Omicron sublineages: a systematic review and metanalysis. BMJ Global Health. 2023;8(7):e012328. doi:10.1136/bmjgh-2023-012328.
    OpenUrlAbstract/FREE Full Text
  48. 48.
    Armstrong RA, Kane AD, Cook TM. Outcomes from intensive care in patients with COVID-19: a systematic review and meta-analysis of observational studies. Anaesthesia. 2020;75(10):1340–1349. doi:10.1111/anae.15201.
    OpenUrlCrossRefPubMed
  49. 49.
    Karagiannidis C, Mostert C, Hentschker C, Voshaar T, Malzahn J, Schillinger G, et al. Case characteristics, resource use, and outcomes of 10 021 patients with COVID-19 admitted to 920 German hospitals: an observational study. The Lancet Respiratory Medicine. 2020;8(9):853–862.doi:10.1016/s2213-2600(20)30316-7.
    OpenUrlCrossRef
  50. 50.
    Ronchini C, Gandini S, Pasqualato S, Mazzarella L, Facciotti F, Mapelli M, et al. Lower probability and shorter duration of infections after COVID-19 vaccine correlate with anti-SARS-CoV-2 circulating IgGs. PLOS ONE. 2022;17(1):e0263014. doi:10.1371/journal.pone.0263014.
    OpenUrlCrossRef
  51. 51.
    Andrews N, Stowe J, Kirsebom F, Toffa S, Rickeard T, Gallagher E, et al. Covid-19 Vaccine Effectiveness against the Omicron (B.1.1.529) Variant. New England Journal of Medicine. 2022;386(16):1532–1546. doi:10.1056/nejmoa2119451.
    OpenUrlCrossRefPubMed
  52. 52.
    UK Health Security Agency. COVID-19 Vaccine Surveillance Report; 2022. Online. Available from: https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/1050721/Vaccine-surveillance-report-week-4.pdf.
  53. 53.↵
    Robert Koch-Institut. SARS-CoV-2 Infektionen in Deutschland; 2023. Zenodo. Available from: doi:10.5281/zenodo.7814948.
    OpenUrlCrossRef
  54. 54.↵
    Bundesregierung. Keine Maskenpflicht mehr im Fernverkehr; 2023. Available from: https://www.bundesregierung.de/breg-de/themen/coronavirus/maskenpflicht-entfaellt-2157682.
  55. 55.↵
    Verkehrsverbund Rhein-Berg. Im Ö PNV gilt weiterhin die Maskenpflicht; 2022. Available from: https://www.vrb-online.de/de/ueber-uns/news/im-%C3%B6pnv-gilt-weiterhin-die-maskenpflicht-296.
  56. 56.↵
    NDR. Maskenpflicht in Bussen und Bahnen: Kleine Ä nderung in Hamburg; 2022. Available from: https://www.ndr.de/nachrichten/hamburg/coronavirus/Hamburg-Kleine-Aenderung-der-Maskenpflicht-in-Bussen-und-Bahnen,corona11074.html.
  57. 57.↵
    Bagheri G, Thiede B, Hejazi B, Schlenczek O, Bodenschatz E. An upper bound on one-to-one exposure to infectious human respiratory particles. Proceedings of the National Academy of Sciences. 2021;118(49). doi:10.1073/pnas.2110117118.
    OpenUrlAbstract/FREE Full Text
  58. 58.↵
    Bundesgesundheitsministerium. Fragen und Antworten zu COVID-19 Tests;. https://www.bundesgesundheitsministerium.de/coronavirus/nationale-teststrategie/faq-covid-19-tests#c25573.
  59. 59.↵
    Wagenhaüser I, Knies K, Hofmann D, Rauschenberger V, Eisenmann M, Reusch J, et al. Virus variant–specific clinical performance of SARS coronavirus two rapid antigen tests in point-of-care use, from November 2020 to January 2022. Clinical Microbiology and Infection. 2023;29(2):225–232. doi:10.1016/j.cmi.2022.08.006.
    OpenUrlCrossRef
  60. 60.↵
    Klee B, Diexer S, Sarajan MH, Glaser N, Binder M, Frese T, et al. Regional Differences in Uptake of Vaccination against COVID-19 and Influenza in Germany: Results from the DigiHero Cohort. Vaccines. 2023;11(11). doi:10.3390/vaccines11111640.
    OpenUrlCrossRef
  61. 61.↵
    Robert Koch-Institut. SARS-CoV-2-PCR-Testungen in Deutschland; 2023. Available from: doi:10.5281/zenodo.8069452.
    OpenUrlCrossRef
Back to top
Posted March 02, 2024.
Download PDF

Supplementary Material

Email
Novel travel time aware metapopulation models: A combination with multi-layer waning immunity to assess late-phase epidemic and endemic scenarios
Henrik Zunker, René Schmieding, David Kerkmann, Alain Schengen, Sophie Diexer, Rafael Mikolajczyk, Michael Meyer-Hermann, Martin J. Kühn
medRxiv 2024.03.01.24303602; doi: https://doi.org/10.1101/2024.03.01.24303602
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area