# An Analysis of Robustness Approaches for the Airport Baggage Sorting Station Assignment Problem

Authors: Amadeo Ascó
Journal: Journal of Optimization
 (2016)

---

## Abstract

In the Airport Baggage Sorting Station Assignment Problem (ABSSAP), the Baggage Sorting Stations (BSSs) are assigned to flights for the period of time necessary to perform their service for a given flights’ schedule. But the flights schedule may change on the day of operation which may deem the original assignment of some flights to BSSs infeasible. These changes may create conflicts between those flights whose schedules have changed and may not be restricted to those flights but propagating to the other flights for different reasons. Conflicts depend on the original assignments for the real arrival and departure flight times on the day of operation. It is therefore desirable to consider potential delays on the day of operation when generating the original flight assignments to BSSs, such that the final flight assignments differ little or do not differ at all from the original assignments on the day of operation. The term robustness is here used to give an indication of the degree to which this has been achieved. Some existing approaches originally presented in the Airport Gate Assignment Problem (AGAP) are adapted to the ABSSAP, other approaches are suggested for generating assignments which take account of potential perturbations on the day of operation for the ABSSAP, and all of them are then compared. It is shown that the suggested approaches by themselves do not perform better than the other considered approaches but when combined they enhance the result further compared to when each approach is used alone.


---

## Body

## 1. Introduction

Flight delays are caused by many factors like airport security, weather conditions, unavailability of required resources (mechanical breakdown), delayed propagation, airport congestion, etc.The ABSSAP corresponds to the assignments of flights to BSSs, whereas the BSSs are the end elements in the baggage system at the ground of the airport terminal where the baggage is collected from and sorted by the baggage handlers or to where the baggage from arrived flights are placed to enter the baggage system. Figure1 shows an overall simple view of the airport baggage handling system.Figure 1
Simple view of an airport baggage handling system with the BSSs.Once the baggage has been loaded into the aircraft hold, the BSSs assigned to that flight are ready for use on the next scheduled flights, and any delay on a flight at the gate will not necessarily have repercussions on these. However, such delays may affect the BSSs assigned to the flight scheduled for that gate next if such a flight is held waiting for the gate to become free. BSSs will then need to hold the baggage longer, thus potentially affecting the following assignment in turn.Delays on the day of operation may render some assignments infeasible which need to be reassigned. It is therefore desirable to account for potential delays on the day of operation when generating the flight assignments to BSSs in the planning stage, such that the final flight assignments differ little or do not differ at all from the original assignments on the day of operation. The degree to which this is achieved is an indication of the solution robustness, so a solution which requires fewer reassignments is said to be more robust than those solutions requiring more reassignments. Robustness is the ability of assignments to resist changes consequence of perturbations by reducing or removing the need to reassign current assignments. Mulvey et al. [1] characterise the desirable properties of a solution to model robustness.There are different ways of increasing robustness depending on the intended effect. One of the most simple and widely used methods is the introduction of a buffer time between assignments which allows the absorption of small disturbances; implementation of this approach for the ABSSAP is presented in [2–5]. The amount of buffer time may take different contributory factors into account, which could perturb the schedules of handlers, airline, airport origin, destination, and flight. The size of the flight is normally related to the travel distance, longer distances presenting a higher probability of disturbances, which may accumulate generating a higher level of delays. Shorter distances present less chance of disturbances. The location of the origin and destination airports has a direct effect on potential disruptions, given that they place a constraint on the permitted routes possible and certain circumstances applicable to them, such as weather patterns during summer and winter or on the equator. Whereas the buffer time approach considers a predefined period of time, where its reduction is considered as a measure of robustness, the “idle time” refers to the time between two consecutive assignments to the same BSS, from the end time of one activity to the base starting time of the following assignment, also called gap. Other approaches for improving the robustness make use of the distribution of “idle time” and the reduction of the number of reassignments of the disrupted schedules.Buffer time has been used in the scheduling of BSSs and airport gates, which may be applied between two consecutive flights to the same BSS in order to absorb small disturbances in the real system behaviour and was studied in [1, 6]. It was used in the ABSSAP in [2, 3, 5, 7] where it was used in the optimization of the aircraft turnaround process. The AGAP has some characteristics similar to the ABSSAP, and buffer times have been commonly considered for the AGAP in [7–11]. Yan and Huo [12] provided a sensitivity analysis for the AGAP buffer time, noting that the length of buffer time significantly influences the gate assignment process, so a reasonable minimum value should be used. Yan et al. [8] looked at the suitability of Flexible Buffer Times (FBTs) where, given low delays, short FBTs usually improve real-time objectives, such as the reassigning of an incoming aircraft at a minimum distance. Wei and Liu [13] showed the feasibility and effectiveness of using a fuzzy model in conjunction with fixed buffer times for the AGAPs. Ascó et al. [4] used buffer times to cope with small perturbations in the ABSSAP, and several constructive algorithms were also studied. Wu and Caves [7, 14] showed the significance of a correct use of scheduled buffer time in maintaining schedule punctuality and performance by balancing the trade-offs between schedule punctuality and aircraft utilisation.


## 2. Problem under Consideration

The ABSSAP considered in this paper may be summarised as the assignment of available BSSs to flights which have already been scheduled. In the ABSSAP, the flights will already have been assigned to stands, which are often grouped along piers around the terminals, and there will usually be some bias in this allocation, such as to airline preferences. This problem was described and studied in [2–5]. In this paper, we concentrate on finding solutions that comply with the constraints and take into account the objectives concentrating on the performance of different robustness approaches described in Section 3.


### 2.1. Overall Problem

Aircraft are usually parked at their allocated stand, around an airport terminal. Two layouts are considered in this paper which were presented in [2].A list of some of the constants for this problem is presented in Table1 and a list of the decision variables is shown in Table 2.Table 1
The constants for the model.
Constant
Description
N
Denotes the number of BSSs.
M
Denotes the number of flights.
P
j
The total number of activities to be serviced by baggage service stations for a given flightj, which also equates to the total number of sorting stations required to fully service flight j, P
j > 0.
i
Denotes a BSS where BSS 0 represents the dummy sorting station to which flights are assigned, if they cannot be assigned to a real BSS,i
∈
{
0
,
…
,
N
}.
j
Denotes a flight,j
∈
{
1
,
…
,
M
}.
T
j
Denotes the required service time for flightj (1 hour for short haul and 1.33 hours for long haul).
B
j
Denotes the desired buffer time for flightj—the time for which its BSS should be idle prior to this flight being serviced. B
j is considered to be 15 minutes for short haul and 30 minutes for long haul flights.
e
j
Denotes the end service time for flightj.
t
j
Denotes the target starting service time for flightj (t
j
=
e
j
-
T
j
-
B
j).
τ
j
Denotes the the base starting service time for flightj (τ
j
=
e
j
-
T
j).Table 2
List of the decision variables used in this ABSSAPs’ model.
Name
Description
y
i
j
p
Specifies the assignment of flights to sorting stations.y
i
j
p
=
1 if baggage sorting station i
∈
[
1
,
…
,
N
] is allocated to flight j
∈
[
1
,
…
,
M
] for p
∈
[
1
,
…
,
P
j
] and 0 otherwise. If each flight only requires one activity, which means that each flight only requires one BSS, then this variable can be expressed as y
i
j.
r
j
p
Specifies the necessary reduction in service time for activityp
∈
[
1
,
…
,
P
j
] of flight j
∈
[
1
,
…
,
M
], given the service starting time allocated, s
j
p.
s
j
p
The service starting time allocated to activityp
∈
[
1
,
…
,
P
j
] of flight j
∈
[
1
,
…
,
M
] and given that a sorting station can only service one flight at a time. s
j
p can be determined from r
j
p since s
j
p
=
t
j
-
r
j
p.The aim is to findy
i
j
p values such that the objective function (3) is maximised, subject to the constraints expressed by (1) and (2) and the different robustness approaches studied in this paper which are presented in Section 3.(1)
∑
i
=
0
N
y
i
j
p
=
P
j
∀
j
∈
1
,
…
,
M
,
∀
p
∈
1
,
…
,
P
j
.For any pair of different flights where service times overlap, if the overlap in service times is greater than the maximum reduction allowed (B
l
q for activity q of flight l), then both flight activities cannot be assigned to the same BSS. Thus, inequality (2) applies to any such pair of flights, j and l (j
≠
l), where t
l
q
<
e
j
≤
e
l and (
e
j
-
t
l
q
)
>
B
l
q. (2)
y
i
j
p
+
y
i
l
p
≤
1
.Equation (1) states that each flight can be assigned either to exactly one real BSS or to the dummy sorting station. Inequality (2) states that flights cannot be assigned to the same real BSS if their service times overlap. If all flights only require one BSS, then the variables are reduced to y
i
j, r
j, and s
j.Although this is inherently a multiobjective problem, the importance of ensuring maximal assignment of flights to BSSs (top priority) and the relative importance of obtained robust solutions (second priority) allow these objectives to be combined into a single compound objective (3) with weights W
1, W
2, and W
3 chosen to implement these priorities. (3)
W
1
∗
∑
i
=
1
N
∑
j
=
1
M
∑
p
=
1
P
j
y
i
j
p
︸
assignment  objective
-
W
2
∗
Robustness
︸
robustness  objective
-
W
3
∗
∑
j
=
1
M
∑
p
=
1
P
j
C
j
·
∑
i
=
1
N
y
i
j
p
·
d
i
j
︸
distance  objective
.The first element in (3) aims to maximise the assignment of flights to BSSs, the second aims to increase the robustness, and the third aims to minimise some distance cost associated with the assignments, where C
j is a factor related to the amount of baggage for flight j (assumed to be 1 in all cases for this paper) and d
i
j is a measure of the distance or cost incurred from assigning BSS i to flight j. This aims to ensure that flights are allocated to appropriate BSSs.The following two points were originally defined in [4] and will be observed to be useful later when defining the new suggested robustness approaches and when interpreting the results for the ABSSAP.TheLower Maximum Assignment Point (LMAP) is the minimum number of resources required to service a certain number of activities when the service starting time (s
j
p) coincides with the base starting service time (τ
j).TheUpper Maximum Assignment Point (UMAP) is the minimum number of resources required to service a certain number of activities when the service starting time (s
j
p) coincides with the target starting service time (t
j
p).Two examples of both points are shown in Figure2, where the maximum values of each diagram correspond to the LMAP and UMAP.Figure 2
Flights service distribution for the two days’ data sets obtained from the website of British Airports Authority (BAA).
(a)
Representation without reduction in service time, Upper Maximum Assignment Point (UMAP)
(b)
Representation with maximum reduction in service time, Lower Maximum Assignment Point (LMAP)


### 2.2. Robustness

A scheduled assignment is said to be in conflict if the completion of its service time is greater than the commencement of the service time of the next assignment at the same BSS. When a delayed flight affects the assignment of subsequent flights to the BSS, then there are two ways it can be corrected: either reassign the conflicting flight or reassign the subsequent flight to the conflicting one. When reassigning a conflicting flight or subsequent flights to another BSS, a situation may arise where the reassigned flight is in conflict with the subsequent flight at the new BSS. Some reassignment may therefore have a downstream effect on the overall schedule, producing further conflicting flights requiring further reassignments, thus potentially increasing the problem difficulty later on.In the ABSSAP, those flights which are late in arriving at their assigned stand are considered to be a perturbation, since baggage cannot be loaded into the aircraft at the scheduled time, and it needs to be held longer at the BSS. Any extension of an aircraft’s stay on its assigned gate should not have an effect on its assigned BSS, as the baggage should have already been loaded into the aircraft, such that the BSS is free for use in its next assignment. This means that not all aircraft delays will affect their assigned BSS.The main objective of this is to reduce the number of BSSs which have to be reassigned on the day the schedule is put into practice. In the ABSSAP,e
j is the end of the service time of flight j, τ
j is the base starting service time of flight j, and the variable y
i
j has a value of 1 if flight j is assigned to BSS i, or y
i
j is zero otherwise, as described in Section 2.1. A new decision variable x
j
k is introduced with a value of 1 if flights j and k are assigned to the same BSS (i.e., if y
i
j
=
y
i
k
=
1 for i
∈
[
1
,
…
,
N
] and j
,
k
∈
[
1
,
…
,
M
]), or x
j
k is 0 otherwise. On the day the schedule is implemented, the real times for the flights j and k correspond to e
j
′, τ
j
′, e
k
′, and τ
k
′ for flights j and k, respectively, where τ
j
′
≥
τ
j and τ
k
′
≥
τ
k, and the flights ordered by their base starting times, for j
<
k and j, k
∈
[
1
,
…
,
M
]; then τ
j
≤
τ
k. A conflict occurs when two flights j and k with j
<
k and j, k
∈
[
1
,
…
,
M
] originally assigned to the same BSS (x
j
k
=
1) have overlapping service times; that is, e
j
′
>
τ
k
′ and τ
j
′
<
e
k
′, as shown in Figure 3. Conflicts depend on the original assignments for the real arrival and departure flight times on the day of operation. It is therefore desirable to consider potential delays on the day of operation when generating the original flight assignments to BSSs, such that the final flight assignments differ little or do not differ at all from the original assignments on the day of operation. B
j refers to the buffer time originally assigned to the flight j, which may be absorbed by perturbations without raising the need for a reassignment. But overlapping, τ
k
′
<
e
,
′, will deem the assignment infeasible and require reassigning it.Figure 3
Example of conflict between two flights originally assigned to the same BSS.


## 2.1. Overall Problem

Aircraft are usually parked at their allocated stand, around an airport terminal. Two layouts are considered in this paper which were presented in [2].A list of some of the constants for this problem is presented in Table1 and a list of the decision variables is shown in Table 2.Table 1
The constants for the model.
Constant
Description
N
Denotes the number of BSSs.
M
Denotes the number of flights.
P
j
The total number of activities to be serviced by baggage service stations for a given flightj, which also equates to the total number of sorting stations required to fully service flight j, P
j > 0.
i
Denotes a BSS where BSS 0 represents the dummy sorting station to which flights are assigned, if they cannot be assigned to a real BSS,i
∈
{
0
,
…
,
N
}.
j
Denotes a flight,j
∈
{
1
,
…
,
M
}.
T
j
Denotes the required service time for flightj (1 hour for short haul and 1.33 hours for long haul).
B
j
Denotes the desired buffer time for flightj—the time for which its BSS should be idle prior to this flight being serviced. B
j is considered to be 15 minutes for short haul and 30 minutes for long haul flights.
e
j
Denotes the end service time for flightj.
t
j
Denotes the target starting service time for flightj (t
j
=
e
j
-
T
j
-
B
j).
τ
j
Denotes the the base starting service time for flightj (τ
j
=
e
j
-
T
j).Table 2
List of the decision variables used in this ABSSAPs’ model.
Name
Description
y
i
j
p
Specifies the assignment of flights to sorting stations.y
i
j
p
=
1 if baggage sorting station i
∈
[
1
,
…
,
N
] is allocated to flight j
∈
[
1
,
…
,
M
] for p
∈
[
1
,
…
,
P
j
] and 0 otherwise. If each flight only requires one activity, which means that each flight only requires one BSS, then this variable can be expressed as y
i
j.
r
j
p
Specifies the necessary reduction in service time for activityp
∈
[
1
,
…
,
P
j
] of flight j
∈
[
1
,
…
,
M
], given the service starting time allocated, s
j
p.
s
j
p
The service starting time allocated to activityp
∈
[
1
,
…
,
P
j
] of flight j
∈
[
1
,
…
,
M
] and given that a sorting station can only service one flight at a time. s
j
p can be determined from r
j
p since s
j
p
=
t
j
-
r
j
p.The aim is to findy
i
j
p values such that the objective function (3) is maximised, subject to the constraints expressed by (1) and (2) and the different robustness approaches studied in this paper which are presented in Section 3.(1)
∑
i
=
0
N
y
i
j
p
=
P
j
∀
j
∈
1
,
…
,
M
,
∀
p
∈
1
,
…
,
P
j
.For any pair of different flights where service times overlap, if the overlap in service times is greater than the maximum reduction allowed (B
l
q for activity q of flight l), then both flight activities cannot be assigned to the same BSS. Thus, inequality (2) applies to any such pair of flights, j and l (j
≠
l), where t
l
q
<
e
j
≤
e
l and (
e
j
-
t
l
q
)
>
B
l
q. (2)
y
i
j
p
+
y
i
l
p
≤
1
.Equation (1) states that each flight can be assigned either to exactly one real BSS or to the dummy sorting station. Inequality (2) states that flights cannot be assigned to the same real BSS if their service times overlap. If all flights only require one BSS, then the variables are reduced to y
i
j, r
j, and s
j.Although this is inherently a multiobjective problem, the importance of ensuring maximal assignment of flights to BSSs (top priority) and the relative importance of obtained robust solutions (second priority) allow these objectives to be combined into a single compound objective (3) with weights W
1, W
2, and W
3 chosen to implement these priorities. (3)
W
1
∗
∑
i
=
1
N
∑
j
=
1
M
∑
p
=
1
P
j
y
i
j
p
︸
assignment  objective
-
W
2
∗
Robustness
︸
robustness  objective
-
W
3
∗
∑
j
=
1
M
∑
p
=
1
P
j
C
j
·
∑
i
=
1
N
y
i
j
p
·
d
i
j
︸
distance  objective
.The first element in (3) aims to maximise the assignment of flights to BSSs, the second aims to increase the robustness, and the third aims to minimise some distance cost associated with the assignments, where C
j is a factor related to the amount of baggage for flight j (assumed to be 1 in all cases for this paper) and d
i
j is a measure of the distance or cost incurred from assigning BSS i to flight j. This aims to ensure that flights are allocated to appropriate BSSs.The following two points were originally defined in [4] and will be observed to be useful later when defining the new suggested robustness approaches and when interpreting the results for the ABSSAP.TheLower Maximum Assignment Point (LMAP) is the minimum number of resources required to service a certain number of activities when the service starting time (s
j
p) coincides with the base starting service time (τ
j).TheUpper Maximum Assignment Point (UMAP) is the minimum number of resources required to service a certain number of activities when the service starting time (s
j
p) coincides with the target starting service time (t
j
p).Two examples of both points are shown in Figure2, where the maximum values of each diagram correspond to the LMAP and UMAP.Figure 2
Flights service distribution for the two days’ data sets obtained from the website of British Airports Authority (BAA).
(a)
Representation without reduction in service time, Upper Maximum Assignment Point (UMAP)
(b)
Representation with maximum reduction in service time, Lower Maximum Assignment Point (LMAP)


## 2.2. Robustness

A scheduled assignment is said to be in conflict if the completion of its service time is greater than the commencement of the service time of the next assignment at the same BSS. When a delayed flight affects the assignment of subsequent flights to the BSS, then there are two ways it can be corrected: either reassign the conflicting flight or reassign the subsequent flight to the conflicting one. When reassigning a conflicting flight or subsequent flights to another BSS, a situation may arise where the reassigned flight is in conflict with the subsequent flight at the new BSS. Some reassignment may therefore have a downstream effect on the overall schedule, producing further conflicting flights requiring further reassignments, thus potentially increasing the problem difficulty later on.In the ABSSAP, those flights which are late in arriving at their assigned stand are considered to be a perturbation, since baggage cannot be loaded into the aircraft at the scheduled time, and it needs to be held longer at the BSS. Any extension of an aircraft’s stay on its assigned gate should not have an effect on its assigned BSS, as the baggage should have already been loaded into the aircraft, such that the BSS is free for use in its next assignment. This means that not all aircraft delays will affect their assigned BSS.The main objective of this is to reduce the number of BSSs which have to be reassigned on the day the schedule is put into practice. In the ABSSAP,e
j is the end of the service time of flight j, τ
j is the base starting service time of flight j, and the variable y
i
j has a value of 1 if flight j is assigned to BSS i, or y
i
j is zero otherwise, as described in Section 2.1. A new decision variable x
j
k is introduced with a value of 1 if flights j and k are assigned to the same BSS (i.e., if y
i
j
=
y
i
k
=
1 for i
∈
[
1
,
…
,
N
] and j
,
k
∈
[
1
,
…
,
M
]), or x
j
k is 0 otherwise. On the day the schedule is implemented, the real times for the flights j and k correspond to e
j
′, τ
j
′, e
k
′, and τ
k
′ for flights j and k, respectively, where τ
j
′
≥
τ
j and τ
k
′
≥
τ
k, and the flights ordered by their base starting times, for j
<
k and j, k
∈
[
1
,
…
,
M
]; then τ
j
≤
τ
k. A conflict occurs when two flights j and k with j
<
k and j, k
∈
[
1
,
…
,
M
] originally assigned to the same BSS (x
j
k
=
1) have overlapping service times; that is, e
j
′
>
τ
k
′ and τ
j
′
<
e
k
′, as shown in Figure 3. Conflicts depend on the original assignments for the real arrival and departure flight times on the day of operation. It is therefore desirable to consider potential delays on the day of operation when generating the original flight assignments to BSSs, such that the final flight assignments differ little or do not differ at all from the original assignments on the day of operation. B
j refers to the buffer time originally assigned to the flight j, which may be absorbed by perturbations without raising the need for a reassignment. But overlapping, τ
k
′
<
e
,
′, will deem the assignment infeasible and require reassigning it.Figure 3
Example of conflict between two flights originally assigned to the same BSS.


## 3. Robustness Approaches

In this section, some robustness approaches are presented and others are suggested. These are then studied and compared in the following section.


### 3.1. Minimise Reduction in Service Time

It may be possible to gain robustness by reordering assignments between BSSs so that “idle time” between flights consecutively assigned to the same BSS is greater, as shown in Figure4.Figure 4
Simple example of two schedules, with the same flights, where one is obviously more robust in respect of perturbations than the other.Figure4 shows two potential solutions with different robustness. While any delay to flight “a” in the “less robust” solution will certainly affect flight “b,” which will in turn have to be reassigned to another BSS, in the “more robust” solution a delay in flight “a” will not affect flight “b.” So the “more robust” solution is preferable to that of the “less robust” solution.Mangoubi and Mathaisel [15] proposed the use of “buffer times” between two flights which are consecutively assigned to the same gate in order to obtain robust assignments, defining the reduction in service as that part of buffer time which overlaps with the previous assignment to the same gate. Given the detrimental effects that the reduction in service time has on the robustness of the assignment as against real-life delays, it is advisable to minimise the total reduction in service time, thus maximising buffer times. This objective, total reduction in service time (TRS), can be expressed by formula (4).(4)
min
⁡
∑
j
=
1
M
∑
p
=
1
P
j
r
j
p
.If a fitness function is defined as a weighted sum of the different objectives, as used here, and the robustness uses the reduction in service time only for assigned flights, then particular weights may compromise the importance of the main objective (maximisation of the number of assignments), where solutions with lower numbers of assignments are favoured over those with higher numbers of assignments because of the robustness objective.The use of formula (4) treats any reduction in service equally, so it does not make a distinction between reducing all the assignment buffer time and allocating it to another flight from another solution where both flights share the available “idle time,” as shown in Figure 5. However, flight 3 with maximum reduction in service in solution “b” (Figure 5) will be unable to absorb any delay on the day of operation, although the same flight in the alternative solution “a” (Figure 5) will be able to do so, making it the preferable choice.Figure 5
Simple example of a more robust schedule of four flights and two sorting stations using the reduction in buffer time.A nonlinear penalty function is required to take account of this, such as formulas (5) and (6).(5)
u
i
j
p
=
arctan
⁡
τ
j
-
e
k
B
j
-
Π
2
i
f
k
<
j
,
y
i
j
p
=
y
i
k
q
=
1
,
∑
l
=
k
+
1
j
∑
p
=
1
P
l
y
i
l
p
=
1
,
0
o
t
h
e
r
w
i
s
e
,
(6)
max
⁡
∑
i
=
1
N
∑
j
=
1
M
-
1
∑
p
=
1
P
j
u
i
j
p
.Formula (5) defines the penalty u
i
j
p for assigning activity p of flight j to sorting station i, with the total penalty being represented by formula (6). There is only a penalty for consecutive assignments to the same sorting station, between the assignment and the previous assignment to the same sorting station. The reason for using arctangent is based on the properties of this function which imposes stronger penalties around the point at which the flights are assigned service time without any buffer time to the point at which all the buffer time is retained (by dividing the member of arctangent by B
k). The effect is reduced as the time separation between the assignments increases, but the contribution of any time separation is always considered negative. The constant π
/
2 is used so that there is always a cost associated with assigning two flights to the same sorting station; otherwise it would be a benefit. Considering the arctangent for each flight increases the individual effect on the overall objective, in contrast to using the arctangent of the sum of all of the flight contributions.The robustness objective represented by formula (6) also takes account of the objective of “minimising the service reduction,” as being the sum of the reduction in buffer time of each flight assigned.


### 3.2. Distribute Idle Time

Bolat [16] proposed the distribution of “idle time” uniformly amongst gates for the AGAP. In the case of the ABSSAP, the distribution of the “idle time” uniformly amongst BSSs can be considered, where “idle time” is calculated as the time between the start of service time on a flight’s assignment to BSS and the end of service time of the flight assigned immediately before the same baggage sorting station, as shown in formula (7). Whereas the buffer time implies preference for a particular gap size between consecutive assignments to the same sorting station, an “idle time” does not. When using the “idle time” it is normally intended that the gap size should be as large as possible. The reason for this is to increase the probability that, even with a delay, completion of service time for flight presently assigned will still be earlier than the start of service time for the next flight assigned to the same BSS. (7)
Δ
t
j
=
s
j
-
e
k
for
j
<
k
,
y
i
j
p
=
y
i
k
q
=
1
where
∄
l
∈
k
,
…
,
M
with
τ
l
<
e
k
;
j
,
k
∈
1
,
…
,
M
,
p
∈
1
,
…
,
P
j
,
q
∈
1
,
…
,
P
k
,
i
∈
1
,
…
,
N
,
,
(8)
v
+
=
max
⁡
Δ
t
j
∣
j
∈
1
,
…
,
M
,
(9)
v
-
=
min
⁡
Δ
t
j
∣
j
∈
1
,
…
,
M
,
(10)
min
⁡
v
+
-
v
-
.Formulas (8) and (9) refer to the maximum and minimum “idle time” for a solution, respectively, while formula (10) represents the objective as the difference between both the maximum (v
+) and the minimum (v
-) “idle times” (Δ
t
j, in (7)) for the same solution. Figure 6 shows a simple example of the robustness of two solutions, where solution “a” is more robust than solution “b.” In solution “a” the delay of flight 3 will not affect flight 4 and needs to be considerably larger to affect flight 5, whereas in solution “b” small delays in flight 3 will affect flight 4.Figure 6
Examples of a range of idle times for different solutions, where solution “a” is clearly better than solution “b” for formula (10).Formula (10) assumes that all flights within a solution have been assigned, which may not be the case in certain circumstances, where the maximum possible assignment is lower than a full assignment (simulation) or where the initial solution(s) has some flights remaining unassigned, and an example of this is shown in Figure 7(a). It should also be observed that this objective may conflict with the maximum assignment objective (first objective), as shown in Figure 7(b), where solution “f” has a smaller Δ
v compared to solution “e,” which is based on formula (10). This means that solution “f” is considered more robust, but solution “e” would be preferable because it achieves more assignments. Thus, it will be necessary to select the objective weights appropriately, where the fitness is a weighted sum of the different objectives, in order to ensure the correct selection of the solution.Figure 7
Examples of a range of idle times for different solutions using formula (10).Bolat [9, 17] extended the model by minimising the variance of the idle times.(11)
min
⁡
∑
j
=
1
M
Δ
t
j
-
∑
j
=
1
M
Δ
t
j
M
︸
mean  “idle  time”
2


### 3.3. Reduce Reassignment on Disruption

The ability to reassign all flights directly affected by a disruption is desirable, without the need to reassign other flights. The intention here is to generate schedules which take account of this objective, allowing such reassignment to be performed more frequently.One way to achieve this objective would be to count the number of assignments between which a reassignment could be placed when necessary. Whether the reassigned flights are on the same pier/side, as well as how many reassignments could be absorbed by a pair of assigned flights, must all be taken into account. Figure8 shows the ability of the “idle time” between the two flights 1 and 8 to accommodate flights 3, 4, and 5, should one of them be delayed. Its reassignment to B
S
S
1 may be sufficient, thus avoiding transfer of any delay to other assignments.Figure 8
An example of the capacity to absorb reassignments.The following model is proposed where the capacity to absorb reassignments may be achieved by weighting each reassignment by the Inverse of 1 plus the distance between the BSSs (d
i
n
′, presented in Section 2.1), given that such distance (d
i
n
′) may be zero, where all flights are ordered by their base start time (τ
j), as shown in (12) and (13).The intention is to use the number of flights which could be reassigned between two flights already consecutively assigned to the same BSS, without the need to reassign either of these flights in order to achieve this. Equation (12) states that φ
j
k
l is equal to 1 if flight k’s base service duration does not overlap with the base service duration of both flights j and l (j
<
l), both of which are assigned to B
S
S
i, and there is no other flight r between these (j
<
r
<
l) already assigned to the same BSS as flight j; otherwise φ
j
k
l is zero, as shown in (12)
φ
j
k
l
=
1
i
f
j
<
k
<
l
,
e
j
≤
τ
k
,
e
k
≤
τ
l
∀
j
,
k
,
l
∈
1
,
…
,
M
∄
r
∈
j
,
…
,
l
,
w
i
t
h
y
i
j
p
=
y
i
r
q
=
1
for  any
p
∈
1
,
…
,
P
j
,
q
∈
1
,
…
,
P
r
,
y
i
k
z
=
0
∀
z
∈
1
,
…
,
P
k
,
0
otherwise
.The objective is to maximise formula (13), which weights the contribution of each potential reassignment, based on which BSS the reassigned flight was originally assigned to. The underlying idea is that closer reassignments are preferred to more distant ones, but the ability to reassign without affecting other assignments is preferable. (13)
∑
i
=
1
N
∑
j
=
1
M
∑
p
=
1
P
j
∑
l
=
j
+
1
M
∑
q
=
1
P
l
∑
k
=
j
+
1
l
-
1
y
i
j
p
∗
y
i
l
q
∗
φ
i
k
l
∗
∑
n
≠
i
N
∑
z
=
1
P
k
y
n
k
z
d
i
n
′
+
1
︸
effect  of  the  BSS  assigned
.The representation in formula (13) also needs to include the border cases relating to the first and last assignments in BSS as shown in Figure 9.Figure 9
An example of border assignments in baggage sorting stations.The previous formula could be extended to cover the border cases by assigning two extra dummy flights to all available BSSs; first,j
=
0, with the end time being the start of the time period studied, and the second flight, j
=
M
+
1, with the start time being the completion of the time period studied, which for our time period would be e
0
=
0 and (
e
M
+
1
-
T
M
+
1
)
=
24hr, with y
i
0
=
y
i
(
M
+
1
)
=
1
∀
i
∈
[
1
,
…
,
N
], as shown in (14)
∑
i
=
1
N
∑
n
≠
i
N
∑
j
=
0
M
∑
p
=
1
P
j
∑
l
=
j
+
1
M
+
1
∑
q
=
1
P
l
∑
k
=
j
+
1
l
-
1
∑
z
=
1
P
k
y
i
j
p
∗
y
i
l
q
∗
y
n
k
z
∗
φ
i
k
l
d
i
n
B
S
S
s
+
1
.


### 3.4. Area of Reduction in Service (ARS)

Bolat [16] examined the minimisation of the range of idle time and the difference between the maximum and minimum idle times for the AGAP, which was later extended in [9] to consider both minimisation of the idle time range and minimisation of the idle time variance (Section 3.2). However, these do not take account of the influence or effect which the disruptions have on the schedule due to the time of their occurrence, as all reductions in service time are treated as being the same, irrespective of the time period considered in the whole “planned schedule” for the given set of departures. It is anticipated that the more congested time periods in the “planned schedule” will also represent periods where disruptions are more likely to occur and propagate, extending their effect and further increasing costs.LMAP and UMAP values are an indication of the difficulty of the problem, and these may be obtained from the distribution of the number of flights requiring service over time, as shown in Figures10 and 11.Figure 10
Flight distributions with LMAP and UMAP for 194 flights on 16th December 2009 at Terminal 1 of London Heathrow airport (identified as H1T091216).Figure 11
Flight distributions with LMAP and UMAP for 163 flights on 1st March 2010 at Terminal 1 of London Heathrow airport (identified as H1T100301).Delay during a high flight density period is more likely to propagate given that fewer resources will be available to absorb any reassignment without repercussions on other flights. At the same time, it is these cases where it is most difficult to keep a sufficiently large gap between assignments to the same BSS. In assessing the importance of the time of day when service reduction is performed, it is suggested that it would be advisable to increase penalisation of flight assignments with reduced service time when there is a higher flight density. This can be accomplished by calculating the required number of BSSs at different times of the day, as shown in Figures10 and 11.Average Assignments Point. The function f
u
(
t
) refers to the distribution of flights over time with t
s being the schedule starting time and t
e being the schedule end time, as shown in Figure 12. The Average Assignment Point (AAP) is here defined as the number of BSSs for which the distribution of flights would be uniform, which can be calculated by (15)
A
A
P
=
∫
t
s
t
e
f
u
t
d
t
t
e
-
t
s
.Figure 12
Distribution of flights over time and Area Reduction in Service Objective.A
j is the density distribution area for the time period from the target service time of flight j (t
j) to the end of service time for the previous flight assigned to the same BSS; for example, in Figure 12  A
18
=
∫
t
18
e
9
‍
f
u
(
t
)
d
t. A
A
P
j is here defined as the mean number of flights over the target start time for flight j and the end of service time for the previous flight assigned to the same BSS as flight j; for example, A
A
P
18
=
A
18
/
(
e
9
-
t
18
). So the contribution to the objective for assignment j is A
A
P
j
/
A
A
P times the reduction in service previously considered in Section 3.1. This corresponds to values greater than one for dense flight regions of the schedule and less than one for underused regions.The approach proposed intends to penalise more those reductions in flight service time during time periods where more flights require servicing compared to those periods with less service load by means of the AAP and the distribution of flights over time, as shown in (16). The idea is that flights which require servicing during congested periods are more likely to have a knock-on effect compared to those in less congested periods. It is therefore preferable not to reduce the service time of flights at more congested times so much, in order to limit the effect of potential delays: (16)
f
2
=
-
1
A
A
P
∗
∑
i
=
1
N
∑
j
=
1
M
∑
p
=
1
P
j
y
i
j
p
∗
A
j
.This approach is compared with some of the other approaches described here in Section4.


### 3.5. Subarea of Reduction in Service (SARS)

In this case, the objective is based on the area between both flight densities “without reduction in service time” and “with reduction in service time” shown in Figure13, being here called a subarea. The flight density subarea corresponds to the area from the flight target start time (described in Section 2.1) to the previously assigned flight’s end service time, which lies between both flight density distributions (f
u and f
l), where a reduction in service time is not permitted and when all the buffer time has been reduced (considering only the base service duration, T). (17)
A
j
′
=
∫
t
j
e
f
u
t
-
f
l
t
d
t
.Figure 13
Subarea Reduction in Service Objective.The approach now proposed has a fitness which covers the subarea divided by the difference between the UMAP and the LMAP for all the assigned flights.(18)
f
2
=
-
∑
i
=
1
N
∑
j
=
1
M
∑
p
=
1
P
j
y
i
j
p
∗
A
j
′
U
M
A
P
-
L
M
A
P
.Similarly to the Area Reduction in Service (Section3.4), the main idea is to penalise more heavily those reductions in flight service time which occur in regions with high flight density, as these are more likely to further disrupt the schedule in case of delays. Adding both approaches together with the TRS will be seen in Section 4 to increase the robustness as compared to using each approach individually.Nevertheless, in cases where the LMAP is equal to the UMAP, the SARS approach cannot be used; an example is shown in Figure14. In these cases, the area between both flight densities could be used, an approach that is here named Base Subarea Reduction in Service (BSARS). As will be seen in Section 4.1, this approach significantly increases the robustness and widens the range of the quantity of BSSs, in which it performs better when compared to the ARS, TRS, and SARS.Figure 14
Example of distributions with the same LMAP and UMAP.These approaches are compared with some of the other approaches introduced in this paper in Section4.


### 3.6. Unsupervised Estimated Stochastic Reduction in Service (UESRS)

Lim and Wang [18] proposed a stochastic programming model for the AGAP with a robustness cost of conflicts, which is estimated by a function, v
(
j
,
k
). Flights are ordered by their base starting service time, so the gap between two flights j and k, l
(
j
,
k
), assigned to the same BSS, where j
<
k and j, k
∈
[
1
,
…
,
M
], is the difference between flight k’s target service time and the prior assigned flight j’s end service time (19), where l
(
j
,
k
)
=
-
r
j for l
(
j
,
k
)
<
0, as shown in Figure 15. v
(
j
,
k
) is used to estimate the mean conflict probability between flights j and k assigned to the same BSS, which is a function of the gap l
(
j
,
k
), where larger gaps between assignments to the same BSSs result in lesser probability of real flight conflicts. v
(
j
,
k
) is normalised in (20). (19)
l
j
,
k
=
t
k
-
e
j
t
k
=
τ
k
-
B
k
,
(20)
E
p
j
,
k
=
∑
i
=
1
N
y
i
j
∗
y
i
k
∗
v
j
,
k
-
v
m
i
n
j
,
k
v
m
a
x
j
,
k
-
v
m
i
n
j
,
k
,
(21)
f
2
=
-
∑
j
=
1
M
-
1
∑
k
=
j
+
1
M
E
p
j
,
k
.Figure 15
Overlap between two flightsj and k assigned to the same BSS.The definition ofv
(
j
,
k
) comes from the application domain, in the absence of historical data; some unsupervised estimation functions were introduced in [18]. Figure 16 shows the penalty (y-axis) incurred for different unsupervised estimation functions as a function of the gap (x-axis). Negative values refer to reductions in service time between two assignments to the same BSS, which are heavily penalised as they may require reassignment should delays occur, whereas positive gaps are penalised less. Wider gaps between two assignments reduce the need to reassign delayed flights, given that the delay has to be larger than the gap in order to affect the following assignment to the same BSS. Similarly, starting the service earlier may not require the flight to be reassigned because the duration of earliness has to be lower than the gap in order to affect the previous assignment. Both earliness and delay probabilities decrease as the gap increases. Sufficiently large gaps may also be used on the day of operation by disrupted flights which need to be reassigned, such that the detrimental effect of disruptions on that day is reduced.Figure 16
Penalty for different unsupervised estimation functions based on the gap between assignments.The unsupervised estimation functions introduced in [18] are presented below and are shown in Figure 16:(1)
Linear estimation is as follows.(22)
v
j
,
k
=
-
l
j
,
k
.
(2)
Exponential estimation is as follows.(23)
v
j
,
k
=
e
-
β
∗
l
j
,
k
.
(3)
Inverse estimation is as follows.(24)
v
j
,
k
=
b
l
j
,
k
+
b
i
f
l
j
,
k
>
0
,
1
otherwise
.
The value of the constant “b” changes the penalisation as shown in Figure 16, so a higher “b” increases the penalisation and a lower “b” decreases it. An appropriate value should be selected to properly weight the influence of the potential conflicts. Lim and Wang [18] used b
=
15 minutes, which proved to provide rather poor results when compared with the exponential estimation function, which may partly be caused by the fixed cost when l
(
j
,
k
)
<
0 (dark red dash line, Figure 16), whereas in the exponential estimation function (purple dash line, Figure 16) this is not the case. The value used for “b” may be too great, and a lower value would make this estimation function provide values closer to those provided by the exponential estimation function which provided fitter solutions in the results presented in [18]. Consequently, a value b
=
6 was seen in the experiments studied in Section 4.1 to provide better results than when b
=
15. In general, an even lower value did appear to perform better in some instances but not as well as b
=
6, as shown in Section 4.1.
The Inverse estimation function as considered in [18] treats all gaps smaller than the buffer time equally, which does not represent a real case since smaller gaps between flights are more likely to result in conflicts than larger ones on the day of operation. Given this and that the exponential estimation function performs best and treats all gaps differently, it is proposed that all of the gaps be treated differently, as shown by the modified version which is herein named “Offset inverse” (25), which is shown in Figure 16 for b
=
15 (green line). (25)
v
j
,
k
=
b
l
j
,
k
-
min
⁡
l
j
,
k
+
b
.
(4)
Sublinear estimation is as follows.(26)
v
j
,
k
=
cos
⁡
π
∗
l
j
,
k
l
m
a
x
i
f
l
j
,
k
>
0
,
1
otherwise
.
This estimation also suffers from the same problem as the Inverse estimation and may be improved by offsetting its value so that the maximum penalisation corresponds tol
m
a
x and the minimum corresponds to l
m
i
n (27), which is shown in Figure 16 for γ
=
0. (27)
v
j
,
k
=
cos
⁡
π
∗
l
j
,
k
-
l
m
i
n
+
γ
l
m
a
x
-
l
m
i
n
+
γ
.The gap definition used takes account of the buffer time, as the target service duration is the base service duration (T
k) plus the buffer time (B
k) for the flight. This makes the estimation functions dependent on the buffer time of each flight, as shown in Figure 17 for two buffer times of 30 and 15 min each. When the buffer time is the same irrespective of the flights, B
k
=
B
∀
k
∈
[
1
,
…
,
M
], as considered in [18], the cost is the same irrespective of the flight, depending only on the separation between consecutive flight assignments, but this is not the situation when the buffer time depends on the flight, namely, long, medium, or short distance flights, which are the cases studied here.Figure 17
Penalty for different unsupervised estimation functions based on the gap between assignments for different buffer times.


### 3.7. Reduction in the Number of Conflicts (RNC)

A stochastic approach for improving schedule performance is described here, when disruptions occur on the day of operation. A similar approach was used in [19] where random delay scenarios are generated in the “planning stage” which are used to account for the potential disruptions in the schedule on the day of implementation by means of calculating the expected semideviation risk measure [20] for all those delay scenarios.This approach is based on reducing the number of conflicts on the day of operation. Given that the real perturbed conditions will not be available until the day the schedule is implemented, these perturbed conditions are simulated by examining a set of perturbed base cases,S, which may be obtained in different ways, such as randomly, for example, from historical data or calculated using known distribution(s) from information available at the time of generating the assignments.It is considered a set of perturbed schedulesS, which simulate the perturbations on the day of operation. A new variable is introduced c
j
s, which for a given solution of assignments has the value of 1 if flight j is in conflict with another flight in the perturbed schedule s
∈
S or zero otherwise. The average number of conflicts in the set of perturbed schedules S is calculated by (28), which is a measure of the solution robustness. (28)
f
2
=
-
1
S
∗
∑
s
∈
S
∑
j
=
1
M
c
j
s
.When reassigning conflicting flights or subsequent flights to other BSSs, a situation can arise where the reassigned flight interferes with the subsequent flight at the new BSS, a so-called secondary conflict. Some reassignment may therefore have a downstream effect on the overall schedule, producing more conflicting flights, in turn requiring further reassignments, thus potentially increasing the difficulty of the problem later on.The above version considers all the conflicts to be of the same importance, but it is preferable to have conflicts which do not have repercussions later, that is, can be reassigned to another BSS without affecting any of the assignments already in existence. To account for this situation, a new variablec
j
s
′ is defined which takes the value of 1 if the reassignment of conflicting flight j in a perturbed schedule “s” affects other assignments already in existence or zero otherwise. The objective is presented as (29) where the constant, α,  0
≤
α
≤
1, denotes the importance of the conflicting flight repercussions on other assignments; α
=
0 corresponds to the case where no account is taken of any repercussion on other assignments, which corresponds in turn to (28), and α
=
1 corresponds to the cases in which both the conflicting flights and their repercussions on other assignments are considered to be of the same importance. α
>
1 refers to the cases where more importance is given to the repercussions of a conflict on assignments other than the conflict itself. (29)
f
2
=
-
1
S
∗
∑
s
∈
S
∑
j
=
1
M
c
j
s
+
α
∗
c
j
s
′
.Calculation of the conflicts is time consuming and even more so if the effect of the conflict repercussions is also calculated, which is further aggravated by the need to use a large number of schedules in the perturbed setS in order to achieve a good representation of all the potential situations. The execution time is one of the disadvantages of using this approach as was experienced when executing the experiments presented in Section 4.


### 3.8. Probability of Conflict Based on the Gap (PCBG)

The previous approach would normally require a large number of perturbed data sets, which makes its application very slow. Given that we are still interested in reducing the number of conflicts, but without the heavy cost in speed imposed by the RNC approach, then it is proposed to use the probability of having a conflict in a given “idle time” for each flight. This can be easily obtained if the delay distribution is known. In the ABSSAP, the early arrival of a flight does not normally affect the assignment to the Baggage Sorting Station Selection (BSSS) as this does not extend the time in which the BSS is required for servicing the flight, but transportation of the baggage already in the BSS to the flight may start earlier, so no earlier arrival is considered here, in the case of a normal folded distribution (no negative numbers and with zero mean) and independent delays for flights. This could also be extended to other distributions and to nonzero means.Independent delays are considered initially, where a conflict between two consecutive flights is independent of other flights assigned to the same BSS. The assignments to different BSSs are independent from the point of view of conflicts. The probability of two consecutive flights having overlapping service times (a conflict) corresponds to the sum of the product between the probability of a sufficiently large delay on the part of the previous flight assigned to that BSS and the probability of the next consecutive assignment to the same BSS not being sufficiently delayed, as shown in Figure18, where η
j
(
t
) is the probability density function for flight j and t
0
=
τ
k
-
e
j.Figure 18
Probability of a conflict between two consecutive flights based on the intervening gap.Equation (30) is the probability of a conflict, which is equal to the probability of exceeding the gap between both assignments, multiplied by the probability of the following assignment not being sufficiently delayed to avoid conflicts. p
(
j
,
k
) is the probability of a conflict between two flights j and k assigned consecutively to the same BSS with a gap between them of t
0
=
τ
k
-
e
j given their respective probability density distribution of delay η
j
(
t
) and η
k
(
t
). (30)
p
j
,
k
=
∫
τ
k
-
e
j
∞
η
j
t
∗
1
-
∫
0
t
-
τ
k
-
e
j
η
k
x
d
x
d
t
.A simplification of the conflict probability is represented in (31), which uses the “Riemann integral” approach for the range of delays between (
τ
k
-
e
j
) and four times the delay distribution standard deviation of flight j, σ
j, and a time increment of Δ
t. (31)
p
j
,
k
=
∑
i
=
0
4
∗
σ
j
-
τ
k
-
e
j
/
Δ
t
η
j
τ
k
-
e
j
+
i
+
1
2
∗
Δ
t
∗
1
-
∑
l
=
0
i
η
k
l
+
1
2
∗
Δ
t
∗
Δ
t
∗
Δ
t
.


## 3.1. Minimise Reduction in Service Time

It may be possible to gain robustness by reordering assignments between BSSs so that “idle time” between flights consecutively assigned to the same BSS is greater, as shown in Figure4.Figure 4
Simple example of two schedules, with the same flights, where one is obviously more robust in respect of perturbations than the other.Figure4 shows two potential solutions with different robustness. While any delay to flight “a” in the “less robust” solution will certainly affect flight “b,” which will in turn have to be reassigned to another BSS, in the “more robust” solution a delay in flight “a” will not affect flight “b.” So the “more robust” solution is preferable to that of the “less robust” solution.Mangoubi and Mathaisel [15] proposed the use of “buffer times” between two flights which are consecutively assigned to the same gate in order to obtain robust assignments, defining the reduction in service as that part of buffer time which overlaps with the previous assignment to the same gate. Given the detrimental effects that the reduction in service time has on the robustness of the assignment as against real-life delays, it is advisable to minimise the total reduction in service time, thus maximising buffer times. This objective, total reduction in service time (TRS), can be expressed by formula (4).(4)
min
⁡
∑
j
=
1
M
∑
p
=
1
P
j
r
j
p
.If a fitness function is defined as a weighted sum of the different objectives, as used here, and the robustness uses the reduction in service time only for assigned flights, then particular weights may compromise the importance of the main objective (maximisation of the number of assignments), where solutions with lower numbers of assignments are favoured over those with higher numbers of assignments because of the robustness objective.The use of formula (4) treats any reduction in service equally, so it does not make a distinction between reducing all the assignment buffer time and allocating it to another flight from another solution where both flights share the available “idle time,” as shown in Figure 5. However, flight 3 with maximum reduction in service in solution “b” (Figure 5) will be unable to absorb any delay on the day of operation, although the same flight in the alternative solution “a” (Figure 5) will be able to do so, making it the preferable choice.Figure 5
Simple example of a more robust schedule of four flights and two sorting stations using the reduction in buffer time.A nonlinear penalty function is required to take account of this, such as formulas (5) and (6).(5)
u
i
j
p
=
arctan
⁡
τ
j
-
e
k
B
j
-
Π
2
i
f
k
<
j
,
y
i
j
p
=
y
i
k
q
=
1
,
∑
l
=
k
+
1
j
∑
p
=
1
P
l
y
i
l
p
=
1
,
0
o
t
h
e
r
w
i
s
e
,
(6)
max
⁡
∑
i
=
1
N
∑
j
=
1
M
-
1
∑
p
=
1
P
j
u
i
j
p
.Formula (5) defines the penalty u
i
j
p for assigning activity p of flight j to sorting station i, with the total penalty being represented by formula (6). There is only a penalty for consecutive assignments to the same sorting station, between the assignment and the previous assignment to the same sorting station. The reason for using arctangent is based on the properties of this function which imposes stronger penalties around the point at which the flights are assigned service time without any buffer time to the point at which all the buffer time is retained (by dividing the member of arctangent by B
k). The effect is reduced as the time separation between the assignments increases, but the contribution of any time separation is always considered negative. The constant π
/
2 is used so that there is always a cost associated with assigning two flights to the same sorting station; otherwise it would be a benefit. Considering the arctangent for each flight increases the individual effect on the overall objective, in contrast to using the arctangent of the sum of all of the flight contributions.The robustness objective represented by formula (6) also takes account of the objective of “minimising the service reduction,” as being the sum of the reduction in buffer time of each flight assigned.


## 3.2. Distribute Idle Time

Bolat [16] proposed the distribution of “idle time” uniformly amongst gates for the AGAP. In the case of the ABSSAP, the distribution of the “idle time” uniformly amongst BSSs can be considered, where “idle time” is calculated as the time between the start of service time on a flight’s assignment to BSS and the end of service time of the flight assigned immediately before the same baggage sorting station, as shown in formula (7). Whereas the buffer time implies preference for a particular gap size between consecutive assignments to the same sorting station, an “idle time” does not. When using the “idle time” it is normally intended that the gap size should be as large as possible. The reason for this is to increase the probability that, even with a delay, completion of service time for flight presently assigned will still be earlier than the start of service time for the next flight assigned to the same BSS. (7)
Δ
t
j
=
s
j
-
e
k
for
j
<
k
,
y
i
j
p
=
y
i
k
q
=
1
where
∄
l
∈
k
,
…
,
M
with
τ
l
<
e
k
;
j
,
k
∈
1
,
…
,
M
,
p
∈
1
,
…
,
P
j
,
q
∈
1
,
…
,
P
k
,
i
∈
1
,
…
,
N
,
,
(8)
v
+
=
max
⁡
Δ
t
j
∣
j
∈
1
,
…
,
M
,
(9)
v
-
=
min
⁡
Δ
t
j
∣
j
∈
1
,
…
,
M
,
(10)
min
⁡
v
+
-
v
-
.Formulas (8) and (9) refer to the maximum and minimum “idle time” for a solution, respectively, while formula (10) represents the objective as the difference between both the maximum (v
+) and the minimum (v
-) “idle times” (Δ
t
j, in (7)) for the same solution. Figure 6 shows a simple example of the robustness of two solutions, where solution “a” is more robust than solution “b.” In solution “a” the delay of flight 3 will not affect flight 4 and needs to be considerably larger to affect flight 5, whereas in solution “b” small delays in flight 3 will affect flight 4.Figure 6
Examples of a range of idle times for different solutions, where solution “a” is clearly better than solution “b” for formula (10).Formula (10) assumes that all flights within a solution have been assigned, which may not be the case in certain circumstances, where the maximum possible assignment is lower than a full assignment (simulation) or where the initial solution(s) has some flights remaining unassigned, and an example of this is shown in Figure 7(a). It should also be observed that this objective may conflict with the maximum assignment objective (first objective), as shown in Figure 7(b), where solution “f” has a smaller Δ
v compared to solution “e,” which is based on formula (10). This means that solution “f” is considered more robust, but solution “e” would be preferable because it achieves more assignments. Thus, it will be necessary to select the objective weights appropriately, where the fitness is a weighted sum of the different objectives, in order to ensure the correct selection of the solution.Figure 7
Examples of a range of idle times for different solutions using formula (10).Bolat [9, 17] extended the model by minimising the variance of the idle times.(11)
min
⁡
∑
j
=
1
M
Δ
t
j
-
∑
j
=
1
M
Δ
t
j
M
︸
mean  “idle  time”
2


## 3.3. Reduce Reassignment on Disruption

The ability to reassign all flights directly affected by a disruption is desirable, without the need to reassign other flights. The intention here is to generate schedules which take account of this objective, allowing such reassignment to be performed more frequently.One way to achieve this objective would be to count the number of assignments between which a reassignment could be placed when necessary. Whether the reassigned flights are on the same pier/side, as well as how many reassignments could be absorbed by a pair of assigned flights, must all be taken into account. Figure8 shows the ability of the “idle time” between the two flights 1 and 8 to accommodate flights 3, 4, and 5, should one of them be delayed. Its reassignment to B
S
S
1 may be sufficient, thus avoiding transfer of any delay to other assignments.Figure 8
An example of the capacity to absorb reassignments.The following model is proposed where the capacity to absorb reassignments may be achieved by weighting each reassignment by the Inverse of 1 plus the distance between the BSSs (d
i
n
′, presented in Section 2.1), given that such distance (d
i
n
′) may be zero, where all flights are ordered by their base start time (τ
j), as shown in (12) and (13).The intention is to use the number of flights which could be reassigned between two flights already consecutively assigned to the same BSS, without the need to reassign either of these flights in order to achieve this. Equation (12) states that φ
j
k
l is equal to 1 if flight k’s base service duration does not overlap with the base service duration of both flights j and l (j
<
l), both of which are assigned to B
S
S
i, and there is no other flight r between these (j
<
r
<
l) already assigned to the same BSS as flight j; otherwise φ
j
k
l is zero, as shown in (12)
φ
j
k
l
=
1
i
f
j
<
k
<
l
,
e
j
≤
τ
k
,
e
k
≤
τ
l
∀
j
,
k
,
l
∈
1
,
…
,
M
∄
r
∈
j
,
…
,
l
,
w
i
t
h
y
i
j
p
=
y
i
r
q
=
1
for  any
p
∈
1
,
…
,
P
j
,
q
∈
1
,
…
,
P
r
,
y
i
k
z
=
0
∀
z
∈
1
,
…
,
P
k
,
0
otherwise
.The objective is to maximise formula (13), which weights the contribution of each potential reassignment, based on which BSS the reassigned flight was originally assigned to. The underlying idea is that closer reassignments are preferred to more distant ones, but the ability to reassign without affecting other assignments is preferable. (13)
∑
i
=
1
N
∑
j
=
1
M
∑
p
=
1
P
j
∑
l
=
j
+
1
M
∑
q
=
1
P
l
∑
k
=
j
+
1
l
-
1
y
i
j
p
∗
y
i
l
q
∗
φ
i
k
l
∗
∑
n
≠
i
N
∑
z
=
1
P
k
y
n
k
z
d
i
n
′
+
1
︸
effect  of  the  BSS  assigned
.The representation in formula (13) also needs to include the border cases relating to the first and last assignments in BSS as shown in Figure 9.Figure 9
An example of border assignments in baggage sorting stations.The previous formula could be extended to cover the border cases by assigning two extra dummy flights to all available BSSs; first,j
=
0, with the end time being the start of the time period studied, and the second flight, j
=
M
+
1, with the start time being the completion of the time period studied, which for our time period would be e
0
=
0 and (
e
M
+
1
-
T
M
+
1
)
=
24hr, with y
i
0
=
y
i
(
M
+
1
)
=
1
∀
i
∈
[
1
,
…
,
N
], as shown in (14)
∑
i
=
1
N
∑
n
≠
i
N
∑
j
=
0
M
∑
p
=
1
P
j
∑
l
=
j
+
1
M
+
1
∑
q
=
1
P
l
∑
k
=
j
+
1
l
-
1
∑
z
=
1
P
k
y
i
j
p
∗
y
i
l
q
∗
y
n
k
z
∗
φ
i
k
l
d
i
n
B
S
S
s
+
1
.


## 3.4. Area of Reduction in Service (ARS)

Bolat [16] examined the minimisation of the range of idle time and the difference between the maximum and minimum idle times for the AGAP, which was later extended in [9] to consider both minimisation of the idle time range and minimisation of the idle time variance (Section 3.2). However, these do not take account of the influence or effect which the disruptions have on the schedule due to the time of their occurrence, as all reductions in service time are treated as being the same, irrespective of the time period considered in the whole “planned schedule” for the given set of departures. It is anticipated that the more congested time periods in the “planned schedule” will also represent periods where disruptions are more likely to occur and propagate, extending their effect and further increasing costs.LMAP and UMAP values are an indication of the difficulty of the problem, and these may be obtained from the distribution of the number of flights requiring service over time, as shown in Figures10 and 11.Figure 10
Flight distributions with LMAP and UMAP for 194 flights on 16th December 2009 at Terminal 1 of London Heathrow airport (identified as H1T091216).Figure 11
Flight distributions with LMAP and UMAP for 163 flights on 1st March 2010 at Terminal 1 of London Heathrow airport (identified as H1T100301).Delay during a high flight density period is more likely to propagate given that fewer resources will be available to absorb any reassignment without repercussions on other flights. At the same time, it is these cases where it is most difficult to keep a sufficiently large gap between assignments to the same BSS. In assessing the importance of the time of day when service reduction is performed, it is suggested that it would be advisable to increase penalisation of flight assignments with reduced service time when there is a higher flight density. This can be accomplished by calculating the required number of BSSs at different times of the day, as shown in Figures10 and 11.Average Assignments Point. The function f
u
(
t
) refers to the distribution of flights over time with t
s being the schedule starting time and t
e being the schedule end time, as shown in Figure 12. The Average Assignment Point (AAP) is here defined as the number of BSSs for which the distribution of flights would be uniform, which can be calculated by (15)
A
A
P
=
∫
t
s
t
e
f
u
t
d
t
t
e
-
t
s
.Figure 12
Distribution of flights over time and Area Reduction in Service Objective.A
j is the density distribution area for the time period from the target service time of flight j (t
j) to the end of service time for the previous flight assigned to the same BSS; for example, in Figure 12  A
18
=
∫
t
18
e
9
‍
f
u
(
t
)
d
t. A
A
P
j is here defined as the mean number of flights over the target start time for flight j and the end of service time for the previous flight assigned to the same BSS as flight j; for example, A
A
P
18
=
A
18
/
(
e
9
-
t
18
). So the contribution to the objective for assignment j is A
A
P
j
/
A
A
P times the reduction in service previously considered in Section 3.1. This corresponds to values greater than one for dense flight regions of the schedule and less than one for underused regions.The approach proposed intends to penalise more those reductions in flight service time during time periods where more flights require servicing compared to those periods with less service load by means of the AAP and the distribution of flights over time, as shown in (16). The idea is that flights which require servicing during congested periods are more likely to have a knock-on effect compared to those in less congested periods. It is therefore preferable not to reduce the service time of flights at more congested times so much, in order to limit the effect of potential delays: (16)
f
2
=
-
1
A
A
P
∗
∑
i
=
1
N
∑
j
=
1
M
∑
p
=
1
P
j
y
i
j
p
∗
A
j
.This approach is compared with some of the other approaches described here in Section4.


## 3.5. Subarea of Reduction in Service (SARS)

In this case, the objective is based on the area between both flight densities “without reduction in service time” and “with reduction in service time” shown in Figure13, being here called a subarea. The flight density subarea corresponds to the area from the flight target start time (described in Section 2.1) to the previously assigned flight’s end service time, which lies between both flight density distributions (f
u and f
l), where a reduction in service time is not permitted and when all the buffer time has been reduced (considering only the base service duration, T). (17)
A
j
′
=
∫
t
j
e
f
u
t
-
f
l
t
d
t
.Figure 13
Subarea Reduction in Service Objective.The approach now proposed has a fitness which covers the subarea divided by the difference between the UMAP and the LMAP for all the assigned flights.(18)
f
2
=
-
∑
i
=
1
N
∑
j
=
1
M
∑
p
=
1
P
j
y
i
j
p
∗
A
j
′
U
M
A
P
-
L
M
A
P
.Similarly to the Area Reduction in Service (Section3.4), the main idea is to penalise more heavily those reductions in flight service time which occur in regions with high flight density, as these are more likely to further disrupt the schedule in case of delays. Adding both approaches together with the TRS will be seen in Section 4 to increase the robustness as compared to using each approach individually.Nevertheless, in cases where the LMAP is equal to the UMAP, the SARS approach cannot be used; an example is shown in Figure14. In these cases, the area between both flight densities could be used, an approach that is here named Base Subarea Reduction in Service (BSARS). As will be seen in Section 4.1, this approach significantly increases the robustness and widens the range of the quantity of BSSs, in which it performs better when compared to the ARS, TRS, and SARS.Figure 14
Example of distributions with the same LMAP and UMAP.These approaches are compared with some of the other approaches introduced in this paper in Section4.


## 3.6. Unsupervised Estimated Stochastic Reduction in Service (UESRS)

Lim and Wang [18] proposed a stochastic programming model for the AGAP with a robustness cost of conflicts, which is estimated by a function, v
(
j
,
k
). Flights are ordered by their base starting service time, so the gap between two flights j and k, l
(
j
,
k
), assigned to the same BSS, where j
<
k and j, k
∈
[
1
,
…
,
M
], is the difference between flight k’s target service time and the prior assigned flight j’s end service time (19), where l
(
j
,
k
)
=
-
r
j for l
(
j
,
k
)
<
0, as shown in Figure 15. v
(
j
,
k
) is used to estimate the mean conflict probability between flights j and k assigned to the same BSS, which is a function of the gap l
(
j
,
k
), where larger gaps between assignments to the same BSSs result in lesser probability of real flight conflicts. v
(
j
,
k
) is normalised in (20). (19)
l
j
,
k
=
t
k
-
e
j
t
k
=
τ
k
-
B
k
,
(20)
E
p
j
,
k
=
∑
i
=
1
N
y
i
j
∗
y
i
k
∗
v
j
,
k
-
v
m
i
n
j
,
k
v
m
a
x
j
,
k
-
v
m
i
n
j
,
k
,
(21)
f
2
=
-
∑
j
=
1
M
-
1
∑
k
=
j
+
1
M
E
p
j
,
k
.Figure 15
Overlap between two flightsj and k assigned to the same BSS.The definition ofv
(
j
,
k
) comes from the application domain, in the absence of historical data; some unsupervised estimation functions were introduced in [18]. Figure 16 shows the penalty (y-axis) incurred for different unsupervised estimation functions as a function of the gap (x-axis). Negative values refer to reductions in service time between two assignments to the same BSS, which are heavily penalised as they may require reassignment should delays occur, whereas positive gaps are penalised less. Wider gaps between two assignments reduce the need to reassign delayed flights, given that the delay has to be larger than the gap in order to affect the following assignment to the same BSS. Similarly, starting the service earlier may not require the flight to be reassigned because the duration of earliness has to be lower than the gap in order to affect the previous assignment. Both earliness and delay probabilities decrease as the gap increases. Sufficiently large gaps may also be used on the day of operation by disrupted flights which need to be reassigned, such that the detrimental effect of disruptions on that day is reduced.Figure 16
Penalty for different unsupervised estimation functions based on the gap between assignments.The unsupervised estimation functions introduced in [18] are presented below and are shown in Figure 16:(1)
Linear estimation is as follows.(22)
v
j
,
k
=
-
l
j
,
k
.
(2)
Exponential estimation is as follows.(23)
v
j
,
k
=
e
-
β
∗
l
j
,
k
.
(3)
Inverse estimation is as follows.(24)
v
j
,
k
=
b
l
j
,
k
+
b
i
f
l
j
,
k
>
0
,
1
otherwise
.
The value of the constant “b” changes the penalisation as shown in Figure 16, so a higher “b” increases the penalisation and a lower “b” decreases it. An appropriate value should be selected to properly weight the influence of the potential conflicts. Lim and Wang [18] used b
=
15 minutes, which proved to provide rather poor results when compared with the exponential estimation function, which may partly be caused by the fixed cost when l
(
j
,
k
)
<
0 (dark red dash line, Figure 16), whereas in the exponential estimation function (purple dash line, Figure 16) this is not the case. The value used for “b” may be too great, and a lower value would make this estimation function provide values closer to those provided by the exponential estimation function which provided fitter solutions in the results presented in [18]. Consequently, a value b
=
6 was seen in the experiments studied in Section 4.1 to provide better results than when b
=
15. In general, an even lower value did appear to perform better in some instances but not as well as b
=
6, as shown in Section 4.1.
The Inverse estimation function as considered in [18] treats all gaps smaller than the buffer time equally, which does not represent a real case since smaller gaps between flights are more likely to result in conflicts than larger ones on the day of operation. Given this and that the exponential estimation function performs best and treats all gaps differently, it is proposed that all of the gaps be treated differently, as shown by the modified version which is herein named “Offset inverse” (25), which is shown in Figure 16 for b
=
15 (green line). (25)
v
j
,
k
=
b
l
j
,
k
-
min
⁡
l
j
,
k
+
b
.
(4)
Sublinear estimation is as follows.(26)
v
j
,
k
=
cos
⁡
π
∗
l
j
,
k
l
m
a
x
i
f
l
j
,
k
>
0
,
1
otherwise
.
This estimation also suffers from the same problem as the Inverse estimation and may be improved by offsetting its value so that the maximum penalisation corresponds tol
m
a
x and the minimum corresponds to l
m
i
n (27), which is shown in Figure 16 for γ
=
0. (27)
v
j
,
k
=
cos
⁡
π
∗
l
j
,
k
-
l
m
i
n
+
γ
l
m
a
x
-
l
m
i
n
+
γ
.The gap definition used takes account of the buffer time, as the target service duration is the base service duration (T
k) plus the buffer time (B
k) for the flight. This makes the estimation functions dependent on the buffer time of each flight, as shown in Figure 17 for two buffer times of 30 and 15 min each. When the buffer time is the same irrespective of the flights, B
k
=
B
∀
k
∈
[
1
,
…
,
M
], as considered in [18], the cost is the same irrespective of the flight, depending only on the separation between consecutive flight assignments, but this is not the situation when the buffer time depends on the flight, namely, long, medium, or short distance flights, which are the cases studied here.Figure 17
Penalty for different unsupervised estimation functions based on the gap between assignments for different buffer times.


## 3.7. Reduction in the Number of Conflicts (RNC)

A stochastic approach for improving schedule performance is described here, when disruptions occur on the day of operation. A similar approach was used in [19] where random delay scenarios are generated in the “planning stage” which are used to account for the potential disruptions in the schedule on the day of implementation by means of calculating the expected semideviation risk measure [20] for all those delay scenarios.This approach is based on reducing the number of conflicts on the day of operation. Given that the real perturbed conditions will not be available until the day the schedule is implemented, these perturbed conditions are simulated by examining a set of perturbed base cases,S, which may be obtained in different ways, such as randomly, for example, from historical data or calculated using known distribution(s) from information available at the time of generating the assignments.It is considered a set of perturbed schedulesS, which simulate the perturbations on the day of operation. A new variable is introduced c
j
s, which for a given solution of assignments has the value of 1 if flight j is in conflict with another flight in the perturbed schedule s
∈
S or zero otherwise. The average number of conflicts in the set of perturbed schedules S is calculated by (28), which is a measure of the solution robustness. (28)
f
2
=
-
1
S
∗
∑
s
∈
S
∑
j
=
1
M
c
j
s
.When reassigning conflicting flights or subsequent flights to other BSSs, a situation can arise where the reassigned flight interferes with the subsequent flight at the new BSS, a so-called secondary conflict. Some reassignment may therefore have a downstream effect on the overall schedule, producing more conflicting flights, in turn requiring further reassignments, thus potentially increasing the difficulty of the problem later on.The above version considers all the conflicts to be of the same importance, but it is preferable to have conflicts which do not have repercussions later, that is, can be reassigned to another BSS without affecting any of the assignments already in existence. To account for this situation, a new variablec
j
s
′ is defined which takes the value of 1 if the reassignment of conflicting flight j in a perturbed schedule “s” affects other assignments already in existence or zero otherwise. The objective is presented as (29) where the constant, α,  0
≤
α
≤
1, denotes the importance of the conflicting flight repercussions on other assignments; α
=
0 corresponds to the case where no account is taken of any repercussion on other assignments, which corresponds in turn to (28), and α
=
1 corresponds to the cases in which both the conflicting flights and their repercussions on other assignments are considered to be of the same importance. α
>
1 refers to the cases where more importance is given to the repercussions of a conflict on assignments other than the conflict itself. (29)
f
2
=
-
1
S
∗
∑
s
∈
S
∑
j
=
1
M
c
j
s
+
α
∗
c
j
s
′
.Calculation of the conflicts is time consuming and even more so if the effect of the conflict repercussions is also calculated, which is further aggravated by the need to use a large number of schedules in the perturbed setS in order to achieve a good representation of all the potential situations. The execution time is one of the disadvantages of using this approach as was experienced when executing the experiments presented in Section 4.


## 3.8. Probability of Conflict Based on the Gap (PCBG)

The previous approach would normally require a large number of perturbed data sets, which makes its application very slow. Given that we are still interested in reducing the number of conflicts, but without the heavy cost in speed imposed by the RNC approach, then it is proposed to use the probability of having a conflict in a given “idle time” for each flight. This can be easily obtained if the delay distribution is known. In the ABSSAP, the early arrival of a flight does not normally affect the assignment to the Baggage Sorting Station Selection (BSSS) as this does not extend the time in which the BSS is required for servicing the flight, but transportation of the baggage already in the BSS to the flight may start earlier, so no earlier arrival is considered here, in the case of a normal folded distribution (no negative numbers and with zero mean) and independent delays for flights. This could also be extended to other distributions and to nonzero means.Independent delays are considered initially, where a conflict between two consecutive flights is independent of other flights assigned to the same BSS. The assignments to different BSSs are independent from the point of view of conflicts. The probability of two consecutive flights having overlapping service times (a conflict) corresponds to the sum of the product between the probability of a sufficiently large delay on the part of the previous flight assigned to that BSS and the probability of the next consecutive assignment to the same BSS not being sufficiently delayed, as shown in Figure18, where η
j
(
t
) is the probability density function for flight j and t
0
=
τ
k
-
e
j.Figure 18
Probability of a conflict between two consecutive flights based on the intervening gap.Equation (30) is the probability of a conflict, which is equal to the probability of exceeding the gap between both assignments, multiplied by the probability of the following assignment not being sufficiently delayed to avoid conflicts. p
(
j
,
k
) is the probability of a conflict between two flights j and k assigned consecutively to the same BSS with a gap between them of t
0
=
τ
k
-
e
j given their respective probability density distribution of delay η
j
(
t
) and η
k
(
t
). (30)
p
j
,
k
=
∫
τ
k
-
e
j
∞
η
j
t
∗
1
-
∫
0
t
-
τ
k
-
e
j
η
k
x
d
x
d
t
.A simplification of the conflict probability is represented in (31), which uses the “Riemann integral” approach for the range of delays between (
τ
k
-
e
j
) and four times the delay distribution standard deviation of flight j, σ
j, and a time increment of Δ
t. (31)
p
j
,
k
=
∑
i
=
0
4
∗
σ
j
-
τ
k
-
e
j
/
Δ
t
η
j
τ
k
-
e
j
+
i
+
1
2
∗
Δ
t
∗
1
-
∑
l
=
0
i
η
k
l
+
1
2
∗
Δ
t
∗
Δ
t
∗
Δ
t
.


## 4. Results

This section looks at the performance from the point of view of the robustness by measuring the number of conflicts for a given set of perturbed schedules, for all the approaches presented in Section3: firstly when they are used alone, secondly when combined with the TRS approach, and finally for the combinations with the approaches which consider the flight densities. The comparison made between the results obtained, when applying the different approaches, uses the Mann-Whitney test to establish the statistical significance of the different approaches, presented in the summary result tables as the number of instances in each range of the number of BSSs which can be said to have no statistically significantly higher numbers of conflicts compared to any of the other approaches compared. Regarding an airport, N referring to the number of BSSs available, three ranges of the number of BSSs have been defined, based on the LMAP and UMAP, where the first is for N
< LMAP, the second is for LMAP ≤
N
< UMAP, and the third is for UMAP ≤
N. These are shown within brackets and separated by a comma in the following tables. Any approach achieving full coverage of a range of the number of BSSs is presented in bold font and the approaches with higher numbers for a range, covering the highest number of BSSs in the range compared, are presented in underlined font to assist in the interpretation of the summary result tables.The robustness approaches described are applied to the ABSSAP using the Steady State Evolutionary Algorithm (SSEA) from [3] and their results are compared and analysed in this section using the data sets obtained from NATS for London Heathrow airport Terminal 1, which were also used in [2, 3].To compare the performance of each of the robustness approaches introduced in the previous sections, three sets of perturbed schedules were generated using a folded normal distribution with a zero mean and 10-, 20-, and 30-minute standard deviations. These sets are used to calculate the average number of conflicts for each robustness approach, where a lower value represents a more robust solution compared to those with higher values. Each experiment is repeated at least 30 times. To calculate the number of conflicts within a solution in the perturbed schedule, each flight in the original solution is assigned to the same BSS as in the original solution where possible; otherwise it is assigned to the dummy. This is repeated until all of the flights are assigned to BSS or the dummy. The number of flights assigned to the dummy represents the number of conflicts. If a solution does not achieve the maximum assignment possible, then the number of extra unassigned flights may be accounted for with a higher contribution, given that they are less desirable solutions. Nevertheless, all of the solutions having the highest fitness in the following experiments achieve maximum assignments, which simplifies comparison of the different robustness approaches, based on the described measure. This is possible since the measure only depends on assignments to the BSSs and not on any other objective, such as those upon which the fitness depends.The number of possible perturbed schedules depends on the number of flights, and this accounts for an extremely large number of possible combinations, making it impossible to consider them all. The number of combinations for 194 flights corresponds to194
!
≈
1.3291
∗
1
0
361 which is far greater than 1
0
82 atoms estimated to exist in the observable universe. The number of perturbed schedules necessary to calculate the quality of a solution should therefore be as large as possible to account for as many potential combinations as possible. However, as the number of perturbed schedules increases, so does the time required to perform the calculations, and the memory requirements also increase likewise: moreover, this is further multiplied by the number of solutions which will be used in the comparison. However, not all combinations are likely to represent a valid schedule. So, 10,000 perturbed schedules were used based on the number of solutions to be processed, the time available, and the memory required.The stochastic robustness approach RNC is time consuming when compared with the other approaches reviewed. It would be desirable to use an approach which provides solutions closer to, or better than, those provided by the stochastic robustness approach, without the heavy cost of the time required. With the aim of assessing the difference in performance when the number of perturbed schedules is reduced, two sets of 1,000 and 25 perturbed schedules used with the RNC were also considered. Initial experiments were conducted using SSEA withl
=
1 and the RNC and Multiexchange between a Fixed Number of 3 Resources (MEFNR3) for 25 perturbed schedules and 800,000 total iterations which required an average execution time of 52 min per instance. These, when extrapolated to 1,000 perturbed schedules, provide an execution time of around 34 hours per instance, whereas the other approaches require no more than two minutes to complete the full set of iterations. Both applications of RNC required too long an execution time for the number of iterations specified, so they were run with a time limit of 30 min, as the other approaches required less than 2 min.The objective importance presented in Section2.1 is used in the following experiment where the most important objective is to achieve maximum assignment, with the second in importance being to maximise robustness and the third objective being minimisation of the distance between flights and their assigned BSS being the last objective considered. The fitness function used to guide the search in the SSEA is a weighted sum of the different objectives evaluated, introduced previously in Section 2.1, whose weights are those calculated in [3]. Both approaches UESRS and PCBG need their robustness weight (W
2) to be recalculated. Thus, given a maximum distance between a flight and its assigned BSS, D
m
a
x, which depends on the airport topology, which for the topologies studied here is D
m
a
x
=
9, a distance of one unit is assumed between different sides of a pier and a distance of two units was assumed between different piers. A new assignment between two previously assigned flights may incur a service reduction for the new assignment and next flight, which is used to obtain the decrease in robustness (the second objective), which for UESRS and PCBG cannot be greater than 1 for each flight, thus totalling 2 in this case. Finally, using the objective priorities the following relations can be established: W
1
∗
1
>
|
W
2
|
∗
2
+
|
W
3
|
∗
D
m
a
x and |
W
2
|
∗
2
>
|
W
3
|
∗
D
m
a
x. The original conditions for the weights, when the “minimise reduction in service time” objective was used (Section 2.1), are W
3
=
-
1 and W
1
>
23.4, which together with the objective priorities obtain |
W
2
|
>
(
1
∗
9
)
/
2
=
4.5 and |
W
2
|
∗
2
+
|
W
3
|
∗
9
<
23.4 giving |
W
2
|
<
14.4
/
2
=
7.2. The value used for W
2 is -
7.2 as W
1
=
90
>
23.4. The value used for W
1 is greater than the value originally used to calculate W
2, so a value of −10 was also used.The unsupervised estimation functions introduced in Section3.6 were used for the same parameter values as those used in [18] and with Offset Inverse b
=
6, Offset Inverse b
=
15, Offset Sublinear γ
=
0, and Offset Sublinear γ
=
1000. A summary of the robustness approaches studied and their parameter values is shown in Table 3. In Table 3, the first column contains the name of the robustness approaches considered, all of which were introduced in Section 3; the second column shows the weights for each approach, and the subsequent columns show the name and value of the parameters for the corresponding robustness approach. The SSEA with l
=
1, operator MEFNR3 from [2, 3], a population size of 10, the replacement strategy being Index Selection with Elitist Selection and a group size of 1 (IS1ES), and a maximum of 800,000 iterations for two data sets from London Heathrow airport Terminal 1 (H1T091216: 16th December 2009 and H1T100301: 1st March 2010) were used to obtain the solutions for comparison. The fitness is the sum of the total number of assignments with weight of 90, less the distance between assigned flights and their assigned BSS with a weight of 1, less the robustness approach with the appropriate weight, all of which are shown in Table 3.Table 3
Robustness approaches used with their parameter values.
Approach
Weight (W
2)
Parameters Name: values
TRS
0.008
Buffer time: 30 min. long-haul and 15 min. others
ATRS
0.008
Buffer time: 30 min. long-haul and 15 min. others
ARS
0.008
Buffer time: 30 min. long-haul and 15 min. others
BSARS
0.008
Buffer time: 30 min. long-haul and 15 min. others
SARS
0.008
Buffer time: 30 min. long-haul and 15 min. others
PCBG
7.2 and 10
Std. deviation: 10, 20 and 30 min.
RNC
10 and 14
Std. deviation: 10, 20 and 30 min.
Num. schedules: 25 and 1000
Max. execution time: 30 min.
UESRS
7.2 and 10
Estimation function: Exp 0.03, Exp 0.05, Inverse 6, Inverse 15, Linear, Offset Inverse 6, Offset Inverse 15, Offset Sublinear 0, Offset Sublinear 1000, and Sublinear
Buffer time: 30 min. long-haul and 15 min. othersThe results presented in the following sections were summarised for simplicity and clarity by considering the average number of times an approach achieves statistically significantly lower conflicts, or at least no worse, than the other approaches, in the different regions of numbers of BSSs, based on the LMAP and UMAP, which divide the range of BSSs studied into three areas corresponding toN
<
L
M
A
P, L
M
A
P
≤
N
<
U
M
A
P, and U
M
A
P
≤
N. The values between brackets correspond to the number of times the approach provides significantly statistically solutions not worse than the other approaches used, for each of the BSS ranges. Table 4 shows that the ARS approach for the data set of H1T091216 as well as a 4-pier topology has the values (
2
,
3
,
3
_
) showing that it achieves a statistically significant number of conflicts not higher in two instances for the range of N
< LMAP, three for the range LMAP ≤
N
< UMAP, and three for N
≥ UMAP. Thus, the larger the number between parenthesis, the better the performances in respect of robustness.Table 4
Number of instances with significantly statistically lower number of conflicts in each range of numbers of BSSs for disruptions with a delayσ
=
10 min for ARS, BSARS, SARS + TRS, and TRS for a significance level of 0.05.
Approach
H1T091216
H1T100301
Max. (9, 5, 3)
Max. (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
ARS
(1, 2, 2)
(2, 3,3
_)
(3, 3, 4)
(1, 2, 1)
BSARS
(9
_, 5
_, 3
_)
(8
_, 2, 3
_)
(4
_, 6
_, 5
_)
(4
_, 5
_, 5
_)
SARS + TRS
(3, 4, 2)
(3,5
_, 3
_)
(3, 3, 4)
(0, 2, 1)
TRS
(0, 1, 1)
(1, 0, 0)
(2, 3, 4)
(1, 2, 0)Bold font is used to identify those cases where the robustness approach achieves good results for all numbers of BSSs in a range; for example, in the previous example the approach performs well for all numbers of BSSs in the range ofN
≥
27 (UMAP). Underlining is used to identify those cases where the robustness approach performs well for greater numbers of BSSs in a range. The maximum quantity of instances of number of BSSs in a range is presented between brackets at the top of the table for each of the ranges discussed, preceded by the word “Max.”; for example, in Table 4 the cell in the second column and second row of the header shows that the first range contains nine instances of numbers of BSSs for N
< LMAP, five for LMAP ≤
N
< UMAP, and three for N
≥ UMAP. For simplicity and clarity, those approaches which do achieve statistically significantly higher number of conflicts on average than any of the other approaches in all three BSSs regions are not shown in the tables which followed.The next sections present the experimental results for the robustness approaches presented.


### 4.1. Results for All Robustness Approaches Considered

In this section, we look at the performance from the point of view of robustness for all the approaches presented in Section3; first when they are used alone and then when combined with the TRS approach.The BSARS provides more robust solutions than those obtained by the other proposed approaches, that is, TRS, ARS, SARS, and TRS + SARS, as shown in Table4. Furthermore, BSARS will be seen to perform even better when combined with TRS in Tables 5, 6, and 7. ARS and SARS were similarly seen to improve robustness when combined with TRS.Table 5
Number of instances with a significantly statistically lower number of conflicts in each range of numbers of BSSs for disruptions with a delayσ
=
10 min and all the approaches alone and combined with TRS for a significance level of 0.05.
Approach
H1T091216
H1T100301
Max. (9, 5, 3)
Max. (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
BSARS
(2, 0, 0)
(1, 0, 0)
(0, 0, 0)
(3, 0, 0)
TRS + BSARS
(3, 0, 0)
(3, 0, 0)
(0, 0, 0)
(4, 0, 0)
PCBG
7.2
(5, 2, 1)
(3, 1, 0)
(1, 1, 0)
(5
_, 2, 0)
10
(7
_, 4
_, 0)
(8
_, 1, 0)
(2,4
_, 0)
(5
_, 4
_, 0)
TRS + UESRS
Exp 0.03
7.2
(2, 1, 0)
(0, 1, 0)
(2, 1, 1)
(1, 2, 0)
10
(2, 1, 0)
(1,4
_, 3
_)
(4
_, 2, 4
_)
(2, 1, 1)
Inverse 6
7.2
(0, 0,2
_)
(0, 1, 0)
(0, 1, 1)
(0, 0, 0)
10
(0, 0,2
_)
(0, 1, 2)
(0, 1,4
_)
(0, 2,4
_)Table 6
Number of instances with a significantly statistically lower number of conflicts in each range of numbers of BSSs for disruptions with a delayσ
=
20 min and all the approaches alone and combined with TRS for a significance level of 0.05.
Approach
H1T091216
H1T1003010
Max. (9, 5, 3)
Max. (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
BSARS
(0, 0, 0)
(1, 0, 0)
(0, 0, 0)
(2, 0, 0)
TRS + BSARS
(0, 0, 0)
(4, 0, 0)
(0, 0, 0)
(3, 0, 0)
PCBG
7.2
(5, 1, 0)
(2, 0, 0)
(3, 2, 0)
(3, 0, 0)
10
(7
_, 2, 1)
(8
_, 5
_, 1)
(1,5
_, 1)
(5
_, 6
_, 1)
TRS + UESRS
Exp 0.03
7.2
(2, 1, 0)
(1, 0, 0)
(2, 1, 0)
(1, 0, 0)
10
(5,3
_, 1)
(2, 0,3
_)
(5
_, 1, 4
_)
(4, 0, 1)
Inverse 6
10
(1, 1,3
_)
(0, 2, 2)
(0, 1, 2)
(2, 2,5
_)Table 7
Number of instances with a significantly statistically lower number of conflicts in each range of numbers of BSSs for disruptions with a delayσ
=
30 min and all the approaches alone and combined with TRS for a significance level of 0.05.
Approach
H1T091216
H1T100301
Max. (9, 5, 3)
Max. (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
BSARS
(0, 0, 0)
(1, 0, 0)
(0, 0, 0)
(1, 0, 0)
TRS + BSARS
(0, 0, 0)
(2, 0, 0)
(0, 0, 0)
(2, 0, 0)
PCBG
7.2
(4, 2, 0)
(0, 0, 0)
(1, 4, 0)
(1, 0, 0)
10
(1, 4,3
_)
(5
_, 4
_, 3
_)
(2,6
_, 5
_)
(3
_, 6
_, 5
_)
TRS + UESRS
Exp 0.03
7.2
(3, 0, 0)
(1, 0, 0)
(3, 0, 0)
(1, 0, 0)
10
(9
_, 1, 0)
(4, 1, 0)
(5
_, 1, 0)
(2, 0, 0)
Inverse 6
7.2
(1, 0, 0)
(1, 0, 0)
(1, 0, 0)
(2, 0, 0)
10
(3, 3, 1)
(2, 3, 0)
(0, 0, 0)
(3
_, 0, 0)The results which are summarised in Tables5, 6, and 7 show that the PCBG does not gain any advantage when combined with the TRS, but it appears to be detrimental, as PCBG alone performs better throughout the whole ranges of the numbers of BSSs compared to when used combined with TRS. Even with the RNC for 25 and 1,000 instances taking account of the distribution of delays, the other approaches provide a statistically significantly lower number of collisions, but with a much lower running time. This could be due to the fact that the number of iterations which it is possible to execute in the 30 minutes is too low to find promising solutions with a lower number of collisions compared to those obtained by the other approaches. Around 33,000 to 70,000 iterations were executed for the experiments conducted here (the numbers depending mainly on the data set and the number of BSSs).The PCBG provides statistically significantly lower number of conflicts through a wider range of BSSs, and such a range also includes the range of BSSs used in the real problems; that is,N
≥
U
M
A
P. Nevertheless, this result could be regarded as biased, given that the PCBG considers a normal folded distribution of the same standard deviation as that from which the perturbed schedules were generated. The PCBG could consider different standard deviations and distributions depending on the aircraft, season, route, destination, and time of the day which should further improve the results in real situations. The running time is also of the same magnitude as that for the other approaches, with the exception of the RNC, which has a much higher running time for the same number of iterations.There is a possibility of combining the different approaches, in particular either ARS or BSARS with UESRS studied below. There is also a question about the preference as to how these approaches are combined, either as a sum of each individual with the TRS, as studied in here, or as a product of their individual contributions.


### 4.2. Other Combinations of These Approaches

In this section we study the combination of the UESRS with each of the ARS, BSARS, and TRS. It may be noted that the SARS is equal to BSARS divided by (UMAP-LMAP), which it is equivalent to thus changing the wight (W
2) in the fitness function. As seen previously in Table 4, BSARS provides more robust solutions compared to when using SARS alone for the same weight (W
2) or when combined with TRS, and these differences between both similar approaches are even more pronounced when BSARS is used in combination with TRS.The combination of TRS and UESRS with exponential estimation function overall performs better forN
< LMAP, as shown in Tables 8, 9, and 10, which correspond to the situation when there are no sufficient resources for the number of flights to service. The results which are summarised in Tables 11, 12, and 13 and which correspond to the cases when two robustness approaches are combined show that the combination of ARS with UESRS for exponential estimation function overall performs better for LMAP ≤
N
< UMAP, which refers to the situation when there are just sufficient resources for the number of flights to service, as, in this case, it is particularly important to consider the distribution of the flights requiring service which ARS takes into account. Also the combination of UESRS with BSARS in general improves the performance more than when using ARS for all the ranges. Similarly, the combination of BSARS with UESRS with Inverse estimation function overall performs better for UMAP ≤
N and also for N lower than UMAP but closer to it. This would be the preferable approach for real problems. Finally, the ARS with UESRS with Inverse estimation function performance decreases for big delays for N
≥ UMAP and improves for N
≤ LMAP.Table 8
Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of0 min average and σ
=
10 min deviation, Mann-Whitney test, and significance level of 0.05 and for the TRS, UESRS, and the combination of TRS with UESRS.
Approach
H091216
H100301
17 (9, 5, 3)
17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008 × TRS + 10 × UESRS (Exp 0.03)
14
_ (9
_, 4
_, 1)
15
_ (7
_, 5
_, 3
_)
14
_ (6
_, 5
_, 3
_)
11
_ (6
_, 4
_, 1)
0.008 × TRS + 10 × UESRS (Inverse 2)
1 (0, 0, 1)
2 (0, 2, 0)
1 (0, 0, 1)
0 (0, 0, 0)
0.008 × TRS + 10 × UESRS (Inverse 4)
2 (0, 0,2
_)
3 (0, 2, 1)
3 (0, 1, 2)
6 (0, 2,4
_)
0.008 × TRS + 10 × UESRS (Inverse 6)
3 (0, 1,2
_)
4 (0, 2, 2)
4 (0, 1,3
_)
6 (0, 2,4
_)
0.008 × TRS + 7.2 × UESRS (Exp 0.03)
8 (7, 1, 0)
5 (4, 1, 0)
9 (4, 4, 1)
8 (4,4
_, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 2)
1 (0, 0, 1)
1 (0, 1, 0)
0 (0, 0, 0)
0 (0, 0, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 4)
1 (0, 0, 1)
1 (0, 1, 0)
2 (0, 1, 1)
1 (0, 1, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 6)
2 (0, 0,2
_)
1 (0, 1, 0)
2 (0, 1, 1)
0 (0, 0, 0)
10 × UESRS (Exp 0.03)
0 (0, 0, 0)
1 (1, 0, 0)
2 (0, 0, 2)
0 (0, 0, 0)Table 9
Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of 0 min average andσ
=
20 min deviation, Mann-Whitney test, and significance level of 0.05 and for the TRS, UESRS, and the combination of TRS with UESRS.
Approach
H091216
H100301
Max. 17 (9, 5, 3)
Max. 17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008× TRS + 10 × UESRS (Exp 0.03)
12
_ (9
_, 3
_, 0)
13
_ (7
_, 3, 3
_)
14
_ (6
_, 4
_, 4
_)
9 (5
_, 3, 1)
0.008 × TRS + 10 × UESRS (Inverse 2)
2 (2, 0, 0)
4 (4, 0, 0)
0 (0, 0, 0)
2 (2, 0, 0)
0.008 × TRS + 10 × UESRS (Inverse 4)
5 (2, 1, 2)
5 (2, 3, 0)
1 (0, 0, 1)
5 (2, 1, 2)
0.008 × TRS + 10 × UESRS (Inverse 6)
5 (0, 2,3)
10 (3,5
_, 2)
4 (0, 1, 3)
12
_ (3, 4
_, 5
_)
0.008 × TRS + 7.2 × UESRS (Exp 0.03)
4 (3, 1, 0)
3 (3, 0, 0)
4 (3, 1, 0)
4 (2, 2, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 2)
1 (0, 1, 0)
1 (1, 0, 0)
0 (0, 0, 0)
1 (1, 0, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 4)
1 (1, 0, 0)
1 (1, 0, 0)
1 (0, 1, 0)
0 (0, 0, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 6)
0 (0, 0, 0)
0 (0, 0, 0)
1 (0, 1, 0)
1 (0, 1, 0)
10 × UESRS (Exp 0.03)
0 (0, 0, 0)
1 (1, 0, 0)
1 (0, 0, 1)
0 (0, 0, 0)Table 10
Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of0 min average and σ
=
30 min deviation, Mann-Whitney test, and significance level of 0.05 and for the TRS, UESRS, and the combination of TRS with UESRS.
Approach
H091216
H100301
Max. 17 (9, 5, 3)
Max. 17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008 × TRS + 10 × UESRS (Exp 0.03)
13
_ (9
_, 3
_, 1)
12
_ (6
_, 3, 3
_)
15
_ (6
_, 5
_, 4
_)
6 (3, 2, 1)
0.008×TRS + 10×UESRS (Inverse 2)
5 (5, 0, 0)
3 (3, 0, 0)
0 (0, 0, 0)
4 (4
_, 0, 0)
0.008×TRS + 10×UESRS (Inverse 4)
7 (5, 1, 1)
8 (5, 3, 0)
3 (1, 1, 1)
3 (3, 0, 0)
0.008 × TRS + 10 × UESRS (Inverse 6)
9 (3,3
_, 3
_)
11 (5,5
_, 1)
10 (2,5
_, 3)
12
_ (3, 4
_, 5
_)
0.008×TRS + 7.2×UESRS (Exp 0.03)
3 (3, 0, 0)
2 (2, 0, 0)
3 (3, 0, 0)
2 (1, 1, 0)
0.008×TRS + 7.2×UESRS (Inverse 2)
1 (1, 0, 0)
3 (3, 0, 0)
0 (0, 0, 0)
2 (2, 0, 0)
0.008×TRS + 7.2×UESRS (Inverse 4)
3 (3, 0, 0)
2 (2, 0, 0)
1 (0, 1, 0)
1 (1, 0, 0)
0.008×TRS + 7.2×UESRS (Inverse 6)
1 (1, 0, 0)
2 (2, 0, 0)
3 (2, 1, 0)
2 (2, 0, 0)
10×UESRS (Exp 0.03)
0 (0, 0, 0)
0 (0, 0, 0)
2 (1, 0, 1)
0 (0, 0, 0)Table 11
Conflicts statistical significance for combined robustness approaches for Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of0 min average and σ
=
10 min deviation, Mann-Whitney test, and significance level of 0.05.
Approach
H091216
H100301
Max. 17 (9, 5, 3)
Max. 17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008 × ARS + 10 × UESRS (Exp 0.03)
3 (1, 2, 0)
5 (0,3
_, 2)
4 (0, 1, 3)
3 (0, 2, 1)
0.008 × ARS + 10 × UESRS (Inverse 6)
1 (0, 0, 1)
4 (0, 1,3
_)
3 (0, 1, 2)
6
_ (0, 1, 5
_)
0.008 × ARS + 7.2 × UESRS (Inverse 6)
3 (0, 2, 1)
0 (0, 0, 0)
1 (0, 0, 1)
0 (0, 0, 0)
0.008 × BSARS + 10 × UESRS (Exp 0.03)
7
_ (4, 2, 1)
6
_ (0, 3
_, 3
_)
3 (0, 0, 3)
6
_ (3, 2, 1)
0.008 × BSARS + 10× UESRS (Inverse 6)
7
_ (1, 3
_, 3
_)
5 (0, 2,3
_)
6
_ (0, 2
_, 4
_)
6
_ (0, 2, 4
_)
0.008 × BSARS + 7.2 × UESRS (Exp 0.03)
3 (2, 1, 0)
2 (2, 0, 0)
1 (0, 0, 1)
5 (4
_, 1, 0)
0.008 × BSARS + 7.2 × UESRS (Inverse 6)
5 (1, 2, 2)
0 (0, 0, 0)
2 (0, 1, 1)
0 (0, 0, 0)
0.008 × TRS + 0.008 × BSARS
0 (0, 0, 0)
2 (2, 0, 0)
0 (0, 0, 0)
1 (1, 0, 0)
0.008 × TRS + 10 × UESRS (Exp 0.03)
0 (0, 0, 0)
2 (0, 0, 2)
4 (2
_, 0, 2)
1 (0, 0, 1)
0.008 × TRS + 10 × UESRS (Inverse 2)
1 (0, 0, 1)
0 (0, 0, 0)
1 (0, 0, 1)
0 (0, 0, 0)
0.008 × TRS + 10 × UESRS (Inverse 4)
1 (0, 0, 1)
0 (0, 0, 0)
2 (0, 0, 2)
4 (0, 1, 3)
0.008 × TRS + 10 × UESRS (Inverse 6)
2 (0, 0, 2)
1 (0, 0, 1)
3 (0, 0, 3)
5 (0, 1, 4)
0.008 × TRS + 7.2 × PCBG (10)
6 (6
_, 0, 0)
2 (2, 0, 0)
4 (2
_, 2, 0)
4 (2, 2, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 4)
1 (0, 0, 1)
0 (0, 0, 0)
1 (0, 0, 1)
0 (0, 0, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 6)
1 (0, 0, 1)
0 (0, 0, 0)
1 (0, 0, 1)
0 (0, 0, 0)
10 × PCBG (10)
5 (4, 1, 0)
5 (5
_, 0, 0)
6 (2
_, 4
_, 0)
5 (2,3
_, 0)
7.2 × PCBG (10)
0 (0, 0, 0)
0 (0, 0, 0)
1 (0, 1, 0)
3 (1, 2, 0)Table 12
Conflicts statistical significance for combined robustness approaches for Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of0 min average and σ
=
20 min deviation, Mann-Whitney test, and significance level of 0.05.
Approach
H091216
H100301
Max. 17 (9, 5, 3)
Max. 17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008 × ARS + 10 × UESRS (Exp 0.03)
3 (1, 2, 0)
0 (0, 0, 0)
3 (0, 1, 2)
1 (0, 0, 1)
0.008 × ARS + 10 × UESRS (Inverse 6)
3 (0, 1, 2)
0 (0, 0, 0)
3 (0, 0, 3)
7 (0, 2,5
_)
0.008 × BSARS + 10 × UESRS (Exp 0.03)
7 (3, 2, 2)
0 (0, 0, 0)
5 (0, 1,4
_)
3 (3, 0, 0)
0.008 × BSARS + 10 × UESRS (Inverse 6)
10
_ (3, 4
_, 3
_)
0 (0, 0, 0)
6 (0, 3, 3)
9
_ (2, 3, 4)
0.008 × BSARS + 7.2 × UESRS (Exp 0.03)
3 (2, 0, 1)
0 (0, 0, 0)
0 (0, 0, 0)
3 (3
_, 0, 0)
0.008 × BSARS + 7.2 × UESRS (Inverse 6)
2 (2, 0, 0)
0 (0, 0, 0)
0 (0, 0, 0)
1 (1, 0, 0)
0.008 × TRS + 0.008 × ARS
0 (0, 0, 0)
17
_ (9
_, 5
_, 3
_)
0 (0, 0, 0)
0 (0, 0, 0)
0.008 × TRS + 10 × UESRS (Exp 0.03)
0 (0, 0, 0)
0 (0, 0, 0)
2 (1, 0, 1)
1 (0, 0, 1)
0.008 × TRS + 10 × UESRS (Inverse 6)
1 (0, 0, 1)
0 (0, 0, 0)
2 (0, 0, 2)
5 (0, 1, 4)
0.008 × TRS + 7.2 × PCBG (20)
7 (7
_, 0, 0)
0 (0, 0, 0)
8
_ (4
_, 4
_, 0)
3 (2, 1, 0)
10 × PCBG (20)
4 (3, 1, 0)
0 (0, 0, 0)
6 (1,4
_, 1)
9
_ (2, 6
_, 1)Table 13
Conflicts statistical significance for combined robustness approaches for Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of0 min average and σ
=
30 min deviation, Mann-Whitney test, and significance level of 0.05.
Approach
H091216
H100301
Max. 17 (9, 5, 3)
Max. 17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008 × ARS + 10 × UESRS (Exp 0.03)
3 (2, 1, 0)
1 (0, 1, 0)
3 (2, 1, 0)
0 (0, 0, 0)
0.008 × ARS + 10 × UESRS (Inverse 6)
2 (0, 1, 1)
1 (0, 1, 0)
0 (0, 0, 0)
1 (1, 0, 0)
0.008 × ARS + 7.2 × UESRS (Exp 0.03)
0 (0, 0, 0)
1 (1, 0, 0)
1 (1, 0, 0)
0 (0, 0, 0)
0.008 × ARS + 7.2 × UESRS (Inverse 6)
0 (0, 0, 0)
0 (0, 0, 0)
1 (1, 0, 0)
1 (1, 0, 0)
0.008 × BSARS + 10 × UESRS (Exp 0.03)
6 (3, 2, 1)
0 (0, 0, 0)
4 (2, 1, 1)
4 (3
_, 1, 0)
0.008 × BSARS + 10 × UESRS (Inverse 6)
11 (6, 4, 1)
1 (0, 1, 0)
3 (1, 2, 0)
5 (3
_, 2, 0)
0.008 × BSARS + 7.2 × UESRS (Exp 0.03)
2 (2, 0, 0)
4 (4, 0, 0)
1 (1, 0, 0)
3 (3
_, 0, 0)
0.008 × BSARS + 7.2 × UESRS (Inverse 6)
2 (2, 0, 0)
5 (5
_, 0, 0)
1 (1, 0, 0)
2 (2, 0, 0)
0.008 × TRS + 10 × UESRS (Exp 0.03)
1 (1, 0, 0)
1 (1, 0, 0)
3 (3, 0, 0)
1 (1, 0, 0)
0.008 × TRS + 10 × UESRS (Inverse 6)
3 (1, 1, 1)
1 (1, 0, 0)
0 (0, 0, 0)
0 (0, 0, 0)
0.008 × TRS + 7.2 × PCBG (30)
13
_ (7
_, 5
_, 1)
5 (3, 2, 0)
11
_ (5
_, 6
_, 0)
5 (2, 3, 0)
0.008 × TRS + 7.2 × UESRS (Exp 0.03)
2 (2, 0, 0)
1 (1, 0, 0)
0 (0, 0, 0)
0 (0, 0, 0)
10 × PCBG (30)
5 (0, 2,3
_)
11
_ (4, 4
_, 3
_)
9 (2, 2,5
_)
13
_ (2, 6
_, 5
_)Future studies could look at the dependency of the robustness objective weight on performance. If other robustness weights are to be considered, all of the weights taking part in the fitness function should be modified accordingly to maintain the order of importance of each objective. When the unassigned flights are not taken into account by the robustness measure used in the fitness function, then particular care has to be taken when selecting the appropriate weights for maximisation of the number of assignments (the first objective), since incorrect selection of this objective weight may sometimes deem solutions with a lower number of assignments to be fitter. This interdependency mainly between the maximum number of assignments’ objective and the robustness could also be decoupled by penalising the unassigned flights in the robustness.


## 4.1. Results for All Robustness Approaches Considered

In this section, we look at the performance from the point of view of robustness for all the approaches presented in Section3; first when they are used alone and then when combined with the TRS approach.The BSARS provides more robust solutions than those obtained by the other proposed approaches, that is, TRS, ARS, SARS, and TRS + SARS, as shown in Table4. Furthermore, BSARS will be seen to perform even better when combined with TRS in Tables 5, 6, and 7. ARS and SARS were similarly seen to improve robustness when combined with TRS.Table 5
Number of instances with a significantly statistically lower number of conflicts in each range of numbers of BSSs for disruptions with a delayσ
=
10 min and all the approaches alone and combined with TRS for a significance level of 0.05.
Approach
H1T091216
H1T100301
Max. (9, 5, 3)
Max. (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
BSARS
(2, 0, 0)
(1, 0, 0)
(0, 0, 0)
(3, 0, 0)
TRS + BSARS
(3, 0, 0)
(3, 0, 0)
(0, 0, 0)
(4, 0, 0)
PCBG
7.2
(5, 2, 1)
(3, 1, 0)
(1, 1, 0)
(5
_, 2, 0)
10
(7
_, 4
_, 0)
(8
_, 1, 0)
(2,4
_, 0)
(5
_, 4
_, 0)
TRS + UESRS
Exp 0.03
7.2
(2, 1, 0)
(0, 1, 0)
(2, 1, 1)
(1, 2, 0)
10
(2, 1, 0)
(1,4
_, 3
_)
(4
_, 2, 4
_)
(2, 1, 1)
Inverse 6
7.2
(0, 0,2
_)
(0, 1, 0)
(0, 1, 1)
(0, 0, 0)
10
(0, 0,2
_)
(0, 1, 2)
(0, 1,4
_)
(0, 2,4
_)Table 6
Number of instances with a significantly statistically lower number of conflicts in each range of numbers of BSSs for disruptions with a delayσ
=
20 min and all the approaches alone and combined with TRS for a significance level of 0.05.
Approach
H1T091216
H1T1003010
Max. (9, 5, 3)
Max. (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
BSARS
(0, 0, 0)
(1, 0, 0)
(0, 0, 0)
(2, 0, 0)
TRS + BSARS
(0, 0, 0)
(4, 0, 0)
(0, 0, 0)
(3, 0, 0)
PCBG
7.2
(5, 1, 0)
(2, 0, 0)
(3, 2, 0)
(3, 0, 0)
10
(7
_, 2, 1)
(8
_, 5
_, 1)
(1,5
_, 1)
(5
_, 6
_, 1)
TRS + UESRS
Exp 0.03
7.2
(2, 1, 0)
(1, 0, 0)
(2, 1, 0)
(1, 0, 0)
10
(5,3
_, 1)
(2, 0,3
_)
(5
_, 1, 4
_)
(4, 0, 1)
Inverse 6
10
(1, 1,3
_)
(0, 2, 2)
(0, 1, 2)
(2, 2,5
_)Table 7
Number of instances with a significantly statistically lower number of conflicts in each range of numbers of BSSs for disruptions with a delayσ
=
30 min and all the approaches alone and combined with TRS for a significance level of 0.05.
Approach
H1T091216
H1T100301
Max. (9, 5, 3)
Max. (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
BSARS
(0, 0, 0)
(1, 0, 0)
(0, 0, 0)
(1, 0, 0)
TRS + BSARS
(0, 0, 0)
(2, 0, 0)
(0, 0, 0)
(2, 0, 0)
PCBG
7.2
(4, 2, 0)
(0, 0, 0)
(1, 4, 0)
(1, 0, 0)
10
(1, 4,3
_)
(5
_, 4
_, 3
_)
(2,6
_, 5
_)
(3
_, 6
_, 5
_)
TRS + UESRS
Exp 0.03
7.2
(3, 0, 0)
(1, 0, 0)
(3, 0, 0)
(1, 0, 0)
10
(9
_, 1, 0)
(4, 1, 0)
(5
_, 1, 0)
(2, 0, 0)
Inverse 6
7.2
(1, 0, 0)
(1, 0, 0)
(1, 0, 0)
(2, 0, 0)
10
(3, 3, 1)
(2, 3, 0)
(0, 0, 0)
(3
_, 0, 0)The results which are summarised in Tables5, 6, and 7 show that the PCBG does not gain any advantage when combined with the TRS, but it appears to be detrimental, as PCBG alone performs better throughout the whole ranges of the numbers of BSSs compared to when used combined with TRS. Even with the RNC for 25 and 1,000 instances taking account of the distribution of delays, the other approaches provide a statistically significantly lower number of collisions, but with a much lower running time. This could be due to the fact that the number of iterations which it is possible to execute in the 30 minutes is too low to find promising solutions with a lower number of collisions compared to those obtained by the other approaches. Around 33,000 to 70,000 iterations were executed for the experiments conducted here (the numbers depending mainly on the data set and the number of BSSs).The PCBG provides statistically significantly lower number of conflicts through a wider range of BSSs, and such a range also includes the range of BSSs used in the real problems; that is,N
≥
U
M
A
P. Nevertheless, this result could be regarded as biased, given that the PCBG considers a normal folded distribution of the same standard deviation as that from which the perturbed schedules were generated. The PCBG could consider different standard deviations and distributions depending on the aircraft, season, route, destination, and time of the day which should further improve the results in real situations. The running time is also of the same magnitude as that for the other approaches, with the exception of the RNC, which has a much higher running time for the same number of iterations.There is a possibility of combining the different approaches, in particular either ARS or BSARS with UESRS studied below. There is also a question about the preference as to how these approaches are combined, either as a sum of each individual with the TRS, as studied in here, or as a product of their individual contributions.


## 4.2. Other Combinations of These Approaches

In this section we study the combination of the UESRS with each of the ARS, BSARS, and TRS. It may be noted that the SARS is equal to BSARS divided by (UMAP-LMAP), which it is equivalent to thus changing the wight (W
2) in the fitness function. As seen previously in Table 4, BSARS provides more robust solutions compared to when using SARS alone for the same weight (W
2) or when combined with TRS, and these differences between both similar approaches are even more pronounced when BSARS is used in combination with TRS.The combination of TRS and UESRS with exponential estimation function overall performs better forN
< LMAP, as shown in Tables 8, 9, and 10, which correspond to the situation when there are no sufficient resources for the number of flights to service. The results which are summarised in Tables 11, 12, and 13 and which correspond to the cases when two robustness approaches are combined show that the combination of ARS with UESRS for exponential estimation function overall performs better for LMAP ≤
N
< UMAP, which refers to the situation when there are just sufficient resources for the number of flights to service, as, in this case, it is particularly important to consider the distribution of the flights requiring service which ARS takes into account. Also the combination of UESRS with BSARS in general improves the performance more than when using ARS for all the ranges. Similarly, the combination of BSARS with UESRS with Inverse estimation function overall performs better for UMAP ≤
N and also for N lower than UMAP but closer to it. This would be the preferable approach for real problems. Finally, the ARS with UESRS with Inverse estimation function performance decreases for big delays for N
≥ UMAP and improves for N
≤ LMAP.Table 8
Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of0 min average and σ
=
10 min deviation, Mann-Whitney test, and significance level of 0.05 and for the TRS, UESRS, and the combination of TRS with UESRS.
Approach
H091216
H100301
17 (9, 5, 3)
17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008 × TRS + 10 × UESRS (Exp 0.03)
14
_ (9
_, 4
_, 1)
15
_ (7
_, 5
_, 3
_)
14
_ (6
_, 5
_, 3
_)
11
_ (6
_, 4
_, 1)
0.008 × TRS + 10 × UESRS (Inverse 2)
1 (0, 0, 1)
2 (0, 2, 0)
1 (0, 0, 1)
0 (0, 0, 0)
0.008 × TRS + 10 × UESRS (Inverse 4)
2 (0, 0,2
_)
3 (0, 2, 1)
3 (0, 1, 2)
6 (0, 2,4
_)
0.008 × TRS + 10 × UESRS (Inverse 6)
3 (0, 1,2
_)
4 (0, 2, 2)
4 (0, 1,3
_)
6 (0, 2,4
_)
0.008 × TRS + 7.2 × UESRS (Exp 0.03)
8 (7, 1, 0)
5 (4, 1, 0)
9 (4, 4, 1)
8 (4,4
_, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 2)
1 (0, 0, 1)
1 (0, 1, 0)
0 (0, 0, 0)
0 (0, 0, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 4)
1 (0, 0, 1)
1 (0, 1, 0)
2 (0, 1, 1)
1 (0, 1, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 6)
2 (0, 0,2
_)
1 (0, 1, 0)
2 (0, 1, 1)
0 (0, 0, 0)
10 × UESRS (Exp 0.03)
0 (0, 0, 0)
1 (1, 0, 0)
2 (0, 0, 2)
0 (0, 0, 0)Table 9
Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of 0 min average andσ
=
20 min deviation, Mann-Whitney test, and significance level of 0.05 and for the TRS, UESRS, and the combination of TRS with UESRS.
Approach
H091216
H100301
Max. 17 (9, 5, 3)
Max. 17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008× TRS + 10 × UESRS (Exp 0.03)
12
_ (9
_, 3
_, 0)
13
_ (7
_, 3, 3
_)
14
_ (6
_, 4
_, 4
_)
9 (5
_, 3, 1)
0.008 × TRS + 10 × UESRS (Inverse 2)
2 (2, 0, 0)
4 (4, 0, 0)
0 (0, 0, 0)
2 (2, 0, 0)
0.008 × TRS + 10 × UESRS (Inverse 4)
5 (2, 1, 2)
5 (2, 3, 0)
1 (0, 0, 1)
5 (2, 1, 2)
0.008 × TRS + 10 × UESRS (Inverse 6)
5 (0, 2,3)
10 (3,5
_, 2)
4 (0, 1, 3)
12
_ (3, 4
_, 5
_)
0.008 × TRS + 7.2 × UESRS (Exp 0.03)
4 (3, 1, 0)
3 (3, 0, 0)
4 (3, 1, 0)
4 (2, 2, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 2)
1 (0, 1, 0)
1 (1, 0, 0)
0 (0, 0, 0)
1 (1, 0, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 4)
1 (1, 0, 0)
1 (1, 0, 0)
1 (0, 1, 0)
0 (0, 0, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 6)
0 (0, 0, 0)
0 (0, 0, 0)
1 (0, 1, 0)
1 (0, 1, 0)
10 × UESRS (Exp 0.03)
0 (0, 0, 0)
1 (1, 0, 0)
1 (0, 0, 1)
0 (0, 0, 0)Table 10
Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of0 min average and σ
=
30 min deviation, Mann-Whitney test, and significance level of 0.05 and for the TRS, UESRS, and the combination of TRS with UESRS.
Approach
H091216
H100301
Max. 17 (9, 5, 3)
Max. 17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008 × TRS + 10 × UESRS (Exp 0.03)
13
_ (9
_, 3
_, 1)
12
_ (6
_, 3, 3
_)
15
_ (6
_, 5
_, 4
_)
6 (3, 2, 1)
0.008×TRS + 10×UESRS (Inverse 2)
5 (5, 0, 0)
3 (3, 0, 0)
0 (0, 0, 0)
4 (4
_, 0, 0)
0.008×TRS + 10×UESRS (Inverse 4)
7 (5, 1, 1)
8 (5, 3, 0)
3 (1, 1, 1)
3 (3, 0, 0)
0.008 × TRS + 10 × UESRS (Inverse 6)
9 (3,3
_, 3
_)
11 (5,5
_, 1)
10 (2,5
_, 3)
12
_ (3, 4
_, 5
_)
0.008×TRS + 7.2×UESRS (Exp 0.03)
3 (3, 0, 0)
2 (2, 0, 0)
3 (3, 0, 0)
2 (1, 1, 0)
0.008×TRS + 7.2×UESRS (Inverse 2)
1 (1, 0, 0)
3 (3, 0, 0)
0 (0, 0, 0)
2 (2, 0, 0)
0.008×TRS + 7.2×UESRS (Inverse 4)
3 (3, 0, 0)
2 (2, 0, 0)
1 (0, 1, 0)
1 (1, 0, 0)
0.008×TRS + 7.2×UESRS (Inverse 6)
1 (1, 0, 0)
2 (2, 0, 0)
3 (2, 1, 0)
2 (2, 0, 0)
10×UESRS (Exp 0.03)
0 (0, 0, 0)
0 (0, 0, 0)
2 (1, 0, 1)
0 (0, 0, 0)Table 11
Conflicts statistical significance for combined robustness approaches for Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of0 min average and σ
=
10 min deviation, Mann-Whitney test, and significance level of 0.05.
Approach
H091216
H100301
Max. 17 (9, 5, 3)
Max. 17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008 × ARS + 10 × UESRS (Exp 0.03)
3 (1, 2, 0)
5 (0,3
_, 2)
4 (0, 1, 3)
3 (0, 2, 1)
0.008 × ARS + 10 × UESRS (Inverse 6)
1 (0, 0, 1)
4 (0, 1,3
_)
3 (0, 1, 2)
6
_ (0, 1, 5
_)
0.008 × ARS + 7.2 × UESRS (Inverse 6)
3 (0, 2, 1)
0 (0, 0, 0)
1 (0, 0, 1)
0 (0, 0, 0)
0.008 × BSARS + 10 × UESRS (Exp 0.03)
7
_ (4, 2, 1)
6
_ (0, 3
_, 3
_)
3 (0, 0, 3)
6
_ (3, 2, 1)
0.008 × BSARS + 10× UESRS (Inverse 6)
7
_ (1, 3
_, 3
_)
5 (0, 2,3
_)
6
_ (0, 2
_, 4
_)
6
_ (0, 2, 4
_)
0.008 × BSARS + 7.2 × UESRS (Exp 0.03)
3 (2, 1, 0)
2 (2, 0, 0)
1 (0, 0, 1)
5 (4
_, 1, 0)
0.008 × BSARS + 7.2 × UESRS (Inverse 6)
5 (1, 2, 2)
0 (0, 0, 0)
2 (0, 1, 1)
0 (0, 0, 0)
0.008 × TRS + 0.008 × BSARS
0 (0, 0, 0)
2 (2, 0, 0)
0 (0, 0, 0)
1 (1, 0, 0)
0.008 × TRS + 10 × UESRS (Exp 0.03)
0 (0, 0, 0)
2 (0, 0, 2)
4 (2
_, 0, 2)
1 (0, 0, 1)
0.008 × TRS + 10 × UESRS (Inverse 2)
1 (0, 0, 1)
0 (0, 0, 0)
1 (0, 0, 1)
0 (0, 0, 0)
0.008 × TRS + 10 × UESRS (Inverse 4)
1 (0, 0, 1)
0 (0, 0, 0)
2 (0, 0, 2)
4 (0, 1, 3)
0.008 × TRS + 10 × UESRS (Inverse 6)
2 (0, 0, 2)
1 (0, 0, 1)
3 (0, 0, 3)
5 (0, 1, 4)
0.008 × TRS + 7.2 × PCBG (10)
6 (6
_, 0, 0)
2 (2, 0, 0)
4 (2
_, 2, 0)
4 (2, 2, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 4)
1 (0, 0, 1)
0 (0, 0, 0)
1 (0, 0, 1)
0 (0, 0, 0)
0.008 × TRS + 7.2 × UESRS (Inverse 6)
1 (0, 0, 1)
0 (0, 0, 0)
1 (0, 0, 1)
0 (0, 0, 0)
10 × PCBG (10)
5 (4, 1, 0)
5 (5
_, 0, 0)
6 (2
_, 4
_, 0)
5 (2,3
_, 0)
7.2 × PCBG (10)
0 (0, 0, 0)
0 (0, 0, 0)
1 (0, 1, 0)
3 (1, 2, 0)Table 12
Conflicts statistical significance for combined robustness approaches for Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of0 min average and σ
=
20 min deviation, Mann-Whitney test, and significance level of 0.05.
Approach
H091216
H100301
Max. 17 (9, 5, 3)
Max. 17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008 × ARS + 10 × UESRS (Exp 0.03)
3 (1, 2, 0)
0 (0, 0, 0)
3 (0, 1, 2)
1 (0, 0, 1)
0.008 × ARS + 10 × UESRS (Inverse 6)
3 (0, 1, 2)
0 (0, 0, 0)
3 (0, 0, 3)
7 (0, 2,5
_)
0.008 × BSARS + 10 × UESRS (Exp 0.03)
7 (3, 2, 2)
0 (0, 0, 0)
5 (0, 1,4
_)
3 (3, 0, 0)
0.008 × BSARS + 10 × UESRS (Inverse 6)
10
_ (3, 4
_, 3
_)
0 (0, 0, 0)
6 (0, 3, 3)
9
_ (2, 3, 4)
0.008 × BSARS + 7.2 × UESRS (Exp 0.03)
3 (2, 0, 1)
0 (0, 0, 0)
0 (0, 0, 0)
3 (3
_, 0, 0)
0.008 × BSARS + 7.2 × UESRS (Inverse 6)
2 (2, 0, 0)
0 (0, 0, 0)
0 (0, 0, 0)
1 (1, 0, 0)
0.008 × TRS + 0.008 × ARS
0 (0, 0, 0)
17
_ (9
_, 5
_, 3
_)
0 (0, 0, 0)
0 (0, 0, 0)
0.008 × TRS + 10 × UESRS (Exp 0.03)
0 (0, 0, 0)
0 (0, 0, 0)
2 (1, 0, 1)
1 (0, 0, 1)
0.008 × TRS + 10 × UESRS (Inverse 6)
1 (0, 0, 1)
0 (0, 0, 0)
2 (0, 0, 2)
5 (0, 1, 4)
0.008 × TRS + 7.2 × PCBG (20)
7 (7
_, 0, 0)
0 (0, 0, 0)
8
_ (4
_, 4
_, 0)
3 (2, 1, 0)
10 × PCBG (20)
4 (3, 1, 0)
0 (0, 0, 0)
6 (1,4
_, 1)
9
_ (2, 6
_, 1)Table 13
Conflicts statistical significance for combined robustness approaches for Terminal 1 with algorithm SSEA1, operator MEFNR3, and population size of 2,000 for perturbations of0 min average and σ
=
30 min deviation, Mann-Whitney test, and significance level of 0.05.
Approach
H091216
H100301
Max. 17 (9, 5, 3)
Max. 17 (6, 6, 5)
3-pier
4-pier
3-pier
4-pier
0.008 × ARS + 10 × UESRS (Exp 0.03)
3 (2, 1, 0)
1 (0, 1, 0)
3 (2, 1, 0)
0 (0, 0, 0)
0.008 × ARS + 10 × UESRS (Inverse 6)
2 (0, 1, 1)
1 (0, 1, 0)
0 (0, 0, 0)
1 (1, 0, 0)
0.008 × ARS + 7.2 × UESRS (Exp 0.03)
0 (0, 0, 0)
1 (1, 0, 0)
1 (1, 0, 0)
0 (0, 0, 0)
0.008 × ARS + 7.2 × UESRS (Inverse 6)
0 (0, 0, 0)
0 (0, 0, 0)
1 (1, 0, 0)
1 (1, 0, 0)
0.008 × BSARS + 10 × UESRS (Exp 0.03)
6 (3, 2, 1)
0 (0, 0, 0)
4 (2, 1, 1)
4 (3
_, 1, 0)
0.008 × BSARS + 10 × UESRS (Inverse 6)
11 (6, 4, 1)
1 (0, 1, 0)
3 (1, 2, 0)
5 (3
_, 2, 0)
0.008 × BSARS + 7.2 × UESRS (Exp 0.03)
2 (2, 0, 0)
4 (4, 0, 0)
1 (1, 0, 0)
3 (3
_, 0, 0)
0.008 × BSARS + 7.2 × UESRS (Inverse 6)
2 (2, 0, 0)
5 (5
_, 0, 0)
1 (1, 0, 0)
2 (2, 0, 0)
0.008 × TRS + 10 × UESRS (Exp 0.03)
1 (1, 0, 0)
1 (1, 0, 0)
3 (3, 0, 0)
1 (1, 0, 0)
0.008 × TRS + 10 × UESRS (Inverse 6)
3 (1, 1, 1)
1 (1, 0, 0)
0 (0, 0, 0)
0 (0, 0, 0)
0.008 × TRS + 7.2 × PCBG (30)
13
_ (7
_, 5
_, 1)
5 (3, 2, 0)
11
_ (5
_, 6
_, 0)
5 (2, 3, 0)
0.008 × TRS + 7.2 × UESRS (Exp 0.03)
2 (2, 0, 0)
1 (1, 0, 0)
0 (0, 0, 0)
0 (0, 0, 0)
10 × PCBG (30)
5 (0, 2,3
_)
11
_ (4, 4
_, 3
_)
9 (2, 2,5
_)
13
_ (2, 6
_, 5
_)Future studies could look at the dependency of the robustness objective weight on performance. If other robustness weights are to be considered, all of the weights taking part in the fitness function should be modified accordingly to maintain the order of importance of each objective. When the unassigned flights are not taken into account by the robustness measure used in the fitness function, then particular care has to be taken when selecting the appropriate weights for maximisation of the number of assignments (the first objective), since incorrect selection of this objective weight may sometimes deem solutions with a lower number of assignments to be fitter. This interdependency mainly between the maximum number of assignments’ objective and the robustness could also be decoupled by penalising the unassigned flights in the robustness.


## 5. Conclusions

Several approaches were presented here to take account of solution robustness when applying the SSEA. The TRS approach provides solutions with a statistically significantly higher number of conflicts than those obtained by many of the other approaches considered, as shown in [21]. TRS does not consider the extra increase in conflict, as the service time is reduced between assignments, thus further penalising those assignments. Where the service reduction is higher the Arc Tangent Reduction in Service (ATRS) was also used, but the results were not very encouraging as some of the other approaches improved on it.When looking at the overall schedule it is evident that the number of flights requiring assignment at each time is not uniformly distributed, as shown in Figure2, such that when the number of flights to be serviced at any one time increases, the potential for conflict also increases, when perturbed, such that there is potential for further improving the results when this is taken into account, so other approaches were proposed and studied. The approaches examined, which consider the distribution of flights in time, were the ARS, BSARS, and SARS. These approaches performed better for a very low number of BSSs but in general provided solutions inferior in respect of conflicts compared to the PCBG, as well as the UESRS with exponential and Inverse estimation functions. The advantage of the TRS approach in conjunction with some of the other approaches was anticipated, as was corroborated by results from the experiments conducted for the combination of TRS with the ARS, BSARS, SARS, and UESRS.Several stochastic approaches were also considered; namely, UESRS was evaluated for different estimation functions and the exponential function withβ
=
0.03 provided the best results in respect of robustness for the measure considered. PCBG provides the best overall performance regarding robustness but did not seem to improve when combined with TRS. RNC appeared to provide solutions with a higher number of conflicts compared to the other approaches, but this could be a consequence of the low number of iterations which it was possible to execute given both the low speed, which is a characteristic of this approach, and the higher memory requirement. The UESRS approach also improved when combined with TRS, showing the Inverse function for b
=
6 with TRS providing good results for high numbers of BSSs (N
≥
L
M
A
P).When UESRS was combined with BSARS the performance improved to when these approaches where used alone or combined with TRS.It would be interesting to apply these approaches to some data sets where the original flight schedules and final real schedules are known, to see if the solutions obtained by these approaches did cope well with the changes, but unfortunately this was impossible at the time, given the unavailability of such data. Furthermore, if the data sets contain the real flight assignments to BSSs on the day of implementation, then it would be possible to quantify the actual improvement which could have been achieved by each of the approaches presented here, if they were implemented. There is also a question about the preference as to how these approaches are combined, either as a sum of each individual with the TRS, BSARS, and ARS, as studied here, or as a product of their individual contributions.Future work should consider the use of multiple distributions, in general one per flight, based on the particular characteristics of each flight, such as aircraft type, airline, destination, route, and season. This could be applied similarly to the generation of the buffer times. This approach was not used, since such information was unavailable at the time this study was conducted. It is envisaged that the use of this information in the PCBG and RNC will improve their performance and may also be used to generate the disrupted schedules used to measure the solution quality provided by all of the approaches, thus assisting in the identification of the best approach for use in the specific problem.More information about the problem studied here, the metaheuristic algorithms used, and the application of these approaches to the ABSSAP and AGAP can be found in [21].



---

*Source: 1213949-2016-09-08.xml*