Inﬂuence of the Apron Parking Stand Management Policy on Aircraft and Ground Support Equipment (GSE) Gaseous Emissions at Airports

: The purpose of this study is to analyze the concept of a hybrid apron with a ﬁxed number of parking positions considering the management model inﬂuence for the average delay per aircraft and the gaseous emissions generated by aircraft and ground support equipment (GSE) altogether. The apron is studied based on two gate management models: in the ﬁrst model, the aircraft are allocated in each gate due to operational factors only; in the second model, the rules of exclusive use of each gate according to the airline are included. The emissions generated by aircraft operations and that of their GSE (produced by the service and movements on the apron) are quantiﬁed and compared in the two gate management models: operation in the standard LTO cycle of the studied aircraft, GSE emissions have a similar relation with the compared gasses (NOx and CO), ranging between 1% and 3%. Further, if it compares the emissions between support vehicles and aircraft taking only into account the in-out taxiway, the relation between both CO sources shows similar values to those of the previous comparison, whereas NOx emissions produced by GSE reach an approximately 20%. The study considers different demand conditions obtained from the average day of the peak month of Aeroparque Jorge Newbery airport. Subsequently, through the SIMMOD PLUS software, the aircraft operations are simulated. The gates assignment and the arrival timetables are used as inputs for the GSE study due to an analytical model developed by us. Once the operational dimension is characterized and evaluated, the necessary data to quantify the gaseous emissions from the sources (Aircraft-GSE), based on the International Civil Aviation Organization (ICAO) guidelines, is obtained.


Introduction
Airports are critical nodes in the air transport system and, as such, they face several challenges, such as sustainable development infrastructures, which must provide a balanced design between the operational and environmental dimensions.
As stated in [1], the aviation industry grows within a complex network of systems. The complexity of this network results in delays produced by many sources, including weather conditions, airspace restrictions, and airport congestion. Of these three, the main cause of delay is associated with weather conditions, which are beyond the control of the industry. In this context, one of the issues to study in depth is ground congestion and its consequences, such as the associated environmental impact, which is an aspect addressed in this study.
Particularly, ground support equipment (GSE) include, among their operations, passenger transportation from terminals to aircraft and vice versa, loading and unloading of goods, aircraft energy and fuel supply, crew transportation, and any maneuver needed to a capacity envelope for an apron with 11 positions, in which the stand configuration (domestic and mixed) is managed and determines the capacity for different market allocations (% of domestic flights), and the turnaround time. Ref. [8] analyses the air level of the airport through the runway system and the apron-gates complex, assuming a fully developed taxiway system. As some authors state, there are few analytical models taking into account the apron capacity. In another study, the same authors [8] pose an airport configuration with 30 stands (22 in contact and 8 remote) with four demand configurations (characterized by the light, medium and heavy aircraft percentage) depending on the type of airport under consideration (HUB, MIXed Hub or O/D), and also propose a demand cyclical behavior, which varies pursuant to the aircraft type and airport type. The results lead to the conclusion that the characteristics associated with the type of traffic and the limiting factors cannot be disregarded.
About the problem of optimization in parking stand assignment, there are several publications focused on minimizing the distance covered by the passenger in the terminal area, summarized in [9]. Following this, Ref. [10] proposes an algorithm of stand and taxiway assignment, and compares the results obtained from SIMMOD simulations. On the other hand, Ref. [11] proposes a simulation model of gate assignment to predict the operational time in congested aprons. The implementation of this type of models in real time permits a 50% reduction of ramp times. As the author states, to identify physical conflicts, the first thing to do is studying the ramp geometry and physical limitations to understand its dynamic. Physical and operational interactions and restrictions include the aircraft taxiway preventing other aircraft to do a pushback, a narrow apron or one limited in the number of entry and exit points, and long aircraft queues waiting in the ramp for the taxiway clearance. Such restrictions depend on the ramp geometry.
The importance of the apron and taxiway system is also pointed out by [12], who explains that, in the most important airports in the United States, aircraft traffic congestion takes place in the apron and surrounding areas, leading to inefficiencies in the movement area, which are then reflected in the reduction of performance time and the increase of taxiing delays. Consequently, there is an increase in the operational costs and the environmental impact due to the emissions produced by fuel combustion. In this context, Ref. [13] analyses the emissions of the Istanbul airport and identifies the taxiway and apron as the components with the most influence over emissions, besides allocating a monetary cost based on European Union (EU) data.
There are studies that have taken into account multiple variables, among which we can mention those of [14], who presents a methodology for the performance analysis, shaping and evaluation of soil use made by airports and states where the interaction between physical, operational, environmental, economic, and social variables play a key role. Physical, operational, and environmental characteristics are directly addressed in this study. Ref. [15] analyses potential problems and limitations regarding the capacity of different apron designs, with the aim of identifying a better arrangement according to the proposed cost-benefit function. Results indicate that it is possible to identify the configuration favoring the reduction of turnaround times and, therefore, the airline costs, and a second configuration which allows strong response conditions and also a higher turnaround time. Afterwards, Ref. [16] uses a model based on a "virtual cycle", which consists of two independent models, feeding each other to optimize apron configuration. Through this "virtual cycle", the groups and number of required vehicles can be determined according to the two proposed apron configurations (independent and supported stands). Following the studies that take into account the operational and environmental aspects, there appears [17], who presents a reference framework for the evaluation of airport project designs based on sustainability performance. As pointed out, environmental and operational capability can be maximized through a long-term planning guaranteeing an effective environmental management, which compensates growth by introducing infrastructure strategies from a technological and operational point of view. This problem is not insignificant, since environmental restrictions affect a huge number of airports. For instance, Ref. [18] indicates that, in the case of Europe, 70% of airports are affected by environmental restrictions. Among the operational problems, it highlights the possibility of answering demand conditions, and the satisfaction of such requirement, which involves taking into account factors such as capability, environmental impact, territorial competitiveness, and operational security conditions.

Objectives of the Research
The aim of this study is to present the results of the application of the methodology developed to quantify the emissions of aircraft on the ground and GSE according to different gate assignation policies, service time discretization and characterization of the circulation process depending on the measures made in the Jorge Newbery (SABE) Airport.

Methodology
This study models different scenarios characterized according to gate assignment policy and demand condition, depending on the number of operations per hour, over which gaseous emissions produced by aircraft and GSE were studied and quantified.
There are three complementary and interrelated stages in the development of this methodology (see Figure 1):

•
In the first stage, SIMMOD Plus tool (called Model 1) is used to simulate operations of aircraft on the ground, particularly on the apron. SIMMOD software is a simulation program of discreet events which allows us to study the flight field dynamic, airspace routes, taxiing operations, and departure queues sequences, among other events related to the system capacity and the delay associated with it. The software allows us to quantify the delay based on the conditional operational rules included in the model, which specify the actions to be carried out by the simulation based on the system state.

•
In this study, delay is measured in minutes and defined as "the time necessary to meet the requirements of a second aircraft when, at the same time, two aircraft request the same service".

•
The second stage corresponds to the implementation of Model 2, aiming at the study and quantification of GSE emissions produced by apron circulation and service. To do this, the covered distances are estimated through the assignment of gates. Aircraft service stage times are discretized, and other factors associated with GSE operations are identified, such as power, loading, damages, and emission factors. • Finally, the third stage follows the guidelines provided by [6] and widely applied for the study of the polluting input, as observed in [5,19] with the purpose of obtaining aircraft emissions.
As observed in the methodological diagram, from the output of Model 1 produced by the 20 iterations made in each scenario, characteristic mean values and standard deviation of taxiway times and delays for the analysis conditions are produced. From the total amount of iterations, one representative of mean values was selected to identify the used stand assignment sequence.
Arrivals, aircraft, and gate assignments are the outgoing data of Model 1 and the entrance to the following one (Model 2), whose objective is quantifying emissions produced by the operation of GSE vehicles apron circulation and service. Both situations are analyzed in parallel to provide the corresponding values to the proposed equation for each case revised by the authors. Values also depend on aircraft support policies determined by the types of flight (point-to-point, with scale or low cost) and the type of service associated with parking stands (use of passenger boarding bridge or remote stands).
On the other hand, Model 3 uses taxiway average and delay aircraft times produced in the arrival and departure stage and provided by Model 1 as input, since they allow us to estimate CO and NOx emissions.
Hereafter, it presents the measuring models and the characteristics of each of the stages above mentioned. As observed in the methodological diagram, from the output of Model 1 produced by the 20 iterations made in each scenario, characteristic mean values and standard deviation of taxiway times and delays for the analysis conditions are produced. From the total amount of iterations, one representative of mean values was selected to identify the used stand assignment sequence.
Arrivals, aircraft, and gate assignments are the outgoing data of Model 1 and the entrance to the following one (Model 2), whose objective is quantifying emissions produced by the operation of GSE vehicles apron circulation and service. Both situations are analyzed in parallel to provide the corresponding values to the proposed equation for each case revised by the authors. Values also depend on aircraft support policies determined by the types of flight (point-to-point, with scale or low cost) and the type of service associated with parking stands (use of passenger boarding bridge or remote stands).
On the other hand, Model 3 uses taxiway average and delay aircraft times produced in the arrival and departure stage and provided by Model 1 as input, since they allow us to estimate CO and NOx emissions.

Infrastructure details
The modelled infrastructure consists of a configuration of a simple runway with a taxiway system adequately developed for the runway system. On the other hand, since the main object of study lays upon the apron complex-parking stands, this is presented including three contact points with the taxiway system (see Figures 2 and 3).

Infrastructure details
The modelled infrastructure consists of a configuration of a simple runway with a taxiway system adequately developed for the runway system. On the other hand, since the main object of study lays upon the apron complex-parking stands, this is presented including three contact points with the taxiway system (see Figures 2 and 3). The taxiway on the apron has positions on both sides, giving rise to a hybrid configuration. The infrastructure has been shaped and characterized with 30 C-type stands. Following the ICAO classification standards, stands G-3 to G-12 are those in contact with independence entrance and supported departure maneuvers for pushback, while the remaining positions are operated with an independent entrance and departure maneuver.
The delay is calculated as the difference between the entry time and the exit time of the queue for every process the aircraft carries out. The construction of the model involves the definition of taxiing rules in the apron entry and exit. For instance, nodes N1 to N6 are specified for the decision-making process. In node N1, it evaluates if the primary taxiway corresponding to that gate is conflict-free.

Model 1-Simulation Tool
Infrastructure details The modelled infrastructure consists of a configuration of a simple runway with a taxiway system adequately developed for the runway system. On the other hand, since the main object of study lays upon the apron complex-parking stands, this is presented including three contact points with the taxiway system (see Figures 2 and 3). The taxiway on the apron has positions on both sides, giving rise to a hybrid configuration. The infrastructure has been shaped and characterized with 30 C-type stands. Following the ICAO classification standards, stands G-3 to G-12 are those in contact with independence entrance and supported departure maneuvers for pushback, while the remaining positions are operated with an independent entrance and departure maneuver.
The delay is calculated as the difference between the entry time and the exit time of the queue for every process the aircraft carries out. The construction of the model involves the definition of taxiing rules in the apron entry and exit. For instance, nodes N1 to N6 are specified for the decision-making process. In node N1, it evaluates if the primary taxiway corresponding to that gate is conflict-free. The taxiway on the apron has positions on both sides, giving rise to a hybrid configuration. The infrastructure has been shaped and characterized with 30 C-type stands. Following the ICAO classification standards, stands G-3 to G-12 are those in contact with independence entrance and supported departure maneuvers for pushback, while the remaining positions are operated with an independent entrance and departure maneuver.
The delay is calculated as the difference between the entry time and the exit time of the queue for every process the aircraft carries out.
The construction of the model involves the definition of taxiing rules in the apron entry and exit. For instance, nodes N1 to N6 are specified for the decision-making process. In node N1, it evaluates if the primary taxiway corresponding to that gate is conflict-free. If it is, the taxiway is the next step; if it is not, conflict resolution rules are applied on the intersections. Some conflict resolution rules are exemplified in the Appendix A.

Model 2-Integrated GSE Model
Measurement consists of two stages: the first one studies vehicle circulation on the apron and the second one consists of the sum of the polluting input due to service time discretization according to GSE equipment waiting time, connection, service (task completion) and disconnection.
From the "sophisticated" model for emission measurement provided by the ICAO, which considers five factors, power, damages, loading and operational time factors were adjusted to strengthen the model and include GSE circulation (a process that was not considered in the base model) and represent in a more precise way the special features of the operations according to the GSE service.
Next, the adapted model corresponding to the emissions produced by GSE service is presented, in which discretization of loading and time factors according to the period (wait, connection, service, and disconnection) for both aircraft unloading and loading are proposed.
Likewise, the standardized calculation of emissions produced by circulation is presented:

Model 3-Aircraft Emission
Measurement of aircraft and GSE emissions is based on the method proposed by [6]. Emission factors associated with the operational mode of the aircraft are obtained from [20] database.
The standardized calculation for aircraft emissions is given by In this study, TIM refers to the taxiway time and delays produced in the apron complex and its associated taxiway. Used times correspond to the average values obtained from the 20 iterations made on each scenario analyzed on Model 1.

Demand description
Several different demand profiles, using SABE peak hour in 2017 as base line, were proposed. These used the IATA methodology to take real conditions from the arrival and departure scheduled time. The identified hour demand was reproduced for the previous and subsequent hours to analyze a representative time frame in a real operational scenario. Afterwards, on the base line (36 arrival movements in three hours), changes were made to generate different demand scenarios (operations/hour) and proportion scenarios of aircraft belonging to an airline with exclusive stand use. Figure 4 shows the average profiles of arrivals and departures obtained using Model 1, resulting from the 20 iterations made on each scenario. As a result, the profile of the demand shows decimal values, since the proposed programming took into account an exact hour for each flight (entry data). However, the influence of random variables included in the model through the probability distributions presented on Appendix B result in some operation being carried out even after 13:00. The main goal of Figure 4 is to represent the level of activity and interaction between arrival and departure movements.
Regarding service times, the occupation times of the aircraft in the parking spot in the airport were registered for a particular time frame. According to the registries, the probability distribution presented on Appendix B and assigned to the flights on Model 1 was drafted. This distribution of turnaround time shows a minimal time of 50 min and a maximum of 150 min. On Appendix B, it presents the cumulative frequency distribution of the process of the turnaround, pushback and arrival and departure latency. in some operation being carried out even after 13:00. The main goal of Figure 4 is to represent the level of activity and interaction between arrival and departure movements. Regarding service times, the occupation times of the aircraft in the parking spot in the airport were registered for a particular time frame. According to the registries, the probability distribution presented on Appendix B and assigned to the flights on Model 1 was drafted. This distribution of turnaround time shows a minimal time of 50 min and a maximum of 150 min. On Appendix B, it presents the cumulative frequency distribution of the process of the turnaround, pushback and arrival and departure latency.

Description of Inputs
Simulated scenarios Scenarios are characterized by the parking stand assignation policy and the demand conditions which depend on the number of operations per hour.
Such policy responds to the exclusive use or not of the stand by an airline. In this study, 7 configurations of the use of stands were modelled and, on each configuration, 4 demand conditions were simulated, taking into account 3 airline participation percentages with parking stand use priority, resulting in a total of 84 scenarios subjected to 20 iterations each.
There follows a list of scenario combinations depending on the variables defining them: demand factor (DF), percentage of stands with exclusive use (%Se), and percentage of aircraft with access to exclusive stands (% Aircraft with priority-Aes) according to Table 1.

Description of Inputs
Simulated scenarios Scenarios are characterized by the parking stand assignation policy and the demand conditions which depend on the number of operations per hour.
Such policy responds to the exclusive use or not of the stand by an airline. In this study, 7 configurations of the use of stands were modelled and, on each configuration, 4 demand conditions were simulated, taking into account 3 airline participation percentages with parking stand use priority, resulting in a total of 84 scenarios subjected to 20 iterations each.
There follows a list of scenario combinations depending on the variables defining them: demand factor (DF), percentage of stands with exclusive use (%Se), and percentage of aircraft with access to exclusive stands (% Aircraft with priority-Aes) according to Table 1.

GSE Characterization
The amount and characteristics of the support vehicles depend on the aircraft type of operation. The following possible configurations were identified: point-to point service with boarding bridge and remote stand; partial service (connecting flight) with boarding bridge and remote stand; and Low Cost (LC). Connecting flight is identified because, in the destination place, a full exchange of passengers, load, and fuel is produced; in a transit flight, passenger, load, and fuel exchanges take place partially and, unlike what occurs in the former case, there is no cabin or wastewater cleaning. On the other hand, on LC flights, speed, efficiency, and precision are important in the operation to minimize response times and ground handling costs. Regarding the case under study of this type of configuration, the ground power unit (GPU), catering, pushback, passenger stair, passenger bus, and cabin cleaning vehicle are not taken into account [21].
It is necessary to identify the characteristics of GSE for the traffic and service stages, which is to say, the power parameters, loading factor, emission and damage factor, traffic distances, and discretized times for the service process. In the Appendix C, the considered factors are presented.
Theoretical service times associated with each aircraft can be found in the corresponding Airport Planning Manual. However, the proposed emission model takes into account the characterization of the aircraft support process according to the times measured in Jorge Newbery Airport. The values have been obtained by conducting measurements for procedures during the service of the proposed GSE. "Waiting", "connection", "service" and "disconnection" are used for both unloading and loading of passengers, food, luggage, among others. As stated in Equation (1) and pointed out in Table 2, the time average is discretized for the waiting, connection, service, and disconnection. Apron traffic: transfer distance Transfer distances depend on the geometrical arrangement of the passenger terminal apron and on the location of parking defined areas for GSE equipment. Apron characterization is not only associated with geometric configuration and the identification of aircraft parking spots, but also with identification of GSE parking spots.
Airport locations have different areas intended for parking and circulation of ground support vehicles. Identifying the arrangement of these parking spots and circulation permitted zones will be of great importance since it will allow the determination of pathways associated with each of the vehicles. Said characterization will be beneficial for estimating circulation times on apron. Therefore, some considerations have to be taken into account for its measurement in apron, according to the type of procedure.
Therefore, some considerations have to be taken into account for its measurement in apron, according to the type of procedure (see Table 3). Catering, cleaning, and water supplying vehicles are parked outside the apron [3] . The tanker has assigned parking stands in the fuel plant on the airport grounds. [1] The traffic sequence in the apron depends on the availability of each group of service vehicles. [2] It is chosen in such a way that it is best if the equipment is always moving between this point and every aircraft parking stand. [3] Vehicles have a larger loading capacity and therefore these exit the parking stand and are able to provide a service to a higher number of aircraft.
Up next, considerations for each procedure are presented: • Type A procedure: According to the equipment availability, a transfer between temporal GSE parking stands (ESA) can be simulated depending on the number of aircraft stands (Type C as they are defined on the apron configuration) that they have to go through for the next service. • Type B procedure: These vehicles always have to return to a particular operation area after providing the service to the aircraft. Therefore, it is easy to measure the transfer distance because the team goes to every aircraft parking stand and returns to its determined fixed area before moving to another aircraft. • Type C procedure: Since this equipment has a fixed parking stand, transfer distances can be estimated depending on the first arrival, along with service sequence for three extra aircraft according to the loading capacity before going back to their fixed parking stand to restock or unload wastes.
To identify the minimal amount of GSE vehicles, some Gantt diagrams have to be developed with the service time corresponding to each type of support for the amount and mixture of aircraft. The obtained service time in the corresponding Gantt diagram shows if the GSE vehicle is available or free to support another aircraft. Diagrams made for each type of service provide, as result, the minimal number of necessary vehicles for supporting aircraft in the apron on the time under study.

Model 1
In addition to generating assignment per gate for aircraft, the simulation performed on Model 1 obtained associated taxiway times and delays, differentiated in the arrival and departure process for each of the 84 scenarios. There were no significant variations among the scenarios on taxiway times, which is why the Figure 5 only shows the total delay (arrival and departure) produced in each scenery. Figure 5 shows that the delay increase occurs due to the effect of three parameters on the system performance: the demand factor (DF), the percentage of stands with exclusive use (%Se) and the percentage of aircraft with access to exclusive stands (%Aes). Further, for the same %Se conditions, the total delay of the system increases along with the DF rise.
In addition to generating assignment per gate for aircraft, the simulation performed on Model 1 obtained associated taxiway times and delays, differentiated in the arrival and departure process for each of the 84 scenarios. There were no significant variations among the scenarios on taxiway times, which is why the Figure 5 only shows the total delay (arrival and departure) produced in each scenery.  Figure 5 shows that the delay increase occurs due to the effect of three parameters on the system performance: the demand factor (DF), the percentage of stands with exclusive use (%Se) and the percentage of aircraft with access to exclusive stands (%Aes). Further, for the same %Se conditions, the total delay of the system increases along with the DF rise.
Simultaneously, it is observed that for a same DF, total delay shows a behavior with two separated characteristic areas for an amount of %Se close to 50%: the first one shows a steady behavior as %Se increases, whereas the second one varies with the DF, in which the total delay presents a strong rise inversely proportional to %Aes.
As a consequence, the apron complex capacity is noted depending on the demand characteristics and the infrastructure management model. Therefore, it seeks to identify the existence of an optimal point associated with the aircraft operation depending on the combination of different variables: DF, %Aes and %Se.
Total produced emissions may rise up to 300% due to the delay increase produced by the combination of variables that characterize the demand. In this sense, emissions increase following a behavior proportional to the relation presented hereafter, defined with the Z variable.
Variable Z is built according to the analysis of data and variables made with R Studio software. The variable is presented with the objective of estimating the relationship between land delays and the characteristics of the simulation model. Equation (8) shows the adjustment, which may be considered valid with a level of confidence of 0.95, given its pvalue and T-Student and Fisher test, as observed on Figure 6, which contrasts the results obtained on the simulation with those provided by the proposed model. Simultaneously, it is observed that for a same DF, total delay shows a behavior with two separated characteristic areas for an amount of %Se close to 50%: the first one shows a steady behavior as %Se increases, whereas the second one varies with the DF, in which the total delay presents a strong rise inversely proportional to %Aes.
As a consequence, the apron complex capacity is noted depending on the demand characteristics and the infrastructure management model. Therefore, it seeks to identify the existence of an optimal point associated with the aircraft operation depending on the combination of different variables: DF, %Aes and %Se.
Total produced emissions may rise up to 300% due to the delay increase produced by the combination of variables that characterize the demand. In this sense, emissions increase following a behavior proportional to the relation presented hereafter, defined with the Z variable. Z = %Se 2 ·(1 − %Aes)·DF Variable Z is built according to the analysis of data and variables made with R Studio software. The variable is presented with the objective of estimating the relationship between land delays and the characteristics of the simulation model. Equation (8) shows the adjustment, which may be considered valid with a level of confidence of 0.95, given its p-value and T-Student and Fisher test, as observed on Figure 6, which contrasts the results obtained on the simulation with those provided by the proposed model.

Models 2 and 3
With the purpose of demonstrating the effect of the gate assignment policy on emissions and the contribution of aircraft and GSE, 8 scenarios were prepared. In order to do this, variable DF is fixed with a value of 0.80 (below the congestion limit according to the reference values set forth by the Federal Aviation Administration (FAA). For variables %Aes and %Se, the values assigned are those presented on Table 4.
Particularly, the emission produced by the processes associated with GSE service and their circulation does not show significant variations due to the gate assignment policy, while that generated by the service is constant for the same demand conditions (aircraft and type of flight).

Models 2 and 3
With the purpose of demonstrating the effect of the gate assignment policy on emissions and the contribution of aircraft and GSE, 8 scenarios were prepared. In order to do this, variable DF is fixed with a value of 0.80 (below the congestion limit according to the reference values set forth by the Federal Aviation Administration (FAA). For variables %Aes and %Se, the values assigned are those presented on Table 4.  Scenario  %Aes  %Se  1  67  57  2  67  67  3  50  57  4  33  57  5  67  77  6  50  67  7  50  77  8 33 67 Particularly, the emission produced by the processes associated with GSE service and their circulation does not show significant variations due to the gate assignment policy, while that generated by the service is constant for the same demand conditions (aircraft and type of flight).
However, if the policy of gate assignment results in aircraft operations of another scale, GSE services will logically be modified. This is because this will substantially represent a variation in the number of vehicles required, which will cause more routes in total km, and service times will be affected. For instance, if we compare an LC flight of a narrow-body aircraft with the flight of a wide-body aircraft point to point, emissions would drastically increase. In this sense, many airports currently work proposing different measures for mitigation, such as GSE electrification, underground passages to provide services and analysis of the location of different base stations for each GSE team, to mention just a few.
Emissions produced by delays show significant variations as the gate assignment policy varies, whilst the contaminant contribution associated with aircraft movement,  However, if the policy of gate assignment results in aircraft operations of another scale, GSE services will logically be modified. This is because this will substantially represent a variation in the number of vehicles required, which will cause more routes in total km, and service times will be affected. For instance, if we compare an LC flight of a narrow-body aircraft with the flight of a wide-body aircraft point to point, emissions would drastically increase. In this sense, many airports currently work proposing different measures for mitigation, such as GSE electrification, underground passages to provide services and analysis of the location of different base stations for each GSE team, to mention just a few.
Emissions produced by delays show significant variations as the gate assignment policy varies, whilst the contaminant contribution associated with aircraft movement, even though it shows variations, is insignificant compared to that produced by delays (see Table 5 and Figures 7 and 8).   Figure 8 shows that the impact of GSE CO emissions (circulation + service) r between 4% and 12% according to the gate assignment policy. Whereas the impact emissions from aircraft associated with the delay-movement process shows a relatio that ranges from 30-59% (scenario 1) to 78-18% (scenario 8) within the emission tot Regarding NOx emissions, the performance of total emission is similar to that o However, it is noted that the impact of GSE vehicles becomes more important. Figures 9 and 10 shows that the impact of GSE NOx emissions (circulation + se ranges between 34 and 66% according to the gate assignment policy. Whereas the im   Figure 8 shows that the impact of GSE CO emissions (circulation + service) r between 4% and 12% according to the gate assignment policy. Whereas the impact emissions from aircraft associated with the delay-movement process shows a relatio that ranges from 30-59% (scenario 1) to 78-18% (scenario 8) within the emission tot Regarding NOx emissions, the performance of total emission is similar to that However, it is noted that the impact of GSE vehicles becomes more important. Figures 9 and 10 shows that the impact of GSE NOx emissions (circulation + se ranges between 34 and 66% according to the gate assignment policy. Whereas the i  Figure 8 shows that the impact of GSE CO emissions (circulation + service) ranges between 4% and 12% according to the gate assignment policy. Whereas the impact of CO emissions from aircraft associated with the delay-movement process shows a relationship that ranges from 30-59% (scenario 1) to 78-18% (scenario 8) within the emission total.
Regarding NOx emissions, the performance of total emission is similar to that of CO. However, it is noted that the impact of GSE vehicles becomes more important. Figures 9 and 10 shows that the impact of GSE NOx emissions (circulation + service) ranges between 34 and 66% according to the gate assignment policy. Whereas the impact of NOx emissions from aircraft associated with the delay-movement process shows a relationship that ranges from 12-24 (scenario 1) to 54-13% (scenario 8) within the emission total.
As mentioned before, the service and number of GSE used in a particular airport depend mainly on the aircraft, type of flight and parking stand. Having established the variables, the corresponding emissions depend on the required vehicles. Furthermore, it is possible to obtain from Model 2 the detailed emission percentage corresponding to all GSE for a point-to-point operation, as is shown in the Figure 11: of NOx emissions from aircraft associated with the delay-movement process shows lationship that ranges from 12-24 (scenario 1) to 54-13% (scenario 8) within the emi total.  As mentioned before, the service and number of GSE used in a particular ai depend mainly on the aircraft, type of flight and parking stand. Having established variables, the corresponding emissions depend on the required vehicles. Furthermo is possible to obtain from Model 2 the detailed emission percentage corresponding GSE for a point-to-point operation, as is shown in the Figure 11: of NOx emissions from aircraft associated with the delay-movement process shows lationship that ranges from 12-24 (scenario 1) to 54-13% (scenario 8) within the emis total.  As mentioned before, the service and number of GSE used in a particular air depend mainly on the aircraft, type of flight and parking stand. Having established variables, the corresponding emissions depend on the required vehicles. Furthermo is possible to obtain from Model 2 the detailed emission percentage corresponding t GSE for a point-to-point operation, as is shown in the Figure 11: Likewise, it is possible to compare GSE emissions during the service according to the type of flight and the parking stand using a particular aircraft. Figure 12 shows GSE emission quantification in the service provided to an A320-200: Likewise, it is possible to compare GSE emissions during the service according to the type of flight and the parking stand using a particular aircraft. Figure 12 shows GSE emission quantification in the service provided to an A320-200: (a) (b) Figure 11. Emissions per GSE, (a) CO, (b) NOx.
Likewise, it is possible to compare GSE emissions during the service according to the type of flight and the parking stand using a particular aircraft. Figure 12 shows GSE emission quantification in the service provided to an A320-200: On the other hand, each scenario and, as a consequence, each stand assignment policy brings forth a different requirement of vehicles and, therefore, of kilometers covered per vehicle. Table A6 shows that the number of vehicles varies in the scenario under analysis. However, it does not present a behavior proportional to the variables used in the characterization of the demand (%Aes y %Se). Furthermore, excluding the GSE associated with passenger transport, the number of vehicles required shows a stable performance with a variation range around its average value of +/− 1 vehicle (approx. +/− 10%). In such cases, it is noted that, given the number of travelling GSE, the inclusion of a new vehicle does not imply a significant increase in the required total (see Table 6). On the other hand, each scenario and, as a consequence, each stand assignment policy brings forth a different requirement of vehicles and, therefore, of kilometers covered per vehicle. Table A3 shows that the number of vehicles varies in the scenario under analysis. However, it does not present a behavior proportional to the variables used in the characterization of the demand (%Aes y %Se). Furthermore, excluding the GSE associated with passenger transport, the number of vehicles required shows a stable performance with a variation range around its average value of +/− 1 vehicle (approx. +/− 10%). In such cases, it is noted that, given the number of travelling GSE, the inclusion of a new vehicle does not imply a significant increase in the required total (see Table 6). Regarding the covered kilometers, no GSE shows a behavior directly dependent on the variables used in the demand characterization. Table 7 shows that the variation of the required kilometers may be twice as high (see GSE case for passenger transport scenario 1 vs. scenario 8), whereas, for those which include services requiring a lower amount of GSE vehicles (GPU, and Tug), it is noted that the total amount of covered kilometers shows a higher balanced performance among scenarios. If the relationship between average covered kilometers per vehicle is analyzed, it is noted that there is a significant variation between scenarios (with the exception of catering, wastewater and cleaning services). Such variation does not show a relationship between the variables used in the characterization of the demand (see Table 8): Due to its particular characteristic, the fuel service appears separately due to the fact that the distance between the fuel plant and the apron varies significantly according to the airport configuration, so in this case the adopted distance was 3.2 km (see Table 9).

Conclusions
The results cannot be extrapolated to all airports, but they serve a guide and reference to compare values and orders of magnitude.
In this context, the capacity of the apron complex-parking stands depends on the demand characteristics and the infrastructure management model. Therefore, the existence of an optimal point depends on the combination of different variables: DF, %Aes, and %Se.
Total emissions (aircraft + GSE) can increase due to the variable combination characterizing the capacity and demand, basically, because total emissions and delay are directly proportional.
CO and NOx quantification associated with aircraft was classified in emissions by movement and delay. For the conditions in which the Equation (7) shows low values, the first main emission input is related to the aircraft movement, whereas, as it increases, those which prevail are the emissions originated by delay due to apron overcrowding either for conditions related to the DF (number of operations per hour) or for conditions associated with the stand assignment policy (%Se) according to the airport operator.
Aircraft contribution to emissions has two constituents: delay and movement. Their impact shows variations according to the analyzed scenario: The delay has a CO impact ranging between 30% and 78%, whereas, in movements, there is a CO impact ranging between 59% and 18%. Regarding NOx, delay input ranges between 12% and 54% and movement input, between 24% and 13%.
As regards GSE emissions (traffic and service) related to the total emissions (including aircraft with delay and in movement), CO impact ranges between 4% and 12% for the different scenarios; whereas those of NOx range between 34% and 66%, according to the gate assignment policy. On the other hand, if only GSE operation is analyzed, the service represents 88% of CO emissions, while NOx emissions range between 88% and 98% according to the gate assignment policy. During the service, the results show that the GPU and baggage tractor emits 70% of the total adding up both on CO and NOx. This is due to the fact that the first one is the one offering the service for most time (4666 s), whereas the high emissions produced by BT are the result of high FC they have in both loading and unloading states. Moreover, in both cases it is noted that the EF is one of the highest after BL.
The point-to-point service is the one producing more emissions, followed by the service in flights with scales and finally by the one in Low Cost (LC) flights. In addition, it is noted that the services for aircraft with boarding bridges are minimally higher both in point-to-point flights and flights with scales, even though they require less vehicles.
Compared to the operation in the standard LTO cycle of the studied aircraft, GSE emissions have a similar relation with the compared gasses (NOx and CO), ranging between 1% and 3%. On the other hand, if it compares the emissions between support vehicles and aircraft taking only into account the in-out taxiway, the conclusion is that, on the one hand, the relation between both CO sources shows similar values to those of the previous comparison, whereas NOx emissions produced by GSE reach an approximately 20%.
Even though time measurements were taken at Jorge Newberry Airport, the times associated to the proposed stages are similar for services with similar features. Furthermore, the methodology can be implemented in other airports with similar geometric characteristics and operations. This is so because of each GSE operation and circulation, which, in turn, depends on the service needed to be provided to the aircraft (commercial or loading flights, number of passengers, etc.).  Figure A1, where you can see the evaluation of the benefits resulting from using a secondary taxiway or waiting for the primary taxiway to be free. The alternative minimizing the waiting times shall be chosen. If the assigned gate is taken, an alternative one shall be assigned, and the decision-making process, repeated.

Appendix A. Conflict Resolution Rules
As it can be observed in Figure A2, in case of conflicts during the taxiway, priority is given to the departing aircraft in order to free the congestion level on the apron. In case of delays in the departure queue, aircraft accommodate the delay in the gate, as long as there are free stands.
Z is a variable built according to the analysis of data and variables using R Studio software. Figure A1, where you can see the evaluation of the benefits resulting from using a secondary taxiway or waiting for the primary taxiway to be free. The alternative minimizing the waiting times shall be chosen. If the assigned gate is taken, an alternative one shall be assigned, and the decision-making process, repeated.  As it can be observed in Figure A2, in case of conflicts during the taxiway, priority is given to the departing aircraft in order to free the congestion level on the apron. In case of delays in the departure queue, aircraft accommodate the delay in the gate, as long as there are free stands. Figure A2. Taxiway conflicts.

Appendix B. Characterization of Processes-Probability Distributions
Turnaround process Figure A3 summarizes that pattern, in which, for instance, it can be observed that approximately 40% of the aircraft show a turnaround time lower than 70 min.

Appendix B. Characterization of Processes-Probability Distributions
Turnaround process Figure A3 summarizes that pattern, in which, for instance, it can be observed that approximately 40% of the aircraft show a turnaround time lower than 70 min. Figure A2. Taxiway conflicts.

Appendix B. Characterization of Processes-Probability Distributions
Turnaround process Figure A3 summarizes that pattern, in which, for instance, it can be observed that approximately 40% of the aircraft show a turnaround time lower than 70 min. Arrival latency and departure latency Probability distributions produce randomness that have an impact on the scheduled hours of the execution process, which, in turn, affect system capacity and performance. In this study, the processes of arrival latency, departure latency, and pushback have been characterized using the cumulative frequency distributions of Figures A4-A7. Arrival latency and departure latency Probability distributions produce randomness that have an impact on the scheduled hours of the execution process, which, in turn, affect system capacity and performance. In this study, the processes of arrival latency, departure latency, and pushback have been characterized using the cumulative frequency distributions of Figures A4-A7.
The arrival latency was modelled by the distribution shown in Figure A4, even though the behavior of this variable was significantly distinct from one airport to the other. The latency distribution is adapted from [22], given that our aim was to add randomness to the process. The arrival latency was modelled by the distribution shown in Figure A4, even though the behavior of this variable was significantly distinct from one airport to the other. The latency distribution is adapted from [22], given that our aim was to add randomness to the process. On the other hand, when considering the scheduled hours for departure flights, a departure latency probability was included ranging from 0 to 20 min with the distribution law observed on Figure A5 based on usual delays that might be produced by processes in the terminal or delays in apron services managed as buffer times. On the other hand, when considering the scheduled hours for departure flights, a departure latency probability was included ranging from 0 to 20 min with the distribution law observed on Figure A5 based on usual delays that might be produced by processes in the terminal or delays in apron services managed as buffer times. On the other hand, when considering the scheduled hours for departure flights, a departure latency probability was included ranging from 0 to 20 min with the distribution law observed on Figure A5 based on usual delays that might be produced by processes in the terminal or delays in apron services managed as buffer times. Figure A5. Departure latency distribution.

Pushback process
The pushback process is divided into two operations. The first one is the towing of aircraft from gate to taxiway line ( Figure A6), and the second one is the start-up of the aircraft ( Figure A7). Both processes were modelled according to established probability

Pushback process
The pushback process is divided into two operations. The first one is the towing of aircraft from gate to taxiway line ( Figure A6), and the second one is the start-up of the aircraft ( Figure A7). Both processes were modelled according to established probability distributions [23]. Later, pushback process was sampled for SABE, with the sampling results compatible with the distribution used. distributions [23]. Later, pushback process was sampled for SABE, with the sampling results compatible with the distribution used. Figure A6. Characterization of pushback time process [24]. Figure A6. Characterization of pushback time process [24]. Figure A6. Characterization of pushback time process [24].

Appendix C. GSE Factor
Maximum power, loading, or power factor Based on the field data gathering, the proposed model times have been measured, and GSE models have been characterized in order to relate their power plant and observe the engine state in the service to estimate loading factors.
Next, maximum powers and load factors used in the model are presented in Tables A1-A3: Figure A7. Characterization of dwell time process.

Appendix C. GSE Factor
Maximum power, loading, or power factor Based on the field data gathering, the proposed model times have been measured, and GSE models have been characterized in order to relate their power plant and observe the engine state in the service to estimate loading factors.
Next, maximum powers and load factors used in the model are presented in Tables A1-A3:   Table A1. Maximum powers studied from observed GSE models and load factors in traffic. Loading factors are those components affecting the maximum power, obtained from the software Aviation Environmental Design Tool (AEDT 2d), and adjusted depending on the field study. Hereafter, values adopted for the model according to the traffic process in apron and also for luggage and goods loading and unloading, and passenger movement can be observed. The importance of such factors has increased due to the characteristic curves of power-RPM according to the mechanization of each vehicle for every operational status. Emission factor Emission factors are obtained from the same database than loading factors. In order to obtain them, tests are performed on engines with no operating hours, at maximum power and according to the type of fuel. Next, Table A4 shows the emission factors according to the pollutant gas for used GSE.

Damages factor
The damages factor enlarges the value of emissions produced by the years of use. Particularly, software EDMS V.5.0 proposes Equation (A1).
For each gaseous pollutant in the study, there are two types of auxiliary coefficients for damages measurement, according to the corresponding power range. Such coefficients were updated to 2010, the proposed year for the study. Therefore, it is presumed by default that the equipment has been used for nine years. The damages factor also depends on the year of service and the average lifetime of the GSE equipment. The following values are then deduced for all analyzed vehicles: for CO, coefficients A and b are 0.150 and 1 respectively, whereas for NOx, A and b are 0.008 and 1, respectively.
From the corresponding coefficients, together with the years of service and lifetime of each equipment, damages factors are obtained, as shown in Table A5.