Optimization of Aircraft Climb Trajectory considering Environmental Impact under RTA Constraints

In order to realize the concept of air traffic sustainable operation, taking the aircraft climbing stage as an example, firstly, we establish the vertical trajectory model of aircraft climbing, analyze the change rule of aircraft performance parameters under different indicated airspeed, and establish the RTA and RHA constraint models according to the waypoint constraints. ,en, considering the fuel economy and the greenhouse effect of pollutant emission, we establish a multiobjective model of aircraft flight parameter optimization, and, based on the multiobjective genetic algorithm, we establish an optimization model. Finally, we use B737-800 aircraft to carry out simulation experiments and find that, with the change of speed, fuel consumption and warming trend are different, and “objective weight, aircraft mass, flight distance, RTA time window, and wind” have different effects on the optimization results. ,e results show that this optimization method has a good compromise between fuel consumption and greenhouse effect by changing the weighting factor. By optimizing the flight parameters of the aircraft, it can effectively reduce the impact on the environment and provide theoretical support for the green flight of the aircraft.


Introduction
e climbing stage is an important stage in the flight of the aircraft, and the speed, flight time, fuel consumption, and climbing mode in the climbing flight performance parameters have important influence on the whole flight stage. erefore, the research on the climbing stage of aircraft can clarify the mechanism of influence on aircraft climbing performance and provide theoretical support for optimizing aircraft fuel efficiency and improving environmental friendliness. In terms of departure trajectory optimization, Ho-Huu et al. [1] introduced a two-stage optimization framework to deal with the optimization and selection of aircraft departure trajectories as well as the flight allocation among these routes, so as to minimize the accumulated noise and fuel consumption. Visser et al. [2] developed an airport abatement tool for aircraft take-off noise based on noise models, geographic information systems, and dynamic trajectory optimization algorithms and validated it at Amsterdam Airport Schiphol using B737-300 aircraft. Prats et al. [3] present a strategy for designing noise abatement procedures aimed at reducing the global annoyance perceived by the population living around the airports. A nonlinear multiobjective optimal control problem is implemented and numerically solved obtaining minimal annoyance trajectories. Torres et al. [4] comprehensively consider the impact of noise and NO x and CO 2 emissions on the environment, establish a nonlinear multiobjective optimization, and solve to get the optimal departure trajectory of the aircraft. Visser et al. [5] proposed a multiphase/ multicriteria trajectory optimization framework for citypair, so that the aircraft can fly to optimum flight paths with the lowest possible noise and emissions. In Ho-Huu et al.'s work [6], in order to mitigate aircrafts negative impacts (e.g., noise and pollutant emission) on near-airport communities, a multiobjective optimization model involving noise and fuel consumption was proposed, and a multiobjective evolutionary algorithm based on decomposition (MOEA/D) was developed to solve the problem and optimize the departure trajectory of the aircraft. en, they present a new multiobjective optimization formulation for the design and allocation of optimal aircraft departure routes. In addition to the two conventional objectives of noise and fuel consumption, a new objective of aircraft flight frequency is introduced and solved by MOEA/D optimization algorithm. Finally, the reliability and applicability of the model are verified by an example [7].
In terms of approach trajectory optimization, Sang et al. [8] adopted the idea of Continuous Descent Operation (CDO) and the strategy of adjusting Top of Descent (TOD) to optimize the vertical flight trajectory of the aircraft for the approach of descent phase of the aircraft and finally realized the goal of reducing fuel consumption of the aircraft. In terms of cruise trajectory optimization, Jensen et al. [9] quantitatively analyzed the influence of different flight levels on flight fuel and optimized the flight level to reduce the fuel consumption of aircraft. Garcia et al. [10,11] used the automatic control optimization algorithm to realize the optimization solution of aircraft fuel-saving trajectory. Four optimization methods were compared, and then two methods with better optimization performance were studied in depth; it was pointed out that the pseudospectral collocation methods could be used to solve the optimized trajectory quickly and effectively. In terms of meteorological factors, Patron et al. [12] studied and modeled the mechanism of the influence of wind on aircraft fuel consumption by combining horizontal and vertical tracks. In terms of combination optimization, Mendoza et al. [13,14] studied the realization method of the optimal trajectory of aircraft. First, cost index was introduced to optimize the vertical fuelsaving trajectory by finding the optimal combination of altitude, surface speed, and Mach number. en the dynamic weight graph is used to find the economic reference trajectories of the aircraft in horizontal and vertical directions. In terms of environmental protection, Rosenow et al. [15] analyzed the factors influencing the generation of contrails and studied the impact of contrails on global warming and then optimized the aircraft trajectory by reducing the generation of contrails. e above studies mainly optimize the flight trajectory of aircraft from the perspective of fuel economy, without considering the impact of airspace operation restrictions or the impact of engine exhaust on the environment. e concept of required time of arrival (RTA) in trajectory based operation (TBO) can well solve the airspace operation limitation. It plans, manages, and optimizes the flight in the whole operation process strategically based on the time management, the information exchange between the air and the ground system, and the ability of the aircraft to accurately plan the flight trajectory in time and space and constructs the collision-free flight path between the aircraft, so as to avoid potential flight conflicts in airspace [16,17]. Ramon and Adrian et al. [18,19] studied CDO in the RTA mode and believed that this mode was a potential solution to reduce fuel consumption, pollutant emission, and noise impact in the airport area without reducing the airport capacity, and it was beneficial to allocate RTA before the aircraft began to descend. Higuchi et al. [20] used a threeparameter model to simulate the flight operation of aircraft and optimized the aircraft in approach according to RTA constraints to improve the operation efficiency of busy airports. Vilardaga et al. [21] carried out quantitative tradeoffs in terms of fuel and time consumption to optimize the aircraft's 4D trajectory under the condition of meeting the constraints of RTA. Alejandro et al. [22] proposed a method that could provide an optimal combination of Mach Numbers for different waypoints to meet the time constraints required for aircraft during flight at a fixed altitude while taking account of aircraft fuel consumption. In the following year, Alejandro et al. [23] proposed an algorithm that satisfies the time of arrival constraint and has an economical vertical reference trajectory. Relative to the flight cost of the reference trajectory, the RTA constraint is satisfied and the flight cost of the optimized trajectory is reduced. However, the study did not take into account the environmental impact of emissions from aircraft engines, which deviates from NextGen's environmental goals. e impact of aircraft on the environment is mainly the greenhouse effect caused by greenhouse gases emitted by engines, which ultimately affects global temperature changes. How to achieve fuel-efficient flight and reduce the impact of emissions on the environment while meeting the requirements of waypoint coordinated operation is an important part of realizing the new generation of air traffic and transportation system and ICAO Aviation System Block Upgrades (ASBU), as well as a problem that the civil aviation industry is concerned about and has to solve. In order to solve the deficiencies in the above studies, this paper considers the requirements of airspace operation restrictions, sets up RTA and required height of arrival (RHA) constraints on waypoints of aircraft in the climbing phase, and considers the temperature changes caused by aircraft emissions of greenhouse gases. e multiobjective optimization model of aircraft climbing stage under single waypoint constraint is established, and the genetic algorithm is designed to optimize aircraft flight parameters, and the sensitivity analysis of relevant influencing factors is carried out, hoping to provide theoretical support for the concept of green development of air transport.

Problem Description.
Route climb generally refers to the climb process in which an aircraft starts from an altitude of 1500 ft above the airport surface and climbs in a certain way and increases to the specified cruising altitude and cruising speed [24]. During the flight, the aircraft shall meet the requirements of waypoint collaborative constraint issued by the air traffic management system, coordinate the RTA constraint window and RHA constraint window of the aircraft at relevant positions, and construct the conflict-free flight path between the aircraft to avoid potential flight conflicts in the airspace. erefore, the aircraft needs to plan the flight distance, flight altitude, and flight time in advance. is section first analyzes the departure vertical climb process of the aircraft, including the end of climb and the beginning of cruise level flight. Assuming that the aircraft climbs according to the same indicated airspeed, the departure physical model of the aircraft is established, and the climb parameters and cruise parameters of the aircraft are calculated.
en, the flight parameters of the aircraft are optimized considering the fuel consumption of the aircraft and the impact of engine exhaust on the environment. Finally, the mathematical model of the optimization objective is described in detail, and the research framework of the problem is shown in Figure 1.

Vertical Section of Aircraft Climb.
Due to the greater flight maneuverability of the aircraft below 10,000 feet, the aircraft is subject to more intervention by air traffic control (ATC), and the flight is less predictable. erefore, this paper only considers the climbing phase of aircraft above 10,000 feet and the beginning of route cruise. e departure flight under the constraint of single waypoint starts from point A to the top of climb (TOC), and the aircraft climbs at the same indicated airspeed and then at the same speed to point B (with the same altitude). In the whole departure process, only point B has RTA and RHA constraints. As shown in Figure 2, the single waypoint constrained departure model of the aircraft, and R is the total flight range of the aircraft, which is composed of climb range r c and cruise range r cr .
Here, we introduce the idea of microelement to divide the departure trajectory into n s segments. e flight parameters of each segment can be determined by the flight parameters of the intermediate trajectory point of the microelement, determined by the horizontal distance and flight height of the aircraft from point A. e horizontal flight distance between the aircraft and point A in the i-th (i � 1, 2, ..., n s ) leg can be expressed as where i represents the serial number of the current segment; Δr and Δh represent the horizontal flight distance and climb height of the aircraft in a microsegment, respectively (unit: m); V T,i is the true airspeed (unit: kt); V W is the effective wind speed (unit: kt); c c,i is the climb rate (unit: m/s). e flight altitude of the aircraft at the middle point of the i-th leg is e flight time for the aircraft to reach the middle point of the i-th leg is where V W,mps,i is the wind speed (unit: m/s). e aircraft uses the maximum continuous thrust to climb a certain altitude. When the flight speed is fixed, the climb range r c and cruise range r cr can be calculated. Suppose that the flight distance from TOC to point B is short; the fuel consumption of the aircraft in this process can be negligible compared with the total mass of the aircraft, and then the total masses of the aircraft can be approximately equal when the aircraft flies from TOC to point B. Note that i cr is the last trajectory point in the cruise phase, and r cr and r c can be calculated as follows: Δr and Δh are calculated as follows: where Δr 0 and Δh 0 are used to represent the preset horizontal flight distance step and climb height step of the aircraft.
In this way, the flight distance, altitude, and time state parameters of the aircraft in the process of climbing under the constraint of single waypoint can be described by the following formulas: Single waypoint constraints include RHA constraints and RTA constraints. e RHA constraint of point B in Figure 2 is represented by the end height H c ; the RTA constraint is presented in the form of time window, which can be represented as where t R, max and t R, min represent the maximum and minimum values of RTA window, respectively.

Aircraft Climb Parameters.
Because the aircraft obeys the law of conservation of energy in the climbing process, we can use the full energy model to model and analyze the climbing process of the aircraft. e work of the external force on the aircraft is equal to the change of the mechanical energy of the aircraft. e forces on the longitudinal axis of the aircraft during the climb mainly include engine thrust T c and drag D. e calculation formula is as follows: e flight path of the aircraft is mainly changed by adjusting the position of the engine thrust throttle and elevator without considering the lifting device. According to (8), it is not difficult to find that any two of the three variables, thrust T c , velocity V T , and altitude change rate dh/dt, can determine another variable. e current research is generally to calculate the rate of altitude change with given speed and thrust. is is also the case in the climbing section. At this time, fix the throttle position of the aircraft and maintain the flight speed at the same indicated airspeed V I or Mach number M as shown in the following equation: By separating dh/dt to the left of the equation, we can get It is a function of Mach number M, which represents the ratio of the residual thrust used for climb to the residual thrust used for acceleration when climbing at a given speed. e specific calculation formula can be expressed as follows: In the process of aircraft climb, for certain climb height and flight distance, different engine thrust T c , flight speed V T , and climb rate c c will lead to different flight time and fuel consumption. For a given climbing altitude ΔH(m), the flight time t c (s) and fuel consumption F c (kg) of the climbing process can be obtained according to the principles of kinematics and dynamics: where F f,c is the climbing fuel flow (kg/s) and the climbing rate c c is calculated as follows: e climbing fuel flow F f,c (kg/s) at maximum thrust is calculated as follows:  where C f,1 and C f,2 are constants related to aircraft type and, for the maximum climb thrust T max climb of aircraft, the reader is referred to [19].

Aircraft Cruise Parameters.
Without considering the turning of the aircraft, the aircraft can be regarded as a process of constant altitude in cruise flight. At this time, the aircraft can be regarded as a point-mass model. e gravity of the aircraft in the vertical axis direction is equal to the lift, and the engine thrust in the longitudinal axis direction is equal to the resistance of the aircraft. We have the following dynamic equations: where m A is the mass of the aircraft (unit: kg); the acceleration of gravity g � 9.8(m/s 2 ); V T,mps is the true speed (unit: m/s); S is the surface area of the wing (unit: m 2 ); C L is the lift coefficient; T cr is the engine thrust (unit: N); D is the resistance (unit: N); C D is the resistance coefficient. C D is a function of C L , and the relationship between them is as follows: where C D0,CR and C D2,CR are parameters related to aircraft type (see [25]). Assuming that the aircraft flies at the same speed V I (kt), V T can be calculated by V I : where δ � P H /P 0 is the pressure ratio; sound speed a 0 � 340.294 m/s. In the process of cruise flight, for the given flight distance segment ΔR(m), according to the kinematics principle, the cruise flight time t cr (s) and fuel consumption F cr (kg) can be calculated. e specific calculation is as follows: where V W,mps is the effective wind speed (unit: m/s); F f,cr is the cruise fuel flow (unit: kg/s):

Establishment of Multiobjective Function.
According to the research problem, we establish a multiobjective optimization model. e optimization objectives of aircraft climb parameters under single-point constraints include fuel consumption and temperature rise. e optimization model can be expressed as follows: where λ 1 and λ 2 represent the impact of fuel consumption objective and temperature rise objective on the optimization results.
(1) Fuel consumption calculation When the flight meets the constraints of RHA and RTA, the fuel consumption F s in the process of aircraft departure can be expressed as follows: where F s,c is the fuel consumption during climbing; F s,cr is the fuel consumption during cruise. (2) e impact of aircraft flight on climate change e development of air transport brings about huge economic benefits to the region and also impacts the global environment. According to the report on air activities issued by IPCC, there are five main substances that affect the global atmosphere: CO 2 , NO x , contrail, water vapor, and cloud. e impact process is shown in Figure 3. e three emissions that have the greatest impact on global temperature change are CO 2 , NO x , and contrail. Because the formation conditions of contrail are very strict, its impact cycle on temperature change is short, and the impact area is small [26], this paper focuses on the effects of CO 2 and NO x emissions on the atmosphere, regardless of the role of contrail.
However, because different gases have different radiation properties and the greenhouse effect is not the same after gas emissions, how to evaluate the impact of engine emissions on climate change is particularly important. At present, the Journal of Advanced Transportation 5 main methods to study the impact of engine emissions on climate change are as follows: global warming potential (GWP), global temperature change potential (GTP), radiation intensity, and surface temperature change [27]. However, GWP only represents the integral effect of greenhouse gas radiative forcing, which cannot reflect the impact of short-lived gases on the climate. Compared with GWP, GTP can directly give the ratio of surface temperature change, which has a better evaluation effect. erefore, it is widely used to characterize the impact of greenhouse gas emissions from aircraft engines on the environment [28]. GTP can be divided into two types: one is absolute pulse global temperature change potential (APGTP); the other is absolute sustained global temperature change potential (ASGTP). Because ASGTP represents the temperature effect of greenhouse gases under the condition of continuous emission, this paper studies the instantaneous emission of greenhouse gases during the climb process of aircraft, so APGTP is selected to characterize the impact of gas emissions on the temperature rise. T AP,x (n/kg) is used for APGTP, and the calculation formula is as follows: where x is the type of emission gas (CO 2 and NO x ); A x is the radiation forcing caused by the emission of 1 kg of greenhouse gas (unit: W · m − 2 · kg − 1 ); T h is the timescale of the influence of exhaust gas on temperature (unit: year); α x is gas life cycle (unit: year); τ is the timescale of climate change response (unit: year); C is the environmental heat capacity of the climate system (unit: year · W · m − 2 · n − 1 ); the absolute values of the three "emissions per kilogram" indicators in different time ranges are given by comparing the climate impacts from well-mixed greenhouse gases and heterogeneous forcing [29,30]. is paper studies the impact of air emissions on the environment over a period of 20 years. e impact of greenhouse gas emissions on the environment can be expressed by the temperature rise, and the calculation formula is as follows: where E x is the amount of emission. E I,CO 2 and E I,NO x represent the reference emission indexes of CO 2 and NO x , respectively, where E I,NO x is obtained by linear interpolation of the emission index provided by the engine manufacturer and the corrected fuel flow in a log-log coordinate system; the actual fuel flow correction method is Boeing Method 2 (BM2). Please refer to [31] for detailed steps. e reference emission index of different emissions is obtained according to the ICAO engine exhaust emission database [32]. Figure 2 is the convergence point in the process of aircraft departure. In order to avoid potential flight conflicts in the airspace, the air traffic management system regulates the trajectory of the aircraft by adjusting the flight speed of the aircraft. At this time, the constraint of single waypoint can be expressed as h s,n s � H cr ,

Constraints. Suppose that point B in
where t R, min t R, max is the RTA constraint time window. When the aircraft flies according to the given RTA window, the flight speed matrix V I,L under the limitation of time window can be obtained, which is expressed as follows: e flight speed under RTA restriction shall be as follows:

Fitness Function.
According to the contents in Section 3.1, the optimization objectives of aircraft flight parameters include (1) minimum fuel consumption of aircraft and (2) minimum temperature rise. e final optimization objective is to make both objective functions as small as possible. erefore, the fitness f s can be set as follows:  where w 1 and w 2 are the weights of the two optimization objectives, meeting w 1 + w 2 � 1; F s and T s,E refer to the individual fuel consumption and temperature changes in the iterative optimization process; F s,r and T s,E,r are the fuel consumption and temperature rise under reference conditions, respectively.

Climbing Trajectory Optimization
Process. Genetic algorithm (GA) is a computational model to simulate the natural selection and genetic mechanism of Darwinian biological evolution. It is a method to search for the optimal solution by simulating the natural evolution. It was proposed by Professor John H. Holland of the University of Michigan in the United States [33]. GA is used to optimize the climb speed V I of the aircraft, and binary encoding method with 8bit encoding length is selected for encoding. Only a gene V I is included on the chromosome, and the chromosome model is shown in Figure 4. Set the population size as N GA � 10; use roulette method to select parents; the selection probability is 0.5; carry out single-point crossing and mutation operation on the parents after the selection operation, the crossover probability is 0.9, and the mutation probability is 0.1. e individuals with the highest fitness value in the current generation are transferred to the next generation with a transfer probability of 1, which ensures that the optimization algorithm converges quickly and the parents will not get worse after cross mutation. Each generation of the optimization process is recorded as i. When the termination generation L E is reached, the optimization ends, and the optimal individual, objective function, and other relevant parameters are retained. e optimization model of genetic algorithm is shown in Figure 5.

Simulation Condition Setting.
is simulation case uses B737-800 aircraft for optimization and uses the optimization algorithm described above to optimize the climb trajectory of the aircraft. Assume that the aircraft is flying under standard atmospheric conditions, with temperature deviation ΔT � 0°C, calm wind, flight distance R � 250 km, initial mass of 65000 kg, initial climb height of 10000 ft, and TOC of 8400 m. e horizontal flight distance and climb altitude in a microsegment are Δr 0 � 1000 m, Δh 0 � 1000 ft, respectively, and the change range of the speed is 200-310 kt. According to this range, the flight time range

Simulation Results.
First of all, in order to get the influence of flight speed on the optimization target, this paper tests the flight conditions at different speed and calculates the flight parameters and trajectories. e parameters corresponding to the lowest fuel consumption F s and the corresponding temperature rise T s,E are shown in Table 1.
e influence of flight speed on flight trajectory is shown in Figure 6, and the step V I is 20 kt. It can be seen from the figure that the larger V I is, the smoother the climbing slope is, and the shorter the level flight distance is.
When the target fuel consumption is the lowest, that is, w 1 � 1, w 2 � 0, the optimization result without RTA time constraint is shown in Figure 7. When the target temperature consumption is the lowest, that is, w 1 � 1, w 2 � 0, the optimization result without RTA time constraint is shown in Figure 8.

Sensitivity Analysis of Influencing Factors.
Due to the volatility of the optimization results of genetic algorithm, in the follow-up study of the influencing factors, each case is carried out five times to eliminate the accidental error and find out the general influence rules of different influencing factors on the optimization results.

5.3.1.
e Influence of Objective Weight on Optimization Results. In order to find the influence of multiobjective weight coefficient on the optimization results, the change of flight parameters is analyzed by increasing w 1 from 0 to 1 according to Step 0.1. According to Figure 9, with the increase of w 1 , the optimal flight speed of the aircraft decreases from 310 kt to 251 kt, and the fitness also shows a decreasing trend. In phase w 1 ∈ [0, 0.6], the decreasing trend of adaptability is more significant than that of phase w 1 ∈ [0.6, 1], and the change trends of fuel consumption and speed of the aircraft are basically the same, but the trend of temperature rise is opposite to that of the other three parameters. is is because the trends of fuel and temperature rise with speed V I are different. When w 1 increases from 0 to 1, the corresponding fuel consumption will be reduced by about 6%, CO 2 emissions will be reduced by about 6%, NO x will be reduced by about 17%, and the temperature rise will be increased by about 4%.

Influence of Aircraft Mass on Optimization Results.
In order to study the influence of aircraft mass on the optimization results, the initial mass of the aircraft is changed from 60000 kg to 70000 kg, and the change step is 1000 kg to analyze the change of each flight parameter. According to Figure 10, with the increase of aircraft mass, the optimal flight speed of the aircraft increases from around 253 kt to around 280 kt, and the fitness also shows a decreasing trend. e fuel consumption of aircraft is basically consistent with the change trend of temperature rise. is is because, with the increase of mass, the aircraft needs more power to reach the same speed, so as to consume more fuel. At the same time, the emissions of CO 2 and NO x also increase correspondingly, resulting in the increase of temperature rise. When the initial mass of the aircraft changes from 60000 kg to 70000 kg, the corresponding fuel

Journal of Advanced Transportation
consumption increases by about 12%, CO 2 emissions increase by about 12%, and NO x increases by about 26%, resulting in an increase of temperature rise by about 5%.

Influence of Flight Distance on Optimization Results.
In order to study the influence of the flight distance on the optimization results, the horizontal flight distance of the aircraft is changed from 200 km to 300 km, and the change step is 10 km to analyze the change of each flight parameter. According to Figure 11, with the increase of flight distance, the optimal flight speed of the aircraft decreases from 274 kt to 265 kt, and the fitness also shows a decreasing trend. e fuel consumption of aircraft is basically consistent with the change trend of temperature rise. is is because, with the increase of flight distance, the aircraft consumes more fuel,    of RTA time, temperature rise shows an upward trend first and then has a downward trend. is is mainly because the increase of t R, max makes the minimum value of the flight speed satisfy the RTA window increase, which further affects the fuel consumption of the aircraft and the change of the temperature rise, explaining that the change trends of the fuel consumption of the aircraft and the global total temperature are different with the change of the flight speed. With the increase of t R, max from 710s to 890 s, the corresponding fuel consumption will be reduced by about 3%, CO 2 emissions will be reduced by about 3%, NO x will be reduced by about 11%, and the temperature rise will be increased by about 1%.
Next, according to the value of t R, min , the impact of optimization results is analyzed. As shown in Figure 13, the minimum value of RTA window has little effect on the optimized flight parameters in interval t R, min � [700, 850]s, but there is a significant reduction in t R, min � 870 s. is is because when the RTA time window is satisfied, the larger the minimum value of the RTA window is, the smaller the corresponding maximum flight speed is. When t R, min � 870 s, the corresponding maximum flight speed is about 263 kt. In conclusion, the influence of t R, min on each index is not obvious.

e Influence of Wind on Optimization
Results. In order to analyze the influence of wind on the optimization results, the wind speed is changed from -30 kt (upwind) to 30 kt (downwind) during the optimization test, and the change step is 5 kt. Set RTA window: [t R, min , t R, max ] � [690, 890]s, and set w 1 � w 2 � 0.5. According to Figure 14, with the increase of V w , the fitness of optimization results shows an increasing trend, while the flight speed shows a decreasing trend. is is because the optimal speed is relatively fixed, and the change of wind speed will lead to the change of true airspeed. In the same way, the fuel consumption of aircraft and the temperature rise show a decreasing trend first and then have an increasing trend.
is is mainly because the downwind or headwind will lead to the change of true airspeed under the same other conditions, which will affect the fuel consumption of the aircraft and the emissions of CO 2 and NO x . When the wind speed changes from -30 kt to 30 kt, the corresponding fuel consumption is reduced by about 4%, CO 2 emissions are reduced by about 4%, NO x is reduced by about 6%, and the temperature rise is reduced by about 2%.

Warming Index of Each Influencing Factor.
Considering the fuel consumption and the corresponding temperature rise, the warming index is defined as the average temperature rise produced by consuming unit fuel. As shown in Figure 15, the warming index corresponding to "objective weight, aircraft mass, horizontal flight distance,      wind speed, and RTA time" is described. With the increase of w 1 from 0 to 1, the warming index increases by about 10%; with the change of aircraft initial mass from 60000 kg to 70000 kg, the warming index decreases by about 7%; with the change of aircraft horizontal flight distance from 200 km to 300 km, the warming index increases by about 6%; with the change of wind speed from -30 kt to 30 kt, the warming index increases by about 2%; with the change of t R, min , the warming index does not change obviously; with the change of t R, max from 710 s at the beginning to 890 s, the warming index increases by 5%.

Conclusion
In this paper, the trajectory optimization problem of aircraft in departure climb phase is studied when the waypoint is constrained. At the same time, the temperature rise caused by the minimum fuel consumption and the emissions of CO 2 and NO x is considered.
Firstly, the climbing trajectory model of the aircraft is established to analyze the flight parameter variation law of the aircraft at different speed; then the multiobjective optimization model is established for the optimization objective and the genetic algorithm is designed to solve it. Finally, the simulation test is carried out according to the optimized parameters. According to the research results, the following conclusions can be obtained: (1) In the process of aircraft climb, with the change of speed, the trend of temperature rise caused by fuel consumption and engine exhaust is different, so the temperature rise tradeoff between fuel consumption and exhaust can be realized by adjusting the size of w 1 and w 2 . When w 1 increases from 0 to 1, the corresponding fuel consumption will be reduced by about 6%, CO 2 emissions will be reduced by about 6%, NO x will be reduced by about 17%, and the temperature rise will be increased by about 4%. (2) With the increase of aircraft mass, the fuel consumption and the temperature rise brought by engine exhaust are increasing. When the initial mass of the aircraft changes from 60000 kg to 70000 kg, the corresponding fuel consumption increases by about 12%, CO 2 emissions increase by about 12%, and NO x increases by about 26%, resulting in an increase of temperature rise by about 5%.
(3) When the horizontal flight distance of the aircraft changes from 200 km to 300 km, the corresponding fuel consumption increases by about 35%, CO 2 emission increases by about 35%, and NO x increases by about 24%, resulting in an increase of temperature rise by about 43%. (4) e increase of the maximum time window of RTA will lead to the decrease of fuel consumption, while the temperature rise brought by exhaust gas shows a trend of first decreasing and then increasing. e minimum value of RTA window has little effect on the optimized flight parameters.
With the increase of t R, max from 710s to 890s, the corresponding fuel consumption will be reduced by about 3%, CO 2 emissions will be reduced by about 3%, NO x will be reduced by about 11%, and the temperature rise will be increased by about 1%. (5) Downwind will reduce the temperature rise caused by fuel consumption and engine exhaust, and upwind will bring adverse effects. When the wind speed changes from -30 kt to 30 kt, the corresponding fuel consumption is reduced by about 4%, CO 2 emissions are reduced by about 4%, NO x is reduced by about 6%, and the temperature rise caused by this is reduced by about 2%. (6) e change rule of temperature rising index corresponding to different influencing factors is not the same. e environmental friendliness of the aircraft can be evaluated according to the warming index.
According to the research content and results of this paper, it can provide theoretical support for the realization of aircraft safety and green and coordinated operation in the future. e next step can be further studied from two aspects. First, the aircraft can be analyzed and studied when it is subject to multiple waypoint constraints. Secondly, the algorithm with higher precision is used to solve the model, which can improve the efficiency and precision of the method.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request. Most of the data of this study are included within the article.  Figure 15: Warming index of each factor.

Conflicts of Interest
e authors declare that they have no conflicts of interest.