Forecasting ultra-early intensive care strain from COVID-19 in England

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\todo{last thing to do}. 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 COVID-19 pandemic has rapidly caused an enormous worldwide medical and socioeconomic impact since the first case emerged on November 16th 2019 [1]. Although a self-limiting illness for most, the percentage of COVID-19 patients with hypoxaemic respiratory failure requiring ICU admission for mechanical ventilation translates into large numbers which may challenge healtcare provision. In Northern Italy, an exponential increase in COVID-19 admissions rapidly overwhelmed normal ICU capacity [2] and surge capacity had to be created quickly. The exact reason for the sudden need for ICU surge capacity in Italy, and whether this will generalise to other countries, is unclear, but both demographic factors and healthcare system structure are likely to be important. Notably, UK availability of ICU beds per capita compares poorly with other highincome countries-including Italy [3]. Whilst standard acute wards may be re-purposed easily, creating ICU capacity is constrained by the need for complex equipment and the delivery of highly specialised medical and nursing care. Nevertheless, there are mechanisms by which ICU capacity could be increased in an emergency, including facilitating ICU discharge of recovering patients through liberating downstream beds by reducing elective work, by stopping elective work likely to require ICU admission, or by changing referral policies. Such changes are likely to be relatively quick to implement but have important repercussions for normal healthcare provision and ideally would not be instituted until absolutely necessary. Furthermore, since ICUs in the UK typically undertake a high proportion of emergency work, much of which will continue despite the pandemic, and since occupancy is typically well above 80%, such strategies are likely to result in only a relatively modest increase in capacity. A greater increase in emergency physical capacity for mechanical ventilation could be achieved by opening new level 3 beds with additional equipment (e.g. operating theatre ventilators). This requires significant changes to infrastructure, processes, or sta ng and is therefore logistically complex, expensive, and most importantly slow to implement. Forecasting was therefore essential in guiding such di cult policy decisions in Italy [2]. The explosion of cases seen in Italy means that an early warning of need for surge capacity is likely to be required in other countries including England.
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 described [6], which suggests an overwhelming demand for critical care, with a peak occurring between May and early June 2020 and lasting 2 to 3 months depending on non-pharmacological intervention (NPI) assumptions. Such models are very useful but do not incorporate up-to-date data. Such projected timescales are unlikely to be reliable and are unsuitable for real-time surveillance 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 data from the literature and under the assumption that the current increase in cases represents the exponential phase of an outbreak rather than a change in ascertainment. Our emphasis is on making an updatable model from the little timeseries data that are available in this ultra-early period, with the understanding that assumptions are necessary where data are unavailable. Our

Methods
We used COVID-19 diagnoses from England as reported by Public Health England (PHE) and matched to NHS commissioning regions [8] as our source data to obtain information on daily cases (reproduced in Figure 1). We started to extract this data feed on 13/03/20 to give daily case data. Figure 1: Model timeline. Our model relies on a n-day window for regression, beginning 13/03 with recent observed data. Themodel predicts two weeks into the future from the time of writing. Figure 2: Regional demographics and expected critical care demand per case of COVID-19, stratified by region and 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.
We assumed that the daily incidence of COVID-19 can be modelled as an exponential growth (in line with what was observed in Italy [2]). Therefore, we forecast the likely distributions of new COVID-19 diagnoses over the next 14 days by using an ordinary least squares fit to linearly extrapolate from the logarithm of the cumulative cases. We di↵erentiated the cumulative model to obtain a daily incidence model i.e. we multiply by the exponent ( d dx ⇥ Ae Bx ⇤ = BAe Bx ), which stabilises the exponent fitting (B).
Using early data from Verity et al. 2020 [9] (reproduced in Table 1), we estimated the ICU admission rate 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 [10]). We adjust by age and sex, which 3 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not peer-reviewed)
The copyright holder for this preprint . https://doi.org/10.1101/2020.03. 19    Overall, Figure 4 shows predictions that will challenge ICU capacity in all areas, particularly in populous regions. London and the Midlands will gain a number of COVID-19 patients equi-valent to 100% standard capacity within 14 days for slightly di↵erent reasons. In London, despite the higher number of ICU beds per capita, the higher case incidence per capita will overwhelm ICU capacity. Whereas in the Midlands, due to an increased average age, a higher percentage of ICU admissions is predicted (see Table 1), which combined with a higher number of cases than other regions will lead to exhausted capacity. As normal occupancy is often above 80%, even the additional minimum capacity due to COVID-19 of 20% in other regions will prove to be di cult to 5 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not peer-reviewed)
The copyright holder for this preprint . https://doi.org/10.1101/2020.03. 19.20039057 doi: medRxiv preprint Figure 4: Projected regional COVID-19 ICU occupancy as a percentage of regional capacity in the seven National Health Service commissioning regions in England. The figure illustrates the large regional variation in projected ICU occupancy, with London and the Midlands reaching near 100% additional capacity due to COVID-19 within 14 days, whereas areas such as the South West and South East will only have approximately 20% of their beds occupied by COVID-19 patients.

Discussion
In this study we demonstrate how publicly available data can be rapidly combined to dynamic-6 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not peer-reviewed)
The copyright holder for this preprint . https://doi.org/10.1101/2020.03. 19.20039057 doi: medRxiv preprint ally model short-term ICU requirements in light of the emerging COVID-19 pandemic. Our data suggests that traditional ICU capacity could be rapidly consumed over a period of approximately 14 days from the time of modelling/writing, such figures hide substantial regional heterogeneity, with London and the Midlands demonstrating the most rapid growth. The key message from our paper varies with the temporal position of the reader. Within the current acute stage of this pandemic, we hope our findings provide a framework to facilitate the rapid development and deployment of additional ICU surge capacity as was required within days of the Italian pandemic [2].
In the future, we feel that our methodological approach (use of emerging epidemic curve data, Monte Carlo simulation, and population standardisation), as well as our simultaneous deployment of an interactive and dynamically updated web-tool have merit in modelling and communicating any future dynamically emerging event.
However, we do recognise that as with all models our is subject to significant assumptions and limitations. Both model and parameter uncertainties are inevitable, particularly when predicting the behaviour of a novel virus in a new population, and this may radically a↵ect our forecasts.
Nevertheless, we set out to provide the earliest possible data-driven forecast and therefore explicitly accept the limitations of the data that were available to us at the time. Our approach has been to keep the model as parsimonious as possible, with parameter estimates from the existing literature, to give a rough guide to early surge needs. With the small amount of available data, there are some specific limitations which need to be discussed. Our use of an updated, online resource we hope is a novel way of allowing interested parties to interact with the mathematical assumptions and formulate their own opinions.
We have used PHE published data for case ascertainment. We recognise that this data is potentially flawed and does not recognise all cases within the wider population. However, being as our explicit aim is to model the most severe cases, the e↵ect of any ascertainment bias on our findings should be minimised due to the roll-out of routine testing for all critical care patients as part of the "COVID Hospitalisation in England Surveillance System-CHESS" [18]. Despite this, case ascertainment and delay between symptom onset and ICU admission is unlikely to be uniform between populations, healthcare settings, or discrete time points. Specifically, at the time of writing, test capacity for COVID-19 in the UK has led to stringent requirements for testing with only cases requiring hospitalisation being routinely tested [16]. By incorporating a 'delay to ICU' parameter within our model, that can take on negative values (thus indicating a diagnosis after ICU admission), we hope to have replicated the anecdotal experience of clinicians in these early phases of the pandemic; that a substantial proportion of patients are being con- forecast numbers as these data inform our standardised estimates of ICU cases from within each population [9]. Similarly, as already intimated, UK case definition has evolved over time. If the upturn in COVID-19 cases seen at the time of writing is largely driven by increased or altered ascertainment rather than a true rise in cases, 7 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not peer-reviewed)
The copyright holder for this preprint .
then this would render our model overly pessimistic. Nevertheless, this will always be a limitation of any early modelling [9].
Although we have resolved ICU admissions to the level of a commissioning region, this assumes that each region behaves as a homogeneously allocated 'pool' of ICU beds, which is not necessarily true because inter-hospital ICU-to-ICU transfers may not be feasible for both operational and clinical reasons. We do not have more granular data available, but it is worth considering the likelihood that individual hospitals will reach critical capacity before the whole region. This has already been seen within the pandemic at certain London hospitals [19]. In this sense our predictions represent a 'best case' scenario and cannot be used for decision making at the level of an individual unit.
Building on this theme, we forecast the percentage of COVID-19 bed requirements in isolation, figures which do not reflect competing ICU burden. Since UK bed occupancy is typically greater than 80% and may frequently exceed 100% [7], clearly not all beds can simply be re-allocated for COVID-19 patients. This may be especially true for specialist ICUs (e.g. neuroscience or cardiac) which may not be able to entirely reconfigure, or high-dependency beds which may not be able to provide mechanical ventilation routinely. Furthermore, we assumed that all adult critical care bed spaces can be used for mechanically ventilated ICU patients, which operationally may not be possible for a variety of reasons including equipment and sta↵ availability. Thus, the precise percentage of additional COVID-19 patients that will actually exhaust routine capacity will vary from unit to unit, particularly in ICUs with a substantial postoperative elective surgical workload. Finally, front-line clinicians reviewing our results in these early days of the pandemic will highlight the disconnect between our estimated 'confirmed' cases and the substantial operational workload that is being created by treating all suspicious cases of respiratory failure as COVID-19 until proven otherwise.
At the time of writing, the potential challenge that faces the health system (including intensive care) is essentially without precedent and has resulted in social and political interventions that have never been seen in peacetime. As such, there is little relevant precedent to judge the utility of approaches such as ours against. We hope that regardless of the assumptions inherent within our model (which we feel are based on best available data) that the overarching message of the need for urgent surge capacity is timely. Finally, we hope that our transparency in releasing our code, model assumptions, and results in a dynamic web format embodies a new paradigm in rapidly developed and reviewed science that persists beyond the current pandemic.

Conclusions
Early warning of an impending need for ICU surge capacity is crucial if there is to be sufficient time to re-configure services. 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 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.
Declarations of interest statement None declared.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
(which was not peer-reviewed) The copyright holder for this preprint . https://doi.org/10.1101/2020.03. 19.20039057 doi: medRxiv preprint thank Chris Fryer for crowd-sourced peer review of the manuscript and code which have been invaluable.