- Research article
- Open Access
- Published:

# Development of the reproduction number from coronavirus SARS-CoV-2 case data in Germany and implications for political measures

*BMC Medicine*
**volume 19**, Article number: 32 (2021)

## Abstract

### Background

SARS-CoV-2 has induced a worldwide pandemic and subsequent non-pharmaceutical interventions (NPIs) to control the spread of the virus. As in many countries, the SARS-CoV-2 pandemic in Germany has led to a consecutive roll-out of different NPIs. As these NPIs have (largely unknown) adverse effects, targeting them precisely and monitoring their effectiveness are essential. We developed a compartmental infection dynamics model with specific features of SARS-CoV-2 that allows daily estimation of a time-varying reproduction number and published this information openly since the beginning of April 2020. Here, we present the transmission dynamics in Germany over time to understand the effect of NPIs and allow adaptive forecasts of the epidemic progression.

### Methods

We used a data-driven estimation of the evolution of the reproduction number for viral spreading in Germany as well as in all its federal states using our model. Using parameter estimates from literature and, alternatively, with parameters derived from a fit to the initial phase of COVID-19 spread in different regions of Italy, the model was optimized to fit data from the Robert Koch Institute.

### Results

The time-varying reproduction number (*R*_{t}) in Germany decreased to <1 in early April 2020, 2–3 weeks after the implementation of NPIs. Partial release of NPIs both nationally and on federal state level correlated with moderate increases in *R*_{t} until August 2020. Implications of state-specific *R*_{t} on other states and on national level are characterized. Retrospective evaluation of the model shows excellent agreement with the data and usage of inpatient facilities well within the healthcare limit. While short-term predictions may work for a few weeks, long-term projections are complicated by unpredictable structural changes.

### Conclusions

The estimated fraction of immunized population by August 2020 warns of a renewed outbreak upon release of measures. A low detection rate prolongs the delay reaching a low case incidence number upon release, showing the importance of an effective testing-quarantine strategy. We show that real-time monitoring of transmission dynamics is important to evaluate the extent of the outbreak, short-term projections for the burden on the healthcare system, and their response to policy changes.

## Background

The outbreak of the novel coronavirus SARS-CoV-2 (CoV/COVID-19) in China has induced a worldwide pandemic. The comparably high lethality in the elderly population and in patients with comorbidities [1, 2], together with a widely absent immunization of the population against the novel virus as well as the limited health system capacity estimated to become overwhelmed by an unlimited virus spreading [3], led to non-pharmaceutical interventions (NPIs) to reduce virus transmission mostly by reducing inter-individual contacts. The aim of these measures was to achieve at least a delay of viral spreading, allowing the healthcare system to extend its capacities and to treat less patients per time or, ideally, achieve a complete stop of viral spreading. The NPIs installed in Germany have been effective in containing viral dissemination [4]. Hence, in the light of economic damage incurred by restrictions [5], a gradual release of NPIs was decided with moderate effects on virus transmission. However, the length of the serial interval, which is in the range of 4–7.5 days (mean values) for CoV [6–8], and an inevitable delay in testing and reporting imply that any sudden outbreaks may be recognized too late and careful continuous monitoring of the infection dynamics on a regional basis is required. Thus, current political decisions need foundational information about current infection dynamics and their response to changes in NPIs such as partial release of contact restrictions or school openings, ideally on a regional basis. In fact, a declining or stable number of daily reported cases despite releasing measures can be misleading if the trend of the achieved reproduction number, the delay between changes in the infection dynamics, and their manifestations in reported case numbers are not taken into account. Furthermore, the high variance of the locally reported new cases adds to this uncertainty. Thus, it is extremely important to construct a model that not only captures the disease dynamics but also has the potential to provide information on the trend of the outbreak by considering the time-dependence of the reproduction number for COVID-19. The situation in Europe was recently analyzed [7]. Here, a systematic analysis of the development of the reproduction number over the time period of the COVID-19 outbreak in Germany and in all federal states of Germany is provided.

A second level of information necessary for political decisions on NPIs is the prospective development of the outbreak under different scenarios. A too early release of NPIs risks to abandon the current level of containment and to initiate a new wave of viral spreading [9]. A too long application of NPIs carries the risk of collateral damage and imposes a strong economic burden [5]. In view of the achieved reproduction number in Germany and its federal states by April 2020, a partial release of NPIs was decided, including partial school re-openings and resumption of catering and hotel business under certain restrictions, including mask obligations in public. The effects of such re-openings are hard to predict and require careful monitoring of local factors governing the infection dynamics and their implications for forecasting the immediate future development of the pandemic.

Whereas existing data-based simple algorithms (e.g., [10]) stand useful for estimating the time-varying reproduction number of an ongoing outbreak for many infectious diseases (e.g., measles, H1N1 swine flu, polio) using symptom onset data and less parameters (e.g., serial interval only), they, in general, are designed considering detection of all cases and are not suitable when the proportion of detected cases may change over time. Furthermore, the absence of reliable symptom onset data and heterogeneity of an effective infectious period among asymptomatic, pre-symptomatic, and symptomatic individuals demand for an alternative method to estimate the time-varying reproduction number by adjusting for such specific features of COVID-19. We have developed an ordinary differential equation (ODE)-based compartmental model specific for COVID-19 transmission dynamics and disease progression, and used it for quantitative evaluation of the time-varying reproduction number under the influence of NPIs in Germany and its federal states. In addition, our contribution retrospectively infers the usage of the healthcare system in Germany and offers short-term predictions based on current developments in terms of infections, number of non-critical hospital beds, and critical/intensive care units (ICUs) needed to treat patients with severe disease progression, as well as fatalities. This analysis provides additional information on when and how strongly to react to potential infection waves in order to avoid unacceptably high mortality and morbidity as well as excessive demands on the healthcare system. As a state-specific estimation of the reproduction numbers and a prospective estimation of the outbreak need to be up-to-date for the purpose of closely monitoring effects of policy changes, we provide daily updates of our analysis results online [11].

## Methods

The implemented SECIR (Susceptible-Exposed-Carrier-Infected-Recovered) model is a deterministic ODE model adapted to the specific properties of SARS-CoV-2 viral infections. It distinguishes healthy individuals without immune memory of COVID-19 (*S*), infected individuals without symptoms but not yet infectious (*E*), infected individuals without symptoms who are infectious (pre-symptomatic (*C*_{I}) and asymptomatic (*C*_{R}) carriers), infected symptomatic individuals who are not yet detected (*I*), and detected (*I*_{H,R}) and undetected (*I*_{X}) symptomatic patients. Further, compartments for hospitalization in non-critical (*H*_{U,R,S}) and critical/intensive care units (*U*_{D,R}) were introduced to monitor the load on the healthcare system. Detected patients recover from different states of the disease (*R*_{Z}) or die (*D*). Undetected individuals who went through the infection and recovered are also taken into account (*R*_{X}). The quantities are defined, and the model is summarized in Fig. 1 with parameters in Table 1. The model equations read

The rates \(R_{2,\dots,12}\) denote the inverse time of transition between the respective states and can be estimated from literature. Parameter *R*_{1} is fitted to the course of reported case numbers in a sliding time window and therefore is a time-varying parameter. Greek letters *α*,*μ*,*ρ*,*𝜗*, and *δ* denote fractions of individuals with a particular fate while other Greek letters, viz., *γ*,*χ*,*ω*, and *β*, reflect the intensity of interaction of corresponding infectious compartments with the susceptible population. The overall case fatality ratio (CFR = *ρ**𝜗**δ*) has a time-varying component modeled with a logistic function

where *t* corresponds to the day of the year starting from January 1, 2020, and *H*=0.139,*L*=0.007,*k*=0.145, and *t*_{0}=87.3 (these values are for *μ*=0.2) are obtained from fitting the curve for cumulative deaths to obtain a time-dependent case fatality rate (CFR), which changed over the course of the epidemic in Germany. This is due to changing testing frequencies [37] and the shifting age structure of the infected over time [38], which we assume to predominantly reflect in a time-varying rate of hospitalization (*ρ*(t)).

The time-dependent *ρ*(*t*) is effectively incorporating the time-varying age distribution of the infected people in the course of the epidemic. An explicit representation of the age distribution [39] was not favored in view of many unknown parameters. Note that hospitalization occurs from a quarantined compartment in the model, thereby having relatively less influence on the *R*_{t} values. Hence, *ρ*(t) as resulted from CFR(t) in the country level was retained for estimating federal state-specific time-varying reproduction numbers. The demographic differences may be more important for the analysis of smaller districts; however, case numbers in smaller districts might not be sufficient for a proper discrimination of age groups.

### Parameterization

The model parameters are critical for the overall behavior of the model and for the quality of the predictions derived from it. For the sake of robustness of the results, we followed two different strategies on how to determine the model parameters. The development for time- varying reproduction numbers in Germany (Fig. 3a) was presented with the parameter sets derived from both strategies.

The first strategy was to derive the estimated values (see Table 1) of the model parameters for Germany based on the available literature (e.g., [26]). Some disease-specific quantities such as the incubation period and infectious period are considered to be independent of a specific country. Uncertainty in the values of the parameters was invoked either using a percentage variation (20% unless otherwise specified) around their estimated mean values or by sampling from an estimated range (e.g., *𝜗* and *δ*) (see Table 1 for details). The resulting ranges were subsequently used to determine the distribution of *R*_{t} values.

In the second approach, we kept model parameters open in a broader range (see Additional file 1: Table 1 [27, 30, 31]) and fitted them to the cases reported in different regions of Italy until March 18, 2020, in a single stretch using our model assuming that the dynamics in this initial phase of the outbreak are not affected by the overwhelmed healthcare system. We optimized the model parameters to result in minimized error over this whole period. As the lockdown was announced on March 9, 2020, in Italy, an additional duration of 9 days (i.e., the sum of the incubation period of 5.2 days and a period of 3.7 days until clinical visit) was considered. Available data for cumulative infected, hospitalized, ICU, and deaths were fitted for Italy and for the regions where the first registered case was on February 28, 2020, or prior [40]. The diversity of resulting parameter values for the different regions in Italy (Fig. 2a) was used to derive a second range of the parameters to determine the distribution of *R*_{t} values (Fig. 3).

### Basic (*R*
_{0}) and time-varying (*R*
_{t}) reproduction number

The basic reproduction number *R*_{0} is defined as the expected number of secondary cases produced by a single infection in a completely susceptible population [44]. It can be calculated from the parameters of a compartmental model [45–48] after fitting the model to data for a given time period during the epidemic. While *R*_{0} provides valuable information on the viral dissemination dynamics in the absence of immunity and awareness of the epidemic, the dynamics of the epidemic over time will be heavily influenced by development of immunity in the population [49], policy changes to minimize infection risk [50], and individual behavior in response to public awareness of a disease [51]. Hence, a practically more useful quantity during an outbreak is the time-dependent reproduction number *R*_{t} describing the expected number of secondary cases per index case at a given time of the epidemic. This quantity has to be derived from incidence data over time and reflects the multifactorial impact of NPIs, behavior changes, seasonal effects, etc. on the dynamics of viral spread.

In epidemic models with multiple compartments, *R*_{0} can be derived with the next generation matrix method [47]. The compartments with infected individuals are divided into two contributions with respect to their dynamics: new infections entering the compartments and transfer of infected into and out of the compartment to other compartments. The Jacobian matrices of these two quantities *F* and *V* describe the generation of new infections and the transfer across compartments, respectively [52, 53]. The elements *G*_{ij} of *G*=*F**V*^{−1} are related to the expected number of secondary infections in compartment *i* caused by a single infected individual of compartment *j*.

The reproduction number *R*_{0} for the present model is given by the dominant eigenvalue of *G*, i.e.:

where *N*_{0} is the total population and *S*_{0} is the susceptible population, both at the start of the outbreak.

For our analysis, we use reported cases (*I*_{H}+*I*_{R}) in the first 4 days and impose the adjusted number of corresponding detectable infected cases (*N*_{i}/[(1−*α*)*μ*], where *N*_{i} represents the number of reported cases at day *i*, *i*=1,2,3,4), as exposed individuals consecutively at an earlier time point given by the sum of one incubation period and the duration until clinical registration following symptom onset (1/*R*_{2}+1/*R*_{3}+1/*R*_{11}=8.9 days). This assumption takes into account that the symptom onset from the first exposed individuals has not happened before the minimum duration of serial interval (around 4 days), and thus, the first reported cases shall represent independent sources of the virus rather than being the result of transmissions. Given the initial conditions and using the parameter sets in Table 1, the transmission parameter *R*_{1}, which mostly contains information on the contact frequency and probability of transmission per contact and, thus, best reflects the individual behavior in the population with respect to social distancing and other measures to minimize the infection risk, is used to optimize the fit of the model dynamics to the observed case data.

In order to assess the impact of political measures and changes in the population response along with other dynamic but unrecognized variables (such as seasonality) onto the development of the time-varying reproduction number *R*_{t}, the cumulative number of registered cases is used. The cumulative registered case number is compared to the sum of infected individuals and all subsequent states in the model, i.e., with *I*_{H}+*I*_{R}+*H*_{U}+*H*_{R}+*H*_{S}+*U*_{R}+*U*_{D}+*D*+*R*_{Z}. A time window of a width of 10 days is defined starting at the day of the first reported case (i) considering the time-difference between viral exposure and clinical registration of a case (which is around 9 days as per the model parameters) and (ii) assuming a delay of 1 day between the announcement of a new measure and changed personal behavior, the combination of which we use to define our default window size. For a window size of *WS* (our default *W**S*=10 days), *R*_{t} shown on a particular date *D*_{t} depicts the reproduction number observed over a period from *D*_{(t−WS)} to *D*_{t}. This allows to determine *R*(*t*_{0}) in the first 10 days and to define the initial conditions for the first sliding time window (Fig. 2b). Then, in repeating cycles, the best *R*(*t*_{k}) (with *k*=1...*N*) for each time window at the starting time *t*_{k} of the *k*th time window is determined by

where *ρ*(*t*_{k}) denotes the average value of the time-varying hospitalization rate in the *k*th time window. A new set of initial conditions is defined a day later, including the reduced fraction of susceptible individuals *S*(*t*_{k})/*N*(*t*_{k}), with *S*(*t*_{k}) and *N*(*t*_{k}) the values of susceptible and total population at the starting time *t*_{k} of the *k*th time window. Note that fatal cases reduce the total population. *R*_{1}(*t*_{k}) is determined by fitting to the data in this time window. In cycles, the time window is shifted 1 day later. The series of *R*_{t} values for each of the sliding time windows at time *t*_{k} is reported at the final date of the time window.

For the prospective study, the state of the model at the last time of *R*_{t} evaluation is kept and used as initial condition for the model.

To assess the impact of different scenarios, a set of constant reproduction numbers (*R*) was imposed based on the history of the epidemic to mimic release, maintenance, or intensification of NPIs (see the “Results” section). The cumulative number of infected individuals and the number of occupied ICUs, hospital beds, and deaths are reported. More observables are found at [11].

The distribution of observables and *R*_{t} values is generated by reiteration of the analysis under varying model parameters randomly drawn from a uniform distribution within the range provided in Table 1. The box plots in the figures show median, 25 and 75 percentiles, and minimum and maximum values from these analyses.

## Results

Based on the classical models of infection epidemics [52], we developed a mathematical model particularly adapted to the specificities of the COVID-19 outbreak (SECIR model, Fig. 1). For the evaluation of effects of NPIs and behavioral changes on viral spreading, a time-varying reproduction number *R*_{t} has to be estimated [10]. We opted for a shifting time window of 10 days (in Fig. 2b, a scheme with a window of 7 days is shown) in each of which *R*_{t} is determined, and developed an automatized algorithm for the fast analysis of the current *R*_{t} (see the “Methods” section, Eq. 19). Importantly, each time window is not analyzed independently but includes the history of the epidemic by starting from the saved state of the simulation at the beginning of each time window.

The cumulative reported cases are reproduced by the model in each time window, giving rise to a time evolution of the reproduction number *R*_{t} (Fig. 3a, b). The large initial value at February 28, 2020, results from a sudden increase of independent first reported infections, possibly related to people coming back from holiday. This leads to an overshoot of the *R*_{t} value in a strength depending on the size of the time window used for analysis (not shown). The initial estimates for *R*_{t} are not reliable because importation is the major contributor to the detected cases instead of local transmission events. Whereas the choice of the serial interval determines the *R*_{t} estimation (median values shown with black line) obtained from the EpiEstim package [7, 54], the *R*_{t} distribution resulting from our model is governed by literature-informed model parameters (blue curve) and detection ratio of the infected cases. Using parameters derived from fitting data from Italy yields similar results (magenta curve) (see Fig. 3a).

The nationwide NPIs imposed in Germany included the recommendation for cancelation of large events on March 10, 2020, followed by recommendation of self-isolation issued on March 12, 2020 [55]. A series of NPIs were implemented subsequently in close spacing, viz., restriction on individual movement, nationwide closure of schools and leisure-related venues on March 16, 2020, and extensive contact restrictions on March 22, 2020 [56, 57] (see Fig. 3c). The Apple mobility trend [58] observed in Germany until March 22, 2020, revealed an altered higher drift since the 7th week of the year (2020) showing a peak on February 22, 2020, and a subsequent declining trend (Additional file 1: Figure S1 [58, 59]). This was followed by a rapid decline in mobility since March 8, 2020 (Additional file 1: Figure S1). As the new cases registered were exposed 9 days earlier, the decline in the reproduction number until March 19–March 21, 2020, observed in the model is unlikely due to NPIs but can be attributed to behavioral changes. Although mobility in Germany showed an upward trend following March 22, 2020, despite extensive contact restrictions (Additional file 1: Figure S1), the reproduction number went downwards until mid-April 2020 achieving a value near unity as of April 6, 2020, in between and after a period of fluctuations, attained a minimum on May 4, 2020. In addition to demonstrating protective awareness among individuals, this illustrates that the NPIs imposed appear to have had a strong effect on the dynamics of the COVID-19 epidemic (see Fig. 3c).

NPIs were released in Germany on April 20, 2020, for the first time. Shops were opened, and a few days later, wearing masks became compulsory. The *R*_{t} value reacted with a delay of 15–19 days (Fig. 3c). On the 19th day following the first release, it increased by roughly 0.38, continued to be in a range of 0.66–0.75, and then decreased again by 0.15 (numbers mentioned for *μ*=0.2), presumably in response to the imposed masks. The second release of measures was widely implemented on May 11, 2020, and involved a cautious re-opening of child care and schools as well as restaurants. However, all of those were opened with imposed social distancing. Again, 19 days later, the *R*_{t} values increased to 0.82 as of May 30, 2020. Following a short span of slight fall in the reproduction numbers, this remained in a range of 0.80–0.88 around June 12, 2020, for Germany (Fig. 3c). This observation illustrates the sensitivity of the viral spreading to NPIs as well as the possibility to partially release NPIs without losing control of the epidemic, provided the population keeps social distancing and hygiene measures in place, and avoids inter-personal contacts [60].

While there was a large diversity of epidemic onset and intermediate developments particular to individual federal states, the overall tendency converged to values below *R*_{t} = 1 around the first week of April 2020 (see Additional file 1: Figure S3 [41–43]). The coherent reduction of the reproduction number after nationwide implementation of several NPIs together with further measures specifically applied in different federal states speaks for the efficiency of the measures and the responsiveness of the population to the NPIs. Re-opening and a continuously increasing trend in mobility induced resurgence of COVID-19 outbreaks in multiple federal states, such as North Rhine-Westphalia, Berlin, Brandenburg, Saxony, and Saxony-Anhalt around the second half of May and the first half of June 2020 (Fig. 3d and Additional file 1: Figure S3). This resulted in a sharp rise in reproduction numbers ranging from 2 to 2.8 in these aforementioned states around the 3rd week of June. It is also reflected in the reproduction numbers of whole Germany around that time, reaching a peak on June 20, 2020. The *R*_{t} values in Germany remained more than unity during June 17, 2020, to June 27, 2020, before maintaining around 0.63–0.82 until July 13, 2020, due to reimposed regional regulations in several federal states. Since then, it had increased and remained mostly in the range of 1.22–1.56 prior to displaying a downward trend since the last week of August. Such variations can be understood as an overall impact resulting from alterations in *R*_{t} development in individual federal states (see Additional file 1: Figure S3). One interesting observation, however, is some level of correlation between a couple of federal states with regard to the evolution of their reproduction numbers, e.g., (i) Berlin and Brandenburg and (ii) Baden-Würtemberg and Bavaria (Fig. 3d). We note that Berlin is encircled by Brandenburg whereas Baden-Würtemberg and Bavaria are neighboring federal states located at similar geographical altitude with Alps in their southern part. Such a correlation suggests that resurgence of outbreak in one region may act like a reservoir of new infections in adjacent regions.

Next, the *R*_{t} distribution resulting from the median values obtained in the last week in the different federal states is compared and ranked in Fig. 3e. As of August 31, 2020, most of the federal states and Germany as a country exhibited a weekly median of *R*_{t} around 1 or higher. Bavaria and Baden-Würtemberg which were hit early on by the COVID-19 outbreak as well as Bremen resurfaced with consistently higher reproduction numbers during the last 3 weeks of August 2020 (see Additional file 1: Figure S3). In contrast, North Rhine-Westphalia which had a significant case load during early weeks in the pandemic and also exhibited local outbreaks and super-spreading events around May–June [61] was relatively doing better in controlling the outbreak by the end of August 2020 (Fig. 3e). Among other federal states which encountered a substantial sharp increase in *R*_{t} at some point after re-opening (see Additional file 1: Figure S3), Mecklenburg-Western Pomerania, Berlin, Saxony-Anhalt, Brandenburg, and Hamburg showed a median *R*_{t} less than or around 1 whereas Saxony displayed a median *R*_{t} higher than 1 for the week ending on August 31, 2020 (Fig. 3e and Additional file 1: Figure S3).

The number of unregistered cases is not well known in Germany. In the model, (1−*α*)*μ* captures the registered fraction of the infected cases (Fig. 1). In order to demonstrate the importance of the number of undetected cases for the interpretation of the results, we compared the results for *μ*=1,*μ*=0.5,*μ*=0.2, and a qualitative time-varying *μ*_{t} informed by the mobility and testing data (see Additional file 1: Figure S1-S2 [32, 41–43, 62, 63]) for a fixed proportion of purely asymptomatic as well as unregistered individuals, i.e., *α* (*α*=0.22 unless otherwise specified). For time-invariant detection ratios [ (1−*α*)*μ*], it turns out that the *R*_{t} value derived from a model with more symptomatic unregistered cases (i.e., a lower *μ*) is slightly enhanced but remains in the same range, the impact being more prominent during initial weeks (Fig. 3c). Temporal evolution of *R*_{t} for a realistic time-varying detection ratio captures the sensitivity towards NPIs better and clearly shows the timeline of induced changes in *R*_{t} due to phased re-opening with a similar delay (Fig. 3c, green curve). A consistent rise in *μ*_{t} prior to June 2020 resulted in lower *R*_{t} values from April 2020 onward. Following a peak detection of infected cases around the end of May 2020, *μ*_{t} fell due to a lower number of tests per confirmed case compared to the increased mobility (see Additional file 1: Figure S1), causing higher *R*_{t} values since June 2020. The resulting reproduction numbers are not significantly sensitive to changes in *α* (see Additional file 1: Figure S4). In our model, the infectious period of individuals who remain asymptomatic throughout is assumed to be shorter than the overall infectious period (including pre-symptomatic stage) of the symptomatic people. This results in a slightly lower reproduction number as we increase *α*.

The model can be used to estimate the dynamics of the load for the healthcare system. Based on the resulting fitting of cumulative detected cases (Fig. 3b) and cumulative deaths (Fig. 4a, see the “Methods” section), we investigated the extent of hospitalization during the epidemic (Fig. 4 and Additional file 1: Figure S5). The number of deaths and new daily reported cases well captured the trend in the data. The highest number (median) of estimated daily reported cases was 5727 (Fig. 4b). The estimated peak (median) for healthcare usage showed 10,690 occupied hospital beds (all types of non-critical care units) on March 28, 2020, and 4938 ICU beds (all types of critical care units) on April 4, 2020 (Fig. 4c, d). These numbers stayed within the capacities of the German healthcare system [64]. The sensitivity of this result to changes in the model parameters is shown for our retrospective analysis in Additional file 1: Figure S6.

We next investigated different prospective scenarios from the final date of the data fitting phase by retaining the state information of the model. We used the hospitalization rate estimated as on August 31, 2020 (*ρ*(*t*_{latest})), and plugged this into Eq. 19 while imposing different reproduction numbers (*R*) for the whole period of prospective analysis. *ρ*(*t*_{latest}) primarily depends on the affected age groups and the extent of an effective testing-quarantine strategy (Fig. 5a). Starting from the last state of the model for Germany, thus, including the complex distribution of individuals onto the different compartments of the model at this time, the simulation was first continued for 28 days with the mean of median *R*_{t} values observed during the last week of August 2020, i.e., with *R* = 1.03 (base scenario). It provides a stable situation without a significant resurgence of cases in short time (Fig. 5b–f, black). The median *R*_{t} value observed over a period from May 5, 2020 (as the first re-opening started to show its first impact from this date), until August 31, 2020, which resulted in *R* = 0.91, contained the epidemic but was not able to stop it in a short time (Fig. 5b–f, magenta).

A release of NPIs to a degree that induces *R*_{t} to be the maximum value estimated during May 5, 2020, to August 31, 2020, i.e., *R* = 1.84, causes a resurgence of around 7000 daily new reported cases at the end of the 4th prospective week but does not cause a significant burden to healthcare in near future (Fig. 5b–f, green). Provided this reproduction number was kept for 1 year, around 19,500 critical care beds and around 40,000 non-critical care beds would be occupied at peak causing an overwhelmed healthcare system in comparison to around 9600 free ICUs on August 31, 2020, as per [64]. This scenario may lead to total of 80,000 deaths (see Additional file 1: Figure S7).

A relatively pessimistic prospective scenario reflecting a complete release of NPIs and the absence of any protective behavioral measure was modeled with the median *R*_{t} observed on March 21, 2020, i.e., with *R* = 4.33. This leads to a major immunization of the population and results in a drastic increase in healthcare usage (e.g., ≈ 7000 occupied critical care beds) and cumulative deaths (≈ 1900 higher than the base scenario) within 4 weeks (Fig. 5, red). Consequently, in the long run, it leads to an overwhelmed healthcare system (see Additional file 1: Figure S7).

An important question is how long NPIs would have to be kept in place until all new cases can be controlled by public health departments. Assuming 300 detected cases per day to be manageable, we calculated the time needed to achieve this number of daily new cases given different levels of the reproduction number (Fig. 6). Given the latest reproduction number of August 31, 2020, in Germany (Fig. 6a), this number could be achieved within 2 months for registered cases, and 5 months for the infectious cases if *R*_{t} remains at the same level throughout (Fig. 6). The detection ratio ((1−*α*)*μ*) influences this time as can be seen for different values of *μ* (*α* = 0.22 fixed). Although the duration to reach 300 new registered cases is comparable for *μ*=1 and *μ*=0.2, the infectious population significantly increases with larger fractions of undetected infections (Fig. 6, black versus red).

## Discussion

The estimated fraction of the immunized population that we calculate as the ratio of recovered to initial susceptible population assuming some form of long-lasting immunity of infected survivors stands at the range of 0.3–3% (0.3% with *μ*=1, 1.7% with *μ*=0.2, and 2.8% for our estimated time-varying *μ*_{t}; all with *α*=0.22 [32]) at the end of August 2020. Hence, the German population seems far from achieving herd immunity and a renewed outbreak or a second wave is extremely likely in the absence of NPIs or continued behavioral changes to prevent a viral infection. While continuation of the epidemic with *R*_{t} values close to unity would avoid a large number of COVID-19 infections and deaths, it may lead to a major economic burden [5], induce unemployment and related collateral damages, and increase risks of suicide in certain individuals due to lasting social distancing [65] and could pose a strong load on the healthcare system due to a delayed/neglected treatment of other diseases [66].

As expected, intensifying contact restriction measures accelerates reaching a low number of daily new cases whereas an increase in *R*_{t} values delays the reduction of the daily new cases. In view of the non-linearity of the delay with larger reproduction numbers, one may speculate that it might be more advantageous in terms of health and economy to maintain a lower *R*_{t} value less than unity through NPIs, social distancing, and protective social culture in order to quickly achieve a controllable and traceable low number of daily new cases [67]. A complete elimination of the virus, as it appears in the model, is hard to achieve in reality because of open borders to European neighbors and unknown viral properties which might allow it to reappear under particular conditions. A combined strategy of rapid identification and isolation of infection as well as efficiently traced contact clusters also worked well in Japan [68]. We showed that a large number of undetected cases increase the delay in reaching a target number of daily new cases. While seroprevalence studies suggest undetected cases substantially less than 10-fold of the detected cases in the studied communities [32, 69], the overall true number of undetected cases in Germany is still not known.

The predictive power of the model was analyzed by comparing a forecast of cumulative new reported cases for 1 or 2 weeks based on the *R*_{t} value in the past week with the corresponding cumulative increase observed in reality (Additional file 1: Figure S8). Regional outbreaks cause a sudden rise in *R*_{t} values globally, such that in this case, the prediction on the country level overestimates cumulative new cases in the upcoming days. One intriguing finding is the tendency of underestimating the cumulative new reported cases on a scale of a couple of weeks prior to a positive overshoot in prediction error (see Additional file 1: Figure S8). It can be inferred as upcoming structural increase in *R*_{t} estimations (see *R*_{t} values in June 2020, Fig. 3c). Apart from the periods of sudden rise in *R*_{t} values, the model prediction works well on a scale of a couple of weeks. Model predictions work assuming a time-invariant detection ratio in the upcoming days, which, in reality, is unlikely the case for all weeks. This would result in an overestimation when the detection ratio falls or an underestimation when the detection ratio rises. Due to having a lower chance of missing an infection induced death in the data and a longer delay to death following viral exposure, the prediction for the death toll is excellent up to 5 weeks (e.g., < 0.5% error while continuing the projection based on the fitting until August 31, 2020). It illustrates the usefulness of the model in determining the burden on the healthcare system in the near future at least on a scale of a month.

Long-term prediction for any pandemic, especially a new one, is challenging due to several less known or unpredictable factors which may impact its transmission dynamics and its effective potency. Examples of these factors include the impact of accumulation of aerosols with viable virus in closed rooms [70, 71], extent of aerosol-mediated transmission and seasonality [72, 73], alterations of behavioral response, future NPIs, and viral mutations. In addition, development of efficient rapid testing methods, extent of reinfection and inherent immunity, and improvement in clinical management would determine and modify the future course of the outbreak. Within these limitations, our model can still guide the government authorities to prepare better by projecting the peak of healthcare usage and estimates of population immunization by the pathogen as well as case fatalities under different circumstances (see Additional file 1: Figure S7).

The analysis of the individual federal states in Germany revealed local differences. The federal states appear to witness different phases of the outbreak, and NPIs exhibit different kinetics of impact. Therefore, state-specific or even district/city-specific development of *R*_{t} provides a better sense to the local authorities to plan the future course of actions to control the epidemic locally, which, if left uncontrolled, may act like a disease hot spot to initiate new clusters of infections across the federal states. One intriguing finding is that a local outbreak results in sharper changes in *R*_{t} while its global impact as observed in the country level *R*_{t} estimates is relatively damped (Fig. 3c, d, Additional file 1: Figure S3). The full analysis for *R*_{t} of all federal states is available at [11] and clearly emphasizes local heterogeneity of the epidemic. Prevalence of super-spreading events, population density, and differences in social structures can be some of the contributing factors in driving the outbreak heterogeneously across different states.

We do not use an age-stratified version of our model for the presented analysis due to incomplete age-stratified data, lack of knowledge on how the model parameters would depend on different age groups, and an undetermined and uneven testing bias across different age groups. The age-dependence is phenomenologically included in the model by using a time-dependent hospitalization rate, which reflects the demography of the infected people. With this, our goal to understand the development of time-varying reproduction number, overall usage of healthcare facilities, and future course of the outbreak can still be achieved using an age-independent mean-field approach. Even though observed case fatality ratios (and indeed infection fatality rates, too) for older age groups are much higher than for young adults and children in COVID-19 [62], we can still fit the cumulative death curve based on the estimated (informed by the literature) fraction of the hospitalized COVID-19 patients who are dying thereby enabling us to capture the time-dependent hospitalization rate eventually resulting in a time-dependent overall CFR for the reported cases.

## Conclusions

In this paper, we developed a compartmental model (SECIR) accounting for the specificities of the recent COVID-19 outbreak. We reported an adaptive methodology to estimate the time-varying reproduction number (*R*_{t}) based on the incidence of reported cases. As parameterization is essential for the quality of our analysis and predictions, two reference parameter sets were determined by thorough analysis of the literature on COVID-19 and an analysis of Italian data. The results discussed are consistent between both parameter sets. Even though both parameter sets are not completely independent, this consistency increases the credibility of the model results. Implementation of NPIs in close spacing, heterogeneity in their application, and withdrawal timings across different federal states as well as cities make it statistically uncertain to disentangle the impact of a particular NPI [7]. Most importantly, the behavior of the people changes over time, examples of which include dynamic and heterogeneous compliance to NPIs and mask usage as well as behavioral exhaustion. Even though the model is constructed by taking into account the biological characteristics of the infection transmission dynamics of COVID-19 such as asymptomatic and pre-symptomatic transmission, and contribution from undetected mild-symptomatic cases, the model presented here can be easily translated to any similar infectious disease with equivalent features while such a methodology can also be applied to common infectious diseases. In addition, our results can capture the qualitative aspects of how the infection incidence, patients admitted in non-critical and critical care units, and deaths change as days progress during the COVID-19 outbreak. Furthermore, it can also guide the authorities in assessing how the pandemic would evolve in the near future and what load on the healthcare facilities to expect under certain scenarios drawn from the history of the outbreak itself. We provide a daily updated evaluation of the reproduction number suitable to support political decisions on NPIs in the course of the COVID-19 outbreak and applied to German federal state data online [11].

## Availability of data and materials

Daily updated results are provided in a version controlled repository [11] including a full revision history. The analysis is based on case data for Germany provided by the Robert-Koch-Institut as on October 27, 2020 [41]. Note that the case numbers in the data might change for newer data sets due to delayed reporting and data corrections.

## Change history

### 23 April 2021

The original article has been corrected to include an OA funding note.

## References

- 1
Istituto Superiore di SanitÃ. Report sulle caratteristiche dei pazienti deceduti positivi a COVID-19 in Italia. Ilpresente report è basato sui dati aggiornati al 20 Marzo 2020. Istituto Superiore di SanitÃ. 2020. https://www.epicentro.iss.it/coronavirus/sars-cov-2-sorveglianza-dati. Accessed 23 Mar 2020.

- 2
Sanyaolu A, Okorie C, Marinkovic A, et al.Comorbidity and its Impact on Patients with COVID-19. SN Compr Clin Med. 2020; 2:1069–76. https://doi.org/10.1007/s42399-020-00363-4.

- 3
an der Heiden M, Buchholz U. Modellierung von Beispielszenarien der SARS-CoV-2-Epidemie 2020 in Deutschland. 2020. https://www.rki.de/DE/Content/InfAZ/N/Neuartiges_Coronavirus/Modellierung_Deutschland.html. Accessed 30 Mar 2020.

- 4
Hartl T, Wälde K, Weber E. Measuring the impact of the German public shutdown on the spread of Covid-19. Covid Econ. 2020; 1:25–32.

- 5
Dorn F, Fuest C, Marcell G, Krolage C, Lautenbacher S, Link S, Peichl A, Reif M, Sauer S, Stöckli M, Wohlrabe K, Wollmershäuser T. Die volkswirtschaftlichen Kosten des Corona-Shutdown für Deutschland: Eine Szenarienrechnung. ifo Schnelldienst. 2020; 73(04):29–35. Accessed 07 Apr 2020.

- 6
Nishiura H, Linton NM, Akhmetzhanov AR. Serial interval of novel coronavirus (COVID-19) infections. Int J Infect Dis. 2020; 93:284–6. https://doi.org/10.1016/j.ijid.2020.02.060. Accessed 07 Apr 2020.

- 7
Flaxman S, Mishra S, Gandy A, Unwin HJT, Coupland H, Mellan TA, Zhu H, Berah T, Eaton JW, Guzman PNP, Schmit N, Callizo L, Ainslie KEC, Baguelin M, Blake I, Boonyasiri A, Boyd O, Cattarino L, Ciavarella C, Laura C, CucunubÃ¡ Z, Cuomo-Dannenburg G, Dighe A, Djaafara B, Dorigatti I, van Elsland S, FitzJohn R, Fu H, Gaythorpe K, Geidelberg L, Grassly N, Green W, Hallett T, Hamlet A, Hinsley W, Jeffrey B, Jorgensen D, Knock E, Laydon D, Nedjati-Gilani G, Nouvellet P, Parag K, Siveroni I, Thompson H, Verity R, Volz E, Walker PGT, Wang H, Wang Y, Watson O, Whittaker C, Winskill P, Xi X, Ghani A, Donnelly CA, Riley S, Okell LC, Vollmer MAC, Ferguson NM, Bhatt S. Estimating the number of infections and the impact of non-pharmaceutical interventions on COVID-19 in 11 European countries. Imperial College London. 2020. https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/covid-19/report-13-europe-npi-impact/. Accessed 07 Apr 2020.

- 8
Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, Ren R, Leung KSM, Lau EHY, Wong JY, Xing X, Xiang N, Wu Y, Li C, Chen Q, Li D, Liu T, Zhao J, Liu M, Tu W, Chen C, Jin L, Yang R, Wang Q, Zhou S, Wang R, Liu H, Luo Y, Liu Y, Shao G, Li H, Tao Z, Yang Y, Deng Z, Liu B, Ma Z, Zhang Y, Shi G, Lam TTY, Wu JT, Gao GF, Cowling BJ, Yang B, Leung GM, Feng Z. Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia. N Engl J Med. 2020; 382(13):1199–207. https://doi.org/10.1056/NEJMoa2001316. Accessed 07 Apr 2020.

- 9
Flaxman S, Mishra S, Gandy A, Unwin H, Mellan T, Coupland H, Whittaker C, Zhu H, Berah T, Eaton J, Monod M, Team ICC-R, Ghani A, Donnelly C, Riley S, Vollmer M, Ferguson N, Okell L, Bhatt S. Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Nature. 2020; 584(7820):257–61. https://doi.org/10.1038/s41586-020-2405-7. Accessed 24 June 2020.

- 10
Cori A, Ferguson NM, Fraser C, Cauchemez S. A new framework and software to estimate time-varying reproduction numbers during epidemics. Am J Epidemiol. 2013; 178(9):1505–12. https://doi.org/10.1093/aje/kwt133. Accessed 07 Apr 2020.

- 11
Systems Immunology Group. Complete and up-to-date analysis of Germany and all federal states. 2020. http://secir.theoretical-biology.de. Accessed 08 June 2020.

- 12
Zhao S, Gao D, Zhuang Z, Chong M, Cai Y, Ran J, Cao P, Wang K, Lou Y, Wang W, Yang L, He D, Wang M. Estimating the serial interval of the novel coronavirus disease (COVID-19): a statistical analysis using the public data in Hongkong from January 16 to February 15, 2020. MedRxiv. 2020. https://doi.org/10.1101/2020.02.21.20026559. Accessed 08 June 2020.

- 13
Wang D, Hu B, Hu C, Zhu F, Liu X, Zhang J, Wang B, Xiang H, Cheng Z, Xiong Y, Zhao Y, Li Y, Wang X, Peng Z. Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan, China. JAMA. 2020; 323(11):1061–9. https://doi.org/10.1001/jama.2020.1585. Accessed 2020-02-09.

- 14
World Health Organization. Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19). 2020. https://www.who.int/publications/i/item/report-of-the-who-china-joint-mission-on-coronavirus-disease-2019-(covid-19). Accessed 24 June 2020.

- 15
Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, Liu L, Shan H, Lei CL, Hui DSC, Du B, Li LJ, Zeng G, Yuen KY, Chen RC, Tang CL, Wang T, Chen PY, Xiang J, Li SY, Wang JL, Liang ZJ, Peng YX, Wei L, Liu Y, Hu YH, Peng P, Wang JM, Liu JY, Chen Z, Li G, Zheng ZJ, Qiu SQ, Luo J, Ye CJ, Zhu SY, Zhong NS, China Medical Treatment Expert Group for Covid-19. Clinical characteristics of coronavirus disease 2019 in China. N Engl J Med. 2020; 382(18):1708–20. https://doi.org/10.1056/NEJMoa2002032. Accessed 07 Apr 2020.

- 16
Cai Q, Huang D, Ou P, Yu H, Zhu Z, Xia Z, Su Y, Ma Z, Zhang Y, Li Z, He Q, Fu Y, Chen J, Liu L. 2019-nCoV pneumonia in a normal work infectious diseases hospital besides Hubei Province, China. SSRN J. 2020. https://doi.org/10.2139/ssrn.3542163. Accessed 07 Apr 2020.

- 17
Wan X, Wang W, Liu J, Tong T. Estimating the sample mean and standard deviation from the sample size, median, range and/or interquartile range. BMC Med Res Methodol. 2014; 14(1):135.

- 18
Mizumoto K, Kagaya K, Zarebski A, Chowell G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the diamond princess cruise ship, Yokohama, Japan, 2020. Euro Surveill. 2020; 25(10):2000180. https://doi.org/10.2807/1560-7917.ES.2020.25.10.20001. Accessed 07 Apr 2020.

- 19
Nishiura H, Kobayashi T, Miyama T, Suzuki A, Jung SM, Hayashi K, Kinoshita R, Yang Y, Yuan B, Akhmetzhanov AR, Linton NM. Estimation of the asymptomatic ratio of novel coronavirus infections (COVID-19). Int J Infect Dis. 2020; 94:154–5. https://doi.org/10.1016/j.ijid.2020.03.020. Accessed 07 Apr 2020.

- 20
Bi Q, Wu Y, Mei S, Ye C, Zou X, Zhang Z, Liu X, Wei L, Truelove SA, Zhang T, Gao W, Cheng C, Tang X, Wu X, Wu Y, Sun B, Huang S, Sun Y, Zhang J, Ma T, Lessler J, Feng T. Epidemiology and transmission of COVID-19 in Shenzhen, China: analysis of 391 cases and 1,286 of their close contacts. medRxiv. 2020. https://doi.org/10.1101/2020.03.03.20028423. Accessed 07 Apr 2020.

- 21
He X, Lau EH, Wu P, Deng X, Wang J, Hao X, Lau YC, Wong JY, Guan Y, Tan X, et al. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nat Med. 2020; 26(5):672–5.

- 22
Hao X, Cheng S, Wu D, Wu T, Lin X, Wang C. Reconstruction of the full transmission dynamics of COVID-19 in Wuhan. Nature. 2020; 584(7821):420–4.

- 23
Woelfel R, Corman V, Guggemos W, Seilmaier M, Zange S, Mueller M, Niemeyer D, Vollmar P, Rothe C, Hoelscher M, Bleicker T, Bruenink S, Schneider J, Ehmann R, Zwirglmaier K, Drosten C, Wendtner C. Clinical presentation and virological assessment of hospitalized cases of coronavirus disease 2019 in a travel-associated transmission cluster. medRxiv. 2020. https://doi.org/10.1101/2020.03.05.20030502. Accessed 12 Mar 2020.

- 24
Clifford S, Quilty BJ, Russell TW, Liu Y, Chan Y-WD, Pearson CA, Eggo RM, Endo A, Flasche S, Edmunds WJ, et al. Strategies to reduce the risk of SARS-CoV-2 re-introduction from international travellers. medRxiv. 2020. https://doi.org/10.1101/2020.07.24.20161281.

- 25
Bullard J, Dust K, Funk D, Strong JE, Alexander D, Garnett L, Boodman C, Bello A, Hedley A, Schiffman Z, et al. Predicting infectious SARS-CoV-2 from diagnostic samples. Clin Infect Dis. 2020;:ciaa638. https://doi.org/10.1093/cid/ciaa638.

- 26
Robert Koch Institut. Steckbrief zur Coronavirus-Krankheit-2019 (COVID-19). 2020. https://www.rki.de/DE/Content/InfAZ/N/Neuartiges_Coronavirus/Steckbrief.html. Accessed 07 Apr 2020.

- 27
ISARIC. International Severe Acute Respiratory and Emerging Infections Consortium COVID-19 Report: 13 July 2020. 2020. https://media.tghn.org/medialibrary/2020/07/ISARIC_Data_Platform_COVID-19_Report_13JUL20.pdf. Accessed 09 Nov 2020.

- 28
Dreher M, Kersten A, Bickenbach J, Balfanz P, Hartmann B, Cornelissen C, Daher A, Stöhr R, Kleines M, Lemmen SW, et al. The characteristics of 50 hospitalized COVID-19 patients with and without ARDS. Deutsches Ärzteblatt Int. 2020; 117(10):271.

- 29
Khalili M, Karamouzian M, Nasiri N, Javadi S, Mirzazadeh A, Sharifi H. Epidemiological characteristics of COVID-19: a systemic review and meta-analysis. MedRxiv. 2020. https://doi.org/10.1101/2020.04.01.20050138.

- 30
Spiteri G, Fielding J, Diercke M, Campese C, Enouf V, Gaymard A, Bella A, Sognamiglio P, Moros MJS, Riutort AN, et al. First cases of coronavirus disease 2019 (COVID-19) in the WHO European region, 24 January to 21 February 2020. Eurosurveillance. 2020; 25(9):2000178.

- 31
Lavezzo E, Franchin E, Ciavarella C, Cuomo-Dannenburg G, Barzon L, Del Vecchio C, Rossi L, Manganelli R, Loregian A, Navarin N, et al. Suppression of a SARS-CoV-2 outbreak in the Italian municipality of Vo. Nature. 2020; 584(7821):425–9.

- 32
Streeck H, Schulte B, Kuemmerer B, Richter E, Hoeller T, Fuhrmann C, Bartok E, Dolscheid R, Berger M, Wessendorf L, Eschbach-Bludau M, Kellings A, Schwaiger A, Coenen M, Hoffmann P, Noethen M, Eis-Huebinger AM, Exner M, Schmithausen R, Schmid M, Kuemmerer B. Infection fatality rate of SARS-CoV-2 infection in a German community with a super-spreading event. MedRxiv. 2020. https://doi.org/10.1101/2020.05.04.20090076. Accessed 17 June 2020.

- 33
Karagiannidis C, Mostert C, Hentschker C, Voshaar T, Malzahn J, Schillinger G, Klauber J, Janssens U, Marx G, Weber-Carstens S, et al. Case characteristics, resource use, and outcomes of 10 021 patients with COVID-19 admitted to 920 German hospitals: an observational study. Lancet Respir Med. 2020; 8(9):853–62.

- 34
Armstrong R, Kane A, Cook T. Outcomes from intensive care in patients with COVID-19: a systematic review and meta-analysis of observational studies. Anaesthesia. 2020; 75(10):1340–9.

- 35
Sabbadini LL. Primi Risultati dell’ Indagine di Sieroprevalenza SARS-CoV-2. Results from ISTAT. 2020. https://www.slideshare.net/slideistat/primi-risultati-dellindagine-di-sieroprevalenza-sarscov2. Accessed 09 Nov 2020.

- 36
Quartagno M. SARS-CoV-2 IFR Italy based on seroprevalence. IFR for Italy. 2020. https://twitter.com/stats_q/status/1290330189810081792. Accessed 09 Nov 2020.

- 37
Robert Koch Institut. Corona virus disease 2019 (COVID-19) daily situation report of the Robert Koch Institute. Robert Koch Institut. 2020. https://www.rki.de/DE/Content/InfAZ/N/Neuartiges_Coronavirus/Situationsberichte/2020-06-17-en.pdf?__blob= publicationFile. Accessed 19 June 2020.

- 38
Vanella P, Wiessner C, Holz A, Krause G, Moehl A, Wiegel S, Lange B, Becher H. The role of age distribution, time lag between reporting and death and healthcare system capacity on case fatality estimates of COVID-19. MedRxiv. 2020. https://doi.org/10.1101/2020.05.16.20104117. Accessed 24 June 2020.

- 39
Davies NG, Klepac P, Liu Y, Prem K, Jit M, CMMID COVID-19 working group, Eggo RM. Age-dependent effects in the transmission and control of COVID-19 epidemics. medRxiv. 2020. https://doi.org/10.1101/2020.03.24.20043018. Accessed 07 Apr 2020.

- 40
Protezione Civile. COVID-19 Italia - Monitoraggio situazione. 2020. https://github.com/Dpcm-dpc/COVID-19. Accessed 09 Nov 2020.

- 41
Nationale Plattform für geographische Daten. RKI COVID19. 2020. https://npgeo-corona-npgeo-de.hub.arcgis.com/datasets/dd4580c810204019a7b8eb3e0b329dd6_0/data. Accessed 04 Apr 2020.

- 42
GENESIS Online. Bevölkerung: Deutschland, Stichtag, Altersjahre. Statistisches Bundesamt. 2020. https://www-genesis.destatis.de. Accessed 01 Apr 2020.

- 43
GENESIS Online. Bevölkerung: Bundesländer, Stichtag, Altersjahre. Statistisches Bundesamt. 2020. https://www-genesis.destatis.de. Accessed 01 Apr 2020.

- 44
Dietz K. The estimation of the basic reproduction number for infectious diseases. Stat Methods Med Res. 1993; 2(1):23–41.

- 45
Wallinga J, Lipsitch M. How generation intervals shape the relationship between growth rates and reproductive numbers. Proc Biol Sci. 2007; 274(1609):599–604. https://doi.org/10.1098/rspb.2006.3754. Accessed 08 June 2020.

- 46
Heffernan JM, Smith RJ, Wahl LM. Perspectives on the basic reproductive ratio. J R Soc Interface. 2005; 2(4):281–93. https://doi.org/10.1098/rsif.2005.0042. Accessed 07 Apr 2020.

- 47
Diekmann O, Heesterbeek JA, Metz JA. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J Math Biol. 1990; 28(4):365–82. https://doi.org/10.1007/BF00178324. Accessed 07 Apr 2020.

- 48
Diekmann O, Heesterbeek JAP, Roberts MG. The construction of next-generation matrices for compartmental epidemic models. J R Soc Interface. 2010; 7(47):873–85. https://doi.org/10.1098/rsif.2009.0386. Accessed 07 Apr 2020.

- 49
Fine P, Eames K, Heymann DL. “Herd immunity”: a rough guide. Clin Infect Dis. 2011; 52(7):911–6. https://doi.org/10.1093/cid/cir007. Accessed 24 June 2020.

- 50
OECD. Coronavirus: the world economy at risk. OECD Interim Economic Assessments. 2020. https://www.oecd.org/berlin/publikationen/Interim-Economic-Assessment-2-March-2020.pdf. Accessed 24 June 2020.

- 51
Weston D, Hauck K, Amlôt R. Infection prevention behaviour and infectious disease modelling: a review of the literature and recommendations for the future. BMC Public Health. 2018; 18(1):336. https://doi.org/10.1186/s12889-018-5223-1. Accessed 24 June 2020.

- 52
van den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002; 180:29–48. https://doi.org/10.1016/S0025-5564(02)00108-6. Accessed 17 June 2020.

- 53
van den Driessche P, Watmough J. Mathematical Epidemiology In: Brauer F, van den Driessche P, Wu J, Morel J-M, Takens F, Teissier B, editors. Lecture notes in mathematics, vol. 1945. Berlin, Heidelberg: Springer: 2008. p. 159–78. https://doi.org/10.1007/978-3-540-78911-6_6. Accessed 17 June 2020.

- 54
Cori A. EpiEstim: estimate time varying reproduction numbers from epidemic curves. R package version 2.2-3. 2020. https://CRAN.R-project.org/package=EpiEstim. Accessed 09 Nov 2020.

- 55
Bundesministerium für Gesundheit. Chronik zum Coronavirus SARS-CoV-2 - Maßnahmen des BMG. 2020. https://www.bundesgesundheitsministerium.de/coronavirus/chronik-coronavirus.html. Accessed 24 June 2020.

- 56
Desvars-Larrive A, Dervic E, Haug N, Niederkrotenthaler T, Chen J, Di Natale A, Lasser J, Gliga DS, Roux A, Chakraborty A, et al. A structured open dataset of government interventions in response to COVID-19. medRxiv. 2020. https://doi.org/10.1101/2020.05.04.20090498.

- 57
Bundesregierung. Besprechung der Bundeskanzlerin mit den Regierungschefinnen und Regierungschefs der Länder. 2020. https://www.bundesregierung.de/breg-de/themen/coronavirus/besprechung-der-bundeskanzlerin-mit-den-regierungschefinnen-und-regierungschefs-der-laender-1733248. Accessed 24 June 2020.

- 58
Apple. COVID-19 - Mobility Trend Reports - Apple. Mobility data from searchs in Apple map. 2020. https://covid19.apple.com/mobility. Accessed 09 Nov 2020.

- 59
Ritchie H. Coronavirus (COVID-19) testing. Statistics on COVID-19 testing: ourworldindata. 2020. https://ourworldindata.org/coronavirus-testing. Accessed 09 Nov 2020.

- 60
Jarvis CI, Van Zandvoort K, Gimma A, Prem K, Klepac P, Rubin GJ, Edmunds WJ, Auzenbergs M, O’Reilly K, Medley G, Emery JC, Houben RMGJ, Davies N, Nightingale ES, Flasche S, Jombart T, Hellewell J, Abbott S, Munday JD, Bosse NI, Funk S, Sun F, Endo A, Rosello A, Procter SR, Kucharski AJ, Russell TW, Knight G, Gibbs H, Leclerc Q, Quilty BJ, Diamond C, Liu Y, Jit M, Clifford S, Pearson CAB, Eggo RM, Deol AK. Quantifying the impact of physical distance measures on the transmission of COVID-19 in the UK. BMC Med. 2020; 18(1):124. https://doi.org/10.1186/s12916-020-01597-8.

- 61
Guenther T, Czech-Sioli M, Indenbirken D, Robitailles A, Tenhaken P, Exner M, Ottinger M, Fischer N, Grundhoff A, Brinkmann M. Investigation of a superspreading event preceding the largest meat processing plant-related sars-coronavirus 2 outbreak in germany. SSRN J. 2020;:3654517. doi:10.2139/ssrn.3654517.

- 62
Levin AT, Cochran KB, Walsh SP. Assessing the age specificity of infection fatality rates for COVID-19: meta-analysis & public policy implications. Technical report, National Bureau of Economic Research. 2020.

- 63
Caglar MU, Teufel AI, Wilke CO. Sicegar: R package for sigmoidal and double-sigmoidal curve fitting. PeerJ. 2018; 6:4251.

- 64
DIVI-IntensivRegister. DIVI IntensivRegister Tagesreport 04.04.2020. Deutsche Interdisziplinäre Vereinigung für Intensiv- und Notfallmedizin (DIVI). 2020. https://www.divi.de/images/Dokumente/{DIVI}-{IntensivRegister_Tagesreport_2020_04_04}.pdf. Accessed 07 Apr 2020.

- 65
Ruiz-Perez I, Rodriguez-Barranco M, Rojas-Garcia A, Mendoza-Garcia O. Economic crisis and suicides in Spain. Socio-demographic and regional variability. Eur J Health Econ. 2017; 18(3):313–20. https://doi.org/10.1007/s10198-016-0774-5.

- 66
Bartlett DL, Howe JR, Chang G, Crago A, Hogg M, Karakousis G, Levine E, Maker A, Mamounas E, McGuire K, Merchant N, Shibata D, Sohn V, Solorzano C, Turaga K, White R, Yang A, Yoon S. Management of cancer surgery cases during the COVID-19 pandemic: considerations. Ann Surg Oncol. 2020; 27(6):1717–20. https://doi.org/10.1245/s10434-020-08461-2.

- 67
Dorn F, Khailaie S, Stöckli M, Binder S, Lange B, Peichl A, Vanella P, Wollmershäuser T, Fuest C, Meyer-Hermann M. Das gemeinsame Interesse von Gesundheit und Wirtschaft: Eine Szenarienrechnung zur Eindämmung der Corona- Pandemie. ifo Schnelldienst. 2020. https://www.ifo.de/en/publikationen/2020/article-journal/das-gemeinsame-interesse-von-gesundheit-und-wirtschaft. Accessed 05 June 2020.

- 68
Iwasaki A, Grubaugh ND. Why does Japan have so few cases of COVID-19?EMBO Mol Med. 2020; 12(5):12481. https://doi.org/10.15252/emmm.202012481.

- 69
Korth J, Wilde B, Dolff S, Anastasiou OE, Krawczyk A, Jahn M, Cordes S, Ross B, Esser S, Lindemann M, Kribben A, Dittmer U, Witzke O, Herrmann A. SARS-CoV-2-specific antibody detection in healthcare workers in Germany with direct contact to COVID-19 patients. J Clin Virol. 2020; 128:104437.

- 70
Morawska L, Cao J. Airborne transmission of SARS-CoV-2: the world should face the reality. Environ Int. 2020; 139:105730. https://doi.org/10.1016/j.envint.2020.105730.

- 71
Setti L, Passarini F, De Gennaro G, Barbieri P, Perrone MG, Borelli M, Palmisani J, Di Gilio A, Piscitelli P, Miani A. Airborne transmission route of COVID-19: why 2 meters/6 feet of inter-personal distance could not be enough. Int J Environ Res Public Health. 2020; 17(8):2932. https://doi.org/10.3390/ijerph17082932.

- 72
Jones RM, Brosseau LM. Aerosol transmission of infectious disease. J Occup Environ Med. 2015; 57(5):501–8. https://doi.org/10.1097/JOM.0000000000000448.

- 73
Asadi S, Wexler AS, Cappa CD, Barreda S, Bouvier NM, D RW. Aerosol emission and superemission during human speech increase with voice loudness. Sci Rep. 2019; 9:2348. https://doi.org/10.1038/s41598-019-38808-z.

## Acknowledgements

We thank Werner Hofmann and Alexander Kuhlmann for insightful suggestions. We thank Rebecca Ludwig for support in data acquisition and for revising the manuscript.

## Funding

This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement no. 101003480 and the Initiative and Networking Fund of the Helmholtz Association. It was supported by German Federal Ministry of Education and Research for the project CoViDec (FKZ: 01KI20102). The funding bodies had no role in the design of the study; collection, analysis, and interpretation of the results; or writing the manuscript. Open Access funding enabled and organized by Projekt DEAL.

## Author information

### Affiliations

### Contributions

TM and MMH developed the SECIR model and the corresponding methodology, and wrote the paper; SK coded the simulation, automated the analysis of daily incoming data, and contributed to the visualization; TM derived the literature-based parameter set, coded the simulation, and contributed to the accuracy assessment and visualization; SK, TM, SB, and MMH organized the simulation analysis and evaluation; AB and MS coded an independent simulation for Italy and derived the corresponding parameter set; TM and PM derived the formula for the reproduction number; PM and SB contributed to the quality assessment of the results and performed the statistical analysis; PV provided the data for the analysis; BL supported the data acquisition and analysis; SB organized the public repository; SB and MMH supervised the study; all revised the paper. The authors read and approved the final manuscript.

### Corresponding authors

## Ethics declarations

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

## Additional information

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Supplementary Information

**Additional file 1**

Figures S1-S8, Parameter description, Additional details of Italy fitting, and Table 1.

## Rights and permissions

Corrected publication. April 2021 **Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

## About this article

### Cite this article

Khailaie, S., Mitra, T., Bandyopadhyay, A. *et al.* Development of the reproduction number from coronavirus SARS-CoV-2 case data in Germany and implications for political measures.
*BMC Med* **19, **32 (2021). https://doi.org/10.1186/s12916-020-01884-4

Received:

Accepted:

Published:

### Keywords

- SARS-CoV-2
- COVID-19
- Epidemiology
- Modeling
- Non-pharmaceutical interventions
- Reproduction number
- Healthcare usage