Multilane Microscopic Modeling to Measure Mobility and Safety Consequences of Mixed Traffic in Freeway Weaving Sections

Weaving sections are components of highway networks that introduce a heightened likelihood for bottlenecks and collisions. Automated vehicle technology could address this as it holds considerable promise for transportation mobility and safety improvements. However, the implications of combining automated vehicles (AuVs) with traditional human-driven vehicles (HuVs) in weaving freeway sections have not been quantitatively measured. To address this gap, this paper objectively experimented with bidirectional (i.e., longitudinal and lateral) motion dynamics in a microscopic modeling framework to measure the mobility and safety implications for mixed traﬃc movement in a freeway weaving section. Our research begins by establishing a multilane microscopic model for studied vehicle types (i.e., AuV and HuV) from model predictive control with the provision to form a CACC platoon of AuV vehicles. The proposed modeling framework was tested ﬁrst with HuV only on a two-lane weaving section and validated using standardized macroscopic parameters from the Highway Capacity Manual. This model was then applied to incrementally expand the AuV share for varying inﬂow rates of traﬃc. Simulation results showed that the maximum ﬂow rate through the weaving section was attained at a 65% AuV share. At the same time, steadiness in the average speed of traﬃc was experienced with increasing AuV share. The results also revealed that a 95% AuV share could reduce potential conﬂicts by 94.28%. Finally, the results of simulated scenarios were consolidated and scaled to report expected mobility and safety outcomes from the prevailing traﬃc state and the optimal AuV share for the current inﬂow rate in weaving sections.


Introduction
e rise of automated vehicle (AuV) technology, as an essential component of a new generation of traffic infrastructure, has been researched by academics and industry who recognize its advantages over the existing transportation framework regarding improved mobility and enhanced safety and reduced environmental impact. However, large-scale transitions to AuV technology-based transportation systems cannot happen overnight. Research on the fusion of such technology with the current, humanoriented transportation system that also considers the restraints of existing roadways is warranted. Certainly, freeway weaving sections are considered restraints since they act as recurrent bottleneck locations due to inherent vehicle trajectory patterns formed by vehicles changing lanes from auxiliary lanes to mainlines and vice versa. As such, substantial research has established that both mobility and safety of the weaving sections are compromised [1][2][3][4][5][6][7][8].
Although numerous studies have established the eminence of mixed traffic over traditional traffic systems from mobility, safety, and environmental perspectives, the exploration of coexistence is primarily limited to partial motion dynamics, most often car-following strategy, of studied vehicle groups [9][10][11][12]. is limitation is significant since both traffic operational and regulatory authorities must base their strategic investment and policy and legislative decision-making on sound, objective facts regarding integrating AVs with conventional traffic systems. Motivated by this gap in knowledge, this research, as part of a broader study, addresses the following question: how can the mobility and safety of varying traffic states in a multilane weaving section be influenced by the shared presence of AuVs and human-driven vehicles (HuVs)? e two-fold objective of this descriptive research includes the following: (i) contriving a comprehensive and realistic modeling framework of mixed traffic with bidirectional motion dynamic and (ii) quantifying and clarifying the causal connection between the presence of AuVs in traffic with a potential shift in mobility and safety benchmarks for a weaving section. e findings of this study aim to contribute to the existing body of knowledge, which currently lacks substantial evidence to support the effect of AuV integration to improve the mobility and safety of weaving freeway sections. Furthermore, the study outcomes can benefit research communities and industries actively committed to intelligent transportation systems and presume that AuVs will play a significant role in overcoming flow efficiency limitations and crash likelihood. e demand for understanding the implications of AuV substantiates the need for more comprehensive research from the mobility and safety perspective. ese perspectives play critical roles in measuring transportation system performance and the effects of mixed traffic on planning decisions. From planning and operational standpoints, traffic mobility is defined as the ability and level of ease of moving goods and services [13]. For instance, freeways can provide designated high occupancy vehicle lanes to increase the overall efficiency of moving people while maintaining the overall number of vehicles. Since a detailed evaluation of mobility shift would require a considerable mobilization of resources, a few mobility parameters are selected in this study to provide a general overview of potential adjustments to traffic mobility. On the other hand, implications on traffic safety can be measured directly from the number of collisions, injuries, and fatalities. Since AuVs are not yet widespread in present traffic, direct measures of potential changes cannot be obtained in most cases. Hence, surrogate safety indicators are adopted to measure positive or negative shifts in overall traffic safety.

Literature Review
Despite real-world pilots of AuVs and the significant advancement in knowledge on this technology, large-scale deployment of AuVs in contemporary traffic streams is not readily achievable, making much of the existing literature that studies the effects of mixed traffic reliant on traffic simulation. Model predictive control (MPC)-based microscopic traffic dynamic has been adopted in few studies to simulate AuV motion. Liu et al. [14] proposed a distributed MPC for cooperative highway driving considering safety, efficiency, and fuel economy. e cooperative lane-changing maneuver was the foundation block for the dynamics in a mixed traffic environment that recognizes the effect of surrounding HuVs on the cooperation procedure. Bertoni et al. [15] took an energy-saving eco-CACC approach to use nonlinear MPC to vary distance and speed between electric vehicles to minimize fuel consumption. Stanger and del Re [16] took a similar approach to demonstrate the fuel and roadway capacity benefits resulting from CACC. Kamal et al. [17] proposed an efficient driving system based on MPC where AuV received vehicle status of surrounding vehicles to generate optimal control decisions. eir study revealed that the proposed system improved the fuel economy and travel efficiency of the AuVs significantly.
ese studies provided valuable insights about roadway capacity changes of mainstream traffic resulting from mixed traffic flow at varying market shares, although few explicitly explores the influence of integrating AuVs into freeway weaving sections. Furthermore, the few studies that included car-following [10,12,28] lack the inclusion of weaving sections and their influence on lane-changing vehicles. In addition to capacity, the impact of AuVs on mobility has been evaluated through traffic speed [29][30][31]. However, the aggregated impact on both of these parameters was underexplored by these studies. Malikopoulos et al. [32] examined environmental implications, travel time, and traffic throughput in mixed traffic scenarios. Rios-Torres and Malikopoulos [33] discussed the environmental and mobility aspects of mixed traffic flow for merging segments in a simulated environment. Even though weaving sections are critical components of a freeway system, scarcely any study examined the possibility of mobility and safety paradigm shifts in these sections due to introducing AuVs in the traffic stream. Fazio et al. [6] used a simulated conflict rates count to identify hazardous locations, comparing this with obtained crash rates. Uno et al. [7] identified the potential conflicts in weaving sections by analyzing vehicle movements from recorded videos. Tilg et al. [34] proposed a mixed traffic model calibrated to replicate the traffic dynamics on a weaving section. Although findings from this study revealed the potential of AuVs to improve the capacity of the weaving section, other aspects of mobility and safety remain uncharted. Ye and Yamamoto [35,36] conducted studies on heterogenous traffic flow and concluded that the resulting capacity improvements depend largely on AuV market penetration and car-following parameters. Both studies discussed the changes in the macroscopic fundamental diagram to determine the changes in traffic flow parameters.
Several studies have been published in recent years that measure the safety impact of AuVs when mixed with human-driven vehicles (HuVs). e study conducted by Hayes [37] reported that fatality rates could be reduced to 1% of current rates once AVs reached a 100% market share. Alternatively, Fagnant and Kockelman [38] predicted that the influence of AuVs could reduce the crash rate by 90% with the elimination of human error possibility. Multiple research papers explored the simulation approach to estimate traffic safety [39][40][41][42]. Fan et al. [39] proposed a two-stage process to use the VISSIM simulation model outputs to calibrate for the surrogate safety assessment model at merging locations. A similar objective was followed by Huang et al. [40] for signalized intersections. On the other hand, Essa and Sayed [41] studied the transferability of the calibrated parameters to different sites for simulation. From the perspective of the safety impact evaluation of mixed traffic, Ye and Yamamoto [11] simulated mixed traffic to study traffic safety under various market shares of connected-automated vehicles (CAVs). ey argued that the cautious and accurate carfollowing strategy from CAVs would significantly contribute to traffic safety. Papadoulis et al. [43] developed a bidirectional decision-making control algorithm for AuVs to evaluate resulting safety implications at different market penetration rates. e results revealed that traffic conflicts could be reduced by up to 90-94% with a full AuV traffic stream. e safety implications of exclusive AuV lanes were identified by Zhang et al. [44] with different market shares and for multiple traffic demand scenarios. is study indicated that a higher number of lanes would be required at high-demand scenarios to attain significant safety improvements. In brief, the research efforts in measuring safety implications of AuVs in mixed traffic environment are summarized in Table 1.
An overall review of the relevant literature solidifies the common assumption that introducing AuVs in traditional transportation systems would generate favorable mobility and safety implications. However, uncertainties regarding the magnitude of improvements as well as potential negative consequences remain equivocal. Furthermore, most studies considered collective implications from different configurations of roadway segments, which diluted the influence of natural impediment and variability resulting from individual segment structure. Hence, this research consciously identified a prevalent bottleneck segment of a freeway and explored both mobility and safety implications from a system-level prospect.

A Framework of Multilane Traffic
e framework applied in this study incorporates longitudinal and lateral motion dynamics of both vehicle types for unimpeded movement along the roadway. In addition to mandatory lane-changing, platoon formation among AuVs is also configured within this framework. More definitive conclusions can be drawn from such a comprehensive modeling architecture, directly examining the mobility and safety implications of AuVs. Agent-based modeling was adapted from the MATLAB library to develop the core model structure of this study. Agents of two different types were defined to represent the two vehicle types (i.e., HuV and AuV). ese two agents enter and exit a predefined roadway segment with average headway input to maintain average inflow rate and AuV share for each simulated scenario. e roadway segments were defined by giving several lanes and length of section values as inputs. e driving strategy and vehicle model were similar for both agents, apart from acceleration/deceleration conditions and platoon formation provisions. Irrespective of the leading vehicle type, HuVs had an average desired headway of 1.4 sec with a standard deviation of 0.3 seconds to account for human driving behavior variability. While the desired headway of an individual vehicle would be fixed (with some exceptions for forced mandatory lane-changing) throughout the simulation period, the parameter would follow a log-normal distribution for the overall HuV proportion of traffic. Instead of a constant desired headway, the provision of varying desired headway distribution of HuVs within the modeling framework presented the opportunity to incorporate the behavioral variation of human drivers.
On the other hand, AuVs inherited desired headways from platoon configurations (i.e., interplatoon headway and intraplatoon headway) depending on their position within a platoon. As mentioned in Seraj et al. [10], the AuVs are programmed to form platoons amongst themselves. If there is an HuV in front of AuV, the subject AuV will pursue ACC with a small desired headway (i.e., 1.25 sec) in the car-following state. However, if the leading vehicle is an AuV, the subject AuV will reduce the desired headway (i.e., 1.0 sec) to form a CACC platoon. e maximum platoon formation length was 3 AuVs in the simulation. Also, intraplatoon distance is taken as 4 sec to accommodate lane-changing vehicles in between platoons. Hence, if the leading AuV is the 3rd AuV in the platoon, then the subject AuV would maintain a 4-sec headway. e agent-based modelling's default driving strategy and vehicle model were customized to develop more realistic vehicle dynamics. Various modeling approaches in the literature established the spatial-temporal anticipation ability of human drivers [51][52][53][54][55][56][57].
e anticipative nature of human drivers compensated for higher reaction time than the ACC/CACC system that can generate quicker reactions to an event [58]. Hence, for an advanced vehicle control system, incorporating such anticipative intelligence naturally improves driving efficiency and opportunities to avoid a collision. In response, the MPC-based motion algorithm is chosen for modeling both vehicle types. MPC utilizes present information to predict the future state by controlling the process through the minimization of objective function under constraints.
is study aimed at developing a control framework that combined connectivity and automation of AuVs to navigate effortlessly in mixed traffic scenarios. Our modeling framework used the MPC framework integrated into the MATLAB Stateflow model to decide between the different states of the driving strategy.
Each vehicle was assigned a vehicle ID to specify their presence on the road and vehicle type to decide on control provisions. e subject vehicle's longitudinal and lateral control was implemented by computing the current traffic state from input data and predicting its future pattern. A P-step prediction horizon and a C-step control horizon were considered for optimizing the control inputs of the subject vehicle. e prediction and control horizon of the MPC was taken as 30 and 3 timesteps, respectively. Since a lanechanging maneuver takes multiple seconds to complete, the following constraint was included: which implies that the subject vehicle was allowed to change lane only once either to the right (δ r ) or to the left (δ l ) lane within the prediction horizon (P). Some additional constraints were defined to ensure safety and comfort in driving. e speed of the subject vehicle at any timestep was bounded by the following upper and lower limit: Here, v max is the maximum achievable speed, which is 10% higher than the speed limit, and v SL � 25 m/s. Acceleration is bounded by the following boundary constraints:

Journal of Advanced Transportation
Previous studies on modeling HuV and AuV chose their limits to account for driving comfort [59,60]. To ensure safe driving in traffic, i.e., avoiding collisions with other vehicles, the minimum gap constraint was imposed as Here, Δg s is the gap between the subject vehicle and leading vehicle, Δg s, min is the minimum allowable gap for the subject vehicle, g 0 is the minimum gap in a vehicle at standstill condition (2.5 m), and h 0,s is the desired headway of the subject vehicle. e following nonlinear constraint was introduced to avoid collision risk during lane changing: or δ (r/l),t−1 (k), or δ (r/l),t (k) � 0.

(6)
Here, the first equation was used to estimate the acceleration of surrounding vehicles. e second equation implies that the subject vehicle assumes no other surrounding vehicle is changing lane at any timestep k. e acceleration of HuVs was updated by an enhanced Intelligent Driver Model (IDM) [61], and the model proposed by Hu et al. [62] was employed to determine the acceleration of AuVs: Reference Safety parameters Findings [45] (i) Time-to-collision (TTC); (ii) nondimensional warning index ACC strategy can prevent vehicle gaps from reducing to an unsafe level in different driving scenarios [46] (i) Time Exposed Time-to-collision (TET); (ii) Time Integrated Time-to-collision (TIT) 90% reduction of rear-end collision risks due to CACC system [47] (i) TTC; (ii) postencroachment time (PET) Low CAV penetration showed an increase in potential conflicts at lower headway and signalized intersections; higher CAV market share showed a global decrease [48] (i) TET; (ii) TIT Degradation from CACC to ACC had a significantly negative impact on longitudinal safety [43] (i) TTC; (ii) PET Estimated conflict rates reduced with a gradual increase in market share of CAV and finally reduced by 90-94% at 100% CAV penetration rate [49] (i) TTC; (ii) TET; (iii) TIT Mixed traffic flow can induce higher longitudinal collision risk at low market penetration of ACC driven vehicle due to individual heterogeneity [50] (i) Number of conflicts; (ii) longitudinal volatility Substantial safety improvement was observed at market penetration higher than 40% 4 Journal of Advanced Transportation rough equation (7), it was assumed that the vehicles surrounding the subject vehicle (i.e., leading vehicle in the current lane x s−1 , leading vehicle in the target lane x t−1 , and the following vehicle in the target lane x t ) are not changing lanes at any timestep k. e predicted states of all vehicles for ensuring the safety in the prediction horizon needed to be estimated using which denotes S as the state of traffic considering S N as the estimated state of the individual vehicle where Since the state vector of the individual vehicle included time-varying continuous variables (e.g., position, velocity, current, and destination lane of the vehicle) and an integer discrete variable, the optimization of both longitudinal and lateral control decisions became computationally demanding and uncertain in obtaining optimal solution within a finite horizon. erefore, based on the control requirement (i.e., longitudinal or lateral) at a timestep, respective control optimization was considered to estimate the state of the vehicles: lateral control: Here, v des is the desired velocity equal to the speed limit (i.e., 25 m/s). e two terms in longitudinal control ensure the vehicles drive close to the desired velocity with little or no acceleration. 1st term in equation (9) compelled the vehicles to attain the desired velocity to minimize the objective function. 2nd term of longitudinal control restricts blunt fluctuation of acceleration while in motion. Both these aspects are important during the car-following state. Equation (10) considered the additional condition of an acceptable gap in the target lane while changing lanes. While the 1st and 2nd terms in equation (10) were the same as equation (9), the third term in lateral control cost function (equation (10)) penalizes for unsafe lane changes in the form of Gaussian function. e value of coefficient α t was measured according to Kamal et al. [17], which defines the shape of the Gaussian function. Furthermore, planning for a lane change near the end of a horizon limit is not preferred. Since it is more likely that the predicted states of the other vehicles vary in the course, this may force the vehicle to give up on a lane change process before completion. Specifically, minimum steps to initiate a lane change can be imposed by suitably choosing δ. Once the vehicle is on course to lane change, the receding horizon approach is applied for successively relaxing δ. Finally, a lane change is executed in a predefined reference trajectory. e subject vehicle was provided with two sets of data at each timestep: (i) vehicle's status memory and (ii) environment inputs. e vehicle's status memory included acceleration, velocity, position, yaw rate, lane position, and destination lane information of the subject vehicle for the final 50 timesteps. Environment inputs included velocity, position, and lane position vehicles surrounding the subject vehicle in the weaving section. With those inputs, the subject vehicle decided between longitudinal and lateral control to update the vehicle status. Regarding the driver's control decisions, both longitudinal and lateral control decisions contained two states. When only the longitudinal control decision was active, the subject vehicles could be either in the free-flow or car-following state. e prerequisites of these states are listed as follows: (i) Free-flow state was triggered when the upstream area of the vehicle is empty, or more than enough safety gap was available for the vehicle to drive at 10% higher than the speed limit. (ii) Car-following state was activated when there was spatial constraint due to the presence of a leading vehicle that restricted the vehicle to drive under the speed limit.
For HuVs, a discretionary lane change was initiated when the vehicle drove at 10% lower than the speed limit for 5 sec. e discretionary lane-changing actions by HuVs were motivated by driving at the speed limit. AuVs sought lanechanging opportunities due to the possibility of forming platoons with the leading vehicle in the target lane. e lateral control state was active when the subject vehicle attempts to change lanes. Discretionary lane changing was omitted for the weaving section in this part of the research.

Journal of Advanced Transportation
Due to the inherent characteristics of weaving section traffic flow patterns, it was assumed that the likelihood of vehicles performing discretionary lane changing in this roadway section would be insignificant. erefore, only mandatory lane-changing occurred based on the destination lane assigned for the vehicle and the vehicle looking for lanechanging opportunities as soon as it enters the weaving section. In the lateral control state, there were two possible scenarios:

Model Validation
According to the Highway Capacity Manual (HCM), the weaving section length varies between 150 m and 750 m. e studied roadway segment was designed as a 600-m long twolane Type-A weaving section of a freeway. e speed limit of the roadway was taken to be 25 m/s (90 kph). e lanechanging vehicle ratio (VR) was assigned to be 10% of the flow rate. e destination lane of individual vehicles was assigned randomly during the entry of the simulated roadway segment. e agent-based modeling approach was taken for simulating different inflow rates and AuV shares, which provides the opportunity to investigate complex interactions between two distinct groups of vehicles. HCM suggests that the weaving flow rate should not exceed 2800 vehicle per hour per lane (vphpl) for a Type-A weaving section as a higher inflow rate would be prone to more frequent operational failure. Although this restriction on inflow rate was imposed specifically for roadways containing only HuVs, the maximum inflow rate of the simulation model was kept at 2800 vphpl since no specific instruction was provided for mixed traffic inflow rates in HCM. Since it was not established before the analysis that AuVs could effectively increase the vehicle movement capability of the roadway segment, increasing inflow rates over the suggested level could induce operational failure due to high inflow rates. Hence, the simulated flow rates were restricted from 1200 vphpl to 2800 vphpl with 100 vphpl increment. e varying inflow rate was simulated by adjusting the average headways of vehicle entrance for both types of vehicles. AuV market share was maintained by regulating the number of AuVs present at a time during the simulation. However, each inflow rate and AuV share scenario was simulated 20 times to offset the randomness of vehicle entrance pattern and human driving behavior in the simulation.
To validate the proposed traffic model, the traffic data generated by the model were compared to the standard used in HCM. Numerous simulations were run, having only HuV scenarios for varying inflow rates; the macroscopic parameters (i.e., flow, density, and space mean speed) were recorded. All these recorded data points were then plotted to develop the fundamental diagram (i.e., flow-density diagram) and determine the segment capacity of the simulated weaving section. As illustrated in Figure 2(a), the fundamental diagram of the model weaving section was generated from numerous passes of the developed model. Figure 2(b) was generated from HCM that illustrated the three-lane Type-A weaving section's confirmed capacity with a 100 kph free-flow speed. e capacity term was defined in HCM as the maximum number of vehicles that can pass a given point under prevailing roadway, traffic, and control conditions. Following the definition, the capacity of the simulated section was identified as 2088 vph from the flow-density diagram. Comparing the simulated weaving section's obtained capacity with Figure 2(b) demonstrated that the capacity value falls near the capacity curve for VR � 0.10 of the three-lane weaving section with 600 m length. Hence, the developed model was deemed to be consistent with realworld roadway and traffic scenarios.
To further consolidate the model validity, the Level of Service (LOS) criteria established for the freeway weaving segment in HCM was matched with obtained model output.
As stated in HCM, the capacity of a weaving segment is the result of flows that causes the density to reach Level of Service (LOS) E/F for freeways. Journal of Advanced Transportation E for freeway weaving sections and therefore complies with the requirement stated in HCM. Since the results from both analyses were consistent with the expected outcome, the developed microscopic model was regarded as an effective representation of real-world traffic and could therefore be applied for future analysis.

Mobility Implications
e mobility consequences resulting from mixed traffic movement were measured considering two key parameters: (i) maximum throughput; (ii) average speed of traffic. e inflow rate varied from 1200 to 2800 vphpl with 100 vphpl   increments in each step. Similarly, the AuV share was increased from 5% to 95%, with a 5% increment in each step. As a result, 323 unique scenarios of mixed traffic movement were generated. To account for the fluctuations resulting from driving behavioral variations of human drivers and AuVs' location and distribution on measured parameters, each simulation scenario (i.e., inflow rate and AuV share) was repeated 20 times. e analyzed outcomes of these scenarios were compared with the base case (0% AuV share) to compute the consequences of mixed traffic movements. To calculate the chosen mobility parameters of each scenario, the simulation period (18000 timesteps � 30 minutes) was divided into six 5-min intervals, and the maximum/ average parameter value was measured for each interval. Finally, the maximum parameter value was chosen by taking the maximum interval values for maximum throughput and average interval values for the average traffic speed.

Maximum roughput.
While much effort has been made to specify the impact of AuVs on the maximum throughput of freeways [9, 20, 24-26, 34, 68-72], the variability in results kept it as an open and debated question. Hence, this study attempts to confront the issue with a precise scope of estimation for Type-A weaving sections with a 90 km/h speed limit and restricted inflow rates. e simulated scenarios varied the AuV shares from 5% to 95% to maintain a mixed traffic environment and compared the performance with a base case to appraise potential changes. As illustrated in Figure 3, the variation in maximum throughput is nonlinear and much greater due to variations in AuV share than inflow rate changes. Maximum detriment in throughput was experienced at 10% AuV share (−3.88%, 1981 vphpl@ 1400 vphpl inflow rate) and maximum gain was attained at 65% AuV share (80.49%, 3720 vphpl@2200 vphpl inflow rate) in comparison to base case (2061 vphpl). We compared our attained results with Tilg et al. [34], where they evaluated the effects of automated vehicles in freeway weaving sections.
eir study found that the maximum throughput of the simulated weaving section increased from 1700 vphpl@0% AuV to 3300 vphpl@100% AuV (94.12%). Although a comparable level of increase was observed in this study, the maximum increase in throughput was observed at a lower AuV share (65%). Recently, Rezaei and Caulfield [73] drew similar conclusions from their study where optimal mobility was achieved for both off-peak and peak traffic scenarios at 60% AuV shares. Attained maximum throughput for increasing AuV shares was further compared with the analysis of Liu et al. [68], which estimated the theoretical capacity of homogenous freeway sections resulting from mixed traffic. e following equation was applied to calculate the theoretical capacity at varying AuV shares: Here, P AuV(l) , P AuV(f) , and P HuV are the probability of any vehicle being AuV platoon leader, follower, and HuV, respectively; HW AuV(f) , HW AuV(l) , and HW HuV are average safety headways of AuV platoon leader (i.e., 1.25 sec), follower (i.e., 1.0 sec), and HuV (i.e., 1.5 sec), respectively; δ is AuV share. Comparison between the theoretical capacity and the highest maximum throughput value for specific AuV share displayed a similar pattern up to 65% AuV share (Figure 3(b)). While theoretical capacity showed an increase, the simulation results demonstrated a steady decline. is phenomenon can be explained by the fact that our inflow rate was bounded at 2800 vphpl based on the HCM guideline restricting potential gains. Additionally, in the simulation, both vehicle types were designed to lane change with the same trajectory with a similar acceptable gap and safety headway conditions that could play a critical role in restricting maximum throughput after a threshold AuV share. At higher AuV shares, platoons formed by AuVs  Analysis of maximum throughput fluctuations due to varying inflow rates at a specific AuV share identified the highest variations at 45% AuV share, which gradually receded with increased AuV shares. In addition to behavioral variations of human drivers, represented by the desired headway parameter of HuVs, the competing lane-changing opportunity searches by both vehicle types reached the maximum level once the traffic achieved equitable shares of AuVs and HuVs, which explains the highest variation at 45% AuV share and lower variations at higher as well as lower AuV shares.

Average Speed of Traffic.
One of the key indicators of mobility is the average speed of traffic. An increase in the average traffic speed can be interpreted as a reduction in travel time and vice versa. Furthermore, an experimental study conducted by Jiang et al. [74], specifically with AuVs forming a 51-vehicle platoon, indicated that traffic instability could be better reported through average traffic speed than traffic density or spacing. In this regard, the average space mean speed of the weaving section was measured for each simulated scenario. Figure 4 illustrates the obtained average speed outcomes for different simulated scenarios. As observed in Figure 4(a), the average traffic speed resulting from various inflow rates was highly dispersed with lower AuV shares (5-20%). Each black dot for a specific AuV share represents the average speed at a particular inflow rate. e average speed value of a specific scenario (i.e., inflow rate and AuV share) was obtained by taking the average speed of all the repeated simulation runs for that scenario.
On the other hand, with increasing AuV shares, the resultant average speed became increasingly concentrated, shifting towards the speed limit. e dispersion of average speed at lower AuV share indicated instability in traffic flow resulting from sparsely spaced AuVs in traffic. e instability of traffic at lower AuV shared was also translated into lower maximum throughput at similar AuV share, as presented in the final section. However, with gradual increase of AuVs brought more stability, which was perceivable from a more compact average speed distribution of traffic. Figure 4(b) provides a more explicit demonstration of the interaction between inflow rate and AuV shares with the average traffic speed. Although the average traffic speed changes at the lower inflow rate were negligible, reductions in the average speed at lower AuV shares were prominent at high inflow rates. Comparing mixed traffic scenarios with base case average speed (82.75 kph) obtained a maximum reduction of 11.53 kph at 20% AuV share and 2700 vphpl inflow rate and a maximum increase of 7.22 kph at 95% AuV share and 2400 vphpl inflow rate. It should be noted that both types of vehicles were simulated to attain the desired speed, which was the speed limit (90 kph). erefore, the maximum possible increase in average speed was 7.25 kph within the design framework. Liu et al. [68] noted a minor improvement at a lower market share as compared to a drastic improvement at a higher market share. However, the study considered market share at 20% increments, limiting the understanding for AuV shares with more nuanced increments. In contrast to Liu et al., the findings from this study showed slower traffic movements at lower AuV shares in high traffic demand. ese events could have been triggered by the dispersed positioning of AuVs in traffic that could not form closely spaced CACC platoons with other AuVs. As a result, these isolated AuVs had to maintain ACC driving principles while accommodating lane-changing vehicles at high inflow rates. Tilg et al. [34] mentioned that these gap-searching vehicles, specifically HuVs, related to speed attenuation until they engaged in successful lane-changing maneuvers. e gradual increase of the AuV share improved platoon forming probability that led to more stable traffic movement and is analogous to the findings by Talebpour and Mahmassani [75]. ose authors demonstrated the instability in traffic flow resulting from lower connected and automated vehicle shares. On the other hand, Spiliopoulou et al. [76] concluded that almost all congestion of the studied corridor could be eliminated at 60% market penetration. Our study also demonstrated a similar pattern since the average speed over the simulation duration of most scenarios beyond 60% market share was close to the speed limit.

Safety Implications
Traffic safety assessment is vital when analyzing unfamiliar transportation system performance.
e developed microscopic simulation framework provides the opportunity to perform a detailed safety implication analysis on simulated mixed traffic scenarios without physically implementing the AuVs in a transportation system. However, two major limitations of the simulation models are (1) the discounting of collision events among simulated vehicles and (2) the inability to simulate human drivers' distraction and misjudgement errors. Hence, the safety analysis of simulated models is dependent on interpreting and comparing surrogate measures of safety. In this research, we have used time-to-collision (TTC) as the primary measure of effectiveness for safety as its use in this capacity has significant precedence in the literature. Safety-critical TTC events were detected from simulated scenarios with maximum throughput and then translated into relatable safety parameters such as potential conflict events and rear-end crash potential.
ere are four types of vehicle interaction scenarios in a mixed traffic environment: HuV with an HuV (potential) leader (HuV-HuV), HuV with an AuV (potential) leader (HuV-AuV), AuV with an HuV (potential) leader (AuV-HuV), and AuV with an HuV (potential) leader (AuV-AuV). However, any commonly accepted threshold of surrogate safety measures has yet to be established for AuV in the literature. Hence, in this study, we only examined the interactions of HuV with other HuV and AuV leaders for potential conflict identification. Since the market share of AuVs changes with different scenarios, the crash risk rate for HuVs and conflict frequencies are used to identify the safety implications of AuVs in the mixed traffic environment.

Potential Conflict
Events. TTC, defined as the expected time before two vehicles collide if they remain on the same path at the same speed, is a widely used surrogate safety measure.
e Surrogate Safety Assessment Model [77] suggests 1.5 sec TTC as the threshold value for identifying potential safety concerns within a simulation environment, adopted in our research to extract the number of possible conflict events from simulated scenarios. e identified possible conflict events were classified into two collision groups: (i) rear-end and (ii) lane changing. e driving state of the subject vehicle was checked to classify the potential conflict type. If the TTC value fell below the threshold during the free-flow or car-following longitudinal control state, it was classified as a rear-end conflict event. On the other hand, if the TTC value fell below the threshold during the initial lateral control state, it was classified as a lanechanging conflict event. e following equations were used to measure and identify a potential conflict event: Here, CE s (k) is the conflict event count of the subject vehicle, s at timestep k. x i , v i , and g i represent i vehicle's position, velocity, and lead gap (current lane), respectively. Since the simulation was repeated several times for specific traffic flow scenarios (i.e., inflow rate and AuV share), the average number of conflict events was reported for each particular scenario. e aim was to reduce the stochastic effects resulting from the behavioral variation of HuVs and placement changes of AuVs. Since the threshold TTC value for identifying the conflict events was lower than the desired headway of AuVs (i.e., ACC � 1.25 sec and CACC platoon � 1.0 sec), the conflict events that emerged from AuVs were disregarded from the total conflict event computation. Hence, the identified conflicts were generated exclusively from HuVs interacting with other HuVs and AuVs in the vicinity. e gradual decline of conflict events experienced an exponential pattern with AuV share expansion ( Figure 5(g)). While larger shares of conflicts were found to be lane-changing (base case: 38.1% rear-end, 61.9% lane-changing), the presence of AuVs was relatively more proficient in mitigating rear-end type conflicts (95% AuV share: 22.6% rear-end, 77.4% lane-changing) by HuVs. e potential conflict events were further explored by the proposed methods of Oh and Kim [78] to estimate the crash risk index (CRI) from the TTC values in the analyzed roadway segment. While it was possible to obtain the deterministic control state of each vehicle under the current modeling framework, the probabilistic estimation of lanechanging decisions was retained for the novelty of the proposed method. e following equations were used to determine the crash risk index of a specific traffic scenario: Journal of Advanced Transportation 11 AuV share = 5%

Conflict event counts
Rear-end Lane-changing (iv) P(C s |TTC s ) is the probability that the subject vehicle would collide with the front vehicle given the current TTC (v) CRI i is the crash risk index of scenario i (vi) K is the total timesteps of analysis (vii) S is the number of vehicles passing through the weaving section Further details and parameter values are available in Oh and Kim [78]. In applying this method, the CRI values for each mixed traffic simulation scenario (i.e., inflow rate and AuV share) were measured from identified potential conflict events. TTC values more than the threshold TTC (i.e., 1.5 sec) were discarded from the CRI calculation. For each traffic state (i.e., inflow rate and AuV share), the simulation instances with maximum numbers of potential conflict events were chosen to measure CRIs by plotting them, as illustrated in Figure 6. As shown in Figure 6(b), the overall CRI showed a nonlinear downward trend with increasing AuV shares. e maximum CRI was experienced at a 5% AuV share (average � 0.0207, standard deviation � 0.0143), which reduced significantly at a 95% AuV share (average � 0.0019, standard deviation � 0.0014). At 5% AuV share, the average reduction in the CRI from the base case (average CRI � 0.0218, standard deviation � 0.0149) was 4.82% that gradually increased to a 91.35% reduction at 95% AuV share. For specific AuV shares, the CRIs experienced an upward trend with increasing inflow rates. Findings from this analysis conform to the findings of Papadoulis et al. [43], Ye and Yamamoto [11], and Rahman and Abdel-Aty [79], who all concluded, based on different sets of parameter analysis, that increasing AuV share could greatly improve traffic safety conditions.
While both analyses discussed so far showed significant improvements in overall traffic safety with a gradual increase of AuV shares, the findings were incomprehensible when determining the implications of AuVs on HuVs' perceived safety. Since the movements of AuVs were regarded as safe, irrespective of TTC values, due to their instantaneous response capability, an overall reduction of potential conflict events and crash risk for all traffic was anticipated with increasing AuV shares. Hence, adjusted CRI values were measured by considering only the number of HuVs passing through the weaving section to obtain more profound insights. e analysis results are presented in Figure 6(c). As portrayed in this figure, the CRI from the human drivers' perspective did not present as straightforward improvements as in the previous analysis. While the pattern of higher CRI with increasing inflow rate persisted for a specific AuV share, the reduction of CRI with increasing AuV share was renounced. e minimum average CRI for HuVs was experienced at 70% AuV share (0.0083) with minor changes until this point. However, the HuVs encountered substantially higher safety issues at high inflow rates from 75% AuV share upwards (at 5% AuV share). e maximum CRI at 95% AUV share was 0.0966 that occurred at a 2800 vphpl inflow rate, more than two times higher than the maximum CRI (0.0401) at a 5% AuV rate. While the overall reduction in CRI up to the 75% AuV share was moderate, the increments in CRI for HuVs were more drastic at high inflow rates with higher AuV shares.  [80,81] introduced TET and TIT, which have been widely used in the traffic safety literature to evaluate rear-end crash potential in traffic scenarios. TET is the summation of instances where the TTC is lower than a threshold value. In the previous section, the conflict events were identified similarly. However, the difference between the identified rear-end conflict events and TET measurements for rear-end crash potential stems from the fact that multiple consecutive TTC values for the subject vehicles lower than threshold TTC were counted as a single potential in the conflict event identification process. On the other hand, individual timesteps are considered in the TET calculation. erefore, higher TET values could be observed with lower rear-end conflict events, which would indicate higher safety concerns. TIT measures the value of TTC lower than threshold TTC. Similar to TET, lower TIT values are expected in safer traffic conditions. e values of these parameters were measured using the following equations: 14 Journal of Advanced Transportation TET s , TIT s , Here, TET i and TIT i are the TET and TIT values for traffic scenario (i.e., inflow rate and AuV share) i, TET s and TIT s are the TIT and TIT values for subject vehicle s, TTC * is threshold TTC (1.5 sec according to [77]), and TTC s (k) is TTC of subject vehicle s at timestep k. Aligned with the previous safety analysis, this analysis measured the TET and TIT values of simulated traffic for all 323 studied scenarios. e simulation instance with maximum throughput was considered for each scenario. Both TET and TIT counted for mixed traffic volume over the duration and adjusted for changing HuV proportions of traffic to make a fair comparison. As observed in Figures 7(a) and 7(c), both TET and TIT of mixed traffic gradually decreased with an increasing share of AuVs. Generally, for a specific AuV share, TET and TIT for mixed traffic increased with increasing inflow rates.
is pattern for specific AuV shares persisted through adjusted TET and TIT. Figures 7(b) and 7(d) account for gradually diminishing shares of HuVs in traffic. Nevertheless, adjusted TET and TIT that had changed due to increasing AuV share were dissimilar to those for the mixed traffic scenario. Although the TET and TIT were reduced considerably from mixed traffic volume, these parameter values increased once considered from only HuVs, particularly at high AuV shares.
is analysis demonstrated the dichotomy of crash potential due to changes in perspective. In plain sight, both TET and TIT showed significant reductions (from average TET �100 and TIT � 0.1318 at 5% AuV share to average TET � 3.3 and TIT � 0.0124 at 95% AuV share) for the mixed traffic passing through the weaving section. However, in-depth observations of these parameters from the HuVs' perspective revealed stark dissimilarities at higher AuV shares (from average adjusted TET �105 and TIT � 0.1388 at 5% AuV share to average adjusted TET � 65.88 and TIT � 0.2472 at 95% AuV share). It is also important to note that rear-end crash potential was nonexistent for the majority (1200-2100 vphpl) of the simulated inflow rate. Hence, the increasing trend in rear-end crash risk by HuVs was pushed by extensive AuV presence at high inflow rates. Moreover, the analysis disclosed the ineffectiveness of AuVs to influence the perceived safety of HuVs at high inflow rates in weaving sections, even with leading market share.

Maximizing Collective Mobility and Safety Implications
Due to the complex correlation of simulated traffic scenarios (i.e., inflow rates and AuV shares) with mobility and safety implications, developing a closed-form objective function to attain combined optimal benefits is considerably complex and time-consuming. Hence, an alternative approach of iterative search with metamodeling was taken in this study.
Metamodeling is a macromodeling method used in the literature [82][83][84][85] that aims to build simple and computationally inexpensive models, replicating the observed correlations when samples of a complex, high-fidelity model or simulation are drawn. Metamodels are often referred to as approximation, surrogate, and response surface models. Different metamodeling techniques include generating analytically inexpensive approximations of computationally intensive true response through various machine learning methods such as polynomial interpolation, support vector regression, kriging, and neural network [86][87][88][89]. is approach has been commonly used for solving simulationbased optimization and analytical dynamic equilibrium problems in transportation applications. e goal of including this modeling in this research is to produce a faster and simpler approximation of simulation-generated results to achieve multiobjective optimization and design space exploration feasible.
Mobility and safety scores were calculated for the range of simulated traffic scenarios from the attained parameter values of chosen mobility and safety parameters. To achieve balanced mobility benefits, both the average travel time and average travel distance were leveraged to calculate mobility scores. Safety scores were achieved by measuring the reduction of CRI from the base case. e following equations were used to measure the mobility and safety scores of each inflow rate and AuV share scenario: Here, MS (i,j) and SS (i,j) are mobility and safety scores for inflow rate I and AuV share j; ρ (i,j) (t) and v i,j (t) are the traffic density and average traffic speed at the weaving section at timestep t for inflow rate I and AuV share j; ρ (i,0) and v (i,0) are the traffic density and average traffic speed at the base case (0% AuV share); v SL is the speed limit of the weaving section. For a better performance by the machine learning algorithm, both scores were scaled within a 0 to 1 range, as they would be supplied as training and testing samples of a neural network (Figures 8(a) and 8(b)). e neural network was developed in MATLAB with inflow rates and AuV shares as input and corresponding mobility and safety scores as output. Seventy percent of samples were used for training, and the remaining samples were used for testing purposes. A Bayesian regularization algorithm was chosen with 20 hidden layers. e performance of each iteration was measured by mean squared error. A final model comparison between predicted and expected values showed R 2 � 0.9884 for training samples and R 2 � 0.9786 for testing samples. e purpose of developing such metamodels was to regress the response surface that characterizes the correlation between decision variable inputs (i.e., inflow rates and AuV shares) and simulation outputs (i.e., mobility and safety score). e finalized training model was then applied to obtain scaled mobility and safety scores for specific inflow rates and AuV shares in the weaving section. Furthermore, the model imparted optimal AuV share information to achieve maximum collective mobility and safety benefits for the specified inflow rate. Suppose the upstream average flow rate of the weaving section 1760 vphpl and AuV share of the traffic is 46%. e model provided the expected mobility and safety scores, which   (Figure 8(c)). To achieve more efficient flow and safer traffic movements through weaving sections that may act as bottlenecks on the freeway network, this information would be crucial for traffic operation and management authorities to impose controls enabling inflow rates and AuV share.

Conclusion
is paper proposed a microscopic modeling framework of mixed traffic that was applied to a freeway weaving section. Multifaceted traffic scenarios of mixed traffic were simulated for the weaving section, a notably underexplored combination in the literature. e modeling framework has been built on MPC-based decision-making principles that account for the anticipative aspects of human drivers. is modeling capability can be attributed to lateral control decisions and control model transitions of both vehicle types. Specifically, the cornerstone of reproducing realistic traffic states for a multilane roadway was anticipatory lanechanging decisions while perceiving the future status of neighboring vehicles and reference lane-changing trajectories built upon empirical observations. In this study, the simulation results contributed to understanding the effects of AuVs on traffic mobility and safety as well as validating the findings from previous relevant studies, such as Tilg et al. [34], Liu et al. [68], Spiliopoulou et al. [76], Papadoulis et al. [43], Ye and Yamamoto [11], and Rahman and Abdel-Aty [79]. While our findings were consistent with this research, our study contributes to the body of knowledge by addressing several underrepresented issues, including dichotomic interaction between mobility and safety in mixed traffic, pseudosafety perception with higher AuV shares, and metamodeling to attain equitable mobility and safety advantage from AuV presence in a weaving section. Analysis results revealed that a low AuV share with high inflow rates could undermine the mobility of the weaving section. is finding was critical for understanding the place of the weaving section in these traffic scenarios due to the distinct and preemptive lane-changing activity experienced. Safety features of traffic in weaving sections were found to be exponentially associated with AuV shares. While increasing traffic demand could raise safety concerns at any level of AuV presence, the range of possible conflict levels was substantially reduced with higher AuV shares. Finally, the simulation results were assembled to develop an assistive application for traffic operation and management authorities, aiding traffic state evaluation from both mobility and safety perspectives. e application can also seek optimal levels of AuV presence to maximize the impact. Altogether, the analysis results were reported to provide clear insight into the implications of AuVs on weaving section traffic mobility and safety and clarify the transformation of mixed traffic flow dynamics with the gradual adoption of AuVs under the current traffic system. One limitation of this study is the counterintuitive result of attaining maximum throughput at 65% AuV share. While the simulation model validation containing HuVs showed adequate resemblance to real-world traffic conditions, mixed traffic conditions may substantially differ. e current version of the developed model could have limited success in replicating the complex interaction between AuVs and HuVs. erefore, further improvements are required prior to policy decision-making and should be based on real-world corroboration of AuV and HuV interactions.
Future research includes further expansion of human driving behavior and the identification of resulting variations. Furthermore, we will continue the research by considering different platoon strategies and distinguishing the impact of different platoon structures and varying human behavior. While this study examined the mobility and safety implications for traffic containing two types of vehicles, the mixed traffic scenario will be enhanced to address HuVs with an advanced driver assistance system that improves driving efficiency without absolute reliance on automation.
Data Availability e simulated data used to support the findings of this study are available from the corresponding author upon request. Disclosure e contents of this paper reflect the views of the authors who are responsible for the facts and the accuracy of the data presented herein. e contents do not necessarily reflect the official views or policies of NSERC, City of Edmonton, and Transport Canada.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.