Abstract

Expressway crashes in interchange areas are a critical concern in China, posing significant economic and social challenges. Utilizing three years of crash data from the Beijing–Shanghai Expressway, this study investigates the transferability and heterogeneity of crash characteristics between interchange and noninterchange areas, as well as the temporal shifts in factors influencing injury severity levels. The research employs four series of random parameters logit models to estimate three potential crash injury severity outcomes of severe injury, minor injury, and no injury (based on the most severely injured individual in each crash) and to identify key determinants, encompassing driver, vehicle, roadway, environmental, temporal, traffic, and crash attributes. Likelihood ratio tests and out-of-sample predictions are utilized to assess the temporal stability and transferability of crash area characteristics. Additionally, the marginal effects of various determinants are calculated to understand their influence across different year periods and crash types. Key variables such as overspeed, single-vehicle, AADT (annual average daily traffic volume), , and other crash type indicators are identified as significant random parameters, demonstrating heterogeneity in means and variances. Notable distinctions are observed between interchange and noninterchange crashes, indicating nontransferability, with most significant indicators revealing temporal instabilities. Furthermore, factors such as multivehicle involvement, buses, and nighttime conditions are identified as risk indicators, notably increasing the likelihood of severe injuries. These insights are invaluable for expressway designers and decision-makers, aiding in understanding the contributing mechanisms of various elements. This study suggests that stricter enforcement measures are crucial to prevent random lane changes, particularly at interchange entrances and exits. Additionally, effective management strategies and enforcement countermeasures should be implemented to mitigate crash injury outcomes in both interchange and noninterchange areas.

1. Introduction

In China, 211,074 motor vehicle crashes occurred in 2020, causing 55,950 fatalities, 214,442 injuries, and property losses of 1228.0 million CNYs [1]. In recent years, the traffic fatalities caused by expressway crashes accounted for nearly 10% of the total traffic fatalities [2], posing enormous economic and social threats. Due to the complexities of infrastructure, traffic, and driving environment, the interchange areas are the bottleneck segments of safety and efficiency in expressways. Hu et al. [3] stated that the crashes that occurred in interchange areas expressways made up 34.5% of all expressway crashes, while the proportion of lengths of interchange areas was only 9.5% from 2014 to 2016 in Jiangsu province. Therefore, it is vital to analyze the injury levels of interchange areas and adopt effective strategies to eliminate the risk levels. However, few studies addressed the injury severity outcomes of crashes that occurred in interchange areas on China expressways.

Recent studies indicated that the determinants affecting injury severity outcomes include driver, vehicle, roadway, environmental, temporal, and crash characteristics. Due to the discrete form of crash data, many discrete-data models were chosen to analyze traffic crashes, among which the logit approach is widely used [4]. However, the traditional probit/logit model assumes that the effects of multiple elements on injury levels are fixed in different crash observations, which is inconsistent with the natural situation. Moreover, traffic crashes are complicated interactions among multiple factors in various crash observations; one specific factor might influence the injury outcomes of various crashes differently. Otherwise, not all the data determining the injury outcomes could be collected, while the missing data or biased information might lead to unobserved heterogeneities [5]. To overcome the biased estimation results and errors of model structure caused by such issues, the random parameters logit models that allow the parameters to vary across observations have been shown to have better performance and utilized in many research efforts [617]. Nowadays, the random parameters logit approaches have been proven to demonstrate statistical superiority and accuracy and capture more flexibility in estimating the unobserved heterogeneity by considering heterogeneity in the means and variances [1822]. For instance, while analyzing the gender differences in injury severities of nonhelmeted motorcyclists, Wang et al. [21] validated that the random parameters logit model with heterogeneity in the means and variances outperformed other models based on tests. Moreover, Hou et al. [23] also indicated the advantages of this model in predictive performance.

Other than capturing the unobserved heterogeneities, temporal instability over year, season, or time of day might be another critical issue [24]. The potential explanation might be the changes in perception, decision, and reaction of drivers in hazardous situations along with the safety attitudes. (During the driving process, the driver’s decision-making behavior can be regarded as a continuous trade-off between the brain’s autonomic response and the controlled process. The driver’s autonomic response is instinctive. For example, when an object or vehicle suddenly appears in front of the vehicle, the driver tends to make an instinctive autonomic response. The driver’s controlled process is the behavior made after thinking for a certain period of time. For example, when he percepts the road ahead in advance, the driver could make an appropriate judgment on braking, lane change, and other behaviors. Both the autonomic response and the controlled process will change with the accumulation of driving experience and the changes in road environments. Therefore, there are different potential factors affecting the driver’s behavior in different time periods, and even when faced with the stimulation of the same factors, the corresponding behavior of the driver will change over time, ultimately resulting in determinants that affect the injury severity of highway crash and the instability of their effects on time), which mainly suffer from variations in various elements including interactions between vehicles, roadway alignment, environment, the influence of macroeconomics, and so on [24]. A growing body of research efforts has confirmed that the effects of the determinant variables change by the year [11, 16, 2529]. For instance, Dabbour et al. [30] established a year-separate model to analyze the effects of driver and vehicle characteristics on the driver injury levels of rear-end crashes in North Carolina from 2004 to 2015. Then, the year-separate model was compared with the aggregated model based on the overall sample. The estimation results show that the accuracy of the aggregated model was less than that of the year-specific model. This might be explained by a structural variation, indicating that the influence of the determinants affecting injury severity will change over the year.

Therefore, this paper intends to reveal the transferability and temporal variation of the significant factors affecting injury severities in interchange crashes (IAC) and noninterchange crashes (NIAC). The transferability of the IA and NIA crashes could make sense for better understanding the time-varying effects of determinants on expressway safety planning, management, and policy in interchange areas in China expressways. Capturing the unobserved heterogeneity, the random parameters logit models with heterogeneity in means and variances are used to reveal the injury severities based on Beijing-Shanghai Expressway crash data in 2017–2019.

Figure 1 presents the research flowchart. First, the two heterogeneity models and collected data are described, followed by the likelihood ratio tests and out-of-sample predictions. Then, the estimation results and corresponding specific recommendations are discussed. Finally, the summaries and conclusions are presented. These findings have significant implications for roadway designers and traffic management to implement appropriate and effective measurements to eliminate the risk levels of IAC and NIAC.

This study aims at (i) investigating the variations in factors determining injury severities in IAC and NIAC crashes; (ii) offering an explicit understanding of how determinants about injury severities in IAC and NIAC vary over different year periods; (iii) distinguishing and interpreting the unobserved heterogonous effects. Not only could the findings of this paper fill the knowledge gaps in the temporal instability and the determinants’ transferability among IAC and NIAC but they also serve as references for roadway designers seeking to mitigate crash occurrences and severities in interchange areas.

2. Data Collection

The crash dataset occurred in the Beijing–Shanghai Expressway being a 4-lane national expressway was collected from the traffic management department, which is designed according to the highway standard with level to rolling terrains. The cross-section width is 28.0 m, and the design speed is 120 km/h. A total of 4,238 traffic crashes were collected from 2017 to 2019. As the number of fatalities is few, the fatal and incapacitating injuries are merged into one classification: severe injury. Then, the injury severity was classified into three levels: severe injury (fatal/incapacitating injury), minor injury (incapacitating/possible injury), and no injury (property damage only) [29]. Specifically, the term “fatal injury” includes not only fatalities that occur at the scene of the accident but also those resulting from crashes within 30 days of being injured. Additionally, the term “possible injury” refers to cases where there was no visible trauma, but individuals reported experiencing pain and potential injuries as a result of the crashes. The statistical results of the severity of crashes in interchange areas and noninterchange areas over the years are shown in Table 1, and Figure 2 shows the injury severity distribution. Table 2 presents the descriptive statistics of significant variables in IAC and NIAC from 2017 to 2019 in terms of the driver, vehicle, roadway, environmental, temporal, traffic, and crash characteristics.

3. Methodology

3.1. Random Parameter Logit Approaches

Capturing heterogeneity in means and variances, the random parameters logit approaches are utilized in this paper to analyze the determinants affecting the injury severity levels of IAC and NIAC. First, the propensity function could be expressed as follows [5]:where represents the function determining the probability of injury severity outcome in crash , is a vector of explanatory variables, denotes the estimated parameter, and is a stochastic error term assumed to follow the generalized extreme-value distribution.

A standard multinomial logit model is defined when follows an assumption of extreme-value distribution [31]:where represents the probability density function of the random vector and denotes the mean or variance vector of the parameters of the probability density function.

While capturing the heterogeneity in the mean and variance, is defined as a vector of estimable parameters [20]:where and represent the vectors capturing heterogeneity in the mean and standard deviation with the corresponding parameter vector for injury severity in crash , respectively. is associated with the estimated parameters, while is a disturbance term. The and characterize the heterogeneity caused by attributes. (If the random parameters logit model shows significance in the vector of and , the model characterizes the unobserved heterogeneity in means and variances. If the model shows only significance in the vector , the model only characterizes heterogeneity in means).

The Halton sampling method is adopted to optimize efficiency and prediction performance [31]. Other than uniform, lognormal, or triangular distribution, the normal distribution is adopted in this study, which could provide a better statistical fit for random parameters logit modeling approaches [32]. After simulating different Halton draws, 500 Halton draws are used based on the maximum likelihood approach considering the trade-off among performance accuracy and estimation efficiency [33].

3.2. Likelihood Ratio Test

To examine the temporal shifts of determinants affecting the injury outcomes of IAC or NIAC across different years, along with the transferability among IAC and NIAC, three series of likelihood ratio tests (LRTs) are estimated. The first two series are used for temporal instability tests, and the third one is the transferability tests. According to the research of Behnood and Mannering [34], there are two kinds of LRTs concerning the temporal instability. The first is the pairwise test, which could propose to determine whether there are significant differences between each year’s model. Another is the global test, which is used to consider all years as a whole and determine if there are significant differences overall. Moreover, the third LRT being the transferability tests is conducted to identify overall group differences, with the group classified by the interchange area and noninterchange area crashes.

Specifically, the first pairwise test is used to compare models in two different periods and explore whether the parameter estimation is stable from one time period to another:where is the log-likelihood (LL) value of the model containing the convergence parameter of using the crash data in and is the log-likelihood value of the model using data and convergence parameters. To obtain two test results for each model comparison, the and subgroups could be reversed in this test, where the degree of freedom of the model is equal to the number of estimated parameters in the model. The LRT results of interchange and noninterchange crashes are, respectively, shown in Tables 3 and 4.

The second LRT being the global test is also used to explore whether the model parameters in two-year periods remain temporally stable, expressed as follows:where is the LL value of the convergence model for the aggregated data of crash classification (IAC and NIAC) in the use period of 2017–2019 and is the log-likelihood value of the model using the crash data and convergence parameter of year for crash classification . The degree of freedom in this test statistic is equal to the number of statistically significant parameters in the separate model minus the number of statistically significant parameters in the combined model [4].

The third series of LRT was used to judge the transferability among IAC and NIAC [14], and the corresponding test statistics are expressed as follows:where is the log-likelihood value of the convergence model combining IAC and NIAC data when period is used and and are log-likelihood values of convergence models for IAC and NIAC at year , respectively.

3.3. Out-of-Sample Prediction

In addition to the standard LRT [4], out-of-sample prediction (OOSP) is another effective method adopted in many studies to analyze the temporal instability and the nontransferability between crash types [21, 25, 28, 35, 36].

According to recent studies [23, 25], OOSP is calculated by adopting one model’s estimated parameters to predict another model’s observed crashes and explore the nontransferability among different crash models. There are two types of OOSP. The first could be used to predict the injury severity levels in another year by using IAC or NIAC in one year to reveal the temporal shifts. Another OOSP is estimated by proposing the estimated parameters of the IAC model to predict the NIAC data and reveal the difference in its prediction results. Corresponding to the estimation model, 500 Halton samples are also taken for OOSP. OOSP can be carried out by calculating the mean probability difference, and the simulation approach implemented to compute the predicted probabilities is specific as follows [21, 23]:where notations and symbols have been represented in equations (1)–(3) and N represents the number of draws used for individual observations.

4. Temporal Stability and Transferability Tests

As shown in Table 3, the LRT statistic value based on equation (4) obtained by using the estimated parameters of the 2017 IAC model using the data of the 2018 IAC is 32.18, with 10 degrees of freedom. Above all, both the two tables show that the assumption that significant parameters generated by test statistics in different years are equal can be rejected in all models, with confidence >99.00%.

Concerning the global test based on equation (5), the test statistics are 149.32 and 570.268, respectively, for IAC and NIAC, and the degrees of freedom are 8 and 30, respectively. The results also indicate that the initial hypothesis that the same significant parameters in the separate IAC and NIAC models is rejected with a confidence of >99.99%. In general, the above two LRTs indicate that crash IAC and NIAC both stay temporally unstable, which is consistent with most relevant research results [29, 33, 37].

The test statistics for 2017, 2018, and 2019 are 176.27, 162.36, and 126.81, respectively, and the degrees of freedom are 18, 17, and 12, respectively. These results also indicate that within the >99.99% confidence interval, the same initial hypothesis of weekday and weekend crash models in different years can be rejected.

Considering the temporal variability, the probability difference (mean value) between the OOSP and the original model is shown in Tables 5 and 6. The IAC and NIAC model’s parameters are respectively used to predict the crash data for the next year. Table 5 indicates that using the IAC model parameters of 2017 to predict the IAC data of 2018 can overestimate the no injury outcome probability (0.0118) and underestimate the probability of minor/severe injury outcomes (−0.0053 and −0.0065, respectively). In addition, as shown in Table 6, the probability of no injury outcome is overpredicted when using the 2017 NIAC model to predict the 2018 NIAC data. Both the two tables also reveal the potential temporal variability between the year periods.

In addition, Table 7 lists the mean values of probability difference by using NIAC model parameters to predict IAC data, while Table 8 lists the mean values of probability difference by using estimated parameters of the IAC model to predict NIAC data. Table 7 shows that the probability of no injury outcomes is overestimated, while the likelihood of minor injury levels is underestimated, with the absolute value greater than 0.01.

5. Out-of-Sample Prediction

Consistent with the results based on the LRT, the results of OOSP also indicate the nontransferability between crashes in IAC and NIAC models. The OOSP of nontransferability further illustrates the difference in the generation mechanism of IAC and NIAC. Therefore, targeted recommendations need to be implemented for these two types of the crash to eliminate the injury severity of IAC and NIAC.

The recent studies indicated high degrees of probability differences [21, 25, 36], while Yan et al. [38] showed a relatively low injury effect of temporal/crash-type shifts. And Wang et al. [21] illustrated that the contradictory phenomenon could be attributed to the overall variation in the integration of the observation size, injury outcomes distribution, and determinants. The OOSP can further verify the transferable issues across the three grades based on the two-unobserved-heterogeneity models. Notably, this can also confirm the necessity of separating the expressway crashes into IAC and NIAC.

Generally, the OOSP results also provide a basis for the nontransferability among IAC and NIAC and temporal variability across different year periods.

6. Results and Discussion

Table 9 shows the goodness-of-fit comparison results of four models, including three basic logit models (fixed parameter logit model, random parameter logit model, and random parameter logit model with mean heterogeneity) and random parameter logit model with heterogeneity in mean and variance (RPLHMV). It should be noted from Table 9 that for different individual models, the random parameter logit model with mean and variance heterogeneity obtains higher . At the same time, there are also lower AIC and BIC, specifying great improvement in the model’s fit. This also demonstrates the superiority of the RPLHMV proposed in this paper.

Then, Table 10 lists the parameter estimation results based on RPLHMV, indicating an overall good fit with the values greater than 0.65. To further study the influence of significant variables in different models on injury severity, Table 11 also demonstrates the marginal effects of determinants of IAC and NIAC models. A specific difference exists in the significant variables among IAC and NIAC models, and the influencing factors of injury severity vary in different periods. Moreover, note that the influence degree and trend of the same determinants will change across the year and by other crash areas.

6.1. Driver Characteristics

The driver’s improper operation is identified to affect minor injury in the 2017 NIAC model, and negative values showed that improper operation reduced the minor injury likelihood. For example, the marginal effect values in Table 11 showed that improper operation reduced the minor injury possibilities by 0.0078, while the no injury likelihood increased by 0.0076. However, overspeed does not show significance in all models.

6.2. Vehicle Characteristics

The single-vehicle indicator significantly affects no injury in the 2019 NIAC model and significantly influences severe injury in the 2017 IAC model. In addition, the estimation parameters of single-vehicle crashes specific to severe injury in the 2018 IAC model are identified to be statistically random. The corresponding normal distribution of N (−6.502, 4.809) means that 91.2% of single-vehicle crashes in 2018 are less likely to involve severe injuries, while another 8.8% tend to have severe injuries.

The two-vehicle crash only has a negative correlation with minor injury outcomes in 2017 NIAC. The marginal effect values in Table 11 show that this variable increases the no injury and severe injury likelihood by 0.0227 and 0.0018, respectively, and decreases the minor injury likelihood results by 0.0245. The difference in the operating speed of two vehicles might be greater, which increases greater impact force and the outcomes of crashes [39]. Moreover, multivehicle crash decreases the minor injury likelihood in the 2017 NIAC and 2018 IAC models, indicating the increased minor injury likelihood in the 2018 IAC and severe injury likelihood in the 2017 IAC. The value of marginal effect shows that this variable presents inconsistent effects on minor injuries. To be specific, the probability of minor injury likelihood decreases in the 2017 IAC model, while it increases in the other three models. However, consistent findings exist in the other two injury outcomes, such as a declined likelihood of no injuries and rising severe injury likelihood. This finding is consistent with recent studies [15, 40], indicating that the number of vehicles in the crashes tends to increase the injury severity levels.

Concerning vehicle types, minibus increases the minor injury likelihood in the 2019 NIAC model, while bus decreases the no injury likelihood in the 2018 NIAC model. In 2018 NIAC, van is positively correlated with serious injury outcomes, and marginal effects show that van improves the probability of severe injury outcomes while decreasing the no/minor injury likelihood. The truck is positively linked to minor and severe injuries in the 2018 IAC and 2017 NIAC models, respectively. The great vehicle sizes, power, and worse braking performance of trucks might produce severe outcomes of the crashes, and impose greater hazards due to the stronger crash tendency of trucks [41]. However, the marginal effect shows that truck makes different impacts on injury severity results in the two models. In contrast, truck increases the minor injury likelihood by 0.0256 in 2018 IAC and decreases the severe injury likelihood by 0.0022. In contrast, truck decreases the minor injury likelihood by 0.0010 and increases the likelihood of severe injuries by 0.0055 in the 2017 NIAC model. This may be because in interchange areas, cars tend to run at a lower speed, and the speed difference between trucks and cars is decreased, thus weakening the crash consequence. Also, the speed variation between cars and trucks in noninterchange areas is larger, increasing the injury outcomes of crashes.

In 2018, heavy trucks in noninterchange areas were positively correlated with minor injuries, and the marginal effect showed that the minor injury likelihood increased by 0.0131. The possible reason might be the fact that the maneuverability and braking performance of heavy trucks are relatively poor. In general, the car driver may be negligent and misestimate the safety distance between the car and the truck, resulting in rear-end collisions [15].

6.3. Roadway Characteristics

For roadway attributes, the significant variables that affect injury severity mainly include , , , and . Among them, is positively correlated with severe injury likelihood in the 2018 IAC model. The marginal effect shows that when the horizontal curve radius of this road section increases by 1%, the severe injury likelihood increases by 0.0533. The possible explanation is that the greater the horizontal curve radius is, the drivers are more inclined to a higher operating speed, and the probability of all kinds of crashes occurring in the interchange area is higher.

not only decreases the no injury likelihood in the 2019 NIAC model but also increases the minor injury likelihood in the 2018 NIAC. The marginal effect shows that this variable is consistent among the two models, specifically as it increases the minor injury possibilities and decreases the no injury/severe injury likelihood.

is negatively correlated with the possibility of minor injuries in noninterchange area crashes in 2017 and 2018. Marginal effects show that with a 1% increase in , the minor injury likelihood in NIAC decreased by 0.0206 and 0.0673, respectively, in 2017 and 2018, while the no/severe injury likelihood increased correspondingly. Previous studies have shown that steeper slopes indicate a shorter line of sight, so the drivers have a shorter time to react appropriately to the upcoming crash [4244]. Another explanation may be the driver’s compensation mechanism for potential risks. Due to the reduced sight distance, drivers tend to operate more cautiously and conservatively and eliminate injury severity outcomes obviously [45].

In addition, is identified as a significant random parameter on the results of severe injury in the 2019 NIAC model. Also, the normal distribution of N (−0.00463, 0.00215) indicates that 98.3% of 2019 NIAC tends to decrease the severe injury likelihood with increased and increase the severe injury likelihood for another 1.7% of observations.

6.4. Environmental Characteristics

For environmental characteristics, sunny, cloudy, rainy, snowy, and icy roads are significant in the models. Among them, all the other weather parameters are effective only in the NIAC models, except for the icy road in the 2019 IAC model. In the 2019 IAC model, icy road increases the severe injury likelihood by 0.0099.

There is a positive correlation between sunny weather and severe injury in the 2017 NIAC model, while the marginal effects demonstrate that the extreme injury likelihood increases by 0.0045 on sunny days. In addition, both cloudy and snowy days are positively relevant to minor injury outcomes in the 2018 NIAC model, and the marginal effects show that cloudy and snowy days increase the minor injury likelihood by 0.0232 and 0.0058, respectively. However, rainy weather is positively correlated with severe injury in the 2019 NIAC model, increasing the severe injury likelihood by 0.0042. Moreover, some recent studies argue more minor and severe injuries during poor weather conditions [14, 46, 47], while other research efforts indicate contradictory findings [17].

6.5. Temporal Characteristics

Monday is positively correlated with severe injuries in the 2018 NIAC model. Tuesday is negatively correlated with no injuries in the 2018 NIAC model and positively correlated with severe injuries in the 2019 IAC model. The results indicate the temporal variation of the two variables to some extent. At the same time, Wednesday shows a positive correlation with the severe injuries in NIAC for the three years, while it was not significant in the NIAC models, revealing the nontransferability of crash occurring among the two different areas. The marginal effect shows that the severe injury likelihood in NIAC increased by 0.0069 to 0.0091 on Wednesday.

As for the time of each day, the severe injury likelihood increases in both early morning and night time. The main reason may be poor visibility in the early morning and night, and easy fatigue or speeding driving, significantly increasing the probability of severe crash consequences [48]. For example, in the 2018 IAC model, the severe injury likelihood during early morning increases by 0.0819, while the probability only increases by 0.0143 in the 2018 NIAC model, indicating that IAC tends to involve more severe injuries in the early morning. At the same time, in the 2017 and 2018 NIAC models, the afternoon time decreases the minor injury likelihood, while the effects on severe injury are relatively small. It is mainly due to the good light and driving environment in the afternoon, so drivers are more likely to percept potential risk scenes and take effective operations to slow down or avoid collisions [45].

However, the influence trend of spring on the injury severity among 2017 and 2019 NIAC models is inconsistent. Specifically, the no injury likelihood increases by 0.0074 in 2017, while the likelihood decreases by 0.0083 in 2018, specifying the temporal variability of the influence trend to some extent. Summer is positively correlated with the probability of minor injuries only in the 2019 NIAC model. In winter, the severe injury likelihood increases in the 2017/2019 IAC models and 2017/2018 NIAC models, indicating more severe injuries during winter. The main explanation might be that the crashes occurred in Jiangsu Province, where the average temperature in winter is close to zero, resulting in frequent icing on the road in winter, aggravating the injury severity outcomes. Moreover, the adverse weather conditions during these months could affect the driver perception reaction time, especially under snow and rain weather [49].

6.6. Traffic Characteristics

AADT (annual average daily traffic volume) shows a normal effect across all models, declining the probability of minor and severe injuries while rising the no injury likelihood. This finding is consistent with previous studies because increasing traffic volume reduces the driving speed of vehicles [42, 45]. However, high-speed vehicles can lead to serious consequences of crashes.

Interestingly, certain differences exist in the degree of effects between IAC and NIAC models. For example, the degree of impact of AADT in the interchange area is greater than that in the noninterchange area. Specifically, when AADT increases by 1%, the decreased probability of severe injuries in IAC is 1.1–2.3 times compared to that in NIAC.

6.7. Crash Characteristic

Sideswipe type significantly affects severe and minor injuries in the 2017 and 2019 NIAC models, respectively. The marginal effect shows that this crash type decreases the no injury likelihood and increases the severe injury likelihood among both models. However, in the 2017 NIAC model, a negative correlation exists between hitting fixed objects and severe injury outcomes. The marginal effect indicated that hitting fixed objects decreases the severe injury likelihood by 0.0031. On the other hand, off-road crashes were positively correlated with minor injury results in 2019 NIAC.

Interestingly, other crash types were negatively and positively correlated with minor injuries among the 2017 NIAC and 2018 IAC models, respectively. However, the marginal effects were utterly opposite in three injury severity levels. Specifically, the severe injury likelihood in other crashes in the IAC decreased by 0.0085 in 2018, while the severe injury likelihood in NIAC increased by 0.0010 in 2017.

For the rescue response time, the rescue response time exceeding 60 minutes increases the probability of severe injury in the 2018 NIAC model, and the marginal effect value is 0.0174.

6.8. Heterogeneity in Mean and Variance

It can be found from the model results that among the crashes in the 2018 NIAC model, the parameter on Tuesday increases the mean value of the random parameter single-vehicle crash by 4.707, indicating an increased likelihood of severe injuries when single-vehicle crashes occur on Tuesday. At the same time, the mean value of the random parameter in 2019 NIAC increased by 0.00215 in the early morning. The severe injury likelihood increases in road sections with large in NIAC. This also reveals the potential relationship between the and early morning, further indicating that the severe injury likelihood increases in the early morning time.

In terms of the heterogeneity of random parameter variance, afternoon decreases the variance of the random parameter single-vehicle in the 2018 NIAC model by 0.343, indicating that, in 2018, the afternoon likelihood of single-vehicle crashes with severe injuries in noninterchange areas decreased. The results could be further evidence of better visibility and driving conditions during afternoon.

The mean and variance of random parameters only exist in the NIAC model. These heterogeneities improve the fit in goodness degree and estimation accuracy of the model to some certain extent. Then, they capture the potential influences between parameters and reveal the degree of influence of the correlation between the two parameters on the injury severity outcomes.

6.9. Practical Implications

To show the differences in the determinants, Table 12 compares the determinants showing the differences in temporal shifts and crash areas. These results have specific guiding recommendations for the safety measures and decision-making of road designers and traffic management departments, and the safety improvement strategies among interchange areas and noninterchange areas could be proposed as follows:(1)During the early morning and night time, it is recommended to set up active luminous warning messages, speed limit signs, and other reasonable measures. These measurements could prevent drivers from speeding or fatigued driving, especially in the interchange entrance and exit areas.(2)In addition, educational programs should be set up to enhance safe driving.(3)Traffic management agencies should intensify efforts to strictly strengthen the punishment and law enforcement of dangerous driving behavior.(4)Necessary exit signs should be set at appropriate locations to ensure that drivers have enough time to change lanes, such as setting solid lane change lines or oscillating lines in front of the on-ramp to control random lane-changing behaviors strictly.(5)Stricter enforcement measures should be raised to prevent random lane changes, especially at interchange entrance and exit areas.(6)For the snow and icy roads in winter, appropriate facilities such as a monitoring system, speedometer, and audio warning system should be added to remind drivers to pay attention to the slippery road surface and reduce the operating speed to improve highway safety.(7)Separate lanes should be implemented for cars and trucks to decrease the interaction among cars and trucks and eliminate the injury severity outcomes of the crashes.(8)During the roadway designing, the horizontal and vertical alignments of the coordination curve should be optimized to have continuous coordination of the road in the three-dimensional geometric parameters.

7. Conclusions

Based on crash data from the Beijing–Shanghai Expressway across three years (2017–2019), this study analyzes the differences between crashes in interchange and noninterchange areas, along with the temporal shifts in the determinants affecting the injury severity levels. Based on this study, the following conclusions are summarized:(1)By comparing with the other three basic logit models (fixed parameter logit model, the random parameter logit model, and the random parameter logit model with heterogeneity in mean), the RPLMHMV is adopted with better fit in goodness and simulation accuracy. Then, severe injury, minor injury, and no injury are considered to identify the significant variables among the driver, vehicle, roadway, environmental, temporal, traffic, and crash characteristics.(2)The LRT and OOSP demonstrate the effects of variables vary significantly across year periods and crash areas. In addition, model results and marginal effects also show that significant parameters, corresponding influence trends, and degrees are different in different individual models. For example, although the early morning and AADT have the same influence trend, specific differences exist in influence degree. Moreover, the models’ results show certain differences in the risk factors among IAC and NIAC models.(3)Future research will be conducted based on more expressways or regions to reveal the general evaluation of crash outcomes in the different areas of China. However, generally, the current finding could help expressway designers and decision-makers understand the contributing mechanism to raise proper management strategies and enforcement countermeasures to eliminate the injury severity outcomes targeted at the interchange and noninterchange areas.(4)Some limitations exit in the current study has limitations, including the missing variables concerning demographic characteristics that might lead to unobserved heterogeneity. Then, the data sample was small, for it was collected from two freeways. The crash data from more freeways could be collected to explore more accurate estimated results. Moreover, it is noted that the crash number for interchange areas is limited compared to that for noninterchange areas. More exploration would be conducted based on more detailed dataset by proposing advanced approaches.

Data Availability

The data in this paper are confidential, and the readers could apply for usage through mail based on reasonable ground.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors reviewed the results and approved the final version of the manuscript.

Acknowledgments

This study was funded by the Project of Ningxia Key Research and Development: Research on Application of Risk Identification and Intelligent Decision Making for Highway Vehicle Operation in Wind and Sand Environment (grant no. 2022BEG02017).