Research Article  Open Access
Integrated Optimization of Bus Line Fare and Operational Strategies Using Elastic Demand
Abstract
An optimization approach for designing a transit service system is proposed. Its objective would be the maximization of total social welfare, by providing a profitable fare structure and tailoring operational strategies to passenger demand. These operational strategies include full route operation (FRO), limited stop, short turn, and a mix of the latter two strategies. The demand function is formulated to reflect the attributes of these strategies, invehicle crowding, and fare effects on demand variation. The fare is either a flat fare or a differential fare structure; the latter is based on trip distance and achieved service levels. This proposed methodology is applied to a case study of Dalian, China. The optimal results indicate that an optimal combination of operational strategies integrated with a differential fare structure results in the highest potential for increasing total social welfare, if the value of parameter related to additional service fee is low. When this value increases up to more than a threshold, strategies with a flat fare show greater benefits. If this value increases beyond yet another threshold, the use of skipped stop strategies is not recommended.
1. Introduction
There has been a considerable increase in numbers for citizen travel, travel frequencies, and travel distances in the past few decades, which arose from the continuous increase in urban interspace and the development of city economies. Aside from the fact that the increased number of private cars could obviously not efficiently cover these increasing trips, they generate urban concerns such as traffic congestion and air pollution, as well as health and safety problems, and so on. Compared with the performances of private cars, public transit clearly can perform in an environmentally friendly, efficient, and sustainable manner [1, 2]. However, in practice, full advantage is not always taken of public transit resources to meet such increasing and imbalanced travel demand, due to the fact that a common full route operation strategy is provided. Such an operation situation could be efficiently improved by the use of other strategies tailored to the greatest degree possible to meet observed passenger demand. This is because the use of strategies opens up an opportunity to further save on vehicles while ensuring that the passenger load on each route segment does not exceed the desired occupancy. Furthermore, by applying strategies, it is expected that not only will passenger travel times and crowding times in vehicles be significantly reduced, but bus operational emissions will even decrease, due to the reduction in the numbers of “stop” and “go” locations which will significantly contribute to air quality and environmental improvement.
Fundamentally, transit operational strategies include short turn, limited stop, deadheading, express, and zonal service all of which were initially illustrated by Furth and Day [3]. For instance, short turn strategy is useful when high passenger demand areas require attention; it can perform shorter round trip time and, thereby, increase service frequencies in specific areas along a bus line (e.g., city center or a residential area). Limited stop serves only a subset of stops on a route; thus stops with high demand are more frequently served than those stops with low demand. This strategy has also demonstrated that it can significantly improve transit system efficiency and better accommodate capabilities and needs with respect to supply and demand. Further studies are of Furth [4] and Ceder [5] who investigated the short turn strategy based on analyzing operational characteristics of full and partial route segments. These studies considered a given timetable and in Ceder [5] also the minimization of the number of vehicles required. Other studies of Delle Site and Filippi [6], Tirachini et al. [7], and Cortés et al. [8] dealt with the minimization of the total social cost for users and operators.
For limited stop strategy, based on the requirement of minimum level service, Fu et al. [9] formulated a model to design a dynamic stopskipping strategy through minimizing total costs for users and operators. In their study, the headway allowed at any stop must not exceed double the dispatch headway. The study was extended by Liu et al. [10] and Chen et al. [11], considering random travel time and inserting deadheading strategy. Based on the same objective, Leiva et al. [12] allowed for two types of vehicles in the model for designing high frequency unscheduled limited stop strategy on a bus corridor. Considering running time variations resulting from the uncertainty of transit system operation, Tétreault and ElGeneidy [13] designed a limited stop strategy using archived automatic passenger counter and automatic vehicle location data. Meanwhile, based on these two types of data, ElGeneidy and SurprenantLegault [14] evaluated the implementation of the limited stop strategy and demonstrated a saving in running time. Recently, a mixed integer nonlinear model, the objective of which was maximization of consumer surplus, was proposed by Chiraphadhanakul and Barnhart [15] for determining limited stop strategy on the basis of given bus trips.
Most of the studies investigated operational strategies with different objective functions; it shows that constructing operational strategies is beneficial to improve the efficiency of the public transit operations system. However, these studies ignore the possible impact on passenger demand. The change of the demand is, therefore, important to consider in the optimization formulation. The use of strategies is linked to service levels changes for two groups of passengers: (a) improved service for passengers whose origins and destinations are served though some of their intermediate stops is skipped and (b) reduced service for those passengers who want to board or alight at skipped stops. Therefore, naturally, potential passengers on segments with high service levels are likely to find use of public transit attractive. Alternatively, it is likely that passengers who have used public transit will discontinue this practice due to the low levels of service offered. As a result, there is a fluctuation in demand and its distribution induced by the use of strategies. Even though Ulusoy et al. [16] and Ulusoy and Chien [17] formulated an estimated demand function to determine the demand share for each strategy, such as all stop, short turn, and express strategies, according to their attributes, it was still based on the assumption of invariable demand. Thus, this paper considers variation of demand in the proposed method to obtain optimal topological characteristics of each strategy.
There are two common function forms used to describe variation of demand in the previous studies. The first form is exponential, presented in such studies as Evans [18], Yang and Huang [19], Lam and Zhou [20], Yang and Kin [21], and Bellei et al. [22], Yao et al. [23], Jørgensen et al. [24], and Klier and Haase [25]. Works on the other form, the linear function, have been done by Chang and Schonfeld [26], Chien and Spasovic [27], Tsai et al. [28], Savage [29], Li et al. [30], Chen et al. [11], and Kim and Schonfeld [31]. In these studies, two function forms have been proved to significantly reflect fluctuation of demand by changes in generalized costs of trips, though they are formulated in different expressions. Elements of generalized cost for a trip have been explored in modal choice [32–34], normally consisting of invehicle time, waiting time, access time, fare, comfort, and so on., but much work in public transport operations planning considered travel time and fare elements, ignoring crowding in vehicles. With improved level of life passengers are prone to increase the value of comfort items. For example, crowding in vehicles increases the value of time of passengers and hence their generalized travel cost [35]. In this paper, a linear demand function is developed to reflect the combined effects of attributes of strategies, crowding, and bus fare on demand.
Fare is found to be a key factor to capture in operational strategies planning context. There are four main elements related to fare: fare policy, strategies, structures, and enforcement technologies. An interaction among these four elements was illustrated in a report by Fleishman et al. [36]. The fare structure is the combination of one or more fare strategies with specific fare levels [36–39]. The different types of fare strategies are grouped into two basic categories: flat and differentiation. The flat one requires all passengers to pay the same fare, as studied by Evans [18], Chien and Spasovic [27], and Kim and Schonfeld [31]. When differential fare strategy is applied, passengers pay fare dependent on one or more factors, including length of trip, time of day, quality of service. Most studies on differential fare strategy, such as, Yang and Kin [21], Tsai et al. [28, 40], and Yook and Heaslip [41], concentrated on distance based factors.
The differences between flat fare and differential fare were discussed in detail by Fleishman et al. [36]. The flat fare was recognized as the simplest and most convenient, though it ignored equity. In contrast, the differentiated fare could display social equity while its implementation required the use of hitech collection systems, for example, smart bus or strip card, because an efficient way was needed to reduce the times by which passengers pay differently, often by a small amount. There is no doubt that adopting a differential fare in implementing public transit services presents a huge obstacle. However, as new technologies become less and less costly, they have gained applicability for use in collecting bus fare onboard. The problem has become how to develop an attractive and profitable fare structure.
These previous studies mainly focused on the optimization of FRO strategy with a flat or a differential fare structure using distance based methods. This paper, as indicated above, aims to examine effects of both fare types on operational strategy optimization, for three basic strategies, FRO, limited stop, and short turn, and a mix of limited stop and short turn strategies. A fare table is constructed which can serve to present either a flat fare or a differential fare structure, in accordance with fare policy objectives. Considering the attributes of operational strategies, the differential fare structure is constructed based on travel distance and obtained service levels. Finally, an optimal combination of fares and strategies is determined by optimizing total benefits consisting of consumer surplus and operator profit, considering the variation of demand.
Following this introductory section, Section 2 systematically describes system characteristics, constructs a fare table, and formulates optimization models. Section 3 presents an indepth analysis of the results generated by applying the proposed models to a real life example of a single bus line. Conclusions and directions for future research are presented in Section 4.
2. System Characteristics and Problem Formulation
Considering a single and bidirectional bus line, stop locations, route zones, and greatest potential passenger demand along this bus line are given. The passenger demand is assumed to be sensitive to fare, waiting time, crowding time, and invehicle time. It is also assumed that passenger boarding and alighting have a significant effect on bus dwell times at stops. Thus, invehicle time depends on buses passing through the number of bus stops and passenger demand at these bus stops, as well as bus speed. Users are assumed to arrive randomly at stops. The distribution of bus arrivals at each stop presents Poisson distribution. We assume that passengers would not transfer between vehicles operated by different strategies on the single bus line because of high transferring cost considered. This study period is assumed to be a single period of one hour, such as the morning peak hour.
2.1. A Fare Structure Table
Fare is of considerable importance to bus operation and management. A low fare may attract more passengers but reduce operator revenue. Alternatively, a high fare may increase the revenue but reduce demand. Therefore, it is necessary to determine an attractive and profitable fare that benefits both users and operators. This paper defines a fare structure for a bus line, as shown in Table 1. This fare table consists of fares between each pair of all stops along a bus line, whether these two stops are successive or not. Suppose that there is a bus line with stops in each direction. A value of fare, , is defined for each pair of all stops along this bus line.
 
Note. , represent passenger origin and destination, respectively. 
A Flat Fare Structure Table. In Table 1, values of a fare table for all , are the same. This means that the bus operator charges all travelers a flat fare, regardless of travel distance, time of day, service quality, or trip route.
A Differential Fare Structure Table. In Table 1, the value of completely or partly differs for each pair of stops along the bus line. If the value of equals the value of , the resulting fare table is symmetrical.
The elements of a differential fare structure table vary in accord with differential fare strategies, such as distance based, time based, quality based, cost based, route based, and patron based strategies, as well as market based fare strategies. They depend on fare policy objectives. For instance, in order to attract more passengers to use public transit in offpeak periods, the transit service planner may consider varied discounts to lower fares during offpeaks, which could cause some users to shift their traveling times from peak to offpeak times. If the objectives of a fare policy are to increase the operator’s revenue and achieve social equity, other differential fare strategies such as cost based, quality based, and distance based strategies are considered when setting bus fares. Accordingly, a differential fare structure table will display hundreds of thousands of possible bus fares for varied purposes of public transit service optimization. However, in consideration of transit system operation characteristics, this work constructs a differential fare structure based on travel distance and obtained service levels.
2.2. Notation
In order to formulate the methodology as a programming model, the parameters and variables are defined in Notation.
2.3. Demand Formulation
In this study, operational strategies are constructed based on passenger demand. Meanwhile, strategy application generates fluctuations in passenger demand because of the beneficial attributes of these strategies, like high speed between successive stops, greater comfort, and so on. Passenger demand on a single bus line is considered variations with respect to waiting time, invehicle time, fare, and invehicle crowding time, when designing strategies to further improve transit system efficiency.
(a) Waiting Time. The waiting time for an OD pair depends on obtained operational strategies and their associated frequencies. The waiting time for pair can be expressed as follows:
The product term between two binary variables, , in (1) equals 1 only if strategy serves stops i and j, otherwise 0.
(b) InVehicle Time. The invehicle time that a user experiences on a bus consists of running time and all dwell times at each intermediate stop visited by the bus. The total dwell time of strategy l between stop i and stop j depends on passenger boarding and alighting times at each stop and a constant consisting of deceleration/acceleration as well as doors opening and closing at each stop. Passenger boarding and alighting times of strategy l at stop , respectively, denoted as and , are expressed as follows:
In addition, it is also assumed in this paper that the processes of boarding and alighting are simultaneous (different doors for boarding and alighting) and that boarding and alighting flows are independent of each other. Therefore, the larger one between passenger boarding and alighting times at a stop is dwell time of strategy serving this stop. Accordingly, dwell time of strategy between stop and stop , , is expressed as follows:Thus, the average invehicle time for pair which also equals the average travel time of a bus from stop i to stop j is formulated as follows:
(c) Flat Fare and Differential Fare
Flat Fare. Clearly, using a flat fare table indicates that all values of pairs are the same , which is expressed as follows:
Differential Fare. In this study, a differential fare structure table is determined based on trip distance and obtained service levels, though additional approaches are mentioned in Section 2.1. A pair with a long travel distance might consume more transit operation resources than a pair with a short travel distance. Subsequently, it is reasonable that a pair with a long distance trip should pay more for additional consumption of resources. In addition, feasible operational strategies are applied in this study, such as FRO, limited stop, short turn, and mixed strategies. Apparently, one of the important attributes characteristic of these strategies is that they provide different service levels for passengers. The greater the obtained service levels by pair , the more it should pay. Therefore, passengers with the same travel distance may need to pay different fares due to their obtained different service levels. Weight factor of fare per unit distance for each pair is used to show this pair obtained service levels. In this study, service level or qualityofservice is measured by travel time savings, especially because of skipped stop strategies used. The number of skipped stops between stops i and j by strategy is calculated as follows:
The weight factor of fare per unit distance for pair on strategy , , is expressed as follows:
Thus, an integrated differential fare for pair based on trip distance and obtained service levels is the product of trip distance, unit fare, and the average associated weight factor. It is formulated as follows:
(d) Crowding Time. The number of onboard passengers using strategy departing from stop , , which consists of passengers boarding at stop p and those passengers passing through stop p, is expressed as follows:
When the number of onboard passengers is beyond a threshold , it leads to crowding discomfort. The greater the rate between the number of onboard passengers and this threshold is, the more the crowding time for onboard passengers is [42]. Thus, the crowding time for pair using strategy , , is formulated as follows:
Consequently, average crowding time for pair , , is formulated as follows:
(e) Demand Function. The resulting number of passengers for pair along the bus line, , is formulated as where , , , and are the sensitivity parameters for the waiting time, invehicle time, fare, and invehicle crowding time, respectively. It is noteworthy that these parameters are not the actual elasticities in this linear function. The ratios, , , and , determine the value of waiting time, the value of invehicle time, and the value of crowding time, respectively [11, 26, 30].
2.4. Fare and Operational Strategies Model with Elastic Demand
An efficient transit system matches the complex relationship between supply and demand from the users’ and operators’ perspectives [1, 2]. Users pursue better service while the operator attempts to maximize operating profit. An attractive and efficient transit system will approach optimization of service from both perspectives.
(a) Consumer Surplus. Consumer surplus, , is the economic benefit to users. It is the price that users are willing to pay minus the price that they actually have to pay. Let be the inverse function of elastic demand function for an OD pair as follows. is a demand variable with the upper boundary .
Consumer surplus can be defined as follows [26, 43]:
Substituting (13) into (14), the final consumer surplus is formulated as follows:
(b) Operator Profit. Operator profit consists of operator revenue minus operator cost. The operator cost, , equals the fleet size (bus round trip time multiplied by the frequency) multiplied by the hourly operating cost per bus. It is expressed as follows:Operator revenue is considered only from passenger fares: Thus, operator profit function is formulated as
(c) Model Formulation. The objective function considered is to maximize the sum of consumer surplus and operator profit. Thus, the operational strategies problem using elastic demand is formulated as follows:
In this formulation, Constraint (20) guarantees that frequencies of strategies are nonnegative. Constraint (21) maintains that the number of onboard passengers departing from each stop is less than vehicle capacity. Constraint (22) shows whether or not strategies serve stops; it equals 1 if strategies serve stops and 0 otherwise. Constraint (23) ensures that all stops, including start and end stops of strategies and intermediate stops, can be visited by FRO and ST strategies. Constraint (24) makes sure that all strategies must visit their start and end stops. Constraint (25) guarantees that strategies involved with skipping stops do not skip more than three successive stops. Constraint (26) ensures that the expected demand for each pair is nonnegative and less than given potential demand.
The above problem is formulated as a mixed integer nonlinear programming (MINLP) problem. It can be handled by the outer approximation with both equality relaxation and augmented penalty (OA/ER/AP) algorithm of Viswanathan and Grossmann [44]. This algorithm starts by relaxing the integer requirement of decision variables and solving the relaxed nonlinear programming (NLP) problem. If an integer solution is not found, a sequence of iterations is solved consisting of NLP subproblems through fixing the integer variables and mixed integer linear program (MIP) master problems generated by augmented penalty function. This algorithm is readily available in the DICOPT solver of GAMS [45, 46]. Conopt and Cplex are adopted to solve the DICOPT subproblems of NLP and MIP.
3. Case Study
In this section, a real bus line in the city of Dalian, China, Line 26, is taken as a numerical example for applying the proposed models. Line 26, as shown in Figure 1 runs from Lingshui Passenger Transport Station to Wuyi Square, visiting 19 stops with 10.9 km in each direction, passing through a school zone (Dalian University of Technology, and the primary and high schools attached to Dalian University of Technology), a shopping center (Xi’an Road), and a business office zone (Software Park Service Center).
3.1. Data and Parameter Settings
In this case, a survey was conducted to collect route and travel demand data of Line 26 in the morning peak hour. The running time and distance between successive stops are shown in Table 2. Based on these collected demand data, the potential passenger demand OD matrix, for Line 26 in Table 3, was estimated. Accordingly, load profiles are also constructed, as shown in Figure 2.


(a) Passenger load profile in Direction 1
(b) Passenger load profile in Direction 2
Boarding time and alighting time per passenger are, respectively, s and s. The value of a parameter consisting of acceleration/deceleration as well as door opening and closing times at each stop is min. The layover time of a strategy at turning point is min. The demand elasticity parameters for wait time, invehicle time, and fare are , , , and , respectively [11, 30]. The vehicle operating cost is /vehiclehour. Capacity of a vehicle with 40 seats is passengers. The maximum acceptable number of passengers inside vehicle is passengers which generally is 125 to 150% of a bus’s seating capacity. Typically loads above 150% of a bus’s seating capacity subject standees and other passengers to unreasonable discomfort [47].
Figure 2 presents passenger load profiles for bus Line 26 in Dalian. The maxload segment of Line 26 is between stops 14 and 15 of direction 1 and between stops 15 and 16 the other direction; it is because both route segments are located in a building material market and the second biggest central business district (Xi’an Road), respectively. Based on the characteristics shown in Figure 2 it is possible to establish the sets of start and end stops of feasible operational strategies as is shown in Table 4.

3.2. Results
The optimal results obtained by the proposed methodology are shown in Table 5, taking a differential fare structure into account. It is evident that the optimal result is to operate a combination of all two strategies, FRO, limited stop, simultaneously with frequencies of 13.0 buses/h and 14.8 bus/h, respectively. Topological characteristics of these two resulting strategies are shown in Figure 3. The value of between two successive stops is in Figure 3. This optimal result is compared with the FRO strategy, without adding other strategies, to illustrate gains of $648.2/h for objective function value and consumer surplus of $771.8/h and attracting more than 124 passengers, despite a small decrease in operator profit of $123.6/h. However, increased objective function value indicates that a gain in consumer surplus compensates for this decreased operator profit. This promising result indicates that optimal strategies for single line bus operation can be applied to real life public transit systems, to attract more passengers, and to increase consumer surplus and total gains.

(a) FRO strategy
(b) Limited strategy
In Section 2, a fare structure table is defined. In order to clearly illustrate optimal fare changes along a bus line, this fare table is converted to a visual figure, as shown in Figure 4. The area in the upper triangle represents optimal fare for each pair in direction 1 of bus Line 26. Accordingly, the area in the lower triangle relates to optimal fares in the other direction. The pairs close to the diagonal have short travel distances, while those pairs far from the diagonal represent long distances. The optimal differential fare for each pair, with and without consideration for strategies, is presented in Figures 4(a) and 4(b), respectively. The optimal differential fare structure shown in Figure 4(b) is symmetrical and represents increasing values of fares from the diagonal toward two sides, while, in Figure 4(a), it is nonsymmetric and not always increasing. This results from the fact that the optimal fare for each pair is determined on the basis of bus service levels provided for this pair and its travel distance, in this study. In Figure 4(b), only FRO strategy is provided along a bus line. This implies that the same bus service level is offered for all pairs in both directions. Therefore, the optimal fares in this transit operation system are attained based on trip distances. That is, pairs and present the same fare due to the same trip distance, and the longer the travel distance for a pair, the higher the fare should be. However, in Figure 4(a), optimal strategies are applied in the transit operation system, which provide different service levels for pairs. According to Figure 4(a), some of the pairs which have achieved high service levels can be identified. For instance, the fare for pair (2, 16) is higher than that for pair (2, 17) in Figure 4(a), though its travel distance is shorter than for pair (2, 17). This is adequate explanation that pair (2, 16) needs to be offered a high service level.
(a) Optimal fare for each pair with strategies
(b) Optimal fare for each pair with FRO strategy and no other strategies added
Interestingly, Figure 5 indicates that OD pairs located in those areas with difference values between their two fares beyond $0 must be offered high quality service in the transit system which considers the use of strategies. Since compared with the unit fare from the model without consideration for strategies, the unit fare from the model considering these strategies is lower, as shown in Table 4. When pairs from the latter pay more than that from the former, it indicates that these pairs must gain better service than from the former.
Figure 6 compares the optimal results from two models with a flat fare and a differential fare structure, respectively. It is found that using a differential fare structure produces greater benefits for passengers and society, though there is a slight decrease in the profit for operators (see Figure 6(e)). These benefits are performed in terms of gains of about $326.9/h for objective function value, an increase of consumer surplus of about $1626.0/h, increased attraction of passengers by approximately 208 passengers, and, in particular, saving of an average fare per passenger by about $0.287, as shown in Figures 6(a)–6(d). It illustrates that using a differential fare structure in optimizing operational strategies shows more profitable, to a degree of considerable significance compared to using a flat fare.
(a) Objective function values ($/h)
(b) Consumer surpluses ($/h)
(c) Demands (passengers/h)
(d) Average optimal fares ($)
(e) Operator profits ($/h)
3.3. Sensitivity Analysis for a Differential Fare Structure
In this study, a differential fare function is developed based on trip distance and obtained service levels. This section will discuss the degree of how much passengers are willing to pay for service levels provided. Weight factor of fare per unit distance for each pair calculated by (7) is used to show their obtained service levels. A parameter is added in (7) as follows:
The second part on the right hand side of (27) represents the proportion of additional service fee in fare per unit distance that passengers from pair on strategy are willing to pay for gaining high service levels. Passengers do not need to pay an additional fee for high quality of service, if equals zero. The greater the value of is, the greater the proportion of additional service fee in fare per unit distance is. Note that the differential fare structure based on trip distance and obtained service levels will be based on travel distances if the value of the second part equals zero.
Figure 7 shows that optimal results change as the value of parameter related to additional service fee increases, in terms of objective function value, operator profit, consumer surplus, demand, and fare per unit distance. As shown in Figure 7(a), it is found that the greatest objective function value is about $20949.0/h when the value of equals 0. It gradually decreases with the value of increasing from 0 to 250. It is invariant when the value of is greater than 250. That is, as is shown by (27), the weight factor of fare per unit distance for each pair must increase with the increase of the value of parameter , thus reducing passenger demand and the objective function value; this is given the use of the limited stop strategy to serve the same set of stops using the service frequency. According to Figure 7(f), it is found that the frequency of the limited stop strategy decreases gradually while the frequency of the FRO strategy increases with the increase of the value of parameter from 0 to 250. This aims to reduce the adverse effects on the objective function value generated by increased weight factor of fare per unit distance. When the value of parameter increases to be greater than 250, the limited stop strategy makes a maximum adverse contribution to the objective function value which is less than that from the model using the FRO strategy. Therefore, the skipped strategy is abandoned, except the FRO strategy. By using only the FRO strategy, the value of weight factor of fare per unit distance equals 1 regardless of the increase of parameter , because of the second term of (27) being zero. That is, the differential fare is constructed only based on trip distances. The objective function value becomes a fixed constant in a FRO strategy problem. Accordingly, this also explains why other components of optimal results, as shown in Figures 7(b)–7(e), are constant when the value of is greater than 250. In addition, the value of objective function using this differential fare structure is less than that by about $20496.3/h using a flat fare, when the value of is greater than 10. This indicates a flat fare is more profitable than a differential fare for the application of strategies, if the additional service fee is high and more than a threshold.
(a) Objective function value ($/h)
(b) Operator profit ($/h)
(c) Consumer surplus ($/h)
(d) Demand (passengers/h)
(e) Fare per unit distance ($/km)
(f) Frequency of operational strategy (vehicles/h)
It is obvious that the effect of an additional service fee on a trip fare will become greater, as the value of rises from 0 to 250, according to (27). In order to lower the increasing adverse effects, of an additional service fee on the value of the objective function, the optimization model also attempts to reduce the optimal fare per unit distance to the extent possible in Figure 7(e); this is except for the method, mentioned above, that reduces frequency of the limited stop strategy and increases frequency of the FRO strategy in Figure 7(f). In Figure 7(b) the operator profit increases to a fixed constant as the value of increases. The demand and consumer surplus present an opposite trend (see Figures 7(c) and 7(d)) to that of operator profit. It is found that there is a drastic change of the operator profit, consumer surplus, demand, and fare per unit distance once the value of increases from 200 to 250. This is directly related to the use of operational strategies. For close to 250, the use of the combined limited stop and FRO strategies is transformed into only using the FRO strategy in Figure 7(f). In the former case, the value of has a great adverse effect on the objective function value, while there is no effect in the latter. Additionally, consumer surplus achieves the maximum value when the value of equals zero, while the operator profit is at its minimum level. This offers transport planners and decision makers a way to make a beneficial service from the perspectives of the society, passengers, and operators. Multiple operational strategies integrated with a distance based fare structure can bring the maximum benefits for the society and passengers, with a minimum profit for operators. The multiple operational strategies integrated with a distance based and service based fare structure can present a better tradeoff between passengers and operators. If decision makers want to obtain the maximum profit for operators, FRO strategies with a distance based fare structure ought to be used.
4. Conclusion
This work focuses on operational strategies and fare problems on a bus line. These operational strategies consist of FRO, limited stop, short turn, and mixed strategies, constructed on the basis of passenger demand. According to characteristics of this transit operation system, a fare table is proposed. This fare table can be presented for a flat fare and a differential fare structure, depending on fare policy objectives. Both types of fare structure are considered in this study, in order to seek a more suitable, profitable fare structure for the application of strategies. The differential fare structure proposed is constructed on the basis of trip distance and obtained service levels. Moreover, passenger demand is sensitive with respect to waiting time, invehicle time, fare, and crowding time in vehicles. This allows for investigating the effects of using operational strategies and two types of fare structures on passengers and the operators, as well as the transit system as a whole. The strategies and fare problems, therefore, are formulated as optimization models, with the objective of maximizing a sum of benefits of users and operators.
This model has been applied to a real life example in Dalian, China. It shows that using strategies can improve results by gains of $648.2/h for objective function value and consumer surplus gains of $771.8/h and attracting more than 124 passengers, despite a small loss in operating profit of $123.6/h, compared with applying FRO strategy exclusively, when a differential fare structure is applied for both optimization models. Clearly, passenger gains can compensate for this loss of operator profit.
In addition, in comparison with optimizing strategies considering a flat fare, a differential structure proves more profitable in terms of objective function value, passenger surplus, and demand, if the value of parameter related to additional service fee is less than 10. This illustrates that using optimal operational strategies integrated with a differential fare structure can further improve a transit operation system.
Moreover, sensitivity analysis is conducted to examine the effect of an additional service fee on optimal results. It is found that operational strategies using a flat fare provide greater objective function value than those using a differential fare, when the value of parameter related to additional service fee is between 10 and 250. Interestingly, when this value increases to the point of exceeding 250, skipped strategy, except for FRO, may be abandoned due to a high additional service fee required for providing high service levels. This indicates that FRO strategy performs in a more beneficial manner when the proportion of the additional service fee for a trip fare is more than a threshold.
Future research could extend the proposed methodology to include provision of information on bus arrival/departure and running, passengers transferring between vehicle trips associated with different operational strategies, and transport emissions considerations. In the modern, multimodel, urban transportation system, the application of an optimal combination of operational strategies integrated into a differential fare structure, may drive more travelers to leave their cars and turn to public transit services in deference to resulting high service levels of the transit service. The online information can be provided via smartphone and Internet and thus allows for reducing waiting time by timing arrival at stops. We also note that some passengers may transfer between vehicle trips associated with different operational strategies when transferring costs are less than their obtained travel savings. Moreover, operational strategies’ attributes can reduce transit emissions by reducing the number of “stops and goes,” though this is not explicitly discussed in this study. This would be inserted into the objective function.
Notations
:  The maximum number of passengers inside a vehicle with an accepted comfortable level (passengers/vehicle) 
:  Passenger alighting times for a bus using strategy at stop (h/vehicle) 
:  Passenger boarding times for a bus using strategy at stop (h/vehicle) 
:  Operating cost per bus hour ($/vehicle h) 
:  Average crowding time for pair (h) 
:  Consumer surplus for trips ($/h) 
:  Operator cost ($/h) 
:  Dwell time of a bus using strategy between stop and stop (h/vehicle) 
:  Trip distance for pair (km) 
:  Set of end stop of strategies 
:  Demand elasticity parameter for waiting time 
:  Demand elasticity parameter for invehicle time 
:  Demand elasticity parameter for fare 
:  Demand elasticity parameter for crowding time 
:  Bus fare for pair ($) 
:  Frequency of strategy (vehicles/h) 
:  Operator revenues ($/h) 
:  Potential demand for pair along a bus line (passengers/h) 
:  Layover time of strategy on a bus line at the turning point (h) 
:  A set of bus strategies serving a bus line, including full route operation (FRO), short turn (ST), limited stop (LS), and a mix of limited stop and short turn (MLS) 
:  Set of stops on a bus line 
:  Set of feasible stops on strategy 
:  The skipped number of stops that a bus using strategy travels from stop to stop 
:  The number of stops from stop to stop on a bus route 
:  Operator profit ($/h) 
:  Fare per unit distance ($/km) 
:  Capacity of vehicle (passengers/vehicle) 
:  An demand variable of the inverse function of elastic demand function with the upper boundary 
:  Running time of a bus using strategy between stop and stop (h/vehicle) 
:  Running time of strategy on a bus line (h) 
:  Running time for a bus from stop to stop (h/vehicle) 
:  Set of start stop of strategies 
:  Average invehicle time for pair (h) 
:  Waiting time for pair (h) 
:  Average boarding time per passenger (h/passenger) 
:  Average alighting time per passenger (h/passenger) 
:  Flat fare for all OD pairs ($) 
:  The resulting number of passengers for pair (passengers/h) 
:  The number of onboard passengers for a bus using strategy departing from stop (passengers/vehicle) 
:  The inverse function of elastic demand function for pair 
:  Crowding time for pair using strategy (h) 
:  The sum of consumer surplus and operator profit ($/h) 
:  Equals 1 if passengers of pair are boarding at stop and 0 otherwise 
:  Equals 1 if passengers of pair are alighting at stop and 0 otherwise 
:  Equals 1 if passengers of pair are passing through stop and 0 otherwise 
:  A constant consisting of deceleration/acceleration as well as doors opening and closing time at each stop (h/vehicle) 
:  A parameter depending on the distribution of bus arrivals at each stop; when bus arrival is Poissondistributed, will be equal to 1 
:  A parameter used to adjust the proportion of obtained additional service fee in fare per unit distance 
:  Weight factor of unit fare for pair using strategy , depending on pair obtained service levels 
:  Binary variable; taking the value of 1 if strategy serves stop and 0 otherwise. 
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The first author would like to acknowledge the support of the China Scholarship Council (CSC) for this study.
References
 A. Ceder, Public Transit Planning and Operation: Theory, Modeling and Practice, Elsevier, Oxford, UK, 2007.
 A. Ceder, Public Transit Planning and Operation: Modeling. Practice and Behavior, CRC Press, Boca Raton, Fla, USA, 2nd edition, 2016.
 P. G. Furth and F. B. Day, “Transit routing and scheduling strategies for heavy demand corridors,” Transportation Research Record, vol. 1101, pp. 23–26, 1985. View at: Google Scholar
 P. G. Furth, “Short turning on transit routes,” Transportation Research Record, vol. 1108, pp. 42–52, 1987. View at: Google Scholar
 A. Ceder, “Optimal design of transit shortturn trips,” Transportation Research Record, vol. 1221, pp. 8–22, 1989. View at: Google Scholar
 P. Delle Site and F. Filippi, “Service optimization for bus corridors with shortturn strategies and variable vehicle size,” Transportation Research Part A: Policy and Practice, vol. 32, no. 1, pp. 19–38, 1998. View at: Publisher Site  Google Scholar
 A. Tirachini, C. E. Cortés, and S. R. JaraDíaz, “Optimal design and benefits of a short turning strategy for a bus corridor,” Transportation, vol. 38, no. 1, pp. 169–189, 2011. View at: Publisher Site  Google Scholar
 C. E. Cortés, S. JaraDíaz, and A. Tirachini, “Integrating short turning and deadheading in the optimization of transit services,” Transportation Research Part A: Policy and Practice, vol. 45, no. 5, pp. 419–434, 2011. View at: Publisher Site  Google Scholar
 L. Fu, Q. Liu, and P. Calamai, “Realtime optimization model for dynamic scheduling of transit operations,” Transportation Research Record, no. 1857, pp. 48–55, 2003. View at: Google Scholar
 Z. Liu, Y. Yan, X. Qu, and Y. Zhang, “Bus stopskipping scheme with random travel time,” Transportation Research Part C: Emerging Technologies, vol. 35, pp. 46–56, 2013. View at: Publisher Site  Google Scholar
 J. Chen, Z. Liu, S. Zhu, and W. Wang, “Design of limitedstop bus service with capacity constraint and stochastic travel time,” Transportation Research Part E: Logistics and Transportation Review, vol. 83, pp. 1–15, 2015. View at: Publisher Site  Google Scholar
 C. Leiva, J. C. Muñoz, R. Giesen, and H. Larrain, “Design of limitedstop services for an urban bus corridor with capacity constraints,” Transportation Research Part B: Methodological, vol. 44, no. 10, pp. 1186–1201, 2010. View at: Publisher Site  Google Scholar
 P. R. Tétreault and A. M. ElGeneidy, “Estimating bus run times for new limitedstop service using archived AVL and APC data,” Transportation Research Part A: Policy and Practice, vol. 44, no. 6, pp. 390–402, 2010. View at: Publisher Site  Google Scholar
 A. M. ElGeneidy and J. SurprenantLegault, “Limitedstop bus service: an evaluation of an implementation strategy,” Public Transport, vol. 2, no. 4, pp. 291–306, 2010. View at: Publisher Site  Google Scholar
 V. Chiraphadhanakul and C. Barnhart, “Incremental bus service design: combining limitedstop and local bus services,” Public Transport, vol. 5, no. 12, pp. 53–78, 2013. View at: Publisher Site  Google Scholar
 Y. Y. Ulusoy, S. Chien, and C.H. Wei, “Optimal allstop, shortturn, and express transit services under heterogeneous demand,” Transportation Research Record, no. 2197, pp. 8–18, 2010. View at: Publisher Site  Google Scholar
 Y. Y. Ulusoy and S. I.J. Chien, “Optimal bus service patterns and frequencies considering transfer demand elasticity with genetic algorithm,” Transportation Planning and Technology, vol. 38, no. 4, pp. 409–424, 2015. View at: Publisher Site  Google Scholar
 A. Evans, “A theoretical comparison of competition with other economic regimes for bus services,” Journal of Transport Economics & Policy, vol. 21, no. 1, pp. 7–36, 1987. View at: Google Scholar
 H. Yang and H.J. Huang, “Analysis of the timevarying pricing of a bottleneck with elastic demand using optimal control theory,” Transportation Research Part B: Methodological, vol. 31, no. 6, pp. 425–440, 1997. View at: Publisher Site  Google Scholar
 W. H. K. Lam and J. Zhou, “Optimal fare structure for transit networks with elastic demand,” Transportation Research Record, no. 1733, pp. 8–14, 2000. View at: Google Scholar
 H. Yang and W. K. Kin, “Modeling bus service under competition and regulation,” Journal of Transportation Engineering, vol. 126, no. 5, pp. 419–425, 2000. View at: Publisher Site  Google Scholar
 G. Bellei, G. Gentile, and N. Papola, “Network pricing optimization in multiuser and multimodal context with elastic demand,” Transportation Research Part B: Methodological, vol. 36, no. 9, pp. 779–798, 2002. View at: Publisher Site  Google Scholar
 J.L. Yao, C.F. Hou, Z.X. Yang, and C. Fu, “Optimization study on road pricing with variable bus frequency,” Journal of Transportation Systems Engineering and Information Technology, vol. 10, no. 2, pp. 93–99, 2010. View at: Google Scholar
 F. Jørgensen, T. A. Mathisen, and B. Larsen, “Evaluating transport user benefits and social surplus in a transport market—the case of the Norwegian ferries,” Transport Policy, vol. 18, no. 1, pp. 76–84, 2011. View at: Publisher Site  Google Scholar
 M. J. Klier and K. Haase, “Urban public transit network optimization with flexible demand,” OR Spectrum. Quantitative Approaches in Management, vol. 37, no. 1, pp. 195–215, 2015. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. K. Chang and P. M. Schonfeld, “Multiple period optimization of bus transit systems,” Transportation Research Part B, vol. 25, no. 6, pp. 453–478, 1991. View at: Publisher Site  Google Scholar
 S. I.J. Chien and L. N. Spasovic, “Optimization of grid bus transit systems with elastic demand,” Journal of Advanced Transportation, vol. 36, no. 1, pp. 63–91, 2002. View at: Publisher Site  Google Scholar
 F.M. Tsai, S. I.J. Chien, and L. N. Spasovic, “Optimizing distancebased fares and headway of an intercity transportation system with elastic demand and trip length differentiation,” Transportation Research Record, no. 2089, pp. 101–109, 2008. View at: Publisher Site  Google Scholar
 I. Savage, “The dynamics of fare and frequency choice in urban transit,” Transportation Research Part A: Policy and Practice, vol. 44, no. 10, pp. 815–829, 2010. View at: Publisher Site  Google Scholar
 Z.C. Li, W. H. K. Lam, S. C. Wong, and A. Sumalee, “Design of a rail transit line for profit maximization in a linear transportation corridor,” Transportation Research Part E: Logistics and Transportation Review, vol. 48, no. 1, pp. 50–70, 2012. View at: Publisher Site  Google Scholar
 M. E. Kim and P. Schonfeld, “Maximizing net benefits for conventional and flexible bus services,” Transportation Research Part A: Policy and Practice, vol. 80, pp. 116–133, 2015. View at: Publisher Site  Google Scholar
 R. Voith, “Fares, service levels, and demographics: what determines commuter rail ridership in the long run?” Journal of Urban Economics, vol. 41, no. 2, pp. 176–197, 1997. View at: Publisher Site  Google Scholar
 R. A. Petitte, “Fare variable construction and rail transit ridership elasticities case study of the Washington, D.C., metrorail system,” Transportation Research Record, no. 1753, pp. 102–110, 2001. View at: Google Scholar
 N. Paulley, R. Balcombe, R. Mackett et al., “The demand for public transport: the effects of fares, quality of service, income and car ownership,” Transport Policy, vol. 13, no. 4, pp. 295–306, 2006. View at: Publisher Site  Google Scholar
 O. Cats, J. West, and J. Eliasson, “A dynamic stochastic model for evaluating congestion and crowding effects in transit systems,” Transportation Research Part B: Methodological, vol. 89, pp. 43–57, 2016. View at: Publisher Site  Google Scholar
 D. Fleishman, N. Shaw, M. Joshi, and R. Freeze, “Fare policies, structures, and technologies,” TCRP Report 10, Transportation Research Board, Washington, DC, USA, 1996. View at: Google Scholar
 S. I.J. Y. Chien and C. F. M. Tsai, “Optimization of fare structure and service frequency for maximum profitability of transit systems,” Transportation Planning and Technology, vol. 30, no. 5, pp. 477–500, 2007. View at: Publisher Site  Google Scholar
 Z.C. Li, W. H. Lam, and S. C. Wong, “The optimal transit fare structure under different market regimes with uncertainty in the network,” Networks and Spatial Economics, vol. 9, no. 2, pp. 191–216, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. Gkritza, M. G. Karlaftis, and F. L. Mannering, “Estimating multimodal transit ridership with a varying fare structure,” Transportation Research Part A: Policy and Practice, vol. 45, no. 2, pp. 148–160, 2011. View at: Publisher Site  Google Scholar
 F.M. Tsai, S. Chien, and C.H. Wei, “Joint optimization of temporal headway and differential fare for transit systems considering heterogeneous demand elasticity,” Journal of Transportation Engineering, vol. 139, no. 1, pp. 30–39, 2013. View at: Publisher Site  Google Scholar
 D. Yook and K. Heaslip, “Determining appropriate fare levels for distancebased fare structure: considering users' behaviors in a timeexpanded network,” Transportation Research Record, vol. 2415, pp. 127–135, 2014. View at: Publisher Site  Google Scholar
 Z. Li and D. A. Hensher, “Crowding and public transport: a review of willingness to pay evidence and its relevance in project appraisal,” Transport Policy, vol. 18, no. 6, pp. 880–887, 2011. View at: Publisher Site  Google Scholar
 G. Kocur and C. Hendrickson, “Design of local bus service with demand equilibration,” Transportation Science, vol. 16, no. 2, pp. 149–170, 1982. View at: Publisher Site  Google Scholar
 J. Viswanathan and I. E. Grossmann, “A combined penalty function and outerapproximation method for MINLP optimization,” Computers and Chemical Engineering, vol. 14, no. 7, pp. 769–782, 1990. View at: Publisher Site  Google Scholar
 A. Brooke, D. Kendrick, A. Meeraus, and R. Raman, GAMS: A User's Guide, GAMS Development Corporation, Washington, DC, USA, 2005.
 I. E. Grossmann, J. Viswanathan, and A. Vecchiette, DICOPT User's Manual, GAMS Development Corporation, Washington, DC, USA, 2008.
 Transportation Research Board, Transit Capacity and Quality of Service Manual, Transit Cooperative Research Program, Report 100, Transportation Research Board, Washington, DC, USA, 2nd edition, 2003.
Copyright
Copyright © 2017 Chunyan Tang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.