MJA
MJA

Population density predicts outcome from out-of-hospital cardiac arrest in Victoria, Australia

Ziad Nehme, Emily Andrew, Peter A Cameron, Janet E Bray, Stephen A Bernard, Ian T Meredith and Karen Smith
Med J Aust 2014; 200 (8): 471-475. || doi: 10.5694/mja13.10856
Published online: 5 May 2014

Abstract

Objective: To examine the impact of population density on incidence and outcome of out-of-hospital cardiac arrest (OHCA).

Design, setting and participants: Data were extracted from the Victorian Ambulance Cardiac Arrest Registry for all adult OHCA cases of presumed cardiac aetiology attended by the emergency medical service (EMS) between 1 January 2003 and 31 December 2011. Cases were allocated into one of five population density groups according to their statistical local area: very low density (≤ 10 people/km2), low density (11–200 people/km2), medium density (201–1000 people/km2), high density (1001–3000 people/km2), and very high density (> 3000 people/km2).

Main outcome measures: Survival to hospital and survival to hospital discharge.

Results: The EMS attended 27 705 adult presumed cardiac OHCA cases across 204 Victorian regions. In 12 007 of these (43.3%), resuscitation was attempted by the EMS. Incidence was lower and arrest characteristics were consistently less favourable for lower population density groups. Survival outcomes, including return of spontaneous circulation, survival to hospital and survival to hospital discharge, were significantly poorer in less densely populated groups (< 0.001 for all comparisons). When compared with very low density populations, the risk-adjusted odds ratios of surviving to hospital discharge were: low density, 1.88 (95% CI, 1.15–3.07); medium density, 2.49 (95% CI, 1.55–4.02); high density, 3.47 (95% CI, 2.20–5.48) and very high density, 4.32 (95% CI, 2.67–6.99).

Conclusion: Population density is independently associated with survival after OHCA, and significant variation in the incidence and characteristics of these events are observed across the state.

Out-of-hospital cardiac arrest (OHCA) remains a significant national and global health problem, and is estimated to affect almost 25 000 Australians annually.1 Timely access to emergency medical service (EMS) treatment remains an integral aspect of the chain of survival, and has an important role in determining the provision and performance of EMSs across both urban and regional communities.2 Regional variation in the characteristics and outcomes of OHCA has been described internationally, although few authors have examined the contribution of population density on their reported outcomes.

Previous reports from Australia have found that population factors contribute to the variation in outcomes observed for conditions such as stroke,3 acute myocardial infarction4 and major trauma.5 Similarly, OHCA outcomes have been examined for spatial relationships across areas of Europe,6,7 North America8,9 and Japan;10 however, methods and findings have been inconsistent. In Australia, several authors have shown poorer survival outcomes for OHCA patients in regional areas, although these reports failed to account for the potential effect of population density on survival.11,12

In this study, we examine the contribution of population density to the reported incidence, characteristics and outcome of OHCA in Victoria, Australia.

Methods

Study design

We retrospectively extracted data from the Victorian Ambulance Cardiac Arrest Registry (VACAR) for all non-EMS witnessed OHCA events between 1 January 2003 and 31 December 2011. Adult patients aged more than 17 years who suffered a cardiac arrest of presumed cardiac aetiology were included in the sample. Ethics approval for this study was granted by the Monash University Human Research Ethics Committee.

Setting

As at 30 June 2007, about 5.2 million people inhabited the state of Victoria, Australia, covering a geographical area of 227 000 km2.13 The state comprises 205 inhabited statistical local areas (SLAs) ranging in population density from 0.2 to 7000 people/km2, and a further four uninhabited SLAs.13 The vast majority of its population, about 3.8 million people, reside in the city of Melbourne.

Emergency medical service

The state operates a single EMS comprising about 2500 paramedics who respond in a two-tier system. Advanced life support (ALS) paramedics are authorised to undertake laryngeal mask airway insertion and administer intravenous adrenaline (epinephrine) during cardiac arrest. Mobile intensive care ambulance (MICA) paramedics provide a higher level of clinical care, performing endotracheal intubation, including rapid sequence intubation, and administer a wider range of cardiotropic medications.

ALS and MICA paramedics are dispatched concurrently to suspected cardiac arrest events in the community. A first responder program for early defibrillation by firefighters and volunteer community teams operates in areas of Melbourne and select rural communities across the state.14 For patients with cardiac arrest, paramedic treatment guidelines follow the recommendations of the Australian Resuscitation Council.15 Resuscitation efforts can be withheld by paramedics in accordance with clinical practice guidelines, and include circumstances where there is strong evidence of prolonged cardiac arrest, or where patient injuries are incompatible with life.

Data sources

The VACAR has been described previously.16 Electronically recorded patient data are synchronised daily with a central organisational clinical database. Potential cases are identified by VACAR using a highly sensitive search strategy, and individually screened for eligibility. All eligible cases are reviewed and entered into the registry according to the Utstein requirements.17 Hospital follow-up data are obtained from hospital medical records in as many as 99% of all transported cases.

Definitions

Cardiac arrests are presumed to be of cardiac aetiology unless otherwise stated on the patient care record. The variable “EMS response time” is defined as the time from emergency call to arrival of the first emergency medical response team to the scene. Bystander cardiopulmonary resuscitation (CPR) is defined as any attempt at chest compressions, with or without ventilations. Survival to hospital is defined as a palpable carotid pulse on arrival at hospital as documented on the patient care record.

Geospatial and statistical analyses

Spatial data of included cases were defined by region according to 209 SLAs. Geospatial mapping software MapInfo Professional version 11.5 (Pitney Bowes Software) was used to assign population characteristics according to the longitude and latitude coordinates of the case scene location. Scene address was used to generate approximate location coordinates in cases with missing data. Age-stratified population characteristics for OHCA cases were obtained from the Australian Bureau of Statistics for population estimates at 30 June 2007.13

For the analysis, cases were stratified into five population density categories according to the population density of their immediate SLA: very low density (≤ 10 people/km2), low density (11–200 people/km2), medium density (201–1000 people/km2), high density (1001–3000 people/km2), and very high density (> 3000 people/km2). Patient characteristics and survival outcomes were compared across density groups using the Pearson χ2 test or the Kruskal–Wallis test as appropriate. The Victorian standard population as at 30 June 2007 was used to calculate the adjusted incidences in each region.

The effect of population density on survival outcome was assessed using a logistic regression model in the overall EMS-treated population and the ventricular fibrillation or pulseless ventricular tachycardia (VF/VT) subgroup. Multivariate models were adjusted for factors with known associations to survival, including: EMS response time, age, sex, public location, the event being witnessed by a bystander, bystander CPR and VF/VT (EMS-treated population only). Independent variables were chosen on the basis of previous reports,16,18 and underwent testing in our models to evaluate potential interaction effects with survival (eg, EMS response time and population density). The very low density group was expressed as the reference category in all models. Effect sizes were presented as risk-adjusted odds ratios (AOR) and 95% CIs. P less than 0.05 was considered statistically significant. The predicted probability of survival was graphically depicted after conditional standardisation. Statistical analyses were undertaken using PASW Statistics 18 (SPSS).

Results

The EMS attended 27 705 adult presumed cardiac OHCA cases across 204 Victorian SLAs (Appendix 1). In 12 007 of these (43.3%), the EMS attempted resuscitation (the EMS-treated sample).

Sample characteristics and incidence rates

The statewide unadjusted incidence in the EMS-attended sample was 76.7 per 100 000 adult person-years and in the EMS-treated sample it was 33.3 per 100 000 adult person-years. The age-adjusted incidence rate, mortality rate and case fatality rate also varied across population densities (Appendix 2).

A lower median age was observed in the lower population density groups when compared with higher density groups (P < 0.001, Box 1). EMS response time and the median distance of the nearest EMS station from the scene were significantly higher in the lower density populations. Arrests in public locations and aged care facilities were more frequent in higher density groups. With the exception of the very low density group, VF/VT accounted for about 40.0% of the first monitored rhythms. While the proportion of bystander-witnessed events was similar between the very high and very low density groups (62.2% v 61.6%, P = 0.962), the proportion of patients receiving bystander CPR varied significantly (45.2% v 59.0%, P < 0.001).

Outcomes

Unadjusted survival outcomes were consistently in favour of high density groups (Box 1). Univariate comparisons of return of spontaneous circulation (40.5% v 15.2%, P < 0.001), survival to hospital (36.9% v 12.3%, P < 0.001), and survival to hospital discharge (15.2% v 2.8%, P < 0.001) were significantly improved in the very high density group when compared with the very low density group. After adjusting for known survival confounders, the AOR of survival to hospital and survival to hospital discharge were strongly associated with density group across both the overall EMS-treated population and the VF/VT subgroup (Box 2). The predicted probability of survival after standardisation (Box 3) and the annual unadjusted case fatality rates (Box 4) were improved with higher density groups.

Discussion

This study identifies that population density is independently associated with survival outcome from OHCA, with trends observed in its incidence and characteristics. The data collected for this study are representative of 99.5% of the inhabited regions of Victoria, and offer a key strength over previous reports, which may have been limited to data collected from select municipalities.6-8

Variation in the incidence of OHCA has been described previously,1,8 with the global variation of EMS-attended OHCA ranging between 24.0 and 186.0 per 100 000 person-years.1 In the current study, the highest incidence of OHCA was found in the very high density region, where the presence of modifiable cardiovascular risk factors is reportedly lower than the Australian average.19 This paradox can be partly explained by a shift in Melbourne's diurnal population where people commuting to areas of work create a sixfold increase in the normal resident population of inner Melbourne.20 This may also explain the increased public nature of arrests in the very high density group. Despite significantly better outcomes in higher density areas, lower density groups were associated with favourable baseline characteristics such as younger age, greater bystander CPR and increased VF/VT (low density group). While these findings are likely to be interdependent, they may also indicate selectivity in the EMS treated population, where poorer response time performance is often mitigated by the presence of younger patients and by bystander-initiated resuscitation.21

While previous authors have coupled poorer survival in regional areas to increased EMS response times, this study shows that the association between population density and survival is independent of known survival confounders. Authors from Japan have described similar findings,10 and suggest that combinations of prehospital and hospital-based factors are contributing to regional disparities in OHCA outcome. For example, our study showed a disparity in the availability of intensive care paramedics in less densely populated areas of Victoria. While the impact of paramedic skill level on OHCA outcome is yet to be well understood, an Australian report conducted in a similar EMS system showed that access to an intensive care paramedic was associated with a significant increase in survival to hospital discharge after OHCA.22 Reduced paramedic exposure to resuscitation practice in lower density areas is also likely to be a contributing factor, although evidence in this area remains scarce and requires further exploration. In addition, the frequency of preventable bystander delays after OHCA has undergone a recent evaluation, which shows that arrests occurring in rural regions were independently associated with emergency call delays and death after OHCA.21

Similarly, hospital-based factors have been associated with regional variation in OHCA outcomes. A recent report from Victoria established that the availability of a 24-hour cardiac interventional facility could partly explain the variation in regional survival from OHCA.23 The greatest survival benefit was observed in patients transported to two centralised trauma hospitals situated in the very high density population of Melbourne, and may partly explain the increased survival observed in this region. While the importance of regionalised trauma centres has been well established in Victoria, the implementation of centralised care pathways for OHCA is less certain.24 In fact, few randomised controlled trials have been conducted to examine the benefit of system pathways that improve access to sophisticated emergency teams and optimal post-arrest practices.2,24 These clinical trials are urgently needed, and will help inform the contribution of regionalised post-arrest care on OHCA outcomes in Australia.

We acknowledge several limitations in this study. First, its retrospective design carries inherent limitations. Higher levels of missing survival data were present in the lower density groups and may have introduced some bias. Second, our population density categories were arbitrarily chosen and may not be generalisable in cross-nation comparisons. Finally, other confounding factors such as transport time to hospital and paramedic skill level were considered, but not included in our multivariate analyses. Cases not involving intensive care paramedics were presumed to be overrepresented by futile attempts at resuscitation, and excluding them reduces the possibility of confounding in our models. Similarly, while transport time offers an appropriate surrogate of time to definitive care, its inclusion in multivariate models provides little prognostic information,25 and excludes patients whose resuscitation efforts are ceased at scene (greater in lower density areas).

Nevertheless, this study provides key insight into the variation in OHCA outcomes observed across a large Australian community, and indicates that population density is independently associated with survival outcome. However, the causative factors associated with the observed variation in outcomes remain speculative and require further investigation. Delineating prehospital and inhospital factors associated with improved outcomes in densely populated areas can be used to inform practice recommendations and service provision in sparse populations.

1 Patient characteristics and survival outcomes across population densities after adult out-of-hospital cardiac arrest of presumed cardiac aetiology treated by the emergency medical service (EMS) in Victoria, 2003–2011*

  

Population density


 

Characteristics and outcomes

Statewide

Very low

Low

Medium

High

Very high

P

Total sample

12 007

842

1709

1993

6017

1446

 

Median age in years (IQR)

70 (58–80)

67 (56–78)

68 (56–77)

68 (56–78)

72 (60–81)

71 (59–80)

< 0.001

Men

8460 (70.5%)

626 (74.3%)

1201 (70.3%)

1411 (70.8%)

4148 (68.9%)

1074 (74.3%)

< 0.001

Median EMS response time in minutes (IQR)

8.0 (6.0–10.5)

13.0 (8.0–21.0)

9.1 (7.0–13.0)

9.0 (7.0–11.7)

7.5 (6.0–9.0)

7.0 (5.0–9.0)

< 0.001

Patients with unknown EMS response time

44 (0.4%)

11 (1.3%)

13 (0.8%)

11 (0.6%)

8 (0.1%)

1 (0.1%)

Median distance between scene and nearest EMS station in km (IQR)

2.7 (1.6–4.0)

7.3 (1.4–18.6)

3.2 (1.8–6.7)

3.0 (1.8–4.5)

2.7 (1.7–3.7)

2.0 (1.4–2.6)

< 0.001

Location of arrest

       

Home

8428 (70.2%)

629 (74.7%)

1244 (72.8%)

1435 (72.0%)

4273 (71.0%)

847 (58.6%)

< 0.001

Nursing home

935 (7.8%)

23 (2.7%)

85 (5.0%)

135 (6.8%)

577 (9.6%)

115 (8.0%)

< 0.001

Public place

2557 (21.3%)

169 (20.1%)

360 (21.1%)

408 (20.5%)

1140 (18.9%)

480 (33.2%)

< 0.001

Other

87 (0.7%)

21 (2.5%)

20 (1.2%)

15 (0.8%)

27 (0.4%)

4 (0.3%)

< 0.001

First monitored rhythm

       

VF/VT

4780 (39.8%)

299 (35.5%)

744 (43.5%)

830 (41.6%)

2332 (38.8%)

575 (39.8%)

< 0.001

Asystole

4567 (38.0%)

385 (45.7%)

646 (37.8%)

742 (37.2%)

2266 (37.7%)

528 (36.5%)

< 0.001

Pulseless electrical activity

2476 (20.6%)

108 (12.8%)

303 (17.7%)

402 (20.2%)

1342 (22.3%)

321 (22.2%)

< 0.001

Not shockable

94 (0.8%)

13 (1.5%)

3 (0.2%)

6 (0.3%)

55 (0.9%)

17 (1.2%)

< 0.001

Unknown — pulseless

72 (0.6%)

28 (3.3%)

10 (0.6%)

9 (0.5%)

20 (0.3%)

5 (0.3%)

< 0.001

Unknown

18 (0.1%)

9 (1.1%)

3 (0.2%)

4 (0.2%)

2 (0.0%)

0

Highest skill level on scene

       

Basic life support

643 (5.4%)

174 (20.7%)

116 (6.8%)

105 (5.3%)

199 (3.3%)

49 (3.4%)

< 0.001

Advanced life support

907 (7.6%)

200 (23.8%)

113 (6.6%)

141 (7.1%)

368 (6.1%)

85 (5.9%)

< 0.001

Intensive care

9091 (75.7%)

361 (42.9%)

1238 (72.4%)

1498 (75.2%)

4818 (80.1%)

1176 (81.3%)

< 0.001

Unknown

1366 (11.4%)

107 (12.7%)

242 (14.2%)

249 (12.5%)

632 (10.5%)

136 (9.4%)

< 0.001

Bystander witnessed

7316 (60.9%)

519 (61.6%)

1070 (62.6%)

1267 (63.6%)

3561 (59.2%)

899 (62.2%)

< 0.001

Unknown

69 (0.6%)

12 (1.4%)

15 (0.9%)

12 (0.6%)

24 (0.4%)

6 (0.4%)

Bystander CPR performed

5803 (48.3%)

497 (59.0%)

924 (54.1%)

1029 (51.6%)

2699 (44.9%)

654 (45.2%)

< 0.001

Unknown

580 (4.8%)

36 (4.3%)

75 (4.4%)

90 (4.5%)

284 (4.7%)

95 (6.6%)

Scene outcomes

       

Efforts ceased at scene

7397 (61.6%)

632 (75.1%)

1089 (63.7%)

1250 (62.7%)

3615 (60.1%)

811 (56.1%)

< 0.001

Transport with ROSC

3625 (30.2%)

104 (12.4%)

401 (23.5%)

554 (27.8%)

2046 (34.0%)

52 (36.0%)

< 0.001

Transport with CPR ongoing

985 (8.2%)

106 (12.6%)

219 (12.8%)

189 (9.5%)

356 (5.9%)

115 (8.0%)

< 0.001

ROSC at any time

4272 (35.6%)

128 (15.2%)

469 (27.4%)

659 (33.1%)

2430 (40.4%)

586 (40.5%)

< 0.001

Survival to hospital

3655 (30.7%)

101 (12.3%)

408 (24.2%)

559 (28.3%)

2058 (34.4%)

529 (36.9%)

< 0.001

Unknown

102 (0.8%)

20 (2.4%)

20 (1.2%)

17 (0.9%)

34 (0.6%)

11 (0.8%)

Survival to hospital discharge

1176 (10.0%)

23 (2.8%)

117 (7.1%)

171 (8.8%)

649 (10.9%)

216 (15.2%)

< 0.001

Unknown

218 (1.8%)

27 (3.2%)

56 (3.3%)

41 (2.1%)

72 (1.2%)

22 (1.5%)


* Values are numbers (%) of patients unless otherwise indicated.  Population density categories were defined as: very low density, ≤ 10 people/km2; low density, 11–200 people/km2; medium density, 201–1000 people/km2; high density, 1001–3000 people/km2; and very high density, > 3000 people/km2.  Survival figures exclude missing data. CPR = cardiopulmonary resuscitation. IQR = interquartile range. ROSC = return of spontaneous circulation. VF/VT = ventricular fibrillation or pulseless ventricular tachycardia.

2 Crude and adjusted odds ratios (and 95% CI) for effect of population density on risk-adjusted odds of survival (95% CI) in population treated by the emergency medical service (EMS) and in ventricular fibrillation or pulseless ventricular tachycardia (VF/VT) subgroup, after adult out-of-hospital cardiac arrest of presumed cardiac aetiology in Victoria, 2003–2011*

 

EMS-treated population


 

VF/VT subgroup


Population density (people/km2)

Crude survival to hospital

Adjusted survival to hospital

Crude survival to hospital discharge

Adjusted survival to hospital discharge

 

Crude survival to hospital

Adjusted survival to hospital

Crude survival to hospital discharge

Adjusted survival to hospital discharge


Very low (≤ 10)

1.00

    

1.00

   

Low (11–200)

2.27 (1.80–2.88)

2.19 (1.69–2.84)

2.62 (1.66–4.14)

1.88 (1.15–3.07)§

 

1.78 (1.30–2.45)

1.87 (1.34–2.60)

1.93 (1.19–3.13)§

1.60 (0.96–2.65)

Medium (201–1000)

2.82 (2.24–3.54)

2.86 (2.21–3.69)

3.31 (2.12–5.20)

2.49 (1.55–4.02)

 

2.48 (1.82–3.39)

2.53 (1.83–3.50)

2.63 (1.64–4.21)

2.12 (1.30–3.45)§

High (1001–3000)

3.74 (3.02–4.64)

4.15 (3.26–5.29)

4.22 (2.77–6.44)

3.47 (2.20–5.48)

 

3.92 (2.94–5.24)

4.10 (3.01–5.56)

3.73 (2.39–5.82)

2.96 (1.86–4.72)

Very high (> 3000)

4.17 (3.30–5.27)

4.37 (3.35–5.69)

6.16 (3.97–9.55)

4.32 (2.67–6.99)

 

4.53 (3.28–6.26)

4.42 (3.14–6.23)

5.94 (3.72–9.50)

3.87 (2.36–6.35)


* Logistic regression analysis adjusted fo EMS response time, age, sex, public location, cardiac arrest witnessed by a bystander, cardiopulmonary resuscitation given by a bystander, and VF/VT (EMS-treated patients only).  Reference category.  P < 0.001. § P < 0.05.

Received 25 June 2013, accepted 26 November 2013

Online responses are no longer available. Please refer to our instructions for authors page for more information.