Next Article in Journal
Research on Evaluation and Influencing Factors of Regional Digital Innovation Ecosystem Resilience—Empirical Research Based on Panel Data of 30 Provinces and Cities in China
Next Article in Special Issue
Blending Efficiency and Resilience in the Performance Assessment of Urban Intersections: A Novel Heuristic Informed by Literature Review
Previous Article in Journal
Research on an Improved Rule-Based Energy Management Strategy Enlightened by the DP Optimization Results
Previous Article in Special Issue
Origin-Destination Traffic Survey—Case Study: Data Analyse for Bacau Municipality
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing the Performance of Highway Safety Manual (HSM) Predictive Models for Brazilian Multilane Highways

by
Olga Beatriz Barbosa Mendes
1,
Ana Paula Camargo Larocca
1,*,
Karla Rodrigues Silva
2 and
Ali Pirdavani
3,4
1
Department of Transportation Engineering (EESC-USP), Sao Carlos School of Engineering, University of Sao Paulo, Sao Carlos 13566-590, Brazil
2
Department of Transportation, RTS Administration Building, Gainesville, FL 32601, USA
3
UHasselt, Faculty of Engineering Technology, Agoralaan, 3590 Diepenbeek, Belgium
4
UHasselt, Transportation Research Institute (IMOB), Martelarenlaan 42, 3500 Hasselt, Belgium
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(13), 10474; https://doi.org/10.3390/su151310474
Submission received: 9 May 2023 / Revised: 28 June 2023 / Accepted: 29 June 2023 / Published: 3 July 2023
(This article belongs to the Special Issue Sustainable Road Transport System Planning and Optimization)

Abstract

:
This paper assesses the performance of Highway Safety Manual (HSM) predictive models when applied to Brazilian highways. The study evaluates five rural multilane highways and calculates calibration factors (Cx) of 2.62 for all types of crashes and 2.35 for Fatal or Injury (FI) crashes. The Goodness of Fit measures show that models for all types of crashes perform better than FI crashes. Additionally, the paper assesses the application of the calibrated prediction model to the atypical year of 2020, in which the COVID-19 pandemic altered traffic patterns worldwide. The HSM method was applied to 2020 using the Cx obtained from the four previous years. Results show that for 2020, the observed counts were about 10% lower than the calibrated predictive model estimate of crash frequency for all types of crashes, while the calibrated prediction of FI crashes was very close to the observed counts. The findings of this study demonstrate the usefulness of HSM predictive models in identifying high-risk areas or situations and improving road safety, contributing to making investment decisions in infrastructure and road safety more sustainable.

1. Introduction

Road safety is a global concern that has prompted nations to implement measures to reduce the fatalities and injuries resulting from road crashes. Despite some success in reducing the number of deaths in road crashes [1], the problem persists, with the proportion of fatal crashes increasing in recent years, causing more than 15 deaths per 100 thousand inhabitants yearly [2]. This number is about three times higher for emerging countries than developed countries [3], which might be related to the rise in motorization across Latin American countries that has led to a significant increase in exposure to traffic risks [4].
Therefore, countries must devise strategies to decrease this figure, including implementing stricter regulations to manage key risk factors and allocating greater resources to initiatives and studies that enhance road safety. By comprehending the factors that significantly influence the likelihood of accidents, it becomes feasible to forecast the probability of their incidence [5,6,7]. Establishing standardized definitions and methodologies for collecting comprehensive data on accidents, risk factors, and exposure occurrence is imperative to facilitate global and regional comparisons [4]. As a result, there is a lack of uniformity in the organization and collection of crash data across different regions and municipalities within the country. Each state and municipality may have its system for collecting crash data, leading to inconsistencies and challenges in data management and analysis [4,8].
Despite a decrease in the total fatalities on federal highways in Brazil over the past ten years, there has been an alarming increase in the proportion of fatal crashes [2,9]. This discrepancy may be attributed to changes in the crash reporting system since 2015, particularly the introduction of self-reporting for non-injury crashes. This could have led to an underreported number of property damage only (PDO) crashes [8,10]. Additionally, Brazil’s technological backwardness resulting from the economic and political crisis that began in 2014 may have contributed to this trend [11]. Thereby, further investigation is needed to address road safety on Brazilian highways, including investments in infrastructure and technology for accident prevention [12].
Developing effective strategies to address road safety requires a comprehensive understanding of contributing factors, which can be achieved through data-driven approaches like safety performance functions (SPF). The Highway Safety Manual (HSM) offers predictive models that integrate SPF with crash modification factors to estimate the crash frequency and identify high-risk areas and scenarios. However, it is crucial to assess the transferability of HSM predictive models when applied to an international context, particularly on Brazilian highways where data availability is limited and local SPFs are lacking. This study aims to bridge this gap by evaluating the performance of HSM predictive models on Brazilian rural multilane highways, thereby contributing to developing effective road safety strategies and advancing the United Nations Sustainable Development Goals.
Additionally, the COVID-19 pandemic had a significant impact on mobility in Brazil, leading to a reduction in the use of public transportation and an increase in individual transport [13,14]. This shift and the overall decrease in mobility during the pandemic resulted in a heterogeneous mobility pattern over time. Therefore, this study aimed to evaluate the performance of the calibrated prediction model under these atypical conditions, offering insights into its resilience and accuracy when confronted with significant changes in traffic patterns and volumes caused by the COVID-19 pandemic.

2. Literature Review

2.1. The Highway Safety Manual Predictive Model

The existing literature on crash prediction models primarily attributes crashes to inadequate driving performance with the demands of the road environment. Factors such as traffic flow, geometric attributes, road signs, and vehicle characteristics have been identified as contributing to this mismatch [15,16,17,18,19,20]. Moreover, SPFs have been developed to estimate crash rates within a specific timeframe or exposure [5,21,22,23,24,25,26]. These SPFs utilize statistical models that analyze risk indicators, including absolute numbers, frequency, and crash rates, as defined by Equation (1).
λ = N × p,
where λ is the expected crash number, N is the exposure, and p is the crash rate. The introduction of the HSM has provided a systematic approach to assessing crashes by employing analytical techniques and tools that quantify the impacts of road network planning, design, operation, and maintenance decisions. In research-based studies, the HSM has played a significant role in evaluating crashes.
The SPFs included in the HSM were developed using negative binomial (NB) regression models. These models were constructed using a generalized linear modeling (GLM) procedure, as outlined by Srinivasan et al. [27]. The SPFs consider both the infrastructure and operational characteristics.
Equation (2) illustrates how the predicted number of crashes (Npredicted) is determined using the SPF [28]. The SPF equation is specific to each facility, considering its base conditions, and adjusted by a calibration factor (Cx) and multiple crash modification factors (CMFs). Each CMF accounts for the operational and geometric characteristics (y) of the facility (x).
Npredicted = NSPFx × Cx × (CMF1x × CMF2x × … × CMFyx)
To determine Cx, Equation (3) provides the necessary calculation. The observed crashes are summed up across all sites and divided by the predicted crashes across all sites. The resulting Cx value is rounded to two decimal places and applied to the predictive model.
Cx = ∑ observed crashes/∑ predicted crashes
Calculating the corresponding Cx value for each facility type and year is advisable to customize the model. By substituting default values with locally derived values, the reliability of the predictive model can be improved. To apply this methodology, the HSM recommends a minimum desirable sample of 30 to 50 sites, representing at least 100 crashes annually [28]. Following the initial calibration, the HSM suggests utilizing the Empirical Bayes (EB) method to enhance the reliability of results and account for the regression-to-the-mean effect.
However, the model has limitations, particularly regarding its failure to consider speed limits. A study by Shirazinejad et al. [29] demonstrated that increasing the speed limit from 70 mph to 75 mph led to a significant 27% increase in total crashes and a notable 35% increase in fatal and injury crashes. Additionally, the HSM methodology fails to account for factors such as road infrastructure damage and unreasonable road design, all of which have been identified as impacting traffic safety [30].

2.2. Previous Studies on the Transferability of the HSM Model

Numerous studies have investigated the transferability and calibration of the HSM predictive model in different countries and regions. Over the past decade, researchers have explored the performance and parameters of the HSM model to assess its applicability and effectiveness in various contexts. The following studies shed light on the transferability and calibration challenges and the practical solutions and results in different countries.
Sun et al. conducted a statewide calibration of the HSM model for rural divided multilane highways in the US [31]. Their findings indicated that the HSM model reasonably predicted crashes in Missouri, with a calibration factor (Cx) of 0.98. In a study on rural two-lane roads in Arizona, Srinivasan et al. identified limitations in applying the HSM predictive models [32]. They emphasized the importance of gathering a larger sample and exploring the estimation of calibration functions to fit local data better. The overall calibration factor in this study was 1.079, indicating the success of the HSM model for US cases. D’Agostino examined the calibration factor for Italian motorways and found that the HSM model underestimated observed crash counts, with a Cx of 1.26 [33]. La Torre et al. concluded that a jurisdiction-specific base model derived from the HSM’s SPF provided a solid and reliable tool for crash prediction on the Italian freeway network [34].
In Brazilian studies, Rodrigues-Silva applied the HSM predictive model to two-lane highways in São Paulo State and found a calibration factor of 3.73 [35]. Barbosa et al. developed SPFs for intersections in Belo Horizonte, Brazil, with a calculated Cx of 2.06 [36]. Another study in Fortaleza city found a calibration factor of 0.65, highlighting the challenges in developing a nationwide SPF. Waihrich & Andrade investigated the calibration of the HSM model for multilane highways in the states of Minas Gerais and Goiás, Brazil. The resulting Cx values were 2.37 and 1.58 for each region, respectively, indicating a lack of transferability of the original HSM model in these scenarios [37]. Rodrigues-Silva compared the transferability between the HSM method and a local SPF for two-lane highways in different regions of Brazil. The calculated calibration factors were 3.67, 3.77, and 2.60 for São Paulo, Minas Gerais, and Paraná, respectively. This study highlighted the need for more parameters and knowledge in models for different facility types [38].
Studies conducted in Egypt by Elagamy et al. and in California, Maine, and Washington by Matarage & Dissanayake found that the HSM model overpredicted crash occurrences on multilane rural roads [39,40]. These studies emphasized the importance of considering local conditions and conducting calibration to improve the accuracy of predictions. Dadvar et al. proposed a method to adjust the HSM crash prediction model to provide a better fit for local data, as misallocating resources due to incorrect calibration factors can be problematic [41]. Al-Ahmadi et al. studied multilane rural highway segments in Saudi Arabia [42]. They found Cx values ranging from 0.63 to 0.78, emphasizing the need for in-depth local calibration and assessment of SPF quality. Researchers agree that the transferability of a model is dependent on the similarity of site characteristics to base conditions, and models must be built by associating regions with similar characteristics. The effectiveness of the local calibration factor as a method for transferring SPFs is widely discussed, considering socio-economic characteristics, traffic safety data distributions, and traffic flow influences on the transferability process. Kronprasert et al. compared different regression models for prediction accuracy, and the calibrated HSM SPF was the most effective model in predicting crashes on horizontal curve segments, underscoring its usefulness [43]. In a comprehensive overview, Heydari S. et al. [44] addressed road safety in low-income countries (LICs). They stressed the importance of accurate and complete road crash data for effective road safety interventions. They acknowledged that traditional sources such as police records suffer from varying levels of under-reporting, especially in LICs. They also emphasized the need to improve the quality and accuracy of road crash data through techniques like combining police and hospital records.
Countries like Brazil, with comprehensive databases integrating crash counts, traffic volume, and infrastructure data, must evaluate the performance of crash prediction models to shape investment planning strategies effectively. Conducting local calibration exercises considering regional peculiarities is crucial to enhance the transferability and precision of the HSM model across diverse countries and regions. These efforts aim to optimize the reliability of crash predictions, facilitating sustainable and informed interventions in transportation systems.

2.3. Goodness of Fit Measures

Assessing the accuracy of crash prediction models is essential in enhancing road safety measures. One approach to improve model performance is incorporating a local calibration factor (Cx) that considers the specific conditions of the target region. However, it is equally important to evaluate the model’s goodness of fit (GOF) and examine how well it aligns with observed data. In this regard, two widely used measures of forecast accuracy, the mean absolute percentage error (MAPE) and the mean absolute deviance (MAD), are commonly employed for comparative analysis.
Table 1 summarizes recent studies applying the Highway Safety Manual (HSM) method to the Brazilian road network. The table provides information on the geographical region, facility type, estimated Cx values, and the GOF tests employed in the prediction models. Various studies have focused on different types of highways, including multilane, rural, and urban roads, while also investigating the influence of road geometry and traffic characteristics on crash frequency.
Moreover, in recent research, the root mean square error (RMSE) has emerged as another evaluation metric for prediction accuracy in studies conducted by Li et al., Yao et al., and Yehia et al. [46,47,48]. However, it should be noted that the effectiveness of CURE plots in assessing model performance may be limited in studies with smaller sample sizes, as highlighted by Dadvar et al. [41].

2.4. The Impact of COVID-19 on Traffic Safety

The COVID-19 pandemic has brought about significant changes in traffic patterns and increased interest in investigating its impact on traffic safety globally. During the period of lockdowns and restrictions, there was a noticeable reduction in traffic flow in many affected countries [49]. However, studies have revealed a concerning increase in the severity of crashes during this period [50].
Research suggests that implementing nonpharmaceutical interventions (NPIs) and the higher percentage of people staying at home have had mixed effects on traffic safety. On the one hand, these measures have been associated with potential improvements in pedestrian and cyclist safety but have also increased crash risk for motor vehicle drivers [51]. Surprisingly, the average number of cyclists killed or injured per crash has tripled compared to previous years [52].
It is important to note that simply reducing traffic volume during the pandemic does not necessarily lead to improved traffic safety. This can be attributed to the homeostasis effect, wherein drivers compensate for reduced traffic by engaging in risky driving behaviors such as speeding and failure to signal [47]. Furthermore, crashes resulting in severe injuries are more likely to occur on highways due to, i.a., increased speeding, reduced law enforcement, lack of seat belt usage, and alcohol and drug abuse [49]. Therefore, effective law enforcement mechanisms should focus on preventing these behaviors [53].
Another significant pandemic effect was the shortened trip lengths and decreased travel frequency as people engaged in more online activities as an alternative to physical travel [50]. These changes in transportation characteristics and reduced traffic intensity on the roads, driven by the rise of e-commerce, have had implications for traffic patterns.
The sudden disruptions in traffic behavior caused by the pandemic offer a valuable opportunity to broaden the understanding of risk factors and the application of SPFs. As such, in this study, the calibrated HSM SPF for 2020 is compared to the crash data count in 2020 to assess its capability in evaluating the impacts of COVID-19 on the studied highways. This analysis can provide valuable insights into the effects of the pandemic on road safety and inform future strategies and interventions.

3. Materials and Methods

3.1. The HSM Crash Prediction Method for Divided Highway Segments

The required and desirable site characteristics for calibrating the SPFs for divided rural multilane roadways are described in Table 2.
The Nspf for rural multilane highways depends on the Annual Average Daily Traffic (AADT) for each year by segment and the segment length (L) in miles, as shown in Equation (4). The regression coefficients a and b are presented in the HSM (Table 3).
N s p f = e a + b × ln A A D T   + ln L
The EB method should be applied to estimate better the expected number of crashes for a single site [54], as described in Equation (5).
k = 1 e c + ln L
Here, k represents the overdispersion parameter associated with the roadway segment, L is the length of the roadway segment (in miles), and c is a regression coefficient used to determine the overdispersion of this model (see Table 3). After determining the k value for each studied segment, the Site-Specific EB Method is applied to obtain the weighted adjustment (w) placed on the predictive model estimate in Equation (6).
w = 1 1 + k × a l l   s t u d y   y e a r s N p r e d i c t e d
The final step is to obtain the Nexpected, as shown in Equation (7). This represents the final calibrated number of crashes for each segment.
N e x p e c t e d = w × N p r e d i c t e d + 1 w × N o b s e r v e d

3.2. Road Network Analysis

Five rural divided highways in São Paulo State managed by toll administration were analyzed. The selected segments are part of the highways SP-255, SP-318, SP-330, SP-334, and SP-345. The sections were chosen based on their geometric aspects and the availability of traffic volume information, as presented in Table 4. The total length of the studied roads is 235.6 km. Traffic volume data was collected through sensors strategically placed along the highways.

3.2.1. Traffic Volume Data

The traffic volume data is detected by sensor devices called “SAT” or “TESC”. The available traffic volume data were verified to match the studied highways. The average annual daily traffic (AADT) data was collected for 2016, 2017, 2018, 2019, and 2020, as presented in Table 5. There are a few cases in which there was a lack of information. For SP 318, the available AADT data corresponds to 2019 and 2020 only. As recommended by HSM, the number has been repeated for previous years (2016, 2017, 2018). For SP330_S01, the AADT for 2016 was missing, completed by linear interpolating the existing data.

3.2.2. Crash Data

The crash data analysis for the study period is presented in Table 6, while Figure 1 depicts the severity-based distribution of crash data. The findings corroborate that the observed KABC data has exhibited lower variability than PDO data since 2015. As anticipated, the number of PDO crashes has been declining since 2015, which may be attributed to regulatory changes in Brazil’s crash reporting system, resulting in underreported crashes.
Table 6 provides key information regarding the observed crash data throughout the study period, encompassing severity types and corresponding totals for each year. The mean and standard deviation values highlight the decreasing trend in crash frequencies, particularly for fatal or injury (FI) crashes. The data ranges from a minimum of zero to a maximum of 39 crashes, with decreasing means and standard deviations over the years.
Figure 1 graphically illustrates the crash frequency distribution by severity type, further emphasizing the decreasing trend in crash occurrences over the study period. Furthermore, Table 7 presents the proportion of crash data categorized by crash type, a crucial aspect for determining crash modification factors (CMFs). The table highlights the distribution of FI and PDO crashes among different collision types, including single-vehicle and multi-vehicle crashes. The proportions provide valuable insights for calculating CMFs, which serve as multiplicative factors in predicting the number of crashes based on specific road features.

3.3. Crash Modification Factor for Divided Roadway Segments (CMFs)

In Equation (2), the predicted number of crashes (Npredicted) is determined by multiplying the corresponding safety performance function (SPF) values (NSPFx) with calibration factors (Cx) and the CMFs specific to each road characteristic (CMF1x × CMF2x × … × CMFyx). The default base conditions for divided roadway segments on rural multilane highways include lane width of 12 feet, right-hand side shoulder width of 8 feet, median width of 30 feet, no lighting, and no automated speed enforcement. CMFs greater than 1.0 indicate an expected increase in crash frequencies due to specific road characteristics, while CMFs less than 1.0 signify a potential reduction in crash numbers.

4. Results and Discussion

The study applies the recommended methodology from Part C, Chapter 11 of the HSM 1st edition. The methodology includes calculating the network screening predicted frequency (Npredicted) using the crash modification factors (CMFs) and the observed crash data. A calibration factor (Cx) is determined by comparing the observed crashes (Nobserved) to the predicted crashes (Npredicted).

4.1. The Local Calibration Factor (Cx)

The results of the calibration process are presented in Table 8, which shows the observed crashes (Nobserved), predicted crashes (Npredicted), expected crashes (EB), and the corresponding calibration factors (Cx) for total and fatal injury (FI) crashes. The values of Cx indicate the similarity between the local road networks and the conditions for which the model was developed. Comparing the Cx values with previous studies, it is observed that the methodology performs closely to existing findings (Cx = 2.37 for the state of Minas Gerais) [37]. However, thoroughly examining the predicted points’ fit to the observed data is necessary to gain better insights.

4.2. The Goodness of Fit (GOF) Measures

The goodness-of-fit measures, including mean absolute deviation (MAD), mean absolute percentage error (MAPE), and root mean square error (RMSE), are presented in Table 9. The smaller values of these measures indicate a better model fit. The results suggest that the total crash model performs better than the FI crash model, indicating the variability in calibrated predicted values. Moreover, the comparison with previous studies shows a good methodology performance regarding these measures.
As anticipated, applying the Empirical Bayes (EB) method yielded estimated values closely aligned with the observed values. The MAD, MAPE, and RMSE metrics indicate that the total crash model outperforms the FI model in the final step after implementing the EB method. Conversely, the MAD and RMSE values suggest that the calibrated predicted values exhibit more substantial variation when including PDO crashes, which is reasonable given that FI crashes represent only a smaller portion (approximately 30%) of the total crash data. The MAD for FI crashes accounted for 43% of the MAD for total crashes, while the RMSE was 39%. Moreover, a study by Waihrich and Andrade reported a MAD value of 5.54 for total crashes in the Minas Gerais state, which is 24% higher than the MAD value found for total crashes in the São Paulo state (MAD = 4.44) [37].
Figure 2 displays the results for the entire study period, comparing the calibrated Npredicted values with the observed total and FI crashes to the centerline. The proximity of the model’s output to the centerline signifies a closer alignment between the predicted and observed data. The graph underscores the dispersion of total crashes in comparison to FI crashes. Furthermore, points below the centerline indicate underestimation by the model, while points above the centerline indicate an overestimation of the observed data. Approximately 56% of the predicted values fall below the centerline trend, indicating that the model has predominantly underestimated the data. Figure 3 also compares the calibrated Npredicted values for total and FI crashes to the centerline for each study year.
To further analyze the difference between Npredicted and Nexpected (results obtained after applying the EB method), Figure 4 and Figure 5 display the estimated data for total crashes, while Figure 6 and Figure 7 illustrate the estimated data for FI crashes. These graphs facilitate the estimation of the R2 values for each severity type (total or FI) and year, as presented in Table 10. Upon applying the EB method, the performance of Nexpected aligns with previous literature studies, as indicated by the proximity of the points to the centerline. This close alignment suggests that Nexpected closely resembles Nobserved. In contrast, Npredicted exhibits a moderate dispersion. The graphs for each study year exhibit a similar pattern to the one depicted in Figure 4. Notably, 2019 displays a denser distribution of Npredicted values, indicating a closer prediction of the actual number of crashes. In contrast, the estimated values for 2016 appear more dispersed, suggesting a less accurate prediction for that particular year.
By adjusting the graph scale to accommodate the smaller sample represented by FI crashes (Figure 6), a clearer understanding of the performance of Npredicted can be achieved. Consistent with previous observations, the predicted values are denser above the centerline, indicating a tendency for underprediction by the model. However, in the case of Nexpected for FI crashes, most values fall below the centerline, indicating that most expected values underestimated the observed FI crashes. This suggests that the EB method did not perform as effectively for all types of crashes. In this context, the expected crashes are more scattered, and the results indicate a general trend of underprediction.
Figure 7 demonstrates a similar performance of Nexpected compared to Figure 6, indicating that the EB method consistently underpredicted the observed crash counts. The graphs reveal a pattern where the model underestimates FI crashes in segments where more than five crashes are observed annually. Lastly, Table 10 provides the R2 estimates obtained from the developed graphs.
Table 10 presents the estimated R2 values for Npredicted and Nexpected, categorized by severity type and year. The R2 value represents the goodness of fit between observed and estimated graphs, providing a correlation measure. A higher R2 value indicates a better fit. As expected, the R2 values for Nexpected, which accounts for the observed crash counts, are significantly higher than those for Npredicted after applying the EB method. Notably, the FI crashes show lower R2 values compared to all crashes.
Table 11 presents the results of various goodness of fit (GOF) tests for calibrated predicted crashes. These tests include MAD, MAPE, RMSE, and R2. It is observed that using FI crashes yields lower MAD and RMSE values, indicating better accuracy, while MAPE and R2 show better performance for all crashes.
On the other hand, Table 12 compares the GOF parameters for expected crashes. Here, improved results are observed for all crash types. This suggests that the prediction of crashes using the HSM model performed better across all crash types.

4.3. Crash Data Analysis for 2020

Due to the COVID-19 pandemic, 2020 witnessed significant disruptions in global traffic patterns. Ongoing studies are exploring the impact of the pandemic on various aspects of human health, including traffic-related fatalities and injuries [55]. In Figure 8, the number of Property Damage Only (PDO) crashes is depicted, while Figure 9 displays the reported traffic-related fatalities on state highways in the years 2019, 2020, and 2021, as documented by the São Paulo State government [56].
Notably, there is a substantial difference in PDO crashes between April 2019 and April 2020, with a 40% reduction in PDO crashes. Similarly, the highest reduction in fatal crashes was observed during the same period, with a decline of approximately 27%. Despite the lockdown measures commencing on 22 March 2020 in São Paulo State, the impact of these measures became more pronounced in April 2020. The graphs also present the moving average of crashes over a twelve-month interval, demonstrating the reduction in crashes during that period.
Table 13 and Figure 10 illustrate the variation in fatal crashes and average AADT concerning the average and counts of the previous year. The data suggest that the decrease in fatal crashes is associated with the reduction in AADT, which can be attributed to the implementation of disease control measures during that period.
To assess the impact of COVID-19 on the analyzed segments, Table 14 presents the crash data specifically from 2020. A comparison is made between the crash counts in 2020 and the average counts from the previous four years. The data indicate a significant reduction in crashes during 2020, with a decrease of approximately 20% for all types of crashes and 11% for FI crashes compared to the average counts. This reduction in crash numbers suggests a notable influence of the COVID-19 pandemic on traffic safety, potentially due to factors such as reduced traffic volume, changes in driver behavior, and altered travel patterns resulting from pandemic-related restrictions and guidelines.
The 2020 AADT data played a crucial role in predicting the crash counts for that year using the HSM prediction model. This model relies on AADT values and segment length to estimate crash counts. In this study, the calibration factors obtained in Section 4.1 (Cx,TOTAL = 2.62 and Cx,FI = 2.35) were utilized to calculate the calibrated N predicted for 2020, as shown in Table 15. These calibration factors were derived from the baseline data of the four previous years. The EB method was also applied using the obtained parameters to calculate the Nexpected, which represents the expected crash counts considering the observed crash data. By comparing the calibrated Npredicted and Nexpected, it becomes possible to evaluate the prediction model’s performance and assess its accuracy in estimating the crash counts for 2020.
The COVID-19 pandemic has affected crash counts and changed the Average Annual Daily Traffic (AADT) values. Consequently, the Nspf (predicted crash counts) reflects the impact of the pandemic on crash frequencies. Despite this, the Nobserved (actual observed crash counts) remains approximately 10% lower than the calibrated Npredicted (predicted crash counts considering calibration factors) for all types of crashes. However, the calibrated prediction of FI crashes aligns closely with the observed counts, indicating high accuracy in predicting FI crashes.
Moreover, applying the EB method, which incorporates the observed number of crashes, brings Nexpected (expected crash counts) closer to Nobserved. This suggests that the prediction model performs well in estimating unseen data.
In Table 16, the evaluation of model performance based on mean absolute deviation (MAD) and root mean square error (RMSE) reveals that FI crashes demonstrate a better model adjustment to the actual inputs compared to all types of crashes. Conversely, the mean absolute percentage error (MAPE) and R2 (coefficient of determination) indicate that using all types of crashes yields successful model adjustment. The high MAPE value for 2020 can be attributed to the sudden reduction in crashes, which the model was not explicitly trained to anticipate as it was based on previous years’ data.
After applying the EB method and considering the 2020 data, the GOF (Goodness of Fit) tests demonstrate that the model performs better when using all types of crashes, as shown in Table 17. However, further investigation is required to fully understand the influence of infection prevention and control procedures on road safety in Brazil.

5. Conclusions

Injury crashes and fatalities have significant economic and social costs, impacting the development of a country by increasing medical expenses, insurance claims, and productivity losses. Additionally, these crashes contribute to a higher carbon footprint and infrastructure repair costs. For developing countries like Brazil, the impact is particularly significant. Predictive models, such as the HSM prediction model, can play a crucial role in identifying high-risk areas and situations, enabling targeted interventions to improve road safety, and contribute to a more sustainable transport system. This aligns with the objectives of the 2030 Agenda, which aims to ensure sustainable transport systems that promote economic growth, social inclusion, and environmental sustainability, including reducing the number of road traffic deaths and injuries by 50% by 2030.
The assessment of the HSM prediction model employment during the atypical year of 2020, marked by the COVID-19 pandemic and the resulting changes in traffic patterns, helped understand the temporal transferability of the model. The calibrated prediction model showed promising results, although there was a slight underestimation of crash counts for all types of crashes compared to the observed values. However, the calibrated prediction of fatal and injury crashes (FI crashes) closely matched the observed counts, demonstrating the model’s capability of capturing severe crash fluctuations.
Using all types of crashes in the model yielded better results in most goodness-of-fit tests, indicating that underreporting crashes did not significantly affect the model’s validity. However, it is essential to acknowledge that additional risk factors not accounted for by the Safety Performance Functions (SPFs) may influence road safety outcomes. The study also highlighted the importance of calibration to local conditions and the need to establish “good-enough” thresholds for other contexts, as Brazilian data and road characteristics may differ from those in the original development of the HSM SPFs.
This research provides valuable insights into applying the HSM prediction model for multilane rural highways in Brazil, serving as a reference for safety assessment and guidance for highway administration, municipalities, and toll agencies. It demonstrates the need for calibration and suggests further investigation into roadway characteristics, driver behavior, and crash patterns specific to the Brazilian context to enhance the accuracy of crash predictions. The findings also contribute to understanding SPF transferability and the model’s performance in atypical years.
Future studies should consider using different calibration methods and functions and explore additional goodness-of-fit tests such as cure plots, chi-square, and coefficient of variation. Developing jurisdiction-specific SPFs for local conditions and addressing questions related to SPF calibration frequency and temporal transferability acceptable thresholds would further enhance the knowledge in this research area.

Author Contributions

Conceptualization, A.P.C.L., A.P., K.R.S. and O.B.B.M.; methodology, O.B.B.M.; software, O.B.B.M.; validation, A.P. and K.R.S.; formal analysis, O.B.B.M.; investigation, O.B.B.M.; resources, A.P.C.L.; data curation, A.P. and K.R.S.; writing—original draft preparation, O.B.B.M. and A.P.C.L.; writing—review and editing, A.P. and K.R.S.; supervision, A.P.C.L.; project administration, A.P.C.L.; funding acquisition, A.P.C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Coordinate Improvement of University Personnel (CAPES code 001); Conselho Nacional de Pesquisa e Desenvolvimento: 407056/2022-0 and São Paulo Research Foundation: 2021/10727-1.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Andrade, F.R.; Antunes, J.L.F. Trends in the Number of Traffic Accident Victims on Brazil’s Federal Highways before and after the Start of the Decade of Action for Road Safety. Rep. Public Health 2019, 35, e00250218. [Google Scholar] [CrossRef]
  2. DATASUS. Observatório Nacional de Segurança Viária Home Page. Available online: https://www.onsv.org.br/ (accessed on 8 December 2020).
  3. World Health Organization. Global Status Report on Road Safety 2018; World Health Organization: Geneva, Switzerland, 2018; p. 424. [Google Scholar]
  4. International Transport Forum. Benchmarking Road Safety in Latin America-Case-Specific Policy Analysis; ITF: Paris, France, 2017. [Google Scholar]
  5. Hauer, E. Statistical Road Safety Modeling. Transp. Res. Rec. J. Transp. Res. Board 2004, 1897, 81–87. [Google Scholar] [CrossRef]
  6. Hauer, E. The Art of Regression Modeling in Road Safety; Springer International Publishing: Cham, Switzerland, 2015; ISBN 978-3-319-12528-2. [Google Scholar]
  7. Hauer, E.; Bonneson, J.A.; Council, F.; Srinivasan, R.; Zegeer, C. Crash Modification Factors. Transp. Res. Rec. 2012, 2279, 67–74. [Google Scholar] [CrossRef]
  8. Bhalla, K.; Job, S.; Mitra, S.; E Harrison, J.; Mbugua, L.W.; Neki, K.; Gutierrez, H.; Balasubramaniyan, R.; Winer, M.; Vos, T.; et al. Assessing discrepancies in estimates of road traffic deaths in Brazil. Inj. Prev. 2023, 2, 044871. [Google Scholar] [CrossRef]
  9. Brazil Ministry of Health. Portaria No 184, de 24 de Junho de 2010; Brazil Ministry of Health: Rio de Janeiro, Brazil, 2010.
  10. Brazilian Federal Highway Police PRF Disponibiliza Servico via Internet Para Registro de Acidentes Sem Vitima. Available online: https://legado.justica.gov.br/news/prf-disponibiliza-servico-via-internet-para-registro-de-acidentes-sem-vitimas (accessed on 24 March 2022).
  11. Bernal-Meza, R. The Rise and Decline of Brazil as a Regional Power (2003–2016). Lat. Am. Perspect. 2022, 49, 51–67. [Google Scholar] [CrossRef]
  12. Brazil Ministry of Justice and Public Security PRF Disponibiliza Serviço via Internet Para Registro de Acidentes Sem Vítimas. Available online: https://www.justica.gov.br/news/prf-disponibiliza-servico-via-internet-para-registro-de-acidentes-sem-vitimas (accessed on 20 March 2021).
  13. Costa, C.S. Análise Da Substituição Do Transporte Público Pelo Serviço de Ridesourcing Durante a Pandemia Da COVID-19 No Brasil. Master’s Thesis, Universidade de São Paulo, Escola de Engenharia de São Carlos, São Carlos, Brazil, 2022. [Google Scholar]
  14. Barbieri, D.M.; Lou, B.; Passavanti, M.; Hui, C.; Lessa, D.A.; Maharaj, B.; Banerjee, A.; Wang, F.; Chang, K.; Naik, B.; et al. A survey dataset to evaluate the changes in mobility and transportation due to COVID-19 travel restrictions in Australia, Brazil, China, Ghana, India, Iran, Italy, Norway, South Africa, United States. Data Brief 2020, 33, 106459. [Google Scholar] [CrossRef]
  15. Aljanahi, A.A.M.; Rhodes, A.H.; Metcalfe, A.V. Speed, Speed Limits and Road Traffic Accidents under Free Flow Conditions. Accid. Anal. Prev. 1999, 31, 161–168. [Google Scholar] [CrossRef]
  16. Chang, L.-Y.; Mannering, F. Analysis of Injury Severity and Vehicle Occupancy in Truck-and Non-Truck-Involved Accidents. Accid. Anal. Prev. 1999, 31, 579–592. [Google Scholar] [CrossRef]
  17. Geurts, K.; Wets, G. Black Spot Analysis Methods: Literature Review. Onderz. Kennis Verkeersonveiligheid 2003, 1, 32. [Google Scholar]
  18. Hydén, C.; Várhelyi, A. The Effects on Safety, Time Consumption and Environment of Large Scale Use of Roundabouts in an Urban Area: A Case Study. Accid. Anal. Prev. 2000, 32, 11–23. [Google Scholar] [CrossRef]
  19. Ivan, J.N.; Pasupathy, R.K.; Ossenbruggen, P.J. Differences in Causality Factors for Single and Multi-Vehicle Crashes on Two-Lane Roads. Accid. Anal. Prev. 1999, 31, 695–704. [Google Scholar] [CrossRef]
  20. Wright, C.; Abbess, C.; Jarrett, D. Estimating the regression-to-mean effect associated with road accident black spot treatment: Towards a more realistic approach. Accid. Anal. Prev. 1986, 20, 199–214. [Google Scholar] [CrossRef] [PubMed]
  21. Abdel-Aty, M.A.; Radwan, A. Modeling traffic accident occurrence and involvement. Accid. Anal. Prev. 2000, 32, 633–642. [Google Scholar] [CrossRef] [PubMed]
  22. Buyco, C.; Saccomanno, F.F. Analysis of truck accident rates using loglinear models. Can. J. Civ. Eng. 1988, 15, 397–408. [Google Scholar] [CrossRef]
  23. Charnes, A.; Cooper, W.; Seiford, L.; Stutz, J. A multiplicative model for efficiency analysis. Socio-Econ. Plan. Sci. 1982, 16, 223–224. [Google Scholar] [CrossRef]
  24. Lord, D.; Miranda-Moreno, L.F. Effects of low sample mean values and small sample size on the estimation of the fixed dispersion parameter of Poisson-gamma models for modeling motor vehicle crashes: A Bayesian perspective. Saf. Sci. 2008, 46, 751–770. [Google Scholar] [CrossRef]
  25. Oppe, S. The use of multiplicative models for analysis of road safety data. Accid. Anal. Prev. 1979, 11, 101–115. [Google Scholar] [CrossRef]
  26. Quddus, M.A. Time series count data models: An empirical application to traffic accidents. Accid. Anal. Prev. 2008, 40, 1732–1741. [Google Scholar] [CrossRef] [Green Version]
  27. Srinivasan, R.; Carter, D.; Bauer, K. Safety Performance Function Decision Guide: SPF Calibration Versus SPF Development; US Department of Transportation: Washington, DC, USA, 2013; pp. 1–31.
  28. American Association of State Highway and Transportation Officials (AASHTO). Highway Safety Manual; American Association of State Highway and Transportation Officials (AASHTO): Washington, DC, USA, 2010; Volume 1, ISBN 978-1-56051-477-0. [Google Scholar]
  29. Shirazinejad, R.S.; Dissanayake, S.; Al-Bayati, A.J.; York, D.D. Evaluating the Safety Impacts of Increased Speed Limits on Freeways in Kansas Using Before-And-After Study Approach. Sustainability 2019, 11, 119. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, Z.; Wu, J.; Yousaf, A.; McIlroy, R.C.; Wang, L.; Liu, M.; Plant, K.L.; Stanton, N.A. A Study of Vulnerable Road Users’ Behaviors Using Schema Theory and the Perceptual Cycle Model. Sustainability 2023, 15, 8339. [Google Scholar] [CrossRef]
  31. Sun, C.; Brown, H.; Edara, P.; Carlos, B.; Nam, K. Calibration of the Highway Safety Manual for Missouri; Final Reports & Technical Briefs; Mid-America Transportation Center: Lincoln, NE, USA, 2013; pp. 1–240, Paper 94. [Google Scholar]
  32. Srinivasan, R.; Colety, M.; Bahar, G.; Crowther, B.; Farmen, M. Estimation of Calibration Functions for Predicting Crashes on Rural Two-Lane Roads in Arizona. Transp. Res. Rec. J. Transp. Res. Board 2016, 2583, 17–24. [Google Scholar] [CrossRef]
  33. D’agostino, C. Investigating Transferability and Goodness of Fit of Two Different Approaches of Segmentation and Model form for estimating safety performance of Motorways. Procedia Eng. 2014, 84, 613–623. [Google Scholar] [CrossRef] [Green Version]
  34. La Torre, F.; Meocci, M.; Domenichini, L.; Branzi, V.; Tanzi, N.; Paliotto, A. Development of an accident prediction model for Italian freeways. Accid. Anal. Prev. 2019, 124, 1–11. [Google Scholar] [CrossRef] [PubMed]
  35. Silva, K.C.R. Aplicação Do Modelo de Previsão de Acidentes Do HSM Em Rodovias de Pista Simples Do Estado de São Paulo; Universidade de São Paulo: São Carlos, Brazil, 2012. [Google Scholar]
  36. Barbosa, H.; Cunto, F.; Bezerra, B.; Nodari, C.; Jacques, M.A. Safety performance models for urban intersections in Brazil. Accid. Anal. Prev. 2014, 70, 258–266. [Google Scholar] [CrossRef] [PubMed]
  37. Waihrich, D.R.L.d.S.; Andrade, M. Calibração Do Método de Previsão de Acidentes Do Highway Safety Manual (HSM) Para Trechos Rodoviários de Pista Dupla No Brasil. In Proceedings of the Xxix Congresso Nacional de Pesquisa em Transporte da Anpet 2015, Ouro Preto, Brazil, 9–13 November 2015; pp. 1434–1437. [Google Scholar]
  38. Silva, K.C.R. Assessing the Transferability of Crash Prediction Models for Two Lane Highways in Brazil; Universidade de São Paulo: São Carlos, Brazil, 2017. [Google Scholar]
  39. Elagamy, S.R.; El-Badawy, S.M.; Shwaly, S.A.; Zidan, Z.M.; Shahdah, U.E. Segmentation Effect on the Transferability of International Safety Performance Functions for Rural Roads in Egypt. Safety 2020, 6, 43. [Google Scholar] [CrossRef]
  40. Matarage, I.C.; Dissanayake, S. Quality assessment between calibrated highway safety manual safety performance functions and calibration functions for predicting crashes on freeway facilities. J. Traffic Transp. Eng. 2020, 7, 76–87. [Google Scholar] [CrossRef]
  41. Dadvar, S.; Lee, Y.-J.; Shin, H.-S. Improving crash predictability of the Highway Safety Manual through optimizing local calibration process. Accid. Anal. Prev. 2020, 136, 105393. [Google Scholar] [CrossRef]
  42. Al-Ahmadi, H.M.; Jamal, A.; Ahmed, T.; Rahman, M.T.; Reza, I.; Farooq, D. Calibrating the Highway Safety Manual Predictive Models for Multilane Rural Highway Segments in Saudi Arabia. Arab. J. Sci. Eng. 2021, 46, 11471–11485. [Google Scholar] [CrossRef]
  43. Kronprasert, N.; Boontan, K.; Kanha, P. Crash Prediction Models for Horizontal Curve Segments on Two-Lane Rural Roads in Thailand. Sustainability 2021, 13, 9011. [Google Scholar] [CrossRef]
  44. Heydari, S.; Hickford, A.; McIlroy, R.; Turner, J.; Bachani, A.M. Road Safety in Low-Income Countries: State of Knowledge and Future Directions. Sustainability 2019, 11, 6249. [Google Scholar] [CrossRef] [Green Version]
  45. Cunto, F.; Sobreira, L.; Ferreira, S. Assessing the Transferability of the Highway Safety Manual Predictive Method for Urban Roads in Fortaleza City, Brazil. J. Transp. Eng. 2015, 141, 734. [Google Scholar] [CrossRef]
  46. Li, L.; Gayah, V.V.; Donnell, E.T. Development of regionalized SPFs for two-lane rural roads in Pennsylvania. Accid. Anal. Prev. 2017, 108, 343–353. [Google Scholar] [CrossRef]
  47. Yao, Y.; Geara, T.G.; Shi, W. Impact of COVID-19 on city-scale transportation and safety: An early experience from Detroit. Smart Health 2021, 22, 100218. [Google Scholar] [CrossRef]
  48. Yehia, A.; Wang, X.; Feng, M.; Yang, X.; Gong, J.; Zhu, Z. Applicability of boosting techniques in calibrating safety performance functions for freeways. Accid. Anal. Prev. 2021, 159, 106193. [Google Scholar] [CrossRef]
  49. Yasin, Y.J.; Grivna, M.; Abu-Zidan, F.M. Global impact of COVID-19 pandemic on road traffic collisions. World J. Emerg. Surg. 2021, 16, 51. [Google Scholar] [CrossRef]
  50. Shaik, E.; Ahmed, S. An overview of the impact of COVID-19 on road traffic safety and travel behavior. Transp. Eng. 2022, 9, 100119. [Google Scholar] [CrossRef]
  51. Dong, N.; Zhang, J.; Liu, X.; Xu, P.; Wu, Y.; Wu, H. Association of human mobility with road crashes for pandemic-ready safer mobility: A New York City case study. Accid. Anal. Prev. 2022, 165, 106478. [Google Scholar] [CrossRef]
  52. Li, J.; Zhao, Z. Impact of COVID-19 travel-restriction policies on road traffic accident patterns with emphasis on cyclists: A case study of New York City. Accid. Anal. Prev. 2022, 167, 106586. [Google Scholar] [CrossRef]
  53. Vanlaar, W.; Woods-Fry, H.; Barrett, H.; Lyon, C.; Brown, S.; Wicklund, C.; Robertson, R. The impact of COVID-19 on road safety in Canada and the United States. Accid. Anal. Prev. 2021, 160, 106324. [Google Scholar] [CrossRef]
  54. Elvik, R. State-of-the-Art Approaches to Road Accident Black Spot Management and Safety Analysis of Road Networks; Transportøkonomisk Institutt: Oslo, Norway, 2007; p. 126. [Google Scholar]
  55. Barnes, S.R.; Beland, L.; Huh, J.; Kim, D. COVID-19 lockdown and traffic accidents: Lessons from the pandemic. Contemp. Econ. Policy 2022, 40, 349–368. [Google Scholar] [CrossRef]
  56. Sao Paulo State Government Statewise System for Data Analysis Foundation State-Level Data on a Variety of COVID-19 Metrics Home Page. Available online: http://200.144.30.103:8081/vdm/Page/Index.aspx (accessed on 18 September 2022).
Figure 1. Crash frequency by severity.
Figure 1. Crash frequency by severity.
Sustainability 15 10474 g001
Figure 2. The correlation of calibrated Npredicted versus Nobserved comparing total and FI crashes for the total period of study.
Figure 2. The correlation of calibrated Npredicted versus Nobserved comparing total and FI crashes for the total period of study.
Sustainability 15 10474 g002
Figure 3. Comparison between Npredicted versus Nobserved for total and FI crashes for 2016, 2017, 2018, and 2019, respectively.
Figure 3. Comparison between Npredicted versus Nobserved for total and FI crashes for 2016, 2017, 2018, and 2019, respectively.
Sustainability 15 10474 g003
Figure 4. The correlation between the observed crash data and the estimated calibrated Npredicted and Nexpected for all crashes for all study years.
Figure 4. The correlation between the observed crash data and the estimated calibrated Npredicted and Nexpected for all crashes for all study years.
Sustainability 15 10474 g004
Figure 5. The correlation between the observed crashes and the estimated calibrated Npredicted and Nexpected for all types of crashes in 2016, 2017, 2018, and 2019.
Figure 5. The correlation between the observed crashes and the estimated calibrated Npredicted and Nexpected for all types of crashes in 2016, 2017, 2018, and 2019.
Sustainability 15 10474 g005
Figure 6. The correlation between the observed crashes and the estimated calibrated Npredicted and Nexpected for FI crashes at the total period of study.
Figure 6. The correlation between the observed crashes and the estimated calibrated Npredicted and Nexpected for FI crashes at the total period of study.
Sustainability 15 10474 g006
Figure 7. The correlation between the observed crashes and the estimated calibrated Npredicted and Nexpected for FI crashes in 2016, 2017, 2018, and 2019.
Figure 7. The correlation between the observed crashes and the estimated calibrated Npredicted and Nexpected for FI crashes in 2016, 2017, 2018, and 2019.
Sustainability 15 10474 g007
Figure 8. PDO crashes on state highways in 2019, 2020, and 2021.
Figure 8. PDO crashes on state highways in 2019, 2020, and 2021.
Sustainability 15 10474 g008
Figure 9. Fatal crashes on state highways in 2019, 2020, and 2021.
Figure 9. Fatal crashes on state highways in 2019, 2020, and 2021.
Sustainability 15 10474 g009
Figure 10. Estimated variation of crashes and AADT for state highways.
Figure 10. Estimated variation of crashes and AADT for state highways.
Sustainability 15 10474 g010
Table 1. Works of HSM method application in Brazil [35,36,37,38,45].
Table 1. Works of HSM method application in Brazil [35,36,37,38,45].
AuthorRegionFacility TypeCxGOF
Rodrigues-Silva (2012) [35]SPTwo-lane Rural Highways3.73Chi square test and Kolmogorov-Smirnov
Barbosa et al. (2014) [36]CEUrban Intersection0.65AIC, R2 statistic, and CURE plots
2.06
Cunto et al. (2015) [45]Fortaleza (CE)Urban Roads0.98 MAD ,   MAPE ,   CURE ,   Pearson   p 2 statistics and z-score
2.15
Waihrich & Andrade (2015) [37]MGMultilane Rural Highways2.37MAD, MAPE and R2Efron
GO/DF1.58
Rodrigues-Silva (2017) [38]SPTwo-lane Rural Highways3.67MAD, MAPE, R2Efron, and CURE plots
PR3.77
MG2.60
Table 2. Data needed to calibrate Part C predictive models by facility type for Rural Multilane Highway Segments [28].
Table 2. Data needed to calibrate Part C predictive models by facility type for Rural Multilane Highway Segments [28].
Data ElementData NeedDefault Assumptions
RequiredDesirable
Segment lengthX Actual data required
Average annual daily traffic (AADT)X Actual data required
Lane widthX Actual data required
Shoulder widthX Actual data required
Presence of LightingX Assume no lighting
Use of automated speed enforcement XBase default on current practice
Median widthX Actual data required
Table 3. Regression coefficients for four-lane highways in HSM [28].
Table 3. Regression coefficients for four-lane highways in HSM [28].
Facility Type/SeverityAbc
4-Lane Total−9.0251.0491.549
4-Lane KABC−8.8370.9581.687
4-Lane KAB−8.5050.8741.740
Table 4. Key aspects of the studied road segments.
Table 4. Key aspects of the studied road segments.
HighwayIDStart Point (km)Endpoint (km)Studied Length (km)Total Number of Crashes (2009–2019)Traffic VolumeData Available
20162017201820192020
SP 255255_S012.848.145.31837XXXXX
255_S0277.183.16.0215XXXXX
SP 318318_S01235.7236.10.424 XX
SP 330330_S01241.0267.326.31865 XXXX
330_S02267.3304.036.73144XXXXX
330_S03304.0318.514.53175XXXXX
SP 334334_S01319.3349.530.22015XXXXX
334_S02349.5396.046.51856XXXXX
334_S03396.0406.010.01237XXXXX
SP 345345_S0119.431.111.7612XXXXX
345_S0231.139.18.0377 XX
Table 5. AADT of the studied period.
Table 5. AADT of the studied period.
Homogeneous SegmentsAADT (veh/day)
SensorHighway IDDirectionStart Point (km)Endpoint (km)20162017201820192020
SAT10255_S01North48.12.831735725604959696374
SAT10255_S01South2.848.134295889616463946144
SAT11255_S02North83.177.166005741681369016855
SAT11255_S02South77.183.179726602835784618498
TESC2318_S01North235.7236.181958195819581951683
TESC2318_S01South236.1235.784888488848884881638
SAT01330_S01North241267.895659573958195758913
SAT01330_S01South267.824183599322943794008570
SAT04330_S02North267.830410,22213,92813,84213,31812,311
SAT04330_S02South304267.827,88913,96913,85113,65612,643
SAT05330_S03North304318.528,12529,50519,36030,58929,153
SAT05330_S03South318.530427,79329,69019,32031,17729,558
SAT06334_S01North319.3349.510,53510,98511,68211,0909764
SAT06334_S01South349.5319.3925110,74611,79511,1619766
SAT08334_S02North349.539648984281416244023910
SAT08334_S02South396349.542574252446143863887
SAT09334_S03North396406928413,75815,52416,28315,655
SAT09334_S03South406396793714,15914,80115,03113,896
SAT13345_S01East31.119.464026750668966926548
SAT13345_S01West19.43154695666554157885748
SAT13345_S02East3631.164026750668966926548
SAT13345_S02West31.13654695666554157885748
Table 6. Key information regarding the observed crash data throughout the study period.
Table 6. Key information regarding the observed crash data throughout the study period.
Severity TypeTotalFI
Year of Study20162017201820192016201720182019
1653159713981301451467406415
Mean2.322.241.961.830.630.660.570.58
Standard deviation3.473.082.943.081.181.131.131.17
Max33272839910715
Min00000000
Table 7. Proportion of crash data categorized by crash type.
Table 7. Proportion of crash data categorized by crash type.
Collision TypeFIPDOTotal
Single vehicle0.6490.7550.724
Multi-vehicle (total)0.3510.2450.276
Angle0.0370.0160.022
Head-on0.0110.0010.004
Rear-end0.2070.1480.165
Sideswipe0.0730.0560.061
Other multi-vehicle0.0220.0250.024
Total Crashes1.0001.0001.000
Table 8. Estimated Npredicted, Nexpected, and Cx.
Table 8. Estimated Npredicted, Nexpected, and Cx.
Severity2016201720182019Four-Years
Observed CrashesTotal16531597139813015949
FI4514674064151739
Predicted Crashes
(Npredicted)
Total5655705455872267
FI182186181191741
Expected Crashes (EB)
(Nexpected)
Total16221581140212985892
FI4574724204221774
CxTotal2.922.802.572.222.62
FI2.472.502.242.172.35
Table 9. The goodness of Fit of the HSM predictive model by MAD, MAPE, and RMSE tests.
Table 9. The goodness of Fit of the HSM predictive model by MAD, MAPE, and RMSE tests.
GOFCalibrated Predicted CrashesExpected Crashes
MADMAPERMSEMADMAPERMSE
Total4.4453%8.590.8010%1.33
FI1.9278%3.380.8635%1.35
Table 10. Estimated R2 for Npredicted and Nexpected by year and by severity type.
Table 10. Estimated R2 for Npredicted and Nexpected by year and by severity type.
Severity TypeTotalFI
Year of Study20162017201820192016–201920162017201820192016–2019
Calibrated Npredicted0.260.430.340.400.450.120.190.100.230.24
Nexpected0.980.980.990.990.990.850.870.890.910.88
Table 11. Result of all the GOF tests applied for calibrated predicted crashes.
Table 11. Result of all the GOF tests applied for calibrated predicted crashes.
GOFMADMAPERMSER2
Total4.4453%8.590.45
FI1.9278%3.380.24
Table 12. Result of all the GOF tests applied for expected crashes.
Table 12. Result of all the GOF tests applied for expected crashes.
GOFMADMAPERMSER2
Total0.8010%1.330.99
FI0.8635%1.350.88
Table 13. Comparison of estimated variance in fatalities and AADT on state highways in recent years.
Table 13. Comparison of estimated variance in fatalities and AADT on state highways in recent years.
YearFatal Crash CountsMean% Change Compared toAADTMean% Change Compared to
Previous YearMeanPrevious YearMean
201518721836No data1.94%No data----
20161853−1.01%0.90%14,59814,1373.26%
201719113.13%4.06%14,149−3.08%0.08%
20181876−1.83%2.15%14,1830.24%0.32%
201919292.83%5.04%15,1286.66%7.00%
20201651−14.41%−10.10%12,092−20.07%−14.47%
202117636.78%−4.00%16,76438.64%18.58%
Table 14. Key aspects related to crash data from 2016, 2017, 2018, 2019, and 2020.
Table 14. Key aspects related to crash data from 2016, 2017, 2018, 2019, and 2020.
Severity TypeTotalFI
Year of Study2016201720182019202020162017201820192020
16531597139813011182451467406415389
Mean2.322.241.961.831.660.630.660.570.580.55
Standard Deviation3.473.082.943.082.731.181.131.131.171.13
Max332728392791071511
Min0000000000
Table 15. HSM prediction model estimation compared to observed crashes.
Table 15. HSM prediction model estimation compared to observed crashes.
Severity2016201720182019Four-Years Median2020
Observed CrashesTotal165315971398130114871182
FI451467406415435389
Predicted Crashes
(Npredicted)
Total565570545587567498
FI182186181191185164
Calibrated Predicted Crashes
(Cal. Npredicted)
Total-----1308
FI-----386
Expected Crashes (EB)
(Nexpected)
Total162215811402129814761205
FI457472420422443394
Table 16. Result of all the GOF tests applied for calibrated predicted crashes, including 2020.
Table 16. Result of all the GOF tests applied for calibrated predicted crashes, including 2020.
GOFMADMAPERMSER2
Total (4 years)4.4453%8.590.45
FI (4 years)1.9278%3.380.24
Total (2020)1.3380%2.300.32
FI (2020)0.62114%1.040.17
Table 17. Result of all the GOF tests applied for expected crashes, including 2020.
Table 17. Result of all the GOF tests applied for expected crashes, including 2020.
GOFMADMAPERMSER2
Total (4 years)0.8010%1.330.99
FI (4 years)0.8635%1.350.88
Total (2020)0.2415%0.370.98
FI (2020)0.2952%0.440.90
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mendes, O.B.B.; Larocca, A.P.C.; Rodrigues Silva, K.; Pirdavani, A. Assessing the Performance of Highway Safety Manual (HSM) Predictive Models for Brazilian Multilane Highways. Sustainability 2023, 15, 10474. https://doi.org/10.3390/su151310474

AMA Style

Mendes OBB, Larocca APC, Rodrigues Silva K, Pirdavani A. Assessing the Performance of Highway Safety Manual (HSM) Predictive Models for Brazilian Multilane Highways. Sustainability. 2023; 15(13):10474. https://doi.org/10.3390/su151310474

Chicago/Turabian Style

Mendes, Olga Beatriz Barbosa, Ana Paula Camargo Larocca, Karla Rodrigues Silva, and Ali Pirdavani. 2023. "Assessing the Performance of Highway Safety Manual (HSM) Predictive Models for Brazilian Multilane Highways" Sustainability 15, no. 13: 10474. https://doi.org/10.3390/su151310474

APA Style

Mendes, O. B. B., Larocca, A. P. C., Rodrigues Silva, K., & Pirdavani, A. (2023). Assessing the Performance of Highway Safety Manual (HSM) Predictive Models for Brazilian Multilane Highways. Sustainability, 15(13), 10474. https://doi.org/10.3390/su151310474

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop