Research Article | Open Access
A New Air Traffic Flow Management User-Driven Prioritisation Process for Low Volume Operator in Constraint: Simulations and Results
This document presents a new potential feature for the User Driven Prioritisation Process (UDPP) concept to give access and flexibility to Airspace Users (AUs) when they operate a low number of flights involved in a particular hotspot, a.k.a., Low Volume Users in Constraint (LVUC). Capacity constraints and congestion in the Air Traffic Management system impose delay to flights that cause large costs on airlines and passengers alike, with no significant capacity increases expected in the near-nor medium-term. Current UDPP features such as Enhanced Slot Swapping can increase flexibility for AUs to adapt their operations during capacity constrained situations. However, AUs are often impacted in their flight schedules by constraints that only affect a reduced number of flights, thus being in a situation of reduced flexibility—or no flexibility at all—to prioritise those flights. Some AUs are more vulnerable to this problem because they typically operate a low number of flights, e.g., business aviation. The new method proposed, named Flexible Credits for LVUC (FCL), is based on the use of “credits”, as a virtual currency, to increase the flexibility of LVUCs irrespective of the number of flights operated or affected by delay. FCL aims at facilitating the smooth coordination between AUs during the optimisation of their operations across multiple constraints and over the time. An initial set of simulations performed under credible conditions are presented to preliminarily analyse the feasibility and limitations of the method and to shed light on future research aspects. A first empirical evidence is given in this paper showing that increasing flexibility for LVUCs is possible without jeopardising equity.
The Air Traffic Flow and Capacity Management (ATFCM) Function of the air traffic management system aims at protecting airport and airspace sector operational capacities from excessive traffic during periods of saturated levels of traffic. For that purpose, if a significant imbalance between forecasted traffic and available capacity is predicted, a.k.a., hotspot, the ATFCM may impose delays or other measures on certain flights before they depart to maintain the operational safety level [1, 2]. The ATFCM measures that consist of applying predeparture ATFCM delays to the flights are denoted in Europe as traffic regulations  and in the US as ground delay programs .
In the current ATFCM paradigm, the delay imposed may cause flight cancellations and delay costs that impact the Airspace Users (AUs) through numerous cost elements (fuel, maintenance, crew, plus passenger “hard” and “soft” costs [4–6]). Since profitability in air transport industry is very sensitive to cost variations (profit margins were in 2017 around 7%, but they might be as low as 1-2% when fuel prices are high) [7, 8], the AUs would like more flexibility, i.e., the ability to accommodate their changing business priorities into the air traffic management system, to reduce the “impact of delay” (cost of delay) during irregular operations. However, under the current paradigm, the AUs are most of the time passive/reactive agents that have little influence on the slot allocation process.
To achieve additional flexibility for AUs to adapt their operations in a more cost-efficient manner, SESAR envisioned the development of the User Driven Prioritisation Process (UDPP) [9, 10]. UDPP is a concept that embraces several mechanisms/features, some of them at different maturity levels of development, to give flexibility to the AU to redistribute the delay across its fleet, through prioritisation of flights with high economic value over flights with relative low cost of delay. With UDPP, the AUs have more opportunities to get their most important flights, i.e., typically the ones with greater impact on their all-day schedule, to be scheduled back on track. Through a reduction of delay in their most important flights, the AUs can often be able to satisfy their most important time-constraints in the schedule (e.g., connecting flights, curfews, crew constraints, among others). Thanks to such flexibility, cancellations and delay-induced cost for freight—including reactionary delay costs—can be reduced substantially. On the other hand, for UDPP, equity is the main constraint; i.e., the actions of one AU shall not generate a direct negative impact on another AU’s flights (e.g., not increasing the delay of another’s flights).
UDPP concept is currently under development and new features are being progressively incorporated aiming to fulfil different operational requirements and implementation constraints. Some of these features have already been proposed and validated with different levels of maturity, such as Enhanced Slot Swapping, which was validated in 2015 in terms of impact on equity and acceptability by the AUs and ATFCM, and have been deployed to operations in May 2017. Other less mature features are still under development and are being validated, such as the ones explained by N. Pilon et al. . However, these UDPP features are especially available for AUs that have several flights/slots in a hotspot, a condition that will only occur at those airports where they have many flights scheduled in short periods, typically at their base airports. Unfortunately, AUs are not always in the ideal situation of having low priority flights (flights with enough margins and/or relatively low economic value) in positions nearby their most impacted flights (the high priority ones), so that they might not be able to swap their positions. A typical situation of lack of flexibility is when an AU has a small number of flights affected in a hotspot (e.g., 3 or less). In that case, the AU is called an LVUC (Low Volume User in Constraint) .
Figure 1 illustrates the lack of flexibility that AUs may experience if they are LVUCs at a given hotspot. In the figure, it is shown that AUs with a large number of flights can use the current UDPP features based on slot swapping between their own fleet, while two LVUCs in the example cannot: one AU has two flights in the sequence but if those are too far from each other they may be non-swappable; the other LVUC has only one flight, which makes it totally impossible for that AU to have access to the UDPP features based on slot swapping.
Based on the analysis of all the European airport regulations in 20 consecutive AIRAC cycles (roughly, one year and a half from Jan 2016 to July 2017, i.e. AIRACs from 1601 to 1707, taken from EUROCONTROL’s Demand Data Repository), we have found that (a) the proportion of LVUCs in daily hotspots is large, being more than 2/3 of the AUs typically affected by regulations; since a typical regulation affects 9 AUs in average, it means that 6 of these AUs would be LVUCs; see Figure 2; and, (b) all AUs are often LVUCs: in 85% of the regulations in which AUs are involved (typically hundreds every day for most commercial airlines), they typically have only a few flights (see Figure 3), meaning that in most regulations the affected AUs will find their flexibility strongly limited.
Such problem is worse for those AUs that often operate just a few flights, since they have little flexibility and may even never be able to prioritise their relatively important flights (e.g., business aviation is especially vulnerable to that problem). This is considered to be inequitable from the access KPA point of view. Therefore, there is a need to explore new features in UDPP to enable more flexibility for all the AUs, in particular for LVUCs.
This paper introduces the problem of lack of flexibility of LVUCs for using UDPP and presents a new potential UDPP feature as a solution: the Flexible Credits for LVUCs (FCL), developed under the SESAR PJ07 project (Optimised Airspace Users Operations). The study in this paper has a special focus on the design phase for reaching operational feasibility of the concept with a potential for future deployment. An early assessment through fast time simulations has been conducted and is also presented in this paper, to preliminary explore the potential operational feasibility and limitations of the proposed mechanism.
The rest of the paper unfolds as follows. Section 2 discusses the state of the art and the research gaps approached by this paper. Section 3 presents the FCL model. Section 4 discusses the preparation of the simulations and scenarios. Section 5 shows the case studies and results. Section 6 contains the conclusions and identifies future work lines. In Appendix, a brief description of the optimisation model used in the simulations is provided.
2. State of the Art
2.1. Classical and Modern Literature about Collaborative Decision-Making in Slot Allocation
ATFCM slots currently follow the First Planned First Served (FPFS) method, i.e., sorting the flights by the estimated time of arrival at or over (ETA/ETO) the constrained airport or sector, according to the information present in the filed flight plans and assigning the slots in such order. This procedure is sometimes called Ration-By-Schedule (RBS) . FPFS is widely accepted by the AUs because it preserves the original sequence of flights (considered fair), and it is well accepted today in ATFCM operations because it minimises the total delay in a regulation [12, 13]. However, although this policy was proven to be effective for the resolution of demand and capacity imbalances (hotspots), it typically does not lead to an optimal allocation of the ATFCM slots in terms of costs for the AUs [14, 15]. Due to the level of saturation in both the European and US airspaces, a massive number of flights receive ATFCM delays every day, which causes important operational disruptions and costs to the Airspace Users (AUs), the airports, and passengers [4, 5]. Similar congestion issues are also faced today—or will be in next years—in other countries due to the worldwide increase of traffic, especially in Asia [4, 16, 17].
In a real operational environment, the AUs are in the best position to take business decisions on their own flights, because they have the required internal (and typically private) information; for this reason, different collaborative decision-making (CDM) processes have been studied in the last decades to allow the AUs to participate in the ATFCM decision loop [2, 15, 18, 19], as an essential strategy to minimise the impacts of deteriorated operations on all such stakeholders [20, 21]. The network and the airports can also be benefited from AUs optimising their own operations, since the reactionary delay is minimised and the adherence to airport slots increased . In Europe, the CDM philosophy is the backbone of the UDPP concept.
Some early developments under CDM include the allocation of arrival capacity to the airlines according to the procedure RBS . See Figure 4. This procedure has improved the collaboration between airlines and ATCFM Function with respect to the previous paradigm in which the flights were sequenced giving priority to the most recent estimated time of arrival. As a result, airlines do not forfeit a slot by reporting a delay or a cancelation, which was the case prior to the CDM improvement; thus the accuracy of information about delays and cancelations has increased significantly .
An important concept introduced with the CDM philosophy is the concept of “property”, by which each AU has total control over the slots allocated to their flights by FPFS/RBS mechanism. It means that AUs have flexibility to exchange/substitute the flights allocated to the slots, without invading the allocations of competing companies (and then equity is preserved) . The trade-off between efficiency and equity is well known and has become an important topic for discussion in the last years [11, 15, 19].
The procedure for inter-airline slot exchange under CDM is called Compression . This procedure seeks to maximize, in a fair and equitable manner, the utilization of the available airport capacity during hotpot situations. If a flight allocated to a slot is cancelled, the concept of “priority” is taken into account by the system, which will try to allocate first those flights of the AU that owns the vacant slot. On this basis, more advanced inter-airline slot exchange mechanisms have been explored to increase the efficiency, but these new approaches are not today in operations, e.g., the ones proposed by Schefers et al. , by Vossen and Ball , or by Ivanov et al. .
In Europe, the slot swapping concept has also been adopted by UDPP, with some particularities to adapt it to the European context and with some operational enhancements developed under the SESAR programme. To increase the flexibility for AUs, the Enhanced Slot Swapping has refined some of the low-level operational and implementation rules, e.g., allowing multislot swapping instead of multiple binary slot swaps . In addition, some functionalities to enhance the level of usability and automation of the optimisation processes that are controlled by the AUs have been developed and partially validated, such as the Fleet Delay Apportionment concept. In this concept, the AU can send a priority list to the ATFCM system so that the system will allocate the flights of the AU according to the list and not according to the FPFS/RBS rule . Further evolutions built on the top of slot swapping concept were developed in the Selective Flight Protection feature of UDPP . With such a feature, the AU can ask ATFCM to reduce the delay of a selected flight by increasing the delay of another flight placed earlier in the sequence. This has as advantage that the AU does not need to request multiple slot swaps to reach the desired result. Another advantage is that the AU does not sacrifice flights if there is the flight protection requested that cannot be implemented. It can happen, for instance, that a hole in the sequence generated by a flight releasing a slot cannot be filled with another flight. If the system cannot reduce the delay of the protected flight, then the extra delays requested for the earlier flights are not implemented. Similar advantages can be achieved with the Slot Credit Substitution implemented in USA .
2.2. Identification of Research Gaps: Lack of Flexibility for Low Volume Operators
In all the above mechanisms, AUs can adjust their schedules by substituting and cancelling flights. The equity is preserved because there is no negative impact to flights of other AUs when slot swapping is used between flights of the same AU. If compression is applied on the other hand (e.g., after a cancellation of a flight), it typically may generate a positive impact to other flights, often to the flights of the owner of the cancelled flight and quite often to flights of other AUs too. However, none of the above mechanisms has been designed to give flexibility to LVUCs; thus all of them fail for that purpose. Indeed, up to the best of our knowledge, no previous approach has been addressing the needs of LVUCs, while it is key—compulsory in Europe indeed—from the point of view of “access and equity” to find a new feature that allows the LVUCs to participate to the UDPP mechanism.
S. Ruiz et al.  showed that due to the nature of the problem (lack of available “owned” slots by LVUCs) it may be necessary to relax the equity constraints (no impact to others) at least momentarily in some situations. However, the impact generated by these new features for LVUCs to other flights might be negligible and thus acceptable by others, or might not be negligible but still accepted by others if they can obtain some direct benefits as a kind of compensation (see Figure 5), e.g., using the same mechanism when they are LVUCs in a hotspot, or if they are compensated with some better slots for some of their high-priority flights in exchange, yet in the same hotspot or in the long-term after multiple hotspots. The compensation should not be necessarily given by exactly the same AU that causes the impact, but rather others can compensate the negative impact. What is important is that in the long-term all the AUs should have no significant deviation from their baseline delay, i.e., no significant deviation from the total FPFS delay they would have gotten at the end of the reference period if UDPP was not enabled (equity should be preserved in the long term, e.g., after 1 year), and that the net economic impact is positive for all AUs.
A survey reported by Crown Consulting Inc. (M. Ball et al.)  concluded that market mechanisms present promising characteristics to achieve the many objectives of NextGen (or SESAR) while giving flexibility to the AUs and improving the economic efficiency of the air transportation industry. Due to the lack of acceptability of the AUs to pay money directly, a system of “credits”—referred to as Priority Points in the report cited—understood as a virtual currency without monetary value could be developed to increase the level of coordination between AUs and facilitate complex slot exchanges between AUs so that efficient and equitable slot sequences can be found.
Many advanced slot allocation mechanisms have been proposed by researchers with the aim of optimising the slot allocation process, such as mechanisms based on game theory, auctions, combinatorial auctions, market mechanisms, and a system based on credits, e.g. [13, 14, 19, 21, 23–25]. However, a large set of these methods have the handicap of being too difficult to implement in real operational environments and/or too complex to be operated by human operators. Therefore, the mechanism developed for LVUCs should be based on simplicity of design and usage; e.g., long iterative negotiation processes should be avoided to avoid generating too much workload to the flight dispatchers. S. Ruiz et al.  showed some early findings based on mathematical analysis suggesting that it might be possible to find a system with such characteristics. This paper builds on the top of such findings, as explained in the following sections, and presents the first empirical evidence that increased access to flexibility might be possible for LVUCs without jeopardising equity.
3. Flexible Credits for LVUCs (FCL): A New UDPP Feature
3.1. Overview of the Mechanism
The Flexible Credits for LVUCs (FCL) is a new UDPP feature built upon the existing collaborative decision-making mechanisms that will allow the AUs to participate in the reallocation of slots. FCL is considered as an extension and complementary to other UDPP features such as the Enhanced Slot Swapping (ESS) or the Selective Flight Protection (SFP).
The new proposed mechanism is based on the use of delay credits (DCs), a virtual currency without monetary value that can be used by LVUCs to enhance their flexibility for redistributing delays among their flights. To make it operationally acceptable for AUs, the FCL feature aims at facilitating the smooth coordination between AUs for the optimisation of their operations during hotspots.
The mechanism follows a ration-by-effort principle, i.e., aiming at keeping a fair balance between the delay saved and given to others. Thus, DCs can be earned and accumulated by AUs when they are LVUCs only if they accept extra delay in their low-priority flights. This will reduce the delay to flights of other AUs. DCs will be used in the right proportion to reduce the delay of their high-priority flights. One important characteristic of the proposed mechanism is that credits can be transferred from one hotspot to another, so LVUCs can increase their chances to prioritise their flights even when they only operate 1 flight in a hotspot. LVUCs can decide at any moment if credits are fully or partially transferred from one hotspot to another or alternatively can be used in the same hotspot (in case that the LVUC has more than one flight in a hotspot).
Business aviation operators usually have few flights involved in a given hotspot, so they will benefit greatly from the introduction of this new mechanism. For instance, an AU operating just one flight may accept more delay in a hotspot when its flight is far from its operational margins (in such situation its cost of delay is relatively low). Giving up part of its operational margins will help other AUs to reduce their delays, and in exchange the AU will earn some credits that he will be able to use in future hotspots to reduce the overall cost when his flight is impacted severely by delay (beyond the operational margins).
Consider the illustrative case in Figure 6 in which an LVUC has only one flight (FL001) in Hotspot 1 and one flight (FL002) in Hotspot 2. In both cases, the AU cannot improve his situation (lower the cost impact of delay) with the current Enhanced Slot Swapping, FDR nor SFP features. However, with FCL the AU could have access to UDPP and could substantially improve his situation. Note that in this example the flight FL001 has a certain amount of delay (D1) that could be increased, due to the operational margins available, with relatively low impact in terms of cost. Accepting the extra delay (D1’ > D1), the AU is giving up his position and reducing the delay (positive impact) of the AUs between the original and the new sequence position. The AU can then be rewarded for this, with an amount of delay credits proportional to the extra delay accepted (in this document 1 DC = 1 minute of delay). In Hotspot 2, the LVUC could use the credits available from Hotspot 1, because flight FL002 has an amount of delay D2 that has an important impact in terms of costs for the LVUC. After using part or all of the credits available, the AU can reduce the delay for that flight from D2 to D2’, which is an amount of delay within the operational margins available for that flight, and therefore with marginal cost of delay. The AUs between the baseline and the new position are impacted negatively, but with just one position in the sequence (typically it means 2 or 3 minutes or less of extra delay).
The consideration of an AU as a LVUC may change in different hotspots according to the circumstances; therefore, even large airlines can often be considered as LVUCs in many hotspots (typically in hotspots at airports in which they operate a few flights). The access to FCL and LVUC rules to any AU—in those hotspots in which the AU can be considered as LVUC—can help to accept some degree of inequity in favour of LVUCs at some moment in time (but with equity compensated over time).
Another example is illustrated in Figure 7. The sequence in the example is dominated in presence by two AUs that are not LVUCs, represented by sticks coloured in green and blue, whereas there is a LVUC with only two flights, identified in the sequence in orange. The example illustrates how the LVUC can reduce the total impact of delay (costs) in the same hotspot by increasing +24 minutes the delay of flight 1, to obtain 24 delay credits, and then reducing the delay of flight 2 by giving the 24 credits gotten previously. In this example, a 20% of cost reduction could be achieved, due to the nonlinear relationship between the fleet cost structure and its delay. The impact on other flights can be considered negligible (3 minutes of extra delay per flight, if impacted), according to the AU experts that have reviewed such figures (further research will be conducted to fully validate such assumption).
3.2. FCL Rules to Preserve Equity under High Flexibility Conditions
In FCL, high flexibility is given to LVUCs to minimise their own global delay costs; i.e., the LVUC has full freedom to transfer its total baseline delay (i.e., initial ATFCM delay) among its flights and to exchange freely flight sequence positions with other AUs while only being subject to two particular equity constraints: AU’s total baseline delay cannot be reduced, and the Maximum Negative Impact of Time (MNIT) for individual flights in a hotspot must be respected.
The main rules of FCL are defined as follows:(1)Any AU with a given number of flights (3 or less is considered as an initial proposal) in a hotspot is considered as LVUC in the context of such a hotspot (note: perhaps other criteria might be used in the future to determine whether an AU is an LVUC in a given hotspot, e.g., the AUs’ share of flights in the hotspot, but such a concept requires further research).(2)Any LVUC in a hotspot can save credits obtained by increasing the delay in some flights (thus giving better sequence positions to other AUs) and use them in that particular hotspot or keep them to use in future hotspots.(3)Any LVUC in a hotspot can use Leftover Operational Credits (LOCs) obtained in the past to protect flights in the current hotspot; LOCs could be weighted up to take into account the differences of the hotspots in terms of duration and severity, so that the equity among AUs can be preserved (out of the scope of this research, though).(4)LVUCs can request any target place in the sequence (e.g., corresponding to the optimal delay allocation for that flight), irrespective of where and when the efforts were done to obtain the credits (note that this rule applies only for LVUCs, which, by definition, are expected to have less flexibility than other AUs; it may occur that in a particular hotspot the only flight available for a suspension could be in positions after the one that is tactically important for the LVUC).(5)The AUs’ total delay at the end of the reference period (e.g., 1 year) must be the same (or near the same) as the baseline delay (delay without UDPP, e.g., FPFS).(6)The Maximum Negative Impact of Time (MNIT) represents the maximum minutes of additional delay that LVUCs can cause to a flight of another AU. LVUCs might be unable to protect and reduce the delay for a particular flight if the MNIT has been reached for any of the flights it affects.(7)All the requests will be sorted by the requested time and will be integrated in the sequence in FIFO order and by minimising the impact to others. No empty positions should be found in the sequence (compression) and requests that could generate impact greater than MNIT will be allocated to the nearest feasible solution.
3.3. Assumptions and Principles of the Mechanism
The proposed method is based on the following assumptions:(i)All transactions of credits and slots exchanges will be initiated through the Airport or ATFCM actor triggering the UDPP mechanism, which will act as a broker and supervisor.(ii)AUs must be consistent and accept any consequence derived from their own decisions.(iii)Regarding the NM and the CFL system, one minute of delay will be considered as having the same “value” for all the flights (i.e., 1 minute delay = 1 DC). (note: to have more control on equity aspects, future research may consider different value of delay depending on the position in the sequence and the level of delay already allocated to such position).(iv)No negative delays are allowed during the slot reallocation (this could be relaxed later, but for the moment this assumption is taken for the sake of simplification).(v)No uncertainty associated with the new slot allocations planned during UDPP is considered (i.e., Confidence Index = 1).(vi)Having credits is useful and positive for the AUs since credits provide flexibility to adapt the operations to changing and unforeseen conditions. AUs must contribute to the network with something in exchange, e.g., assuming extra delays for some flights under a regulation, to give their sequence positions to other AUs.(vii)No loss of value or expiry of credits is currently considered (although this may be reviewed later).(viii)AUs are able to and have the necessary information (e.g., historical records) to make decisions involving different hotspots over the time, thus making a stochastic/probabilistic management of delay and management of impact of delay.(ix)The FCL decisions should be made and supervised carefully by expert human operators that will ensure the stability of the system and the application of good practice (i.e., automation shall only be provided to support human decision-making and control).(x)If an AU sends no sequence position request for a given participating flight, it is assumed that it is willing to take advantage of any delay reduction opportunity for such flight (this assumption is consistent with the assumptions behind the compression algorithm detailed by T. bossen and M. Ball ).(xi)FCL can coexist with the FDR and SFP features that will be used by non-LVUC operators.(xii)Special care should be taken when extending the validity of the credits through different hotspots, as the severity, duration, or the flights involved are different in each situation. Hence, the credits should not have the same value in all hotspots. This is out of the scope of this research, but in future research an equivalence factor could be developed to update the number of credits to transfer credits from one hotspot to another, being equitable for the AUs involved, e.g., associating the number of credits given or requested according to the amount of delay associated with the positions released or taken.
4. Preparation of the Simulations
4.1. User Delay Optimisation Model (UDOM) Used in the Simulations
In order to generate trustworthy evidence through simulations, it is necessary to approximate the expected decisions potentially made by LVUCs using FCL with regard to the management of delay. For that purpose, the User Delay Optimisation Model (UDOM) presented by S. Ruiz et al.  has been used. A summary of the model together with an extension to include some implementation aspects can be found in Appendix.
4.2. Description of the Traffic Scenario
The scenario used for the case study has been previously used in some UDPP validation exercises and has been slightly adapted for the purpose of this new research. It is based on historical traffic demand of 96 flights at a coordinated airport from 12:15 to 17:00 approximately.
Figure 8 shows the schedule list in detail, showing the Expected Time of Arrival (ETA) planned for each flight, the new Controlled Time of Arrival (CTA) assigned by FPFS policy (CTA FPFS), and the calculated delay, also the new user-preferred CTA after the application of FCL prioritizations (UCTA FCL) and the calculated delay. The last column shows, for each flight, the difference in the delay allocated by FPFS and ESFP mechanisms. Note that ESFP takes as a baseline the same policy applied at a given airport, in this case FPFS; thus at the beginning there is no difference in the baseline delay.
Figure 9 shows the evolution of the delay in the hotspot, which is directly correlated with the number of positions, , that a flight i is far from its original position in the sequence (delay is proportional to the number of positions and the size in time of each position in the sequence).
Table 1 shows the distribution of flights and delay per Airspace User. Note that delay is distributed among AUs, after aggregating the delay of all their flights, as a direct proportion of the share of flights. A certain degree of concentration of some flights has been identified after analysing the flight positions, especially for the user HUB and LC2. In this particular case, HUB has many of its flights in positions with low delay and LC2 has the flights concentrated in positions subject to relatively long delay.
Note that some AUs have three or less flights in the sequence and therefore they can be considered as LVUCs.
4.3. Description of the Cost Model Used and Its Parameterisation
A simplified cost model has been taken into account in the simulations to quantify the impact of delay on flights. The costs of delay for flights have been modelled with a quadratic function as follows:Using a quadratic model instead of more realistic nonsmooth and nonconvex utility curves is a strong simplification (see Figure 16) that cannot be made in all the contexts to generate reliable evidence. However, for the purpose of this research, such simplification brings more advantages than disadvantages. For instance, it eases the analytical study—with UDOM—of the dominant strategies of the AUs followed during the optimal redistribution of their delay, and thus to anticipate the expected dynamics in the system, which is the main purpose of this research. The emergent dynamics of the system will be analysed under the hypothetical what-if scenario in which some AUs are allowed to optimise the allocation of their delay by redistributing, under high-flexibility conditions, their baseline delay among their constrained flights. It should be noted that the underlying system dynamics should be similar either with the simplified model or with more realistic utility curves, since these dynamics depend mostly on the trends driven by the AUs’ dominant strategies. Such simplified approach also allows calculating a first-order approximation to the cost reduction potential in relative terms (absolute realistic cost figures are hardly difficult to find, due to their privacy nature and their strategic value for airlines). For a more refined study of the potential benefits for the AUs, it may be necessary to adopt realistic nonsmooth curves and the approached with nonsmooth optimisation techniques (see, e.g., G. Stojkovic et al. , or R. Hoffman ).
Each flight of the scenario has been randomly parameterised with different sensitivities to delay, i.e., with different (a.k.a., elasticity). The assignation of epsilons/elasticities for each flight has been done randomly following a uniform probabilistic distribution bounded within the range . All the flights had the same probability to get any value of elasticity within such range. In practice, it means that all the flights were assigned costs as a function of their delay within the maximum and minimum range illustrated in Table 2.
Such approximation of AUs costs is not fully representative of the high level of complexity of the actual cost structures (e.g., different types and sizes of knock-on delay impacts). However, the approach is considered valid at this maturity level to observe the potential benefits and/or impacts of FCL when flights and AUs with heterogeneous cost curves coexist in the same hotspot. The purpose is to early assess the potential feasibility of the FCL concept, to enable the observation of potential emergent dynamics, and to start quantifying the potential benefits of flexibility for the LVUCs and for the rest of the AUs.
Table 1 shows the baseline delay and baseline costs per AU. Note that since the elasticities were randomly assigned to flights it was already expected that AUs had a proportion of costs similar to the share of delay.
5. Case Studies and Results
The following three case studies were simulated and analysed with the purpose of generating early empirical evidence about the potential feasibility and limitations of the FCL method. Case study 1 answers the question of whether the FCL feature can be used by all the AUs and for all the flights in a hotspot; case study 2 sheds light on whether the FCL mechanism could generate acceptable outcomes with a reduced number of flights using the mechanism in a single hotspot (3 flights was considered appropriate by a panel of AUs experts as a starting point for the analysis); case study 3 illustrates an example similar to case study 2 but including delay credits that are assumed to be transferred from past hotspots.
5.1. Case Study 1: All the AUs Using the FCL System
A first question that seems relevant to answer is whether the high flexibility conditions given by FCL could be used by all the AUs and not only by those considered LVUCs in a given hotspot. This question was preliminary discussed by S. Ruiz et al. in . After that, several AU experts showed interest in exploring the possibility of using FCL by all the AUs, i.e., not only for LVUCs. This case study explores this idea.
Figure 10 shows the requested optimal delay changes (taus) calculated with UDOM, as well as the impact of the delay change on each flight. For convenience, it has been assumed that only the flights in and after the position of flight LC1 2612, a flight that has 9 minutes of delay, could be included in the FCL optimisation; i.e., the AUs could not use FCL to change the position of those flights earlier than flight LC1 2612.
Table 3 shows the aggregated impact on the costs per AU. Note that in relative terms all the AUs can optimise their costs in similar proportions, and it is worthy to point out that the AU with more potential for cost reduction has only 3 flights (an LVUC). The potential for cost reduction in this scenario for the AUs has been found within the range of approximately 7% up to 36% (average cost reduction of around 17%). The potential for cost reduction depends on the cost structure of the flights in the hotspot, and on the amount of delay allocated by the baseline sequence to each of the flights. Also it must be observed that some AUs—OA2, OA3, OA4, and OA7—cannot minimise their costs. In the case of OA2, OA3, and OA4, they have indeed just one of their flights in the hotspot (OA2 and OA3 have another flight but were not candidates for using FCL, since they are considered to be out of the hotspot); thus they cannot exchange delay between their flights. The use of credits with validity over the time could be used for giving flexibility to those AUs.
Figure 11 shows the slots demanded by the AUs with the aim of optimising their own costs. It can be observed in the figure that the demand is more or less evenly distributed. The presence of holes (slots not demanded), which are also spread more or less evenly along the sequence, is a consequence of the random distribution of cost structures (with different epsilons randomly assigned), together with the effect of the equity rule that forces the AUs to release slots in exchange of better slots for their relatively more expensive flights. Nevertheless, due to the need to allocate one flight to each available sequence position, the flights have been sorted by requested slot time and then they have been compressed to avoid empty slots and to optimise capacity.
Figure 12 is an evidence of complex dynamics and inter-AUs demand incompatibilities emerging after trying to join all the individual requests in a single sequence. The figure shows at each sequence position the difference in minutes between the allocated position and the one requested. In general, it can be expected that less delay than the one requested (negative values in the chart) will have positive impact on AUs, which is advantageous for those flights. In the middle of the sequence, due to the limited capacity with respect to new demand for slots, some flights can be found that requested less delay but that will receive more delay than requested (they will be placed at later positions). Therefore, due to the incompatibilities of the many AUs requests, some flights are getting less delay (and less cost) and some more delay (and more cost). Table 4 shows the total effect in the total cost supported by the AUs, while Table 5 shows the new distribution of delay per AU. It can be noted that some AUs are significantly worse-off in terms of delays and costs after trying to optimise their costs through the FCL mechanism.
5.2. Case Study 2: LVUCs with up to 3 Flights in the Hotspot (Single Hotspot Optimisation)
In this case study, the purpose is to assess the benefits of enabling FLC mechanism only for LVUCs that have up to three flights. The impact to others is also assessed. As shown in Figure 13, only OA2, OA3, O4, and OA5 are LVUCs in this hotspot scenario. Note that OA3 and OA4 have 1 flight each in the hotspot (OA3 has another flight in the scenario, but is placed before the hotspot), OA2 has 2 flights, and OA5 has 3 flights. In the figure, they can be found highlighted with different colours. The optimal delay requests found with UDOM for those AUs can also be found in the figure, together with the cost variation for those flights.
In this case study, optimisation is only allowed within the hotspot; i.e., AUs cannot optimise their delay throughout multiple hotspots. Thus, the LVUCs with one single flight, OA3 and OA4, will not have flexibility to minimise the impact of their baseline delay.
Figure 14 shows the changes in the sequence after the implementation of the LVUCs requests. It can be observed that the maximum impact on individual flights of other AUs has been 3 minutes (corresponding in this scenario to one sequence position).
Table 6 shows the aggregated variation of delay for the different AUs, which in the light of the figures and after consultation with a panel of operational experts (which included AUs representatives) it can be considered as negligible impact. Note that LC2 is the AU receiving more delay, i.e., 9 minutes, with no more than 3 minutes per flight. Note as well that OA2 and OA5 have the same delay as the baseline (no impact).
In terms of costs, Table 7 shows that OA2 and OA5 could potentially optimise their costs to figures that are quite close to what UDOM found, i.e., 22.34% of cost reduction for OA2 (who expected a 25.61% of cost reduction) and 38.38% of cost reduction for OA6 (who aimed at reducing the cost in 36.21%).
Regarding the impact on the rest of the AUs, some relatively small variations can be observed in Table 7. It must be pointed out that variation of 1-2% of cost could be considered as very huge impacts in reality; however to interpret the results of this simulation exercises it must be reminded that flights typically operate with some tolerances that makes the flight costs to be typically not so sensitive to an extra minute of delay (the quadratic cost structure is most likely not representing well the impact to others). In addition, it is assumed that the AUs that are not LVUCs may participate in other UDPP features (e.g., through Enhanced Slot Swapping, FDR/SFP, or others), thus optimising their costs and most likely compensating the relatively small impacts received from LVUCs that used FCL.
5.3. Case Study 3: LVUCs with up to 3 Flights in the Hotspot (Multiple Hotspot Optimisation)
To illustrate the impact of LVUCs optimising over the time (multiple hotspots) let us assume that OA3 and OA4 can reduce the delay of their flights because they made efforts in other hotspots (i.e., they increased the delay for some of their flights, thus reducing the delay to flights of other AUs). In the paper of S. Ruiz et al. , the reader can find examples of how the UDOM can be adapted to optimise the costs of operations by managing delay over the time across multiple hotspots.
For the purpose of this case study, let us assume that in the same situation as shown in previous case study 2, the LVUCs OA3 and OA4 want to reduce the delay of their flights in 18 minutes each (six sequence positions in this hotspot scenario). It is assumed that proportional efforts were done by such AUs in the past in order to have the right for reducing their delay in this occasion.
Figure 15 shows the new sequence situation after the implementation of all the LVUCs requests. It can be noted that in this case a maximum delay impact of 6 minutes has been found in some flights (the flight of OA3 and the flight OA5 108 overcome the same three flights, thus generating an impact of two sequence positions, i.e., 6 minutes, for some flights).
Tables 8 and 9 show similar figures as in case study 2 for delay and cost impacts. After consultation in an expert panel workshop (which included AUs representatives), it was concluded the impact on delay caused by LVUCs using FCL could again be considered negligible, whereas the impact to other in terms of costs is not conclusive due to the simplifications done in the cost curves of flights. OA2 and OA5 could optimise their costs in the same order of magnitude as calculated by UDOM. The figures shown for OA3 and OA4 are not taking into account the cost of the efforts done in the past hotspots to have flexibility in the studied hotspot. It is expected however that they should be able to reduce their total—over the time—costs in the same order of magnitude as OA2 and OA5.
6. Discussion of the Simulation Results
With the current FCL rules, and due to the complex interactions and emergent dynamics between the different AUs prioritisations, it cannot be guaranteed that the resulting FCL sequence will have a positive or neutral impact on the AUs costs. In the light of the results (see case study 1), it seems quite unlikely that all the AUs can use FCL in its current form. Further research needs to be done to investigate new strategies and rules that will allow to develop solutions with positive cost reduction for all the AUs. The FCL feature seems to be a good mechanism to give access for LVUCs with 3 flights or less. In addition, it deserves special mention that, for those LVUCs with only 1 flight, the FCL mechanism was able to give access to UDPP (i.e., flexibility) through multiple hotspots, thus guaranteeing access and equity of UDPP.
For a reduced number of LVUCs using the FCL, the expected impact to other AUs can in general be considered negligible (to be further confirmed with operational experts after the analysis of many more cases), while the benefits for the LVUCs using FCL may be notably large. In addition, it can be assumed that the—negligible—impact that LVUCs in a hotspot could generate to other non-LVUCs flights could be easily compensated by the important cost reductions that such non-LVUCs could achieve through access to other UDPP features (e.g., Enhanced Slot Swapping). In addition, according to the results of the statistical analysis (Figure 3), most of the AUs are LVUCs quite often: in 85% of the hotspots that affect their flights; therefore, the AUs that are not LVUCs in a specific hotspot scenario could potentially be benefited from the FCLs feature in other hotspots in which they might be LVUCs.
The hotspot scenario used in the simulated case studies is quite typical in terms of duration, severity, and number of AUs and LVUCs affected; therefore, it is expected that a certain degree of extrapolation of the results to other hotspot situations is possible.
Finally, the case studies considered that LVUCs are the AUs with a maximum number of flights caught in the hotspot equal to three. Such definition of LVUC was preliminary agreed through discussion with a panel of AUs and operational experts. However, in the light of this early results, a sensitivity analysis should be conducted with multiple realistic scenarios in different conditions and with different AUs involved to determine which should be the number of flights that could use FCL in each hotspot while preserving equity over time.
7. Conclusions and Future Work
The FCL mechanism has shown good properties to give effective access for LVUCs to UDPP. In particular, it has been found that (a) a LVUC can find for each of his flights the optimal amount of delay minutes that should be increased or decreased in order to minimise his total costs; (b) a LVUC can obtain delay credits by accepting extra delay (amount of extra delay controlled/decided by the LVUC) and use the credits in the same hotspot (if the LVUC has more than one flight) or save the credits and use them in future hotspots (i.e., LVUCs are able to optimise their operations over the time); and (c) the FCL rules are simple and the level of coordination required is very efficient, so the LVUCs and the rest of the AUs can concentrate on optimising their own operations.
FCL cannot be used—in its current form of baseline rules and principles—by AUs in a hotspot that are not LVUCs due to the likeliness of emerging inequities. However, early evidence showed that equity could be preserved over the time if the number of flights being prioritised with FCL and the maximum negative impact allowed to other flights is limited at each hotspot. The precise definition of these parameters to control the trade-off between flexibility and equity requires further research. Three flights could be recommended as a starting point to determine whether an AU is an LVUC in a hotspot, while six minutes of additional negative impact at some flights could be considered in principle negligible by the AUs.
For the flights that are not allowed to use FCL in a hotspot, the AUs could use the other UDPP features based on slot swapping, since such mechanisms can guarantee full equity over the time (“equity” understood as no direct negative delay impact to others) and may likely provide enough flexibility to the AUs when they are not LVUCs. Nevertheless, all the AUs are often LVUCs—even the largest airlines—in many regulations that occur every day. Therefore the possibility to include all the AUs—irrespective of the number of flights operated per day—using FCL when they are LVUCs in a hotspot is recommended, since it may contribute to reduce the cost for all the AUs and possibly increase the predictability and robustness of fleet schedules (e.g., for a hub carrier a single flight delayed in a remote/non-hub airport could have a knock-on impact on the hub connections and cause disruption on the full-day’s operations). Further research is needed to determine whether special rules may apply for AUs that are always LVUCs, e.g., Business Aviation.
Further research and validation is required to determine the exact number of flights that could use FCL in any hotspot without jeopardising equity over the time and to bring the FCL concept towards higher levels of maturity.
A. User Delay Optimisation Model (UDOM)
A.1. Definitions and Assumptions
One of the main assumptions in UDOM is that the AUs taking part in the system can be modelled as utility maximising; i.e., the major objective of each AU is to maximise its utility function. Therefore, AUs are assumed to have a utility function, which is depending on several variables, such as the delay cost structure of each flight.
Utility is an important concept in economics and game theory, because it represents satisfaction experienced by the consumer of a good (see H. Varian ). In the context of this document, the concept of “utility” will be understood as the value perceived by a particular AU if a given slot is allocated to a particular flight operated. Without loss of generality, in this document it is assumed that utility is directly related to economic profits obtained by the AUs for operating their flights; however the concept of utility may also include any type of operational constraints known by the flight dispatcher, and any indirect economic or non-economic type of benefits or costs.
Table 10 describes all the mathematical symbols that will be used in the model formulation.
A.2. Mathematical Representation of the Utility or Delay Cost Functions for Flights
A utility function for a single flight, U(t), can be represented analytically as a continuous quadratic function, as it can be observed in Figure 16. If negative delay is not considered (simplification), the utility as a function of the delay, d, assigned to a flight can be expressed aswhere is the elasticity of the utility function, the marginal utility, and the maximum utility in case of no delay allocated to that flight. It is assumed for the sake of simplicity that every single flight has a maximum utility when the delay, d, is zero, at a certain slot, and the utility is then progressively decreasing as far as the delay is increasing. Since the elasticity is negative, any delay incurred by a flight will cause a reduction on the utility perceived by the AU for this flight. Utility will be negative if the cost of delay has become greater than the maximum economic value expected for that flight if operated on time. Different revenue and cost structures can be modelled by changing the parameters; thus they can be adapted to different types of AUs’ type of activity (low cost, HUB, business aviation, or others).
In reality, utility functions are unknown and may be nonlinear and nonconvex (the cost functions as well). However, in this document it is assumed for the sake of simplicity that every single flight has a maximum utility when the delay, d, is zero, at a certain slot, and the utility is then progressively decreasing as far as the delay is increasing. Utility will be negative if the cost of delay has become greater than the maximum economic value expected for that flight if operated on time.
A.3. UDOM for N Flights and Multiple Hotspots
The used delay optimisation model can be expressed for N flights of the same AU, each of those flights characterised by a different elasticity and affected by delays with different probabilities and with a particular and different expected delay for each of them:After some mathematical development (e.g., using multipliers of Lagrange), the optimal delay shift for each flight i can be expressed byFurther explanations about how this equation has been developed can be found in Section A.4 where the UDOM model will be extended with inequality constraints to consider the duration of the hotspot.
Figure 17 illustrates an example in which an LVUC has 3 flights, F1, F2, and F3, each one with a maximum utility of 500 monetary units if they were not delayed. The three flights have been originally allocated to some ATFCM slots that generate to them 5, 12, and 20 minutes of delay, respectively. Each of the flights has different cost structure and thus different sensitivity/elasticity to delay, i.e., , , and . The baseline delay utility can be calculated asAfter computing the optimal delay allocation with UDOM, the LVUC would request a slot with 21 minutes of extra delay for flight F1 and a delay reduction of 7 and 14 minutes for flights F2 and F3, respectively. Note that the total delay for the LVUC is the same. The optimised utility is illustrating the change in utility for the AU, and comparing the maximum utility affordable (no delay), Figure 18 expected utility if the baseline FPFS sequence is applied and the optimised sequence if the LVUC is allowed to use FCL.
A.4. Implementation Aspects and Generalisation of UDOM with Inequality Constraints (KKT Conditions)
The method of Lagrange multipliers allows solving convex optimisation problems subject to equality constraints, as presented in (A.3). However, for implementing and running simulations with realistic hotspot scenarios some extra inequality constraints are needed, in particular the hotspot time-window boundaries. If these constraints are not included in the model, in some cases, the AUs with flights close to the boundaries of the hotspot could request, after optimising their delay with UDOM, slots that were beyond these boundaries, thus leading to nonfeasible results in practice. In order to incorporate the hotspot boundary constraints, a new set of inequality constraints is needed.
The mathematical model to optimise becomes where is the utility function of flight i , , and are, respectively, the delay allocated to flight i, the delay shift (variable of control) that the AU can decide to optimise his costs, and the probability of that flight for being delayed. and are the new parameters with respect to the basic UDOM shown in the paper, which are , the distance in time between the slot allocated to flight i, and the initial time of the hotspot, and , the distance in time between the slot allocated to flight i, and the final time of the hotspot.
The Karush-Kuhn-Tucker (KKT) conditions allow generalising the method of Lagrange multipliers to include inequality constraints. Let a generic optimisation problem be expressed as follows (vector form): Any optimal solution for the above problem must fulfil the following four KKT conditions:
(2) Primal Feasibility
(3) Dual Feasibility
(4) Complementary SlacknessIn practice, the resolution of problems with KKT conditions are often solved with recursive algorithms. The fourth condition forces the solution to be at the boundary of an inequality constraint when an optimal solution would be beyond such boundary in case that such inequality constraint would not be present in the problem.
Taking the above into account, the UDOM model can be generalised with the KKT conditions by applying the following recursive algorithm:(1)Let be the total excess delay of the AU out of the hotspot boundaries calculated as follows (note the inequality constraints are incorporated in the calculation of ): The recursive algorithm must be initialised with (assume no excess before the first trial/iteration).(2)Find the optimal solution for the general UDOM-KKT problem, which can be expressed as (3)Calculate with (A.11).(4)For any flight i with (i.e., for each flight it checks if the optimal solution found in previous steps requires slots that are out of the hotspot boundaries), do the following:(a)Set the solution/position for flights with excess at the boundary (i.e., they will request the first or the last slots in the hotspot).(b)Remove these flights from the optimisation in the next iterations (they already have been assigned a solution).(5)Repeat steps , , and until ( must be zero for any feasible solution).
Note that the Lagrange condition for optimality, , applied to the general UDOM-KKT problem in (A.12) generates the following system of equations (dual problem):Once the system of equations is solved, the optimal delay per flight can be found:If in the first iteration (i.e., all the optimal/requested slots are within the hotspot boundaries), the solution is exactly the same as in (A.3), i.e., the normal UDOM (no UDOM-KKT). And if the delay excess generated by flights that have been assigned to the boundaries of the hotspot will be shared among the flights of the AU according to the relative importance/cost structure of each flight (indeed in the same proportion as the total delay is reallocated by the AU among his flights). Therefore, once the recursive algorithm converges to , the KKT conditions can be satisfied, meaning that an optimal solution for the AU has been found.
The flight list data and cost parameterisations used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
- EUROCONTROL, “ATFCM Users Manual,” Ed. 22.1, 2018.
- M. O. Ball, C. Chen, R. Hoffman, and T. Vossen, “Collaborative Decision Making in Air Traffic Management: Current and Future Research Directions,” in New Concepts and Methods in Air Traffic Management, Transportation Analysis, pp. 17–30, Springer Berlin Heidelberg, Berlin, Germany, 2001.
- FAA, “Interface control document for substitutions during ground delay programs, ground stops, and airspace flow programs,” Version 3.1, 2007.
- EUROCONTROL, “Performance Review Report, an assessment of air traffic management in Europe during the calendar year 2017,” 2017.
- Univerity of Westminster, “European airline delay cost reference values - updated and extended values,” Ed. 4.1. 2015.
- D. Xiaoxu, H. Minghua, W. Tian, and H. Liu, “Modeling congestion propagation in multistage schedule within an airport network,” Journal of Advanced Transportation, vol. 2018, Article ID 6814348, 11 pages, 2018.
- IATA ECONOMICS, “Profitability and the air transport value chain,” IATA ECONOMICS N°10, 2013.
- IATA ECONOMICS, “Annual review 2018,” 2018.
- “SESAR P.07.06.04. UDPP Step 2 V1 OSED. Partners: EUROCONTROL, AA, AF-Group, DLH-Group, EBAA, ELFAA, IAOPA, SAS, IATA,” 2013.
- N. Pilon, S. Ruiz, A. Bujor, A. Cook, and L. Castelli, Improved Flexibility and Equity for Airspace Users during Demand-Capacity Imbalance - An Introduction to The User-Driven Prioritisation Process, Sixth SESAR Innovation Days, Delft, Netherlands, 2016.
- S. Ruiz, L. Guichard, and N. Pilon, Optimal Delay Allocation under High Flexibility Conditions during Demand-Capacity Imbalance, SESAR Innovation Days, 2017.
- T. Vossen and M. Ball, “Optimisation and mediated bartering models for ground delay programs,” Naval Research Logistics, vol. 53, no. 1, pp. 75–90, 2006.
- L. Castelli, R. Pesenti, and A. Ranieri, “The design of a market mechanism to allocate air traffic flow management slots,” Transportation Research Part C: Emerging Technologies, vol. 19, no. 5, pp. 931–943, 2011.
- N. Schefers, J. J. Ramos González, P. Folch, and J. L. Munoz-Gamarra, “A constraint programming model with time uncertainty for cooperative flight departures,” Transportation Research Part C: Emerging Technologies, vol. 96, pp. 170–191, 2018.
- R. Jovanovic, O. Babic, V. Toic, and V. Tošic, “Pricing to reconcile predictability, efficiency and equity in ATM,” in Proceedings of the 11th USA/Europe ATM R&D Seminar, vol. 2015, Lisbon.
- A. Harada, T. Ezaki, T. Wakayama, and K. Oka, “Air traffic efficiency analysis of airliner scheduled flights using collaborative actions for renovation of air traffic systems open data,” Journal of Advanced Transportation, vol. 2018, Article ID 2734763, 14 pages, 2018.
- T. Kistan, A. Gardi, R. Sabatini, S. Ramasamy, and E. Batuwangala, “An evolutionary outlook of air traffic flow management techniques,” Progress in Aerospace Sciences, vol. 88, pp. 15–42, 2017.
- A. R. Odoni, L. Bianco, and G. Szegö, “The flow management problem in air traffic control,” in Flow Control of Congested Networks, Springer, Berlin, Germany, 1987.
- M. Ball, G. Donohue, and K. Hoffman, “Auctions for the safe, efficient and equitable allocation of airspace system resources,” in Combinatorial Auctions, P. Cramton, Y. Shoham, and R. Steinberg, Eds., pp. 507–538, MIT Press, Cambridge, 2005.
- R. H. Coase, The Problem of Social Cost, 1960.
- N. Ivanov, F. Netjasov, R. Jovanović, S. Starita, and A. Strauss, “Air Traffic Flow Management slot allocation to minimize propagated delay and improve airport slot adherence,” Transportation Research Part A: Policy and Practice, vol. 95, pp. 183–197, 2017.
- “Crown Consluting, Inc. (for NextGen). Flight Priority Deep Dive. Final Report,” 2011.
- V. Krishna, Auction Theory, Academic Press, 2nd edition, 2009.
- R. Gibbons, Game Theory for Applied Economists, Princeton University Press, Princeton, NJ, USA, 1992.
- S. J. Rassenti, V. L. Smith, and R. L. Bulfin, “A Combinatorial Auction Mechanism for Airport Time Slot Allocation,” RAND Journal of Economics, vol. 13, no. 2, p. 402, 1982.
- G. Stojković, F. Soumis, J. Desrosiers, and M. M. Solomon, “An optimization model for a real-time flight scheduling problem,” Transportation Research Part A: Policy and Practice, vol. 36, no. 9, pp. 779–788, 2002.
- R. L. Hoffman, Integer Programming Models for Ground- Holding in Air Traffic Flow Management, Department ISR, 1997.
- H. Varian, Intermediate Microeconomics: A Modern Approach, W. W. Norton & Company, 2009.
Copyright © 2019 Sergio Ruiz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.