Abstract
The COVID-19 pandemic has led to unprecedented strain on intensive care unit (ICU) admission in parts of the world. Strategies to create surge ICU capacity requires complex local and national service reconfiguration and reduction or cancellation of elective activity. Theses measures require time to implement and have an inevitable lag before additional capacity comes on-line. An accurate short-range forecast would be helpful in guiding such difficult, costly and ethically challenging decisions.
At the time this work began, cases in England were starting to increase. Here we present an attempt at an agile short-range forecast based on published real-time COVID-19 case data from the seven National Health Service commissioning regions in England (East of England, London, Midlands, North East and Yorkshire, North West, South East and South West). We use a Monte Carlo approach to model the likely impact of current diagnoses on regional ICU capacity over a 14 day horizon. Our model is designed to be parsimonious and based on plausible epidemiological data from the literature available.
On the basis of the modelling assumptions made, ICU occupancy is likely to increase dramatically in the the days following the time of modelling. If the current exponential growth continues, 5 out of 7 commissioning regions will have more critically ill COVID-19 patients than there are ICU beds within two weeks. Despite variable growth in absolute patients, all commissioning regions are forecast to be heavily burdened under the assumptions used.
Whilst, like any forecast model, there remain uncertainties both in terms of model specification and robust epidemiological data in this early prospective phase, it would seem that surge capacity will be required in the very near future. We hope that our model will help policy decision makers with their preparations. The uncertainties in the data highlight the urgent need for ongoing real-time surveillance to allow forecasts to be constantly updated using high quality local patient-facing data as it emerges.
Introduction
The emergence of a novel COVID-19 coronavirus pandemic has rapidly caused an enormous world-wide medical and socioeconomic impact since the first case emerged on November 16th 2019 [1]. There is a potential for cases of respiratory complications from such infectious outbreaks to over-whelm intensive care unit (ICU) capacity and this was seen previously with the influenza-A [H1N1] pandemic of 2009. The relatively large numbers of COVID-19 patients with hypoxaemic respiratory failure requiring ICU admission [1] for mechanical ventilation and high mortality (compared to seasonal influenza) is of particular concern in the current circumstances. In Northern Italy, an exponential increase in COVID-19 admissions rapidly overwhelmed normal ICU capacity [2] and surge capacity had to be created quickly. It is currently unclear exactly why the ICUs in Italy were so badly affected and whether this will generalise to other countries: Both demographic factors and healthcare system structure are likely to be important. However, it is noteworthy that the UK compares poorly with other high-income countries—including Italy—in terms of available ICU bed capacity per capita [3] and so it is far from inconceivable that a similar pattern may emerge.
Whilst normal wards may be relatively straightforwardly re-purposed from a basic care perspective, this is more difficult for ICUs which are naturally constrained by the need for complex equipment and for the delivery of highly specialised care. Nevertheless, whilst difficult, there are a number of mechanisms by which ICU capacity can be increased in an emergency situation. These range from improving patient flow through a reduction in elective work through changing referral networks to the provision of additional emergency physical capacity for mechanical ventilation—for example using operating theatre ventilators—at the extreme. None of these measures are simple, having complex knock-on effects and requiring changes to infrastructure and staffing which are expensive and cannot be instantaneous to accommodate. As a result all such measures share a common feature of lag time between implementation and realisation of the additional capacity. Furthermore, the implications in staff redeployment and resourcing are significant. Forecasting was therefore essential in guiding such difficult policy decisions in Italy [2]. The explosion in cases seen in Italy means that an early warning of when surge capacity is likely to be required is urgently needed elsewhere.
Epidemiological simulation has previously been successful in predicting the need for surge H1N1 ICU capacity in 2009 [4, 5]. In recent days, a similar simulation model for COVID-19 has been circulating [6] which suggests an overwhelming demand for critical care, with peak demand occurring between May and early June 2020 and lasting 2–3 months depending on various non-pharmacological intervention (NPI) assumptions. Crucially, however, such simulation models do not incorporate up to date data on outbreak start time. The rapid dynamics of an outbreak can be very variable and this means that such simulations are unsuitable for real time short-range forecasting and early warning.
In this paper we use published COVID-19 diagnosis data for England to generate the earliest possible estimates of additional ICU demand due to infections in the coming days based on cautious epidemiological assumptions from the literature. Our emphasis is making an updatable model from what little time-series is available in this ultra-early period, cognisant that this will inevitably involve a number of assumptions where data is not available. Our model predictions account for the latest results from the rapidly developing COVID-19 literature, account for English demographics, and are stratified by the National Health Service (NHS) commissioning regions across the country.
Methods
We begin with COVID-19 diagnoses from England as reported by PHE and matched to NHS commissioning regions [8] as our source data to obtain information on daily cases. We began extracting this data feed on 13/03/20 to give daily case data. Since this data was only publicly available in machine readable form recently, backward extrapolation was necessary for the small number of cases in the week before 13/03/20. This was achieved by scaling the previous country-wide new case estimates– which are available as a time series from the publicly available dashboard [9]– by the average probability of a case falling in a given region over the observed days.
Since the daily incidence of COVID-19 appears exponential (in line with what was observed in Italy [2]), we forecast the likely distributions of new COVID-19 diagnoses by extrapolating from the log of the daily incidence using ordinary least squares over the next 14 days. We chose this simple model which neglects pseudoreplication because, at the time of modelling, we do not have sufficient data to look for the underlying autocorrelation structure in this time series robustly.
Using early data from Verity et al. 2020 [10] (reproduced in Table 1), we estimate the ICU mortality and ICU admission rate per case by standardising to the local population in each NHS commissioning region in England (see Figure 2, obtained from the Clinical Commissioning Group population estimates for mid-2018 [11]). We adjust by age and sex, which have both been shown to strongly correlate with mortality [10]. There is a higher male:female ratio in China than in England, so we use data from the 2017 Chinese census [12] to adjust for this under the assumption that this holds on a univariate basis. We use the ICU admission rate estimates in each region to transform the expected new cases per day into ICU admissions.
Model timeline. Our model relies on a 10-day window for regression, constructed from recent observed data and estimated prior data. The model predicts two weeks into the future from the time of writing.
Regional demographics and expected critical care demand per case of COVID-19 stratified by region, each compared to ICU bed capacity per 100,000 people. Population is divided into age categories, and percentage of cases requiring ICU is divided into expected percentage of survivors and non-survivors. The numerical data can be seen in Table 1.
As new cases will take time to arrive in the ICU, we model the delay between symptoms (which we assume to be concurrent with diagnosis) and ICU admission with a normal dis-tribution (µ = 10, σ = 3.5). This is in line with multiple early reports from China [13, 14,15, 16], which indicate a consistent pattern of disease progression with hospital admission occurring around 8 days (concurrent with the onset of dyspnoea), Acute Respiratory Distress Syndrome (ARDS) following on day 9 when it occurs, and patients admitted to the ICU for mechanical ventilation around day 10. Ten days was also used in the models published in Ferguson et al. [6]. We model the ICU admission date according to this delay distribution which is considered to be normally distributed.
Bed occupancy also requires estimates of the length of stay for patients currently in the ICU. The most up-to-date information indicates a median length of stay of approximately 8 days, with a wide interquartile range and positive skew [14, 15]. We model this with a gamma dis-tribution (α = 8, β = 1).
We use a Monte Carlo simulation with 100 samples to link the expected daily incidence distribution to the expected excess ICU bed occupancy due to COVID-19. This involves sampling from the delay and length of stay distributions as well as the daily incidence distribution. We represent bed occupancy as a fraction of the total number of ICU beds in the commissioning region based on data from a snapshot in December 2019 [7].
Results
Figure 3 shows the expected number of new COVID-19 patients per day assuming the extra-polated current exponential trend in England. We note that all R2 values are, on average, above 0.9 implying that, for now, new case growth in England is well modelled by exponential growth. In London, the R2 value for the fit is particularly high and the 95% confidence interval is tight demonstrating this region is likely to already be in a period of exponential growth of new cases. However, the fit is less reliable in the North East and Yorkshire, a region wheres we have separately confirmed the trend of observed new cases to have the weakest correlation with the other regions. The log-linear fit produces a strong R2 value in the non-London populous areas (the Midlands and the North West), strengthening the argument for exponential growth in areas where the virus has gained a foothold. Irrespective of how many of these cases translate into ICU patients, even the lower confidence bounds indicate more than 1000 new cases per day in all of the NHS commissioning regions in England within a week of writing.
Projected new COVID-19 patients in the seven National Health Service commissioning regions in England. Exponential projections were calculated by an ordinary least squares fit on both estimates of previous values and reported recent values in log-space.
Figure 4 shows the projected additional ICU occupancy due to COVID-19 from our model over a 14 day horizon for each of the NHS commissioning regions as a percentage of total ICU beds normally commissioned.
Projected regional COVID-19 ICU occupancy as a percentage of regional capacity the seven National Health Service commissioning regions in England.
All source code has been made available at https://github.com/ariercole/Cambridge_COVID-19_ICU [17]. In addition to our analysis and open-source code, an interactive version of our model is available at http://covid19icu.cl.cam.ac.uk. As documented cases in the literature evolve, we hope clinical and policy decision makers will be able to experiment based on their region or the statistics demonstrated by their cohort. Given the highly dynamic situation at the time this work was carried out, with model data changing on a daily basis, minimising model development time was crucial. Using agile project management methodologies, we were able to develop a working model, documentation and web-implementation in less than one week.
Discussion
It is crucial to appreciate that both model and parameter uncertainties are inevitable, particularly when predicting the behaviour of a novel virus in a new population, and this may radically affect our forecasts. However, we set out to provide the earliest possible data-driven forecast and must therefore explicitly accept the limita-tions of the data available at the time. Our approach has been to keep the model as parsimonious as possible with what we hope are plausible parameter estimates from the existing literature to give a ‘ball park’ estimate of early surge needs. Whilst we believe this approach has merit in that our findings suggest that waiting for better data may mean forecasting is too late, there are some specific limitations which need to be discussed.
Crucially, we have used Public Health England published data for case ascertainment as this is publicly available and seems reasonable as the earliest available starting point. We recognise that this data is potentially flawed and does not recognise all cases within the wider population. However, for our needs we are specifically focused on the severe end of cases where, arguably, ascertainment bias is likely minimal due to the roll-out of routine testing for all critical care patients by PHE as part of the ‘COVID Hospitalisation in England Surveillance System - CHESS’ [18]. A persistent issue however is that case ascertainment and timing after symptom onset is unlikely to be uniform between healthcare settings or even over time in a single setting as the surveyed population varies. Thus the compatibility between our data and the literature is uncertain. Nevertheless, this will always be a limitation of any early modelling including that used in reference [10].
The ICU bed estimates are assumed to be fixed at the December 2019 baseline. Whilst this may vary and modify our results in detail, we do not believe that this is likely to change so greatly as to alter our conclusions. Furthermore, we assume that each region behaves as a perfect ‘pool’ of beds in a way that they may be optimally allocated. Of course this is not necessarily the case, since inter-hospital ICU-to-ICU transfers are not always feasible for both operational and clinical reasons, so the patient load is not necessarily uniformly distributed. We do not have more granular data available to model this, but it is worth considering that it is very possible for an individual hospital to reach a locally critical capacity even before the region as a whole does. In this sense our predictions represent a ‘best case’ scenario.
We have forecast the percentage COVID-19 bed requirements in isolation. In reality, the ICU must continue to provide ‘business as usual’ care for other types of patients. Since UK bed occupancy is typically greater than 80% and may frequently exceed 100%, it is clearly not the case that all open beds can simply be re-allocated for COVID-19 patients. Furthermore, we have assumed that all adult critical care beds can be used for level 3 or mechanically ventilated ICU patients and, operationally, this may also not be possible. Some specialist ICUs may not be able to reconfigure at all. Thus the precise percentage of additional COVID-19 patients which will actually exhaust routine capacity will vary from unit to unit, particularly in ICUs with a substantial post-operative elective surgical workload. Overall, it seems unlikely that any ICU could routinely accommodate the numbers forecast by our model without provision of surge capacity.
Conclusions
Early warning of an impending need for ICU surge capacity is crucial if there is to be sufficient time to re-configure services to provide additional capacity. We have shown that ultra-early data can be used to make time-sensitive forecasts of ICU occupancy. We show, subject to our assumptions, that it is credible that ICU requirements may become challenging within weeks. There remains a significant degree of uncertainty in the predictions due both to limitations of the reporting data and also modelling assumptions. This emphasises the need for the collection of real-time patient-facing local data by initiatives such as CHESS [18] and a dynamic approach to improving models as new data becomes available.
Data Availability
Source data is already publicly available.
Declarations of interest statement
None declared.
Acknowledgments
The authors would like to thank Ronan O’Leary and Isobel Ramsay for useful discusions as well as Mark Cresham for rapidly procuring compute facilities for our online model.
