Intelligent Intersection Control for Delay Optimization: Using Meta-Heuristic Search Algorithms

: Traffic signal control is an integral component of an intelligent transportation system (ITS) that play a vital role in alleviating traffic congestion. Poor traffic management and inefficient operations at signalized intersections cause numerous problems as excessive vehicle delays, increased fuel consumption


Introduction
Traffic congestion has become a threatening social concern in major urban cities around the world. It has severe negative consequences for smooth mobility, urban economy as well as increased environmental issues [1][2][3]. For example, in 2014 alone, the US economy suffered $160 billion loss, and over 560 billion pounds of harmful CO2 emissions to the environment in 2011 because of road traffic congestion [4,5]. Functional biodiversity is considered one of the primary index for environment monitoring [6]. In the European Union, traffic congestion costs accounted for approximately 1% of entire gross domestic product (GDP) [7]. In Europe, the concept of green revolution has become quite popular to achieve the goals of sustainable transport infrastructures in recent years [8]. Likewise, it was reported that influential Chinese cities cumulatively suffered a huge daily economic loss worth $1 billion because of traffic congestion [9]. Kingdom of Saudi Arabia (KSA) along with neighboring Gulf states have witnessed rapid rise in economy, increased population growth, and motorization rate since the oil boom in 1970s, which has led to less unsustainable travel patterns such as extreme traffic congestion, social costs, and environmental emissions [10]. According to a study, traffic congestion costs in the city of Riyadh, KSA were estimated to be around $8 billion [11], which is an alarming situation. The primary cause of increased traffic congestion in KSA is lack of public transport and infrastructure except few main cities (Jeddah, Riyadh, Mecca, and Madinah), which has consequently resulted in increased auto-ownership. In 2013, the auto-ownership surged to 220 vehicles per 1000 persons [12]. As per statistics of US Department of Commerce, KSA imported over one million vehicles in 2016 [13]. Rapid motorization in the country has also brought miserable safety concerns including increased road traffic crashes and fatalities [14]. Thus, there is an urgent need to take appropriate steps to make transport infrastructure in KSA more user friendly and environmentally sustainable.
Mobility, safety, and environmental impacts are critical concerns in urban transport networks, particularly at intersections. Intersections are essential component of urban road network, and aimed to ensure safety and smooth flow through urban corridors [15]. However, intersection management has become a global challenge in the present era of rapid motorization. Statistics from previous research conducted by fatality analysis reporting system (FARS) suggest that approximately 40% of total crashes occur near intersections, and about 22% of crash fatalities are intersection related [16]. A recent study conducted in Eastern province, KSA suggested that a large proportion of severe crashes at intersection are associated with driver distractions [17]. In the conventional traffic control paradigm, traffic at the intersection is regulated either by signs or through traffic signals. The vehicles are bound to remain idle and make frequent stops, thus restraining the overall capacity of intersections. In other words, intersection traffic lights are unable to adapt to varying traffic situations leading to extreme congestions, fuel waste, crashes, and exhaust emissions. In recent past, cluster analysis coupled with unsupervised classification has been successfully employed to lower CO2 emissions and to boost circular economy [18]. Signalized intersections are important points along the highway and streets. To evaluate the optimal performance and quality of operations is a challenging task. Various performance measures are adopted by traffic engineers to evaluate and quantify the performance of signalized intersections. Some of these performance measures include average queue length (number of vehicle waiting to cross the intersection), the average delay per vehicle (s/veh.) number of stops, speed and average travel through intersection, etc. It is estimated that, on average, delay at intersection accounts for 15-30% of all traffic delays [19]. Length of a queue is an important measure to determine when the discharge from an adjacent upstream intersection will begin to impede through a given intersection. Delay is a measure of extra time needed to traverse an intersection. The number of stops made during a journey is also an important input variable, particularly in air quality models.
Traffic light settings play a vital role in determining the capacity and quality of operating conditions (level of service, LOS) at signalized intersections. Aiming to minimize average vehicle delay or some other measures of effectiveness (MOEs like emissions, fuel consumptions) through intersections, various models for optimizing signal cycle length have been proposed in previous research studies. The first class of methods in this regard is based on linear and non-linear regressions, and other probability-based methods [20][21][22]. However, the performance of these methods has been frequently questioned under unforeseen traffic conditions (such as traffic incidents, special events etc.,), and oversaturated conditions. A second class of methods deals with optimizing cycle length under both saturated and oversaturated intersection conditions [23,24]. The third group of methods are utilized to optimize other performance measures like pollutant emissions, fuel consumptions by adjusting the cycle length [25][26][27]. The past decade has witnessed rapid development in information and communication technologies such as fourth generation (4G), dedicated short range communication (DSRC), Internet of Things (IOT) sensors that provide unprecedented opportunities for detailed data collection. Consequently, the fourth class of methods are based on simulations and intelligent algorithms to enhance operations at signalized intersections by optimizing cycle length, green splits and offsets [28][29][30][31]. Methods based on simulation and heuristics are more realistic, robust, and efficient in capturing the non-linear and stochastic traffic conditions [32,33]. Various methods discussed above for traffic signal control were mainly focused on lane-based homogenous traffic conditions. However, real-time traffic characteristics are highly stochastic, non-linear, and heterogenous which limits the large-scale adoption of these methodologies. Hence, current study is dedicated to develop intelligent traffic control methods based on meta-heuristics under heterogenous transient traffic conditions to optimize operations at signalized intersections.
During past three decades, KSA has witnessed rapid increase in motorization and auto ownership due to lack of public transportation, except few major cities resulting in huge annual billions of dollars loss to economy due to traffic congestion. However, heuristic-based traffic control has not been thoroughly studied so far, although some signalized intersections in urban metropolitans are operating under actuated traffic control. Thus, heuristic-based intersection traffic control is minimal because of the lack of research in local conditions. Without concrete and authentic research to demonstrate that how heuristic-based intelligent intersection control could significantly improve the prevailing miserable traffic conditions, the government officials, policy and decisionmakers have yet not been convinced to city-wide adoption of traffic control using such methods. Therefore, this study is conducted to assess the applicability and performance of heuristic-based intelligent intersection control to mitigate congestion through a case study in the city Dhahran, Eastern Province, KSA. Meta-heuristics were applied and validated to solve an important real-world traffic control problem. Specifically, GA and DE based demand-responsive traffic control schemes were proposed to optimize delay by rationally adjusting signal timings plan at two isolated signalized intersections. The proposed methods were compared in terms of effectiveness, robustness, and convergence to optimal objective function i.e., minimum delay prediction. Validation results in TRANSYT 7F proved the adequacy and superior performance of both methods compared to existing fixed-time traffic control. Although, from scientific view point this study may seem slightly less innovative, however, we firmly believe that it will likely to have significant societal impact and add value to scientific community in local context. Since, to the best of our knowledge, it is among the very first scientific research on current topic in the region.
The remainder of this paper is organized into different sections as follows. Section 2 provides an overview of related work including: strategies for signalized intersection control, detailed concepts of delay through intersection, and brief summary of previous studies focused on delay optimization. Section 3 describes the study area and data description. Section 4 first introduces basics concepts of GA and DE-based optimization, followed by step-wise procedure of search algorithms, and constraints setting. Section 5 provides results and discussions and validation of proposed methods. Finally, Section 6 presents a brief summary of key findings, study, and outlook for future research.

Signal Control Strategies
Controlling traffic signals at intersections needs a system for managing the duration of different phases of traffic signals in order to decrease the total waiting and total delay, and to maximize the flow efficiently. There are many methods proposed in the existing literature to overcome the problems associated with poor traffic management [34]. For example, a fixed time control method can be applied by using a fixed cycle length and fixed split of the green time duration for all the time plans of the signal. This method is appropriate for signalized intersection with stable and regular traffic arrivals. However, the traffic system is not a fixed system; it is changing from one time period to another, so the fixed traffic signal time method cannot efficiently fit actual dynamic traffic conditions. Real-time traffic-responsive control was introduced to overcome the difficulties of the fixed method [35]. This method uses some sensing technologies that collect the data of traffic at each approach. The concept here is that the traffic signal control managing is based on a specific control strategy according to actual-time traffic data. This method is appropriate and efficient for the conditions where traffic is not stable. In recent years, with advances in vehicle communication technologies, autonomous driving has become very popular for pro-active traffic management [36].

Delay at Signalized Intersections
Among all the performance measures, delay is most frequently used to evaluate the operational performance of signalized intersections [37]. Intersection level of service (LOS), a criterion used to indicate operating conditions, is established using average delay per vehicle. Beckmann et al. were the pioneers first to present an intersection delay equation based on queuing theory [38]. Webster et al. later on, extended their work to estimate delay and signal timings based on the rational assumption [39]. The assumption made for deriving delay formulae was that queue discharge rate increases until it attains a maximum value (termed as saturation flow rate) when the signal becomes green, and it remains stable until the entire queue is dissipated or the green phase is exhausted. This assumption made it easy to calculate lane group capacity, which is the product green to cycle time ratio (g/C) and the corresponding saturation flow rate. The formulae developed by Webster et al. is given by Equation (1).
where d is the average delay per vehicle in seconds, C is the cycle length or cycle time, λ is the proportion of cycle length which is effectively green (g/C), x is the degree of saturation, and q is the arrival rate in vehicle per second at the intersection. It may be noted from Equation (1) that total delay is composed of three terms. The first term (d1) is delay due to uniform vehicle arrivals, second term (d2) is the delay component due to random vehicle arrivals. In contrast, the third component (d3) is delay component in oversaturated flow condition, and can reasonably be replaced by 0.9, as suggested by one previous research study [40]. For the undersaturated flow condition, this component can be assumed zero. Webster formula provides delay estimates based on traffic flow phase-related valuation, while Akçelik later proposed movement-related model for delay estimation [41]. The delay equation proposed by Akçelik has been widely used as objective function for signal timing plan optimization is presented below in Equation (2): where D indicate the total delay (sec/veh.), qC is the average number of vehicles arriving per cycle (veh./sec), c is same as λ in Webster formula and represent green-time ratio (g/C), y indicate flow ratio (q/s), s is the saturation flow rate, No is the average overflow queue (number of vehicles waiting to cross the intersection), and x is the degree of saturation. The term No rises continuously till average arrival flow rate descends below capacity, and is predicted for both under-saturated and oversaturated conditions using relation given by Akçelik in his study [41]. After determining the total vehicle delay (D), consequently average vehicle delay was calculated using Equation (3): where d is the average vehicle delay (s/veh.), D is the total delay, q is the arrival flow (veh./s), and m is the total number of lanes at intersection.
Delay can be classified as uniform or random. However, methods based on the random delay concept are more rational and have been adopted by several previous studies [42,43]. Intersection delay can also be classified as stopped delay, approach delay, travel time delay, time in queue delay, and control delay. A stopped time delay may be defined as the time spent in the queue from the instant a vehicle is stopped by signal control and wait to pass through the intersection when the signal is green. Approach delay is the summation of stopped-time delay plus the time lost in deceleration while approaching the intersection and time lost for acceleration to gain the desired speed. A travel time delay may be defined as the difference in drivers' perceived travel time through an intersection and actual time spent while crossing the intersection. Time in queue delay is the total time lost from the instant a vehicle joins the queue until it clears the stop bar at an intersection. While control delay is the delay caused by either stop sign or traffic signal. These different delay measures may be significantly different from each other depending upon prevailing conditions at the intersection. Figure 1 presents schematic for different delay measures for single vehicles traversing an intersection.

Previous Studies
During the past few decades, numerous studies have been conducted to address different issues prevailing to traffic signals cycle length. Most of these studies were aimed to minimize the delay at intersections by optimizing cycle length [39,[44][45][46][47]. The most popular models used for signal time setting are Australian road research board (ARRB) and transport research laboratory (TRRL) models. The optimum cycle length from TRRL is calculated using the following relation: where; c0 is the optimal cycle length, L is the total lost time at the intersection, and Y is the summation of critical lane flow ratios. The highway capacity manual (HCM) also proposed a cycle length model and corresponding delay model [48]. However, the cycle length from the HCM model is based on expected saturation and thus does not provide a minimum delay. HCM cycle length is given by: where is the cycle length (sec); is total lost time per cycle, CS is the summation of the critical phase traffic-volumes (veh. /h); and RS is the reference sum flow rate, and can be found by relation (1710*PHF*fa). PHF in the above relation indicates the peak-hour factor, and fa indicates the area type adjustment factor, and its value is 0.9 for central business district (CBD) and 1.0 otherwise. Researches illustrated that there were some deficiencies in these analytical methods, and some novel algorithms and AI-based techniques should be utilized to solve equations that were proposed for the delay. Most of the existing traffic control schemes deploy fixed time program are based on historical traffic information without considering real-time traffic information [49]. Accurate short-term traffic sate prediction have been reported to have numerous applications for intelligent traffic control [50,51]. Numerous research studies have focused on minimizing delay through signalized intersections by rationally optimizing the cycle length. Studies prior to 2000, mainly concentrated on utilizing linear programming, and fuzzy logic [52,53], because of limited computational efficiency and powerful simulation tools. For example, Yu et al. and Mihaita et al. adopted a nonlinear programming procedure to address the signal timing optimization problem. With advent of advanced artificial intelligence methods and deep learning, researchers in recent years are more inclined to use these methods to mimic and control at signalized intersections. In conventional traffic control paradigm, traffic conditions are influenced by series of sequential actions, thus reinforcement learning (RL) has been widely opted since 1990s. El-Tantawy et al. presented a brief summary of evolution of RL methods for traffic light control during 1997 to 2010 [54]. Recently, several studies have explored the efficiency of RL and DRL for flow optimization through signalized intersections [55][56][57]. Lertworawanich et al. suggested novel methodology to optimize signal timing plan for oversaturated network. Maximum capacity, delay, and overflow at intersection were chosen as optimization objectives. Bargegol et al. proposed a novel method based on the kinematic wave theory to optimize the signal timing plan at the signalized intersection [31]. The study conducted that the proposed algorithm was more efficient in reducing delays and enhance the level of service for signalized intersections. Some of the widely used approaches for multi-objective signal timing optimization are; GA, simulated annealing (SA), and DE.
Recently, GA and ACO (ant colony optimization) have been investigated in several previous studies to optimize cycle length and other traffic performance measures [58][59][60][61]. Foy et al. used GA to identify minimum delay estimates by optimizing the cycle length under multiple scenarios [62]. Five GA traffic runs were executed for 60 GA generations with different initial populations. The authors found that the population was converging to optimum fitness between 20th to 30th generations. The average minimum waiting time found was around 40 s. It was concluded that GA could produce a rational traffic signal timing plan, and it can enhance existing traffic control techniques. Similarly, for optimizing traffic performance of real-time traffic signal control, Singh et al. used GA. Traffic emulator in JAVA was generated to represent dynamic traffic conditions. The traffic emulator was programmed to perform surveillance periodically at fixed intervals, and send the output data to GA, which provided green time extensions for each phase and optimized cycle length. Factors considered for GA were weights assigned to each road, total stops at each incoming lane, fixed cycle timings, and fixed maximum and minimum green timings. Their results showed that system application appeared to very useful. Under the given constraints, the developed model was shown to perform better than conventional fixed timed signals. Rahbari et al. in their study, attempted to reduce traffic congestion using GA [63]. The researchers developed binary-coded GA for a specified initial population, followed by the generation of a new population through mutation and cross-over, and finally, members with optimal fitness values were selected as a solution set. Putha et al. solved traffic signals coordination problem in oversaturated network by utilizing ACO [64]. ACO is based on swarm intelligence. The researchers found that ACO provide more reliable solutions compared to GA, particularly at large number of trials. Renfrew and Yu also proposed ACO to reduce vehicle waiting time at signalized intersection [60]. ACO outperformed fully actuated control signal, particularly under high traffic demand.
DE and AIS (artificial immune system) are some other heuristic search optimization techniques that have been recently used successfully for traffic engineering applications [65][66][67][68][69]. AIS is inspired by human biological immune system. Liu et al. proposed a bacterial foraging optimization (BFO) algorithm based on DE to obtain optimized estimates of delay at signalized intersections [70]. The competition and collaboration formed the basis for optimization in BFO. DE was used to update the bacteria located in the chemotax process to enhance convergence precision. The authors found that delay estimates from DE-based BFO were reduced by 28.3% from fixed time signals and 5.6% from GA delay values. DEBFO was also found to converge faster than GA, and its results were validated. Korkmaz et al. studied three types of differential evolution delay estimation models (DEDEM): linear, exponential, and quadratic [67]. Their results showed that the quadratic model yielded robust outcomes. The researchers further found that proposed DEDEM was capable of predicting the vehicle delay more closely in terms of relative errors between simulated and estimated values. It was concluded that the quadratic form of the DEDEM model can be used as an alternative estimation model for the delay. Ceylan also applied the DE algorithm for the optimal design of signal-controlled road intersections [71]. Results from his study showed that the DE is better than GA and Harmony-Search (HS)-based models in terms of the network performance. Yunrui also examined traffic signal control with DE. The authors found that the DE was more efficient to decide system parameters, and had a good model performance in terms of reduction of the delay and the queue length [72]. Louati et al. utilized immune network algorithm (INA) to optimize delay, queue, and throughput at signalized intersections under different traffic scenarios [69]. The authors found that proposed INA outperform conventional pre-timed and adaptive traffic control, and validated the results on microscopic traffic simulation platform VISSIM. Trabelsi et al. also investigated the performance of AIS for detecting and controlling the anomalous traffic conditions at signalized intersections [73]. Simulation results demonstrated the pertinence and robustness of proposed AIS approach.

Study Area and Data Description
The data utilized in this study were collected from two signalized intersections in the city of Dhahran, Eastern Province, KSA. The province has witnessed a rapid increase in motorization since the oil boom in early 1970, particularly during recent years, with the number of registered motor vehicles increasing from 2.9 million in 2011 to 5.2 million in 2017 [74]. The city of Dhahran has an area of about 100 km 2 , with a total population of approximately 0.3 million [75]. A significant proportion of the city's population is dominated by expatriates belonging to diverse driving and ethnic, cultural backgrounds producing heterogeneous road user populations. The city hosts Saudi ARAMCO (Arabian American Oil Company, the world's largest oil processing facility), academic institutions, and commercial activities, which are the primary sources of traffic productions/attractions in the area.
Intersection-I is located at the junction of Faisal Bin Fahd road and Dhahran Techno Valley road near the US consulate, while Intersection-II is stationed at the junction of Faisal Bin Fahd road with King Saud road (as shown in Figure 2). Both the intersections are located in mixed residential and commercial zones. In addition, both are four-leg intersections having at least two lanes from each approach. Traffic flow at these intersections is controlled and regulated by fixed/pre-timed signals. Two types of data were mainly acquired i.e., traffic volume data and signals cycle lengths, green splits, including the phasing sequence. Standard sated procedures were opted in collecting the field data. Traffic volume data were collected manually using a mechanical counting board for two hours during evening peak periods from 3:30 pm to 5:30 pm. Signal control data for both intersections were also collected manually. Cycle length for Intersection-I was 220 s, while Intersection-II had a cycle length of 160 s. Besides, traffic volume and phasing, data on vehicle queues and saturation flow was also collected. Figure 3 show the phasing scheme of traffic movements at both intersections in the study area. Intersections traffic and signal control data collected from the field is shown in Table 1. Average vehicle delay was 102.5 and 86.5 (s/veh.) for Intersection I and II respectively in existing condition, both of which correspond to congested regime conditions.

Methodology
The word "optimization" came from the same root "optimal" which indicate best. Thus, optimization may be defined as the process of achieving the best possible solution under given constraints. A typical optimization problem can be mathematically formulated as: where f(x) is the output/objective function to be optimized, gi(x) denote set of inequality constraints, hi(x) represent equality constraints. Constraints are limits within which the variables may be varied. The variables x1, x2, x3…., xn € R n (solution space) are the control or decision variables. By convention, optimization means a minimization problem, however it can also be designed as maximization problem by negating the sign of objective function. Optimization problem given above could be considered as decision problem that involves finding the best vector x of all control variables from the solution space. In this study GA and DE optimization techniques were employed to minimize the average vehicle delay (objective function) through isolated signalized intersections as a function of several input variables (i.e., traffic demand, existing phasing scheme, g/C ratio, saturation flow etc.), and constraints conditions (cycle length and individual green splits). Green splits in each phase (g1, g2, g3, and g4) were decision or control variables. The primary objective was to minimize average vehicle delay in response to approaching traffic demand by selecting best possible combinations of signal timings plan. GA and DE optimization program for current study were developed on MATLAB interface (version R2018b). The following sections provide a brief summary of architecture and working principle of both algorithms, step-wise respective algorithms procedures for current optimization problem, and constraints conditions for obtaining optimized objective function.

Genetic Algorithm (GA)
GA is a robust machine learning search method that is based on the process of natural selection for finding the optimal solution of the objective function [76]. A random population of solutions is initialized, which converges to a nearly optimal solution over a set of iterations or generations. The intrinsic characteristics of GA to mimic nature's evolutionary theory is based on "survival of fittest," [77]. In GA algorithm and other search optimization methods, the solution of the problem is stored by genes or chromosomes or genes in a population. While the population indicates the boundary or area of interest within which the number of genes or chromosomes can be determined. Fitness function is the rule of GA to select the gene with a fit solution. Every solution in the process is screened by the fitness function to calculate the fitness of the gene. The fitness of genes so obtained is then used for a relative ranking process that lists the genes according to their fitness. Ranking of genes allows the algorithm to select the best individual for next stage iteration by eliminating the genes having unfit solutions [78]. The population of genes in offspring is maintained the same as the parent population during reproduction. Transition to the next population generation may utilize either random selection or tournament selection procedure. In random selection, individual genes have equal chances of selection with no reference to fitness at all. While for the tournament, a group of "n" individual genes are first randomly selected to take part in a tournament, with genes having the best fitness values continued to preceding generation [79]. Next iteration offspring population is generated through the application of evolutionary operators such as cross-over, mutation, and elitism. The objective of these evolutionary operators is to ensure a better solution to the next population generation. The evolutionary operators are alternatively applied until the stopping criterion is satisfied.

Differential Evolution (DE)
Differential evolution (DE), like GA, is a population-based optimization evolutionary algorithm characterized by its simplicity, robustness, and fast convergence to the objective function. In recent years, DE has been successfully used for signalized intersection management [66,80]. DE is capable of solving non-differentiable and non-linear optimization problems more efficiently [81]. The search mechanism of DE algorithm is reinforced by a differential mutation operator, which utilizes vectors difference generated by the set of vectors of the original parent population. The resulting vector difference is combined with a randomly selected new individual to yield a new mutant population. This new mutant population is recombined with the individual from the parent population producing a new offspring population. The fitness function evaluates each solution from offspring and corresponding parent solutions, and the best solutions are carried to the next generation. There are several selection variants DE optimization algorithm [82]. However, widely used in existing literature are: DE/Rand; which represent traditional DE version that is based on base vector random selection with uniform probability; and DE/Best; that select base vector of best individuals for the next-generation population. In general, the production of successive offspring population from the parent population is achieved through the application of the same genetic operators i.e., cross-over, mutation, and elitism, just like in GA. Similarly, the process is repeated until either of the stopping criteria is reached.

Algorithm Procedure and Parmeters Setting
In general, a typical GA and DE follows a sequence of decisions summarized in the following steps. Figure 4 shows the schematic flowchart for the sequential decisions faced during the execution of both algorithms molded in terms of current optimization problem.
Step 1-Problem Encoding: It is the process of representing a solution in the form of bits that communicate essential information. Each bit in the string indicates characteristics of the solution just like each gene-associated chromosome control particular characteristic of an individual. Binary encoding is mostly used as the encoding method. Other encoding methods that have been used include, permutation encoding and value encoding.
Step 2-Initialization/Generation of Initial Population: It refers to the random generation of a population of potential/candidate individuals that belong to the solution of the optimization problem.
Step 3-Evaluation of Fitness Function: Individual candidate solutions are decoded and evaluated based on pre-defined quality criteria, commonly referred to as fitness or objective function.
Step 4-Selection for reproduction: These candidate solutions are then ranked according to their fitness values, and prepared for an application of evolutionary operators for next stage iteration.
Step 5-Cross-over: It is a genetic operator that varies programming of parent chromosomes aiming to improve the solution quality of the offspring population. It is a form of sexual reproduction with two strings randomly picked from the mating pool to yield superior offspring. Cross-over could be done single point or multiple points to identify points for the exchange of genetic materials at these points.
Step 6-Mutation: It is again a genetic operator inspired by a biological mutation that maintains genetic diversity from one population to the next. Mutation represents a small random tweak or distortion in parent chromosomes with a very low probability to diversify the search in global space. Mutation could be done through different approaches, such as bit-flip, random resetting, scramble, and inversion mutation methods.
Step 7-Elitism: It represents a genetic operator that copies a small proportion parent chromosome having the best fit value to next-generation, thus ensuring maximum fitness value does not drop from one iteration to the next.
Step 8-Testing for stopping criterion: Once a new generation is produced, the global best is updated. This criterion then decides on the termination of the algorithm. If convergence to objective function at a fixed number of generations/iterations is satisfied, otherwise, the procedure is repeated from step 3 onward. Table 2 summarizes the sequential stages of GA and DE algorithms along with values for individual parameters used in this study. Values for different parameters for both algorithms were kept the same to avoid any bias while comparing optimized delay estimates from both methods. It is worth to note, that value for parameters "selection parameters and alpha (α)" are not shown for DE since it does not require these parameters. Table 3 presents optimization objective functions, control/decision variables, and constraints. The constraint limits were deemed desirable and recommended by a previous research study [66].  Green split for phase III g4 Green split for phase IV   Figure 5 plots convergence to minimum objective function for GA and DE against number of generation/iterations. The objective function represents the minimum average vehicle delay (given in units of s/veh.) through the intersection. To comply with the central limit theorem, the convergence curves shown in Figure 5 indicate the mean rating of the best candidate solutions in the population at each iteration for 30 executions set for both algorithms [82]. Analyzing the convergence pattern from the plots (Figure 5), it is worth to note that both the curves converge rapidly to smaller and smaller values of the objective function. Both the curves become flat onwards after the objective function reaches a steady value as the number of iterations is increased. It can be observed from the plots that DE converges much faster than GA. However, GA solution quality from GA is slightly better as both the algorithms ultimately proceed to complete the fixed number of iterations. The fact that DE outperforms GA in terms of convergence to fitness function is consistent with several previous studies [83,84]. For Intersection-I (shown in Figure 5a Table 4 presents optimized signal cycle lengths and corresponding green splits for each bound. Green splits were the input decision variables for both algorithms and were programmed and controlled to have values between 8 to 45 s, as recommended by a previous study [66]. The objective was to select the best possible combination of green splits in response approaching traffic volume to minimize the objective function i.e., average vehicle delay at the intersection. It may be noted from Table 4 that the optimized green splits distribution from both GA and DE are intuitively based on the proportion of traffic volume from the respective approach. Intersection clearance interval consisting of "yellow and all red" was the same from existing conditions to ensure safety and smooth operations at intersections. There was a 15 to 42% reduction in green splits from an existing pashing scheme using GA for both intersections, whereas DE yielded percent difference ranging between 26 to 52%. Although percent decrease in individual green splits was more pronounced from DE, overall, reduction in average vehicle delay was less compared to GA.  Note: Values in parenthesis indicate percent difference from existing conditions. Figure 6 provides a comparison of average vehicle delay for existing conditions as well as the optimum delay estimates obtained from GA and DE at both intersections. Intersection-I had an existing cycle length of hefty 220 s with corresponding average vehicle delay through intersection being 102.5 s/veh. Signal control at this intersection is very poorly designed, and during rush hours, vehicles have to wait in long queues waiting to pass the intersection. Unless some preventive measures are adopted, in the near future, it is likely that the current prevailing level of service (LOS) D and E operating conditions could rapidly turn to breakdown LOS conditions, particularly during peak periods. This will cause chaos in the vicinity of the intersection. Optimization using GA estimated a cycle length of 148, having a reasonable average delay value of 64.3 s. While DE yielded an optimized signal cycle length of 125 s for Intersection-I with slightly high average vehicle delay value of about 78 s. Intersection-II had slightly better operating conditions having an existing cycle length of 160 s and associated delay of 86.5 s/veh. Again, reduction in optimized cycle length from DE (92 s) was more pronounced compared to GA (105 s). However, GA was found more efficient and robust, yielding significantly low estimates for optimum/minimum delay (55.9 s/veh.) than DE (72.8 s/veh.) algorithm. Table 5 presents the percentage reduction in cycle length and average vehicle delay for GA and DE algorithms. Optimized cycle length from GA was approximately 32% of existing cycle length at Intersection-I and about 34% of original cycle length at Intersection-II. DE, on the other hand, provided a cycle length reduction of a little over 40% for both intersections. Comparing the percent reduction in average vehicle delay, GA was capable of achieving a decrease of approximately 37% and 35% at Intersection-I and Intersection-II, respectively. While using DE, average vehicle delay was reduced by 23% at Intersection-I, and around 16% at Intersection-II.

Methods Validation
To assess the performance, and validate the adequacy and of proposed methods (GA and DE), TRANSYT 7F was used. TRNASYT 7F is the state-of-the-art macroscopic, deterministic simulation and optimization software tool, widely used for the analysis of signalized intersections. The software offers detailed information on various measures of effectiveness (MOEs) for signalized intersection like total vehicle delay, average delay, emissions, number of stops, etc. In TRANSYT 7F, a step-wise multi-cycle hill-climbing search optimization approach was employed to find global optimum in solution space. Figure 7 plots average vehicle delay as a function of cycle length for GA, DE, and using TRANSYT 7F. The relationship plotted in Figure 7 is well-intuitive since a traffic light with shorter cycle length in unable to provide desired green splits on individual approaches to clear entire vehicle queues in a single cycle. Thus, the resulting average vehicle delay increases. In contrast, an excessively large cycle length provides relatively longer green splits for one or all approaches, which again increases average delay for vehicles waiting idly in queues in other directions. This general pattern of cycle length vs. delay is consistent with previous studies [85,86]. Figure 7a,b present cycle lengths vs. average vehicle relationship for Intersection-I and Intersection-II, respectively. It is evident (from Figure 7) that GA and TRANSYT 7F curves closely follow a similar pattern of relationship, while the cycle length-delay curve for DE is slightly off. To be specific, GA provides a minimum delay of 64.3 (s/veh.) at an optimized cycle length of 148 s for Intersection-I, whereas TRANSYT 7F offers a minimum delay of 68.5 (s/veh.) that correspond to cycle length of 164 s.  It is evident from the above discussion and results in Figure 7 that both GA and DE algorithms are powerful tools for the analysis of signalized intersections. Simulation and optimization results from TRANSYT further reinforce and demonstrate the robustness of the proposed methods. Although GA and TRNASYT 7F have a very similar pattern of cycle length-delay relationship, however, global optimum (minimum) delay estimates obtained GA among all scenarios.

Conclusions and Outlook
This study aimed to optimize delay and achieve intelligent traffic control at isolated signalized intersections through application of meta-heuristic search optimization methods. A critical analysis of existing literature on the specific topic revealed that most of these studies have focused on lanebased homogenous traffic conditions. However, real-time traffic is usually heterogeneous, having non-linear, stochastic, and intricate characteristics. Thus, we proposed a couple of meta-heuristicbased including GA and DE methods for efficient traffic control at signalized intersections. Traffic volume and phasing scheme data were collected for two isolated signalized intersections in the study area, the city of Dhahran, KSA. Average vehicle delay through the intersection was chosen as the primary performance index. Optimum delay estimated were obtained by optimizing the signal timing plan (cycle length and corresponding green splits) in response to existing traffic demand. Both GA and DE yielded rational signal timing plans. The results indicated that both methods effectively reduced the average travel time delay ranging from 15 to 35% compared to existing conditions. Although DE converged much faster to the objective function, GA slightly outperforms DE in terms of solution quality i.e., minimum delay value. To assess the validity of proposed methods, cycle length and delay curves were compared with outputs from state-of-the-art signals traffic simulation and optimization package, TRANSYT 7F. Validation results demonstrated the adequacy of the proposed methods. Again, the average adaptive traffic signal timings from GA follow optimized outputs timings from TRANSYT 7F more closely. To conclude, findings from the current study suggest that traffic control at signalized intersections should be based on adaptive heuristics approaches rather than conventional fixed and actuated signals to ensure more efficient and sustainable operations.
The findings of this study could serve as essential guidance for policy and decision-makers, and other concerned, particularly in the study area to initiate concrete steps to improve the miserable traffic situations. This study could be extended in the future from various perspectives. Recently network optimizations have become very popular, so future studies could focus on network-wide optimization problems. Further, the current study did not consider the impact on non-motorized modes, which would be interesting to explore in the forthcoming studies. Lastly, studies could utilize multi-objective heuristic modeling to concurrently consider different performance measures such as number of stops, energy consumption, and emissions.