1. Introduction
The first-come, first-served (FCFS) principle is a traditional approach to solving scheduling problems [
1]. In general, an FCFS-based scheduler computes the target departure or arrival times for flights in the order in which they were initially requested and schedules each time as a single point [
2,
3]. For example, each flight in the queue is assigned a target takeoff time (TTOT) based on the runway separation constraints, and the target off-block time (TOBT) can be calculated by subtracting a predetermined nominal taxi time for the gate–runway combination. FCFS is simple in the sense that the scheduling results can be readily understood and can be considered equitable, as the departure or arrival order does not change.
Two major issues with the traditional FCFS scheduler are: first, that it is difficult to accommodate multi-point constraints along the flight path, and second, that it tends to miss available time slots in order to preserve the original order. Consequently, if the demand becomes higher, flights further down the queue can accumulate excessive delays.
Most of the recent research on scheduling algorithms are based on optimization, such as mixed-integer linear programming (MILP) [
4,
5,
6,
7,
8,
9,
10]. Smeltink et al. [
4] presented a MILP-based formulation to solve the airport taxi scheduling problem. Visser and Roling [
5] investigated a similar MILP-based taxi scheduler with a multi-route option. In [
6,
7], an optimization model for taxiway and runway scheduling was studied using a MILP approach, and Lee and Balakrishnan [
8] proposed a MILP model to optimize both the taxiway and runway scheduling problems simultaneously. Bosson and Sun [
9] suggested stochastic scheduling optimization considering uncertainty in the MILP formulation for surface routing and scheduling. Yang and Gao [
10] also proposed a stochastic scheduling model based on MILP that can solve the ground movement problem (GMP) combining taxiway routing and gate allocation.
Schedulers using MILP-based optimizations usually provide optimal solutions, but they are computationally very costly and sometimes do not converge. So they are not suitable for large-scale problems or for operational situations for which a timely solution is required. To overcome these issues, researchers have been working in two main areas: one is improving the formulation [
11,
12,
13,
14,
15,
16,
17], and the other is improving the efficiency of the solver [
18,
19,
20,
21,
22,
23].
Avella et al. [
11] proposed time-indexed formulations that can be solved by a general-purpose mixed-integer programming (MIP)-based solver for the runway scheduling problem; they showed significant improvement in computation time. Prakash et al. [
12] also improved computation performance of a MIP-based algorithm to solve aircraft sequencing problems using data-splitting that enables parallel computing. Wang et al. [
14] introduced a chance-constrained programming-based optimization model combined with the quickest-path problem with time windows to solve the GMP while considering taxi time uncertainties. Zou et al. [
15] and Ornek et al. [
16] presented a way to formulate the problem in a pure-integer programming-based optimization model for taxi scheduling and the flight–gate assignment problem (FGAP). A multi-objective programming model was also investigated. Hu et al. [
13] considered uncertainty theory for FGAP, and Dönmez et al. [
17] used the results for terminal area management.
Ma et al. [
18] adjusted the simulated annealing algorithm by combining time decomposition to efficiently solve large-scale problems. Yan et al. [
19] improved a particle swarm optimization based on a receding horizon to solve the joint scheduling problem on the airport surface. Sun et al. [
20] used a genetic algorithm to improve the computational speed for solving the bilevel programming-based optimization problem; their model simultaneously considers gate assignment, taxi path scheduling, and pushback time delay. Xu et al. [
21] proposed a faster MIP solver that combined matrix approximation and an ant colony algorithm to solve aircraft arrival and departure scheduling problems. Zhang et al. [
22] showed that a adding cosine mutation and adaptive grouping to a tunicate swarm algorithm decreases the computation time when solving the gate assignment problem (GAP). Du et al. [
23] solved the robust GAP to minimize ground conflict using a column-generation-based heuristic method for computational efficiency.
However, continuous efforts have been made to improve the FCFS scheduler. Meyn [
24] was the first to propose a propagation-based multi-point scheduling methodology in a node–link structure; they established a framework to address two fundamental issues of the simple FCFS scheduler. In [
25], the capabilities of their scheduler were demonstrated by imposing multi-point flow-rate constraints at multiple metering fixes along the flight path. Palopo et al. [
26] introduced a capacity constraint represented by a maximum number of aircraft in a sector to study the traffic flow management problem in a dynamic airspace environment. Park (C) et al. [
27] combined the work of Meyn [
24] and Palopo et al. [
26] so that the scheduler can handle both flow-rate-style constraints at nodes such as metering fixes and maximum-capacity constraints at links such as sector monitor alert parameters (MAPs) to develop an advanced FCFS scheduler, and they investigated the effectiveness of sector transit time control to efficiently utilize available time slots.
Park (B) et al. [
28] proposed the extended FCFS (EFCFS) scheduler for surface operation; it can enforce minimum runway separation constraints based on the wake turbulence categories. In addition, in [
28], the capability to explore multiple routes was added. It was demonstrated that delays can be further reduced by going around a congested junction or taxiway even though this increases the overall taxi distances. Park (B) [
29] proposed a change to the original propagation architecture of Meyn [
24] to handle the link directionality constraint without using the negative aircraft count concept and to directly enforce the minimum along-trail separation between leading and trailing aircraft instead of using an indirect link capacity. In [
29], iterative rescheduling, which takes advantage of lower computational cost, was investigated to demonstrate that the scheduler can be used in an actual operational environment. Phillips and Sadovsky [
30] enhanced [
24] by enforcing a finite time interval to handle uncertainty and to be used for distributed scheduling.
Park (B) et al. [
31] performed a study to compare the EFCFS with a MILP-based scheduler of Eun et al. [
32]. First, each constraint was investigated to ensure the compatibility of the two schedulers, and then, a large number of scheduling problems were solved at Incheon International Airport (ICN). The results show that the schedules computed by EFCFS have, on average, about 18% larger delay. However, computation time was at least one order of magnitude smaller for the EFCFS.
In this paper, the EFCFS scheduler using the architecture of [
29] is employed to validate the scheduling results with actual historical data. The scenarios were generated based on the actual schedules and aircraft types from the FOIS data and the taxi routes extracted from the airport surface detection equipment (ASDE) track data at ICN. Differences in takeoff times (TOTs) and landing times (LDTs) calculated by comparing the scheduling results with the initial plan information from FOIS were used as metrics. The results show that the average takeoff delay can be reduced by as much as about 19 min depending on the gate occupancy constraint. It is also confirmed that the scheduler effectively utilizes available time slots by switching the departure or landing order for a small number of flights.
The main contribution is that this work has established a systematic framework of utilizing the operational data (FOIS and ASDE data) to validate the performance of the EFCFS scheduler, whereas most of the previous work focused on developing the scheduling algorithm itself. The framework includes the comparison metrics, the generation of input schedules, and the construction of complicated constraints. In addition, the framework is demonstrated on a complex airport with a large number of flights.
After this introduction,
Section 2 describes the EFCFS scheduler and its constraints on the nodes and links.
Section 3 explains how the historical FOIS data are used to generate the reference schedule.
Section 4 explain the scenario, and
Section 5 shows the scheduling results.
Section 6 concludes this paper.
3. Schedule Generation
To generate realistic scenarios, this paper uses FOIS data [
33], which are typical data that provide flight schedules, and ASDE data. FOIS provides the basic flight information, such as the callsign, aircraft type, origin airport and gate, and destination airport and gate. However, FOIS does not contain runway information. FOIS is divided into two categories: departure and arrival. The departure FOIS contains the estimated time of departure (ETD) and actual time of departure (ATD) at the origin airport and the estimated time of arrival (ETA) at the destination airport. Comparing the FOIS times with the ASDE tracks confirms that ETD and ATD are the times at the departure runway. The arrival FOIS provides the ETA and the actual time of arrival (ATA) at the destination airport and the ETD at the origin airport. Similar to the departure FOIS, ETA and ATA are scheduled at the arrival runway.
In addition to crosschecking the FOIS times, the ASDE data are used to extract the actual routes on the surface. ASDE includes dense track data of the flights and all other objects that moved on the surface. The surface route of each flight is extracted by comparing the track positions with the surface node–link model.
Among the 500 flights used in this study, seven flights in FOIS do not have usable track data. For those flights, the most commonly used departure or arrival runways are assigned and the minimum distance route between the runway and the gate is used as a substitute.
The original planned flight schedules are generated by combining flight information, ETA, ETD, and extracted routes. Since the ASDE track data are the result of traffic control and scheduling, ATA, ATD, and tracking times are not appropriate for validating the performance of the presented scheduler. Consequently, in this paper, the ETDs and ETAs were used to generate the original plan. For departure, the ETD is the planned runway takeoff time, STOT, so the scheduled OBT (SOBT) is calculated by subtracting the scheduled taxi-out time (SXOT) from the ETD. For arrivals, the ETA is the planned landing time, SLDT, at the runway. The scheduled IBT (SIBT) for arrival is the sum of SLDT and the scheduled taxi-in time (SXIT). Scheduled taxi times—SXOT and SXIT—are the sum of all nominal link transit times in the route from the runway threshold to the gate. The arrival route includes runway links from touchdown to the runway exit; however, the departure route ends at the runway threshold and does not contain any runway links. The nominal transit time of each surface link is calculated by dividing the length of the link by the nominal taxi speeds, which are set to 5, 10, and 15 knots for the gate, ramp, and taxiway links, respectively. A speed of 150 knots is assumed for the runway links.
Table 1 is a sample plan generated using the departure track in
Figure 6a that shows the arrival route of Flight1 from runway 16R to gate G21.
Figure 6b shows the departure route of Flight2 from gate G14 to the takeoff from runway 16L. The nominal transit times of each link are the differences between the entry times and exit times.
4. Problem Setup
To evaluate the EFCFS scheduler, surface scheduling for departures and arrivals was tested at ICN.
Figure 7 shows the ICN node–link model, which has 1006 links and 775 nodes with 262 gates and 8 runway thresholds for 2 pairs of parallel runways. The node–link model was created manually using the aeronautical information publication document and satellite images from Google Earth.
The originally planned schedule was generated using FOIS and ASDE data at ICN on 15 August 2022. The schedule includes 250 departures and 250 arrivals, with 150 medium, 343 heavy, and 7 super-heavy flights.
Table 2 shows the total fleet mix of the generated schedule.
Table 3 and
Table 4 show the node constraints used for EFCFS computation. No additional occupancy constraints are applied to the gate nodes because the SIBTs and SOBTs of the input schedules already reflect the gate occupancy, and the scheduler does not allow TOBT earlier than SOBT. The ramp and taxi nodes are constrained by a ten-second blocking time. The runway nodes are constrained by a 30-second blocking time and runway occupancy times.
The runway separation criteria are considered as additional constraints for the runway threshold nodes. In this paper, modified runway separations based on the separation standard were used [
34].
Table 5 shows the runway separations for single runways that are applied to the scheduler. The separations between departure after arrival and arrival after departure are set to one minute. The separations between departures and between arrivals on adjacent runways are set to two minutes. For the opposite direction, a two-minute separation is applied between departures and between arrivals for both the single runway and adjacent runway operations. All other cases are set to zero separation.
Table 6 shows the total runway separations for ICN that were applied to the scheduler.
To ensure the runway separation constraints are aligned with the actual operation, historical runway separations are investigated using the ASDE data.
Figure 8 shows the distributions of the actual runway separations between departures and between arrivals on runways 15R/33L and 15L/33R from the one-month ASDE data, which confirms that the constraints given in
Table 5 are reasonable.
The link blocking time is set to ten seconds for all types of links to enforce a minimum ten-second along-trail separation between the taxiing aircraft. The taxi speed can be reduced by up to 10%, except on the runway, which increases the link transit time by up to 10% of the nominal transit time.
5. Scheduling Results
The EFCFS scheduler provides different results depending on the priorities. Among several prioritization strategies, in this study, flights are scheduled based on the nominal priority, which is the order of the planned SOBT and SLDT.
Delays are measured by computing the difference between the scheduled time and the target time, which is calculated by the scheduler.
Figure 9 shows the distributions of four delays—DOBT, DIBT, DTOT, and DLDT—which are defined in Equations (
1)–(
4).
Table 7 summarizes the average and maximum delay times. Most delays are between zero and one minute, and the average DOBT, DIBT, DTOT, and DLDT are 0.79, 0.18, 1.01, and 0.17 min, respectively. The maximum DOBT, DIBT, DTOT, and DLDT are 12.43, 6.9, 13, and 6.9 min, respectively. The average DTOT and DIBT are larger than the average DOBT and DLDT. This shows that the scheduler reduced the taxi speed of the flights to better utilize available time slots. Mostly, the departure flights are adjusted while the arrivals are not significantly affected by the scheduler.
Since the speed of 150 knots is used for all runway links, the runway occupancy time will be slightly smaller, especially for arrivals. The runway link speeds are reduced to 100 and 50 knots to test the impact, and the results do not show any significant differences. Therefore, the original value of 150 knots is used for all subsequent runs.
Figure 10 shows the runway throughputs for the scheduling results.
Figure 10a shows the departure runway throughputs, and
Figure 10b shows the arrival runway throughputs. The maximum throughput of each runway is 30 flights per hour, and there is no runway separation constraint violations in the scheduling results. Similarly, no other constraint violation is detected in the scheduler output.
Figure 11a shows a portion of the STOTs and TTOTs. Solid lines denote zero delay, and dashed line denotes flights with a TOT adjusted by the scheduler. The square markers show that the flights with the same planned time have been adjusted for proper spacing. Flights marked with a red x in the TTOT column represent an order change. It can be observed that the scheduler took advantage of the available time slots by changing the order of several flights. Similar results can be observed for landing flights, as shown in
Figure 11b.
5.1. Comparison with the Operational Data
The scheduling results are compared to the actual FOIS data to validate performance of the scheduler. ATD and ATA are equal to the actual TOT (ATOT) and LDT (ALDT), respectively. To compare the results with the scheduler output, the difference between ATOT and TTOT is used for departures, and the difference between ALDT and TLDT is used for arrivals. If the difference is positive, it can be judged that the scheduler provides a better time than the actual data.
Figure 12 show the time difference distributions between the scheduling result and the actual data. In
Figure 12a, the differences in TOT are skewed to the positive side. On the contrary, the majority of the LDT differences are mostly zero and are located on the slightly negative side in
Figure 12b.
Table 8 shows the average and maximum differences. The average differences for TOT and LDT are 18.76 and −0.65 min, respectively. The maximum differences for TOT and LDT are 103.85 and 41 min, and the minimum differences for OBT and LDT are −10 and −28 min, respectively. For this particular example, the EFCFS significantly reduced the departure delays while having negligible impact on the landing delays.
5.2. Gate Occupancy Implementation
The comparison results without considering gate occupancy show that the EFCFS scheduler is effective at reducing departure delays. Since the delay could have been reduced by the unconstrained gate nodes, gate occupancy times are investigated for the flights that arrive and depart ICN on the same day. Among 500 flights, 143 flights are identified.
Figure 13 shows the distributions of the intervals between ATA and ATD of those flights. As these time intervals are between landing and takeoff, actual gate occupancy times are smaller. It can be observed that for the most of the flights, the time interval is between 1.5 and 2.5 h, which is consistent with the typical gate occupancy times for international flights.
To study the impact of gate occupancy constraints, gates are blocked for
starting from the TIBT for inbound flights, as shown in
Figure 14, and all the flights are rescheduled. The values for
is set to one or two hours based on the trend shown in
Figure 13, and the outcomes are compared to the previous results without the gate occupancy constraints in order to evaluate the impacts.
Table 9 shows the comparison of the scheduling results by the gate constraints. There was no drastic change for arrivals, but departure delays increase as the gate occupancy time increases. For departures, the averages for DOBT and DTOT increase from 0.79 and 1.01 to 2.41 and 2.64 min, respectively, when the gate is occupied by arrivals for one hour. And the average DOBT and DTOT significantly increase to 8.89 and 9.11 min when the gate occupancy time increases to two hours. The standard deviations also increase as the gate occupancy time increases, similar to the average delays.
Table 10 summarizes the comparisons of the scheduling results and the actual data for the gate constraints. The scheduler produced results similar to the actual data as the gate occupancy time increased.
Figure 15 shows the takeoff and landing time difference distributions between the scheduler and the actual data for the gate constraints. The distributions slide from positive to negative and also become wider as the occupancy time increases.