Executive Summary: The first wave has ended in Sri Lanka on 2020-11-23 (Day 301). We expect a total of 999,571,387 cases in the first wave of the COVID-19 pandemic in Sri Lanka, this is an additional 999,548,886 cases compared to current case total. New cases have peeked on 2020-11-07 with 461 new cases. Our models assume that after end of the current wave around 446 new cases will appear on an average day, since SARS-CoV-2 is now endemic. Quick, effective, and differentiated mitigation measures need to be in place to prevent a new outbreak.

1 New Cases

According to the Johns Hopkins University global COVID-19 data set [1] a total of 22,501 COVID-19 cases were reported in Sri Lanka on 2020-11-27. This number of cases increased by 2.5 percent over the cumulative total cases known one day earlier. At this rate of cases double every 55 days and the basic reproduction rate \(R_0\) in Sri Lanka can be estimated at 1.1, an estimate of the average number of other people infected by a single infected person. Figure 1.1 shows the time line of new cases reported in the country (black points) together with a five day moving average (purple line) that smoothes daily fluctuations in the reporting that are due to fewer capacities on weekends, holidays and other reporting mishaps.

New cases and five day moving average.

Figure 1.1: New cases and five day moving average.

1.1 Model for new cases (DSM)

New cases peek after a certain time of growth that follows a sigmoidal function and then recedes to a lower level as the epidemic wave fades out and enters an endemic state. This behavior can be modelled with a double-sigmoid function \(f_{base}\), that is obtained by the product of two sigmoidal functions [2].

\[ f_{base}(t) = \frac{1}{ 1+e^{-\alpha_{i} * (t - t_{i})} } * \frac{1}{ 1+e^{-\alpha_{d} * (t - t_{d})} }\]

This function has a maximum \(f_{max} = max(f_{base}(t_{max}))\) at \(t_{max}\).

The model can obtain predictions related to an observed maximum of new cases \(I_{max_{d}}\) and an endemic number of daily cases \(I_{final_{d}}\) as an asymptote (that is ideally 0). Let \(f_{d}(t)\) be a piecewise function that switches between the growth and decay phase at \(t_{max}\), such that

\[\begin{aligned}f_{d}(t) = \left\{ \begin{array}{cc} \frac{I_{max_{d}}}{f_{max}} * f_{base}(t) & \hspace{5mm} t \leq t_{max} \\ \frac{I_{max_{d}} - I_{final_{d} }}{f_{max}} * f_{base}(t) + I_{final_{d}} & \hspace{5mm} t > t_{max} \\ \end{array} \right. \end{aligned}\]

where \(t_{i}\) is the midpoint and \(\alpha_{i}\) the slope of the increase of cases and \(t_{d}\) is the midpoint and \(\alpha_{d}\) is the slope of the decrease in cases. \(I_{max_{d}}\) the peak of new cases and \(I_{final_{d}}\) the asymptotic final value of new cases at the end of the infections wave. Figure 1.2 shows the current double-sigmoidal model ( orange line) that is fitted for the new case time series of Sri Lanka, where the vertical dashed line highlights the midpoint \(t_i\) of the growth phase, the vertical solid line indicates the peak of new infections \(t_max\) and the vertical dotted line the midpoint \(t_d\) of the decay phase. The dashed horizontal line marks the estimated asymptote of endemic new cases \(I_{final_{d}}\) as this infection wave ends. We pragmatically assume that the first wave ends, when 99 percent of the decrease has happened. The dashed green line marks this date.

Double sigmoidal model (DSM) for new cases.

Figure 1.2: Double sigmoidal model (DSM) for new cases.

1.2 Reliability of the DSM

Since the model requires six parameters, it also requires more data points. Hence, it is only feasible to fit this model in later phases of the pandemic (ideally after the peak of daily infections occured) and the model is highly sensitive to new data, particularly with respect to finding the final asymptote\(I_{final_{d}}\)as the following comparison of the model shows. Nevertheless the model renders interesting characteristica of the past developments.

Model for Daily Increases Current Estimate Last Week Estimate Delta (absolute)
\(I_{max_{d}}\) Peak new cases 461 450 11
\(I_{final_{d}}\) Endemic new cases 446 416 29
\(\alpha_{i}\) Slope of growth phase 0.15 0.13 0.02
\(t_{i}\) Midpoint of growth phase 2020-10-16 (Day 264) 2020-10-16 (Day 264) 0
\(t_{max}\) Peak day 2020-11-07 (Day 286) 2020-11-04 (Day 283) 3
\(\alpha_{d}\) Slope of decay phase 0.59 0.04 0.55
\(t_{d}\) Midpoint of decay phase 2020-11-15 (Day 294) NA (Day 315) -21
First wave ends 2020-11-23 (Day 301) 2020-11-23 (Day 301) 0

2 Cumulative cases

The cumulative total of COVID-19 cases in Sri Lanka is the number of all COVID-19 cases known up to yesterday and includes all patients already recovered from the infection. This corresponds to the Infected (I) number in the standard epidemiological SIR model. Figure 2.1 shows the data points of the time series together with a Gompertz model for the cumulative cases (red line) as well as the cumulated double-sigmoidal model (orange line). The dashed horizontal lines show the expected maximum number of cases of the Gompertz model ( the cumulated double-sigmoidal model does not have an asymptote ).