Skip to main content
  • Research article
  • Open access
  • Published:

The association between mechanical ventilator compatible bed occupancy and mortality risk in intensive care patients with COVID-19: a national retrospective cohort study

Abstract

Background

The literature paints a complex picture of the association between mortality risk and ICU strain.

In this study, we sought to determine if there is an association between mortality risk in intensive care units (ICU) and occupancy of beds compatible with mechanical ventilation, as a proxy for strain.

Methods

A national retrospective observational cohort study of 89 English hospital trusts (i.e. groups of hospitals functioning as single operational units). Seven thousand one hundred thirty-three adults admitted to an ICU in England between 2 April and 1 December, 2020 (inclusive), with presumed or confirmed COVID-19, for whom data was submitted to the national surveillance programme and met study inclusion criteria. A Bayesian hierarchical approach was used to model the association between hospital trust level (mechanical ventilation compatible), bed occupancy, and in-hospital all-cause mortality. Results were adjusted for unit characteristics (pre-pandemic size), individual patient-level demographic characteristics (age, sex, ethnicity, deprivation index, time-to-ICU admission), and recorded chronic comorbidities (obesity, diabetes, respiratory disease, liver disease, heart disease, hypertension, immunosuppression, neurological disease, renal disease).

Results

One hundred thirty-five thousand six hundred patient days were observed, with a mortality rate of 19.4 per 1000 patient days. Adjusting for patient-level factors, mortality was higher for admissions during periods of high occupancy (> 85% occupancy versus the baseline of 45 to 85%) [OR 1.23 (95% posterior credible interval (PCI): 1.08 to 1.39)]. In contrast, mortality was decreased for admissions during periods of low occupancy (< 45% relative to the baseline) [OR 0.83 (95% PCI 0.75 to 0.94)].

Conclusion

Increasing occupancy of beds compatible with mechanical ventilation, a proxy for operational strain, is associated with a higher mortality risk for individuals admitted to ICU. Further research is required to establish if this is a causal relationship or whether it reflects strain on other operational factors such as staff. If causal, the result highlights the importance of strategies to keep ICU occupancy low to mitigate the impact of this type of resource saturation.

Peer Review reports

Background

From the first reports of a novel coronavirus (SARS-CoV-2) in late 2019, to date, global mortality associated with the resultant disease (COVID-19) has exceeded 1.7 million people [1]. The virulence of the pathogen has prompted persistent concern about the ability of health services around the world to effectively care for the vast numbers of people affected [2]. These concerns are most relevant in the context of scarce resources (e.g. mechanical ventilation) required by patients in need of high-acuity support, which is relatively common in patients with COVID-19 [3]. Notably, even with the introduction of non-pharmacological interventions such as stay-at-home orders, almost a third of all hospitals in England reached 100% occupancy of their ‘surge’ mechanical ventilation capacity (i.e. including the additional beds that were created through the re-allocation of resources) during the first wave of the pandemic [4]. The second wave in England appears to have been worse than the first, with 40% more people in hospital, many hospitals overwhelmed and staff wellbeing suffering due to the prolonged onslaught [5]. What remains unclear is whether and to what extent operating at these extremes of critical care occupancy impacted patient outcomes.

Availability of certain resources, such as bed occupancy/availability, is well-established proxies for operational strain [6]. Prior to COVID-19, the literature painted a complex picture of the association between morbidity and mortality and ICU strain [7]. Some studies suggest that higher strain in intensive care was associated with increased morbidity and mortality risk [8,9,10], whereas others observed no association at all [11, 12]. Data from the natural experiment caused by the COVID-19 pandemic is accumulating; however, it is similarly inconsistent. A recent study from Belgium reported 42% higher mortality during periods of ICU surge capacity deployment, although in the analysis utilisation of surge capacity was only evaluated as a binary variable [13]. Similarly, a study based on the US Veterans Affairs database evaluating quartile-based bins for occupancy found a 94% increased risk of mortality when occupancy exceeded 75% [14]. The aforementioned observations both directly contradict earlier results from smaller studies in Australia and Wales, where no association between ICU occupancy and mortality was identified [15, 16].

A better understanding of how operating under such extreme circumstances affects outcomes is crucial for two reasons: firstly, to allow hospitals to adapt practice to improve outcomes, and secondly, to provide policy makers with more accurate information about the potential consequences of allowing health services to be overwhelmed. As such, in this study, we sought to evaluate the extent to which mortality risk in intensive care units (ICUs) during the COVID-19 pandemic in England could be explained by differences in occupancy.

Methods

Data sources

Data on all intensive care unit (ICU) admissions across England were extracted from the COVID-19 Hospitalisation in England Surveillance System (CHESS) [17]. Information on occupancy rates were extracted from the daily situation reports (i.e. ‘SitRep’) [4]. Both datasets are mandatory regulatory submissions for all National Health Service (NHS) acute care providers in England, and further details about them can be found in additional file 1: eMethods.

Study population

All ICU admissions reported to CHESS between 2 April and 1 December (see additional file 1: eMethods for details on date selection), with presumed/confirmed COVID-19 (100% tested positive during their admission; see additional file 1: eMethods for details on diagnosis), aged 18–99, non-pregnant, and with valid admission and occupancy data, were eligible for inclusion (additional file 1: Fig. S1).

Recorded clinical features

Patient-level data

Information extracted from CHESS about each patient comprised: administrative features (admitting trust, admission date, and first segment of postcode—used to identify the relevant indices of multiple deprivation for the corresponding areas in England)), demographic characteristics (age, sex, ethnicity), recorded comorbidities (obesity, diabetes, asthma, other chronic respiratory disease, chronic heart disease, hypertension, immunosuppression due to disease or treatment, chronic neurological disease, chronic renal disease, chronic liver disease), and a binary indicator for use of mechanical ventilation. Ethnicity was coded as white, Asian (Subcontinent and other), black, mixed, and other; comorbidities were recorded as binary indicator variables, with missing entries assumed to reflect the absence of a comorbidity. The appropriateness of this assumption in CHESS has been previously explored [18].

Occupancy data

Trusts are groups of geographically co-located hospitals that function as a single organisational unit within the UK’s national healthcare system. Information extracted from SitRep about each trust comprised: pre-pandemic (January–March 2020) number of beds compatible with mechanical ventilation, the proportion of beds compatible with mechanical ventilation occupied on each day of the study period, and each trust’s geographical region [4]. Linkage was carried out based on the trust that an individual was admitted to and the date of ICU admission in CHESS; patients in CHESS were matched via their admission date to the relevant occupancy information from the corresponding date in SitRep.

Outcome

The primary outcome was in-hospital all-cause mortality; patients were followed up until death, discharge, transfer, or a minimum of 90 days from admission. Discharge and transfer were both treated as suggesting that the patient survived.

Statistical analysis

Descriptive summaries were generated as median and interquartile ranges for continuous variables and frequency and percentage incidence for categorical factors. Exploratory analyses included: the relationship of the COVID-19 epidemic curve to bed occupancy at a national level (additional file 1: Fig. S2); the distribution of missingness amongst patients and trusts (additional file 1: Fig. S3); variation in age and comorbidity burden over the first wave (additional file 1: Fig. S4); the impact of modelling continuous variables either linearly, through the use of threshold functions, or via (standard cubic) splines/smooths (additional file 1: eMethods).

A Bayesian hierarchical approach was used to model the association between the trust, group and patient-level factors and mortality risk. Specifically, a generalised additive mixed model was utilised, with intercept and slope coefficients for population and group level effects, and a Bernoulli likelihood with logit function to link to mortality outcome. Coefficients were inferred by Markov chain Monte Carlo sampling, using Stan (CmdStan V2.25.0), with a model specified using BRMS (V2.14.4) in R (V4.0.3) [19,20,21]. Further information on the Bayesian prior specification and modelling methods are reported in additional file 1: eMethods. The rationale for the different thresholds utilised for bed occupancy is also detailed in additional file 1: eMethods; in short, the 85% occupancy threshold utilised is based on guidance from the Royal College of Emergency Medicine [22], whereas the 45% (low occupancy threshold) was derived from the exploratory analysis.

As secondary analyses, two potential interactions were assessed: (1) baseline trust size and occupancy and (2) patient age and occupancy (results not reported due to insignificance). We also assessed the association of occupancy on the recorded outcome date with mortality and occupancy expressed in terms of pre-pandemic ICU size. Several sensitivity analyses were carried out (see additional file 1: eMethods for justifications): (1) filtering for different degrees of missingness of patient-level comorbidity data at trust-level; (2) including all individuals still on the unit as of 2 March (and assuming they all survive or die; n = 250); (3) including all individuals with a known outcome but no date whom were excluded in the data cleaning process (n = 10); (4) a random effect for all patient post 16 June 2020 (i.e. the date of the Recovery Trials press-release regarding the effectiveness of dexamethasone); (5) adjusting for week of admission; (6) adjusting for trust and region as random effects; (7) additional patient-level factors: time from hospital admission to ICU admission, chronic liver disease and obesity; (8) the residual bed availability in the sustainability and transformation partnership (STPs) to which each trust belongs—STPs are collections of geographically co-located hospital trusts that are meant to operate collaboratively to support one another, and include ambulance providers and community services so as to allow end-to-end integration of care; (9) the weighted average of the indices of multiple deprivation (IMD) corresponding to each patient’s given postcode district (i.e. the first half of their postcode), as a proxy for socio-economic status; and (10) separated into the different waves of the pandemic in the UK (i.e. 1st, inter-wave/summer, and 2nd).

Results

Seven thousand one hundred thirty-three individuals were included in this study following application of the inclusion/exclusion criteria (see additional file 1: Fig. S1), of whom 2637 (37.0%) died. In total, 135,600 (median 12 days per patient; IQR 6–25) patient days were observed, equating to a mortality rate of 19.4 per 1000 patient days. A full summary of the recorded patient-level characteristics is reported in Table 1. Occupancy of beds compatible with mechanical ventilation over the study observation period is illustrated in additional file 1: Fig. S5, expanding on the results presented elsewhere which are restricted to the first wave [4].

Table 1 Characteristics of the study cohort stratified by occupancy on the day of admission

Mechanical ventilator occupancy rate on the day of admission is associated with mortality risk

For high occupancy rates (85–100% of total beds, including surge capacity), the unadjusted odds ratio (OR) for mortality based on the mechanical ventilator occupancy rate on the day of admission was 1.22 (95% posterior credibility interval (PCI) 1.06–1.40), relative to the baseline of 45–85%. For low occupancy rates (0–45%), the unadjusted odds ratio was 0.82 (95% PCI 0.73–0.93), relative to the baseline of 45–85%. Following adjustment for patient-level factors (age, sex, ethnicity, and comorbidities), the OR was 1.23 (95% PCI 1.05–1.43) for high occupancy rates and 0.83 (95% PCI 0.73–0.96) for low occupancy rates. Figure 1 illustrates the posterior probabilities for the fully adjusted ORs (see additional file 1: Fig. S6 & S7 for other model coefficients). To aid interpretation, the difference in risk for a 70-year-old man with no comorbidities being admitted during a period of high versus low occupancy is equivalent to the risk of them being over a decade older (see full relationship in Fig. 2). Sensitivity analyses as detailed in additional file 1: Table S1 were all concordant with the primary analysis, and the results do not appear to be sensitive to the choice of prior as illustrated in additional file 1: Table S2.

Fig. 1
figure 1

The adjusted odds ratios for the risk of mortality based on different ICU bed occupancy rates (treated as a three-level categorical variable). The full posterior distribution of the odds ratio (OR) for mortality given low occupancy 0–45% (top; green), and high occupancy 85–100% (bottom; red) are presented. The PCIs presented are equally tailed credibility intervals for the posterior odds ratio distributions. The outer (less saturated) interval is the 95% PCI, and the inner (more saturated) interval shows the 90% PCI. The reference category is 45–85% occupancy. Exact values are tabulated

Fig. 2
figure 2

The adjusted odds ratios for the risk of mortality based on ICU bed occupancy (treated as a linear continuous variable) on the day of admission (top) and each individual’s recorded outcome date (bottom). The full posterior distribution of the odds ratio (OR) for mortality given occupancy on the date of ICU admission (top; purple), mean occupancy during ICU stay (middle: pink), and occupancy on the date of each individual’s recorded outcome (bottom; blue) are presented. The PCIs presented are equally tailed credibility intervals for the posterior odds ratio distributions. Occupancy was specified without multiplying out by 100 (i.e. 20% = 0.20); therefore, the odds ratio is for a step from 0% to 100% (i.e. 0.0 to 1.0). Exact values are tabulated

Mortality risk increases linearly with admission date specific occupancy and average exposure over the full length of stay

The fully adjusted OR for mortality (Fig. 3), using occupancy on the day of admission coded as a continuous linear variable ranging from 0 to 1 (i.e. specified without multiplying out by 100 (e.g. 20% = 0.20), therefore the resulting odds ratio is for a step from 0% to 100%), was 1.59 (95% PCI 1.22–2.03). Furthermore, using the occupancy from each individuals’ outcome date identified an even larger association (full model specification in additional file 1: Table S3), OR 2.29 (95% PCI 1.80–2.92). Finally, the fully adjusted OR for the mean occupancy over a patient’s total length of stay in ICU was 2.69 (95% PCI 2.05–3.50); see additional file 1: Table S4.

Fig. 3
figure 3

The increase in mortality risk associated with admission to intensive care during periods of different occupancy rates, expressed in terms of the equivalent increase in years of age. (Left) The predicted mortality curves arising from predictions made by the primary model across a range of age values for a white male patient is shown alongside 95% credible intervals in a ribbon either side of the median line. The black dotted line intersects all three curves; the 0–45% and 85–100% occupancy curve y value probabilities can then be used to solve back onto the reference curve to determine effective ages of equal risk to the chosen age under reference 45–85% occupancy, shown by the corresponding green and red dotted lines respectively [23]. (Right) The plot illustrates the number of years of additional age that ICU admission on a day with each different mechanical ventilation bed occupancy rate equates to. For example, an individual with a chronological age of 40 has an effective age of 31 in a low occupancy setting (green) and 45 in a high occupancy setting (orange). Both of the aforementioned comparisons are relative to the baseline occupancy of 45–85%). A comparison of the difference in risk between being admitted to the highest occupancy range relative to the lowest is shown in (red) and for a 40-year-old is equivalent to an increase in their age by 11 years

Mortality risk is evident regardless of comparison to surge or baseline capacity as the reference

Using the alternative definition of occupancy, relative to baseline (i.e. pre-pandemic) capacity instead of the real-time (surge) capacity on the day of admission illustrates that mortality risk given a linear continuous factor was 1.31 (95% PCI 1.16–1.48; additional file 1: Table S5). In other words, every 100% increase in occupancy above baseline capacity was associated with a 31% increase in risk of mortality.

Larger ICUs experience exaggerated impacts of extremes of mechanical ventilator occupancy rates

Pre-pandemic number of beds did not substantially alter the OR of occupancy as a sensitivity analysis (additional file 1: Table S1). Moreover, for every additional 25 beds compatible with mechanical ventilation (pre-pandemic), there was no observed (significant at the 5% threshold) decrease in mortality (additional file 1: Table S6). However, introduction of an interaction term between pre-pandemic size and real-time occupancy identified that larger ICUs experience exaggerated impacts of extremes of mechanical ventilator occupancy rates (additional file 1: Fig. S8).

Discussion

The results of this study, including the individual-level experience of over 7000 patients admitted to intensive care units across England, suggest that survival rates for patients with COVID-19 in intensive care settings appears to deteriorate as the occupancy of (surge capacity) beds compatible with mechanical ventilation (a proxy for operational pressure [6]), increases. Moreover, this risk does not occur above a specific threshold, but rather appears linear, whereby going from 0% occupancy to 100% occupancy increases the odds of mortality by 59% (after adjusting for relevant individual-level factors). Furthermore, and expectedly, greater exposure to periods of high occupancy over an individual’s total length of stay is also associated with an increased risk of mortality. This lends credibility (via a potential dose-response relationship) to the argument that this is a causal relationship that is being observed; however, we cannot make a definite claim based on our results. An additional piece of context for interpreting these results is that the relevant authorities in the UK made clear that no explicit triaging unique to the pandemic took place [24], meaning that everyone whom would have previously been deemed suitable for ICU treatment was treated similarly to prior to the pandemic.

The remaining result for which the interpretation is not immediately apparent is that the risk of mortality based on occupancy on the date of recorded outcome is also high; OR 2.29 (95% posterior credible interval 1.80–2.92). There are several plausible explanations, but for the sake of brevity, we will mention just two. The first is based on the same principle as the other results, i.e. as occupancy increases the availability of staff and resources decreases which can indirect contribute to an individual’s mortality risk, and/or there may be a direct effect due to rationing of care or treatment limiting orders. Alternatively, this result could also reflect a specific type of observation bias wherein patients who are likely to survive are kept in ICU and in hospital longer [25], thus resulting in them being eventually discharged during periods of lower occupancy. In both instances further research is needed to clarify the underlying causal mechanism which is driving these observations.

In context of the literature

Our observations are consistent with the aforementioned Belgian and US VA study [13, 14], except that they extend the results to suggest a linear association rather than stepped increase at specific thresholds. Although previous COVID-19 specific studies report an association between larger ICUs and lower COVID-19 mortality [26], our results are equivocal. However, we did observe an interaction with (pre-pandemic) unit size, whereby larger ICUs experience more exaggerated impacts from both higher and lower (surge) occupancy rates, although it is unclear from our data what is driving this heterogeneity. Finally, these results might partially explain the decreased mortality rate seen in the latter half of the first wave in the UK [27], where occupancy rates were much lower than at the peak [4]. Importantly, the average number of co-morbidities per patient appears to rise when occupancy decreases (additional file 1: Fig. S4), suggesting if there was implicit/soft triaging it would not have been in a manner that would explain the observed effect (i.e. more co-morbidities are generally thought to be associated with a worse outcome). However, more research is necessary to definitively exclude this potential explanation for the results.

Strengths and limitations

The strengths of this study are the national cohort of patient-level data with extensive capture of admissions [25], coupled with a rigorous modelling method (additional file 1: eMethods). However, there are also several limitations. For example, we lack 30-day outcome data for discharged and transferred individuals and thus were forced to model under a naïve assumption that these individuals survived, which may have impacted our estimates of the risk of mortality. Below is a more detailed discussion on three key issues that should be accounted for when interpreting the associations reported in this study:

Availability of physiological and clinical data

The lack of physiological data and minimal clinical intervention data available in the CHESS dataset limits our ability to adjust for individual-level differences in severity and treatment. Thus, we cannot exclude that the observed associations are a result of unmeasured confounding due to differences in severity on admission or variation in natural history due to the evolving treatment strategies that have been used over the course of the pandemic. At present, there is no single research environment that links primary care, secondary care, administrative mortality records, surveillance datasets (e.g. CHESS), routine COVID testing (e.g. Public Health England’s testing data), and national audit data such as that collected by the UK Intensive Care National Audit and Research Centre (ICNARC) which includes some of the necessary information at the patient-level. Notably, there are on-going attempts to link the aforementioned datasets, e.g. [28]. However, this exercise is still not complete, and at the time of writing this, the (linked) ICNARC dataset is not available. In essence, a comprehensive analysis including physiological and treatment data alongside all the other ancillary datasets is not possible. To our knowledge, the UK consortium is the first attempt to make such a national dataset and thus replication of our analysis with these data once available is a critical future task.

Although direct analysis to address the above issues is not possible for the reasons outlined previously, it is still worth noting that several studies from the UK suggest that any potential differences in severity on admission are unlikely to explain the observed association between occupancy and mortality. For example, a study using linked CHESS data from the first wave of the pandemic did not find between centre variation in severity scores (e.g. mean APACHE-II) to be associated with mortality risk [29]. Moreover, the UK national audit (ICNARC) did not identify any significant temporal variation in mean severity scores of patients admitted to ICU over the course of this study [30]; both results in combination suggest that unmeasured confounding due to differences in severity of admission are unlikely to fully explain the observed association.

Incomplete characterisation of operational strain

Characterisation of operational strain as a function of surge occupancy likely fails to fully reflect the complexity of using non-specialist staff and other resource allocation issues present when ‘creating’ new ICU beds, which could be the underlying causal explanation for the increased mortality risk observed during periods of high occupancy. Notably, previous research has identified a detrimental impact of insufficient (adequately trained) nursing and medical consultant staff on patient mortality risk in intensive care units [31], and staff absence rates were raised 3-fold from the baseline of 4% at the peak of the first wave of the pandemic [32]. As such, it is not just functional strain due to the creation of new beds, but also absence that could underlie a staff-specific explanation for the observed association. Again, this is a key area for future research to explore, especially with mounting evidence of operational strain being associated with increased mortality risk [33].

The prevalence of variants of concern

Despite the more recent identification of several effective treatments for severe COVID-19, as bed saturation in ICUs again became an issue from November 2020 to January 2021, mortality rates increased across England to rates similar to that seen during the 1st wave [34]. Notably, B117/Alpha was the most prevalent strain of COVID-19 in the UK during November 2020 [35] and is known to have increased transmission and mortality risk [36]. Although we attempted to limit the potential impact of varying prevalence of different variants of concern (VOC), using specific date windows for our analysis and assessing week of admission as a sensitivity analysis, without knowledge of the specific variant COVID-19 for each individual diagnosis, it is impossible to exclude this as a cause of residual confounding.

Implications for policy makers and clinicians

The results presented in this study are especially relevant given the identification of variants of concern with increased risk of transmission and mortality [36], as well as vaccine/immunity escape potential [37], which suggest that the risk of further epidemics driven by Sars-CoV-2 are not implausible. From a mitigation perspective, our study highlights the importance of public health interventions (such as expeditious vaccination programmes and non-pharmacological interventions), to control both incidence and prevalence of COVID-19, and therefore prevent ICU saturation, as there is evidence of direct harm to patients as a consequence. In terms of adaptation, during the second wave over 1000 patients in critical care units across the UK were transferred due to capacity reasons. Whilst there are pre-existing mechanisms to support these action [38], there is need for more investment in digital infrastructure to support better coordinated action at the national level. Again, the natural experiment produced by the pandemic has produced substantial amounts of data on transferred patients and load-sharing between trusts. The community would benefit greatly from analysis of this real world data to inform future best practices, from whom to consider transferring, how far, as well as a host of other factors that influence decision making, rather than relying solely on model-based inferences which suggest load-sharing is possible and probably effective [39]. Moreover, there are number of other adaptation interventions that might be worth considering for future research, such as modified treatment strategies with lower sedation doses and more aggressive extubation strategies to increase discharge rates during periods of high occupancy and therefore an increased baseline risk of mortality independent of individual-level factors.

Conclusion

The results of this study suggest that increasing occupancy of beds compatible with mechanical ventilation, a proxy for operational strain, is associated with a higher mortality risk for individuals admitted to ICU. Further research is required to establish if this is a causal relationship or whether it reflects strain on other operational factors such as staff/residual confounding not addressed due to the lack of physiological severity data. If causal, the result highlights the importance of strategies to keep ICU occupancy low to mitigate the impact of this type of resource saturation, rather than waiting for arbitrary thresholds to be triggered, especially given the observation that average exposure (over an individual’s total length of stay) to this operational risk factor is also associated with mortality risk.

Availability of data and materials

Data, even a de-identified version, cannot be shared publicly as it was collected by Public Health England (PHE) as part of their statutory responsibilities, which allows them to process patient confidential data without explicit patient consent. Given this particular set of circumstances, it was deemed by PHE to not be appropriate to share the data publically as then they would then be unable to fulfil their duty of care to ensure that the data is used in a manner consistent with the statutory remit under which it was collected. Moreover, to retain the fidelity of the data and given that there are numerous unique cases due to low prevalence in certain geographic areas over the time period we analyse, de-identification would be insufficient to comply with the minimum Office for National Statistics’ disclosure principles which would come with significant risk of re-identification. Data utilised in this study were made available through an agreement between the University of Warwick and PHE. Individual requests for access to the raw data are considered directly by PHE (contact via covid19surv@phe.gov.uk).

Abbreviations

ICU:

Intensive care units

CHESS:

COVID-19 Hospitalisation in England Surveillance System

SitRep:

Daily situation reports

NHS:

National Health Service

STPs:

Sustainability and transformation partnership

IMD:

Indices of multiple deprivation

OR:

Odds ratio

PCI:

Posterior credibility interval

IQR:

Interquartile range

VOC:

Variant of concern

ICNARC:

Intensive Care National Audit and Research Centre

References

  1. World Health Organisation. COVID-19 Weekly Epidemiological Update - 29 December 2020; 2020. https://www.who.int/publications/m/item/weekly-epidemiological-update%2D%2D-29-december-2020. Accessed 1 Jan 2021.

  2. Ferguson N, Laydon D, Nedjati Gilani G, et al. Report 9: Impact of non-pharmaceutical interventions (NPIs) to reduce COVID19 mortality and healthcare demand. https://www.imperial.ac.uk/media/imperial-college/medicine/mrc-gida/2020-03-16-COVID19-Report-9.pdf . Accessed 1 Jan 2021.

  3. Docherty AB, Harrison EM, Green CA, Hardwick HE, Pius R, Norman L, et al. Features of 20 133 UK patients in hospital with covid-19 using the ISARIC WHO Clinical Characterisation Protocol: prospective observational cohort study. BMJ. 2020;22:369.

    Google Scholar 

  4. Mateen BA, Wilde H, Dennis JM, Duncan A, Thomas N, McGovern A, et al. Hospital bed capacity and usage across secondary healthcare providers in England during the first wave of the COVID-19 pandemic: a descriptive analysis. BMJ Open. 2021;11(1):e042945. https://doi.org/10.1136/bmjopen-2020-042945.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Mahase E. COVID-19: hospitals in crisis as ambulances queue and staff are asked to cancel leave. BMJ. 2020;371:m4980. https://doi.org/10.1136/bmj.m4980.

    Article  PubMed  Google Scholar 

  6. Bagshaw SM, Wang X, Zygun DA, Zuege D, Dodek P, Garland A, et al. Association between strained capacity and mortality among patients admitted to intensive care: a path-analysis modeling strategy. J Crit Care. 2018;43:81–7. https://doi.org/10.1016/j.jcrc.2017.08.032.

    Article  PubMed  Google Scholar 

  7. Eriksson CO, Stoner RC, Eden KB, Newgard CD, Guise JM. The association between hospital capacity strain and inpatient outcomes in highly developed countries: a systematic review. J Gen Intern Med. 2017;32(6):686–96. https://doi.org/10.1007/s11606-016-3936-3.

    Article  PubMed  Google Scholar 

  8. Chrusch CA, Olafson KP, McMillan PM, Roberts DE, Gray PR. High occupancy increases the risk of early death or readmission after transfer from intensive care. Crit Care Med. 2009;37(10):2753–8. https://doi.org/10.1097/CCM.0b013e3181a57b0c.

    Article  PubMed  Google Scholar 

  9. Howie AJ, Ridley SA. Bed occupancy and incidence of Methicillin-resistant Staphylococcus aureus infection in an intensive care unit. Anaesthesia. 2008;63(10):1070–3. https://doi.org/10.1111/j.1365-2044.2008.05575.x.

    Article  CAS  PubMed  Google Scholar 

  10. Gattinoni L, Radrizzani D, Simini B, Bertolini G, Ferla L, Mistraletti G, et al. Volume of activity and occupancy rate in intensive care units. Association with mortality. Intensive Care Med. 2004;30(2):290–7.

    Article  Google Scholar 

  11. Wagner J, Gabler NB, Ratcliffe SJ, Brown SE, Strom BL, Halpern SD. Outcomes among patients discharged from busy intensive care units. Ann Intern Med. 2013;159(7):447–55. https://doi.org/10.7326/0003-4819-159-7-201310010-00004.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Iwashyna TJ, Kramer AA, Kahn JM. Intensive care unit occupancy and patient outcomes. Crit Care Med. 2009;37(5):1545–57. https://doi.org/10.1097/CCM.0b013e31819fe8f8.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Taccone FS, Van Goethem N, De Pauw R, Wittebole X, Blot K, Van Oyen H, et al. The role of organizational characteristics on the outcome of COVID-19 patients admitted to the ICU in Belgium. Lancet Reg Health-Europe. 2021;2:100019. https://doi.org/10.1016/j.lanepe.2020.100019.

    Article  Google Scholar 

  14. Bravata DM, Perkins AJ, Myers LJ, Arling G, Zhang Y, Zillich AJ, et al. Association of intensive care unit patient load and demand with mortality rates in US Department of Veterans Affairs hospitals during the COVID-19 pandemic. JAMA Netw Open. 2021;4(1):e2034266.

    Article  Google Scholar 

  15. Burrell AJ, Pellegrini B, Salimi F, Begum H, Broadley T, Campbell LT, et al. Outcomes of COVID-19 patients admitted to Australian intensive care units during the early phase of the pandemic. Med J Aust. 2020;16:1.

    Google Scholar 

  16. Baumer T, Phillips E, Dhadda A, Szakmany T. Epidemiology of the first wave of COVID-19 ICU admissions in South Wales—the interplay between ethnicity and deprivation. Front Med. 2020;7. https://doi.org/10.3389/fmed.2020.569714.

  17. Public_Health_England. COVID-19 Hospitalisation in England Surveillance System (CHESS) 2020. https://www.england.nhs.uk/coronavirus/wp-content/uploads/sites/52/2020/03/phe-letter-to-trusts-re-daily-covid-19-hospital-surveillance-11-march-2020.pdf. Accessed 1 Jan 2021.

  18. Dennis JM, Mateen BA, Sonabend R, Thomas NJ, Patel KA, Hattersley AT, et al. Type 2 diabetes and COVID-19–related mortality in the critical care setting: a national cohort study in England, March–July 2020. Diabetes Care. 2021;44(1):50–7. https://doi.org/10.2337/dc20-1444.

    Article  CAS  PubMed  Google Scholar 

  19. Carpenter B, Gelman A, Hoffman MD, Lee D, Goodrich B, Betancourt M, et al. Stan: a probabilistic programming language. J Stat Softw. 2017;76(1).

  20. Bürkner PC. brms: an R package for Bayesian multilevel models using Stan. J Stat Softw. 2017;80(1):1–28.

    Article  Google Scholar 

  21. Development Core Team RR. R: a language and environment for statistical computing. 2011.

    Google Scholar 

  22. The Royal College of Emergency Medicine. Fixing emergency department staffing, systems & support to deliver excellent patient centered care. https://www.rcem.ac.uk/docs/Policy/RCEM%20Vision%202020.pdf. Accessed 3 Jan 2021.

  23. Spiegelhalter D. How old are you, really? Communicating chronic risk through ‘effective age’ of your body and organs. BMC Med Inform Decis Mak. 2016;16(1):104. https://doi.org/10.1186/s12911-016-0342-z.

    Article  PubMed  PubMed Central  Google Scholar 

  24. England N. NHS England » NHS and other professional bodies’ response to the Sunday Times, 25 October 2020 [Internet]. England.nhs.uk. 2020 [cited 14 July 2021]. Available from: https://www.england.nhs.uk/2020/10/nhs-and-other-professional-bodies-response-to-sunday-times/.

  25. Doidge JC, Gould DW, Ferrando-Vivas P, Mouncey PR, Thomas K, Shankar-Hari M, et al. Trends in intensive care for patients with COVID-19 in England, Wales, and Northern Ireland. Am J Respir Crit Care Med. 2021;203(5):565–74. https://doi.org/10.1164/rccm.202008-3212OC.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Gupta S, Hayek SS, Wang W, Chan L, Mathews KS, Melamed ML, et al. Factors associated with death in critically ill patients with coronavirus disease 2019 in the US. JAMA Intern Med. 2020;180(11):1436–46. https://doi.org/10.1001/jamainternmed.2020.3596.

    Article  CAS  PubMed  Google Scholar 

  27. Dennis JM, McGovern AP, Vollmer SJ, Mateen BA. Improving survival of critical care patients with coronavirus disease 2019 in England: a national cohort study, March to June 2020. Crit Care Med. 2020.

  28. Wood A, Denholm R, Hollings S, Cooper J, Ip S, Walker V, et al. Linked electronic health records for research on a nationwide cohort of more than 54 million people in England: data resource. BMJ. 2021;373.

  29. Qian Z, Alaa AM, van der Schaar M, Ercole A. Between-centre differences for COVID-19 ICU mortality from early data in England. Intensive Care Med. 2020;46(9):1779–80. https://doi.org/10.1007/s00134-020-06150-y.

    Article  CAS  PubMed  Google Scholar 

  30. Intensive Care National Audit and Research Centre. Report on COVID-19 in critical care: England, Wales and Northern Ireland (31 December 2020). 2020.

    Google Scholar 

  31. West E, Barron DN, Harrison D, Rafferty AM, Rowan K, Sanderson C. Nurse staffing, medical staffing and mortality in intensive care: an observational study. Int J Nurs Stud. 2014;51(5):781–94. https://doi.org/10.1016/j.ijnurstu.2014.02.007.

    Article  PubMed  Google Scholar 

  32. NHS hospital staff illness or absence up to three times usual level. Guardian. 2021 [cited 7 March 2021]. Available from: https://www.theguardian.com/uk-news/2020/dec/30/nhs-hospital-staff-illness-or-absence-up-to-three-times-usual-level.

  33. Wilde H, Dennis JM, McGovern AP, Vollmer SJ, Mateen BA. A national retrospective study of the association between serious operational problems and COVID-19 specific intensive care mortality risk. PLoS One. 2021;16(7):e0255377. https://doi.org/10.1371/journal.pone.0255377.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Dennis JM, McGovern AP, Thomas NJ, Wilde H, Vollmer SJ, Mateen BA. Trends in 28-day mortality of critical care patients with coronavirus disease 2019 in United Kingdom: a national cohort study, March 2020 to January 2021. Crit Care Med. 2021.

  35. Steel K, Donnarumma H. COVID-19 infection survey: UK Office for National Statistics; 2021.

  36. Challen R, Brooks-Pollock E, Read JM, Dyson L, Tsaneva-Atanasova K, Danon L. Risk of mortality in patients infected with SARS-CoV-2 variant of concern 202012/1: matched cohort study. BMJ. 2021;372.

  37. Wibmer CK, Ayres F, Hermanus T, Madzivhandila M, Kgagudi P, Oosthuysen B, et al. SARS-CoV-2 501Y. V2 escapes neutralization by South African COVID-19 donor plasma. Nat Med. 2021;27(4):622–5. https://doi.org/10.1038/s41591-021-01285-x.

    Article  CAS  PubMed  Google Scholar 

  38. Campbell D, Marsh S. NHS moving dozens of Covid patients a day in search for ICU beds. the Guardian. 2021 [cited 17 April 2021]. Available from: https://www.theguardian.com/world/2021/feb/03/dozens-covid-patients-move-daily-search-nhs-icu-beds.

  39. Lacasa L, Challen R, Brooks-Pollock E, Danon L. A flexible method for optimising sharing of healthcare resources and demand in the context of the COVID-19 pandemic. PLoS One. 2020;15(10):e0241027. https://doi.org/10.1371/journal.pone.0241027.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

HW is supported by the Feuer International Scholarship in Artificial Intelligence. TAM, IH, and SB acknowledge funding from the Medical Research Council (MRC) Centre for Global Infectious Disease Analysis (reference no. MR/R015600/1), jointly funded by the UK MRC and the UK Foreign, Commonwealth & Development Office (FCDO), under the MRC/FCDO Concordat agreement and is also part of the EDCTP2 programme supported by the European Union. SB and SF were support by the MRC to undertake research on COVID-19 (MR/V038109/1). SB acknowledges funding from The Academy of Medical Sciences (SBF004/1080), Bill & Melinda Gates Foundation (CRR00280) and Imperial College Healthcare NHS Trust - BRC Funding (RDA02). JMD is supported by an Independent Fellowship funded by Research England’s Expanding Excellence in England (E3) fund. SJV and BAM are supported by The Alan Turing Institute (Engineering and Physical Sciences Research Council (EPSRC) grant EP/N510129/). SJV is also supported by the University of Warwick IAA (Impact Acceleration Account) funding. SF is supported by EPSRC (EP/V002910/1). AD is supported by Wave 1 of the UKRI Strategic Priorities Fund under the EPSRC Grant EP/T001569/1.

Funding

This study received no direct funding. The individual authors’ funder had no role in the design of the study, the analysis, or the formulation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

BAM conceived the study, with input from SJV and JMD. HW, JMD, and SJV prepared the data. HW, TM, and IH undertook the analysis, with input from SJV, SF, and BAM. All authors (HW, TM, IW, JMD, SD, CP, AD, SB, SF, BAM, SJV) contributed to the development of the visualisations. All authors (HW, TM, IW, JMD, SD, CP, AD, SB, SF, BAM, SJV) contributed to the interpretation of results. BAM, HW, and TM drafted the manuscript. All authors (HW, TM, IW, JMD, SD, CP, AD, SB, SF, BAM, SJV) critically revised the article and approved the final version for submission.

Corresponding author

Correspondence to Bilal A. Mateen.

Ethics declarations

Ethics approval and consent to participate

The study was reviewed and approved by the Warwick Biosciences Research Ethics Committee (BSREC 120/19-20-V1.1) and sponsorship is being provided by University of Warwick (SOC.28/19-20). As a retrospective study of data collected by PHE under their statutory authority, the requirement for individual consent was unnecessary for this project by the BSREC.

Consent for publication

Not applicable.

Competing interests

SJV declares funding from IQVIA and Microsoft. BAM is an employee of the Wellcome Trust and holds a Wellcome funded honorary post at University College London for the purposes of carrying out independent research; the views expressed in this manuscript do not necessarily reflect the views or position of the Wellcome Trust. All other authors declare 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:. eMethods, Figures S1-8

, & Tables S1-6. eMethods – Additional detail on data handling and modelling methods. Fig S1. - Application of inclusion and exclusion criteria to raw data. Fig S2. - Unadjusted mortality rate in the ICU, alongside the total, confirmed COVID and non-COVID national-level occupancies (Top) and percentages of the study cohort in each of the three main ICU occupancy bins used in the primary model (Bottom) across the duration of the study. Fig S3. - Trust-level missingness at varying inclusion thresholds. Fig S4. - Seven-day rolling averages for Age (Top) and Number of Comorbidities (Bottom) across the duration of the study. Fig S5. - Trust-Level Ventilator Bed Occupancy (Based on Surge Capacities) Across England. Fig S6. - Marginal posterior densities for linear population level coefficients in the model, for mechanically ventilated bed occupancy binned as described in the main text, comorbidities, sex and week of admission. Fig S7. - Marginal posterior densities for group level ethnicity intercepts. Fig S8. - The interaction between the baseline availability of beds supporting mechanical ventilation, and occupancy on the day of ICU admission. Table S1. - Marginal posterior densities for occupancy under sensitivity analyses. Table S2. - Comparison of prior choices via their effect on primary model OR estimates. Table S3. - Marginal posterior densities for mortality ORs when occupancy is utilised at each individual’s final outcome date, rather than their date of admission to ICU. Table S4. - Marginal posterior densities for mortality ORs when occupancy is taken as the mean or median occupancy achieved during each patient’s stay. Table S5. - Marginal posterior densities for mortality ORs when occupancy is defined as a ratio relative to baseline capacity. Table S6. - Marginal posterior densities for mortality ORs including baseline bed availability and interaction with occupancy.

Rights and permissions

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wilde, H., Mellan, T., Hawryluk, I. et al. The association between mechanical ventilator compatible bed occupancy and mortality risk in intensive care patients with COVID-19: a national retrospective cohort study. BMC Med 19, 213 (2021). https://doi.org/10.1186/s12916-021-02096-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-021-02096-0

Keywords