Truck Scheduling Problem Considering Carbon Emissions under Truck Appointment System

: Aiming at the truck scheduling problem between the outer yard and multi-terminals, the appointment optimization model of truck is established. In this model, the queue time and the operation time of truck during the appointment period of di ﬀ erent terminals are di ﬀ erent. Under the restriction of given appointment quotas of each appointment period, determine the arrival amount of trucks in each appointment period. The goal is to reduce carbon emissions and total costs, improve the e ﬃ ciency of truck scheduling. To solve this model, hybrid genetic algorithm with variable neighborhood search was designed. Firstly, generate chromosomes, and the front part of the chromosome represents the demand for 40 ft containers and the back part represents the demand for 20 ft containers. Then, the route is generated according to the time constraint and appointment quotas of each appointment period. Finally, the neighborhood search strategy is adopted to improve the solution quality. The validity of the model and algorithm were veriﬁed by an example. A low-carbon scheduling scheme was obtained under truck appointment system. The results show that the scheduling scheme under truck appointment system uses fewer trucks, improves the e ﬃ ciency of delivery, reduces the total costs, and it takes into account the requirements of low carbon.


Introduction
In recent years, the rapid growth of container throughput has put forward higher requirements on the efficiency and service capacity of container terminals. On the one hand, the amount and intensity of operation in the terminals have increased. On the other hand, many trucks arrive at the terminal in the peak period causing terminal gate congestion. Terminal gate congestion leads to longer queue time, which reduces truck turnover efficiency, thereby trucks perform fewer tasks in a day. It is important to note that the truck is idling during the queueing at the terminal gate, they will consume fuel, resulting in fuel costs. The way for micro and small enterprises to reduce total costs is correct to supply chain management [1], while the main cost structure of larger enterprises is fuel cost [2]. Reducing fuel costs can effectively reduce total costs of larger enterprises. If idling time is too long, trucks will emit a greater amount of emissions compared to when they are moving. As is known to all, these emissions from trucks increase the risk of heart disease, asthma attacks, strokes and untimely deaths [3], and massive carbon emissions contribute to global warming. In recent years, the hot issue of greenhouse gas emissions reduction has attracted close attention from all over the world [4].
At present, there are two approaches to alleviate terminal gate congestion, one approach is to increase infrastructures, such as the number of gates and the area of the terminal. Another approach is to manage truck arrival and adopt a truck appointment system (TAS). In the TAS, the terminal operator announces the opening time and the maximum number of containers that can be accepted (appointment quotas) in each appointment period through the network information system. The truck driver can choose the arrival time according to his preference. TAS is an effective approach to alleviate the terminal gate congestion. TAS was first proposed by the local government of California [5]. TAS saves a lot of time and reduces air pollution and alleviates terminal gate congestion by spreading the demands of terminals over different appointment periods of the day. Since 1999, in practical application, TAS has achieved successful experience in Vancouver port [6] and Tianjin port [7].
Truck scheduling has always been a key link in the process of container terminal operation. The terminal divides the day into multiple appointment periods and sets the appointment quotas for each appointment period. The appointment quotas are determined by the yard crane availability. They are also set to avoid potential conflicts with other operations in a certain yard block or zones, such as vessel operations, warehouse operations, rail operations, and customs inspections [8]. Some scholars have studied truck scheduling under TAS. Chen G, Chen X et al. [9,10] established the appointment optimization model of trucks with the goal of simultaneously minimizing the waiting time of trucks and adjusting the number of trucks in each appointment period before and after the appointment. However, most researches under TAS is that truck is scheduled between single outer yard and single terminal at present. Assuming that terminal only has the demand for one type of container. However, in reality, the outer yard has a strong service capacity and generally performs the delivery service for multiple terminals at the same time. Terminal has the demand for both 40 ft containers and 20 ft containers. Thus, an effective truck scheduling scheme under TAS must service not only a single terminal, but also multiple terminals, and trucks can deliver different types of containers.
Based on the above analysis, this paper studies the truck scheduling problem about a single outer yard serving multiple terminals from the perspective of carbon emissions and TAS. The mixed-integer programming formulation of truck scheduling is proposed. A hybrid genetic algorithm with a variable neighborhood search (HGAVNS) is designed. Finally, an example is given to show that the HGAVNS can effectively reduce the total costs of scheduling, thereby reduce the carbon emissions of trucks and achieve the goal of truck scheduling between a single outer yard and multiple terminals.
The remainder of the paper is organized as follows: Section 2 reviews the related literature. Section 3 presents a description of the problem and proposed a mixed-integer programming formulation. Section 4 gives the details of the hybrid genetic algorithm with a variable neighborhood search. Computational results on the example of a port in China are reported in Section 5. Conclusions and future research directions are suggested in Section 6.

Literature Review
As mentioned previously, a number of scholars have studied the truck scheduling under TAS. Huynh [5] took into account both the turnover time of trucks and the utilization rate of yard crane and calculated the appointment quotas of each appointment period with the method of planning model and simulation joint optimization. On this basis, Huynh et al. [11] also considered how to determine the upper limit of the number of containers that can be accommodated in each container area in each time window in the case of missing the appointment. Zhang et al. [12] constructed a second-level queueing network of gate and yard and established appointment optimization model of trucks in order to minimize the average turnover time of trucks. In order to further improve the flexibility of TAS, Chen G et al. [13] relaxed the assumption that the expected arrival mode of trucks is known and proposed a dynamic TAS. Torkjazi et al. [8] designed a TAS to minimize the impact on the terminal and the operation of trucks, which not only alleviated terminal gate congestion but also made the truck deviate from its original arrival time as little as possible. Zehendner et al. [14] proposed a mixed-integer linear programming model of appointment quotas and cross-truck scheduling optimization, which not only could improve the service quality of trucks, but also improve the service quality of trains and ships. Obviously, delivery under TAS effectively reduces the queue time of trucks in terminal gate and the turnover efficiency of truck is greatly improved. However, appointment periods and appointment quotas also add new constraints to the operation and management of truck company. Namboothiri et al. [15] constructed a truck scheduling optimization model of a single truck company considering the constraint of appointment periods to reduce the transportation costs of truck fleet. Islam and Olsen [16] emphasized that the reducing of the truck's flexibility and the possibility of sharing transportation capacity between truck companies were major challenges faced by the implementation of the TAS. Therefore, the operation strategy directly affected the benefits generated by the implementation of TAS. Schulte et al. [17] established a cooperative scheduling model of multiple trucks under the TAS, which reduced the travel distance of trucks when truck is empty, thus reduced the costs and carbon emissions. Phan et al. [18] studied the trucks' appointment optimization problem under the joint decision of terminal and truck company by establishing a mathematical model that included two subproblems, namely, prediction of truck turnover time and truck scheduling optimization plan. Chen et al. [19] proposed a "ship dependent time window" method to control the arrival of the truck and alleviate the terminal gate congestion. Traditional genetic algorithms, multi-social genetic algorithms, hybrid genetic algorithms, and simulated annealing algorithms were used to solve the problem. Niu et al. [20] studied a problem with truck scheduling and container area allocation. He proposed a particle swarm optimization algorithm and colony optimization algorithm based on swarm intelligence technology to solve the problem. Zhang et al. [21] proposed a new truck system optimization method to reduce the waiting time of trucks at the terminal gate and the operation time of trucks in the terminal, and estimated the waiting time of trucks with non-stationary queuing theory and verified the validity of the algorithm and model with numerical experiments. However, these studies are the truck scheduling between single outer yard and single terminal. Terminal only has a demand for one type of container. They haven't studied the problem of trucks service multiply terminals and each terminal has a demand for two types of containers under TAS.
For research on calculating carbon emissions has been addressed by Sim [22], he proposed a system dynamic method to evaluate the total carbon emissions of trucks in container terminals, including loading and unloading containers, container transportation, container terminal operation, and calculated the carbon emissions reduction of trucks from 2017 to 2030. Heilig et al. [23] put forward a multi-objective delivery problem considering carbon emissions, using the archive of simulated annealing method and a kind of visualization technology to solve. On this basis, Heilig et al. [24] deepen the earlier study by building a multi-objective model. The goal of this model is to minimize fixed costs, variable transportation costs, and carbon emissions. Yu et al. [25] proposed a method to evaluate the carbon emissions of trucks in the loading process. Li et al. [26] used the total waiting time and the carbon emissions when trucks are idling to balance the service quality and the green performance of the whole system. He et al. [27] proposed a mixed-integer programming model on the truck scheduling problem. The goal is to minimize all trucks delay deviation and the total transportation energy consumption of trucks. A synthetic optimization method based on simulation is proposed to solve the problem. The method is evaluated by simulation and spatial optimization algorithm. Wang et al. [28] predicted the CO 2 emissions using the Motor Vehicle Emission Simulator (MOVES). The results show that good road conditions can reduce CO 2 emissions. Tolliver et al. [29] used two methods to estimate fuel consumption both from weight and distances based on different truck types. It pointed out that different truck configurations have a great impact on fuel consumption and route selection.
At present, most of the researches on truck scheduling are focused on the inner trucks and few studies are focused on the outer trucks. Compared with the inner trucks, the outer trucks serve multiple terminals and have a longer transportation distance, so the truck scheduling scheme has a greater impact on the waiting time. As the number of outer trucks is also limited, its waiting time will decrease the overall turnover efficiency of trucks. In view of these problems, we study the problem of trucks delivering containers from outer yard to terminals. The truck appointment optimization model is established from the perspective of carbon emissions and TAS. The hybrid genetic algorithm with variable neighborhood search (HGAVNS) was designed based on the advantages of genetic algorithm (GA) and variable neighborhood search (VNS) to solve the model. The goal of this model is to minimize the total costs including fixed costs, opportunity loss costs, carbon emissions costs, and fuel costs.

Problem Description
This paper considers truck scheduling between an outer yard and multiple terminals under TAS form the perspective of carbon emissions. There are a number of trucks can be used to deliver containers. Trucks must meet the demand of terminals within the time window.
Under the TAS, an effective link has been established between each terminal and the outer yard. The outer yard can predict (generally 24 h in advance) the demand of each terminal. The demand of each terminal is determined by the number of vessels need to load. The terminal divides the day into multiple appointment periods and gives the appointment quotas of each appointment period. The appointment quotas are different in different appointment periods of each terminal. The waiting time at the terminal gate and the operation time inside terminal are also different. The terminal needs to ensure a certain length of each appointment period so that trucks can reach the terminal during appointment period. The outer yard selects appropriate appointment periods according to its own transportation capacity and gives the plan of the number of different types of containers will be carried by trucks during this period. Terminals determine the final plan and notice outer yard. The outer yard arranges the transportation capacity according to the specific plan of each terminal. The number of containers transported to the terminal must meet the demand of each terminal. The scheduling scheme of this paper is to use limited truck resources to arrange a certain number of trucks delivering containers to different terminals in different appointment periods. A truck can carry one 40 ft container or two 20 ft containers at the same time. While meeting the limitation of appointment quotas in each appointment period and the demand of each terminal, the scheduling scheme should save energy consumption and improve the benefit and efficiency of the outer yard. Figure 1 shows the operation flow of trucks from a single outer yard to multiple container terminals. Truck 1 and truck 3 both deliver containers going to one terminal respectively. Trucks enter the terminal gate and unload containers. After finishing the scheduling task, trucks return to the outer yard and wait until the next scheduling task. In reality, if the demand for 20 ft containers is odd, this can happen a truck delivers two 20 ft containers going to two terminals. For example, truck 2 delivers two 20 ft containers going to terminal 1 and terminal 2 respectively. Firstly, truck 2 arrives at terminal 1 and unloads the container. Then, truck 2 goes to terminal 2. If the arrival time is within the appointment period of terminal 2, truck 2 enter the terminal 2 and unload the container immediately. Otherwise, truck 2 need to wait outside the terminal 2 until the appointment period. After truck 2 unload the container, truck 2 return to the outer yard and wait until the next scheduling task. The trucks are idling during queuing at the terminal gate and operating inside terminal. At this point, truck has energy consumption. When truck waits outside the terminal, there will be opportunity loss costs.
Outside yard

Assumptions
This study is based on the following assumptions: (1) Under the TAS, each terminal divides each day into 12 appointment periods (the length of each appointment period is 2h ). The period division of each terminal is consistent. If each appointment period is too long, the queue time of the truck at the terminal gate and the operation time inside the terminal are inaccurate. If each appointment period is too short, the queue time of the truck at the terminal gate and the operation time inside the terminal are accurate, but the flexibility of the truck scheduling is too poor. This paper assumes each appointment period is

Assumptions
This study is based on the following assumptions: (1) Under the TAS, each terminal divides each day into 12 appointment periods (the length of each appointment period is 2 h). The period division of each terminal is consistent. If each appointment period is too long, the queue time of the truck at the terminal gate and the operation time inside the terminal are inaccurate. If each appointment period is too short, the queue time of the truck at the terminal gate and the operation time inside the terminal are accurate, but the flexibility of the truck scheduling is too poor. This paper assumes each appointment period is 2 h, which ensures the accuracy of queue time and operation time and the flexibility of the truck scheduling. (3) A terminal has the demand for both 20 ft containers and 40 ft containers and the number is constant. Since non-standard containers are rarely used in reality and there are no trucks of the same size as the non-standard containers. This paper assumes that each terminal only has the demand for 20 ft containers and 40 ft containers. (4) After unloading the container, the truck should return to the outer yard for the next scheduling.
There is no energy consumption when the truck is waiting outside the terminal, but there are opportunity loss costs (in the comparison experiment, this can happen that truck is waiting outside the terminal). There are no costs when the truck is waiting at the outer yard for the next scheduling. There is energy consumption when the truck queues at the terminal gate, operates inside the terminal and loads containers at the outer yard. (5) This paper assumes only has one type of truck. Driving speed, the coefficient values of CO 2 emissions and other performance of each truck are the same. All trucks use diesel and the coefficient values of CO 2 emissions are only related to the type of fuel. (6) Trucks can make multiple trips between the outer yard and terminals to complete the scheduling task. A truck can carry one 40 ft container or two 20 ft containers. (7) This paper doesn't consider other interference factors on the truck. Road conditions, weather, and other factors will affect the speed of the truck. This paper doesn't consider these factors and assumes the speed of the truck is fixed. (8) Under the TAS, the terminal will allocate the time window and appointment quotas for the outer yard. Containers must be delivered to the terminal within the time window of it. If containers aren't delivered to the terminal within the time window of it, it will affect the vessel's delivery plan and cause great costs.

Parameters and Variables
This paper sets the following variables: (1) Input variables: pi,k : The number of 20 ft container delivered by truck k to the terminal i during the appointment period p, T pin,k : The moment of truck k arriving at terminal i during the appointment period p, this is the n th delivery of truck k, T d pn,k : The moment of truck k departing outer yard during the appointment period p, this is the n th delivery of truck k, t w in,k : The waiting time of truck k outside the terminal i, this is the n th delivery of truck k, T m in,k : The moment of truck k departing terminal i, this is the n th delivery of truck k, T r pn,k : The moment of truck k returning to the outer yard during the appointment period p, this is the n th delivery of truck k, t y pn,k : The waiting time of truck k at the outer yard before departing outer yard during the appointment period p, this is the n th delivery of truck k, (3) Decision variables: y k : if truck k is used, y k = 1, otherwise, y k = 0, y i,k : if truck k delivers containers to terminal i, y i,k = 1, otherwise, y i,k = 0, y in,k : if truck k delivers containers to terminal i, this is the n th delivery of truck k, y in,k = 1, otherwise, y in,k = 0, y ij,k : if truck k delivers containers form node i to node j, y ij,k = 1, otherwise, y ij,k = 0, y pi,k : if truck k delivers containers to terminal i during the appointment period p, y pi,k = 1, otherwise, y pi,k = 0.

The Optimization Goal Setting
(1) The fixed costs of trucks The fixed costs of trucks include the purchase costs of trucks, driver's salaries, maintenance costs and so on. The trucks that incur fixed costs mainly refer to those perform scheduling tasks, except idle trucks. Additionally, it is not related to the mileage and the number of terminals. The fixed costs C 1 can be expressed as: (2) Opportunity loss costs of trucks waiting outside the terminal The trucks wait outside the terminal without energy consumption but there are opportunity loss costs. The opportunity loss costs C 2 can be expressed as: (

3) The carbon emissions costs and fuel costs
Trucks emit CO 2 when trucks are idling or moving. According to literature [30], the distance traveled per volume unit of fuel is closely related to the truck's load. Referring to the calculation of fuel consumption costs in literature [30], this paper assumes that the weight of a truck is Q 0 when the truck is empty. One 40 ft container is heavier as two 20 ft containers. The weight of each 20 ft container is same. The maximum weight that a truck can carry is Q. The fuel consumption per unit distance is β 1 when truck is empty. The fuel consumption per unit distance is β 2 when truck is at full load. The fuel consumption of per unit distance can be expressed as: are the coefficients of unitary linear expressions. Q x is the weight that the truck carried in the truck transportation process. The fuel consumption per unit distance in the truck transportation process is β(Q x ). The fuel consumption per unit distance can be expressed as β 1 = aQ 0 + b when truck is empty. The fuel consumption per unit distance can be expressed as β 2 = a(Q 0 + Q) + b when a truck is at full load. It follows that a = Then the fuel consumption of per unit distance in the truck transportation process can be expressed as: The calculation formula of EM 1 (the number of CO 2 emissions that trucks emit when trucks deliver containers form node i to node j) can e expressed as: The load of trucks in the transport process from node i to node j is Q ij and the distance from node i to node j is d ij and the fuel consumption of per unit distance in the truck transport process from node i to node j is β(Q ij ).
The truck is idling when truck loads containers and queues at the terminal gate and operates inside the terminal. Truck emits CO 2 when it is idling. The calculation formula of EM 2 (the number of CO 2 emissions that trucks emit when they are idling) can be expressed as: This paper introduces carbon tax mechanism to calculate carbon emissions costs. If the carbon tax is c o , the carbon emissions costs C 3 can be expressed as: The fuel costs C 4 can be an expression as: Sustainability 2019, 11, 6256 9 of 23

Optimize Model Settings
In Section 3.4.1, fixed costs, opportunity loss costs, carbon emissions costs, and fuel costs are comprehensively considered. The corresponding optimization model is given by the following: T pin,k = T d pn,k + d 0i /v, ∀p ∈ P, ∀i ∈ I, ∀k ∈ K, ∀n ∈ N T m in,k = T pin,k + t g pi + t o pi , ∀p ∈ P, ∀i ∈ I, ∀k ∈ K, ∀n ∈ N T r pn,k = T m in,k + d 0i /v, ∀p ∈ P, ∀i ∈ I, ∀k ∈ K, ∀n ∈ N T pjn,k = T m in,k + d ij /v, ∀p ∈ P, ∀i, j ∈ I, ∀k ∈ K, ∀n ∈ N (25) T e pj − T pjn,k , T e pj > T pjn,k 0, T e pj ≤ T pjn,k , ∀p ∈ P, ∀j ∈ I, ∀k ∈ K, ∀n ∈ N T m jn,k = T pjn,k + t w jn,k + t g pj + t o pj , ∀p ∈ P, ∀ j ∈ I, ∀k ∈ K, ∀n ∈ N T pi1,k ≥ T e i , ∀p ∈ P, ∀i ∈ I, ∀k ∈ K T pin,k ≤ T l i , ∀p ∈ P, ∀i ∈ I, ∀k ∈ K, ∀n ∈ N (29) y in,k ∈ {0, 1}, ∀i ∈ I, ∀k ∈ K, ∀n ∈ N (32) The objective of our problem is to minimize the total costs in objective function (8). The total costs include fixed costs, opportunity loss costs, carbon emissions costs, and fuel costs. Constraints (9)- (15) represent the relationship between variables. Constraints (16) and (17) indicate that the number of containers delivered to the terminal i by all trucks during the appointment period p should less the appointment quotas of the appointment period p of terminal i. Constraints (18) and (19) indicate the moment when truck k departures the outer yard during the appointment period p, it is the first delivery of truck k. T r (p−1)n,k indicates the moment when truck returns to the outer yard during the appointment period p − 1, it is the n th delivery of truck k. t y p1,k indicates the time that truck k waits at the outer yard before it departures the outer yard during the appointment period p, it is the first delivery of truck k. The difference between constraint (18) and (19) is that constraint (18) indicates that truck k only loads one container, and constraint (19) indicates that truck k loads two containers. Constraints (20) and (21) indicate the moment when truck k departures the outer yard during the appointment period p, it is the n th delivery of truck k. T r p(n−1),k indicates the moment when truck k returns to the outer yard during the appointment period p, it is the n − 1 th delivery of truck k. The difference between constraint (20) and (21) is that constraint (20) indicates that truck k only loads one container, and constraint (21) indicates that truck k loads two containers. Constraints (22) indicates the moment when truck k arrives at the terminal i during the appointment period p, it is the n th delivery of truck k. Constraints (23) indicates the moment when truck k departures the terminal i, it is the n th delivery of truck k. Constraints (24) indicates the moment that truck k returns to the outer yard form terminal i during the appointment period p, it is the n th delivery of truck k. Constraints (25) indicate the moment that truck k arrives at the terminal j from terminal i during the appointment period p, it is the n th delivery of truck k. Constraints (26) indicates the waiting time of truck k outside the terminal j, it is the n th delivery of truck k. If the moment of truck k arrived at terminal j earlier than the moment of terminal j begins to accept containers during the appointment period p, the waiting time t w jn,k = T e pj − T pjn,k , otherwise, t w jn,k = 0. Constraints (27) indicate the moment that the truck k departures terminal j, it is the n th delivery of truck k. Constraints (28) and (29)

Hybrid Genetic Algorithm with Variable Neighborhood Search
Genetic algorithm (GA), also known as evolutionary algorithm, was first implemented in the year 1984 [31]. The GA is a kind of highly parallel and strongly robust global search algorithm based on the principles of biological evolution in nature. It has disadvantages such as prematurity and slow convergence speed. Variable neighborhood search (VNS) contains a variety of neighborhood structures, which can effectively combine with other heuristic algorithms to conduct targeted local search. In this paper, the hybrid genetic algorithm with a variable neighborhood search (HGAVNS) was designed based on the advantages of GA and VNS. The algorithm flow is shown in Figure 2.

Coding and Generation Initial Population
This paper studies the scheduling problem of trucks between the outer yard and multi-terminals.

Fitness Function
The fitness function of chromosomes can be expressed as:

Coding and Generation Initial Population
This paper studies the scheduling problem of trucks between the outer yard and multi-terminals.

Coding and Generation Initial Population
This paper studies the scheduling problem of trucks between the outer yard and multi-terminals.

Fitness Function
The fitness function of chromosomes can be expressed as:

Fitness Function
The fitness function of chromosomes can be expressed as: C R is the objective function value of chromosome R, f R is the fitness function value of the chromosome R.

Selection Strategy
The individuals are arranged according to the fitness value. The individual with the highest fitness was marked 0. The individual with the lowest fitness was marked pop_size − 1. It generates a random number ε following a normal distribution. It has a 99.9% possibility of greater than -3 and smaller than 3. If ε is greater than 3, a new random number ε is generated once again. Make ε * = ε 3 , ξ = ε * (pop_size − 1). Select ξ individual to enter the new population for the next iteration.

Variable Neighborhood Search Strategy
In the traditional genetic algorithm, individuals are disturbed by crossover and mutation operations. The algorithm's ability to search the solution depth is insufficient and the computational efficiency is low. This paper introduces variable neighborhood search algorithm to conduct depth search on each individual and design adaptive evolutionary pressure, adaptive neighborhood search range and adaptive neighborhood search times (see Section 4.5 for setting method). It can reduce the neighborhood search scope and improve search efficiency and balance the breadth and depth of search. The pseudo-code of the variable neighborhood search algorithm is given below:

Algorithm 1: Variable Neighborhood Search Algorithm
Begin 1: set the neighborhood structure N k = {N 1 , N 2 , . . . , N l }, N l is the l th neighborhood structure, 2: loop max_N times, 2.1: Individual R disturbed from the first neighborhood structure iter ← iter + 1 }, 2.4 until (iter = max_iter), 2.5 Individual R continues to be disturbed by the next neighborhood structure N l , iter ← 1 , 2.6 repeat, 2.7 until (iter = max_iter), 2.8 until (the optimal solution hasn't update when individuals were disturbed by the last neighborhood structure or neighborhood search has loop max_N times), End.

Adaptive Search Strategy
Due to the different pressures required in different stages of population evolution, adaptive evolutionary pressures, adaptive neighborhood search times strategy and adaptive neighborhood search range strategy are introduced in this paper to balance the breadth and depth by changing the population search range and neighborhood search times. The adaptive search strategies in this paper are as follows: (1) Adaptive evolutionary pressure The adaptive evolutionary pressure in this paper is as follows: ( ) (2) Adaptive neighborhood search times strategy The adaptive search times strategy in this paper is as follows: The adaptive search range strategy in this paper is as follows:

Adaptive Search Strategy
Due to the different pressures required in different stages of population evolution, adaptive evolutionary pressures, adaptive neighborhood search times strategy and adaptive neighborhood search range strategy are introduced in this paper to balance the breadth and depth by changing the population search range and neighborhood search times. The adaptive search strategies in this paper are as follows: (1) Adaptive evolutionary pressure The adaptive evolutionary pressure in this paper is as follows: λ is the evolutionary pressure, ∂ is the coefficient set to prevent the evolutionary pressure at the later stage of evolution from being 1(stop the evolution) and max_gen is the preset maximum number of iterations. gen indicates algorithm has been iterated gen times. This paper set ∂= 0.8.
(2) Adaptive neighborhood search times strategy The adaptive search times strategy in this paper is as follows: S t is the neighborhood search time. α 1 represents the minimum search times and α 2 is the adaptive search times. max_gen is the preset maximum number of iterations. gen indicates algorithm has been iterated gen times. represents round down. It can be seen from the equation (37) that in the initial stage of evolution, the population has a low number of search times α 1 . With the continuous iteration of the population, the search times gradually increase. The maximum search times is α 1 + α 2 . On the basis of maintaining the diversity of the population, the depth of exploration should be increased.

(3) Adaptive neighborhood search range strategy
The adaptive search range strategy in this paper is as follows: S r is the neighborhood search range. max_gen is the preset maximum number of iterations. gen indicates algorithm has been iterated gen times.
represents round down. Q 1 is the demand for 40 ft containers of all terminals. Q 2 is the demand for 20 ft containers of all terminals. In the initial stage of evolution, the minimum neighborhood search range is 1 (Q 1 + Q 2 ). With the continuous iteration of the population, the neighborhood search range increases and the maximum search range is ( 1 + 2 )(Q 1 + Q 2 ). Obviously, with the increase in the number of iterations, the search scope gradually expands from 1 (Q 1 + Q 2 ) to ( 1 + 2 )(Q 1 + Q 2 ). In order to increase the diversity of the population and make it difficult for the population to fall into the local optimum, the poor solution is accepted with a certain bad solution acceptance criterion.

Iteration Stop Criteria
In the population iteration process, if it hasn't updated the optimal solution γ · max_gen times or the number of iterations has reached the preset maximum iteration number, the iteration will be stopped. This paper set γ = 0.1.

Time Complexity Analysis
For the HGAVNS, assuming the demand for 40 ft containers of all terminals is Q 1 , the demand for 20 ft containers of all terminals is Q 2 , the population size is pop_size, the maximum number of iterations is max_gen, the number of neighborhood searches is S n , the number of neighborhood cycles is max_N. According to the algorithm flow in Figure 2, calculate the complexity O(T i ) of each link i.

Description of Experimental Problems
A port in China is selected to study the truck scheduling problem. One outer yard has 100 trucks, each of which can carry one 40 ft container or two 20 ft containers. Seven terminals adjacent outer yard have established TAS, all of which have the demand for 20 ft containers and 40 ft containers. Each terminal divides a day into 12 appointment periods. The length of each appointment period is 2 h. The average driving speed of each truck is 60 km/h. The fuel consumption per unit distance is 0.8 L when the truck is empty. The fuel consumption per unit distance is 1.2 L when truck is at full load. Fuel consumption is 2.5 L/h when truck is idling. The CO 2 emissions coefficient is 2.65 kg/L. The carbon tax is 0.04 Yuan/kg. The fuel cost is 6.41 Yuan/L. The unit cost of the truck waiting outside the terminal is 40 Yuan/h. The fixed cost of each truck is 300 Yuan. The average loading time of a container at the outer yard is 3 min. The distance between the outer yard and each terminal is shown in Table 1. The delivery information of each terminal is shown in Table 2. The appointment quotas in different appointment periods of each terminal are shown in Table 3. The average waiting time of truck at the terminal gate and the average operation time inside each terminal are different in different appointment periods of each terminal. In this paper, time interval is divided by 1 h and the changing trend is shown in Figures 5 and 6. For ease of calculation, the unit of time is expressed in minutes.  0  1  2  3  4  5  6  7   0  0  20  40  30  25  40  25  35  1  20  0  25  15  10  25  15  15  2  40  25  0  15  15  20  10  15  3  30  15  15  0  25  10  15  10  4  25  10  15  25  0  15  25  10  5  40  25  20  10  15  0  15  20  6  25  15  10  15  25  15  0  25  7  35  15  15  10  10  20 25 0 Table 2. Delivery information of each terminal.

Parameter Setting and Results Analysis
This paper uses MATLAB R2018b to implement the proposed algorithm and model, and all experiments in this paper are evaluated on PCs with Intel(R) Core (TM) i7-7700 CPU @3.60GHz and the memory of 16.0GB. The parameters setting of the algorithm is shown in Table 4.

Parameter Setting and Results Analysis
This paper uses MATLAB R2018b to implement the proposed algorithm and model, and all experiments in this paper are evaluated on PCs with Intel(R) Core (TM) i7-7700 CPU @3.60 GHz and the memory of 16.0 GB. The parameters setting of the algorithm is shown in Table 4. The model is solved by HGAVNS, and the optimal solution is 66642.83. In the optimal scheme, the scheduling sequence is shown in Appendix A Table A1. In Table A1, the order 1-14 in the first row represents the delivery order of relevant truck. This paper assumes a represents the appointment period of the truck for this delivery. b represents the number of containers loaded by the truck when departures from the outer yard. c represents the serial number of the terminal. d represents the time that returns to the outer yard after the previous delivery. e represents the waiting time before this delivery at the outer yard. The first nine trucks deliver 40 ft containers. The delivery and operation conditions of these trucks are marked by 1 × 4 array. The 1 × 4 array's form is (a, c, d, e). The last five trucks deliver 20 ft containers. The delivery and operation conditions of these trucks are marked by a 1 × 5 array. The 1 × 5 array's form is (a, b, c, d, e).
The data in Table 5 presents the transportation mileage, working time, carbon emissions and the returning time in the last delivery of each truck under TAS. d represents the transportation mileage of each truck. T l represents the total loading time of each truck. T w represents the total waiting time at outer yard of each truck. T q represents the total queue time at terminal gate of each truck. T o represents the total operation time at terminal of each truck, including unloading time and transportation time at terminal. EM represents the number of carbon emissions of each truck. T r represents the returning time in the last delivery of each truck. C represents the total costs.

Sensitivity analysis
The goal of this paper is to reduce carbon emissions and total costs. The speed and fixed cost of each truck is set in Section 5.1. To discuss the effect of these variables on the optimal solution, the following sensitivity analysis is done in this paper. The data in Table 6 is the average of 10 runs of the program. The data in Table 6 presents the number of trucks were used, the number of CO 2 emissions, the total costs respectively. N represents the number of trucks were used, W represents the number of CO 2 emissions, C represents the total costs. Figure 7 presents the changing trend of total costs under different speeds and fixed costs of trucks. According Table 6 and Figure 7, we can find that as the fixed costs increase, so does the total costs. Under a fixed cost of the truck, as the speed of the truck increases, the number of trucks used decreases and the total costs decrease. In order to ensure the safety of transportation, the speed of each truck is usually controlled between 55 and 65(unit: km/h). This paper takes the intermediate value and assumes the speed of the truck is 60 km/h.  According Table 6 and Figure 7, we can find that as the fixed costs increase, so does the total costs. Under a fixed cost of the truck, as the speed of the truck increases, the number of trucks used decreases and the total costs decrease. In order to ensure the safety of transportation, the speed of each truck is usually controlled between 55 and 65(unit: / km h ). This paper takes the intermediate value and assumes the speed of the truck is 60 / km h .

Comparative Analysis of Calculation Examples
Based on the HGAVNS solves the scheduling scheme of trucks between outer yard and multiple terminals considering carbon emissions under TAS. To verify the effectiveness of TAS, the following comparative experiment was conducted in this paper. Assume that the terminal establishes TAS, the outer yard doesn't adopt TAS. Appendix Table A2 shows the scheduling scheme without TAS. In Table A2, the order 1-13 in the first row represents the delivery service order of relevant trucks. This paper assumes f represents the waiting time of the truck outside the first terminal. g represents the waiting time of the truck outside the second terminal. The first nine trucks and the first five delivery services of the tenth truck diver 40 ft containers. The delivery and operation conditions of these trucks are marked by 1×4 array. The 1×4 array's form is ( ) a c d f . The last delivery service of the tenth truck and the last five trucks diver 20 ft containers. Which were marked with shadow are the operating conditions of a truck loading two 20 ft containers going to two terminals at the same

Comparative Analysis of Calculation Examples
Based on the HGAVNS solves the scheduling scheme of trucks between outer yard and multiple terminals considering carbon emissions under TAS. To verify the effectiveness of TAS, the following comparative experiment was conducted in this paper. Assume that the terminal establishes TAS, the outer yard doesn't adopt TAS. Appendix A Table A2 shows the scheduling scheme without TAS. In Table A2, the order 1-13 in the first row represents the delivery service order of relevant trucks. This paper assumes f represents the waiting time of the truck outside the first terminal. g represents the waiting time of the truck outside the second terminal. The first nine trucks and the first five delivery services of the tenth truck diver 40 ft containers. The delivery and operation conditions of these trucks are marked by 1 × 4 array. The 1 × 4 array's form is (a, c, d, f ). The last delivery service of the tenth truck and the last five trucks diver 20 ft containers. Which were marked with shadow are the operating conditions of a truck loading two 20 ft containers going to two terminals at the same time. These are marked by 1 × 6 array. The 1 × 6 array's form is ([a, a], b, [c, c], d, f , g). Which weren't marked with shadow are the operating conditions of a truck loading one or two 20 ft containers going to one terminal. These are marked by a 1 × 5 array. The 1 × 5 array's form is (a, b, c, d, f ).
The data in Table 7 presents the transportation mileage, working time, carbon emissions and the returning time in the last delivery of each truck without TAS. T w1 represents the total waiting time outside the first terminal. T w2 represents the total waiting time outside the second terminal.
The  In terms of the number of trucks used in the scheduling scheme, the new scheduling model under TAS is 14. Use one less truck than the traditional scheduling scheme without TAS. Under TAS, the total costs were reduced by 2.73%, and the carbon emissions were reduced by 0.22%.
In terms of the working time, that the total queue time of the trucks decreases from 775.5 min to 699 min after TAS is adopted. The maximum queue time decreased from 84.7 min to 72 min. The average queue time decreases from 51.7 min to 49.93 min. The total operation time of the trucks decreases from 737.3 min to 558.1 min after TAS is adopted. The maximum operation time decreased from 68.2 min to 59.4 min. The average operation time decreases from 49.15 min to 39.86 min. It can be seen that the TAS can effectively reduce the queue time and the operation time, improve work efficiency. Figures 8 and 9 respectively compare the queue time of trucks and the operation time at terminal under TAS and without TAS.     In terms of the waiting time, under TAS, the waiting time of a truck is mostly generated at the outer yard. When trucks wait at the outer yard, there are no costs. Without TAS, trucks haven't appointment information communication, trucks will depart from the outer yard immediately after the last delivery completed. However, often it arrived at the terminal due to terminal in a variety of reasons need to wait for a long time outside the terminal, there will be opportunity loss costs.
Compared with the traditional scheduling scheme, the scheduling scheme under TAS need one less truck, improves the delivery efficiency of trucks, effectively reduces energy consumption, and finally reduces the total costs of the truck company.

Conclusion
In this paper, the scheduling problem of trucks between the outer yard and multiple terminals under the TAS is solved by a hybrid genetic algorithm with a variable neighborhood search based on the corresponding optimization model and the consideration of carbon emissions. The numerical example and comparative analysis show the validity of the proposed model and algorithm. The truck fleet can be reasonably scheduled to complete the delivery tasks during each appointment period. The number of trucks used under TAS is one less than the traditional scheduling scheme, it improves the work efficiency of trucks. Fuel consumption and carbon emissions can be reduced greatly. Finally, the total costs of a truck company can be reduced. The research results expand the theoretical research on truck scheduling between the outer yard and multi-terminals. In terms of the waiting time, under TAS, the waiting time of a truck is mostly generated at the outer yard. When trucks wait at the outer yard, there are no costs. Without TAS, trucks haven't appointment information communication, trucks will depart from the outer yard immediately after the last delivery completed. However, often it arrived at the terminal due to terminal in a variety of reasons need to wait for a long time outside the terminal, there will be opportunity loss costs.
Compared with the traditional scheduling scheme, the scheduling scheme under TAS need one less truck, improves the delivery efficiency of trucks, effectively reduces energy consumption, and finally reduces the total costs of the truck company.

Conclusions
In this paper, the scheduling problem of trucks between the outer yard and multiple terminals under the TAS is solved by a hybrid genetic algorithm with a variable neighborhood search based on the corresponding optimization model and the consideration of carbon emissions. The numerical example and comparative analysis show the validity of the proposed model and algorithm. The truck fleet can be reasonably scheduled to complete the delivery tasks during each appointment period. The number of trucks used under TAS is one less than the traditional scheduling scheme, it improves the work efficiency of trucks. Fuel consumption and carbon emissions can be reduced greatly. Finally, the total costs of a truck company can be reduced. The research results expand the theoretical research on truck scheduling between the outer yard and multi-terminals.
In terms of practical application, the case analysis of terminal in China can provide some basis for the implementation of TAS and provide some references for the optimization decision of the truck company on truck scheduling.
It should be pointed out that because of the complexity of various links and influencing factors involved in the delivery, the queue rules of trucks at the terminal gate, relevant interference factors and trucks missing the appointment in the process of delivery were not considered in the solution and analysis, which should be further studied in the future.

Conflicts of Interest:
The authors declare no conflict of interest.