Research Article  Open Access
Emre Tokgöz, Iddrisu Awudu, Saravanan Kuppusamy, "Advancing Transportation Routing Decisions Using Riemannian Manifold Surfaces", Journal of Advanced Transportation, vol. 2020, Article ID 2098495, 17 pages, 2020. https://doi.org/10.1155/2020/2098495
Advancing Transportation Routing Decisions Using Riemannian Manifold Surfaces
Abstract
We consider several realworld driving factors such as the time spent at traffic signs (e.g., yield signs and stop signs), speed limits, and the topology of the surface to develop realistic and accurate routing solutions. Though these factors increase the complexity of modeling, they provide the flexibility to evaluate the routing solutions from different perspectives: cost, distance, and time, to name a few. First, we develop a set of algorithms based on the Riemannian manifold surface (RMS) to factor in the Earth’s curvature to calculate distances. Second, we present a multiobjective, nonlinear, mixedinteger model (MINLP) that minimizes the distance traveled, time traveled, traveling costs, and time spent on traffic signs to design and evaluate the routes where the waiting times associated with traffic lights, stop signs, and yield signs are stochastic. Finally, we apply MINLP and RMSbased algorithms to a set of reallife and short and longdistance transportation problems and analyze the results from computational experiments and discrete event simulations. We show that our approaches are on par with the stateoftheart application, Google Maps, and yield realistic routing solutions that generate significant cost savings.
1. Introduction
Transportation is one of the essential activities in a supply chain. For example, goods manufactured at plant locations must be distributed to customer locations through warehouses, distribution centers, or crossdocking facilities. Typically, a fleet of trucks (owned or outsourced) performs the task of distributing goods. In the view of fuel prices, the complexity of the supply chain, and the uncertainties in the market place, the trucking operations must be efficient in order to keep costs under control. For instance, in 2019, US companies spent 668 billion dollars on trucking costs alone [1]. As a result, businesses are continually looking for a new approach or technique to reduce transportation costs. In product distribution, the sequence (or route) in which goods are delivered by a truck determines the distance traveled by the truck, which, in turn, determines the cost of transportation. The distance between locations is often used as a surrogate measure for the transportation cost. It is therefore vital to understand how distances are calculated.
The assumptions related to distance calculations play a key role in determining the accuracy of the measurement. One of the assumptions is the surface on which the distance is measured. Often, the Euclidean surface is used because it is less complex to calculate distances [2]. However, the Euclidean surface generates inaccuracies in distance calculations because it does not take into account the curvature of the Earth. As a result, the use of the Euclidean surface has implications such as suboptimal routing decisions, resulting in increased transportation costs. Motivated by the above, we develop an approach based on the Riemannian manifold surface (RMS) that factors in the curvature of the Earth when calculating distances. For instance, the curvature of a Riemannian manifold surface indicates how much Earth’s surface is curved in a local neighborhood. As the surface changes, the curvature changes. We use this curvature change to calculate the shortest path distances. Subsequently, we demonstrate that the RMSbased approach improves the distance calculation, yields realistic distances, and generates significant transportation cost savings.
In this work, to factor in the curvature of the Earth, we use the generalized notion of a “straight line” to “curved spaces,” i.e., a geodesic. The calculation of geodesic distances involves the use of local topological regions on Earth’s surface. For example, while sprinting on multilane running tracks, a runner on the inner lane starts at the back, while a runner on the outer lane starts at the front. The starting points of the runners are determined by factoring in the geodesic distances, which are two dimensional as the running tracks are flat. However, driving a vehicle involves threedimensional geodesic calculations that we use in this research.
Although companies use distance as a critical measure to evaluate the routing decisions, end consumers use various indicators (criteria) to determine routes for work and leisure travel. For instance, travel time, travel costs, travel distance, and waiting time spent at traffic signs (e.g., stop signs, yield signs, and traffic lights) are often used. Nevertheless, existing global positioning and global information systems do not provide, for example, the impact of traffic signals on travel cost and time and the flexibility to evaluate routing decisions using various criteria. Moreover, these systems do not consider local shortest geodesic distances. Inspired by the above, we develop a multiobjective, nonlinear, mixedinteger model (MINLP) where the end users can evaluate the routes with multiple criteria based on travel distance, cost, time, and traffic controls such as stop signs, traffic lights, and yield signs. Next, we illustrate the RMSbased approach and the MINLP model, and later, we perform simulation experiments using reallife instances.
1.1. An Illustration of Riemannian Manifold Surfaces and Multiple Criteria Route Evaluation
In this section, we illustrate two key concepts that are used in our methodology: (i) Riemannian manifold surface (RMS) based approach to calculate the distance; (ii) multiple criteria to evaluate the routes. For illustration purposes, we choose two points, Origin, , and Destination, . Let and be Hamden and New Haven, two cities in the state of Connecticut, for the ease of illustration as the authors are familiar with the geography. There are three routes between Hamden and New Haven: (A) via I91 South—using highway; (B) via Whitney Avenue—using local roads; (C) via Ridge Road—using local roads.
To illustrate the concept of RMS, we focus on Route A. Figure 1 is a schematic representation of Route A. In this figure, the solid lines represent the boundaries of the road, and the road is divided into two lanes: Lane 1 and Lane 2. For instance, a vehicle traveling from the origin to the destination can either start from Lane 1 or 2. If the vehicle starts from Lane 1, switches to Lane 2 when the road curves, and stays in Lane 2 until it reaches the destination, the path traced by the vehicle is referred as the innermost lane. Conversely, if the vehicle starts at Lane 2, switches to Lane 1 when the road curves, and stays in Lane 1 until it reaches the destination, the path traced by the vehicle is referred as the outermost lane. The dashed lines in Figure 1 represent the inner and outermost lanes of the road. The distance traveled by the vehicle using inner and outermost lanes is 10.84 miles and 11.38 miles, respectively. Note that both the inner and outer lanes are created based on the shortest and the longest geodesic paths, respectively, that takes the Earth’s curvature into account which is a central theme of the RMSbased approach (the details are provided later in the paper). The approach based on RMS yields savings of 11.38 − 10.84 = 0.53 miles for Route A. Similarly, savings on Routes B and C are identified using RMS, and Table 1 summarizes the results. Note that when two locations are only a few miles apart, the impact of the Earth’s curvature on the distance calculated is very minimal. This can be seen in 1, as the difference between the outermost and innermost lanes is less than a mile for all routes. Later, in the numerical section, we show that the difference between the outermost and innermost lanes (i.e., distance savings) can be significant for longdistance routes. In such cases, it is feasible for a traveler to devise a driving plan to realize the savings in the distance.

Next, using the above realworld instance, we illustrate how a decision maker can use multiple criteria to evaluate routes A, B, and C. We consider the following criteria: distance traveled, waiting time at traffic signs, travel cost, the number of traffic signs, and travel time. Table 2 provides a summary of the evaluation of the routes, and the associated details are in Section 6. Note that for each criterion, the best values are marked in bold.

For instance, if the distance traveled or the traveling cost is used as a sole criterion, then Route B is the best. However, Route B is the worst when the time traveled is considered. In this case, both A and C are better than B. When the waiting time at traffic signs is considered, Route C is the best. For instance, Route A has 1 stop sign, 1 yield sign, and 16 traffic lights which totals to 18 traffic signs. Route B and C has a total of 29 and 21 traffic signs, respectively. If the number of traffic signs is used as a criterion, Route A is the best. The above illustration highlights the importance of using multiple criteria in route evaluation as the route selection varies with the criteria considerably.
In summary, we develop a set of algorithms that are based on Riemannian manifold surfaces to address the need to calculate the distance accurately. In addition, we formulate a mathematical programming model that considers the realworld driving factors such as the time spent at traffic signs (e.g., yield signs and stop signs) and speed limits to evaluate routing decisions from multiple perspectives: cost, distance, and time, to name a few. Furthermore, we apply our approaches to reallife transportation problems and analyze the results from computational experiments.
The rest of this paper is organized as follows: Section 2 summarizes the relevant literature and the research contributions of this paper. Sections 3 and 4 offer details on the RMSbased approach to calculate the distance. Section 5 provides details on mathematical programming formulations that use multiple criteria to evaluate the routes. Section 6 summarizes computational experiments and simulation results. Section 7 concludes with a general discussion of managerial insights and avenues for future research.
2. Literature Review
In this section, we first present the literature related to the distance calculation since one of the main contributions of our research is in developing a new methodology to calculate distances. Next, we look into relevant literature where multiple criteria are used to evaluate routing solutions. Finally, we present our contributions in relation to the past work.
2.1. Distance Calculation Methodologies
In practice, there have been several software applications developed to calculate distances using travel information such as oncoming traffic situation, resting areas, and time to destination. However, the visibility into mathematical methodologies that are used in such applications is very minimal. Therefore, it is a challenging task to subject the methodologies used by software to rigorous academic discussion. On the contrary, in the academia, research related to distance calculation methodologies is limited. In the papers where the distance calculation is addressed, typically, rectangular or Euclidean distance functions are used. In their seminal work, Love and Morris [3] provide mathematical proofs for calculating transportation distances using convex functions that extend beyond the rectangular or Euclidean distance calculations. Brimberg and Love [4] consider a new distance measure, the weighted onetwo norm, which is a positive linear combination of the preceding norms as compared to the Euclidean or rectangular distance calculations. Fernández et al. [5] compare the weighted norm, which has been proven to be an accurate distance calculation function, to the norm which uses different parameter estimates for the distance prediction. This methodology is extended by Brimberg and Love [4] and Üster and Love [6]. Other papers such as Von Hohenbalken and West [7]; Miyagawa [8]; Alpaydin et al. [9]; Griffith et al. [10]; Herica and Servio [11]; Havelock et al. [12] have all used different methodologies to calculate distances, from the Manhattan and Euclidean distances to Lobachevskian and norm methodologies.
2.2. Multiple Objective Routing Solutions
An extensive amount of research has been done in formulating methodologies to aid the decision maker with optimal routing decisions; see, Toth and Vigo [13] for a review on vehicle routing problems (VRP). We review only the research that considers multiple criteria/objectives to evaluate the routing solutions. For instance, Bowerman et al. [14] develop an approach to model the urban school bus routing problem that considers multiple criteria such as route length, number of routes, load balancing, length balancing, and student walking distance. Zografos and Androutsopoulos [15] identify routing solutions that minimize both transportation time and risk as a part of a decision support system to transport hazardous materials and address emergency situations. Demir et al. [16] consider a variant of the pollution routing problem that minimizes both the fuel consumption and driving time. Kovacs et al. [17] extend the VRP from a service perspective by considering several objectives related to improving driver consistency, arrival time consistency, and cost minimization. Kumar et al. [18] consider a production and pollution routing problem with time windows that minimizes the total operational cost and emissions. See Jozefowiez et al. [19] for an extensive review on research related to multiobjective optimization in routing problems. Though our research uses criteria that are related to the total cost and the distance, it also considers the number of traffic signs and the time spent on these signs as criteria that have not been considered in the past research.
As mentioned above, one of our primary goals is to highlight the lack of research in using criteria related to traffic signs and to stress the significance of using multiple criteria when making routing decisions. Therefore, we restrict our attention to modeling multiple criteria and demonstrating the model’s usefulness by applying our techniques to a few reallife instances. We note that our research does not involve designing heuristics for the routing problem. As future research, our work can be extended by employing approaches such as genetic algorithmbased techniques (e.g., [20, 21]), greedy heuristics with variable neighborhood search (e.g., [22]), learningbased heuristics (e.g., [23]), and fuzzy programming (e.g., [24]), to name a few, to solve largescale, multiobjective, routing problems.
2.3. Research Contributions of This Paper
Our work is closely related to the study by Alpaydin et al. [9], Griffith et al. [10], and Havelock et al. [12]. Alpaydin et al. [9] propose a nonparametric approach using neural networks for estimating actual distances compared to the known distance estimators which are parameterized functions of the coordinates of the points. Griffith et al. [10] explore the nonEuclidean nature of network spaces, including Lobachevskian space (via the Poncaŕe disc), and compare the Euclidean, Manhattan, Minkowskian, and Lobachevskian distance calculation predicating parameters. The authors conclude that the skewed transportation networks in largescale regional landscapes are characterized by a Lobachevskian geometry, whereas most smaller scale landscapes appear to be best characterized by a Minkowskian space whose parameters are between those of a Manhattan and Euclidean space. Havelock et al. [12] propose the use of learningbased strategies to solve the problem of distance calculation by using a stochastic method called the adaptive tertiary search (ATS) strategy. This methodology works by utilizing the information provided in the coordinates of the nodes to calculate the actual distances.
While we build on the previous work, our research differs in the following ways:(1)We develop a novel methodology based on the Riemannian manifold surface (RMS), a concept that uses the road network topology and factors Earth’s curvature into account. As part of our methodology, we propose a set of algorithms where we adapt and modify the linked chain methodology (LCM) developed by Tokgöz et al. [25] to calculate distances.(2)We formulate a multiobjective, nonlinear, mixedinteger model (MINLP) to design and evaluate the routes where distance traveled, time traveled, traveling costs, and time spent on traffic signs are minimized. In this model, the waiting times associated with traffic lights, stop signs, and yield signs are stochastic.(3)We apply our approaches to a set of reallife and short and longdistance transportation problems and analyze the results from computational experiments and discrete event simulations. We show that our approach yields realistic solutions and generates significant cost savings compared to using the stateoftheart applications (e.g., Google Maps).
3. Riemannian Manifold Surfaces and Linked Chain Method
Carl Friedrich Gauss proved Theorema Egregium (meaning remarkable theorem in Latin) in 1828 by establishing an important property of surfaces. The theorem states that the Gaussian curvature of a surface can be determined entirely by measuring distances along the paths on the surface, and the curvature does not depend on how the surface is embedded in a threedimensional space. Later, Bernhard Riemann extended Gaussian theory to higherdimensional spaces that allows distances and angles to be measured and the notion of curvature to be defined. The Riemannian approach is an intrinsic way to the manifold (a topological space that is locally Euclidean) and does not depend upon how the manifold is embedded in higherdimensional spaces. A Riemann manifold is a metric space on a smooth manifold with a smooth section of the positivedefinite quadratic forms on the tangent bundle when the length of a continuously differentiable curve is defined by the following equation:
The definition of length in (1) allows every connected Riemann manifold to become a metric space on which the distance between and on is defined as the infimum of when is a continuously differentiable curve joining and . Later, we use (1) to calculate the geodesic length which is instrumental in calculating route distances. One of the most popular use of the Riemannian manifold is in the development of general theory of relativity by Albert Einstein; in particular, his equations for gravitation are constraints on the curvature of spacetime.
The curvature of a Riemannian manifold surface (RMS) indicates how much the surface is curved in a local neighborhood. For instance, the Riemann manifold is a plane in the Euclidean space when the curvature is zero; however, when the curvature changes from zero to one, the surface changes from a plane to a 3dimensional sphere. This change in the curvature is important when calculating shortest path distances on Earth’s surface. This surface is a natural manifold, and the pathways on this manifold are natural geodesics.
The research on RMSbased approaches in transportation is almost nonexistent except by Tokgöz et al. [25]. For instance, Tokgöz et al. [25] solved a location routing problem (LRP) on RMS and showed that their solution methodology generated results that are more realistic than using Euclidean surfaces since the Earth’s surface has changing local RMS curvatures. The authors developed an approach called the linked chain method (LCM) to calculate geodesic distances. We adapt and modify the LCM approach to our problem. Typically, in the LCM approach, given an origin () and destination (), circles are generated to cover the focus area that contains the OD pair. Unlike Tokgöz et al. [25], these circles are not randomly generated, and in addition, the circles represent the underlying road structure in our case. For instance, using the LCM, we generate inner and outermost lanes for highways that provide a range of traveling distance which has not been addressed in the past research. Subsequently, the circles that cover geodesics are identified, and the length of the geodesics is calculated. The overlapping circles that contain geodesics are identified, and they are linked (hence, this approach is called the “linked chain”) to create a path that connects and . Next, we provide details on the linked chain method.
4. Algorithmic Solution for Shortest Path Geodesic Distance Calculation
In this section, we introduce four algorithms to determine the shortest path between two points: origin, , and destination, , that can be used for routing vehicles. In addition, there are two theorems and a corollary stated and proven by using isomorphism and LCMrelated routing results. The set of notations used in the development of algorithms is given in Table 3.

We assume geodesics are determined based on the pathway intersections and illustrate this assumption by using a Tintersection. For instance, a Tintersection, see Figure 2, has three arms, and each arm could be a geodesic. Note that geodesics can be curved lines, but for illustration purposes, geodesics are shown as straight lines, see Figure 2. We assume that there exist nonoverlapping circles containing all the geodesics; these circles form the topology on the manifold. As a result, each circle can have only one geodesic. For instance, in Figure 2, three geodesics (three arms) are covered by three different circles. The radius of these circles is different due to varying lengths of the roads that form the Tintersection. Moreover, if a path uses the Tintersection, two circles of varying radii that cover two geodesics will be a part of this path. For instance, geodesics 1 and 2 could be one path and geodesics 1 and 3 can be another. Our first objective is to minimize the total distance traveled which is calculated using (1) by replacing and with and ; i.e.,
4.1. Subroutine 1
Algorithm 1 shows the pseudocode for the first subroutine. First, a surrounding area within of and is selected. Next, is set to 0.1, and circles with radius are drawn so that the entire surrounding area is covered with circles. Some of these circles may have geodesics in them and some may not. Once geodesics are identified for a radius , i.e., ’s, they are added to a set . This process of identifying geodesics is repeated by incrementing , i.e., where is set to 0.001. Note that a geodesic that is identified by multiple radii will be added only once to . As a result, geodesics in are not duplicated. Each one of these geodesics could have a minimum and a maximum length based on its curvature, and these lengths are calculated using Algorithm 4.

Figure 3 illustrates the steps involved in Algorithm 1. For illustration purposes, a focus area that is similar to a busy neighborhood in Manhattan, New York, is selected. In this illustration, the focus area is represented by a rectangle that has 8 intersections: four on the top row and four on the bottom row, see Figure 3(a). First, the entire is covered with the circles of initial radius . For the purposes of clarity, the fully covered in Figure 3(a) is not shown; therefore, some of the areas in the figure are not covered. Next, the circles with geodesics are identified. In Figure 3(a), four geodesics are identified (on the top row) by four circles that are indicated by the dashed lines. We assign circles , , , and to geodesics , , , and , respectively. Finally, the above geodesics are added to set . The above steps are repeated with a new radius, . If geodesics are identified with the circles of new radii, they will be added to , otherwise, not.
(a)
(b)
(c)
Figure 3(b) illustrates the geodesics that are identified when the radius is equal to . In this figure, the circles that identify geodesics are represented by the dashed lines. We assign circles , , , and to geodesics , , , and , respectively. We add the above geodesics to set . Figure 3(c) is an illustration of the focus area when all geodesics are identified. In this figure, 22 geodesics are identified, which are numbered from 1 to 22 by 22 different circles which are in .
4.2. Subroutine 2
Algorithm 2 shows the pseudocode for the second subroutine. We use this subroutine to calculate the length of geodesics, , that are in . Initially, two empty sets are created, and . First, in set , we select a radius and examine the geodesics, , identified with that radius. Second, the length, , of every geodesic is calculated using the following function:
Third, we define to be a vector with three elements: the geodesic , the length of the geodesic , and the circle that covers the geodesic. Subsequently, we add the geodesic to set and delete the same geodesic from set . Finally, is added to set . The steps explained above are repeated until the distances for all the geodesics are calculated, i.e., until is empty. As an output from this subroutine, we identify a set that consists of vectors that contain the geodesics, length of geodesics, and unique circles that contain the geodesics.
4.3. Subroutine 3
In topology, geodesics are connected and form a path if circles that contain them are overlapping. This subroutine searches for overlapping circles to determine the paths within the focus area, . Algorithm 3 shows the pseudocode for the subroutine. Initially, an empty set, , is created to store all possible paths, i.e., linked list vectors of geodesics. Then, two elements and of set are chosen such that . If circles and that contain and are overlapping, the link is added to set ; otherwise, a different pair of is chosen. This process is repeated until all pairwise comparisons in are exhausted. As a final output, a set that consists of all possible paths is generated.

4.4. Subroutine 4
Algorithm 4 shows the pseudocode for the fourth subroutine. First, is examined to identify the circles that contain origin, , and destination, . Once the circles are identified, they are split into nonoverlapping circles and geodesics. This is done to improve the accuracy of calculation of the distance between and . For instance, if is the circle that contains and is the geodesic in , then is split into and and is split into and . Subsequently, the corresponding vectors, and , are added to . Similarly, , the circle that contains , is split into and and is split into and and the corresponding vectors, and , are added to . We assume number of pathways exist between and , and the set , a subset of , contains all those pathways. The length, , of all pathways is calculated to determine all possible distances from to . Finally, the path with the least distance, , is selected to be the shortest path between and .

Note that steps 14–16 in Algorithm 4 are used to calculate minimum and maximum geodesic distances for all from intersection to intersection which are then used to calculate the minimum and maximum length route between and . Algorithms 1–4 yield the following set of theorems.
Theorem 1. and are isomorphic.
Proof. There exists a set of unique circles, , that covers the set of geodesics within the given bounded region. One circle can cover only one geodesic, and one geodesic can be covered by only one circle. It is evident from Algorithm 2 that there exists a unique circle per geodesic; therefore, there exists a map that is onetoone. To show the converse, assume there are two circles in that correspond to the same geodesic. In this case, two circles must have the same overlapping intersections which would conflict with Algorithm 1 implementation. Therefore, is an onto map; there exists a 11 and onto map that projects to . The inverse map maps each in to in .
Theorem 2. LCM determines the distance minimizing shortest path geodesic between two points chosen on a compact pathconnected manifold region on Earth’s surface.
Proof. A compact pathconnected manifold region on Earth’s surface indicates the existence of a path that allows the existence of geodesics between the two points chosen in a compact region. Algorithm 2 is used to determine all possible pathways between the two points, while Algorithm 4 is used to determine the distance minimizing shortest path on the compact pathconnected predetermined region. The set in Algorithm 4 is the set of all pathways on the manifold’s surface.
Corollary 1. LCM applied on Earth’s surface does not always yield the shortest path geodesics between the two points.
Proof. Assume there exists a geodesic that yields the shortest path between the two chosen points on Earth’s surface. Such a geodesic would require to identify predetermined pathways with finite lengths and nonoverlapping loops; however, the looping nature of the manifold region could result in infinite distances. For instance, consider a region that only contains ocean. In this case, the chosen manifold’s curvature is always 1; therefore, there exist infinitely many geodesics on this surface due to the looping nature on a spherical surface. Such a compact region of Earth’s surface can generate infinite geodesics and infinite pathways. In this case, Algorithms 1–4 are not applicable.
5. A Multiobjective, MixedInteger, Nonlinear Mathematical Model
In this section, in addition to (2), we introduce four objective functions that minimizes (i) total time spent on traffic signs, (ii) total travel cost, (iii) total number of traffic signs, and (iv) total travel time. As the additional objective functions involve traffic signs such as stop signs, traffic lights, yield signs, and the speed limits, we introduce the assumptions and parameters that are associated with them. In some cases, objective functions such as travel distance and time might be correlated, but having them as individual functions provide the flexibility for the decision maker. For instance, in the densely populated areas (e.g., northeastern parts of US), mostly travel time and distance are not correlated. Note that we do not model road congestion in our work which could be an avenue for future research; however, the impact of the network congestion is captured by the objective functions such as travel time and cost. The set of notations used in the mathematical model is given in Table 4.

5.1. Assumptions and Parameters
(i)Stop Signs. There are four types of stop signs used for mathematical formulation: 1way, 2way, 3way, and 4way. Let be the number of stop signs for the stop on route for all , where is the total number of routes that connect to . The expected duration of a stop sign cannot be zero for lawful reasons; therefore, we assume that the wait time of a vehicle on the stop sign, , a stochastic parameter, varies between waiting times and . The set of all stop signs is where represents the set of stop signs on route for all and . For instance, a specific choice of a vehicle’s wait time on a stop sign can vary between waiting times and seconds.(ii)Traffic Lights. As one may have to wait at a traffic light longer than another (e.g., busy vs. nonbusy intersections), the maximum time spent can vary among traffic lights. Hence, it is reasonable to use the maximum time spent at traffic lights to categorize them. We categorize the traffic lights based on their maximum durations. Let the set represent the categoryi traffic lights based on and the maximum amount of time that a vehicle can spend waiting on redlight duration on route such that . For instance, represents the category of traffic lights that has a maximum duration of 25 seconds on the route. The set of traffic lights for different paths from to takes the form for all and where is the number of traffic light categories. The number of traffic lights on the pathway in category is represented by , and the total number of traffic lights on each path from to is for all . For the example given above, if the set has 5 elements, then indicates that there are 5 traffic lights with 25 seconds of maximum wait time on route 100. Let represent the stochastic parameter that varies between 0 and for waiting at the categoryi traffic light on route .(iii)Yield Signs. There are two types of traffic yields considered for a vehicle to join traffic: the yield signs to enter the highway traffic and the city traffic. Let represent the number of yield signs to enter the highway and represent the number of yield signs to enter the city traffic on route for all . Let for each , and let be the stochastic parameter representing the wait time at for , 2. For instance, the average time it takes for a vehicle to enter a highway when there is a yield sign can be assumed to vary between 4 and 8 seconds, and the time it takes to enter the city traffic when there is a yield can be assumed to vary between 5 and 10 seconds on average.(iv)Speed Limit. The distance calculations depend on speed limit, geodesic distance, and traffic sign durations. The LCM (linked chain method) circles are used to determine the geodesic distances from one intersection to the other on each route. Each LCM circle contains the corresponding road’s speed limit. We assume the set of speed limits to be where there are different number of speed limits that could occur from to . If there are more than one in a LCM circle, then we consider an average of speed limits for timerelated calculations. For instance, using geodesics from Algorithm 2 (subscript is modified as ) and the corresponding length , the distance traveled on the LCM circle can be calculated as follows: , where is the speed limit value on the road covered by . Later, we use this formula in the MIP formulations to determine the travel time from to .
5.2. Objective Functions and Constraints
In this section, we introduce five objective functions and develop three different formulations. Let represent the geodesic on route that connects and and is a set of all possible geodesics between and with the index set . Let and be the fixed permile fuel cost and permile operational cost, respectively. The constrained, mixedinteger, nonlinear program formulation for the path problem to travel from to that includes can be formulated as follows:subject to
Recall that algorithms from Section 4 are used to determine the shortest geodesic path, and the corresponding objective function minimizes the travel distance. We reintroduce in (10) for clarity. In this formulation, is an integer variable that denotes the route. In addition, in (10) is a nonlinear, continuous variable that is used to calculate the geodesic distance on nonlinear curves. The waiting times associated with traffic lights, stop signs, and yield signs are stochastic parameters. Since the number of traffic signs on each route is wellknown and do not vary, they are constant values. Equation (11) minimizes the total waiting time on traffic signals while traveling from to . Equation (12) minimizes the cost incurred (in dollars) to travel from to . Equation (13) minimizes the total number of traffic signs, and equation (14) minimizes the total travel time. Equations (15) and (16) determine the waiting times at the city and highway yield signs, respectively. Equation (17) determines the waiting time at stop signs, and equation (18) provides the waiting times at the traffic lights. The set element specifications are given in (19)–(21). Equations (10)–(14) determine the minimum values of for all and .
5.2.1. Weighted Sum Formulation (WSF)
The weighted sum formulation (WSF) allows the decision maker to choose a route based on the weights assigned to the objective functions through . Function introduced below normalizes the objective function values – for each route, and it is used as a part of minimizing the weighted sum. This option of formulation requires weighted sum formulation of equations (15)–(19) by introducing the weight function with the corresponding objective function:where is the finite index set and are the normalized values of the objective functions calculated for each route by using the formula:for each and each such that . The values of are determined for all and that forms a set; this set of values is used for determining the minimum value of . The weights are known constants that are assigned by the decision maker, and the sum of the weights must be equal to 100%.
5.2.2. Best Vector Set Solution (BV)
In this section, we propose a best vector set solution approach where the vector set consists of the best solutions determined by minimizing each . Since the optimal route can vary with , the corresponding vector set solution can be different. In this case,where , , b is the total number of the routes, and the best vector set which is a subset of all solutions. For instance, assume b = 4. Let route 1 (i.e., u = 1) be optimal for and , route 2 (u = 2) be optimal for and , route 3 (u = 3) be optimal for , , and , and route 4 (u = 4) be not optimal for any . The best vector (BV) solution is given as follows:
Note in BV that there is no minimum objective function value for route 4 as it is not optimal for any . The routes in the BV solution may or may not be displayed in the corresponding vectors listed in the set. Based on the end user’s preference, a specific route could be picked from the vector set. For instance, if the user is interested in saving time and money, a route (s) that minimizes both travel time and travel cost can be chosen from the vector set. Later, we illustrate the BV solution for two reallife routing problems in Section 6.
5.3. An Algorithmic Solution
In this section, we develop an algorithmic solution for the mathematical models presented in Section 5.2. Recall that equations (10), (12), and (14) are nonlinear since the distance calculations of geodesics depend on nonlinear curves. The numbers of traffic controls such as stop signs, yield signs, and traffic lights are known constants for each route. Therefore, the objective function in equation (13) is a constant for a route. Algorithm 5 uses the outputs of Algorithms 1–4 from Section 4.

In this algorithm, constraints (15)–(18) introduced in Section 5.2 are randomized by using the minimum and maximum values (that are constant values and assumed to be known) for traffic signs. These constraint values are then used for calculating . Values of are the inputs to the algorithm for calculating the weighted sum value . The values attained from are then used for calculating the best vector set solution ; the decision maker can choose a route from a variety of routes by using that depends on the minimum values attained for . A reallife application of Algorithm 5 with specific entries for minimum and maximum values is given in Section 6.
For instance, the set which is an output of Algorithm 2 is a critical input to Algorithm 5. The set consists of vectors which has the following information: (i) the geodesic , (ii) the length of the geodesic , and (iii) the LCM circle that has . The stochastic parameters are calculated by using a uniform distribution.
Note that running Algorithm 5 for sufficient number of times and averaging the outcomes will yield the desired objective function values. The highest and lowest values of the parameters presented in constraints (15)–(18) determine the minimum and maximum values for the objective function . The locations and numbers of traffic signs along with the speed limits and maximum durations of traffic lights are assumed to be known.
6. Numerical Analysis and Simulation
In this section, we apply the methodologies from Sections 4 and 5 to two different reallife routing cases, short and long distances, and analyze the results from computational experiments and discrete event simulations. The first case includes three different shortdistance routes that connect a location in Hamden () and another location in New Haven (), two small cities located in the state of Connecticut. One of the routes is via a highway, while the other two routes are via local roads. For this OD pair, we assume a travel pattern of two trips per day for all year. The second case includes three different longdistance routes that connect an origin and a destination in the states of Connecticut (CT) and California (CA), respectively. Note that the set of longdistance routes is chosen to mimic the actual truckload routes of a multinational stainless steel and specialty metal manufacturer located in the state of Connecticut. For this OD pair, we assume a travel pattern of oneround trip per two weeks for all year. In both short and longdistance routing cases, though there are many routes existing between OD pairs, we restrict our attention to only three routes to demonstrate the theoretical results and the benchmark with the stateoftheart application, Google Maps™. A summary of parameter values that are used in the computational experiments (unless otherwise stated, the parameter values are identical for both short and longdistance routing cases) is given in Table 5.

Several types of traffic signals are present in each route: traffic lights (TL), yield signs (YS), and stop signs (SS). Furthermore, under the SS category, four different types of stop signs are considered. The time spent in waiting at the traffic signals varies with the type. We assume minimum and maximum waiting times for each type of the traffic signals: traffic lights (15 seconds, 75 seconds), yield signs (3 seconds, 9 seconds), and jway stop signs (3j seconds, 3j + 3 seconds) where j is a type of a stop sign for all [26]. The average speed limits are considered for every route. All three highways that connect CT and CA are assumed to have 70 miles per hour speed limit, while the city roads are assumed to be 35 miles per hour, and the highway that connects Hamden and New Haven is assumed to be 55 miles per hour.
In addition, fuel costs () and operational costs () that include maintenance, repair, and depreciation are used to calculate the traveling costs for each route. From the American Automobile Association (AAA), we assume the fuel price to be $2.376 per gallon (AAA, 2019) and fuel consumption of 29 city/41 highway miles per gallon. We set the fuel consumption parameter, , equal to $0.0813 per mile (i.e., $2.376/29). If a vehicle travels mostly on highways, then is $0.0579 per mile (i.e., $2.376/41). To calculate the maintenance cost per mile, , we use two data points: (i) the national, maintenance cost average for a new Sedan which is $242 per year (AAA, 2019); (ii) annual miles traveled on an average which is 13,476 [27]. Using (i) and (ii), we set equal to $0.0179 (i.e., 242/13,476).
There are three scenarios considered for each evaluation criteria of the objective functions , , , and : best, worst, and average. The best and worst scenarios correspond to the objective function values when the minimum and maximum duration are used for waiting times at the traffic signals, respectively. The average scenario corresponds to the mean of the objective function values when the waiting times are allowed to vary randomly between the minimum and maximum. For instance, for each objective function, we perform 1000 replications of simulation of the waiting times and collect the average statistics.
The best, worst, and average scenarios are different for ; Algorithms 1–4 are used to calculate the total distance values for . This calculation has two phases: for instance, in the first phase, for every geodesic and linked chain circle combination, the minimum length of the geodesic (i.e., the minimum distance to travel from intersection to intersection) is determined; in the second phase, using the minimum length geodesics, the shortest path from to is identified. The best case value of is calculated using the shortest path geodesic (i.e., inner lane) distance, while the worstcase value is calculated using the longest geodesic distance (i.e., outer lane) on each route.
Since includes the cost incurred due to waiting at traffic signals, the best and worstcase values of correspond to the shortest and longest duration of waiting times at the traffic signals. Similarly, the waiting times at the traffic signals are included in the total travel time; as a result, best and worstcase values of correspond to the shortest and longest duration of waiting times. Note that there are no best and worstcase values for as the number of traffic signals does not vary on a given route. Next, we analyze the computational results from the shortdistance case.
6.1. ShortDistance Routing: A Case Study
For demonstration, we choose two locations, a coffee shop in Hamden and the train station in New Haven, where Hamden and New Haven are two cities in the state of Connecticut. There are three routes between the two places: (A) via I91 South—using highway; (B) via Whitney Avenue—using local roads; (C) via Ridge Road—using local roads. On each route, there is a certain number of traffic lights, yield signs, and stop signs.
Table 6 summarizes the computational results for the routes A, B, and C concerning criteria to . The best values are given in bold. We compare the results with Google Maps when the relevant information is available. We make the following observations:(1)When is considered, Route B is a good choice in all scenarios: best, worst, and average. The difference between the best and the worst cases is the highest for Route A which is 0.54 miles (11.38–10.84), about 5% per trip. In other words, if Route A is selected for traveling, using the innermost lanes can generate savings on the distance of about 5.(2)When is considered, Route C is a good choice in all scenarios. Though the simulated average values for all routes are close to each other, Route C has the minimum. Moreover, Route C has the minimum value in the best and worstcase scenarios. The waiting time information is not available from Google Maps. Note that the simulated average values obtained for has a significant impact on total travel time.(3)Route B is the most cost effective when is considered. The simulated average value for Route B is $0.727. When compared to Routes A and C, Route B generates an average cost savings of 0.185 and 0.166, respectively. With respect to , Route A is the best as it has a minimum number of traffic signs.(4)When is considered, Route C is a good choice in all scenarios. If the innermost lanes (i.e., bestcase scenario) are used to travel, both Routes A and C can save travel time when compared to Google Maps of the order of 1 minute (i.e., 6%). In the average scenario, Route C is a good choice; however, Google Maps performs better with a travel time of 17 minutes. The difference between our results and Google Maps is due to the variability in the waiting times that impact the total travel time. Further, we note that when the total distance is considered, Route A (highway) is the worst choice; however, it does well when travel time is considered due to the minimum number of traffic signals and increased speed limits.(5)Route A is a good choice with a corresponding value of 0.18 when equal weights are assumed for all the criteria (i.e., for all where i = 1, …, 5). The priority vector solution, , for Routes A, B, and C is the set BV = (10.84, 4.1, 0.868, 18, 15.93), (7.21, 5.65, 0.715, 29, 18.01), and (7.366, 3.3, 0.0.729, 22, 15.93). This set displays all the solution vectors and assist the decision maker(s) to examine the best options for all routes.

Overall, the above computational experiments for shortdistance routes highlight the complexity of decisionmaking in evaluating routes when multiple criteria are considered. With criteria to , a commuter can pick and choose from several travel options based on individual preference. The geodesic distances calculated based on the RMS approach are on par with the stateoftheart application, Google Maps. By including traffic sign duration and geodesic distances, our approach enhances the Google Maps’ results. Next, we analyze the computational results from longdistance problem instances.
6.2. LongDistance Routing: A Case Study
We choose an origin () and destination () located in the states of Connecticut (CT) and California (CA), respectively, and restrict our attention to three different longdistance routes that connect and : (A) via I80 West; (B) via I40 West; (C) via I70 West. Note that all the above pathways are highway routes as they connect the east and west coasts of the US.
Table 7 summarizes the computational results for the routes A, B, and C with respect to criteria to . The best values are given in bold. We compare the results with Google Maps when the relevant information is available. The following are the observations extracted from Table 7:(1)Route A is a good choice in the bestcase scenario (using innermost lanes) when is considered. On the other hand, in the worst (using outermost lanes) and average (simulation) scenarios, Route B is a good choice. When compared to Google’s best, Route A in the bestcase scenario which generates savings of 24 miles (i.e., 3003 − 2978.88) per trip. Though the savings per trip is minimal, when a large fleet of trucks and hundreds of deliveries are considered, the total savings can be a considerable amount. Further, the difference between the best and worstcase scenarios is at the maximum with Route A (3029.88 − 2978.88 = 51 miles). This difference impacts the total travel time, cost of traveling an additional 51 miles, and all the costs associated with the transportation (e.g., driver, truck maintenance, and fuel consumption costs).(2)Route C is the best choice in all scenarios when the total waiting time at traffic signs, , is considered. Google Maps does not provide this information; therefore, no comparisons are made. When is considered, the difference between the simulation average values for all routes is minimal. However, in the long run, a small difference can make a big impact. For instance, Route B generates a savings of $0.44 per trip (i.e., 295.96 − 295.52) when compared to Route A. For a fleet of 1000 trucks and 50 trips per year, the total savings is equal to 0.44 1000 50 = $22,000. A similar comparison of Route B and Route C generates a savings of $208,500. With respect to , both Routes A and C are the good choices as they have the minimum number of traffic signals.(3)Route A is the good choice in the bestcase scenario when the total travel time, , is considered. However, in both worst and average cases, Route B is a good choice. Furthermore, Google Maps suggest Route A or C, while simulation averages indicate Route B is the best when the total travel time is considered. Note that the results from Google Maps are higher than our time estimates due to the difference in the speed limits that is used to calculate the travel times.(4)When equal weights are assumed for all the criteria (i.e., for all where i = 1, …, 5), any route is a good choice since all of them have a corresponding value of 0.16. Following is the set of priority vector solutions, , for Routes A, B, and C: BV = (2979.8, 3.45, 293.489, 17, 3254.23), (2986.54, 3.95, 294.18, 19, 3261.99), and (3027.41, 3.25, 298.17, 17, 3305.88). This set displays all the solution vectors and assist the decision maker(s) to examine the best options for all routes.

In summary, the geodesic distances calculated based on the RMS approach with traffic signals considered for the routing solution outperform the stateoftheart application, Google Maps.
Though our model does not consider the realtime traffic conditions and congestion effects, the results show that the geodesic distances calculated based on the RMS approach can be helpful in enhancing the accuracy of the results from Google Maps, see Figure 4. For instance, using the innermost lanes (best case) on Route A, one can generate distance savings of 25 miles per trip when compared to the best result from Google Maps. Google Maps displays a single distance for traveling from the origin to destination, while the nature of the traffic roads suggests a distance to be traveled between two values: minimum and maximum geodesic distances. Though the savings may seem minimal, when the scope is increased from one trip to hundreds of trips per year and from one truck to the entire fleet of trucks, the magnitude of the savings can be significant.
(a)
(b)
7. Conclusion and Future Research
Transportation is critical for both end consumers and the business world. One of the key components of the transportation is to identify the best route to travel that connects an origin and a destination. For instance, a route can be evaluated using multiple criteria such as travel cost and travel time. So, for some, the “best” route is cost effective, and for others, it is time sensitive. The approach to identify the routes, the assumptions, and the underlying calculations is important as they determine the accuracy of the routing solutions. In this research, we consider several realworld driving factors such as the time spent at traffic signs (e.g., yield signs and stop signs), speed limits, and the topology of the surface to develop realistic and accurate routing solutions. Though these factors increase the complexity of modeling, they provide the flexibility to evaluate the routing solutions from different perspectives: cost, distance, and time, to name a few.
In this research, we develop a nonlinear, multiobjective, mixedinteger mathematical model (MINLP) with five objectives that minimizes travel distance, travel time, travel cost, time spent on traffic signs, and the number of traffic signs to design and evaluate the routes between an origin and a destination where the waiting times associated with traffic lights, stop signs, and yield signs are stochastic. In addition, we present a weighted sum formulation (WSF) and a best vector set (BV) approach to assist in evaluating routing solutions. A set of algorithms are developed to solve the MINLP that are based on Riemannian manifold surfaces (RMS) to factor in the Earth’s curvature when calculating distances. A modified version of the linked chain method (LCM) is used to determine geodesic distances that are important in calculating the distances for possible routes between an origin and a destination. The geodesic distances determined in the first four algorithms based on RMS are used in the MINLP model to ensure that the routing solutions are realistic and accurate. For instance, using the LCM, we generate inner and outermost lanes for highways that provide a range of traveling distances (instead of one distance) which has not been addressed in the past research.
We demonstrate the complexity of our algorithms and apply them in reallife instances. Two different reallife routing cases, short and long distances, for computational experiments and discrete event simulations are used for solving the MINLP and WSF and determining BV. The first case includes three different shortdistance routes, and the second includes three longdistance routes. In both cases, we show that the geodesic distances calculated based on the RMS approach are on par with the stateoftheart application, Google Maps. For instance, using the innermost lanes (best case) can generate savings on the distance traveled when compared to the best result from Google Maps. When a large fleet of trucks and hundreds of deliveries are considered, the total savings can be significant. In addition, we show that the MINLP, WSF, and BV assist decisionmaking in evaluating routes when multiple criteria are considered. We believe that there are significant research opportunities for our work to be applied in the areas of traffic planning/management for routing decision strategies.
The results obtained in this work highlight the impact of geodesic distance calculations and the role of traffic signs in evaluating the routes which cannot be overlooked. We show that the waiting times at the traffic signs can be significant and play a critical role in evaluating the routes for shortdistance transportation problems. For instance, Route B which uses the local roads to connect a coffee shop and a train station minimizes the total travel distance. But, the same route is a bad choice when the travel time is considered since the route has many traffic signs that create travel interruptions. On the other hand, for longdistance transportation problems, geodesic distance calculations play a major role. For instance, Route A, an actual truckload route of a multinational stainless steel and specialty metal manufacturer located in the state of Connecticut that follows the interstate highway I80, has the minimum travel distance of 2,978.88 miles based on the shortest geodesic. When compared with the longest geodesic distance of 3029.33 miles, Route A generates a savings of 50 miles per trip per truck (1.7%) which can be significant when the total number of trips and the fleet size are considered.
Our research has the following practical implications: (i) it is difficult to obtain the realistic and accurate routing solutions from the available mapping software when they do not consider the nature of the Earth’s surface. Our approach that is based on the Riemannian manifold surface addresses this issue by considering the shortest path geodesics. (ii) Given the magnitude of the transportation spending, i.e., 0.91 trillion dollars in 2014, and the pervasiveness of the distance calculations, our approach can be applied in many different areas such as the supply chain (e.g., network design), transportation (e.g., vehicle routing), renewable energy (e.g., market distribution), and alternative technology (e.g., locationallocation of charging stations for electric vehicles), to name a few. (iii) Furthermore, our multicriteria mathematical model where we include the waiting times at the traffic signs to evaluate the routes based on distance, time, and cost addresses the lack of consideration of the traffic signs in the stateoftheart applications (e.g., Google Maps).
Data Availability
A set of data is collected by using the distance calculator of Google Maps and locations readily available in Google Maps; therefore, the numerical results attained in this work as a part of the tables can be attained by other researchers by using Google Maps. The data set and the corresponding analysis can be shared by Dr. Emre Tokgoz via email (Emre.Tokgoz@qu.edu) if requested.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
Professor Tokgoz was funded by the Quinnipiac School of Engineering for completing this research.
References
 A. T. Kearney, State of Logistics Report, Penske Logistics, Washington, DC, USA, 2019.
 J.P. Rodrigue, The Geography of Transport Systems, Routledge, Taylor Francis, Abingdon, UK, 3rd edition, 2013.
 R. F. Love and J. G. Morris, “Mathematical models of road travel distances,” Management Science, vol. 25, no. 2, pp. 130–139, 1979. View at: Publisher Site  Google Scholar
 J. Brimberg and R. F. Love, “A new distance function for modeling travel distances in a transportation network,” Transportation Science, vol. 26, no. 2, pp. 129–137, 1992. View at: Publisher Site  Google Scholar
 J. Fernández, P. Fernández, and B. Pelegrín, “Estimating actual distances by norm functions: a comparison between the norm and the norm and a study about the selection of the data set,” Computers & Operations Research, vol. 29, no. 6, pp. 609–623, 2002. View at: Publisher Site  Google Scholar
 H. Üster and R. F. Love, “Application of a weighted sum of order p to distance estimation,” IIE Transactions, vol. 33, no. 8, pp. 675–684, 2001. View at: Publisher Site  Google Scholar
 B. Von Hohenbalken and D. S. West, “Manhattan versus Euclid: market areas computed and compared,” Regional Science and Urban Economics, vol. 14, no. 1, pp. 19–35, 1984. View at: Publisher Site  Google Scholar
 M. Miyagawa, “First and second nearest distances in Archimedean tilings,” Forma, vol. 29, pp. 5–11, 2014. View at: Google Scholar
 E. Alpaydin, İ. K. Altinel, and N. Aras, “Parametric distance functions vs. nonparametric neural networks for estimating road travel distances,” European Journal of Operational Research, vol. 93, no. 2, pp. 230–243, 1996. View at: Publisher Site  Google Scholar
 D. A. Griffith, I. Vojnovic, and J. Messina, “Distances in residential space: implications from estimated metric functions for minimum path distances,” GIScience & Remote Sensing, vol. 49, no. 1, pp. 1–30, 2012. View at: Publisher Site  Google Scholar
 S. L. Herica and G. B. Servio, “A method for modeling generalized distance functions,” in Proceedings of the IEEE International Conference on Service Operations and Logistics, and Informatics, vol. 2, pp. 3075–3079, Beijing, China, October 2008. View at: Publisher Site  Google Scholar
 J. Havelock, B. J. Oommen, and O.C. Granmo, “Novel distance estimation methods using “stochastic learning on the line” strategies,” IEEE Access, vol. 6, pp. 48438–48454, 2018. View at: Publisher Site  Google Scholar
 P. Toth and D. Vigo, “An overview of vehicle routing problems,” The Vehicle Routing Problem, SIAM Monographs on Discrete Mathematics and Applications, Society for Industrial and Applied Mathematics, Philadelphia, PA, USA, 2002. View at: Google Scholar
 R. Bowerman, B. Hall, and P. Calamai, “A multiobjective optimization approach to urban school bus routing: formulation and solution method,” Transportation Research Part A: Policy and Practice, vol. 29, no. 2, pp. 107–123, 1995. View at: Publisher Site  Google Scholar
 K. G. Zografos and K. N. Androutsopoulos, “A decision support system for integrated hazardous materials routing and emergency response decisions,” Transportation Research Part C: Emerging Technologies, vol. 16, no. 6, pp. 684–703, 2008. View at: Publisher Site  Google Scholar
 E. Demir, T. Bektaş, and G. Laporte, “The biobjective pollutionrouting problem,” European Journal of Operational Research, vol. 232, no. 3, pp. 464–478, 2014. View at: Publisher Site  Google Scholar
 A. A. Kovacs, S. N. Parragh, and R. F. Hartl, “The multiobjective generalized consistent vehicle routing problem,” European Journal of Operational Research, vol. 247, no. 2, pp. 441–458, 2015. View at: Publisher Site  Google Scholar
 R. S. Kumar, K. Kondapaneni, V. Dixit, A. Goswami, L. S. Thakur, and M. K. Tiwari, “Multiobjective modeling of production and pollution routing problem with time window: a selflearning particle swarm optimization approach,” Computers & Industrial Engineering, vol. 99, pp. 29–40, 2016.