Modeling and Control of Network Macroscopic Fundamental Diagram during Holidays: A Case Study of Qingming Festival in Tianjin

: In this paper, the macroscopic trafﬁc states and network trafﬁc dynamics during the Qingming Festival holiday are explored with the Macroscopic Fundament Diagram (MFD), including the weekday before the holiday (WBH), the day of Qingming Festival (DQF)


Introduction
During holidays (such as New Year's Day, National Day), the travel behavior and pattern are different from those of ordinary weekdays, and the distribution of traffic on the urban network is heterogeneous [1,2].It is essential to identify the congested area of the network and clarify the dynamic evolution characteristics.On this basis, control strategies for regional traffic flow could be implemented to relieve large-scale congestion.In recent years, the widespread collection and application of traffic data enables urban traffic flow to be analyzed on a large scale.Geroliminis and Daganzo [3] used the field data from Yokohama (Japan) to verify the existence of the Macroscopic Fundamental Diagram (MFD) in a real network, which provided a novel way for the study of network traffic flow.Later, different methods have been proposed to investigate the existence, characteristics, influencing factors of the MFD based on field data, theoretical analysis, and numerical simulations [4,5].Moreover, several studies have been conducted to develop traffic management strategies based on the MFD [6].The MFD describes the stable and reproducible relationship between the average flow, average density, and average speed of a network in a specific area.Take the advantage of the MFD, the states of network traffic flow can be predicted based on concise data.Furthermore, an important conclusion is derived that there is a constant ratio between the unit production and the rate at which cars reach their destinations, including cars leave the region along the perimeter and internally [3].Based on this conclusion, the MFD could be explored for improving management and planning control strategies at the regional level.With the aid of the MFD, the macroscopic characteristics and dynamic evolutions of traffic flow could be clarified [7], and perimeter control strategies could be designed to alleviate regional traffic congestions [8][9][10].In terms of network dynamic evolution and equilibrium analysis, many studies constructed two-region [11,12] or multi-region [13,14] dynamic evolution models, theoretically analyzing the equilibrium and stability and revealing the impacts of heterogeneity on the MFD.Furthermore, in order to deal with traffic congestions among different regions, some macroscopic traffic models were proposed to study regional perimeter control [15][16][17].The controllers were designed to manipulate the ratio of transfer flows on the borders of two neighbor regions.Studying two-region models, Geroliminis et al. [17] and Yang et al. [18] used the MFD to design the perimeter control and signal control strategy in the two-region system, which solved the traffic congestion problem in the core area and improved the mobility of the traffic flow in the region.Studying multi-region (more than two regions) models, Guo et al. [19] and Ding et al. [10] proposed a perimeter control for multiple regions to solve the problem of regional congestion and reduce the total delay of the network.
The aforementioned literature regarding the MFD mainly concentrated on ordinary weekdays or weekends, whereas studies focusing on specific holidays with some traditional activities have not been fully developed.In terms of data analysis, Shim et al. [20] observed distinct characteristics of the MFD on weekdays and weekends with the empirical data in the CBD of Daegu.The MFD showed a highly aggregated relationship on weekends, but exhibited bifurcations in a high-density state on weekdays.Wang et al. [21] analyzed the network data of Sendai City throughout one year and found that the characteristics of the MFD on weekdays and weekends were different.On weekdays, hysteresis loops appeared only in the morning and evening peaks, whereas a larger hysteresis loop appeared throughout the whole day on Saturday.On Sunday, there would be no hysteresis.Regarding the research on the network traffic flow characteristics during holidays, Shi et al. [22] investigated the MFD on expressways in Shanghai for one week (26 September to 1 October 2009, including the Chinese National Day) of different days.They also found two hysteresis loops in the morning and evening peaks on weekdays, while during the National Day, the hysteresis of the MFD had larger loops or different types.
Literature shows that the traffic flow characteristics and hysteresis types on public holidays are significantly different from those on weekdays and weekends.It is necessary to analyze the traffic characteristics from the macroscopic perspective of networks and develop control strategies to alleviate congestion during holidays.In many studies, the heterogeneity of density distribution is perceived as a critical factor that causes the hysteresis pattern [23,24].For the heterogeneous network, Haddad et al. [15] introduced different control strategies with different levels of coordination for metropolitan transportation networks that have a hierarchical structure which consists of freeways and urban roads.Furthermore, some researchers proposed hierarchical control frameworks to investigate the dynamics of heterogeneity in urban regions, aiming at decreasing system delays and congestion heterogeneity [25,26].
In response to the issues mentioned above, this study analyzes the characteristics of regional traffic flow and the MFD on the network of Tianjin during the public holiday for Qingming Festival (1st to 4th April 2017), which is a traditional Chinese festival about ancestor worship and spring outings.The contribution of this paper is threefold.Firstly, the regional traffic flow characteristics are clarified based on the empirical data.The analysis involves different days, including the weekday before the holiday, during the holiday, and the day of Qingming Festival (the last day of the holiday).Secondly, the whole network is partitioned into several homogeneous subregions with different network densities and congestion levels.The calibrated dynamic models are applied to analyze the evolution process of traffic flow during the holidays.Finally, some control strategies are proposed to mitigate regional congestion during holidays based on the results of network state evolution under different model parameters.

Data Description
The study area is the urban network of Tianjin, as shown in Figure 1, and the data are obtained from about 400 cameras on the network within the study area.The cameras are installed at the stop lines in each direction of most intersections, covering more than 80% of the network, and the data can represent the overall traffic states.The cameras continuously take pictures of vehicles passing through the stop lines of the intersection, records the ID and GPS coordinates of the camera, the license plate number, photographing time, and driving direction.
proposed to mitigate regional congestion during holidays based on the re state evolution under different model parameters.

Data Description
The study area is the urban network of Tianjin, as shown in Figure are obtained from about 400 cameras on the network within the study ar are installed at the stop lines in each direction of most intersections, cove 80% of the network, and the data can represent the overall traffic states.T tinuously take pictures of vehicles passing through the stop lines of the i ords the ID and GPS coordinates of the camera, the license plate number, time, and driving direction.During the public holiday for Qingming Festival, there are three day The first two days are generally leisure time for residents to travel to the spring.And the last day is the Qingming Festival.On the day, residents cemeteries and worship ancestors in the morning.Several major cemeteri distributed around the city center, and the connecting roads with the c pressways or arterial roads, as shown in Figure 1.On the day of Qingm department of traffic management in Tianjin opened special bus lines from to cemeteries.The departure time was between 6:00 and 7:30, and the r around 9:00.The traffic characteristics will be very different from the ord due to multiple factors, such as the tomb-sweeping activities, spring outin the data of ordinary weekdays are also analyzed for comparisons to ex ences in traffic characteristics between holidays and ordinary weekdays.

Calculations of Network Traffic Variables
The characteristics and disciplines of the regional traffic flow are inve work traffic variables, including average flow, average speed, and avera network traffic variables are calculated using the average values over all la area.The flow is obtained directly through the original data from the cam speed and the density are obtained by calculation.
Figure 2 is a schematic illustration of the position of the cameras.Th from two adjacent cameras along a link are used to calculate the traffic var 2, a pair of neighboring intersections ( , ) i j is defined as a link with the l the number of lanes ( , ) i j n .There are three upstream cameras (straight-d During the public holiday for Qingming Festival, there are three days with no work.The first two days are generally leisure time for residents to travel to the countryside for spring.And the last day is the Qingming Festival.On the day, residents usually sweep cemeteries and worship ancestors in the morning.Several major cemeteries in Tianjin are distributed around the city center, and the connecting roads with the center are all expressways or arterial roads, as shown in Figure 1.On the day of Qingming Festival, the department of traffic management in Tianjin opened special bus lines from the city center to cemeteries.The departure time was between 6:00 and 7:30, and the return time was around 9:00.The traffic characteristics will be very different from the ordinary weekends due to multiple factors, such as the tomb-sweeping activities, spring outings, etc. Further, the data of ordinary weekdays are also analyzed for comparisons to explore the differences in traffic characteristics between holidays and ordinary weekdays.

Calculations of Network Traffic Variables
The characteristics and disciplines of the regional traffic flow are investigated by network traffic variables, including average flow, average speed, and average density.The network traffic variables are calculated using the average values over all lanes in the study area.The flow is obtained directly through the original data from the cameras, while the speed and the density are obtained by calculation.
Figure 2 is a schematic illustration of the position of the cameras.The data collected from two adjacent cameras along a link are used to calculate the traffic variables.In Figure 2, a pair of neighboring intersections (i, j) is defined as a link with the length l (i,j) and the number of lanes n (i,j) .There are three upstream cameras (straight-direction camera i c , leftdirection camera i l , and right-direction camera i r ) and one downstream camera j c .Define Ω be the set of all links, (i, j) ∈ Ω. Link flows are counted at the downstream stop line of links and aggregated every 5 min as a period T. Denote q (i,j) ((i, j) ∈ Ω) as the cumulative flow through the link (i, j) in the period T, and the network average flow q u is the mean value of all lanes covered by cameras in period T [3], given as: ( , ) (( , ) ) i j q i j  as the cumulative flow through the link ( , ) i j in the per network average flow u q is the mean value of all lanes covered by camer [3], given as: The network average density is denoted as u k , which is the ratio of t and the average density:

Network Partitioning
The related research reveals [27] that the heterogeneous distribution traffic flow has a significant impact on the MFD.Here the network is grid 1km*1km, and the density in each grid is calculated.Three days are cons weekday before the holiday (WBH), the day of Qingming Festival (DQF), a weekday (OW).Figure 3 shows the comparisons in the network density these days.The traffic densities are heterogeneously distributed on the ne congested areas are different in location and size.On the WBH, the den center is low, and the links with high densities are distributed on the su (Figure 3a).The traffic flow on the DQF is mainly on the roads to cemete rounding area (the west of the network in Figure 3b).On the OW, conges in the central area (Figure 3c).Therefore, it is necessary to divide the who subregions with similar congestion levels, which provides a basis for furth fic control.Among graph segmentation methods, the normalized cut par gorithm is commonly utilized due to its advantages effectiveness and ro research of heterogeneous networks [28,29].Therefore, the whole networ The travel time t s (i,j) for vehicle s passed link (i, j) is the time difference between the two adjacent records on the link by license plate matching.Suppose there are m T (i,j) vehicles passed through the link (i, j) in period T. The network average speed v u is calculated by the total distance L T and the total travel time t T [3], given as: The network average density is denoted as k u , which is the ratio of the average flow and the average density:

Network Partitioning
The related research reveals [27] that the heterogeneous distribution of the network traffic flow has a significant impact on the MFD.Here the network is gridded to a size of 1 km × 1 km, and the density in each grid is calculated.Three days are considered, i.e., the weekday before the holiday (WBH), the day of Qingming Festival (DQF), and the ordinary weekday (OW).Figure 3 shows the comparisons in the network density distribution of these days.The traffic densities are heterogeneously distributed on the network, and the congested areas are different in location and size.On the WBH, the density in the city center is low, and the links with high densities are distributed on the suburb of the city (Figure 3a).The traffic flow on the DQF is mainly on the roads to cemeteries and its surrounding area (the west of the network in Figure 3b).On the OW, congestion is localized in the central area (Figure 3c).Therefore, it is necessary to divide the whole network into subregions with similar congestion levels, which provides a basis for further regional traffic control.Among graph segmentation methods, the normalized cut partition (Ncut) algorithm is commonly utilized due to its advantages effectiveness and robustness in the research of heterogeneous networks [28,29].Therefore, the whole network is partitioned based on the link densities during 7:00-8:00 with congestion by Ncut algorithm.The basic idea is to partition the whole network into some subregions first and then merge them one by one; the details of the Ncut algorithm are given in Appendix A. The network densities during the peak hour are analyzed for the three days, i.e., the OW, the WBH, and the DQF. Figure 4 shows subregions by using the network partitioning on these days.The network is distinguished with different densities by partitioning.In Figure 4, the high-density area is labeled by red, the middle-density area is labeled by yellow, and the low-density area is labeled by green.The traffic density is similar within each subregion, while much distinct among different subregions.Table 1 lists the average density of each subregion and the evaluation indicators of partitioning.The network densities during the peak hour are analyzed for the three days, i.e., the OW, the WBH, and the DQF. Figure 4 shows subregions by using the network partitioning on these days.The network is distinguished with different densities by partitioning.In Figure 4, the high-density area is labeled by red, the middle-density area is labeled by yellow, and the low-density area is labeled by green.The traffic density is similar within each subregion, while much distinct among different subregions.On the WBH and the OW, the network is partitioned into three annular subregions.The regions from the inner to the outer are labeled as subregion 1, 2, and 3. On the WBH, the high-density areas are located in the suburb (subregion 3 in Figure 4a), with an average density of 15.85 veh/km/lane, much smaller than the density (41.28 veh/km/lane) in the city center on the OW.And the average density of the two inner subregions is below 10 veh/km/lane.In contrast, on the DQF, the network is partitioned into three adjacent subregions.On the DQF, the high-density area is located on the connecting roads to the cemeteries (subregion 2 in Figure 4b), and the average density of subregion 2 is more than 40 veh/km/lane.The partitioning results are consistent with the spatial distribution of the average density on the network (Figure 3).
After partitioning, the MFD on each subregion can be studied.Figures 5-7 show the MFD of each subregion on the WBH, the DQF, and the OW, respectively.In the figures, the red square scatters represent subregion 1, the blue triangle scatters ∆ represent subregion 2, and the black circles • represent subregion 3. The differences in the MFD of subregions for these days can be concluded as follows.(1) The MFDs on the WBH and the OW are different (Figures 5 and 7).Both days are weekdays, and the network is partitioned into three annular subregions.For subregions 1 and 2, the MFD on the WBH has no congested branch, while the congested branch appears on the MFD of the OW.On the WBH, the average densities of subre- (1) The MFDs on the WBH and the OW are different (Figures 5 and 7).Both days are weekdays, and the network is partitioned into three annular subregions.For subregions 1 and 2, the MFD on the WBH has no congested branch, while the congested branch appears on the MFD of the OW.On the WBH, the average densities of subregions 1 and 2 are much lower than those on the OW, reaching 50 and 37 veh/km/lane, respectively.In addition, the average speeds of all subregions on the WBH are above 16 km/h ( and ∆ in Figure 5b), while on the OW, the average speeds of subregions 1 and 2 decline significantly with the average density, and the minimal speeds are about 5.1 and 9.4 km/h (Figure 7b).On the OW, commuting trips are centralized during peak hours, resulting in congestion in the city center.The trips on the WBH are dispersed on the expressways in the suburb.Therefore, compared to that on the OW, the network average density is reduced on the WBH.(2) On the WBH (Figure 5), traffic states are all in free-flow branch with density less than 28 veh/km/lane, while a congested branch appears on the DQF (∆ in Figure 6a), where the average density approaches 50 veh/km/lane in subregion 2. The reason is that cemeteries are located in the west of the city, and a large amount of traffic flow on the DQF is squeezed into the roads to cemeteries in subregion 2.
(3) The location and duration of congestion on the DQF are different from those on the OW.On the DQF, only a small number of scatters during 7:00-8:00 are distributed in the congested branches of the MFD, which correspond to the congestion period on the network.The average density increases and speed drops significantly in subregion 2 in a short time, resulting in congestion on the network about one hour.Due to the activities regarding tomb sweeping and ancestor worship on the DQF, the high-density area is localized in the west of the city.On the OW, more scatters during peak hours are distributed on the congested branch (Figure 7a) and the congestion lasts about 2-3 h.The congested areas are located in the central of the city.
In summary, due to differences in travel characteristics, the congestion level, location, and duration on the network during the holidays vary dramatically from those on weekdays.Therefore, it is necessary to design proper regional traffic control according to the different characteristics of each subregion on the network during the holidays.

Network Dynamic Evolution and Empirical Analysis
Based on the network partitioning in Section 2, two dynamic models are calibrated to analyze the regional traffic flow evolution based on the positional relationship of the subregions.The simulation results are compared with empirical data to confirm the validity of the models, and then regional control strategies are proposed based on the calibrated model and ground truth.Figure 8 depicts the construction of the models and regional control.

The Model on the WBH
Based on the dynamic model [10,30], the data in this paper are used to cali dynamic models for different days first, and then the models are applied to com effects of different control strategies.

Model Descriptions
Based on the network partitioning in Section 2, the network for the WBH c three annular subregions, and then is built as illustrated in Figure 9.Each subre its own MFD, and the average density of each region is MFD is expressed as Equation ( 4) [10]:

The Model on the WBH
Based on the dynamic model [10,30], the data in this paper are used to calibrate the dynamic models for different days first, and then the models are applied to compare the effects of different control strategies.

Model Descriptions
Based on the network partitioning in Section 2, the network for the WBH consist of three annular subregions, and then is built as illustrated in Figure 9.Each subregion has its own MFD, and the average density of each region is k i (i = 1, 2, 3) (veh/km/lane).The MFD is expressed as Equation (4) [10]: where, v i is the free-flow speed of subregion i, k ic (veh/km/lane) is the critical density of subregion i, k b (veh/km/lane) is the jam density of the network, and w i is the backward speed of subregion i. Assume that f ij (veh/h/lane) is the transfer flow from subregion i (i = 1, 2, 3) to its neighbor subregion j, f i,in (veh/h/lane) is the external inflow of subregion i, and f i,out (veh/h/lane) is the outflow from subregion i.There are linear relationships for the variables: f i,in = P iA q ic (5b) where, P ij is the transfer ratio from subregion i to subregion j, q i = Q i (k i ) is the average flow of subregion i, P iA is the external inflow ratio of subregion i, and P iE is the outflow ratio from subregion i. q ic is the maximal average flow of subregion i.Typically, network state is defined by the average density i k of each subregion.I subregion is not completely jammed, network state is determined by the external inflo outflow, and transfer flow with other subregions.As the annular model has no extern inflow and outflow of internal subregion 1 and 2, the inflow and the outflow of subregi 3 are the total.Subregions 1 and 3 do not lie adjacent, so there is no transfer flow betwe them.That is, . The network state evolution process can be described as follows: where, i L is the average travel distance of subregion i .Substituting Equation ( 5) into Equation ( 6): Typically, network state is defined by the average density i of each subregion.If a subregion is not completely jammed, network state is determined by the external inflow, outflow, and transfer flow with other subregions.As the annular model has no external inflow and outflow of internal subregion 1 and 2, the inflow and the outflow of subregion 3 are the total.Subregions 1 and 3 do not lie adjacent, so there is no transfer flow between them.That is, The network state evolution process can be described as follows: where, L i is the average travel distance of subregion i.

Model Calibration and Empirical Analysis
The MFD of each subregion on the two days in Figures 5 and 7 is calibrated by empirical data.Generally, a congested branch will appear in the MFD where the average flow decreases as the density increases, and the density corresponding to the maximum flow q ic is the critical density k ic .For each network subregion of the WBH (Figure 5a) and subregion 3 of the OW (Figure 7a), no congested branch appears.The scatters on the MFD are all distributed within the critical density.Therefore, according to previous literature [31], the critical density of the MFD without a congested branch is calibrated as 30 veh/km/lane in this paper.The proportional relationship between the exit flow from the boundary and the flow inside the area appears to be a constant linear proportional trend after 6:00.The calibration results are listed in Table 2.The jam density k j is calculated as the average vehicle length of 5.2 m and the minimum headway distance of 1.95 m, k b = 140 veh/km/lane.The free-flow speed is the average value of the space-mean speed at night (0:00-5:00), that is v i = 20 km/h.The reason for the low free-flow speed is that the travel time included delays such as signal waiting and queuing at intersections.Moreover, the network of the study area includes different grades of roads, and the design speeds of sub-arterial roads and branches are lower than 40 km/h.There are also disturbances from non-motor vehicles and pedestrians that cause an increase in travel time.The space-mean average speed is the average of all vehicles on the network, and all these factors cause lower average speed.The calibrated model is used to simulate the evolution of the network state.Taking the average density at 7:00 of the network as the initial density, the initial densities of the network evolution are (4,6,10) veh/km/lane on the WBH, and are (24,19,4) veh/km/lane on the OW. Figure 10 shows the network evolution process over time on the WBH and the OW. Figure 10(a1) is the evolutions of subregions 1 and 2, and Figure 10(a2) is subregions 2 and 3 on the WBH, and Figure 10(b1,b2) are the evolutions on the OW.The blue line is the simulation result of the model, and the scatters are the empirical data.Points A and B represent the initial state and final state of the simulation, respectively.On the WBH, the state converges to a free-flow state in subregions 1 and 2, while it evolves to congestion in subregion 3. On the OW, both states in subregions 1 and 2 evolve to congested state over time.The simulation results are consistent with the ground truth, indicating that the model can reflect the actual evolution of the network.Compared with the OW, the traffic flow in the city center decreases on the WBH, the transfer ratio from subregion 2 to subregion 1 (P 21 ) decreases, and the congestion has not occurred in subregion 1.The transfer ratio P 23 and the external inflow ratio P 3A increase, leading to an increase in traffic flow in subregion 3, and the network becomes congested.

Model Descriptions
The network is partitioned into three adjacent subregions on the DQF.Each one connects to the external system, which means each subregion has external inflow and outflow.

  P P
), as shown in Figure 11.

The Model on the DQF 3.2.1. Model Descriptions
The network is partitioned into three adjacent subregions on the DQF.Each one connects to the external system, which means each subregion has external inflow and outflow.The transfer ratio between subregion 1 and 3 is 0 (P 13 = P 31 = 0), as shown in Figure 11.The network state evolution process can be described as follows The network state evolution process can be described as follows: Substituting Equation ( 5) into Equation ( 8):

Model Calibration and Empirical Analysis
The MFD on the DQF in Figure 6 is also calibrated by empirical data.For subregion 2, k 2c = 34.17veh/km/lane, and the scatters on the MFD of subregion 1 and 3 are distributed within the critical density, then k 1c = 30 veh/km/lane, k 3c = 30 veh/km/lane, q 1c = 364 veh/h/lane, q 2c = 450 veh/h/lane, and q 3c = 328 veh/h/lane.The calibration results of the transfer ratios, the external inflow ratio, and the outflow ratio are listed in Table 3.In simulation, the initial densities on the DQF are (4,3,2) veh/km/lane at 7:00, and the evolutions of the system over time are shown in Figure 12. Figure 12a is the evolutions of subregions 1 and 2, and Figure 12b is that of subregions 2 and 3.The blue line is the model simulation result, and the scatters are the empirical data.Points A and B represent the initial state and final state, respectively.Since the traffic flow is localized to the west of the urban network during the morning peak, the transfer ratio P 32 is significantly higher than those of other subregions.This leads to the network density of subregion 2 quickly exceeding the critical density and evolving toward congestion.The network states in subregions 1 and 3 are in free-flow state, and the model simulation results are consistent with the ground truth.
urban network during the morning peak, the transfer ratio 32 P is significantly higher than those of other subregions.leads to the network density of subregion 2 quickly exceeding the critical density and evolving toward congestion.The network states in subregions 1 and 3 are in free-flow state, and the model simulation results are consistent with the ground truth.

Regional Control Strategies Analysis
Current literature has studied traffic control measures to relieve congestion, due to large amounts of trips and heterogeneous distribution of densities during holidays.In this study, the control strategies for a congested subregion mainly include two aspects: (1) to restrict vehicles from other regions entering the congested subregion, and (2) to guide vehicles from the congested subregion to other regions.It can be achieved by restricting vehicles from entering the congested subregion during peak hours or suggesting travelers avoid peak hours and guide navigation users to use other less crowded routes in realtime.
During the Qingming Festival holiday, the traffic flow of each subregion is in a freeflow state on the WBH, while the network state of subregion 2 is congested on the DQF.Therefore, the control strategies are designed for the DQF.Based on the calibrated models, some feasible control strategies are compared by simulations, and then proper control strategies are proposed to avoid large-scale network congestion.

Control Strategies for the DQF
The peak time for cemetery-sweeping of Qingming Festival is 7:00-8:00 during the morning, and the traffic flows from the east to west on the urban network.The model calibration results indicate that the transfer ratio from subregion 3 to subregion 2 in 7:00-8:00 is 32 0.6 P  , and is 32 0.18 P  during other periods.That is to say, the amount of traffic transfer from subregion 3 to subregion 2 during peak hours is much higher than that in other periods, causing that subregion 2 evolves to congestion.Therefore, the network density should be reduced by decreasing the inflow ratio or increasing the outflow ratio of subregion 2.
On the day of Qingming Festival, residents went to cemeteries from the urban area.Most of the demand for entering sub-region 2 is rigid, so it is difficult to reduce the total inflow to sub-region 2. In practical control, residents can be reminded to travel early or delay their travel time.By extending the control time, the total demand can be dispersed over time, avoiding concentrated travel in a short period to reduce the inflow ratio to subregion 2. On the DQF, the control strategies for subregion 2 are designed as follows: (1) Strategy 1: decreasing 12 P or 2 A P , which is to decrease the transfer ratio from subregion

Regional Control Strategies Analysis
Current literature has studied traffic control measures to relieve congestion, due to large amounts of trips and heterogeneous distribution of densities during holidays.In this study, the control strategies for a congested subregion mainly include two aspects: (1) to restrict vehicles from other regions entering the congested subregion, and (2) to guide vehicles from the congested subregion to other regions.It can be achieved by restricting vehicles from entering the congested subregion during peak hours or suggesting travelers avoid peak hours and guide navigation users to use other less crowded routes in real-time.
During the Qingming Festival holiday, the traffic flow of each subregion is in a freeflow state on the WBH, while the network state of subregion 2 is congested on the DQF.Therefore, the control strategies are designed for the DQF.Based on the calibrated models, some feasible control strategies are compared by simulations, and then proper control strategies are proposed to avoid large-scale network congestion.

Control Strategies for the DQF
The peak time for cemetery-sweeping of Qingming Festival is 7:00-8:00 during the morning, and the traffic flows from the east to west on the urban network.The model calibration results indicate that the transfer ratio from subregion 3 to subregion 2 in 7:00-8:00 is P 32 = 0.6, and is P 32 = 0.18 during other periods.That is to say, the amount of traffic transfer from subregion 3 to subregion 2 during peak hours is much higher than that in other periods, causing that subregion 2 evolves to congestion.Therefore, the network density should be reduced by decreasing the inflow ratio or increasing the outflow ratio of subregion 2.
On the day of Qingming Festival, residents went to cemeteries from the urban area.Most of the demand for entering sub-region 2 is rigid, so it is difficult to reduce the total inflow to sub-region 2. In practical control, residents can be reminded to travel early or delay their travel time.By extending the control time, the total demand can be dispersed over time, avoiding concentrated travel in a short period to reduce the inflow ratio to sub-region 2. On the DQF, the control strategies for subregion 2 are designed as follows: (1) Strategy 1: decreasing P 12 or P 2A , which is to decrease the transfer ratio from subregion 1 to subregion 2 or decrease the external inflow ratio.(2) Strategy 2: increasing P 21 or P 2E , which is to increase the transfer ratio from subregion 2 to subregion 1 or increase the outflow ratio.The two control strategies limit vehicles entering subregion 2 and guide internal traffic flow away from the subregion, thus reducing the traffic density in subregion 2. The simulations under the two strategies are compared to analyze the impacts of different strategies on the network state.

Simulation Results
(1) Strategy 1: decreasing P 12 or decreasing P 2A Strategy 1 can be implemented restricting vehicles entering subregion 2 from subregion 1 or outside.In simulations, the parameters P 12 and P 2A are decreased by 5% and 10%, respectively, to compare the impacts of different control strategies.Figure 13(a1,a2) display the network state evolution of the ground truth and after decreasing P 12 .The solid blue line is the evolution with no control, and the red and green dashed lines are the evolutions when P 12 being decreased by 5% and 10%.In the same simulation period, the evolution final points are B1, B2, and B3.The average density of subregion 2 has reduced, but the average density in subregion 1 has increased.The traffic state of subregion 1 evolves to congestion under P 12 decreasing by 10% (green line in Figure 13(a1)).It indicates that this strategy can alleviate the congestion of subregion 2 but cause subregion 1 in congestion.The simulation results of network state evolution under decreasing P 2A are shown in Figure 13(b1,b2).In this strategy, traffic state of each subregion is in free-flow state under P 2A decreasing by 10% (green line in Figure 13b).Thus, among the two strategies, decreasing P 2A strategy is better., which is to increase the transfer ratio from subregion 2 to subregion 1 or increase the outflow ratio.The two control strategies limit vehicles entering subregion 2 and guide internal traffic flow away from the subregion, thus reducing the traffic density in subregion 2. The simulations under the two strategies are compared to analyze the impacts of different strategies on the network state.

Simulation Results
(1) Strategy 1: decreasing 12 P or decreasing 2 A P Strategy 1 can be implemented by restricting vehicles entering subregion 2 from subregion 1 or outside.In simulations, the parameters 12 P and 2 A P are decreased by 5% and 10%, respectively, to compare the impacts of different control strategies.Figure 13a1,a2 display the network state evolution of the ground truth and after decreasing 12 P .The solid blue line is the evolution with no control, and the red and green dashed lines are the evolutions when 12 P being decreased by 5% and 10%.In the same simulation period, the evolution final points are B1, B2, and B3.The average density of subregion 2 has reduced, but the average density in subregion 1 has increased.The traffic state of subregion 1 evolves to congestion under 12 P decreasing by 10% (green line in Figure 13a1).It indicates that this strategy can alleviate the congestion of subregion 2 but cause subregion 1 in congestion.The simulation results of network state evolution under decreasing 2 A P are shown in Figure 13b1,b2.In this strategy, traffic state of each subregion is in free-flow state under 2 A P decreasing by 10% (green line in Figure 13b).Thus, among the two strategies, decreasing 2 A P strategy is better.(2) Strategy 2: increasing P 21 or increasing P 2E For the morning peak hours on the DQF, the outflow from subregion 2 can be increased by directing vehicles away from the region in time.In simulations, the parameters P 21 and P 2E are increased by 5% and 10%, respectively.The evolutions of the network are shown in Figure 14, where the solid blue line is the evolution with no control, and the red and green dashed lines are the evolutions under the two strategies.Points B1, B2, and B3 are final points in the same simulation period.Under increasing P 21 strategy, the average density of subregion 2 reduces with the increase in P 21 , but the traffic state of subregions 1 evolves to congestion.While under increasing P 2E strategy, the average density of subregion 2 decreases, and traffic states of subregions 1 and 2 are both in free-flow.Therefore, increasing P 2E strategy is better.are shown in Figure 14, where the solid blue line is the evolution with no control, and the red and green lines are the evolutions under the two strategies.Points B1, B2, and B3 are final points in the same simulation period.Under increasing 21 P strategy, the average density of subregion 2 reduces with the increase in 21 P , but the traffic state of subregions 1 evolves to congestion.While under increasing 2 E P strategy, the average density of subregion 2 decreases, and traffic states of subregions 1 and 2 are both in free-flow.Therefore, increasing 2 E P strategy is better.In summary, when 12 P and 2 A P are decreased or 21 P and 2 E P are increased in the same proportion, controlling the external inflow or outflow ratio of subregion 2 is more effective in reducing the density of the regional network.To avoid a large amount of traffic gathering in this region in a short period, the total traffic flow can be dispersed in time, or plan detour routes for vehicles entering from outside, and implement temporary control upstream of the entrance.Since the roads connected to the cemeteries are mostly highgrade roads such as expressways and highways, the traffic flow can be limited at the entrances of expressways in actual traffic control.

Conclusions
This study analyzes the macroscopic traffic states during the Qingming Festival holiday, including the weekday before the holiday (WBH), the day of Qingming Festival (DQF), and the ordinary weekday (OW).The network densities distribution of these days In summary, when P 12 and P 2A are decreased or P 21 and P 2E are increased in the same proportion, controlling the external inflow or outflow ratio of subregion 2 is more effective in reducing the density of the regional network.To avoid a large amount of traffic gathering in this region in a short period, the total traffic flow can be dispersed in time, or plan detour routes for vehicles entering from outside, and implement temporary control upstream of the entrance.Since the roads connected to the cemeteries are mostly high-grade roads such as expressways and highways, the traffic flow can be limited at the entrances of expressways in actual traffic control.

Conclusions
This study analyzes the macroscopic traffic states during the Qingming Festival holiday, including the weekday before the holiday (WBH), the day of Qingming Festival (DQF), and the ordinary weekday (OW).The network densities distribution of these days are heterogeneous, and the congested areas are different in location and size.Based on the normalized cut (Ncut) algorithm, the network is partitioned into homogeneous subregions on different days and identifies the network congestion area on holidays.On the weekday before the holiday (WBH), the area with high density is in the suburb of the network, which is different from the ordinary weekday (OW).On the day of Qingming Festival (DQF), the congested areas distributed west on the network and subregion 2 is regional congestion during peak hours.Two dynamic models are calibrated to analyze the network evolution based on the partitioning results, and the simulation results of the calibrated model are consistent with the ground truth.
Based on the calibrated models, some proper control strategies are proposed to avoid or alleviate the regional network congestion on holidays.The results show that congestion on the DQF can be alleviated by controlling the external inflow ratio of subregion 2 or limiting the amount of traffic entering subregion 2 from the outside.Therefore, the Ncut algorithm is used to minimize Ncut(A, B), which can simultaneously maximize the total disassociation between different subregions and association within a subregion.Solving the problem is an NP-complete problem.Shi and Malik [29] proposed the equivalent eigenvalue system to partition the eigenvector corresponding to the second smallest eigenvalue.
(2) Merging and evaluation indicators for partitioning After the initial partitioning, the network is partitioned into several subregions with similar densities.A region with a higher degree of similarity is possibly divided since each subregion always is divided into two parts at each iteration step.Therefore, merging is necessary to improve the rationality of partitioning further.In the merging process, the average density differences in all two adjacent subregions are calculated.The two subregions with the smallest difference are merged until the evaluation indicators introduced in the next step are optimized.
Suppose that the whole network is partitioned into k subregions.The indicator NS k is defined to evaluate the difference between the two subregions: where N A and N B are the number of grids in subregion A and B. To make sure the highest internal similarity and low similarity with other subregions, the following self-similarity can be used: where NSN k (A, B) = min{NS k (A, X)|X ∈ Neighbor(A)},Neighbor(A) is the set of the subregions adjacent to subregion A in space.If NSS k (A) < 1, subregion A is properly partitioned.The overall partitioning can be evaluated by the value of the average similarity ANS k : where K is the set of subregions and k is the number of subregions partitioned.The optimal partitioning result is obtained when ANS k reaches the minimum.

Figure 1 .
Figure 1.Distribution of cemeteries and connecting roads in Tianjin.

Figure 1 .
Figure 1.Distribution of cemeteries and connecting roads in Tianjin.

Figure 2 .
Figure 2. Schematic illustration of the locations of cameras on a link.

Figure 4 . 1 .
Figure 4. Network subregions during peak periods on: (a) the weekday before the holid day of Qingming Festival, (c) the ordinary weekday.The numbers 1, 2 and 3 in each fig labels of the subregion, representing subregion1, subregion 2, and subregion 3, respectiv

Figure 4 .
Figure 4. Network subregions during peak periods on: (a) the weekday before the holiday, (b) the day of Qingming Festival, (c) the ordinary weekday.The numbers 1, 2 and 3 in each figure are the labels of the subregion, representing subregion1, subregion 2, and subregion 3, respectively.

21 Figure 5 .Figure 5 .
Figure 5. MFD of subregions on the weekday before the holiday: (a) average flow-average density, (b) average speed-average density.

Figure 7 .
Figure 7. MFD of subregions on the ordinary weekday: (a) average flow-average density, (b) average speed-average density.

Figure 8 .
Figure 8.The flow diagram of model construction and regional control.

Figure 8 .
Figure 8.The flow diagram of model construction and regional control.

1 Figure 9 .
Figure 9.The illustration of the model on the weekday before the holiday.

Figure 9 .
Figure 9.The illustration of the model on the weekday before the holiday.

21 Figure 10 .
Figure 10.The network evolution of: (a1) subregions 1 and 2 on the weekday before the holiday, (a2) subregions 2 and 3 on the weekday before the holiday, (b1) subregions 1 and 2 on the ordinary weekday, (b2) subregions 2 and 3 on the ordinary weekday.The blue line is the simulation result of the model, and the arrow indicates the direction of system evolution.The blue scatters are the empirical data.Points A and B represent the initial state and final state of the simulation, respectively.

Figure 10 .
Figure 10.The network evolution of: (a1) subregions 1 and 2 on the weekday before the holiday, (a2) subregions 2 and 3 on the weekday before the holiday, (b1) subregions 1 and 2 on the ordinary weekday, (b2) subregions 2 and 3 on the ordinary weekday.The blue line is the simulation result of the model, and the arrow indicates the direction of system evolution.The blue scatters are the empirical data.Points A and B represent the initial state and final state of the simulation, respectively.

Figure 11 .
Figure 11.The illustration of the model on the day of Qingming Festival.

Figure 11 .
Figure 11.The illustration of the model on the day of Qingming Festival.

Figure 12 .
Figure 12.The network evolution on the day of Qingming Festival: (a) subregions 1 and 2, (b) subregions 2 and 3.The blue line is the simulation result of the model, and the arrow indicates the direction of system evolution.The blue scatters are the empirical data.Points A and B represent the initial state and final state of the simulation, respectively.

Figure 12 .
Figure 12.The network evolution on the day of Qingming Festival: (a) subregions 1 and 2, (b) subregions 2 and 3.The blue line is the simulation result of the model, and the arrow indicates the direction of system evolution.The blue scatters are the empirical data.Points A and B represent the initial state and final state of the simulation, respectively.

Figure 13 .
Figure 13.The network evolution on the day of Qingming Festival: (a1) subregions 1 and 2, (a2) subregions 2 and 3 under 12P is decreased, (b1) subregions 1 and 2, (b2) subregions 2 and 3 under2 AP is decreased.The line is the simulation result of the model, and the arrow indicates the direction of system evolution.Point A represents the initial state and points B1, B2, and B3 represent final state of the simulation under different control strategies.

Figure 13 .
Figure 13.The network evolution on the day of Qingming Festival: (a1) subregions 1 and 2, (a2) subregions 2 and 3 under P 12 is decreased, (b1) subregions 1 and 2, (b2) subregions 2 and 3 under P 2A is decreased.The line is the simulation result of the model, and the arrow indicates the direction of system evolution.Point A represents the initial state and points B1, B2, and B3 represent final state of the simulation under different control strategies.

Figure 14 .
Figure 14.The network evolution on the day of Qingming Festival: (a1) subregions 1 and 2, (a2) subregions 2 and 3 under 21 P is increased; (b1) subregions 1 and 2, (b2) subregions 2 and 3 under 2E P is increased.The line is the simulation result of the model, and the arrow indicates the direction of system evolution.Point A represents the initial state and points B1, B2, and B3 represent final state of the simulation under different control strategies.

Figure 14 .
Figure 14.The network evolution on the day of Qingming Festival: (a1) subregions 1 and 2, (a2) subregions 2 and 3 under P 21 is increased; (b1) subregions 1 and 2, (b2) subregions 2 and 3 under P 2E is increased.The line is the simulation result of the model, and the arrow indicates the direction of system evolution.Point A represents the initial state and points B1, B2, and B3 represent final state of the simulation under different control strategies.
illustration of the locations of cameras on a link.

Table 1
lists the average density of each subregion and the evaluation indicators of partitioning.Appl.Sci.2023, 13, x FOR PEER REVIEW

Table 1 .
Evaluation indicators of partitioning (The evaluation indicators of partitioning is defined as the average similarity ANS k , it is used to the evaluate the result of partitioning and calculated by Equation (A18) in Appendix A).

Table 2 .
Parameters calibration of the model on the weekday before the holiday and the ordinary weekday.

Table 3 .
Parameters calibration of the model on the day of Qingming Festival.