Advertisement
Original Research| Volume 42, ISSUE 1, P48-53, January 2023

Effect of Air Transport Delay on Mortality in Critical Illness: A Population-Based Cohort Study

Published:October 13, 2022DOI:https://doi.org/10.1016/j.amj.2022.09.001

      Highlights

      • Air transport delay from remote areas may affect mortality in critical illness.
      • We evaluated fixed wing air ambulance transport of unselected patients to the nearest intensive care unit.
      • Thirty-day mortality ranged from 11% better to 22% worse for each 1 hour of delay.
      • There is a need for additional, larger studies on this topic.

      Abstract

      Objective

      For critically ill patients in remote areas, we assessed the association of transport delay via fixed wing air ambulance on 30-day mortality, excluding interhospital transports.

      Methods

      This was a retrospective cohort analysis of all such adult transports in Manitoba, Canada, over 5.4 years. Causal mediation analysis was used, with the Acute Physiology and Chronic Health Evaluation II Acute Physiology Score at the destination intensive care unit as the mediator. The covariates were age, sex, comorbidities, socioeconomic status, and physiologic variables from the sending site.

      Results

      The primary cohort was composed of 554 patients; 113 (20.4%) died within 30 days. The total transport delay (mean ± standard deviation) was 5.1 ± 1.7 hours. Compared with no delay, the average 5-hour transport delay was associated with an odds ratio for mortality of 1.34 with a 95% confidence interval from 40% lower to 270% higher, with 60% of the influence of total travel time attributable to worsening of patients’ acute physiologic status during the delay in intensive care unit admission due to transport.

      Conclusions

      Although these findings provide insufficient evidence for an effect of fixed wing air transport delay on mortality among critically ill patients, they underscore the need for additional and larger studies on this topic.
      Critically ill patients require highly specialized care that is generally provided in specialized locations, such as the intensive care units (ICUs) of acute care hospitals. Most critically ill patients experience some delay before ICU admission, and there has long been concern that delays result in worse outcomes.
      • Kim J-T
      • Fonarow GC
      • Smith EE
      • et al.
      Treatment with TPA in the golden hour and the shape of the 4.5 hour time-benefit curve in the National United States Get With The Guidelines-Stroke Population.
      This concept is best known as the “golden hour” in trauma care.
      • Tien HC
      • Jung V
      • Pinto R
      • Mainprize T
      • Scales DC
      • Rizoli SB.
      Reducing time-to-treatment decreases mortality of trauma patients with acute subdural hematoma.
      A small amount of evidence suggests that even delays in ICU admission from emergency departments and wards of the same hospital increase mortality.
      • Chalfin D
      • Trzeciak S
      • Likourezos A
      • Baumann B
      • Dellinger R
      • Group DELAY-ED Study
      Impact of delayed transfer of critically ill patients from the emergency department to the intensive care unit.
      ,
      • Young M
      • Gooder V
      • McBride K
      • James B
      • Fisher E.
      Inpatient transfers to the intensive care unit: delays are associated with increased mortality and morbidity.
      This is a complicated issue, subtending different types of critical illness and numerous diagnostic and therapeutic interventions, each combination of which arguably deserves to be evaluated for the nature of its delay-outcomes relationship. For example, in sepsis, the literature is inconsistent regarding the relationship between survival and delayed initial antibiotics
      • Wheeler DS.
      Is the "golden age" of the "golden hour" in sepsis over?.
      ,
      • Sterling SA
      • Miller WR
      • Pryor J
      • Puskarich MA
      • Jones AE.
      The impact of timing of antibiotics on outcomes in severe sepsis and septic shock: a systematic review and meta-analysis.
      and resuscitation.
      • Rivers E
      • Nguyen B
      • Havstad S
      • et al.
      Early goal-directed therapy in the treatment of severe sepsis and septic shock.
      ,
      • Rowan KM
      • Angus DC
      • Bailey M
      • et al.
      Early, goal-directed therapy for septic shock - a patient-level meta-analysis.
      Indeed, whether 1 hour is a critical threshold is uncertain even in trauma care.
      • Hsieh SL
      • Hsiao CH
      • Chiang WC
      • et al.
      Association between the time to definitive care and trauma patient outcomes: every minute in the golden hour matters.
      ,
      • Okada K
      • Matsumoto H
      • Saito N
      • Yagi T
      • Lee M.
      Revision of 'golden hour' for hemodynamically unstable trauma patients: an analysis of nationwide hospital-based registry in Japan.
      The longest treatment delays occur for patients who become critically ill far from the nearest ICU, necessitating critical care air transport. This is not only necessary in large, sparsely populated countries such as Canada and Australia,
      • Powell DG
      • McCallum AL.
      Canadian transport medicine: waypoints and destinations.
      ,
      • Garwood J
      • Wilkinson B
      • Bartholomew H
      • Langford SA
      • O'Connell A
      Air ambulance and emergency retrieval services in Western Australia: caseload analysis over 5 years.
      but even in more modestly sized countries with remote areas, such as Norway, Denmark, and Azerbaijan.
      • Jagtenberg CJ
      • Vollebergh MA
      • Uleberg O
      • Roislien J.
      Introducing fairness in Norwegian air ambulance base location planning.
      • Sorensen LM
      • Rognas L
      • Alstrup K.
      Trauma versus nontrauma patients treated by the Danish helicopter emergency medical service: a register-based study of patient characteristics, hospitalization, and mortality.
      • Ghaffarzad A
      • Ghalandarzadeh A
      • Rahmani F
      • Rajaei Ghafouri R
      • Dorosti F
      • Morteza-Bagi HR
      Helicopter emergency medical services in East Azerbaijan province: assessment of patients' outcome.
      A large and increasing number of regions in the United States lack any ICU beds,

      Kaiser Health News. Millions of older Americans live in counties with no ICU beds as pandemic intensifies. Available at:https://khn.org/news/as-coronavirus-spreads-widely-millions-of-older-americans-live-in-counties-with-no-icu-beds/. Accessed January 9, 2022.

      requiring such transport.
      Air ambulance transport may be by helicopter or fixed wing aircraft. The latter have far longer range, whereas the former can land close to the location of patients, even at the scenes of accidents. Most of the published literature on air transport is restricted to helicopters; a European Union report on emergency transport found so little information about fixed wing transport that it was excluded from consideration.
      Emergency Medical Services Systems in the European Union.
      In addition, most of the literature has been focused on the transport of trauma patients,
      • Tien HC
      • Jung V
      • Pinto R
      • Mainprize T
      • Scales DC
      • Rizoli SB.
      Reducing time-to-treatment decreases mortality of trauma patients with acute subdural hematoma.
      ,
      • Ghaffarzad A
      • Ghalandarzadeh A
      • Rahmani F
      • Rajaei Ghafouri R
      • Dorosti F
      • Morteza-Bagi HR
      Helicopter emergency medical services in East Azerbaijan province: assessment of patients' outcome.
      ,
      • Neagle G
      • Curatolo L
      • Ferris J
      • Donald M
      • Hearns S
      • Corfield AR.
      Epidemiology and location of primary retrieval missions in a Scottish aeromedical service.
      • Schneider AM
      • Ewing JA
      • Cull JD.
      Helicopter transport of trauma patients improves survival irrespective of transport time.
      • Jarman MP
      • Curriero FC
      • Haut ER
      • Pollack Porter K
      • Castillo RC
      Associations of distance to trauma care, community income, and neighborhood median age with rates of injury mortality.
      who compared with air-transported nontrauma patients generally have lower severity of illness and lower burdens of comorbidity.
      • Sorensen LM
      • Rognas L
      • Alstrup K.
      Trauma versus nontrauma patients treated by the Danish helicopter emergency medical service: a register-based study of patient characteristics, hospitalization, and mortality.
      There are multiple reasons to be concerned about the effects on the outcomes of air transport of critically ill patients. The travel and transfer time itself, even with the critical care capabilities afforded in many air ambulances, may delay some diagnostic and treatment modalities and lead to deterioration of the patient's physiologic status. Also, intubation and deep sedation are often performed routinely in air transport
      • George BP
      • Vakkalanka JP
      • Harland KK
      • et al.
      Sedation depth is associated with increased hospital length of stay in mechanically ventilated air medical transport patients: a cohort study.
      despite strong evidence that deep sedation, in the admittedly different setting of the ICU, is harmful and usually unnecessary.
      • Devlin JW
      • Skrobik Y
      • Gelinas C
      • et al.
      Clinical practice guidelines for the prevention and management of pain, agitation/sedation, delirium, immobility, and sleep disruption in adult patients in the ICU.
      This study was performed in the Canadian province of Manitoba. Compared with the large US state of Texas, Manitoba's land area is 82%, and the population density is 14% (population of 1.3 million). The only ICUs in the province are in only 2 population centers, both in the south of the province (ie, Winnipeg with a population of 768,000 and Brandon with a population of 49,000). Most of the province is very sparsely populated and far from those ICUs; in some parts, travel times are substantial even to get to a nursing station or rural hospital lacking an ICU. The straight-line distance is 1,000 km from Winnipeg to the most significant northern population center, Churchill, on the Arctic Ocean. During this study, Manitoba's medical air ambulance transport included a helicopter service (Stars, privately run and primarily used for trauma patients; range = 268 km;) and a fixed wing service (LifeFlight, governmentally run; range 1,000 km). Our analysis only included the fixed wing air ambulance transports; every LifeFlight transport included a nurse and a physician trained in emergency medicine, critical care, or anesthesia.
      We evaluated the influence of fixed wing air transport delays in ICU admission on mortality among unselected critically ill outpatients with initial presentations to remote acute care sites (hospital emergency departments or regional nursing stations). We hypothesized that among such patients who were transported to an ICU, longer transport delay is associated with higher risk-adjusted, 30-day mortality. We chose this cohort, excluding interhospital transfers and helicopter transport for 2 reasons. First, it is the least studied of all air ambulance transports.
      • Alstrup K
      • Rognas L
      • Sollid S
      • Johnsen SP
      • Valentin JB
      • Petersen JA.
      Association of helicopter vs ground emergency medical transportation with 1-year mortality in Denmark.
      More importantly, if our hypothesis is correct, it has important implications for future research. Specifically, it would mandate further studies seeking ways to ameliorate the elevated rates of death resulting from transport delays and, subsequently, the organization of critical care services.

      Methods

      This retrospective cohort study used 5 linked, existing databases (Table 1) from the Canadian province of Manitoba derived from the province's universal, government-funded, single-payer health care system and held at the Manitoba Centre for Health Policy.

      Manitoba Population Research Data Repository University of Manitoba. 2021. Available at:https://umanitoba.ca/manitoba-centre-for-health-policy/data-repository. Accessed January 9, 2022.

      The LifeFlight database included sending and accepting locations and times of initial request for transport from the sending site, aircraft dispatch, departure for sending airport, arrival at sending airport, departure for destination airport, arrival at destination airport, patient arrival at destination hospital, and patient care transferred to personnel at destination hospital. The Winnipeg ICU Database contains detailed clinical data for all ICU admissions, including admission date/time, the total Acute Physiology and Chronic Health Evaluation II (APACHE II) score, and its Acute Physiology Subscore (APS).
      • Knaus WA
      • Draper EA
      • Zimmerman JE.
      APACHE II: a severity of disease classification system.
      It identified the time of the initial admission to any of the 11 adult ICUs in the 6 hospitals in the Winnipeg Health Region during the study period.
      Table 1Databases Used for the Study
      NameContent
      Patient Registry (A)Dates of provincial health coverage and date of death
      Discharge Abstract Database (A)Details of all provincial hospitalizations
      Medical Claims (A)Details of all provincial outpatient visits
      Winnipeg ICU Database (C)Clinical data on all admissions to adult intensive care units in the Winnipeg Health Region
      LifeFlight (A)Details of all provincial, fixed wing medical air transports
      A = administrative database; C = clinical database.
      The following additional data were manually abstracted from photocopied medical records sent with patients from the sending to the receiving sites (Supplemental Appendix 1): facility type (hospital emergency department vs. others), initial triage date/time, the use of life support (as invasive mechanical ventilation and/or continuous infusions of intravenous vasoactive medications), and the worst value before air ambulance arrival of 11 aspects of acute physiologic status (ie, vital signs [heart rate, mean arterial pressure, respiratory rate, and temperature], laboratory values [white blood cell (WBC) count, hemoglobin, sodium, potassium, creatinine, and Pao2/fraction of inspired oxygen [Fio2] ratio), and the Glasgow Coma Scale (GCS) score. The worst values were per APACHE II definitions,
      • Knaus WA
      • Draper EA
      • Zimmerman JE.
      APACHE II: a severity of disease classification system.
      except for the Pao2/FIO2 ratio. For each of these 11 parameters, the worst value in the interval after triage at the sending site and aircraft arrival was chosen and converted to the 0 to 4 ordinal scale per APACHE II definitions,
      • Knaus WA
      • Draper EA
      • Zimmerman JE.
      APACHE II: a severity of disease classification system.
      except for the Pao2/Fio2 ratio. For that parameter, we used the categories as described by Marshall et al
      • Marshall JC
      • Cook DJ
      • Christou NV
      • Bernard GR
      • Sprung CL
      • Sibbald WJ.
      Multiple organ dysfunction score: a reliable descriptor of a complex clinical outcome.
      ; furthermore, when arterial blood gases had not been performed, we imputed it from the oxygen saturation/Fio2 ratio.
      • Brown SM
      • Duggal A
      • Hou PC
      • et al.
      Nonlinear imputation of PaO2/FIO2 from SpO2/FIO2 among mechanically ventilated patients in the ICU: a prospective, observational study.
      Individuals were included if they were Manitoba residents 18 years or older; transported by fixed wing air ambulance with “emergent” priority between November 1, 2010, and March 31, 2016; and their first inpatient location was a high-intensity, adult ICU in the Winnipeg Health Region. Patients were excluded if they were transported from outside of Manitoba, they were an inpatient in any acute care hospital during the 96 hours before destination arrival (thus excluding interhospital transfers), there was > 6 hours between arrival at the destination hospital and ICU admission, or they were obstetric patients. All eligible subjects were included.
      For patients who after ICU admission experienced inter-ICU and/or interhospital transfers, we used the Winnipeg ICU Database and hospitalization data as previously described to reconstruct the entire episodes of posttransport, ICU-containing hospital care,
      • Fransoo R
      • Yogendran M
      • Olafson K
      • Ramsey C
      • McGowan K-L
      • Garland A.
      Constructing episodes of inpatient care: data infrastructure for population-based research.
      which were the unit of measure for this project. For individuals with multiple such episodes, we only included the first one.

      Analysis

      Our outcome was mortality at 30 days after admission to the destination hospital ICU. The exposure of interest was the total transport delay (TTD), which was defined as the interval from the initial request for transport from the sending site until admission to the destination hospital ICU. The regression model covariates included age, sex, Charlson Comorbidity Index
      • Charlson M
      • Pompei P
      • Ales K
      • MacKenzie C.
      A new method of classifying prognostic comorbidity in longitudinal studies: development and validation.
      (calculated using diagnoses listed in the index hospitalization and all inpatient and outpatient claims during the prior 12 months
      • Quan H
      • Sundararajan V
      • Halfon P
      • et al.
      Coding algorithms for defining comorbidities in ICD-9-CM and ICD-10 administrative data.
      and categorized as 0, 1, or ≥ 2), and socioeconomic status (determined by the SEFI-2 [Socio-economic Factor Index] score with lower status indicated by higher scores).
      • Chateau D
      • Metge C
      • Prior H
      • Soodeen R.
      Learning from the census: the Socio-economic Factor Index (SEFI) and health outcomes in Manitoba.
      A key consideration was to appropriately account for the severity of acute illness, 1 of the most influential determinants of short-term mortality in critical illness.
      • Garland A
      • Olafson K
      • Ramsey C
      • Yogendren M
      • Fransoo R.
      Distinct determinants of long-term and short-term survival in critical illness.
      In our ICU database, it is assessed as the APACHE II APS, which is calculated from data from the initial ICU day.
      • Knaus WA
      • Draper EA
      • Zimmerman JE.
      APACHE II: a severity of disease classification system.
      However, a plausible way that transport delays might cause higher mortality is indirectly via physiologic deterioration during transport. Because simply including APS as a regression model covariate could obscure the effect of transport delay, we used causal mediation analysis;
      • Yung Y-F
      • Lamm M
      • Zhang W.
      Paper SAS1991-2018: causal mediation analysis with the CAUSALMED Procedure.
      with APS as the mediator variable it estimates not only the total effect of TTD on the outcome but also partitions it into indirect and direct portions (Supplemental Fig. 1).
      Another analytic challenge was created by the possibility that, on average, patients with longer transport delays (generally from more remote areas) would have more severe acute illness. This would confound the relationship between TTD and mortality, biasing it away from the null. For this reason, we included as covariates the following information derived from the sending facility: type (hospital emergency department vs. nursing station), the sending site triage interval (ie, the interval from the initial triage at the sending site until the air ambulance was requested), the use of life support (invasive mechanical ventilation or vasoactive drugs) before air ambulance arrival, and the 11 physiologic measures listed earlier. Using multiple imputation for missing values of sending site WBC count, creatinine, temperature, and GCS, we generated 100 complete data sets (Supplemental Appendices 1 and 2, Supplemental Table 1, Supplemental Figure 2). The 26 variables used in imputation were 1) the actual values for variables from the sending site; 2) age, sex, and socioeconomic status (determined by the SEFI-2 score), and 3) APACHE II–related parameters from the first day in the ICU at the receiving hospital (ie, heart rate, mean arterial pressure, respiratory rate, temperature, WBC count, hemoglobin, sodium, potassium, creatinine, and GCS). After imputation, we converted the missing values to the 0 to 4 ordinal scale as described previously. Following the study by Lall,
      • Lall R.
      How multiple imputation makes a difference.
      we evaluated the appropriateness of imputation by assessing whether missingness was strongly related to 1 or a just a few of the observed variables.
      A final analytic issue was that our analysis was at risk for model overfitting from including too many independent variables for the number of outcomes.
      • Harrell Jr., FE
      Regression Modeling Strategies: Overfitting.
      Although the well-known rule of at least 10 events/variable has support,
      • Pedruzzi P
      • Concato J
      • Kemper E
      • Holford TR
      • Feinstein AR.
      A simulation study of the number of events per variable in logistic regression analysis.
      other work indicates that 5 to 9 events/variable may be sufficient.
      • Vittinghoff E
      • McCulloch CE.
      Relaxing the rule of ten events per variable in logistic and Cox regression.
      Accordingly, we performed variable reduction on the 11 sending site physiologic variables via principal component analysis
      • Harrell Jr., FE
      Regression Modeling Strategies: Data Reduction.
      using the mean correlation method of van Ginkel et al
      • van Ginkel JR
      • Kroonenberg PM.
      Using generalized procrustes analysis for multiple imputation in principal component analysis.
      as recommended for combining it with multiple imputation (Supplemental Appendix 1). The goal was to include as much of the information content of those 11 variables as possible while avoiding model overfitting for the 113 events in our primary analysis. In line with the Kaiser-Guttman rule,
      • Dunteman GH.
      Principal Components Analysis.
      we retained the 4 principal components with eigenvalues > 1, which for our 113 outcomes resulted in 14 model covariates, providing 8.1 outcomes/variable. In a sensitivity analysis, we included only the single most influential component, providing 10.3 outcomes/variable.
      In a sensitivity analysis, we restricted consideration to the 486 records without any missing covariate data elements. We separately performed principal component analysis on this smaller data set.
      Statistical analysis was performed using SAS Version 9.4 (SAS Institute, Cary, NC). The primary interpretation used 95% confidence intervals (CIs).
      • Simon R.
      Confidence intervals for reporting results of clinical trials.

      Results

      Over the 65 months of the study, there were 737 unique adults transported by fixed wing air ambulance, 630 of whom met all the inclusion and exclusion criteria; 85 were then excluded for missing physiologic parameters from the sending site other than temperature, WBC count, creatinine, and GCS (Supplemental Appendix 1, Fig. 1). Of the remaining 554 records, 468 had no additional missing parameters, whereas 68 (12.3%) required imputation of 1 or more of sending site temperature, WBC count, creatinine, and GCS. GCS had the most missing cases (32 representing 5.8% of records). The results of the logistic regression model for missingness of any of temperature, WBC count, creatinine, or GCS was highly significant (χ2 = 41.3, P = .0075), with the dominant predictor being that the patient was sent from a medical site other than the emergency department of a hospital (Supplemental Table 3; odds ratio = 3.3, P = .001).
      Figure 1
      Figure 1The Consolidated Standards of Reporting Trials diagram.
      The baseline characteristics of this cohort (Table 2) showed that most were males transported from hospital emergency departments who received life support at the sending site. The mean time from the initial triage to requesting air transport was 13.2 hours. The mean TTD was 5.1 hours, ranging from 1 to 12 hours. The air travel time and the time on the ground at the sending site averaged 1.7 and 1.9 hours, respectively, each approximately one third of the TTD.
      Table 2Baseline Characteristics and Outcomes of the Final Cohort
      VariableValueN
      Sociodemographic variables
       Age56.0 ± 16.1554
       Female sex227 (41.0)554
       Charlson Comorbidity Index554
       0221 (39.9)
       1143 (25.8)
       ≥ 2190 (34.3)
      SEFI-2 (Socioeconomic Index)
      Population distribution has a mean of 0 and a standard deviation of 1.0; higher values indicate lower status.
      1.00 ± 1.36554
      Variables from the sending site
       Facility type554
        Hospital emergency department469 (84.7)
        Regional nursing station85 (15.3)
       Life support in the sending site
      Invasive mechanical ventilation and/or continuous infusions of intravenous vasoactive agent(s).
      373 (67.3)554
       Sending site triage interval (h)
      The interval from the initial presentation at the sending site until air transport was requested.
      13.2 ± 18.8554
       Time of day of initial arrival at sending site (24-h day)12.8 ± 6.2554
       Heart rate (beats/min)104 ± 42554
       Mean arterial blood pressure (mm Hg)72.3 ± 30.3554
       Respiratory rate (breaths/min)27 ± 10554
       Temperature (°C)36.4 ± 1.7538
       Serum sodium (mmol/L)138 ± 10554
       Serum potassium (mmol/L)4.3 ± 1.1554
       Serum creatinine (μmol/L)199.2 ± 242.6533
       WBC count (1,000/μL)14.4 ± 8.1531
       Hematocrit (%)37.2 ± 8.8554
       Pao2/FIO2 ratio332 ± 158554
       Glasgow Coma Scale score8.9 ± 5.6522
      Air transport-related variables
       Straight-line travel distance (km)454 ± 185554
       Total transport time/delay (h)5.1 ± 1.7554
       Air travel time (h)1.7 ± 0.7554
       Time on ground during pickup (h)1.9 ± 0.9554
       Time from landing at destination airport to ICU admission (h)1.6 ± 1.7554
      Variables from admitting ICU
       Most responsible hospital diagnosis category554
        Cardiovascular257 (46.4)
        Infection82 (14.8)
        Respiratory62 (11.2)
        Trauma46 (8.3)
        Neuropsychiatric42 (7.6)
        Poisoning/overdose26 (4.7)
        Others39 (7.0)
       APACHE II Acute Physiology Score11.1 ± 6.8554
       Total APACHE II score18.3 ± 9.8554
       Heart rate (beats/min)90 ± 35554
       Mean arterial blood pressure (mm Hg)76.1 ± 26.8554
       Respiratory rate (breaths/min)23 ± 9554
       Temperature (°C)37.1 ± 1.2554
       Serum sodium (mmol/L)138.3 ± 6.4554
       Serum potassium (mmol/L)4.1 ± 1.0554
       Serum creatinine (μmol/L)170.7 ± 213.2554
       WBC count (1,000/μL)14.3 ± 8.2554
       Hematocrit (%)33.9 ± 7.9554
       Glasgow Coma Scale score10.2 ± 2.4554
      Outcomes
       Hospital length of stay (h)23.9 ± 42.8554
       30-day mortality113 (20.4)554
      APACHE II = Acute Physiology and Chronic Health Evaluation II; Fio2 = fraction of inspired oxygen; ICU = intensive care unit.
      The values are mean ± standard deviation or n (%).
      a Population distribution has a mean of 0 and a standard deviation of 1.0; higher values indicate lower status.
      b Invasive mechanical ventilation and/or continuous infusions of intravenous vasoactive agent(s).
      c The interval from the initial presentation at the sending site until air transport was requested.
      The 4 retained principal components included 57% of the variation (information content) of the original 11 sending site physiologic variables (Supplemental Tables 1 and 2).
      Primary analysis on the primary cohort of 554 air-transported patients (Tables 3 and 4, Supplemental Table 4) gave a point estimate indicating higher 30-day mortality with an odds ratio of 1.06 per hour of transport delay (95% CI, 0.89-1.22). On average, 60% of the influence of the total travel time was attributed to worsening of patients’ acute physiologic status during the delay in ICU admission due to transport; this is also consistent with a systematically higher APS with a longer total transport time (Table 4). The results were essentially unchanged including only the most influential principal component (Supplemental Table 5). Sensitivity analysis including only the 486 complete records showed similar results (Table 3); the odds ratio for 30-day mortality was 1.05 (95% CI, 0.88-1.24) for the 486 complete records and 1.06 (95% CI, 0.89-1.22) for the 554 records including imputed data.
      Table 3Primary Mediation Analysis Results: The Effect of Each 1 Hour of Total Travel Time on 30-Day Mortality
      N = 554, Including Imputed Data (113 With the Outcome)N = 486, Complete Records Only (91 With the Outcome)
      Odds Ratio95% CIOdds Ratio95% CI
      Including sending site covariates
       Total effect1.060.89-1.221.050.88-1.24
       Direct effect1.020.87-1.181.020.85-1.19
       Indirect effect1.030.99-1.071.031.00-1.08
       Percentage mediated59.867.2
      Excluding sending site covariates
       Total effect1.31
      The CI did not include the null result.
      1.12-1.531.26
      The CI did not include the null result.
      1.08-1.46
       Direct effect1.130.97-1.321.110.96-1.30
       Indirect effect1.16
      The CI did not include the null result.
      1.10-1.241.14
      The CI did not include the null result.
      1.07-1.21
       Percentage mediated58.857.2
      CI = confidence interval.
      a The CI did not include the null result.
      Table 4Ordinary Least Squares Regression Model of the Mediator Variable (Acute Physiology and Chronic Health Evaluation II Acute Physiology Score) at the Destination Site Intensive Care Unit (N = 554)
      ParameterIncluding Sending Site DataExcluding Sending Site Data
      Coefficient95% CIP ValueCoefficient95% CIP Value
      Age (per year)0.006−0.021 to 0.033.66−0.029−0.063 to 0.006.11
      Female sex (compared with male)0.543−0.293 to 1.378.201.1650.101-2.229.03
      Charlson Comorbidity Index
       0ReferenceReference
       10.809−0.207 to 1.825.121.9270.618-3.237.0039
       ≥ 21.5500.542-2.561.0034.3913.145-5.636<.0001
      SEFI-2 (Socioeconomic Index)
      Population distribution has a mean of 0 and standard deviation of 1.0; higher values indicate a lower status.
      0.196−0.166 to 0.558.290.6870.267-1.107.0014
      Sending site type
       Hospital emergency departmentReference
       Regional nursing station0.145−1.114 to 1.405.82
      Life support at sending facility
      Invasive mechanical ventilation and/or continuous infusions of intravenous vasoactive agent(s).
      1.7060.505-2.907.005
      Sending site triage interval
      The interval from the initial presentation at the sending site until air transport was requested.
      (per hour)
      0.001−0.021 to 0.023.92
      Principal component 1
      Derived from 11 variables per Supplemental Appendix 1.
      1.9261.687-2.166< .0001
      Principal component 2
      Derived from 11 variables per Supplemental Appendix 1.
      −1.053−1.290 to −0.816< .0001
      Principal component 3
      Derived from 11 variables per Supplemental Appendix 1.
      0.058−0.212 to 0.328.67
      Principal component 4
      Derived from 11 variables per Supplemental Appendix 1.
      −0.330−0.709 to 0.049.088
      Total transport time (per hour)0.2660.012-0.520.040.9720.664-1.280< .0001
      CI = confidence interval.
      a Population distribution has a mean of 0 and standard deviation of 1.0; higher values indicate a lower status.
      b Invasive mechanical ventilation and/or continuous infusions of intravenous vasoactive agent(s).
      c The interval from the initial presentation at the sending site until air transport was requested.
      d Derived from 11 variables per Supplemental Appendix 1.
      The relevance of including the severity of illness at the sending site in the analysis is demonstrated by performing analysis without the sending site variables (Table 3). Those results showed a much greater detrimental impact of travel time, with the total effect of the odds of 30-day mortality being 1.31-fold higher for each hour of transport delay (95% CI, 1.12-1.53).

      Discussion

      We assessed the association of 30-day mortality to delays in transporting unselected critically ill patients from remote areas to ICUs using fixed wing air ambulances. For the average travel delay of 5 hours, the point effect was 33.8% higher odds of death compared with no delay, with a 95% CI ranging from 40% lower to 270% higher; approximately 60% of the point effect is attributable to physiologic deterioration during the travel delay.
      Because these findings did not satisfy the usual criteria for statistical significance, our study provides insufficient evidence of an effect of delay on mortality. However, because our study is not inconsistent with a substantial detrimental effect of longer transport times on mortality, it underscores the need for additional and larger studies on this topic. Consistent findings that such air ambulance transports do not lead to worse outcomes would provide reassurance for the care of critically ill patients transported long distances to the nearest ICU. However, if such a transport does worsen outcomes, it would have important implications for the organization of critical care in remote areas. Specifically, it would indicate the need for some combination of: shortening air transport (eg, using a faster aircraft), providing more stabilization and treatment at the sending site, or reducing the need for long transports by establishing ICUs in remote, underserved areas.
      The relevance of this topic derives from 1) the increasing need for distant transport to the nearest ICU as indicated earlier and 2) an extreme paucity of prior evidence. Only 1 prior study has assessed the influence of air ambulance transport of critically ill patients not limited to trauma patients and not limited to interhospital transports; the Danish study compared the time to death over 1 year among critically ill patients transported to the hospital via helicopter versus ground ambulance when the helicopter was called but unavailable.
      • Alstrup K
      • Rognas L
      • Sollid S
      • Johnsen SP
      • Valentin JB
      • Petersen JA.
      Association of helicopter vs ground emergency medical transportation with 1-year mortality in Denmark.
      Although the reported hazard ratio for helicopter transport was 0.94 (95% C.I. 0.84-1.06), the analysis did not account for actual travel times and included no adjustment for the severity of acute illness. The latter is relevant because among helicopter transports nontrauma patients have a systematically lower severity of illness.
      • Sorensen LM
      • Rognas L
      • Alstrup K.
      Trauma versus nontrauma patients treated by the Danish helicopter emergency medical service: a register-based study of patient characteristics, hospitalization, and mortality.
      Our study has notable strengths. First, it is a population-based study of all such fixed wing transports in a Canadian province over 5 years. By including all types of ICU admissions rather than being restricted to trauma patients, our results are applicable to a wide range of patient presentations. Second, our study used a rigorous statistical approach that avoids incorrect conclusions by performing causal mediation analysis and including measures of severity of acute illness at both the sending and destination sites. A comparison of the top and bottom halves of Table 3 illustrates that failing to include sending site information greatly overestimates the effect of travel delay on mortality.
      Our study also has limitations. Despite including a wide range of covariates in our analysis, residual confounding is always possible, which could be negative confounding.
      • Choi AL
      • Cordier S
      • Weihe P
      • Grandjean P.
      Negative confounding in the evaluation of toxicity: the case of methylmercury in fish and seafood.
      Despite using all available data from an entire province over a number of years, the relatively small number of mortality events (113) could have resulted in a type 2 error. Only considering subjects who reached the ICU alive could have conceivably resulted in survival bias, which for our outcome would be expected to bias toward the null.
      • Shintani AK
      • Girard TD
      • Eden SK
      • Arbogast PG
      • Moons KG
      • Ely EW.
      Immortal time bias in critical care research: application of time-varying Cox regression for observational cohort studies.
      However, we consider this unlikely because during our study interval the LifeFlight database contained only 10 patients who died after the initial request for transport from the sending site but before they could be admitted to the destination ICU. Finally, the results from a single Canadian province might not be applicable to other jurisdictions. In particular, with substantial variation in the aircraft personnel transporting patients,
      • Powell DG
      • McCallum AL.
      Canadian transport medicine: waypoints and destinations.
      it is possible that the Manitoba practice of including a physician and nurse could have resulted in lower mortality compared with other staffing paradigms.

      Appendix. Supplementary materials

      References

        • Kim J-T
        • Fonarow GC
        • Smith EE
        • et al.
        Treatment with TPA in the golden hour and the shape of the 4.5 hour time-benefit curve in the National United States Get With The Guidelines-Stroke Population.
        Circulation. 2017; 135: 128-139
        • Tien HC
        • Jung V
        • Pinto R
        • Mainprize T
        • Scales DC
        • Rizoli SB.
        Reducing time-to-treatment decreases mortality of trauma patients with acute subdural hematoma.
        Ann Surg. 2011; 253: 1178-1183
        • Chalfin D
        • Trzeciak S
        • Likourezos A
        • Baumann B
        • Dellinger R
        • Group DELAY-ED Study
        Impact of delayed transfer of critically ill patients from the emergency department to the intensive care unit.
        Crit Care Med. 2007; 35: 1477-1483
        • Young M
        • Gooder V
        • McBride K
        • James B
        • Fisher E.
        Inpatient transfers to the intensive care unit: delays are associated with increased mortality and morbidity.
        J Gen Intern Med. 2003; 18: 77-83
        • Wheeler DS.
        Is the "golden age" of the "golden hour" in sepsis over?.
        Crit Care. 2015; 19: 447
        • Sterling SA
        • Miller WR
        • Pryor J
        • Puskarich MA
        • Jones AE.
        The impact of timing of antibiotics on outcomes in severe sepsis and septic shock: a systematic review and meta-analysis.
        Crit Care Med. 2015; 43: 1907-1915
        • Rivers E
        • Nguyen B
        • Havstad S
        • et al.
        Early goal-directed therapy in the treatment of severe sepsis and septic shock.
        N Engl J Med. 2001; 345: 1368-1377
        • Rowan KM
        • Angus DC
        • Bailey M
        • et al.
        Early, goal-directed therapy for septic shock - a patient-level meta-analysis.
        N Engl J Med. 2017; 376: 2223-2234
        • Hsieh SL
        • Hsiao CH
        • Chiang WC
        • et al.
        Association between the time to definitive care and trauma patient outcomes: every minute in the golden hour matters.
        Eur J Trauma Emerg Surg. 2022; 48: 2709-2716
        • Okada K
        • Matsumoto H
        • Saito N
        • Yagi T
        • Lee M.
        Revision of 'golden hour' for hemodynamically unstable trauma patients: an analysis of nationwide hospital-based registry in Japan.
        Trauma Surg Acute Care Open. 2020; 5e000405
        • Powell DG
        • McCallum AL.
        Canadian transport medicine: waypoints and destinations.
        CJEM. 2020; 22: S1-S3
        • Garwood J
        • Wilkinson B
        • Bartholomew H
        • Langford SA
        • O'Connell A
        Air ambulance and emergency retrieval services in Western Australia: caseload analysis over 5 years.
        Air Med J. 2020; 39: 35-43
        • Jagtenberg CJ
        • Vollebergh MA
        • Uleberg O
        • Roislien J.
        Introducing fairness in Norwegian air ambulance base location planning.
        Scand J Trauma Resusc Emerg Med. 2021; 29: 50
        • Sorensen LM
        • Rognas L
        • Alstrup K.
        Trauma versus nontrauma patients treated by the Danish helicopter emergency medical service: a register-based study of patient characteristics, hospitalization, and mortality.
        Air Med J. 2021; 40: 325-330
        • Ghaffarzad A
        • Ghalandarzadeh A
        • Rahmani F
        • Rajaei Ghafouri R
        • Dorosti F
        • Morteza-Bagi HR
        Helicopter emergency medical services in East Azerbaijan province: assessment of patients' outcome.
        Ulus Travma Acil Cerrahi Derg. 2021; 27: 427-433
      1. Kaiser Health News. Millions of older Americans live in counties with no ICU beds as pandemic intensifies. Available at:https://khn.org/news/as-coronavirus-spreads-widely-millions-of-older-americans-live-in-counties-with-no-icu-beds/. Accessed January 9, 2022.

      2. Emergency Medical Services Systems in the European Union.
        World Health Organization, Copenhagen, Denmark2008
        • Neagle G
        • Curatolo L
        • Ferris J
        • Donald M
        • Hearns S
        • Corfield AR.
        Epidemiology and location of primary retrieval missions in a Scottish aeromedical service.
        Eur J Emerg Med. 2019; 26: 123-127
        • Schneider AM
        • Ewing JA
        • Cull JD.
        Helicopter transport of trauma patients improves survival irrespective of transport time.
        Am Surg. 2021; 87: 538-542
        • Jarman MP
        • Curriero FC
        • Haut ER
        • Pollack Porter K
        • Castillo RC
        Associations of distance to trauma care, community income, and neighborhood median age with rates of injury mortality.
        JAMA Surg. 2018; 153: 535-543
        • George BP
        • Vakkalanka JP
        • Harland KK
        • et al.
        Sedation depth is associated with increased hospital length of stay in mechanically ventilated air medical transport patients: a cohort study.
        Prehosp Emerg Care. 2020; 24: 783-792
        • Devlin JW
        • Skrobik Y
        • Gelinas C
        • et al.
        Clinical practice guidelines for the prevention and management of pain, agitation/sedation, delirium, immobility, and sleep disruption in adult patients in the ICU.
        Crit Care Med. 2018; 46: e825-e873
        • Alstrup K
        • Rognas L
        • Sollid S
        • Johnsen SP
        • Valentin JB
        • Petersen JA.
        Association of helicopter vs ground emergency medical transportation with 1-year mortality in Denmark.
        JAMA Netw Open. 2021; 4e2033318
      3. Manitoba Population Research Data Repository University of Manitoba. 2021. Available at:https://umanitoba.ca/manitoba-centre-for-health-policy/data-repository. Accessed January 9, 2022.

        • Knaus WA
        • Draper EA
        • Zimmerman JE.
        APACHE II: a severity of disease classification system.
        Crit Care Med. 1985; 13: 818-829
        • Marshall JC
        • Cook DJ
        • Christou NV
        • Bernard GR
        • Sprung CL
        • Sibbald WJ.
        Multiple organ dysfunction score: a reliable descriptor of a complex clinical outcome.
        Crit Care Med. 1995; 23: 1638-1652
        • Brown SM
        • Duggal A
        • Hou PC
        • et al.
        Nonlinear imputation of PaO2/FIO2 from SpO2/FIO2 among mechanically ventilated patients in the ICU: a prospective, observational study.
        Crit Care Med. 2017; 45: 1317-1324
        • Fransoo R
        • Yogendran M
        • Olafson K
        • Ramsey C
        • McGowan K-L
        • Garland A.
        Constructing episodes of inpatient care: data infrastructure for population-based research.
        BMC Med Res Methodol. 2012; 12: 133
        • Charlson M
        • Pompei P
        • Ales K
        • MacKenzie C.
        A new method of classifying prognostic comorbidity in longitudinal studies: development and validation.
        J Chronic Dis. 1987; 40: 373-383
        • Quan H
        • Sundararajan V
        • Halfon P
        • et al.
        Coding algorithms for defining comorbidities in ICD-9-CM and ICD-10 administrative data.
        Med Care. 2005; 43: 1130-1139
        • Chateau D
        • Metge C
        • Prior H
        • Soodeen R.
        Learning from the census: the Socio-economic Factor Index (SEFI) and health outcomes in Manitoba.
        Can J Public Health. 2012; 103: S23-S27
        • Garland A
        • Olafson K
        • Ramsey C
        • Yogendren M
        • Fransoo R.
        Distinct determinants of long-term and short-term survival in critical illness.
        Intensive Care Med. 2014; 40: 1097-1105
        • Yung Y-F
        • Lamm M
        • Zhang W.
        Paper SAS1991-2018: causal mediation analysis with the CAUSALMED Procedure.
        SAS Institute Inc, 2019 (Available at) (Accessed January 29)
        • Lall R.
        How multiple imputation makes a difference.
        Polit Anal. 2016; 24: 414-433
        • Harrell Jr., FE
        Regression Modeling Strategies: Overfitting.
        Springer, New York, NY2001: 60-61
        • Pedruzzi P
        • Concato J
        • Kemper E
        • Holford TR
        • Feinstein AR.
        A simulation study of the number of events per variable in logistic regression analysis.
        J Clin Epidemiol. 1996; 49: 1373-1379
        • Vittinghoff E
        • McCulloch CE.
        Relaxing the rule of ten events per variable in logistic and Cox regression.
        Am J Epidemiol. 2006; 165: 710-718
        • Harrell Jr., FE
        Regression Modeling Strategies: Data Reduction.
        Springer, New York2001: 66-74
        • van Ginkel JR
        • Kroonenberg PM.
        Using generalized procrustes analysis for multiple imputation in principal component analysis.
        J Classif. 2014; 31: 242-269
        • Dunteman GH.
        Principal Components Analysis.
        Sage Publications, Inc, 1989: 69
        • Simon R.
        Confidence intervals for reporting results of clinical trials.
        Ann Intern Med. 1986; 105: 429-435
        • Choi AL
        • Cordier S
        • Weihe P
        • Grandjean P.
        Negative confounding in the evaluation of toxicity: the case of methylmercury in fish and seafood.
        Crit Rev Toxicol. 2008; 38: 877-893
        • Shintani AK
        • Girard TD
        • Eden SK
        • Arbogast PG
        • Moons KG
        • Ely EW.
        Immortal time bias in critical care research: application of time-varying Cox regression for observational cohort studies.
        Crit Care Med. 2009; 37: 2939-2945