Abstract

Since crosswalk width and pedestrian green time directly affect the safety of signalized crosswalks, modeling an exact total crossing time model to estimate those two variables is imperative. The total crossing time needed by a group of pedestrians to cross a signalized crosswalk contains the discharge time and the crossing time. The discharge time depends primarily on the maximum queue length, which is determined by pedestrian arrival rate, red interval, waiting position distribution, and the crosswalk width. Crossing time increases when interactions between bidirectional pedestrian flows become more serious. Thus, quantifying the impacts of the start-up process on the discharge time and the effects of the interactions on the crossing time is a prerequisite for optimizing the design of signalized crosswalks. This paper establishes a modified total crossing time model consisting of modified pedestrian discharge and crossing time. Discharge time is modeled by applying traffic wave theory, and crossing time is modeled based on drag force theory. The proposed models provide guidance for the design of crosswalk width and pedestrian green intervals.

1. Introduction

An accurate total crossing time model of bidirectional pedestrian flow plays a significant role in improving pedestrian safety and crosswalk design. So, modeling the relationship between total crossing time, number of pedestrians, and crosswalk width becomes a research focus. Over the past decades, many kinds of methods have been adopted to model the total crossing time at signalized crosswalks. However, few studies have focused on the interaction process of bidirectional pedestrians. Most of the existing studies tried to explore the effect of bidirectional pedestrians on the crossing time based on several assumptions, such as the motion of pedestrians exhibiting uniform deceleration when they meet pedestrians moving in the opposite direction. However, the interaction between bidirectional pedestrians is more complex than that.

As described in the Pedestrian chapter of the Highway Capacity Manual [1], the total crossing time can be estimated by (1). In (1), 3.2 represents the initial start-up lost time, and it is contained within the discharge time in this paper. The total crossing time is a function of the bidirectional pedestrian flow, the bidirectional effect, and the crosswalk width:

Based on the method of the HCM (1), Golani and Damti put forward a modified model to estimate total crossing time [2]. They took the effect of bidirectional pedestrian flow on the total crossing time into consideration and developed a model including start-up lost time, average walking speed, and pedestrian headway. However, it was difficult to model the interaction in the proposed model. Lam et al. [3] concluded that bidirectional flow ratios have remarkable effects on the average walking speed and the maximum pedestrian capacity. However, their study ignored the impacts of different walkway widths on the average walking speed and the maximum pedestrian flow. Alhajyaseen et al. [4] investigated the effects of bidirectional flow ratios on the maximum pedestrian flow on crosswalks. Marisamynathan and Vedagiri [5] divided the total crossing delay into a waiting delay and a crossing delay, according to the crossing process, and finally modeled the total crossing delay and the total crossing time. Zhang et al. [6] evaluated pedestrian crossing time by considering the effects of crosswalk width and signal timing on crossing speed. Xu and Wang [7] used the conflict theory and the critical-gap accepted theory to model pedestrian crossing time when considering the effects of left-turning cars on pedestrian crossing. Zhao and Guo [8] built up an estimating model of crossing times. Their model is divided into a time model for the pedestrian crossing waiting area based on wave theory and a time model for the pedestrian crossing at crosswalks. A pedestrian total crossing time, equal to the sum of discharge time and crossing time, was established by Alhajyaseen et al. [911]. They proposed a macroscopic method. For modeling discharge time, the traffic wave theory was adopted. Then, they applied drag force theory to model crossing time. However, their models contain two shortcomings. Firstly, pedestrians were assumed to be uniformly distributed in the waiting area. Secondly, the interaction between bidirectional pedestrians was assumed to be fixed. Both assumptions cannot conform to reality. There is no doubt that these models are useful but they still contain some disadvantages when referring to bidirectional pedestrians.

This paper represents an improvement over the previous methods of Alhajyaseen et al. [911] and aims to improve the total crossing time model from two aspects. The purpose of this study is to establish modified discharge time model and modified crossing time. Firstly, the pedestrian waiting position distribution model, taking the crosswalk geometry and pedestrian density into consideration, is established based on actual survey data from seven survey sites. Secondly, the maximum queue length of pedestrians is calculated using the traffic wave model, and then the discharge time model of pedestrians is proposed. Thirdly, considering the nature of bidirectional pedestrian flow, a modified crossing time model based on the drag force theory is formulated as a function of the number of pedestrians, the directional ratio, and the crosswalk width. Finally, the parameters of the models are calibrated, the sensitivity of the proposed models to bidirectional pedestrian factors and crosswalk width are analyzed, and the established models are validated in a comparison with the HCM and Alhajyaseen models.

2. Data Collection and Methods

2.1. Parameters Definition

Pedestrians modify their behavior according to the traffic environment and interactions with other pedestrians. In order to ensure their safety, pedestrians tend to cross in the middle of crosswalks, and the edges of crosswalks are used less frequently. According to data collected on pedestrians’ waiting positions in the crosswalk waiting area during the red light period, the frequencies of pedestrians in different positions were analyzed. It is abundantly clear that the crosswalk geometry and volume of pedestrians have a significant influence on the waiting positions. The geometric characteristics of the crosswalk, pedestrian density, and waiting positions are defined in Figure 1.

2.1.1. Geometric Characteristics of Crosswalk

As shown in Figure 1(a), the geometric characteristics include crosswalk width and crosswalk length.

2.1.2. Waiting Positon and Average Waiting Position

The coordinate system is established in Figure 1(a), where is the origin. As shown in Figure 1(b), the coordinates of pedestrians in the waiting area are defined as the waiting position (e.g., and ). The average waiting position represents the average -coordinate across all the observed waiting pedestrians at a single survey site.

2.1.3. Waiting Pedestrian Probability

is used to represent the probability that a pedestrian will stand within a particular interval in the waiting area. represents the total number of observed waiting pedestrians at one survey site. represents the number of waiting pedestrians. For example, means that the -coordinate is located in the range and means that the -coordinate is located in the [] range. equals the ratio of to .

2.1.4. Waiting Pedestrian Density

When pedestrians enter the waiting area, they choose the appropriate waiting position according to the waiting pedestrian density. This equals the ratio of waiting pedestrians arriving at the waiting area during one green interval to the waiting area.

2.2. Study Sites

In order to study the influences of different factors on pedestrians’ waiting positions in the waiting area, seven crosswalks at four signalized intersections were chosen as the study sites. All are located in Chang Chun in China, and all are in commercial centers, ensuring sufficient pedestrian samples were obtained to enable a thorough analysis of pedestrians’ waiting behavior.

The geometric characteristics and sample sizes of pedestrians at the seven surveyed crosswalks are shown in Table 1. Series of traffic data were collected using video-processing technology at the seven crosswalks, which had different geometric characteristics and pedestrian volumes. Software based on a video image processing system was applied to extract pedestrian coordinates at every moment, and on the basis of the pedestrian coordinates, the distributions of pedestrian waiting positions in the waiting area were analyzed. Besides this, acceleration, velocity, discharge time , crossing time , and direction split ratio could all be recorded. represents the discharge time required for a pedestrian group to walk from the waiting area and step onto the crosswalk; represents the time needed to cross the crosswalk. All these data will be used to estimate parameters such as the drag coefficient ().

3. Development of Modified Discharge Time Model

This paper is a modified model of Alhajyaseen models [10], so the development process of method established by Alhajyaseen is briefly introduced. The total crossing time was divided into two parts: discharge time and crossing time . As shown in (2),where represents the total crossing time, is the discharge time needed for a pedestrian group to walk from the waiting area and step onto the crosswalk, and is the time needed for them to cross the crosswalk.

Discharge time represents the time from the first waiting pedestrian stepping into the crosswalk to the last waiting pedestrian stepping into the crosswalk. Crossing time is the necessary time to cross the crosswalk, because the discharge time is similar to the discharge time of cars, which is frequently modeled by shockwave theory. So, the shockwave theory also be used to model the discharge time. Start-up lost time is included in discharge time (Figure 2). By utilizing shockwave theory, the discharge time needed for the pedestrians in the longest queue to step into crosswalk can be expressed as where represents pedestrian arrival rate (ped/s); represents lateral distance (m); represents cycle length (s); represents pedestrian green interval (s); represents pedestrian discharge rate (ped/s); represents jam density (ped/m2); represents average speed of arriving pedestrians (m/s); represents free-flow speed (m/s); represents the crosswalk width (m).

Crossing time is determined by pedestrian speed, which is influenced by the number of opposite pedestrians and crosswalk width. Because the interaction between bidirectional pedestrians is analogous to the fluid, the drag force theory is adopted to model the crossing time . is expressed as where represents average trajectory length of the subject pedestrian demand (m); is crosswalk length (m); is adjusted drag coefficient; represents average speed of discharge pedestrians (m/s); represents the crosswalk width (m); represents the number of subject pedestrians.

However, Alhajyaseen models had several shortcomings: () Pedestrian arrival rate is assumed to be uniform. Moreover, pedestrians are assumed to arrive in a unit of “pedestrian row” per second. () The distribution of pedestrian waiting positions is assumed to be uniform. () All pedestrians occupy the same fixed space. () The deceleration of the subject pedestrian flow is assumed to be invariable. The development of their model is based on this assumption, so Newton’s Second Law is selected to model, which ignored the complex interactions between bidirectional pedestrians.

Frequently, Alhajyaseen models are proposed to define the required crosswalk width under different pedestrian demand volumes and estimate pedestrian crossing and discharge times to pedestrian signal timing optimization. However, these above shortcomings result into unreasonable computation, such as too long signal timing or too narrow pedestrian crosswalk.

Besides, this is because pedestrians prefer to obey the traffic rules and many crosswalks with directions are set to separate bidirectional pedestrians to reduce the interactions between bidirectional pedestrians’ flow in Japan. In China, pedestrians prefer to wait in the middle region of waiting area; that is, the edges of the waiting area and crosswalk have lower utilization. Besides, there are larger pedestrian volume and complex interactions between bidirectional pedestrians’ flows, and the design manual is very rough in China. All of these differences result in unreasonable crosswalk design and signal timing in China. So it is necessary to make some improvement of the method to overcome these shortcomings.

3.1. Factors Affecting Distribution of Pedestrian Waiting Positions

According to the collected pedestrian waiting positions, the pedestrian waiting frequencies in different waiting positions were modeled. When a pedestrian arrives in the waiting area, he/she will choose a waiting position according to the pedestrian density in the waiting area. To model the distribution of pedestrian waiting positions, firstly, the effects of crosswalk geometry and waiting pedestrian density are analyzed.

As shown in Figure 3(a), the relationship between crosswalk width and average waiting position is approximately linear. High values indicate a significant relationship between crosswalk width and average waiting position. The average waiting position in every waiting area increases with the crosswalk width. Similarly, Figure 3(b) shows that there is a significant relationship between crosswalk length and average waiting position. is 0.89, which indicates that the average waiting position in the waiting area increases with the crosswalk length.

In the regression analysis, the value is larger than 0.9, which indicates that the curve fit is good. As shown in Figure 3(c), the relationship between pedestrian density and average waiting position conforms to a linear relationship in one waiting area, which means that the average waiting position increases with pedestrian density. The same method is applied to analyze the relationship between pedestrian density and average waiting position at the other study sites. The minimum adjusted value at the other study sites is 0.88, which suggests that the average waiting position indeed increases with the waiting pedestrian density.

In conclusion, the crosswalk geometry and waiting pedestrian density can be selected as variables for modeling the distribution of waiting positions.

3.2. Development and Validation of Waiting Position Distribution Model

Figure 4 shows the raw distributions of pedestrian waiting positions at each waiting area of the seven observed crosswalks. As shown in Figure 4, the peaks of the histograms in the waiting areas shift to the right, which means that the waiting positions increase. To capture the different characteristics of the shapes of these waiting position distributions, an appropriate and flexible model is needed.

Walck [12] indicated that the Weibull distribution could be used to represent various random distributions through the adjustment of the shape and scale parameters. Therefore, the Weibull distribution is selected to model the waiting position distributions in the waiting area. Based on the established waiting position distributions, the discharge time model is obtained. The density probability function of the Weibull distribution is shown in where is the probability function of the Weibull distribution, is the shape parameter which controls the shape of the distribution, and is the scale parameter which controls the mean and standard deviation of the distribution.

The distributions of pedestrian waiting positions obey the Weibull distribution. According to the above analyses, the crosswalk geometry and waiting pedestrians density are the key factors influencing the waiting position distribution. Thus, they are selected as the variables for modeling the distribution, with the shape parameter and scale parameter of the Weibull distribution modeled as functions of these three variables, as shown inwhere and are the model coefficients, is the crosswalk width (m), is the crosswalk length (m), and is the density of waiting pedestrians (ped/m2).

The coefficients of the shape parameter and scale parameter are estimated by the maximum likelihood method. The resulting functions are shown in (7). The adjusted of and are 0.76 and 0.71, respectively:

Model validation is next conducted to confirm whether the established model adequately represents the distribution of waiting positions in the waiting area. In comparing the observed distributions of pedestrian waiting positions and the estimated results, Hongqi Street’s east crosswalk is selected as the validation site, because of its longer crosswalk length, width, and pedestrian volumes. The crosswalk geometry was shown in Table 1. As shown in Figure 5, based on two sample -tests for the difference in waiting positions, at a 95% confidence level there is no significant difference.

3.3. Modeling Discharge Time

The discharge time relies primarily on the maximum queue length, and then the maximum queue length is determined by pedestrian arrival rate, red interval, distribution of waiting positions, and crosswalk width. The longest queue is located in the position with the biggest waiting probability, as shown in Figure 6. The discharge time is the time needed for all the waiting pedestrians to step into the crosswalk. As analyzed in the previous section, the distribution of waiting positions of the pedestrians arriving in one-minute interval obeys the Weibull distribution. When all the pedestrians in the longest queue have stepped into the crosswalk, all the waiting pedestrians in the entire waiting area will have discharged entirely. In other words, the discharge time of the pedestrians in the longest queue represents the discharge time of all the waiting pedestrians. Thus, how to calculate the discharge time of the pedestrians in the longest queue according to the waiting positions distribution model is the key step.

Traffic wave theory is adopted for modeling the discharge time. Start-up lost time is included in the discharge time . represents the pedestrian arrival rate and is assumed to be uniform. Moreover, pedestrians are assumed to arrive in units of “pedestrian rows” per second. Thus, the pedestrian arrival rate at the position with the biggest waiting probability is , where is the biggest waiting probability and can be calculated using (5). Using traffic wave theory, the discharge time needed for the pedestrians in the longest queue to step into the crosswalk can be expressed as (9), whose derivation is shown in (8):where represents the pedestrian arrival rate (ped/s), is the biggest waiting probability in the waiting area, is the cycle length (s), is the pedestrian green interval (s), is the pedestrian discharge rate (ped/s), is the jam density (ped/m2), is the average speed of arriving pedestrians (m/s), and is the average speed of discharge of pedestrians (m/s).

4. Development of Crossing Time Model

4.1. Analysis of Speed and Acceleration

To analyze the interaction between bidirectional pedestrians, their speed and acceleration first need to be analyzed. Note that the average speeds and acceleration of pedestrians moving in one direction are recorded at every moment in every green interval. Because we need to analyze the interaction between pedestrians moving in opposing directions, their volume needs to be large enough to form “lanes.” At every study site, the same analytical method was used. However, since presenting all the analysis results in this paper would be difficult, we choose just a representative example to explain the interaction between bidirectional pedestrians. Figure 7 shows the variations in the subject pedestrians’ speeds and acceleration in one green interval at Ziyou Road’s north crosswalk. As shown in Figure 7(a), the subject pedestrians maintained the same speed before they met pedestrians moving in the opposite direction (opposing pedestrians). Then, they began to slow down to a minimum speed. After the first subject pedestrian crossed the interaction zone, all subject pedestrians began to accelerate until the last one crossed the interaction zone. The change in acceleration is a reflection of the level of interaction. Figure 7(b) shows that the deceleration of the subject pedestrians increases gradually as the interaction zone expands. That is, the interaction between bidirectional pedestrians becomes more and more significant. After the first of the focal pedestrians crosses the interaction zone, they all begin to accelerate, and the acceleration increases gradually with time, indicating that the interaction between the bidirectional pedestrians reduces.

4.2. Modeling Crossing Time

The drag force is a force that acts on an object and constrains its movement, causing resistance. When an object moves in a fluid, such as water, the resistance acting on the object is called hydrodynamic drag. When the focal and opposing pedestrians meet each other, resistance from the opposing pedestrian flow will act on the focal pedestrians. The interaction of bidirectional pedestrians is a complicated phenomenon and one that is hard to model and explain. Drag force theory is a well-developed method of modeling interaction. A drag force model was proposed by Pugh [13, 14] in which the drag force was expressed as in where is the drag force (kg m/s2), the drag coefficient (a dimensionless number), dynamic pressure (force per unit area), and cross-sectional area (m2).

The dynamic pressure can be expressed as in (11) to (14):where is the density of the fluid (kg/m3) and the speed of the subject relative to the fluid (m/s).

By substituting (11) into (10), the final drag force formula is expressed as in

Drag force is related to the shape of the subject and includes form drag and skin friction. The relation between the shape of the subject and the drag force is shown in Table 2.

The subject pedestrians are assumed to be spherical by Alhajyaseen and Nakamura [10, 11], and the drag force acting upon them is shown in Figure 9(a). This assumption has one shortcoming. The drag force will stay the same because it is determined by the cross-sectional area . As introduced in those studies, is estimated by where is the cross-sectional area of the focal pedestrian flow (m2), represents the average width of one pedestrian, and is the number of pedestrians in the focal pedestrian flow.

As shown in Figure 8, lane formation is a frequently occurring phenomenon with multiple pedestrians. The subject pedestrians pass the opposing pedestrians in a queue. In Alhajyaseen and Nakamura’s model, pedestrians interact with each other one-to-one. In fact, when lane formation occurs in bidirectional pedestrian flow, the pedestrians in one “lane” walk like a team. Thus, the shape of the subject pedestrians looks more like cylindrical, and the drag force is as shown in Figure 9(b). With cylindrical shapes, the drag force is equal to the skin friction, and the skin friction increases progressively as the contact area increases. in (12) becomes the contact area. Figure 9(b) better represents this phenomenon (where the pedestrians in one “lane” walk like a team).

As analyzed in Section 4.1, as the interaction between the subject and opposing pedestrians becomes stronger, the interference received by the subject pedestrians will become more serious. In other words, the drag force from the opposing pedestrian flow will become larger over the interaction time. The deceleration of the focal pedestrian changes over time rather than staying the same (constant deceleration). Thus, the crossing time model needs to be modified so as to reconcile with the reality.

To apply the drag force theory to model the interactions between the focal and opposing pedestrian flows, three assumptions are made:(i)Before the focal and opposing pedestrians meet each other, the speeds of both are assumed to be the free-flow speed, , so that the speed of the focal pedestrians relative to the fluid becomes (ii) in (12) represents the contact area and, because the contact area is proportional to the interaction length , it is reasonable to replace with .(iii)The opposing pedestrian flow is regarded as a fluid and the focal pedestrian flow as a moving body. Thus, the density of the fluid can be expressed as in where is the crosswalk width (m), is the number of opposing pedestrians, and is the maximum interaction length (m), as shown in Figure 10(d).

After substituting (14) and (15) into (12), the drag force formula becomeswhere is the adjusted drag coefficient, expressed as in

Figure 10 illustrates the process for a platoon of pedestrians crossing a crosswalk. In Figure 10, the time from moment (a) to moment (b) equals the discharge time, and the time from moment (a) to moment (f) equals the total crossing time.

Figure 11 shows the time-space diagram of the bidirectional pedestrians from both sides. Before the subject and opposite pedestrians meet each other, their speeds are assumed to be the free-flow speed . In order to reflect the real interaction time, the time from moment (c) to moment (e), as shown in Figure 10, is named the interaction time () here. From moment (c) to moment (d), the resulting deceleration increases with the assumed interaction time (). From moment (d) to moment (e), the resulting deceleration decreases as the assumed interaction time () increases. equals the sum of and . At moment (e), the speeds of the subject and the opposite pedestrians return to the free-flow speed .

For variable deceleration motion, the theorem of momentum and the work-energy principle are the most effective solutions. From moment (c) to moment (d), they can be expressed as in where is the mass of the group of focal pedestrians, equal to the number of focal pedestrians, is the speed of the focal pedestrians at moment (d) (m/s), is the time from moment (c) to moment (d), and is the maximum interaction length (m).

and can be expressed as in (19) by solving (18):

Similarly, from moment (d) to moment (e), the theorem of momentum and the work-energy principle can be expressed as in

and can be expressed as in (21) by solving (20):

The drag force caused by the opposing pedestrian flow is equivalent to the resistance force that causes the deceleration. By solving the equations for the crossing time , can be calculated as in (23):

4.3. Parameters Estimation

To calculate the crossing time by applying (23), the adjusted drag coefficient and maximum interaction length first need to be estimated according to observed data. 35 hours of video data from the seven surveyed crosswalks was analyzed. The bidirectional pedestrian flow in every cycle, the maximum interaction length, and the pedestrian crossing time in the same cycle were acquired. equals 1.45 m/s [10]. Then, using (23), was estimated for 187 samples. Alhajyaseen et al. [911] drew the conclusion that was proportional to . After analyzing the collected data, the relationship between and was found to be as shown in Figure 12(a). can be expressed as a function of the split ratio , as shown in (24), and can be defined as in (25):

Similarly, the maximum interaction length of the focal pedestrian platoon can be estimated by (26). The relationship between and is shown in Figure 12(b). As increases, also increases, and the coefficient of correlation (adjusted ) is 0.84.

5. Discussion and Validation

After modeling the drag coefficient and maximum interaction length under different crosswalk widths, (23) can be applied directly to estimate the crossing time for different bidirectional pedestrian demands. Figure 13(a) shows the relation between crossing time and crosswalk width. When the pedestrian split ratios take different fixed values, the crossing time decreases as the crosswalk width increases, and finally the crossing time approaches a constant. In other words, the pedestrian flow approaches free-flow when the crosswalk width is large enough. However, as the crosswalk width decreases, the crossing time becomes larger, because the interactions between the subject and opposing pedestrians become greater. When the crosswalk width is smaller than a particular value, the opposing pedestrians will definitely hinder the movement of the subject pedestrians.

To validate the established total crossing time model, the discharge times and crossing times of 187 cycles were collected at seven crosswalks. The observed total crossing times were compared with the calculated total crossing time from the established model. The differences between the observed and calculated total crossing times are shown in Figure 13(b). Based on two sample -tests for the difference in waiting positions in the waiting area, at a 95% confidence level, there is no significant difference. The total crossing time calculated using the proposed model was also compared with the total crossing time calculated using the existing model in the HCM [1], as shown in Figure 13(b). The mean absolute value of relative error of the calculated results of the proposed model was 7.38%, that of the calculated results from the Alhajyaseen model was 15.26%, and that of the calculated results from the HCM model was 12.42%. This indicates that the proposed model is more accurate and reliable. The HCM model gives a time that is less than the observed crossing time, and the cause might be that the interactions between bidirectional pedestrians were ignored. To the contrary, the interactions between bidirectional pedestrians were overestimated, and the calculated decelerations based on the Alhajyaseen’s model of the subject pedestrians were larger than the actual values. The calculated decelerations based on our model were also larger than the observed crossing times, but we tried to analyze the interaction process in detail and model different interaction stages. As we can see in Figure 13(b), when the observed crossing time is less than 30 s, the total crossing time calculated by the proposed model is consistent with the empirical data and the calculated total crossing time from the existing model in the HCM. However, when the observed crossing times are more than 30 s, the proposed model is more consistent with the observed data. In other words, the proposed model is more effective at high-occupancy crosswalks.

6. Conclusion

This paper introduces a modified method to model the discharge time and crossing time considering bidirectional pedestrian flow at signalized crosswalks. Traffic wave and drag force theories are used to model the discharge time and the crossing time, respectively. The discharge time takes the characteristics of the waiting position distribution into consideration, and the proposed model satisfies the need for the longest queue to discharge. Based on the drag force theory, the model of the crossing time is a function of the crosswalk width, the number of bidirectional pedestrians, and maximum interaction length and speed. The maximum interaction length is also modeled based on empirical data. But the effects of the traffic signal control of pedestrian on the pedestrian speed were ignored. Therefore, a continuation of this study will quantify that factor. The key step in modeling an accurate crossing time is to estimate parameters accurately. Thus, it is necessary that a large amount of data be gathered and analyzed in order that the accuracy of the drag coefficient can be improved. The proposed models provide guidance for the design of crosswalk width and pedestrian green intervals.

Conflicts of Interest

The authors declare that they have no conflicts of interest.