Optimization of the Waterbus Operation Plan Considering Carbon Emissions: the Case of Zhoushan City

Recently, as more people are concerned with the issues around environment protection, research about how to reduce carbon emissions has drawn increasing attention. Encouraging public transportation is an effective measure to reduce carbon emissions. However, overland public transportation does less to lower carbon because of the gradually increasing pressure of the urban road traffic. Therefore, the waterbus along the coast becomes a new direction of the urban public transport development. In order to optimize the operation plan of the waterbus, a bi-level model considering carbon emissions is proposed in this paper. In the upper-level model, a multiple objective model is established, which considers both the interests of the passengers and the operator while considering the carbon emissions. The lower-level model is a traffic model split by using a Nested Logit model. A NSGA-II (Non-dominated Sorting Genetic Algorithm-II) algorithm is proposed to solve the model. Finally, the city of Zhoushan is chosen as an example to prove the feasibility of the model and the algorithm. The result shows that the proposed model for waterbus operation optimization can efficiently reduce transportation carbon emissions and satisfy passenger demand at the same time.


Introduction
Climate change resulting from the emission of carbon dioxide associated with human economic activities has the potential to create an unstable global climatic future [1].Nowadays, carbon emission has become a global issue.According to the IEA (International Energy Agency), 23% of the global energy-related carbon emission comes from the transport sector [2].Additionally, that figure will grow to 41% by 2030 [3].Vehicles are the biggest carbon emitter in urban traffic.A study by Barth [4] has shown that carbon emissions of vehicles are the highest in cases of stop-and-go traffic and high-speed situations.It is clear that when road congestion occurs, all the vehicles will be in a stop-and-go situation, and the carbon emissions of vehicles will be higher accordingly.However, developing waterbuses in coastal and riverside cities can both alleviate road pressure and reduce the carbon emissions of the transport sector.Furthermore, the waterbus needs little investment because it can operate using existing water resources without building roads.In addition, though the waterbus has opened in some cities, the poor level of waterbus service due to bad arrangements leads to declining traffic.Therefore, the existing waterbus cannot compete with overland transport to share the traffic of vehicles.
Because of the reasons above, this paper proposes an optimization of the waterbus operation plan in terms of carbon emissions.Part of overland traffic is transferred to the waterbus by improving the waterbus service level, which can alleviate the road pressure and reduce the carbon emissions of the transport sector.However, massive expenditures need to be spent to improve the waterbus service level.In order to improve the feasibility of the waterbus operation plan, the toll revenue of the road is used to subsidize the waterbus in order to alleviate the pressure of funding.At the same time, raising the traffic tolls increases the trip cost of cars and somewhat curbs the usage of cars, therefore encouraging people to make greater use of the waterbus.Thus, the target to cut emissions of the transport sector would be reached by the above method.
Compared with the optimization of overland bus operation plans, the waterbus is much simpler because of its smaller site scale and lack of interference from other transportation modes.Furthermore, unlike the road, the sea in which waterbuses sail has no capacity limitations [5].However, on account of the narrow application range of the waterbus (limited to coastal and riverside cities), there are few studies about the waterbus.Ye et al. [6] used a disaggregate model to analyze the attractiveness of the waterbus with Shanghai and Huai'an as examples.Sato and Takadama [7] studied the waterbus line planning problem using a Pittsburgh-style learning classifier system (LCS) and researched a way that ensured the robustness of the waterbus system to make the waterbus continue to serve passengers even in bad weather.Li and Huang [8] studied the characteristics of Guangzhou waterbus passengers based on IC (Integrated Circuit) cards, and they drew a conclusion that the waterbus is more attractive to non-commuters than commuters.Keiki and Keiji [9] studied the classifier system for multiple environments: towards a robust waterbus route for several situations, and in the same year, he found a new method to optimize the waterbus route through the Pittsburgh-style learning classifier system [10].
The study of the overland bus is more sophisticated than the waterbus, so it can be a reference for the study of the waterbus.Ortuzar and Willumson [11] presented the most important transport modeling techniques and they have approached the subject from the point of view of a modeling exercise, discussing the role of theory, data, model specification in its widest sense, model estimation, validation, and forecasting.Cascetta [12] provided a comprehensive and systematic presentation of the mathematical models for the simulation of transportation systems.Methodologies for the analysis and design of these systems were provided as well.Lam and Bell [13] addressed the current problem of improving public transit systems by taking advantage of new technologies and advanced modeling techniques.The key areas open to improvement were service planning and operations management.Russo et al. [14] estimated the target time of trucks that transport agriculture and food products subject to commercial and production constraints.Vitetta et al. [15] investigated the shortcomings of the frequency approach using optimal strategy compared with the schedule-based approach in simulating ex-urban transit systems in terms of user behavior simulation and subsequently verifying a numerical application.Bloom and Mathew [16] thought that the main emphasis of bus transit route network optimization was bus routes and schedules.The passenger demand was satisfied and the cost of the bus company was minimized by optimizing bus routes and schedules.Dubois et al. [17] optimized a bus transit route network considering factors such as geography, capacity, direct demands, and transfers.Baaj and Mahmassani [18] designed a route-generation algorithm (RCA) to optimize bus lines and frequency.Agrawal and Mathew [19] set up an optimization model of a bus transit route network, in which the objective was to minimize the transportation costs, and calculated the model with the genetic algorithm (GA).Yu et al. [20] proposed models to predict bus arrival times at the same bus stop but with different routes and several methods, which included support vector machine (SVM), artificial neural network (ANN), k nearest neighbors algorithm (k-NN), and linear regression (LR).Shih [21] and Chakroborty et al. [22] proposed the optimization models to optimize the fleet size distribution and scheduling.Cordeau [23] presented recent optimization models for the most commonly studied rail transportation problems, which concentrated on routing and scheduling problems.Similar research problems have been studied by Chakroborty [24], Yao et al. [25,26], and Yu et al. [27].
Referring to existing research and studies, this paper studied the waterbus operation plan problem in three aspects: environment, passengers, and waterbus operators.A bi-level model is proposed in this paper.In the upper-level model, a multiple objective model is established, which considers the benefits of both the passengers and the operator while considering the carbon emissions.The lower-level model is a traffic model split by using a Nested Logit model.After calculation, we expect to find a balance point among passenger demand, a low cost for the operators, and low carbon emissions of the transport sector.Finally, the trip between Zhoushan and Ningbo cities is chosen as an example to test the model.This paper is organized as follows: the second section describes the problems and introduces the background of the city of Zhoushan.A bi-level model is proposed in the third section.The fourth section introduces the NSGA-II algorithm and, in the fifth part, the city of Zhoushan is chosen as a setting to test the model.The conclusion can be found in the last part.

Problem Description
As is shown in Figure 1, there is both road-based and marine transportation between the cities of Zhoushan and Ningbo.People can drive across the cross-sea bridge or take the waterbus, which includes two types of ships, the passenger ship and the ro-ro passenger ship.This paper aims at transferring part of the overland traffic to the waterbus by improving the waterbus service level.The transfer traffic can alleviate the road pressure and reduce the carbon emissions of the transport sector.However, massive expenditures need to be spent to improve the waterbus service level.In order to improve the feasibility of the waterbus operation plan, the toll revenue of the cross-sea bridge is used to subsidize the waterbus to alleviate the pressure of funding.At the same time, raising the traffic tolls increases the trip cost of cars and somewhat curbs the usage of cars, which encourages people to make greater use of the waterbus.Due to the improvement of the waterbus service level and the rising toll of the cross-sea bridge, an increasing number of people enjoy traveling by waterbus while fewer people drive across the cross-sea bridge.It is expected to strike a balance between the waterbus and the cross-sea bridge which meet the passenger demand, provide a low cost for the operators, and lower the carbon emissions of the transport sector.

Background
The city of Zhoushan is located on an archipelago in China, which has 1390 islands and reefs.Therefore, Zhoushan is known as "the city of a thousand islands".Zhoushan Island, with an area of 502 square kilometers, is the biggest island in the city of Zhoushan.The permanent population of Zhoushan Island is 700,000, which accounts for 62.5 percent of the population of Zhoushan.The city of Ningbo is one of the most important cities that connect Zhoushan with the outside world.As Figure 2 shows, there are two waterbus lines between Ningbo and Zhoushan Island proper: the Baifeng passenger terminal of Ningbo to Yadanshan passenger terminal of Zhoushan Island is Line 1; the Baifeng passenger terminal of Ningbo to Shenjiamen passenger terminal of Zhoushan Island is Line 2. Additionally, there is also a cross-sea bridge connecting Ningbo to Zhoushan Island.
The data of this paper comes from the survey undertaken in Zhoushan via Dalian Maritime University in September 2012.According to the survey, origin-destination (OD) traffic volume was 33,928 from Ningbo to Zhoushan Island, in which there were 15,448 cars across the cross-sea bridge.Saturation of the cross-sea bridge was about 0.7.The average saturation of the expressway is 0.4, so the traffic congestion occurred on the cross-sea bridge.The service capacity of the cross-sea bridge cannot meet traffic demands as vehicle ownership continues to increase.Although the waterbus has opened between Ningbo and Zhoushan Island, the waterbus lost its attraction to passengers due to the poor service level.Therefore, the existing waterbus between Ningbo and Zhoushan Island cannot compete with the cross-sea bridge to share the traffic of vehicles.
There are two kinds of people traveling between Ningbo and Zhoushan Island.The first are local residents who work or study in Ningbo, and the trip modes of these people are taking the waterbus (Line 1 and Line 2) and driving across the cross-sea bridge.The other people are tourists who travel to the Putuo District of Zhoushan Island, and they usually drive across the cross-sea bridge or take the waterbus (Line 2).This paper studies the waterbus operation plans of Line 1 and Line 2. As most people tour on the non-workday, the waterbus operation plan of Line 2 is optimized on the workday and non-workday, respectively.

Optimization Model for the Waterbus Operation Plan
A bi-level model is proposed in this paper.The upper-level model is a multiple objective model [28], and the decision variables are the numbers of waterbuses, the departure frequency, and the fares.Objective equation 1 is the minimum of the total trip costs of passengers.Objective equation 2 is the minimum of the operation costs of operators.Minimizing the carbon emission costs of both waterbuses and cars on the cross-sea bridge is objective equation 3. The lower-level model is a traffic model split by using a Nested Logit model [29].The attributes under examination are the trip time costs and fares associated with each available mode.We first present a formal definition of the problem.Please see Table 1 for a summary of the notation used in this paper.

Upper-Level Model
Objective Equation (1): where t d is the unit time cost.As the values of t d and emission costs are difficult to measure and it is also time-consuming to measure unit time cost and emission cost, we referred to the paper of Yu et al. [30].μ is the conversion coefficient for the cars.μ is 1.5 in this paper, which denotes that the average people in a car is 1.5; and T represents the base period of the calculation.In this paper, one year is the base period.
The penalty cost is considered because passengers should transfer over land traffic modes to go to the destination after going ashore from the waterbus.This paper assumes that passengers take a taxi to go to the destination.
The objective Equation ( 1) is the minimum of the total travel cost of people.The first part of the equation is the travel costs of passengers who take the waterbus, which consists of the waterbus traveling time cost (waterbus traveling time cost and the passengers' waiting time cost), the fares of the waterbus, and the fares of taxi.There are parking lots and riding facilities in the passenger terminals of the waterbuses.If passengers go to the terminal by car and do not want to take the ro-ro passenger ship, they can park in the parking lot.After these passengers go ashore, they can use other modes of transportation, such as taxis, buses, and so on.If passengers drive to the terminal and want to arrive at the destination in their cars, they can take the ro-ro passenger ship.The second part of the equation is the travel costs of people who drive across the cross-sea bridge, which consists of the traveling time cost (drive), the toll of the cross-sea bridge, and the costs of fuel.
Objective Equation ( 2): ( ) The objective Equation ( 2) is the minimum of the operation cost of the operator, which consists of the depreciation of the waterbus and the cost of fuel consumption minus the subsidy from the toll revenue of the cross-sea bridge.
Objective Equation (3): (3a) ( ) According to IPCC (Intergovernmental Panel on Climate Change), the value of w Eq is 3.17 tons of carbon dioxide generated by one ton of fuel [31].
Eq is a equation of the speed [32], where 1 η , 2 η , 3 η , 4 η are constant terms and ε is the conversion coefficient of the emission cost.
The objective Equation ( 3) is the minimum of the carbon emission cost, which includes the carbon emission cost of passengers taking the waterbus and people driving across the cross-sea bridge.Additionally, the carbon emission cost of passengers is made up of the waterbus carbon emission cost and the taxi carbon emission cost.

Lower-Level Model
When choosing travel modes, passengers are usually affected by many factors, such as traveling time, cost, security, and comfort.After weighing all the influence factors, passengers will select the most advantageous mode for their own travel [33].The traffic flow of each route for each travel mode can be forecasted by traffic model split and traffic assignment, which are the last two steps of the fourstage transportation forecasting model.The four-stage transportation forecasting model, which is based on the person trip survey and consists of trip generation/attraction, trip distribution, traffic model split, and traffic assignment, is a method to forecast traffic volume.
The trip time cost and the fares associated with each available mode are the attributes in this paper, which have impacts upon individuals' transport behaviors.There are two travel modes for people to choose: waterbus and driving across cross-sea bridge.If people take the waterbus, there are three types of ships for them to choose from: the passenger ship, the ro-ro passenger ship for passengers, and the ro-ro passenger ship for cars.As alternatives are not independent, a Nested Logit model [29] is used to solve the traffic assignment problem.As is shown in Figure 3, the alternatives are the two mode choices (waterbus and drive across the cross-sea bridge).Additionally, there are two kinds of waterbuses (passenger ship and ro-ro passenger ship).Passenger ships can only carry passengers and ro-ro passenger ships can carry passengers as well as cars at the same time.For the convenience of calculation, the ro-ro passenger ship is divided into two kinds (for people and for cars) because it can both carry passengers and cars.The specific equations is as follows: (5) The utility is an equation of the attributes associated with that alternative, which is as follows: ( ) where α , β , γ , θ, ' α , ' β are unknown constants, among which α and ' α are the utility coefficients for traveling time costs and β and '  β are the utility coefficients for fares, respectively; γ and θ are the utility coefficients of comfort.Constraint conditions are as follows: Equation ( 12) ensures that the number of waterbus passengers must be no greater than the capacity limit.Equation ( 13) means that the carbon emissions after optimizing cannot be greater than before.

Solution Algorithm
The bi-level model is a NP (Non-deterministic Polynomial)-hard problem that is mostly solved by heuristic algorithm.In this paper, the upper-level model is a multiple objective model optimizing the waterbus operation plan.The lower-level model is a traffic model split, which can be solved by general software.In this paper, we use the Lingo program to solve the parameters of the lower-level model.The multiple objective model of the upper-level model is solved using a NSGA-II algorithm.
The NSGA-II (Non-dominated Sorting Genetic Algorithm-II) algorithm proposed by Deb et al. [29,34] is one of the most efficient and famous multiple objective algorithms.The fast non-dominated sorting technique and a crowding distance are used to rank and select the population fronts in the algorithm.After that, the algorithm uses the standard bimodal crossover and polynomial operators to combine the current population and its offspring generated as the next generation.Lastly, the best individuals in terms of non-dominance and diversity are selected as the solutions.The specific process of the NSGA-II algorithm is shown as the following: Step 1. Generating the initial population Pt.Let t = 0.The decision variables are the waterbus number, the fares, and the toll of the cross-sea bridge.Therefore, the coding of the algorithm consists of the waterbus number, the fares, and the toll of the cross-sea bridge.The coding is shown in Figure 4.The whole coding is encoded by binary.The first part of the coding represents the waterbus number, the second part represents the waterbus fares, and the third part represents the toll of the cross-sea bridge.Step ( ) To explain the NSGA-II algorithm, a small numerical example is proposed in this paper.It is assumed that there are two objective equations: The values of 1 x and 2 x range from 1 to 5. When 1 x and 2 x get the minimal value at the same time, the objective equation 1 reaches the minimum.At the moment, the value of 1 x is 1 and the value of 2 x is also 1.When 1 x gets the minimal value and 2 x gets the maximal value, the objective equation 2 reaches the minimum.At this time, . The two objective equations cannot reach the minimum at the same time, so we use the NSGA-II algorithm to solve the problem.The solving process of the numerical example is as follows.
First, we generate the initial population 5 and transfer the values into binary format: 101 . After that, carry on crossover and mutation and generate new offspring: 101 100 . Place new offspring in the new population: 101 . Then, find the non-dominated individuals from the population: 100 and add them to the Pareto frontier.Finally, repeat the process and the maximum iteration number is set to 200 generations.The result of the small numerical example is shown in Figure 5.As shown in Figure 5, each point represents a set of decision variables and each point is an optimal solution.

Figure 5.
Pareto frontier of the small numerical example.

Case Study
In order to test the above model and algorithm, the trip between Ningbo and Zhoushan Island has been chosen as a research example.As Figure 2 shows, this paper optimizes the waterbus operation plan of Line 1 and Line 2.

Parameters of the Upper-Level Model
For Line 1 (from Baifeng of Ningbo to Yadanshan of Zhoushan Island) and Line 2 (from Baifeng of Ningbo to Shenjiamen of Zhoushan Island), the path length of the waterbus is 14.8 km and 27.4 km, respectively.The path length of the cross-sea bridge is 34.47 km.This paper assumes that the destination of the passengers is the center of Zhoushan Island for convenience of calculations.There is 9.1 km from Yadanshan to the center (Line 1) and 14.1 km from Shenjiamen to the center (Line 2).The distance from the cross-sea bridge to the center of Zhoushan Island is 13.5 km.
The vessel type and sailing speed are set in accordance with Zhoushan's existing ships: the rated number of passengers is 200, and the speed of the ships is 15 km/h.The existing operational ships in Zhoushan are used as waterbuses to utilize existing resources, and only the depreciation costs of these ships are considered.The fuel cost is 7.5 RMB/liter, or $7500/ton.Other relevant data were shown in Table 2.In the Nested Logit model, the unknown parameters α , β , γ , θ , ' α , ' β are determined by the individual point of view.By fitting and analyzing the existing data using the least squares method, the coefficient associated with each mode is calculated: α , and θ 0.47 = .Both the values of the coefficients associated with fares and travel time costs are negative, which indicates that the utility associated with the alternative decreases with respect to rises in fares or travel time costs.

Calculation Results
Figure 6 shows the Pareto frontier for triple-objective optimization (the upper-level model) obtained using the NSGA-II algorithm.The population size of the NSGA-II algorithm is 100 and the crossover rate is set as 0.7, the mutation rate is 0.1, and the target M is 3, which means that the number of the objective equation is three, and the maximum iteration number is 3000 generations.Each point in Figure 6 represents a set of decision variables, which include the number and fare of waterbuses as well as the toll of cross-sea bridge.The three coordinate values of each point are the values of the three objective equations, which are the travel costs of people, the operation costs of operators, and carbon emission costs of the transportation system.The x-axis represents the value of objective Equation (1).The value of objective Equation ( 2) is on the y-axis, and the value of objective Equation ( 3) is on the z-axis.All the points of the Pareto frontier are optimal solutions which can minimize the travel costs of passengers, the waterbus operation costs of operators, and the carbon emissions of the system, respectively.In addition, the optimal solutions are widely distributed, which guarantees the diversity of the optimal solutions.Thus, the policymakers are provided with many solutions to choose from.Considering the carbon emissions, the ideal solution to minimize the carbon emissions is selected.The results are shown in Tables 3-5.As was shown in Table 3, the best operation plan of Line 1 is 23 waterbuses, which includes 15 passenger ships and eight ro-ro passenger ships.The departure frequency and fares of the passenger ship are 8 min and 15 RMB for each person, respectively.The departure frequency of the ro-ro passenger ship is 20 min and the fare is 15 RMB for each person or 70 RMB for each car.
Table 4 shows the waterbus operation plan of Line 2 during the non-workday.The best operation plan of Line 2 during the non-workday is 28 waterbuses, including 18 passenger ships and 10 ro-ro passenger ships, and the departure frequency of the two kinds of waterbus is 12 min and 30 min, respectively.In addition, the fare of the passenger ship is 16 RMB for each person and the fare of ro-ro passenger ship is 16 RMB for each person or 80 RMB for each car.
The waterbus operation plan of Line 2 during the workday is shown in Table 5.The best operation plan of Line 2 during workday is 14 waterbuses, which includes 10 passenger ships and four ro-ro passenger ships.The departure frequency of the two kinds of waterbus is 28 min and 96 min, respectively.The fare of the passenger ship is 16 RMB for each person and the fare of the ro-ro passenger ship is 16 RMB for each person or 80 RMB for each car, which is the same as on the workday.After optimizing, the toll of the cross-sea bridge increases to 140 RMB for each car from the original 120 RMB for each car.
Figure 7 compares the traffic of each mode after optimizing and the traffic of each mode before optimizing.It can be seen in Figure 7 that the traffic of the cross-sea bridge decreases 32.7% after optimizing.It decreases to 15,468 people, which is equivalent to 10,312 cars.Additionally, the saturation of the cross-sea bridge is about 0.5 at that moment.Though the saturation after optimizing is still higher than the average of the expressway, it is observably lower than before.The traffic congestion of the cross-sea bridge is alleviated after optimizing.The traffic of the waterbus increases both on Line 1 and Line 2. The growth of the traffic on Line 1 is bigger than on Line 2, whether during the non-workday or the workday.The reason for this is that Line 1 is closer to the center of Zhoushan Island than Line 2 and the traffic of the cross-sea bridge is more likely to transfer to Line 1 as compared to Line 2. In addition, because Line 2 connects with the Putuo scenic zone, some travelers may choose to take the waterbus instead of driving across the cross-sea bridge after optimizing.Thus, the traffic of Line 2 during the non-workday obviously increases.
The Figure 8 compares the traffic of each mode after optimizing and the traffic of each mode before optimizing.After optimizing the waterbus operation plan, the carbon emission cost of the system decreases from 464,265.4 to 346,969.There is about a 25.3% decrease in the carbon emission cost and a 30.6% decrease in the travel cost of people.The reason for the decrease in the travel cost is that the fares of the waterbus are lower than the cross-sea bridge, and the more people that transfer to take waterbus, the less travel costs are paid from the system.At the same time, the operation cost of the operators increases 26.8% due to the increase of the waterbuses.

Conclusions
In order to cut the carbon emissions of the transport sector, this paper proposes a bi-level model that considers carbon emissions to optimize the waterbus operation plan.Meanwhile, the interests of both the passengers and the operators are also considered.The waterbus operation plan consists of departure frequency, fares, and the number of waterbuses on each line.Additionally, two kinds of ships are considered in this paper: passenger ships and ro-ro passenger ships.Moreover, the NSGA-II algorithm is used to solve the model.Finally, through the case study of Zhoushan Island, the results show that the proposed model in this paper can optimize the waterbus operation plan with the consideration of carbon emissions.After optimizing, 32.7% of overland traffic is transferred to the waterbus, which can alleviate the cross-sea bridge pressure and reduce the carbon emissions of the transport sector at the same time.In addition, the travel costs of passengers are lower than before, thus the plan is more likely to be accepted by the public.Though the operating costs increase after optimizing, toll revenue from the cross-sea bridge is used to subsidize the waterbus, and the income of the companies also increases as the traffic of the waterbus increases.Thus, the result of optimizing involves no harm for operators.
Considering the interests of passengers, operators, and the environment, the plan we propose in this paper is feasible.Nevertheless, funds are one of the biggest hurdles in this plan.In regard to environmental protection and sustainable development, the government can take some measures to improve the feasibility of the waterbus operation plan; for example, the government can increase the subsidy of the waterbus and make some preferential policies for waterbus operators.This paper has significant practical meaning for Zhoushan Island.

Figure 1 .
Figure 1.Trip modes of people.

Qd
ro-ro passenger ship (for people) 3 = i denotes ro-ro passenger ship (for cars) The line between Ningbo and Zhoushan k Average running speed of waterbus w i v Waterbus distance traveled through line k w ik s Cross-sea bridge distance traveled through line k c k s Waterbus wait time of passengers through line k w ik T Unit time cost t d Fares of the three kinds of waterbus through line k w ik C Cross-sea bridge toll of each car through line k c k C Penalty cost punish C Number of passengers of waterbus through line k w ik q Number of people across cross-sea bridge through line k c k q Average running speed of cars c v Number of people across cross-sea bridge through line k c k q Conversion coefficient for the cars μ Unit fuel consumption of cars c Unit cost of gasoline c Base period of the calculation T Departure frequency of waterbus of line k w ik f Number of waterbuses on line k w ik n Travel distance of passengers through line k k S Fare of taxi taken by passenger after going ashore from the waterbus Carbon emissions of waterbus on line k w ik EQ Carbon emission coefficient for waterbus or cars w Eq , c Eq Conversion coefficient of emission cost ε Traffic assignment of cross-sea bridge or waterbus of line k c k P , w k P Traffic assignment of the three kinds of waterbus of line k w ik P Utilities of waterbus or cross-sea bridge of line k w k V , c k V Utilities of three kinds of waterbus of line k w ik V Traveling time costs of the three kinds of waterbus of line k w ik c Traveling time costs of cross-sea bridge c k c Capacity limit of each kind of waterbus i E

Figure 3 .
Figure 3.The structure of the Nested Logit model.

5 )
crowding selection operator x > y if y x r r = and x d > y d ; (6) use new generated population for a further run of the algorithm; (7) if the established number of generations is reached, stop and return to the Pareto frontier [35] of the best solution in the current population.Otherwise, go to Step 2.

Figure 7 .
Figure 7. Contrast of traffic after optimizing with before optimizing.

Figure 8 .
Figure 8. Contrast of costs after optimizing with before optimizing.

Table 1 .
Table of notations.

Table 2 .
The related parameters for solving the model.

Table 4 .
Waterbus operation plan of Line 2 during the non-workday.

Table 5 .
Waterbus operation plan of Line 2 during the workday.