Summary
Background Hospitals need to plan for the surge in demand in each state or region in the United States and the European Economic Area (EEA) due to the COVID-19 pandemic. Planners need forecasts of the most likely trajectory in the coming weeks and will want to plan for the higher values in the range of those forecasts. To date, forecasts of what is most likely to occur in the weeks ahead are not available for states in the USA or for all countries in the EEA.
Methods This study used data on confirmed COVID-19 deaths by day from local and national government websites and WHO. Data on hospital capacity and utilisation and observed COVID-19 utilisation data from select locations were obtained from publicly available sources and direct contributions of data from select local governments. We develop a mixed effects non-linear regression framework to estimate the trajectory of the cumulative and daily death rate as a function of the implementation of social distancing measures, supported by additional evidence from mobile phone data. An extended mixture model was used in data rich settings to capture asymmetric daily death patterns. Health service needs were forecast using a micro-simulation model that estimates hospital admissions, ICU admissions, length of stay, and ventilator need using available data on clinical practices in COVID-19 patients. We assume that those jurisdictions that have not implemented school closures, non-essential business closures, and stay at home orders will do so within twenty-one days.
Findings Compared to licensed capacity and average annual occupancy rates, excess demand in the USA from COVID-19 at the estimated peak of the epidemic (the end of the second week of April) is predicted to be 9,079 (95% UI 253–61,937) total beds and 9,356 (3,526–29,714) ICU beds. At the peak of the epidemic, ventilator use is predicted to be 16,545 (8,083–41,991). The corresponding numbers for EEA countries are 120,080 (119,183–121,107), 32,291 (32,157– 32,425) and 28,973 (28,868–29,085) at a peak of April 6. The date of peak daily deaths varies from March 30 through May 12 by state in the USA and March 27 through May 4 by country in the EEA. We estimate that through the end of July, there will be 60,308 (34,063–140,381) deaths from COVID-19 in the USA and 143,088 (101,131–253,163) deaths in the EEA. Deaths from COVID-19 are estimated to drop below 0.3 per million between May 4 and June 29 by state in the USA and between May 4 and July 13 by country in the EEA. Timing of the peak need for hospital resource requirements varies considerably across states in the USA and across regions of Europe.
Interpretation In addition to a large number of deaths from COVID-19, the epidemic will place a load on health system resources well beyond the current capacity of hospitals in the USA and EEA to manage, especially for ICU care and ventilator use. These estimates can help inform the development and implementation of strategies to mitigate this gap, including reducing non-COVID-19 demand for services and temporarily increasing system capacity. The estimated excess demand on hospital systems is predicated on the enactment of social distancing measures within three weeks in all locations that have not done so already and maintenance of these measures throughout the epidemic, emphasising the importance of implementing, enforcing, and maintaining these measures to mitigate hospital system overload and prevent deaths.
Funding Bill & Melinda Gates Foundation and the state of Washington
Introduction
The Coronavirus Disease 2019 (COVID-19) pandemic started in Wuhan, China, in December 20191 and has since spread to the vast majority of countries.2 As of April 16, twelve countries have recorded more than a thousand deaths: Italy, USA, Spain, France, UK, Iran, China, Netherlands, Germany, Belgium, Canada, and Switzerland. COVID-19 is not only causing mortality but is also putting considerable stress on health systems, with large case numbers and many patients needing critical care including mechanical ventilation. Estimates of the potential magnitude of COVID-19 patient volume – particularly at the local peak of the epidemic –are urgently needed for USA and European hospitals still early in the epidemic to effectively manage the rising case load and provide the highest quality of care possible.
COVID-19 scenarios and forecasts have largely been based on mathematical compartmental models that capture the probability of moving between susceptible, exposed, and infected states, and then to a recovered state or death (SEIR models). Many SEIR or SIR models have been published or posted online.3–20 In general, these models assume random mixing between all individuals in a given population. While results of these models are sensitive to starting assumptions and thus differ between models considerably, they generally suggest that given current estimates of the basic reproductive rate (the number of cases caused by each case in a susceptible population), 25% to 90% of the population could eventually become infected unless mitigation measures are put in place and maintained.6,20 Based on reported case-fatality rates, these projections imply that there would be millions of deaths in the USA and Europe due to COVID-19. Individual behavioural responses and government-mandated social distancing (school closures, non-essential service closures, and shelter-in-place orders), however, can dramatically influence the course of the epidemic. As of April 14, 2020, for Wuhan City in China – and also for at least 12 additional regions in Italy (Liguria, Lombardia, Emilia-Romagna, Marche, Lazio, Campania), Spain (Community of Madrid, Castile and Leon, Catalonia, Navarre), and the USA (King County, Snohomish County) – strict social distancing has led to the peak of the first wave of the epidemic, implying that the effective reproduction number (Reffective) has dropped below unity in these settings. Planning tools based on SEIR models provide high-level information across populations. Few of these planning models have forecasted peaks in deaths or cases and subsequent declines. Using reported case numbers and models based on those for health service planning is also not ideal because of widely varying COVID-19 testing rates and strategies. For example, countries such as Germany, Iceland, and South Korea have undertaken widespread testing, while in the USA and elsewhere, limited test availability has led to largely restricting testing, particularly early in the epidemic, to those with more severe disease or those who are at risk of serious complications.
An alternative strategy is to focus on modelling the empirically observed COVID-19 population death rate curves, which directly reflect both the transmission of the virus and the infection-fatality rates in each specific community. Deaths are likely more accurately reported than cases in settings with limited testing capacity, where tests are usually prioritised for the more severely ill patients. Hospital service need is likely to be highly correlated with deaths, given predictable disease progression probabilities by age for severe cases. In this study, we use statistical modelling to implement this approach and derive state-specific and country-specific forecasts with uncertainty for deaths and for health service resource needs and compare these to available resources in the USA and countries in the European Economic Area (EEA). This model is regularly updated to incorporate new data for the location of interest as well as data from other locations.
Methods
The modelling approach in this study is divided into four components: (i) identification and processing of COVID-19 data; (ii) statistical model estimation for population death rates as a function of time since the death rate exceeds a threshold in a location; (iii) predicting time to exceed a given population death threshold in locations early in the pandemic; and (iv) modelling health service utilisation as a function of deaths. Additional information on the determination of hospital resource utilisation and capacity is provided in Appendix A; details on curve fitting methods, quantification of uncertainty, and a full specification of the statistical model are available in Appendix B. This study complies with the Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) statement.21
Data identification and processing
Local government, national government, and WHO websites, and third-party aggregators22–26 were used to identify data on confirmed COVID-19 deaths by day of death at the first administrative level (state or province, hereafter “admin 1”). Data on licensed bed and ICU capacity and average annual utilisation by location were obtained from a variety of sources for most countries to estimate baseline capacities; observed COVID-19 utilisation data were obtained for a range of countries and USA states providing information on inpatient and ICU use or were imputed from available resources (Appendix A). Other parameters were sourced from the scientific literature and an analysis of available patient-level data. Age-specific data on the relative population death rate by age are available from China,28 Italy,29 South Korea,30 the USA,31,32 Netherlands,33 Sweden,34 and Germany23 and show a strong relationship with age (Figure 1).
Using the average observed relationship between the population death rate and age, data from different locations can be standardised to the age structure using indirect standardisation (Appendix B). For the estimation of statistical models for the population death rate, only admin 1 locations with an observed death rate greater than 0.31 per million (exp(−15)) were used. This threshold was selected by testing which threshold minimised the variance of the slope of the death rate across locations in subsequent days.
Government declarations were used to identify the day that different jurisdictions implemented various social distancing policies (school closures, closures of non-essential services focused on bars and restaurants, stay-at-home or shelter-in-place orders, and the deployment of severe travel restrictions) following the New Zealand government COVID-19 alert schema.35 Data on timings of interventions were compiled by checking national and state governmental websites, executive orders, and newly initiated COVID-19 laws, and cross-referencing other policy compilation resources (see Supplementary Information). Covariates of days with expected exponential growth in the cumulative death rate were created using information on the number of days after the death rate exceeded 0.31 per million that six different social distancing measures were mandated by local and national governments: school closures, partial non-essential business closures, complete non-essential business closures, restricting group gatherings, stay-at-home recommendations, and severe local travel restrictions including public transport closures. To derive weighting schemes for each of the social distancing mandates, we determined the effect of social distancing measures on mobility data published by Google (average of retail, workplace, and transit mobility dimensions),36 Descartes Lab (distance travelled)37 and Safegraph (time spent at home)38 using random effects regression where the dependent variable was the log of mobility measures with social distancing measures as a series of dummy variables. The three different weighting schemes were used to create covariates for an ensemble of three models (Appendix B, section 5). For locations that have not yet implemented all of the closure measures, we assumed that the remaining measures will be put in place within 3 weeks. This lag between reaching a threshold death rate and implementing more aggressive social distancing was combined with the observed period of exponential growth in multiple locations that reached their peak after Level 4 social distancing from the New Zealand alert schema35 was implemented, adjusted for the median time from incidence to death. For ease of interpretation of statistical coefficients, this covariate was normalised so the value for Wuhan was 1.
Statistical model for the cumulative death rate
We developed a curve-fitting tool to fit a nonlinear mixed effects model to the available admin 1 cumulative death data. See Appendix B: Curvefit Tool and Analyses for greater detail. The cumulative death rate for each location is assumed to follow a parametrised Gaussian error function: where the function Ψ is the Gaussian error function (written explicitly above), p controls the maximum cumulative death rate at each location, t is the time since death rate exceeded exp(−15), ß (beta) is a location-specific inflection point (time at which rate of increase of the daily death rate is maximum), and α (alpha) is a location-specific growth parameter. Other sigmoidal functional forms (alternatives to Ψ) were considered but did not fit the data as well. Data were fit to the log of the death rate in the available data, using an optimisation framework described in Appendix B. For data-rich cases, we also developed linear curve fitting extension, where after a Gaussian curve in daily death is obtained, we fit the data to a weighted combination (with constraints on weights) of such curves propagated forward and backward in time. The resulting models can capture more complex behavior in the data.
An ensemble of three models was used to produce the estimates. In all models, we parametrised the time-axis shift parameter beta to depend on a covariate based on time from when the initial ln(death rate) exceeds exp(−15) to the implementation of social distancing. The models differed by the definition of the social distancing covariate. In each model, the value of the covariate multiplier was obtained by fitting a joint model on all the locations that were considered to have peaked; that is, the generalisable information from these locations was the impact that social distancing had on the time to reach the inflection point. Using 13 locations where peak deaths had occurred as of April 14, 2020 – China (Wuhan City), Italy (Liguria, Lombardia, Emilia-Romagna, Marche, Lazio, Campania), Spain (Community of Madrid, Castile and Leon, Catalonia, Navarre), and the USA (King County, Snohomish County) – we fit mixed effects models to get the mean and variance of the relationship between the social distancing covariates and the peak time, and used this information to build priors for location-specific estimates.
We use hospitalization data to generate additional short-term predicted deaths (pseudo-data). On average, the time between hospitalization and death is 8 days. Using location-specific hospitalization data which has more than 10 deaths, we estimate the ratio of cumulative deaths to cumulative hospitalizations up to 8 days in the past. We use this ratio to generate pseudo-data for 8 days, and incorporate this pseudo-data into the CurveFit model. Details are given in Section 11 of Appendix B.
For locations with fewer than 18 days, we use the following analysis. For each type of model (based on definition of the covariate), we considered both “short-range” and “long-range” variants, to explain existing data and forecast long-term trends, respectively. In the former case, covariate multipliers could deviate from those estimated using peaked locations, while in the latter, the joint model fit from peaked locations had a larger impact on the final covariate multiplier. The two remaining parameters (not modelled using covariates) were allowed to vary among locations to fit location-specific data. Uncertainty for every model was obtained using the predictive validity framework that analyses errors in predicting out-of-sample observations. Using these methods, we obtain model realisations using draws, for both short- and long-term models across the forecast horizon. We then obtain forecasts that linearly interpolate between short-term and long-term models, with next days closely following short-term models and long-term forecasts following long-term models. Finally, we ensemble these draws across the model types (based on the definition of the social distancing covariate).
For locations with 18 or more days, we first fit a long-term model, borrowing strength from peaked locations and obtaining location-specific representative daily deaths Gaussian curves. We then fit a linear combination of 13 of the inferred Gaussian curves from the long-term model, placed two days apart (12 days back from the inferred peak to 12 days forward of the inferred peak). We then ensemble across draws for different model types. See Appendix B (section 11) for full details.
The dataset age-standardised to the age-structure of California is shown in Figure 2.
Time to threshold death rate
All states except Wyoming have deaths greater than 0.31 per million (e-15) and more than 2 deaths and were included in the model estimation along with data on 66 other admin 1 locations. For other USA states or locations in the EEA, we estimated the expected time from the current case count to reach the threshold level for the population death rate model. Using the observed distribution of the time from each level of case count to the threshold death rate for all admin 1 locations with data, we estimated this distribution. We used the mean and standard deviation of days from a given case count to the future threshold death rate to develop the probability distribution for the day each state will cross over the threshold death rate, and then we applied the death rate epidemic curve after crossing the threshold.
Hospital service utilisation microsimulation model
From the projected death rates, we estimated hospital service utilisation using an individual-level microsimulation model – additional details are provided in Appendix A. We simulated deaths by age using the average age pattern (Figure 1). For each simulated death, we estimated the date of admission using the median length of stay for deaths from available data (six days). Simulated individuals requiring admission who were discharged alive were generated using the location-specific ratios of admissions to deaths; where location-specific ratios were not available we used the EEA pooled estimate for other EEA countries and the USA pooled estimate for other USA states. An age pattern of the ratio was based on available data (Appendix A). The age-specific fraction of admissions requiring ICU care was based on data from the USA. The fraction of ICU admissions requiring invasive ventilation was estimated as 85%. To determine daily bed and ICU occupancy and ventilator use, we applied median lengths of stay of eight days for those not requiring ICU care and discharged alive and 20 days for those admissions with ICU care, with 13 of those days in the ICU.39
Role of the funding source
The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The authors had access to the data in the study and the final responsibility to submit the paper.
Results
By aggregating forecasts across location, we determined the overall trajectory of expected health-care needs in different categories and deaths, as shown in Figure 3 for the USA (Panel A) and for EEA countries (Panel B). These figures highlight the earlier beginning of the epidemic in EEA countries compared to the USA. The USA projected peak was reached on April 15 with almost 3,500 deaths daily. In EEA the peak was on April 6 with more than 4,000 deaths daily but with a flatter peak, reflecting the considerable variability in the timing of the epidemic by country. Our estimated peak hospital demand was 68,884 (95% UI 34,599–175,312) beds, 18,269 (9,621–44,223) ICU beds and 16,545 (8,083–41,991) ventilators in the USA; for EEA nations the corresponding numbers were 120,080 (119,183–121,107) hospital beds, 32,291 (32,157–32,425) ICU beds, and 28,973 (28,868–29,085) ventilators.
The peak of daily deaths varies considerably between EEA countries and subnational locations (Figure 4, Panel A) and USA states (Panel B). Several regions in Italy reached their peak by the end of March, with parts of Spain, France, Netherlands, Norway, Denmark, Greece, and Estonia following suit by the beginning of April. Other countries such as the UK, Germany, and Sweden are at the peak or are approaching the peak. In the USA, states with earlier peaks include Washington, Nevada, Arizona, Montana and Florida. States at the peak or just approaching the peak include Texas, California and parts of New England. States in the middle of the country, including North Dakota, South Dakota, Iowa and Wyoming are expected to peak later.
Figure 5 shows total excess demand for the USA (panel A) and EEA countries (panel B) overall. In the USA, peak excess demand for hospitalisation above usual capacity was estimated as 9,079 (95% UI 253–61,937); ICU bed excess demand was 9,356 (3,526–29,714). We estimated that EEA countries experienced a peak excess demand above usual capacity for total beds of 28,270 (0 to 126,788) at peak; the ICU bed shortfall was 16,090 (15,973–16,211). Excess demand is concentrated in particular countries and USA states as shown in Figure 6, which shows the percentage excess demand for ICU beds by location: in the USA (panel A), excess demand for ICU beds is concentrated in New York, New Jersey, Connecticut, Wyoming, Michigan, Rhode Island, and Massachusetts; in the EEA (panel B), ICU excess demand above usual capacity is particularly high in Sweden, Spain, Northern Ireland, Italy, France, and Belgium. We have not been able to estimate current ventilator capacity; however, the number of ventilators per person implied by the peak (Figure 3) also suggests potentially large gaps in availability of ventilators.
Figure 7 shows the expected cumulative death numbers with 95% uncertainty intervals for the USA (Panel A) and EEA (Panel B). In the USA, the average forecast suggests 60,308 deaths, but the range is large, from 34,063 to 140,381 deaths. The figure shows that uncertainty widens markedly as the peak of the epidemic approaches, given that the exact timing of the peak is uncertain. In the EEA, 91,972 (95% UI 91,212–93,620) deaths have already been recorded so far, with the majority of these coming from Italy, Spain, and France. Our forecast suggests a cumulative total of 143,088 (101,131–253,163) deaths in the EEA. A large number of these deaths are projected to occur in the UK (13,759 observed to date; 37,521 [17,625–89,385] total), Sweden (1,333 observed to date; 5,890 [1,965–16,883] total), Germany (3,570 observed to date; 4,957 [3,697–9,379] total) and France (18,485 observed to date; 22,555 [19,455–29,314] total).
Figure 8 shows a map of the cumulative number of deaths per capita by location for the USA and EEA. In Europe, the highest number of estimated cumulative deaths per capita are in Italy – particularly the northern regions – Spain, Belgium, Sweden and the UK. In the USA, states with the highest per capita deaths are New York, Rhode Island, New Jersey, Connecticut, Massachusetts, Wyoming, Louisiana, and Michigan.
Figure 9 shows the date by location by which projected daily deaths drop below 0.3 per million. As expected, there is a strong correlation between the timing of the peak daily death and when the daily death rate will drop below this threshold. In Europe, countries where this will happen later include the UK, Norway, Denmark, Sweden and the Netherlands. In the USA, states that will not cross this threshold until the end of May include South Dakota, North Dakota, Iowa, Oklahoma, Arkansas and Utah.
Results for each location are accessible through a visualisation tool at http://covid19.healthdata.org/projections – the estimates presented in this tool will be continually updated as new data are incorporated and ultimately will supersede the results in this paper. Summary information on cumulative deaths, the date of peak demand, the peak demand, peak excess demand, and aggregate demand are provided for each location in Table 1.
Discussion
This study has generated estimates of predicted health service utilisation and deaths due to COVID-19 by day through the end of July for all USA states and EEA countries, assuming that social distancing efforts will continue until deaths reach a very low level. The analysis shows large gaps between need for hospital services and usual capacity, especially for inpatient and ICU beds. A similar or perhaps even greater gap for ventilators is also likely, but detailed state or country data on ventilator capacity are not available to directly estimate that gap. Uncertainty in the time course of the epidemic, its duration, and the peak of utilisation and deaths is large, particularly for when locations are early in the epidemic and where there are few deaths. Given this, it is critical to update these projections as the pandemic progresses and new data are collected. Uncertainty will also be reduced as we gain more knowledge about the epidemic peak and subsequent decline in daily deaths across more than 13 locations. A critical aspect to the size of the peak is when aggressive measures for social distancing are implemented in each state, region, or country and for how long they are maintained. Delays in implementing government-mandated social distancing and relaxing policies will have an important effect on the resource gaps that health systems will be required to manage.
Our estimates of excess demand show that hospital systems have already or will face difficult choices to continue providing high-quality care to their patients in need. This model was first developed for use by the UW Medicine system in Washington state, and the practical experience of that system provides insight into how it has been useful for planning purposes. From the perspective of planning for the UW Medicine system, these projections immediately made apparent the need to rapidly build available capacity. Strategies to do so included suspending elective and non-urgent surgeries and procedures, while supporting surge planning efforts and reconfiguration of medical/surgical and ICU beds across the system. These targets also supported a proactive discussion regarding the potential shift from current standards of care to crisis standards of care, with the goal to do the most good for the greatest number in the setting of limited resources.
There are a variety of options available to deal with the situation, some of which have already been implemented or are being implemented. One option is to reduce non-COVID-19 patient use. In the USA and in many EEA countries, local, state, or national governments have cancelled elective procedures40–45 and many, but not all, hospitals have complied. This decision has significant financial implications for USA health systems, however, as elective procedures are a major source of revenue.46 Also, aggressive social distancing policies reduce not only the transmission of COVID-19 but will likely have the added benefit of reducing health-care utilisation due to other causes such as injuries.47 Reducing non-COVID-19 demand alone will not be sufficient, and strategies to increase capacity are also clearly needed. This includes setting up additional beds by repurposing unused operating rooms, pre- and post-recovery rooms, procedural areas, medical and nursing staff quarters, and hallways.
Currently, one of the largest constraints on effective care may be the lack of ventilators. One supplement to ventilator capacity is using anesthesia machines freed up by deferring or cancelling elective surgeries. Other options go beyond the capacity or control of specific hospitals. The use of mobile military resources has the potential to address some capacity limitations, particularly in the USA given the differently timed epidemics across states. Other innovative strategies will need to be found, including the construction of temporary hospital facilities as has been done in Wuhan,48 Washington state,49 New York,50,51 Italy,52 France,53 and Spain.54
In this study, we have quantified the potential gap in physical resources, but there is an even larger potential gap in human resources (HR). Expanding bed capacity beyond licensed bed capacity may require an even larger increase in the HR to provide care. The average annual bed-day utilisation rate in the US is 66% and ranges from 46%to 92% among EEA countries. Many hospital systems are staffed appropriately at their usual capacity utilisation rate, and expanding even up to, but then potentially well beyond, licensed capacity will require finding substantial additional HR. Strategies include increasing overtime, training operating room and community clinic staff in inpatient care or physician specialties in COVID-19 patient care, rehiring recently separated workers, and the use of volunteers. In academic health systems such as UW Medicine, clinical faculty time can be redirected from research and teaching to clinical care during the COVID-19 surge. A more concerning HR bottleneck identified, given the need for ICU care for COVID-19 patients, is for ICU nurses, for which there are very limited options for increasing capacity. In addition to HR, what should not be overlooked is the increased demand for supplies ranging from personal protective equipment (PPE), medication, and ventilator supplies to basics such as bed linen. Add to these the need to expand other infrastructure required to meet the COVID-19 surge, such as information technology (IT) for electronic medical records. The overall financial cost over a short period of time is likely to be enormous, particularly when juxtaposed against the substantial reductions in revenue for many hospitals due to the cancellation of elective procedures and the broader economic consequences of social distancing mandates.
Our model suggests that the timing of the implementation of social distancing mandates is a critical determinant of peak demand and cumulative deaths. Mobility data derived from cell phone use has provided the basis for evaluating the importance of the different social distancing mandates included in the social distancing covariate. It is important to note the social distancing mandates do not capture all variation in mobility, and that the data in some locations suggests behavioural change prior to the introduction of these mandates. Understanding what drives individual change, e.g. levels of awareness or fear of the pandemic or the private sector implementing remote work policies prior to government mandates, will be important for understanding what may drive the change in behaviour after official social distancing mandates are relaxed, which is now beginning in some European countries and US States.
Based on our experience thus far, we have derived important insights into the epidemic trajectories and health service demand as data have accumulated. These have led to improved forecasts reflecting both new data and method refinements.55,56 For this reason, we are continuing to revise the model as new data are available, providing an updated forecast for health service providers, governments, and the public. In some regions that have peaked, such as regions in Italy like Liguria, or New York, the duration of the peak is much longer than in other places such as Madrid. The mixture model we use accommodates this longer peak but it remains unclear why some communities have the prolonged peak and others do not. The prolonged peak leads to substantially increased total mortality. There is also marked variation across locations in how steeply the epidemic curve rises, captured by the alpha parameter in our model. Understanding why some locations have an epidemic like New York and others like Washington State will be important to make robust forecasts in other regions of the world.
Any attempt to forecast the COVID-19 epidemic has many limitations. Only a limited number of locations with generalised epidemics have reached the peak in terms of daily deaths, and only one has currently brought new cases to near zero, namely Wuhan. Many other locations, including all other provinces in China, have so far successfully contained transmission, preventing a general outbreak. Modelling based on one completed epidemic, at least for the first wave, and many incomplete epidemics is intrinsically challenging. The main limitation of our study is that observed epidemic curves for COVID-19 deaths define the likely trajectory. In this study, we do include a covariate meant to capture the timing of social distancing measures and their effect on various measures of population mobility. Our model also relies on the accuracy of reporting of deaths due to COVID-19; reports suggest that in some locations not all deaths may be included in country reported totals.57,58 Our models explicitly take into account variation in age-structure, which is a key driver of all-age mortality. But these efforts at quantification do not take into account many other factors that may influence the epidemic trajectory: sex, the prevalence of co-morbidity, population density, individual behavior change not captured by mobility metrics, and a host of other individual factors that may potentially influence the immune response. We also have not explicitly incorporated the effect of reduced quality of care due to stressed and overloaded health systems beyond what is captured in the data. For example, the higher mortality rate in Italy may in part be due to policies around restricting invasive ventilation in the elderly. The model ensemble used does suggest that locations with faster increases in the death rate are likely to have greater peak caseload and cumulative deaths, but our uncertainty intervals are appropriately large. Finally, it is critical to note that we restrict our projections to the first wave of the pandemic under a scenario of continued implementation of social distancing mandates and do not yet incorporate the possibility of a resurgence or subsequent waves. This is an essential area for future work.
Conclusion
COVID-19 is an extraordinary challenge to health and the health-care system. In this study, we forecast a large excess of demand for hospital bed-days and ICU bed-days and our estimate of 1,584,737 (95% UI 1,050,954–3,082,999) deaths in the USA and EEA from the first wave of pandemic is an alarming number. This number could be substantially higher if excess demand for health system resources is not addressed and if social distancing policies are not continued, vigorously implemented, and enforced. This planning model will hopefully provide an up-to-date tool for improved hospital resource allocation.
Data Availability
Additional information on the determination of hospital resource utilisation and capacity is provided in Appendix A; details on curve fitting methods, quantification of uncertainty, and a full specification of the statistical model are available in Appendix B. This study complies with the Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) statement.