Organizations: Institute for Sustainable Industries and Liveable Cities (ISILC), The Centre of Excellence in Paramedicine (CoEP), Victoria University, Melbourne, Australia · Center for Research and Evaluation Ambulance Victoria, Blackburn North, Victoria, Australia
Ambulance ramping, the delay between hospital arrival and patient handover, is a critical operational bottleneck in Emergency Medical Services (EMS), yet its systemic magnitude and dynamics remain inadequately characterised at scale. This paper quantifies the scale, trajectory, and operational correlates of ramping across an entire statewide EMS system, analysing 2,850,575 ambulance attendances in Victoria, Australia from January 2020 to March 2024 using an Exploratory Data Analysis (EDA). After systematic preprocessing, an analytical cohort of 2,026,569 Emergency Department (ED) transports across 59 hospitals with ED and 79 Local Government Areas (LGA) was examined through interval decomposition, Pareto concentration, hourly cross-correlation, hospital arrival concurrency and priority-stratified operational comparisons. Cumulative Ambulance Hours Lost (AHL) totalled 1,491,127 hours, equivalent to approximately 96 ten-hour ambulance shift lost every day of the study window. Ten of 59 hospitals account for 57.8% of lost hours from 50.9% of cases. Annual losses rose 57% to a 2022 peak while transported demand fell 3.7%, indicating deterioration in per-case handover rather than growth in demand. Handover duration varies little with patient acuity, but rises monotonically with the number of ambulances arriving at the same hospital in the preceding hour, an effect persisting within every hour of the day. Hourly demand is moderately associated with ramping two to four hours later (r = 0.365). These findings establish the empirical preconditions for hospital-state aware ambulance routing.
Figures & tables
Fig. 1: Seven AV regions: AV’s services are divided into seven regions across the state : the metropolitan (East,West) and five regional and rural regions (Barwon South West, Gippsland, Grampians, Hume and Loddon Mallee).
Fig. 2: Timestamped events in the EMS operational workflow. Seven key events recorded in the dataset (call received, ambulance dispatch, arrival at scene, begin transport, arrival at the medical facility, patient transfer, and back to service) are mapped to the corresponding operational intervals.
Fig. 3: Derivation of the analytical cohorts, with the sample size for each results subsection. Cases not originating from public emergency calls were excluded throughout. Each cohort applies only the requirements of the analyses it supports.
Fig. 4: LGA wise average monthly case loads over the full study period. first 15 LGAs with the most ambulance demand cases are labelled. Shape colors depict LGA regional mapping.
Fig. 5: Duration of service time interval by patient acuity, shown as grouped bars. Intervals completing before the patient’s on-scene assessment, stratified by dispatch code (left panel bars). The intervals occurring after it, stratified by transport priority. Components are not stacked, as percentiles are not additive; the two panels use different case groupings and are not to be read across.
Priority 0
Code 1
Code 2
Code 3
Interval
38,616
1,123,700
630,583
358,998
Med (p90)
Med (p90)
Med (p90)
Med (p90)
Call to dispatch
2.0 (11.7)
2.1 (15.2)
11.6 (59.4)
29.4 (94.5)
Dispatch to scene
8.5 (23.1)
9.5 (21.1)
14.3 (29.6)
20.7 (55.9)
On-scene
25.9 (60.5)
20.8 (38.0)
20.5 (41.7)
23.9 (45.8)
Transport 1
Transport 2
Transport 3
Transport 4
TABLE I: Duration (minutes) of each service-time interval by patient acuity, reported as median and 90th percentile (January 2020-March 2024). Intervals completing before the paramedic’s transport assessment are stratified by dispatch code; those occurring after it by transport priority. Denominators differ because each interval is computed on the subset of cases with a plausible recorded value for that interval.
Fig. 6: Yearly trajectory of Ambulance Hours Lost (AHL) all the transports to ED between Jan 2020-March 2024. Total AHL hours per calendar year, annotated with the equivalent number of 10-hour ambulance shifts lost.
Region
Cases
AHL (h)
Mean
Shift
Ramp (min)
Lost/day
Metro
1,465,470
1,143,320.9
46.8
73.6
Barwon SW
149,332
91,513.9
36.8
5.9
Hume
105,094
72,315.4
41.3
4.7
Loddon Mallee
121,119
71,660.9
35.5
4.6
Gippsland
102,497
60,640.7
35.5
3.9
TABLE II: Regional ambulance hours lost (AHL) by hospital region across the study window.
Year
n_cases
AHL
Mean
Avg Shifts
Ramp(h)
Ramp/day
Lost/day
2020
479,203
255,522.6
32.0
69.8
2021
509,647
373,298.1
43.9
102.3
2022
461,418
401,980.7
52.3
110.1
2023
471,688
382,261.1
48.6
104.7
2024
104,641
78,064.1
44.8
84.9
TABLE III: Yearly trajectory of Ambulance Hours Lost (AHL) across the study window. The cohort comprises all transports taken place during 4.25 years to the 59 hospitals with ED. Shifts lost are expressed in 10-hour shift equivalents. The 2024 row reflects partial year data (January–March, 91 days) and is reported with its daily average rate for comparability with full calendar years.
Fig. 7: Yearly ramping trajectory by patient’s transport priority code: median (left) and 90th percentile (right). Ramping times increased across all transport priorities from 2020 to a peak in 2022, with partial easing in 2023; Transport Priority 4 carries the heaviest tail throughout. 2024 (shaded) covers January–March only
Fig. 8: Hospital wise ambulance hours lost across the study period Jan 2020 - Mar 2024, demonstrates that the first 15 medical facilities bear the 75% AHL while their case share is 68.6%. Out of those 15 facilities, 13 are in the Metro (Refer Fig. 9 ) region, Victoria.
Year
Priority 1
Priority 2
Priority 3
Priority 4
Total
2020
47,957
353,460
49,880
27,906
479,203
2021
53,378
371,670
54,366
30,233
509,647
2022
57,565
338,261
45,904
19,688
461,418
2023
62,745
341,457
46,178
21,308
471,688
2024
14,026
73,365
11,490
5,760
104,641
TABLE IV: Yearly distribution of cases by transport priority(1-4) from Jan 2020-March 2024.
Fig. 9: Pareto analysis: The medical facilities ranked 41-59 contribute a combined 2-3% of the total AHL, their bars are near-invisible. The facilities showing cumulative AHL share (solid line) and cumulative case share (dashed line) as functions of rank.
Rank
Cum. AHL (%)
Cum. Cases (%)
Difference (pp)
5
32.9
28.2
+4.7
10
57.8
50.9
+6.9
15
75.0
68.6
+6.4
20
86.0
80.3
+5.7
25
91.0
86.6
+4.3
30
94.5
91.0
+3.5
TABLE V: Cumulative hospital contributions to AHL across 59 EDs, assessed at key ranking positions. Hospitals are ordered by decreasing total AHL contribution. Difference (percentage points, pp) quantifies the excess of cumulative AHL share over cumulative case share.
Lag
Demand
Arrivals
Lag
Demand
Arrivals
0
0.2372
0.3168
7
0.2031
0.1165
1
0.3130
0.3525
8
0.1274
0.0349
2
0.3526
0.3639
9
0.0452
-0.0437
3
0.3652
0.3524
10
-0.0335
-0.1289
4
0.3532
0.3181
11
-0.1171
-0.2086
5
0.3209
0.2643
12
-0.1974
-0.2674
TABLE VI: Cross-correlation between hourly demand and arrivals with hourly mean ramping.
Fig. 10: Hour of day profile of ambulance demand, ED arrivals and mean ramping duration, averaged across the study window. Demand peak 130,620 cases per hour of day, arrivals peak 125,583, mean ramping peaks 50.65 minutes.
Fig. 11: Hourly 90th percentile (p90) ambulance ramping time (minutes) for the 15 hospitals with the highest overall p90 ramping. Hospitals are anonymised and ranked in descending order of their overall p90 ramping across the study period. Cell values and colour intensity represent the hourly p90 ramping duration, facilitating comparison of temporal congestion patterns among the most severely affected hospitals.
Fig. 12: Hourly 90th percentile (p90) ambulance ramping time (minutes) for the 15 hospitals with the highest ambulance case loads. Hospitals (H01-H15) are anonymised and ordered by decreasing ambulance case loads. Cell values and colour intensity represent the hourly p90 ramping duration, highlighting temporal variation in severe ramping across high demand hospitals.
Fig. 13: Median and the 90th percentile ramping duration stratified by hour of day (rows) and concurrency bin (columns), across the ED cohort. Concurrency is defined as the number of other ambulance arrivals at the same hospital in the 60 minutes preceding each case’s own arrival. Colour intensity reflects their ramping magnitude. The joint gradient, rising ramping from left to right within any row AND from top to bottom within any column, indicates that hospital arrival concurrency and hour of day each contribute independently to ramping outcomes; neither is an artifact of the other. Cells with fewer than 50 cases are masked in white.
Concurrency
Cases
Median
p75
p90
bin
( n )
(min)
(min)
(min)
0
388,362
23.6
38.0
66.1
1
406,521
28.4
47.6
82.4
2
365,989
31.6
53.7
91.8
3
302,032
34.1
58.1
97.4
4-5
380,616
36.6
62.4
103.5
TABLE VII: Relationship between prior ambulance arrivals within the preceding hour (hospital concurrency) and ramping duration.
Fig. 14: Median and 90th-percentile ramping durations by concurrency bin for the nine highest–case load ED hospitals. Each panel corresponds to a single hospital. Blue circles indicate median ramping; red triangles indicate the p90 value.
Ambulance response is time-critical in out-of-hospital cardiac arrest (OHCA), where dispatchers must balance timely arrivals with limited fleet capacity. Static territories and deterministic travel-time estimates are vulnerable to dynamic congestion, while always-dual dispatch adds redundancy but consumes fleet capacity. We propose IDEAL (Intelligent Dual dispatch of Emergency AmbuLances), a selective dual-dispatch framework that sends a second ambulance only when the optimistic gap between primary and secondary paths exceeds a threshold. IDEAL learns context-specific edge travel times from trip-level dispatch records, including unobserved routes, using a weakly supervised bilevel representation network. We train the nonsmooth model with mini-batch conservative gradients and prove an asymptotic convergence guarantee. IDEAL models uncertainty via Burg-divergence perturbations to a shared metric in the learned representation space, thereby inducing correlated changes in edge travel times and learning context-specific radii from historical underprediction errors. For real-time decisions, IDEAL casts optimistic-gap computation as a difference-of-convex program and derives an efficient oracle with complexity guarantees. In collaboration with the Hong Kong Fire Services Department, we evaluate IDEAL using historical OHCA records and real-time adaptive simulations. The results achieve a stronger response-time/resource trade-off relative to all region-based and Google-based baselines.
Zikun Lin, Daniel Zhuoyu Long, Viet Anh Nguyen
Department of Systems Engineering and Engineering Management, The Chinese University of Hong Kong
The ambulance service is the main transport for diseased or injured people which suffers the same acceleration forces as regular vehicles. These accelerations, caused by the movement of the vehicle, impact the performance of tasks executed by sanitary personnel, which can affect patient survival or recovery time. In this paper, we have trained, validated, and tested a system to assess driving in ambulance services. The proposed system is composed of a sensor node which measures the vehicle vibrations using an accelerometer. It also includes a GPS sensor, a battery, a display, and a speaker. When two possible routes reach the same destination point, the system compares the two routes based on previously classified data and calculates an index and a score. Thus, the index balances the possible routes in terms of time to reach the destination and the vibrations suffered in the patient cabin to recommend the route that minimises those vibrations. Three datasets are used to train, validate, and test the system. Based on an Artificial Neural network (ANN), the classification model is trained with tagged data classified as low, medium, and high vibrations, and 97% accuracy is achieved. Then, the obtained model is validated using data from three routes of another region. Finally, the system is tested in two new scenarios with two possible routes to reach the destination. The results indicate that the route with less vibration is preferred when there are low time differences (less than 6%) between the two possible routes. Nonetheless, with the current weighting factors, the shortest route is preferred when time differences between routes are higher than 20%, regardless of the higher vibrations in the shortest route.
Abdulaziz Aldegheishem, Nabil Alrajeh, Lorena Parra +2
Emergency Departments (EDs) are critical access points in healthcare systems, yet they face persistent pressure from unpredictable patient demand, seasonal surges, and non-urgent visits. Effective ED planning requires forecasts at multiple decision-making levels: hospitals need local demand estimates for staffing and bed management, regions require forecasts to coordinate healthcare units, and national authorities need system-wide projections for capacity planning. However, most existing approaches forecast ED demand independently at a single level, ignoring the hierarchy linking hospitals, regions, and national systems. This can produce incoherent predictions, where hospital-level forecasts do not aggregate consistently to regional or national demand. We propose HierSTT, a hierarchical Transformer-based framework for coherent multi-level ED forecasting. HierSTT jointly predicts hospital, regional, and national level demand in a single end-to-end model. A Temporal Fusion Transformer captures national dynamics, while spatio-temporal Transformer encoder-decoder modules model regional and hospital demand conditioned on higher-level forecasts. A coherence-aware loss penalizes cross-level inconsistencies during training. We further introduce a nationwide Portuguese ED dataset covering 81 hospitals across 5 regional health administrations, with heterogeneous covariates at each level. Experiments show that HierSTT reduces average WAPE by 32% relative to the best non-hierarchical deep learning baseline and outperforms all classical hierarchical reconciliation methods, while producing near-coherent predictions across levels. Additional resources associated with this work are available at https://github.com/FilipaLino/HierSTT.
Filipa Lino, Bárbara Tavares, Carlos Santiago +2
Institute for Systems and Robotics, LARSyS, Instituto Superior Técnico, Portugal · NOVA School of Science and Technology, Caparica, Portugal