Research Article  Open Access
Yu Jiang, Zhaolong Xu, Xinxing Xu, Zhihua Liao, Yuxiao Luo, "A Schedule Optimization Model on Multirunway Based on Ant Colony Algorithm", Mathematical Problems in Engineering, vol. 2014, Article ID 368208, 11 pages, 2014. https://doi.org/10.1155/2014/368208
A Schedule Optimization Model on Multirunway Based on Ant Colony Algorithm
Abstract
In order to make full use of the slot of runway, reduce flight delay, and ensure fairness among airlines, a schedule optimization model for arrivaldeparture flights is established in the paper. The total delay cost and fairness among airlines are two objective functions. The ant colony algorithm is adopted to solve this problem and the result is more efficient and reasonable when compared with FCFS (first come first served) strategy. Optimization results show that the flight delay and fair deviation are decreased by 42.22% and 38.64%, respectively. Therefore, the optimization model makes great significance in reducing flight delay and improving the fairness among all airlines.
1. Introduction
With the rapid development of the Chinese civil aviation industry, the number of flights increases sharply and hub airports change their single runway to multirunway. The air traffic control managers put first come first served (FCFS) to use to schedule the arrivaldeparture flights, which can lead to the waste of air resources and make the terminal area congestion more serious. Therefore, the conflict between demand and supply is more and more sharp Hu and Paolo et al [1]. Meanwhile, with the implementation of collaborative decision making (CDM) mechanism in airport resources management, flight scheduling problem should account for not only flight delay but also equity among airlines. A more efficient and reasonable model or algorithm is needed to approach the problem. Therefore, the air traffic congestion will be alleviated and the total operation cost of airlines will be reduced by approaching the scheduling and optimization of arrivaldeparture flights.
In recent years, the capacity of most airports cannot satisfy the rapid increase demand because of the increase in number of flights and severe weather. At present, ground delay procedure (GDP) is the main method to be used to approach the contradiction between capacity supply and demand [2]. With the improvement of airport resource management technique, a modified GDP strategy, collaborative ground delay program (CDMGDP) has been implemented in some hub airports. The modified program can improve the utilization efficiency of airport resources observably [3]. The basic process of CDMGD is that the air traffic control departments allocate landing slots to airlines; then the airlines can merge and cancel the flights according to the allocated information of slots and feed the adjustment information back to the air traffic control management, and then the air traffic control departments make the final decision after receiving the feedback information from airlines [4]. Many scholars at home and abroad have done some research in CDMGDP. Hoffman et al. [5] discussed several tools applied to CDM system and a new algorithm rationbyschedule (RBS) was proposed. Vossen et al. [6] described a new allocation procedure based on FCFS in CDM strategy to schedule the arrivaldeparture flights. They established a equity allocation mechanism, but the efficiency of algorithm should be improved. Mukherjee and Hansen [7] put forward a dynamic stochastic integer programming (IP) model for the airport ground holding problem, but the equity was ignored for small airlines. Ball et al. [8] presented a new rationbydistance (RBD) algorithm showing that the equity and efficiency were improved at a certain extent. Hu and Su [9] modeled the groundholding management system, which provides the theoretical basis and method for the actual traffic management, but they ignored the influence of limited capacity for taking off flights. Ma et al. [10] designed a scheme of approaching queue and optimization schedule of flights. But the satisfaction function is just a local optimization and it is too complicated to make an adjustment in the simulation. Zhou et al. [11] proposed an effectivenessfairness (EE) standard aimed at arrival slot time allocated method by analyzing and simulating the traffic flow model in CDMGDP. The simulation results showed that the single priority decreased the total delay cost more efficiently than the double priority and the fair factors were also taken into account. Zhou et al. [12] proposed an evaluation function which was used to evaluate the priority of delay cost coefficient on the basis of existing slot allocation algorithm. The method used in this paper was more flexible and effective compared with the traditional method on the total delay cost and equity, but it lacked the flexibility in slot time allocation for arrival flights. Zhang and Hu [13] proposed a multiobjective optimization model based on the principles of effectivenessefficiencyequity tradeoffs. But the research lacked further research into the slot time reassignment in CDM GDP based on the real information of aircraft. Zhan et al. [14] adopted the Ant Colony Algorithm (ACA) and receding horizon control (RHC) to optimize the scheduling from the robustness and effectiveness of the queue model on arrival flights. But the instability of the algorithm should be improved and the real scheduling of multirunway should be taken into account. Andrea D’Ariano et al. [15] studied the problem of flight sorting at congested airports. The research regarded the flight scheduling problem as an extension of workshop scheduling problem, but the optimized results were suboptimal feasible solutions. Helmke et al. [16] presented an integrated approach to solve mixedmode runway scheduling problem by using mixedinteger program techniques. However, further researches into the flight scheduling problem under multirunway with mixed operation were required. Samà et al. [17] presented the realtime scheduling flight in order to reschedule the flight based on receding horizon control strategy and conflict detection. However, the optimization approaches requiring frequent retiming and rerouting in consecutive time horizons decreased the scheduling robustness. Hancerliogullari et al. [18] researched into the aircraft sequencing problem (ASP) under multirunway with mixed operation mode. They put greedy algorithm to simulate the model. However, the fairness among airlines was not taken into consideration.
All in all, though the optimization results of most researches at home and abroad satisfied the scheduling of flight, the studies pay little attention to realtime flight information. Most of the studies consider collaborative ground delays program of approach flight without analyzing departure flight. In fact, if the delay of departure flight is dealt with unreasonably, it can lead to unfairness among arrivaldeparture flights and the increase of flight delay. In the paper, a multiobjective optimization model is established based on multirunway arrivaldeparture flight. The maximum delay of departure flight is limited to ensure the fairness between arrivaldeparture flights according to realtime flight information. Meanwhile, a fair runway slot allocation mechanism is established with the objective of minimizing the cost caused by airline delays. As Ant Colony Algorithm has unique advantages in continuous dynamic optimization, the paper introduces it to simulate and validate the model with the expectation of reducing the loss of airline delays, improving runway utilization, and ensuring the fairness among airlines.
2. Model
2.1. Description
The queue of arrivaldeparture flights in multirunway airport is a continuous dynamic process, and it changes with the realtime information of flights. The schedule optimization of arrivaldeparture flight in multirunway airport can be described as follows: within a time window, a number of flights belonging to different airlines are waiting for landing or taking off. The managers in airport should make a reasonable allocation schedule (such as arrivaldeparture time, sequence, and operation runway) for all flights to minimize the total delay cost in the study period under the condition of safe operation of flight and airport resources and to balance the cost of total delay among airlines. The paper selects the research time period in rush hour in hub airport to study the airport surface operation. After modeling and simulation, the results can be applied in the management of airport surface operation in any type of airports.
2.2. Assumptions
(1)The parallel runways studied in the paper operate independently.(2)In the research time period, the runway capacity cannot meet the demands of flights.(3)All the arrival flights do not delay when they are in the takeoff airport, and they can arrive at the destination terminal aerospace studied on time.(4)The basic information (such as flight plans and other information of all flights) within the studied period is known.(5)Each arrivaldeparture flight can only be assigned to one time slot in the studied period.
2.3. Definition
: Set of arrival flight, : Set of departure flight, : Set of arrivaldeparture flights, : Set of airlines, : Set of arrivaldeparture flights slots, : Set of runways in the airport, Flight belongs to airline : Total flight delay cost : Total flight delay cost of airline : The sum of absolute deviation of flight delay cost : The end time of flight taking off from or landing on runway after optimization : The original time of flight taking off from or landing on at runway after optimization : The minimum safety interval of continuous landing on runway : The minimum safety interval of continuous taking off on runway : The minimum time interval when a departure flight follows an arrival flight on runway : The minimum time interval when an arrival flight follows a departure flight on runway : The maximum delay time when an arrival flight lands in advance compared to scheduled time : The maximum delay time when an arrival flight lands later than the scheduled time : The maximum delay time when a departure flight takes off in advance compared to the scheduled time : The maximum delay time when a departure flight takes off later than the scheduled time : The estimated arrival time of flight which belongs to : The estimated departure time of flight which belongs to : The actual arrival time of flight belonging to after optimization schedule : The actual departure time of flight belonging to after optimization schedule : The unit time delay cost of arrival flight belonging to after optimization schedule : The unit time delay cost of departure flight belonging to after optimization schedule.
2.4. Objective Function
Two objectives are taken into consideration in the paper: total delay cost and the fairness among airlines. The total delay cost of all flights is used to reduce the light delay and improve the utilization of runway. The fairness is used to balance the equity among all the airlines and to protect the benefit of small airlines. So a multiobjective function based on it is modeled.
2.4.1. The Objective Function of Delay Cost
Different wake vortex separations between different arrivaldeparture flights are different according to the types of aircraft. Therefore, we can improve the capacity of runway and reduce total delay time by adjusting the arrivaldeparture order of all flights. ICAO aircraft wake turbulence separation criteria are specified in Table 1.

The optimized target of delay cost in the paper is to minimize the total delay of all arrivaldeparture flights, which is based on improving the capacity of runway. The objective function of delay cost can be described as follows:
2.4.2. The Objective Function of Fairness
Delay cost is related to the aircraft type. In general, large airlines are preferred to small aircraft types. If the research only takes the delay cost as a single function, it is likely to lead to serious unfairness among airlines, especially to small airlines. Therefore, absolute deviation of delay cost is introduced to ensure the fairness among airlines.
Definition of standard flight: assume some type of flight to be a standard flight and all other flights can be transformed into it according to aircraft type and delay cost. For example, if we take a large aircraft as a standard flight and the number equals 1, then a light aircraft may be transformed as 0.6 and a heavy aircraft as 1.8. If a standard flight is denoted by , is defined as the number of standard flights after transformation from flight ; then the relation expression can be demonstrated as follows: In order to ensure fairness among airlines, the researchers first transform all flights belonging to different airlines into standard flights. Then the researchers can calculate the average delay cost of standard flight by using the total delay cost of all flights. In the same way, the researcher can get average delay cost of each airline and the total absolute deviation of delay cost. It is obvious that the lower the total absolute deviation is, the more fairness we can balance among airlines. The fairness optimization objective function of airlines is as follows: where is the average delay cost of flights belonging to airline ; is the average delay cost of all flights.
2.5. Constraints
Consider An arrival flight occupies one slot resource and each slot resource can be assigned to a certain arrival flight. Constraint (4) is the constraint of slot resource allocation for arrival flights: Similarly, a departure flight occupies one slot resource and each slot resource can be assigned to a certain departure flight. Constraint (5) is the constraint of slot resource allocation for departure flights: For the safety of flight operation, each flight can only occupy one runway and only one or none aircraft can occupy the runway at the same time. A constraint of slot resource allocation of runways is expressed as constraint (6): According to constraint (6), at a certain time, the number of flights occupying runways will not be greater than the number of runways. Meanwhile, in order to ensure safety in a certain slot, only one flight is arriving or leaving at runway . The runway resource controlling constraint is shown as constraint (7): To ensure the safety of arrival flights, the actual arrival time should meet the maximum delay time ( and ) after optimization. Constraint (8) is to ensure the time of arrival flight: Constraint (9) is to ensure the time of departure flight. Extending the departure flight delay can decrease the service level in most airports and it is necessary to limit the amount of delay time. Like the constraint for arrival flights, the actual departure time should meet the maximum delay time ( and ) after optimization. Consider The minimum safety interval in different conditions which is related to the order of arrivaldeparture flight should be taken into account. Constraint (10) is the minimum flight safety interval constraint of runway .
3. Ant Colony Algorithm Design
Ant Colony Algorithm (ACA) is a metaheuristic algorithm, which uses a heuristic method to search for the space that may be related to feasible solutions. In ACA, the ant can select the path comprehensively based on pheromones and heuristic factors of the environment. It can release the pheromones after traveling the path of the network. In the algorithm, an individual ant can identify and release all the pheromones. All pheromones from the ant colony are used to complete the whole and complex optimization process. ACA has the characteristics of selforganization and distributed computing. Therefore it is able to make a global search. It can effectively avoid local solutions to an extent. Meanwhile, ACA can get the optimized solution faster than other traditional algorithms.
3.1. Algorithm Description
3.1.1. Single Runway Flight Scheduling of ACA
Single runway scheduling problem can be transformed to a TSP which takes each flight as a node and the interval time between flights as the path length of nodes. The problem can be solved by traditional ACA, and the results are satisfactory. The model is built as follows: the node in the network is the element of flight set and the distance between nodes and is the time interval. When the algorithm begins, the ant heads from a virtual starting node . The starting node is set up to ensure that all the ants in ACA can start at the same node. The ant traverses all the nodes of network, so a flight sequence is obtained. We can calculate wait time and delay cost in the queue according to the sequence. The ACA for single runway flight scheduling is shown in Figure 1.
3.1.2. Multirunway Flight Scheduling of ACA
In order to adapt to the multirunway flight scheduling model, the ACA for single runway model should be modified. In the multirunway flight scheduling model, a node may contain several subnodes . The distance between the subnode of and the subnode of node is the minimum safety interval. The ant heads from a virtual starting node and travels all nodes of the network. When the ant arrives at a node, it selects a subnode to get a sequence contained runway number. The ACA for multirunway flight scheduling is shown in Figure 2.
3.1.3. The State Transition Equation
The amount of information on each path and the heuristic information can decide the transition direction of ant . It records the traveled nodes by search table . is denoted as the state transition probability of ant changing direction from the subnode to the subnode at time . Formula (11) is as follows:
is the pheromone concentration on path at time . means that the ant can choose the node which is the node that never traversed next step. is the pheromone heuristic factor, which decides how the pheromones have impact on path choosing; is the expected heuristic factor, which decides the degree of attention of visibility when ants make a choice.
is an expected factor, which is evaluated as where is the minimum safety interval of the flight and former flight .
3.1.4. The Updating Strategy of Pheromone
Updating pheromone of all nodes is needed when all the iterations are completed. With the increasing of pheromone concentration, the residual pheromone evaporates in proportion. In order to get better optimization results, only the best ant can release pheromone of iteration. Therefore, the pheromone updating can be adjusted as the following rules: where is the volatilization coefficient of pheromone; is the amount of pheromone; is the total incremental of pheromone in this circulation of node .
3.2. The Design of ACA
The design of ACA in the simulation is as follows.
Step 1. Set parameters. We set , , , and .
Step 2. Get the flight information. We get flight information (including the type of flight, the estimated time of arrival or departure, and so forth) and other known data by reading the files.
Step 3. Initialize the pheromone and expectations of paths of solution space and empty the tabu list.
Step 4. Set ( is the iteration). Generate the initial ants on the virtual node .
Step 5. Ants select a node orderly according to Formula (11).
Step 6. If all ants complete a traversal, turn to Step 7, otherwise turn to Step 5.
Step 7. If the searching results of all ants meet the constraints, then reduce the pheromone increment when updating pheromone.
Step 8. Calculate the target value of all ants and record the best ant solutions.
Step 9. Update the pheromone of each node according to Formula (13).
Step 10. If , without stagnating, delete the ant and set ; then reset the data and turn to Step 5; otherwise output the optimal results and the calculation is over.
4. Simulation and Vertification
In the paper, we put C program to use to simulate flight schedule problem on multirunway with mixed operation mode. The core algorithm of the program is the ACA design. In order to achieve the objectives of delay cost and fairness, we first take the delay cost as the objective and we take fairness as a constraint and we can get the flight sequence with minimal delay cost. After that, we set fairness as the objective and the delay cost as a constraint, and we get a flight sequence with the best fairness. Finally, the initial flight sequence and these two flight sequences are compared. A peak hour is selected from a certain large airport of China and the authors chose flights in the busiest 15 minutes from that peak hour. Two parallel runways run independently. There are 38 flights (belonging to 7 airlines) to be scheduled. After optimization, the initial flight information and optimized results are shown in Table 2.

In the paper, some real operation data is selected from a certain hub airport and ACA is designed to solve the model. As shown in Table 2, flight delay is serious due to unreasonable slot assignment; it can lead to unfair competition among airlines before optimization. The total delay cost of three types of flight sequences is listed in Table 3 and Figure 3.

The initial sequence is the initial flight sequence before optimization. The minimal delay cost and fairness among airlines are solved according to object function (1) and function (2). The minimal delay sequence means that we take the objective function (1) as the main objective function and objective function (2) as a constraint in optimization. The best fairness sequence means that we set objective function (2) as the main objective function and the objective function (1) as a constraint in optimization. Then the paper contrasts the three flight data among airlines.
The detailed information of delay cost of standard flight is listed in Table 4 and Figure 4. After transforming flights into standard flights, we can compare delay cost and fairness directly.
 
Note: H_{1 }to H_{7} refer to standard flights in Table 4. 
From Figures 3 and 4, we can draw the conclusion that the sequence of minimal delay cost can decrease the delay cost obviously after optimization. The sequence of best fairness can improve the fairness among all the airlines obviously. But from Figure 4, we can see that the fairness among 7 airlines decreases when delay cost is minimal; the delay cost of 7 airlines improves obviously when the fairness is best.
In order to reduce the delay cost of airlines and increase the fairness among airlines, we view the minimum delay cost as objective and control the range of flight delay cost variation. The simulation steps are as follows.(1)First, the fairness is not taken into account and we get the minimum delay cost. We statistic the total delay and the delay cost deviation.(2)Then limit the range of airline flight delay cost deviation by a large number of data. We get the simulated data of total delay and delay cost deviation of flights in the cases , , , , , respectively.(3)Finally, we make an analysis of simulation data and study the relationship between the delay cost and fairness. The relation curve of delay cost and delay cost deviation is shown in Figure 5.
As shown in Figure 5, there is a trend relationship between delay cost and fairness. When the delay cost decreases, the fairness among airlines is not satisfied. Reducing the delay deviation of flights may lead to the increase in total delay cost. In order to make balance of the relationship between the delay cost and fairness and make a better sequence of flights, we can make flight sequence by limiting the delay deviation of flights. So it not only reduces the total delay cost of airlines but also takes care of the fairness among airlines.
From the simulation results, we can get a number of flight sequences by controlling the sum of absolute deviation of flight delay cost. Five sequences (minimal delay cost, optimization 1, optimization 2, optimization 3, and the best fairness) to make a comparison of total delay cost of airlines are shown in Table 5 and Figure 6. Table 6 and Figure 7 show the same comparison by transforming flights into standard flights.

 
Note: H_{1 }to H_{7} refer to standard flights in Table 6. 
The total delay cost and the fairness among airlines have been improved obviously after optimization. In actual operation, the decision makers can get several optimized flight sequences by controlling the range of flight delay cost deviation and by selecting a preferred one according to realtime information.
In Table 7, we list the optimal flight sequence in which both the delay cost and fairness are acceptable. The optimized results show that the total delay cost reduces greatly and the fairness is also acceptable when compared with the initial flight sequence. The standard flight delay cost deviation of initial flight sequence and optimized flight sequence is calculated based on Table 7. The results are shown in Table 8 and Figure 8. The histogram shows the contrast of flight delay deviation between the initial flight sequence and the optimized flight sequence. From the histogram we can find that the delay cost of airlines declines 42.22% at least after optimization. The sum of delay deviation declines 38.64%. So the schedule model and solution have not only reduced the total delay cost significantly but also ensured the fairness among all the airlines.
 
Total: the total cost of the initial flight sequence delay is 343446.5/CNY; the total cost of optimized flight sequence delay is 102681.0/CNY. 

5. Conclusions
In the paper, a mixed multirunway operation flight scheduling optimization model based on multiobjective is proposed. Two objectives are considered: the total delay cost and fairness among airlines are two objective functions. The ACA is introduced to solve the model. The simulation results show that the total delay cost decreases significantly and the fairness among airlines is also acceptable. Meanwhile, ACA used in the paper solves the model with great efficiency.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China, the Civil Aviation Administration of China (no. U1333117), China Postdoctoral Science Foundation (no. 2012M511275), and the Fundamental Research Funds for the Central Universities (no. NS2013067).
References
 X. Hu and E. di Paolo, “Binaryrepresentationbased genetic algorithm for aircraft arrival sequencing and scheduling,” IEEE Transactions on Intelligent Transportation Systems, vol. 9, no. 2, pp. 301–310, 2008. View at: Publisher Site  Google Scholar
 M. C. Wambsganss, “Collaborative decision making in air traffic management,” in New Concepts and Methods in Air Traffic Management, pp. 1–15, Springer, Berlin, Germany, 2001. View at: Google Scholar
 M. Wambsganss, “Collaborative decision making through dynamic information transfer,” Air Traffic Control Quarterly, vol. 4, no. 2, pp. 107–123, 1997. View at: Google Scholar
 K. Chang, K. Howard, R. Oiesen, L. Shisler, M. Tanino, and M. C. Wambsganss, “Enhancements to the FAA grounddelay program under collaborative decision making,” Interfaces, vol. 31, no. 1, pp. 57–76, 2001. View at: Google Scholar
 R. Hoffman, W. Hall, M. Ball, A. R. Odoni, and M. Wambsganss, “Collaborative decision making in air traffic flow management,” NEXTOR Research Report RR992, UC Berkley, 1999. View at: Google Scholar
 T. Vossen, M. Ball, R. Hoffman et al., “A general approach to equity in traffic flow management and its application to mitigating exemption bias in ground delay programs,” Air Traffic Control Quarterly, vol. 11, no. 4, pp. 277–292, 2003. View at: Google Scholar
 A. Mukherjee and M. Hansen, “A dynamic stochastic model for the single airport ground holding problem,” Transportation Science, vol. 41, no. 4, pp. 444–456, 2007. View at: Publisher Site  Google Scholar
 M. O. Ball, R. Hoffman, and A. Mukherjee, “Ground delay program planning under uncertainty based on the rationbydistance principle,” Transportation Science, vol. 44, no. 1, pp. 1–14, 2010. View at: Publisher Site  Google Scholar
 M. H. Hu and L. G. Su, “Modeling research in air traffic flow control system,” Journal of Nanjing University Aeronautics & Astronautics, vol. 32, no. 5, pp. 586–590, 2000. View at: Google Scholar
 Z. Ma, D. Cui, and C. Chen, “Sequencing and optimal scheduling for approach control of air traffic,” Journal of Tsinghua University, vol. 44, no. 1, pp. 122–125, 2004. View at: Google Scholar
 Q. Zhou, J. Zhang, and X. J. Zhang, “Equity and effectiveness study of airport flow management model,” China Science and Technology Information, vol. 4, pp. 126–116, 2005. View at: Google Scholar
 Q. Zhou, X. Zhang, and Z. Liu, “Slots allocation in CDM GDP,” Journal of Beijing University of Aeronautics and Astronautics, vol. 32, no. 9, pp. 1043–1045, 2006. View at: Google Scholar
 H. H. Zhang and M. H. Hu, “Multiobjection optimization allocation of aircraft landing slot in CDM GDP,” Journal of Systems & Management, vol. 18, no. 3, pp. 302–308, 2009. View at: Google Scholar
 Z. Zhan, J. Zhang, Y. Li et al., “An efficient ant colony system based on receding horizon control for the aircraft arrival sequencing and scheduling problem,” IEEE Transactions on Intelligent Transportation Systems, vol. 11, no. 2, pp. 399–412, 2010. View at: Publisher Site  Google Scholar
 A. D'Ariano, P. D'Urgolo, D. Pacciarelli, and M. Pranzo, “Optimal sequencing of aircrafts takeoff and landing at a busy airport,” in Proceeding of the 13th International IEEE Conference on Intelligent Transportation Systems (ITSC 10), pp. 1569–1574, Funchal, Portugal, September 2010. View at: Publisher Site  Google Scholar
 H. Helmke, O. Gluchshenko, A. Martin, A. Peter, S. Pokutta, and U. Siebert, “Optimal mixedmode runway scheduling—mixedinteger programming for ATC scheduling,” in Proceedings of the IEEE/AIAA 30th Digital Avionics Systems Conference (DASC '11), pp. C41–C413, IEEE, Seattle, Wash, USA, October 2011. View at: Publisher Site  Google Scholar
 M. Samà, A. D'Ariano, and D. Pacciarelli, “Optimal aircraft traffic flow management at a terminal control area during disturbances,” Procedia—Social and Behavioral Sciences, vol. 54, pp. 460–469, 2012. View at: Publisher Site  Google Scholar
 G. Hancerliogullari, G. Rabadi, A. H. AlSalem, and M. Kharbeche, “Greedy algorithms and metaheuristics for a multiple runway combined arrivaldeparture aircraft sequencing problem,” Journal of Air Transport Management, vol. 32, pp. 39–48, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Yu Jiang 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.