[go: up one dir, main page]
Next Article in Journal
UWB and MEMS IMU Integrated Positioning Algorithm for a Work-Tool Tracking System
Previous Article in Journal
Mercury Uptake and Transport by Plants in Aquatic Environments: A Meta-Analysis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Applied Sciences
Volume 11
Issue 19
10.3390/app11198828
applsci-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Severity, Spatial Pattern and Statistical Analysis of Road Traffic Crash Hot Spots in Ethiopia

by
Alamirew Mulugeta Tola
1,2,*,
Tamene Adugna Demissie
2,
Fokke Saathoff
1 and
Alemayehu Gebissa
1
1
Faculty of Agricultural and Environmental Sciences, Geotechnics and Coastal Engineering, Rostock University, 18051 Rostock, Germany
2
Faculty of Civil & Environmental Engineering, Jimma University, Jimma 378, Ethiopia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(19), 8828; https://doi.org/10.3390/app11198828
Submission received: 15 August 2021 / Revised: 17 September 2021 / Accepted: 21 September 2021 / Published: 23 September 2021
Browse Figures
Figure 1
<p>Map of the study area in Ethiopia.</p> "> Figure 2
<p>Number of deaths recorded in Ethiopia from 2010–2017.</p> "> Figure 3
<p>(<b>a</b>) Result of spatial autocorrelation; (<b>b</b>) Attribute table result of spatial autocorrelation.</p> "> Figure 4
<p>RTC Hot spot map of Ethiopia regions and federal territories.</p> "> Figure 5
<p>(<b>a</b>) Total lives lost due to RTCs; (<b>b</b>) Percentage deaths (Female versus Male).</p> "> Figure 6
<p>(<b>a</b>) Crashes occurred in a day of the week: from 2014–2017; (<b>b</b>) Crashes occurred in the time of the day: from 2014–2017.</p> "> Figure 7
<p>Crashes analyzed by collision types: 2014–2017.</p> "> Figure 8
<p>The spatial autocorrelation of Oromia in (<b>a</b>) Graph and (<b>b</b>) Attribute table.</p> "> Figure 9
<p>Crash hot spot map of Getis Ord <span class="html-italic">Gi*</span> statistics of Oromia region.</p> ">
Versions Notes

Abstract

:
The reduction of traffic crashes, as well as their socio-economic consequences, has captivated the attention of safety professionals and transportation agencies. The most important activity for an effective road safety practice is to identify hazardous roadway areas based on a spatial pattern analysis of crashes and an evaluation of crash spatial relations with neighboring areas and other relevant factors. For decades, safety researchers have adopted several techniques to analyze historical road traffic crash (RTC) information using the advanced GIS-based hot spot analysis. The objective of this study is to present a GIS technique for identifying crash hot spots based on spatial autocorrelation analysis using a four-year (2014–2017) crash data across Ethiopian regions, as well as zones and towns in the Oromia region. The study considered the corresponding severity values of RTCs for the analysis and ranking of crash hot spot areas. The spatial autocorrelation tool in ArcGIS 10.5 was used to analyze the spatial patterns of RTCs and then the Getis Ord Gi* statistics tool was used to identify high and low crash severity cluster zones. The results showed that the methods used in this analysis, which incorporated Moran’s I spatial autocorrelation of crash incidents, Getis Ord Gi* and crash severity index, proved to be a fruitful strategy for identifying and ranking crash hot spots. The identified crash hot spot areas are along the entrance to and exit from Addis Ababa, Ethiopia’s capital city, so the responsible bodies and traffic management agencies should give top priority attention and conduct a thorough study to reduce the socio-economic effect of RTCs.

1. Introduction

Road traffic safety has become a major apprehension for human beings since the emerging of roadway transport and motor vehicles. According to the World Health Organization (WHO), RTC (road traffic crash) is the 8th most common means of death for all age groups. Nowadays, it is the top reason for the death of children and young adults aged 5–29 years [1]. Road traffic safety analysis has been employed to save the loss of lives by understanding the cause of traffic crashes and coming up with safety mitigations. The analysis aims to investigate pieces of information needed by decision makers to apply suitable safety measures to eliminate and minimize the occurrence of traffic crashes [2].
Because of the presence of distance in the road networks, the spatial pattern of crashes must be examined in traffic safety studies. Spatial analysis is the inspection of crash occurrence patterns by considering their relative locations or zones. Traffic crashes meet the main characteristics of spatial heterogeneity and spatial dependence of point data. Spatial dependence belongs to the influence of events at a location by neighboring events, while spatial heterogeneity happens when the spatial relationships among observed incidents and random parameters in the developed model are not established spatially [3].
The emergence of GIS has provided a vital tool for community health study [4,5,6]. Moellering used a Geographical Information System (GIS) to analyze traffic collisions for the first time in 1976 in his work of “a computer-animated film in the analysis of geographical patterns of traffic crashes” [7]. Since then, the geographical information system has been used extensively in road traffic safety studies over the past five decades [3,8]. Its application varies from simple mapping and visualization roles to further advanced methods such as spatial statistical models and the analysis of large data methods. Nowadays, the precise and exact location of RTCs and their attributes are stored in the GIS database. GIS software allows us to gather spatial data in which we can store, manipulate, analyze and visualize it with ease [9,10]. Even though a comprehensive review and evaluation of analytical approaches have been conducted [11], numerous research works have demonstrated that spatially enriched crash analysis shows potential in establishing a better insight into road safety by revealing areas with safety issues [12,13].
In terms of crash prevention programs, there are minimal indications regarding road safety measures being implemented in low and middle-income countries [14], where the majority of crashes occur. Ethiopia (a low-income country), along with Uganda, Bangladesh and Vietnam, is among the four nations in the world with an exceptionally high health concern, with more than 1000 fatalities per 100,000 motorized vehicles [15]. Therefore, it is important to map crash hot spots by analyzing the severity and spatial pattern of crash incidences in Ethiopia in order to provide information on where to invest a limited budget to improve road safety efficiently.
In Ethiopia, the majority of road safety analysis methodologies employed by transportation authorities and safety specialists are classic descriptive approaches focused on quantifying and summarizing crash data. While crashes are random incidents that occur in space and time [10], they also reflect spatial dependency and spatial autocorrelation, which should be taken into account while analyzing them [16]. Furthermore, at present, in Ethiopia, GIS is not used widely in crash data recording and identifying RTC hot spot locations. To identify RTC hot spots, a GIS-based spatial analysis has been a promising tool and widely used [17,18]. The GIS-based hot spot analysis result is presented accurately on a single map with the related attribute data of each RTC which may advance the understanding of road safety researchers regarding the reasons for each crash occurrence [19].
With the most up-to-date analytical techniques of the GIS’s Getis Ord Gi* statistics, this study aimed to analyze the severity and spatial pattern of crash incidents, as well as map crash hot spots in Ethiopia. While some scholars have studied crash events in Ethiopia [20,21,22,23], none have utilized spatially enriched and statistically integrated analyses like Getis Ord Gi* statistics. Hence, the current study is significant in that it illustrates the benefit of utilizing spatial autocorrelation and statistical techniques to identify crash-intensive prone zones, as well as demonstrating its effectiveness in Ethiopia. Moreover, the study will provide guidance to decision-makers on where to best invest or implement safety measures.

2. Literature Review

There is no universally acknowledged definition of Road Traffic Crash (RTC) hot spot identification [10]. Some safety researchers prioritize sites by crash rate (i.e., crashes per vehicle-kilometers for segments and crashes per entering vehicles for intersections), others use crash frequency (i.e., crashes per year), some use spatial analysis of crashes and others adopt the integration of the above and the rest uses regression modeling (i.e., the predictive approach of Crash Prediction Model, CPM) [24]. RTC hot spot analysis, in general, attempted to identify and prioritize road segments in need of immediate safety improvements to achieve valuable crash reduction through effective safety mitigation. In addition to crash analysis, local expertise and expert judgment are used to identify RTC hot spots. This work (road safety management) comprises the following tasks [25,26]:
  • Targeting or identifying crash hot spots on the road network;
  • Studying safety issues in each hot spot;
  • Identifying contributing factors and design mitigations;
  • Evaluating the safety effects of the possible mitigations;
  • Prioritizing the hot spots to apply cost-effective safety mitigations; and
  • Evaluating the effectiveness of applied treatments.
This review has mainly focused on the first task of identifying crash hot spots on the road network. Since RTC analysis is such an important activity for improving roadway safety, it has been extensively researched in the academic press. Despite the lack of consensus on crash hot spot analysis, researchers and experts have developed various methods for crash event analysis. The simplest and most straightforward method for crash hot spot analysis is detecting where the crash rate per unit exposure exceeds a given normal threshold [27].
Although most of the traditional approaches of crash analysis focus on the time dimension, nowadays, the spatial dimension of traffic crashes has got more attention from researchers [16,28]. This facilitates the application of GIS into crash hot spot analysis. In recent years, several RTC databases recorded the precise locations of crashes through the application of GPS devices; hence, it is no longer necessary to identify the crash hot spot section of a road from aggregated data [2,29]. By having these accurate locations of RTC, road safety analysts can focus on the highly clustered crash locations.
Generally, the most frequently used methods of crash hot spot analysis can be categorized into two. These are:
  • Non-Spatial Model Analysis (NSMA): which uses the traditional approaches of statistics such as regression models [30], Empirical Bayes [31], and full Bayesian [32];
  • Geo-Statistical Analysis (GSA): by analyzing the spatial units of crashes (i.e., Density Estimation, DE) [33,34] or spatial arrangement of each crash attribute value (i.e., Spatial-Autocorrelation, SA) [35,36,37].
As compared to NSMA, Geo-Statistical crash hot spot analysis needs fewer data and is easier to apply due to the simplified mathematical calculations [38] and its integration of crash incidence with spatial factors and, also, the presentation of clear visualization in the result. The adoption of the spatial arrangement of attribute values from each unit is more advantageous than using a spatial unit of crashes due to its consideration of spatial dependence of attribute values and the geographical location of incident points [39]. For a thorough understanding of the Geo-Statistical analysis of crashes, readers are recommended to a review paper [40] dedicated to reviewing the scientific literature on spatial approaches and spatial analysis in road safety.
GSA can be further divided into two groups such as Global Indexes (GI) for instance Global Moran’s I (Spatial-Autocorrelation), Getis-Ord G statistic and Geary’s C; and Local Indexes (LI) these are Local Anselin Moran’s I (Cluster and Outlier Analysis), planar Kernel Density Estimation (KDE), Getis Ord Gi* and kriging. Except for kriging and KDE, the rest incorporates a procedure for testing the statistical significance of clustered incidents. In the global index analysis, the spatial pattern incidence of the entire network is assessed and it analyzes whether data is clustered in a group, dispersed, or randomly distributed, whereas the local indexes are used to study the microscopic pattern of incidents to determine the spatial location and extension of these clusters [41,42].
Several studies have used planar KDE [17,33], which works by creating a continuous surface of total event density within a search bandwidth and network space KDE [34,43], which is generalized to calculate the crash density over a distance unit rather than an area unit for crash hot spot identification. The major limitations of both planar and network space KDE’s are the lack of statistical significance test in the analysis [17,34,44] and also, there are no criteria for prioritizing crash hot spots [34]. Research has developed and used KDE+ as an effective tool for crash hot spot identification and site prioritization to overcome the limitation of the network and planar KDE approach [45,46]. The KDE+ analysis works based on the standard KDE approach with the addition of statistical significance testing of clustered points.
Compared to KDE, for both Moran’s I and Geary’s C the statistical significance of clustered crashes is determined using Z-Score [47,48]. Indeed, Moran’s I and Geary’s C follow the global statistics approach which is the measure of the whole study network. The local statistics approach (such as Getis-Ord Gi* statistics [37] and local Moran’s I [35]) is superior to the global statistics system for studying spatial variance and spatial dependency. The Getis-Ord Gi* statistic, in particular, has been validated to identify statistically significant values of crash hot spots or cold spots and has proven to be useful [10,49].

3. Methods and Materials

3.1. Study Area and Data Collection

This research was carried out in Ethiopia to identify and prioritize Road Traffic Crash (RTC) hotspot regions, with a particular emphasis on the Oromia region (see Figure 1). RTC data of Ethiopia was obtained from the Ethiopian federal police commission, traffic police division. The commission records crash data for all regions and federal territories. Crash data of four consecutive years (from 2014 to 2017) were used for crash severity analysis and hot spot identification. All collision types acquired from the police commission (i.e., single-vehicle collision, multiple vehicle collision, collision with pedestrians, collision with animals and so on) were utilized for analysis. The Ethiopian Geospatial Information Agency (EGIA) provided a map of Ethiopia divided into regions and federal territories, as well as a map of the Oromia region, divided into zones and towns. In these maps, roadway networks and other important features were included as a shapefile.

3.2. Methods

For RTC hot spot identification, this study used the spatial statistics toolbox of ArcGIS 10.5. Supporting tools such as Google Earth and Google Maps were also utilized. For RTC hot spot identification, the procedure used in this study can be divided into the following steps: (i) projecting the map, (ii) computing the severity value of each crash, (iii) analyzing the spatial pattern of RTCs to measure the fluctuating threshold distance of RTC clustering by using Moran’s I index of incremental spatial autocorrelation tool and (iv) identifying and prioritizing highly clustered RTCs or hot spot zones. The detailed steps and the theoretical basis are discussed as follows.
Instead of using a geographical coordinate system, the map used to identify crash hot spots was a projected coordinate system. The map dataset obtained from Ethiopian Geospatial Information Agency (EGIA) was in a geographic coordinate system. Since distance is among parameters to be considered in crash hot spot analysis [50], the analyzed maps were projected using the UTM (Universal Transverse Mercator) coordinate system in World Geodetic System WGS84.

3.2.1. Crash Severity

In addition to crash rates, crash hot spot zones are classified based on their severity index. The crash severity index determines the weight of a single crash. So, to accurately identify high or low clustered zones, it is important to incorporate crash severity index in hot spot analysis. The concept of crash severity index is accrediting the higher values to severer crashes based on the expenditures. Many past studies associated crash severity weights in crash hot spot analysis, for instance, Geurts et al. [51] used severity index developed by the Belgian government by giving values of 5, 3 and 1 for fatal, serious injury and slight injury, respectively. Because Ethiopia lacks a traffic crash costing platform for particular crash severity levels, the crash severity index developed by the Roads and Traffic Authority of New South Wales [52] was utilized in this study. In this system, each crash incident is provided a value of 3.0 for fatal, 1.8 for serious injury, 1.3 for slight injury and 1.0 for property damage only crashes. The crash severity index of each zone can be computed by Equation (1) [52].
S I = 3.0 ⋅ X 1 + 1.8 ⋅ X 2 + 1.3 ⋅ X 3 + 1.0 ⋅ X 4
where X1 is fatal crashes, X2 is serious injury crashes, X3 is slight injury crashes, X4 is Property-damage-only crashes.

3.2.2. Getis Ord Gi*

For RTC hot spot identification, spatial analysis using local spatial autocorrelation is preferable. Local Moran’s I is among the well-known local spatial autocorrelation approach used commonly in motor vehicle crash hot spot analysis [53]. However, the family of Moran indices does not differentiate between hot or cold spots. Getis Ord Gi* is, therefore, more appropriate since it distinguishes clusters with high and low feature attribute values among local events. In this study, Getis Ord Gi* statistics method was used to identify RTC hot spots.
Gi* statistics aims to investigate the presence of a spatial pattern for an arbitrary variable X, where a selected event (with a value of xi) is autonomously connected to the field. As a result, if xi is analogous to adjacent regions, it can be seen that variable X has spatial autocorrelation over area i. A simple form of Gi* statistics is defined as Equation (2), which is derived by dividing the study zone into n indefinite extent of regions, each with accurate Cartesian coordinates and a central point i (i = 1, 2, 3, … n) [37].
G i * = ∑ j = n n w ij x i ∑ j = 1 n x j
where Gi* is a statistic that describes the spatial dependency of feature i, xj is the value of variable X at feature location j. wij is the spatial weight between feature i and j.
The conceptualized-spatial relationship around distance d (for instance Cartesian distance) is used to calculate wij. The result of Gi* statistics may vary based on the choice of d. The value of d is a user-defined threshold. The most straightforward way to think about wij is in binary form, with 1 indicating inclusion and 0 indicating exclusion of the association between i and j events. However, in practice, wij can have non-binary values and the total weights (Wi) are expressed in Equation (3).
W i = ∑ j = 1 n w i j
Readers are referred to a paper on Getis Ord Gi* statistics [54] to learn more about the statistical elements such as Gi* expectation, Gi* variance and sample mean and variance for the variable X. Finally, as shown in Equation (4), Gi* statistics are typically standardized using the sample mean and variance for a normal asymptotical condition.
Z G i * = ∑ j = 1 n w i j x j − x ¯ ∑ j = 1 n w i j 2 s n ∑ j = 1 n w i j 2 − ∑ j = 1 n w i j 2 n − 1
This standardized Gi* statistics is a Z-score, which is associated with the statistical significance of each target region. The null hypothesis for the Getis Ord Gi* hot spot analysis tool is that crashes are caused by Complete Spatial Randomness (CSR) or are randomly distributed. The null hypothesis would be rejected or accepted based on the Z-score and p-values obtained from the Getis Ord Gi* analysis tool. The p-values for the spatial pattern analysis is the probability that the observed crash incidents occurred in some arbitrary manner. When the returned p-values from the hot spot analysis tool is very small, it is very unlikely that the observed crash incident is to be distributed randomly; thus, the null hypothesis can be rejected. The tails of the standard normal distribution curve have very small p-values and very high absolute values of Z-scores. When the value of the Z-score is closer to zero it implies that there is a random distribution of spatial events in the region. The maximum absolute values of Gi* statistics, on the other hand, represent clusters of low-valued events for negative or high-valued events for positive.
ArcGIS’ Getis Ord Gi* hot spot analysis tool uses the Gi* statistics index to identify a significant hot/cold spot based on neighboring attribute values. Regions with statistically significant positive Z-scores (hot spots) are surrounded by neighbors with high feature values, whereas regions with statistically significant negative Z-scores (cold spots) are surrounded by neighbors with low feature values. If the local aggregate of the target area and its neighboring values differ significantly from the likely local value to be random distribution, a statistically significant Z-score is produced.
The final statistical significance for hot spot identification has been adjusted by taking multiple testing and spatial dependency into account and has taken the form of Equation (5). This final equation is incorporated in ArcGIS’ Getis Ord Gi* hot spot analysis tool to identify and prioritize crash hot spots. Thus, the final form of statistical significance for hot spot identification (in Equation (5)) was used in this study by running ArcGIS’ Getis Ord Gi* hot spot analysis tool.
G i * = ∑ j = 1 n w i j x j − ∑ j = 1 n x j n ∗ ∑ j = 1 n w i j ∑ j = 1 n x j 2 − ∑ j = 1 n x j n ∗ n ∑ j = 1 n w i j 2 − ∑ j = 1 n w i j 2 n − 1

3.2.3. Spatial Autocorrelation

In the spatial planar analysis, the distance methods used are Manhattan or Euclidean. In the analysis of spatial autocorrelation, the choice of conceptualization of the spatial association among crash events could be based on knowledge of the interactions among the features to be analyzed. There are several ways to conceptualize the spatial association between events, including fixed distance, inverse distance, inverse squared distance, K-nearest neighbors, the zone of indifference, the space-time window method and contiguity edges and corners. The result of the crash hot spot from global spatial autocorrelation and Getis Ord Gi* can be improved by applying Anselin Local Moran’s I to the aggregated traffic crashes [55]. The selection of appropriate spatial relationships among features for the use of local spatial autocorrelation analysis may help to reflect the distributional and spatial situation of definite target features [56]. The fixed distance approach was used in this study to conceptualize the spatial association among events.
In this study, the incremental spatial autocorrelation tool of ArcGIS 10.5 was used to estimate the fluctuating value of spatial autocorrelation as a distance threshold. In this tool, Global Moran’s I index approach is employed to quantify the distance bandwidth (i.e., a distance in which maximum crash events are clustered [57]) across the entire area. The equations of Moran’s I index (I), the expected value (E[I]), variance (V[I]) and Z-score (Z_Score) are presented in Equation (6) to Equation (9) respectively:
I = n ∑ i = 1 n ∑ i = 1 n w i , j ∗ ∑ i = 1 n ∑ i = 1 n w i , j x i − X ¯ x j − X ¯ ∑ i = 1 n x i − X ¯ 2
E I = − 1 n − 1
V I = E I 2 − E I 2
Z _ S c o r e = I − E I V I
where xi is an attribute value of target feature at location i, xj is an attribute value of neighboring feature at location j, wij is the spatial weight between features at location i and j.
A statistically significant Z-score (peaked) shows a critical distance threshold at which spatial autocorrelations are highly clustered. The first peak was used as the distance threshold in this study because it tends to best describe the spatial variation analysis and is recommended for use when there are multiple peaks [58].

4. Results and Discussions

This section presents statistical trends of crash fatalities as well as the findings of crash hot spot analyses for Ethiopian regions and federal territories, as well as Oromia zones and special towns, followed by discussions.

4.1. Death Rate and Trend of Road Safety in Ethiopia

The analysis to figure out the trends of fatalities that occurred due to road traffic crashes (RTCs) in Ethiopia was done based on eight consecutive years (2010–2017). The crash history from Ethiopian Traffic Police Commission has shown that in 2010, 2541 fatalities occurred as a result of traffic collisions (see Figure 2). There was a slight increase in fatalities until 2013, then a significant increase in 2014, when the number of people killed in crashes shockingly reached 4883. Even if a slight decrease was seen in 2015, it increased again and finally, 5118 deaths were recorded in 2017. According to the data from the Ethiopian Roads Authority (ERA) [59], from 2009 to 2017, 77,980 kilometers of road networks were newly constructed total highways (asphalt and gravel), with an annual growth rate of 13.08%. Despite the Ethiopian government’s efforts to improve access and mobility for road users, the significant growth of motor vehicles was a contributing factor to the increase in traffic fatalities. The motor vehicle registration record of the Ethiopian Federal Transport Authority disclosed that in 2009 the number of vehicles in Ethiopia was about 276,794 in which had exceeded 831,000 motor vehicles in 2017 with an annual average growth rate of 15.28% [60]. The Ethiopian Police Commission-Department of Traffic Police reported that the day-to-day increase of RTC fatalities caused by alcohol use has become their main concern. As previously stated, the reasons for the increase of fatalities may be attributed to a significant growth rate of motor vehicles relative to road networks (the road network’s growth rate is 13.08%, while motor vehicles’ growth rate is 15.28%) and alcohol consumptions of drivers. The detail of these trends is presented in Figure 2.

4.2. Crash Hot Spot Analysis of Ethiopian Regions and Towns

Crash hot spot analysis was done at a regional level as a polygon feature dataset by applying the computed crash severity values of each region in the attribute table. The incremental spatial autocorrelation toolbox of ArcGIS was used to choose the distance bandwidth and the results are shown in Figure 3. The maximum and first peak distance obtained from incremental spatial autocorrelation analysis was 352.441 kilometers, indicating that RTC incidents are highly clustered at a distance band of 352.441 kilometers with a statistical significance of 0.01 (the p-values obtained is 0.402% < 1%; thus, the confidence level is 99%). The bigger threshold distance values can be explained by the fact that the polygons in the analysis (i.e., the regions and administrative cities) have larger area coverages. For example, Oromia, Ethiopia’s largest region, has an area of around 355,423 km2. In distance computations for polygon features, feature centroids are considered. The weighted mean center of all feature components is used to compute the centroid for various portions of the polygon feature. The weighted mean center for point features is 1, length for line features and area for polygon features [61]. Thus, in the case of this analysis, features such as polygons with wider regions were expected to have a larger threshold distance. Furthermore, the default Beginning Distance was applied in this study while executing the incremental spatial autocorrelation. The minimum distance for which each feature in the dataset has at least one neighbor is set as the default value for Beginning Distance. As a result, the beginning distance for incremental spatial autocorrelation analysis increases since larger regions (polygons) demands longer distances to reach at least one neighbor from their centroids. In fact, the higher beginning distance used to determine the first peak has a direct influence on the value of the threshold distance. This value was used as the threshold distance in Getis Ord Gi* crash hot spot identification. Accordingly, crash hot spot identification and ranking of Ethiopian regions and federal territories were performed by considering the spatial patterns and spatial dependence of crashes. The result obtained from Getis Ord Gi* crash hot spot analysis is presented in Figure 4 and Table 1.
With confidence intervals of 99% and 95%, Addis Ababa and the Oromia region were identified as crash hotspots. Furthermore, with a 90% confidence interval, two federal territories, ‘Dire-Dawa’ and ‘Harari,’ were revealed as low clusters of crash severity areas or cold spots. Crash hot spot zones are primarily characterized as highly densified areas that are associated with the presence of numerous activities that result in a large number of pedestrians, traffic volumes, populations and intersections. Addis Ababa, the capital city of Ethiopia was revealed as the first crash hot spot. The city is the home of AU (African Union), the Economic Commission for Africa (ECA) and over one hundred (100) embassies of several countries. The financial and administrative capital city of Ethiopia, Addis Ababa, is facing persistent growth and change. The city is transforming from its current administrative center to the central business district and an industrial center. Due to this change and rapid growth, there is an excessive transportation demand for the mobility of goods and people and, also, there is a higher attraction of people to the city in which a population of 3.6 million in 2013 is forecasted to be nearly 10 million in 2037 [62]. Whereas the city’s efficient share of roadway infrastructures in land use requires 20% to 25%, the current share is 7%, indicating that the city’s road transport system is inadequate [63].
The city’s modernization, combined with an increase in motor vehicle ownership, increases the number of vehicles in the city; for example, of all registered motor vehicles in Ethiopia, 77% are found in Addis Ababa, with a yearly growth rate of 5.8% [64]. Walking (pedestrians) is a substantial portion of Addis Ababa’s transportation modes, accounting for even more than 60% of daily journeys [65]. This high proportion of pedestrian road users in the city, combined with limited sidewalk facilities (such as crowded sidewalks, unpaved or inaccessible sidewalks), increases the number of pedestrian-related crashes. This is demonstrated by the fact that about 89% of the city’s recorded crash causalities were pedestrian-related [66].
Addis Ababa City, the first-ranked crash hot spot zone, has a radial-shaped road network with five major roadways radiating in and out of the city’s central business and administrative district. All traffic flows entering and leaving the city via these five roadways pass through the Oromia region, which was prioritized as the second crash hot spot area. High traffic volumes entering and leaving Addis Ababa have an impact on the Oromia region, which shares a common periphery on all sides with Addis Ababa (Finfinnee). Oromia was prioritized as the second hot spot due to higher traffic volumes, spatial dependence with Addis Ababa, larger size (inland coverage, roadway length and population number) and unsustainable traffic management in the region.

4.3. Crash Analysis of Oromia Region

4.3.1. Number of Deaths Due to RTCs in the Oromia Region

According to the Oromia traffic police bureau’s eight-year (2010–2017) crash database record, fatalities that occurred due to RTCs had shown an increasing trend. Despite population growth, the number of crash fatalities increased by more than double from 906 in 2010 to 1882 in 2017 (see Figure 5a). The Ethiopian Federal Transport Authority has stated that there has been significant motor vehicle growth in the Oromia region, which has a strong correlation with the increased fatalities [60]. According to the Oromia Traffic Police Bureau crash data, denying priority for pedestrians pass, unsafe vehicle passing, driver’s alcohol and drug consumptions, over-speeding, driver workload, insufficient vehicle following time and vehicle problems are the recognized significant factors for the increasing number of deaths in Oromia [67]. This increase in fatalities in the region demonstrates the impact of road traffic crashes on the socio-economic well-being of the country as a whole and the Oromia region in particular. To minimize road traffic crashes and their consequences, the Oromia transport bureau, traffic management, traffic police, road authority and concerned governmental and non-governmental organizations should all work together to address this serious problem. Figure 5a depicts the number of RTC fatalities in the Oromia region over eight consecutive years.
Male fatalities were 2.7 times higher than females (Figure 5). According to the population census commission’s 2007 survey, the percentage of males and females in Oromia was roughly equal; moreover, males accounted for 73% of the region’s road traffic crash victims. This can be explained by the fact that, in comparison to females, males have greater exposure to road transportation as drivers and passengers for a variety of reasons. For example, males are over-represented in professional driving careers (i.e., jobs in remote areas, operators, taxi, nighttime driving and so on) and females are under-represented in self/paid employment or outdoor activities in Ethiopia [68], which could reduce female vulnerability to RTCs. Furthermore, males are more likely than females to engage in risky driving behaviors such as careless overtaking, the use of alcohol, impatience, exceeding posted speed limits and having a lower risk perception [69,70]. As compared to females, males’ higher use of two-wheeled motor vehicles and bicycles for transport can also affect their susceptibility to traffic crashes [71].

4.3.2. Crash Severity

Table 2 shows that there were 17,259 traffic crashes in the Oromia region during the last four years, from 2014 to 2017.

4.3.3. Crashes in the Day of Week and Time of Day

RTCs occur more frequently on Thursday and then on Sunday in a week and from 12:00 to 13:00, 18:00 to 19:00 and 09:00 to 10:00 in a day. We recommend that the Oromia region’s traffic management, transportation authority and traffic police bureau intervene and monitor traffic flows during peak crash days and times when a higher number of crashes are likely to occur. The details are demonstrated in Figure 6. Thursday is more likely associated with a market day in the Oromia region, when a large number of people may travel from place to place and Sunday is an off-day, but most of the society and road users may prefer to consume alcohol, resulting in more traffic crashes. Crash rates were higher from 12:00 to 13:00 and 09:00 to 10:00 during the day and these are the times when higher traffic flows are observed (peak hours), so the risk of RTC could be increased. Another critical time revealed was 18:00 to 19:00. This was most often associated with overnight collisions when there would be insufficient visibility for vision.

4.3.4. Crashes by Collision Type

In Oromia, crashes were more likely to occur with a collision type of ‘Ran-Off-Road (ROR)’ see Figure 7. A run-off-road (ROR) collision happens when a moving car overturns on impassable ground or goes out of the roadway and strikes an object. From the critical reasons for ROR collisions, over 95% were driver-related factors [72]. Among the significant driver-related factors that contribute to ROR collisions were the internal problem (i.e., heart attack or other health deficiency), too fast for over-compensation, speeding on curves and sleeping or ‘fatigued driving’, respectively, in a descending manner [72]. However, a thorough investigation into the factors that contribute to the frequency of higher ROR collisions in the Oromia region is needed.

4.3.5. Crash Hot Spot Analysis

The threshold distance for identifying crash hot spots was determined using the results of Moran’s I spatial autocorrelation analysis, which is shown in Figure 8. The threshold distance used in Getis Ord Gi* crash hot spot analysis is the distance at which crashes are highly clustered. As a result, the incremental spatial autocorrelation analysis of the Oromia region yielded a threshold distance of 161.764 kilometers with a statistical significance level of 0.10 (such that the p-values obtained is 5.14% < 10% thus, the confidence level is 90%). Even though the current analysis’s first peak has been cut by more than half due to the reduced analysis coverage from regions to zones within a single region (Oromia), the threshold distance remains higher. As explained earlier (in Section 4.2), larger zones (e.g., Borena with 45,463 km2 area coverage and Bale with 44,912 km2 area coverage) necessitate longer distances to reach at least one neighbor from their centroids, resulting in a higher default Beginning Distance and threshold distance. Figure 9 and Table 3 show the results of Getis Ord Gi* statistics for identifying crash hot spots in the Oromia region. Getis Ord Gi* statistics’ crash hot spot analysis identified six hot spot areas, which are listed in the attribute table’s result (see Table 3). The rankings were based on each zona’s computed Z-score.
According to the findings, the East ‘Shewa’ (E. Shewa) zone and ‘Burayu’ town (Burayu T.) have the highest values of crash severity clusters and are ranked as priority crash hot spot areas across the Oromia region with a confidence interval of 99%. According to Section 3.1 of this document, Addis Ababa is encircled by five major roadways. These are locally named as; (1) ‘Ambo-Ber’, (2) ‘Tarma-ber’ (3), ‘Jimma-Ber’, (4) ‘Gojjam-Ber’ and, (5) ‘Kality-Ber’ or ‘Bishoftu-Ber’. The AADT for freight vehicles entering and leaving Addis Ababa is estimated to be 10,725 and 12,890, respectively [73]. More than 70% of these entering and exiting freight vehicle shares have gone through Gate-5 (Kality/Bishoftu-Ber) [73], which is located in the East Shewa zone, Oromia’s first-ranked crash hot spot. ‘Bishoftu-Ber’ is also the key entry point for import and export freight vehicles between Addis Ababa and the port of Djibouti. As a result, the East ‘Shewa’ zone was expected to experience higher traffic volumes and freight vehicles, which are directly related to the increased crash frequencies. When compared to others, issues related to the second crash hot spot zone, ‘Burayu town,’ were primarily associated with poor roadway facilities (such as narrow lane width, deteriorated pavement, lack of pedestrian walkway and so on).

5. Limitation of the Study

The availability of crash data recorded with precise incident locations in coordinates (using Global Coordinate System, GPS) or mileposts that could be linked to roadway geometric characteristics drives a well-organized RTC analysis as well as results for specific roadway segments and intersections. The Ethiopian ministry of transport, regional transport authorities and federal police commission-traffic police departments are in charge of collecting and reporting the essential crash data to safety researchers in need of it. However, crash locations that occur along a road alignment (milepost) are not recorded in Ethiopia. This study was unable to investigate the spatial dependence of crashes with roadway parameters (such as road length, degree of curvature, AADT, gradient, stopping sight distance and others) as a contributing factor to crash occurrences due to a lack of crash locations in the given data. Moreover, the crash severity index of each severity level is derived from the respective crash costs in which these costs vary from region to region. However, due to the lack of an agreed-upon or defined crash costing platform for specific crash severity levels in Ethiopia, the severity index developed by the Roads and Traffic Authority of New South Wales [52] was used in this study.

6. Conclusions

The primary goal of crash hot spot analysis is to answer a critical question in road safety practice: where are the hazardous street areas located? This inquiry is addressed scientifically by studying the spatial pattern of crashes. The advantage of using advanced GIS-based hot spot analysis is not limited to the simple presentation of hot spot roadway areas; it also provides the ability to investigate the spatial dependency of crash incidents and spatial connections with other factors. For the past five decades, a GIS application in the field of traffic safety has aided in the advanced realization of crash characteristics, which is then used as a piece of information to improve traffic safety.
This study aimed to present a GIS application to identify and quantify statistically significant spatial patterns based on crash numbers and severity. Based on the severity of the crashes, crash hot spots were identified using the integration of spatial autocorrelation of crashes and Getis Ord Gi*. The use of spatial autocorrelation has the advantage of allowing statistical analysis of crash spatial patterns. The ability of Gi* statistics to differentiate high crash clusters from low crash clusters makes it a better approach for identifying crash hot spots than Moran’s I index. In this study, rather than the total number of collisions, crash counts and severity values were used.
Following the analysis of spatial autocorrelation of crashes, Gi* statistics were used to identify high and low crash severity clusters. The first analysis result of crash hot spot identification for all Ethiopian regions and federal territories identified Addis Ababa and Oromia region as high crash severity clusters with confidence intervals of 99 percent and 95 percent, respectively. Second, six statistically significant crash hot spots in Oromia were detected through the integration of spatial autocorrelation and Gi* statistics. With a confidence level of 99%, the quantitative results of Z-Scores prioritized East ‘Shewa’ zone and ‘Burayu’ town as the most crash hot spots. All of the six identified crash hot spot zones and towns are near the entrance to Addis Ababa, Ethiopia’s capital city. Addis Ababa is linked to other towns, zones and regions by five major roads and all of these connecting roads pass through the identified crash hot spots in Oromia. It can be concluded that the identified crash hot spot locations are along the entrance and exit of Addis Ababa city; therefore, the concerned bodies and traffic management agencies should give these areas top priority and conduct a thorough study in order to reduce the socio-economic impact of traffic collisions.
The results of the study acknowledged that the applied method of crash hot spot identification (such as the integration of spatial autocorrelation of crashes and Getis Ord Gi*) has the ability to analyze spatial patterns of crashes and identify highly clustered crash severity with a statistical significance. The study demonstrated that using GIS in Ethiopia has a significant benefit when it comes to prioritizing a promising site for safety improvements. Thus, the application of GIS in crash hot spot analysis must have to be considered and need to be used as a tool for road safety study in Ethiopia in the future. The recommendations for future road safety improvements include the development of a national crash hot spot (blackspot) identification manual, the frequent spatial and temporal mapping of crash patterns and the provision of training for safety engineers on the methodologies for identifying crash hot spots, as well as on the advancements of GIS technology to be implemented in the field of road safety. Furthermore, the Ethiopian Ministry of Transport (MoT) must support the crash recording and data-sharing platform with advanced technologies such as developing a computerized crash database system and centralized mobile applications that enable GPS, which improves the identification and analysis of the crash spatial pattern to apply effective safety measures.

Author Contributions

Conceptualization, A.M.T. and A.G.; methodology, A.M.T. and A.G.; software, A.M.T.; data curation, A.M.T.; writing—original draft preparation, A.M.T.; writing—review and editing, A.G., T.A.D. and F.S.; supervision, A.G., T.A.D. and F.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was part of the DAAD-EECBP Homegrown PhD Program-2018. The APC was funded by the Open Access Publication Fund of the University of Rostock and German Research Foundation (DFG).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this study are available from the corresponding author upon reasonable request.

Acknowledgments

We would like to acknowledge the Oromia Traffic Police Bureau and Ethiopian Federal Police Commission-Traffic Police Department for providing crash data. It is our pleasure to thank the German Academic Exchange Service (DAAD) for issuing a Ph.D. scholarship to the first author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WHO (World Health Organization). Global Status Report on Road Safety; World Health Organization: Geneva, Switzerland, 2018. [Google Scholar]
  2. Li, L.; Zhu, L.; Sui, D.Z. A GIS-based Bayesian approach for analyzing spatial-temporal patterns of intra-city motor vehicle crashes. J. Transp. Geogr. 2007, 15, 274–285. [Google Scholar] [CrossRef]
  3. Dereli, M.A.; Erdogan, S. A new model for determining the traffic accident black spots using GIS-aided spatial statistical methods. Transp. Res. Part A Policy Pract. 2017, 103, 106–117. [Google Scholar] [CrossRef]
  4. Kwan, M.P. How GIS can help address the uncertain geographic context problem in social science research. Ann. GIS 2012, 18, 245–255. [Google Scholar] [CrossRef]
  5. Neutens, T. Accessibility, equity and health care: Review and research directions for transport geographers. J. Transp. Geogr. 2015, 43, 14–27. [Google Scholar] [CrossRef]
  6. Goodchild, M.F. Space, place and health. Ann. GIS 2015, 21, 97–100. [Google Scholar] [CrossRef] [Green Version]
  7. Moellering, H. The potential uses of a computer animated film in the analysis of geographical patterns of traffic crashes. Accid. Anal. Prev. 1976, 8, 215–227. [Google Scholar] [CrossRef]
  8. Mohaymany, A.S.; Shahri, M.; Mirbagheri, B. GIS-based method for detecting high-crash-risk road segments using network kernel density estimation. Geo-Spatial Inf. Sci. 2013, 16, 113–119. [Google Scholar] [CrossRef]
  9. Lloyd, C.D. Spatial Data Analysis: An Introduction for GIS Users; Oxford University Press: New York, NY, USA, 2010; ISBN 9780199554324. [Google Scholar]
  10. Loo, B.P.Y.; Anderson, T.K. Spatial Analysis Methods of Road Traffic Collisions; CRC Press: New York, NY, USA, 2015; ISBN 9781315266398. [Google Scholar]
  11. Lord, D.; Mannering, F. The statistical analysis of crash-frequency data: A review and assessment of methodological alternatives. Transp. Res. Part A Policy Pract. 2010, 44, 291–305. [Google Scholar] [CrossRef] [Green Version]
  12. Mehta, G.S. Analyzing Crash Frequency and Severity Data Using Novel Techniques. Ph.D. Thesis, University of Alabama, Tuscaloosa, AL, USA, 2014. [Google Scholar]
  13. Mannering, F.L.; Bhat, C.R. Analytic methods in accident research: Methodological frontier and future directions. Anal. Methods Accid. Res. 2014, 1, 1–22. [Google Scholar] [CrossRef]
  14. Lagarde, E. Road Traffic Injury Is an Escalating Burden in Africa and Deserves Proportionate Research Efforts. PLoS Med. 2007, 4, 170. [Google Scholar] [CrossRef] [Green Version]
  15. Elvik, R.; Høye, A.; Vaa, T.; Sørensen, M. The Handbook of Road Safety Measures, 2nd ed.; Emerald Group Publishing Limited: Bingley, WA, UK, 2009; ISBN 9781848552500. [Google Scholar]
  16. Yao, S.; Loo, B.P.Y.; Yang, B.Z. Traffic collisions in space: Four decades of advancement in applied GIS. Ann. GIS 2016, 22, 1–14. [Google Scholar] [CrossRef]
  17. Anderson, T.K. Kernel density estimation and K-means clustering to profile road accident hotspots. Accid. Anal. Prev. 2009, 41, 359–364. [Google Scholar] [CrossRef] [PubMed]
  18. Vemulapalli, S.S.; Ulak, M.B.; Ozguven, E.E.; Sando, T.; Horner, M.W.; Abdelrazig, Y.; Moses, R. GIS-based Spatial and Temporal Analysis of Aging-Involved Accidents: A Case Study of Three Counties in Florida. Appl. Spat. Anal. Policy 2017, 10, 537–563. [Google Scholar] [CrossRef]
  19. Le, K.G.; Liu, P.; Lin, L. Determining the road traffic accident hotspots using GIS-based temporal-spatial statistical analytic techniques in Hanoi, Vietnam. Geo-Spat. Inf. Sci. 2019, 23, 1–12. [Google Scholar] [CrossRef] [Green Version]
  20. Abegaz, T.; Gebremedhin, S. Magnitude of road traffic accident related injuries and fatalities in Ethiopia. PLoS ONE 2019, 14, e0202240. [Google Scholar] [CrossRef] [Green Version]
  21. Guyu, F. Identifying major urban road traffic accident black-spots (RTABSs): A sub-city based analysis of evidences from the city of Addis Ababa, Ethiopia. J. Sustain. Dev. Afr. 2013, 15, 110–130. [Google Scholar]
  22. Anteneh, A.; Endris, B.S. Injury related adult deaths in Addis Ababa, Ethiopia: Analysis of data from verbal autopsy. BMC Public Health 2020, 20, 926. [Google Scholar] [CrossRef] [PubMed]
  23. Hayidso, T.H.; Gemeda, D.O.; Abraham, A.M. Identifying Road Traffic Accidents Hotspots Areas Using GIS in Ethiopia: A Case Study of Hosanna Town. Transp. Telecommun. J. 2019, 20, 123–132. [Google Scholar] [CrossRef] [Green Version]
  24. Geurts, K.; Wets, G. Black Spot Analysis Methods: Literature Review; Steunpunt Verkeersveiligheid: Diepenbeek, Belgium, 2003; Volume RA-2003-07. [Google Scholar]
  25. Vistisen, D. Models and Methods for Hot Spot Safety Work. Ph.D. Thesis, Technical University of Denmark DTU, Lyngby, Denmark, 2002. [Google Scholar]
  26. Sweroad. General Directorate of Highways, A. Road Improvement and Traffic Safety Project: Black Spot Manual; Sweroad: Ankara, Turkey, 2001; pp. 1–82. [Google Scholar]
  27. Taylor, M.A.P.; Bonsall, P.W.; Young, W. Understanding Traffic Systems: Data Analysis and Presentation, 2nd ed.; Taylor and Francis Group: Oxford, UK, 2017; ISBN 9781315235370. [Google Scholar]
  28. Yuan, T.; Zeng, X.; Shi, T. Identifying Urban Road Black Spots with a Novel Method Based on the Firefly Clustering Algorithm and a Geographic Information System. Sustainability 2020, 12, 2091. [Google Scholar] [CrossRef] [Green Version]
  29. Gundogdu, I.B. Applying linear analysis methods to GIS-supported procedures for preventing traffic accidents: Case study of Konya. Saf. Sci. 2010, 48, 763–769. [Google Scholar] [CrossRef]
  30. Zhang, C.; Ivan, J.N. Effects of geometric characteristics on head-on crash incidence on two-lane roads in Connecticut. Transp. Res. Rec. 2005, 159–164. [Google Scholar] [CrossRef]
  31. Elvik, R. Comparative analysis of techniques for identifying locations of hazardous roads. Transp. Res. Rec. 2008, 72–75. [Google Scholar] [CrossRef]
  32. Sacchi, E.; Sayed, T.; El-Basyouny, K. Multivariate full Bayesian hot spot identification and ranking new technique. Transp. Res. Rec. 2015, 1–9. [Google Scholar] [CrossRef]
  33. Erdogan, S.; Yilmaz, I.; Baybura, T.; Gullu, M. Geographical information systems aided traffic accident analysis system case study: City of Afyonkarahisar. Accid. Anal. Prev. 2008, 40, 174–181. [Google Scholar] [CrossRef]
  34. Xie, Z.; Yan, J. Kernel Density Estimation of traffic accidents in a network space. Comput. Environ. Urban Syst. 2008, 32, 396–406. [Google Scholar] [CrossRef] [Green Version]
  35. Anselin, L. Local Indicators of Spatial Association—LISA. Geogr. Anal. 1995, 27, 93–115. [Google Scholar] [CrossRef]
  36. Ord, J.K.; Getis, A. Local Spatial Autocorrelation Statistics: Distributional Issues and an Application. Geogr. Anal. 1995, 27, 286–306. [Google Scholar] [CrossRef]
  37. Getis, A.; Ord, J.K. The analysis of spatial association by use of distance statistics. Adv. Spat. Sci. 2010, 61, 127–145. [Google Scholar] [CrossRef]
  38. Lee, M.; Khattak, A.J. Case Study of Crash Severity Spatial Pattern Identification in Hot Spot Analysis. Transp. Res. Rec. 2019. [Google Scholar] [CrossRef]
  39. Levine, N. CrimeStat II: A spatial Statistics Program for the Analysis of Crime Incident Locations, Part I.; The National Institute of Justice: Washington, DC, USA, 2002.
  40. Ziakopoulos, A.; Yannis, G. A review of spatial approaches in road safety. Accid. Anal. Prev. 2020, 135, 105323. [Google Scholar] [CrossRef]
  41. Fu, W.; Zhao, K.; Zhang, C.; Tunney, H. Using Moran’s I and geostatistics to identify spatial patterns of soil nutrients in two different long-term phosphorus-application plots. J. Plant Nutr. Soil Sci. 2011, 174, 785–798. [Google Scholar] [CrossRef]
  42. Cheng, Z.; Zu, Z.; Lu, J. Traffic Crash Evolution Characteristic Analysis and Spatiotemporal Hotspot Identification of Urban Road Intersections. Sustainability 2019, 11, 160. [Google Scholar] [CrossRef] [Green Version]
  43. Yamada, I.; Thill, J.C. Comparison of planar and network K-functions in traffic accident analysis. J. Transp. Geogr. 2004, 12, 149–158. [Google Scholar] [CrossRef]
  44. Plug, C.; Xia, J.; Caulfield, C. Spatial and temporal visualisation techniques for crash analysis. Accid. Anal. Prev. 2011, 43, 1937–1946. [Google Scholar] [CrossRef] [PubMed]
  45. Bíl, M.; Andrášik, R.; Janoška, Z. Identification of hazardous road locations of traffic accidents by means of kernel density estimation and cluster significance evaluation. Accid. Anal. Prev. 2013, 55, 265–273. [Google Scholar] [CrossRef] [PubMed]
  46. Bíl, M.; Andrášik, R.; Svoboda, T.; Sedoník, J. The KDE+ software: A tool for effective identification and ranking of animal-vehicle collision hotspots along networks. Landsc. Ecol. 2016, 31, 231–237. [Google Scholar] [CrossRef]
  47. Wong, D.W.S.; Lee, J. Statistical Analysis of Geographic Information with ArcView GIS and ArcGIS; Geographic Information Science; Wiley: Hoboken, NJ, USA, 2005; ISBN 9780471468998. [Google Scholar]
  48. Erdogan, S. Explorative spatial analysis of traffic accident statistics and road mortality among the provinces of Turkey. J. Safety Res. 2009, 40, 341–351. [Google Scholar] [CrossRef] [PubMed]
  49. Khan, G.; Qin, X.; Noyce, D.A. Spatial analysis of weather crash patterns. J. Transp. Eng. 2008, 134, 191–202. [Google Scholar] [CrossRef]
  50. Tola, A.M.; Gebissa, A. Identifying Black Spot Accident Zones using a Geographical Information System on Kombolcha-Dessie Road in Ethiopia. Int. J. Sci. Basic Appl. Res. 2019, 48, 1–12. [Google Scholar]
  51. Geurts, K.; Wets, G.; Brijs, T.; Vanhoof, K. Identification and ranking of black spots: Sensitivity analysis. Transp. Res. Rec. 2004, 1897, 34–42. [Google Scholar] [CrossRef]
  52. NSW Road Safety and Traffic Management Directorate. Road Traffic Accidents in New South Wales–1997-Statistical Statement: Year Ended 31 December 1997; Roads and Traffic Authority of NSW: Surry Hills, NSW, Australia, 1999. [Google Scholar]
  53. Mitra, S. Spatial autocorrelation and Bayesian spatial statistical method for analyzing intersections Prone to injury crashes. Transp. Res. Rec. 2009, 2136, 92–100. [Google Scholar] [CrossRef]
  54. Songchitruksa, P.; Zeng, X. Getis–Ord Spatial Statistics to Identify Hot Spots by Using Incident Management Data. Transp. Res. Rec. J. Transp. Res. Board 2010, 2165, 42–51. [Google Scholar] [CrossRef]
  55. Erdogan, S.; Ilçi, V.; Soysal, O.M.; Kormaz, A. A model suggestion for the determination of the traffic accident hotspots on the Turkish highway road network: A pilot study. Bol. Ciênc. Geodésicas 2015, 21, 169–188. [Google Scholar] [CrossRef]
  56. O’Sullivan, D.; Unwin, D.J. Geographic Information Analysis, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2010; ISBN 9780470288573. [Google Scholar]
  57. Maingi, J.K.; Mukeka, J.M.; Kyale, D.M.; Muasya, R.M. Spatiotemporal patterns of elephant poaching in South-Eastern Kenya. Wildl. Res. 2012, 39, 234–249. [Google Scholar] [CrossRef]
  58. Blank, L.; Cohen, Y.; Borenstein, M.; Shulhani, R.; Lofthouse, M.; Sofer, M.; Shtienberg, D. Variables Associated with Severity of Bacterial Canker and Wilt Caused by Clavibacter Michiganensis Subsp. Michiganensis in Tomato Greenhouses. Phytopathology 2016, 106, 254–261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. ERA (Ethiopian Roads Authority). Road Sector Development Program 21 Years Performance Assessment; ERA (Ethiopian Roads Authority): Addis Ababa, Ethiopia, 2019. [Google Scholar]
  60. Federal Transport Authority. Registered Vehicle Statistics; Federal Transport Authority: Addis Ababa, Ethiopia, 2019. [Google Scholar]
  61. ESRI Multi-Distance Spatial Cluster Analysis (Ripley’s k-Function). Available online: http://resources.esri.com/help/9.3/arcgisdesktop/com/gp_toolref/spatial_statistics_tools/multi_distance_spatial_cluster_analysis_ripley_s_k_function_spatial_statistics_.htm (accessed on 15 September 2021).
  62. World Bank. Addis Ababa, Ethiopia—Enhancing Urban Resilience (City Strength: Resilient Cities Program); World Bank: Washington, DC, USA, 2015. [Google Scholar]
  63. Ethiopian Roads Authority (ERA). Urban Transport Study and Preparation of Pilot Project of Addis Ababa Final Report; Ethiopian Roads Authority (ERA): Addis Ababa, Ethiopia, 2005. [Google Scholar]
  64. Samson, F. Analysis of Traffic Accident in Addis Ababa: Traffic Simulation. Master’s Thesis, Addis Ababa University, Addis Ababa, Ethiopia, 2006. [Google Scholar]
  65. Tulu, G.S.; Washington, S.; Haque, M.M.; King, M.J. Injury severity of pedestrians involved in road traffic crashes in Addis Ababa, Ethiopia. J. Transp. Saf. Secur. 2017, 9, 47–66. [Google Scholar] [CrossRef]
  66. World Bank Group. Scoping Study—Urban Mobility in Three Cities: Addis Ababa, Dar es Salaam, and Nairobi (English); World Bank Group: Washington, DC, USA, 2002. [Google Scholar]
  67. Oromia Traffic Police Bureau. Annual Crash Data Report; Oromia Traffic Police Bureau: Addis Ababa, Ethiopia, 2019. [Google Scholar]
  68. Quisumbing, A.R.; Yohannes, Y. How Fair is Workfare? Gender, Public Works, and Employment in Rural Ethiopia; Policy Research Working Papers; The World Bank: Washington, DC, USA, 2005. [Google Scholar]
  69. Al-Balbissi, A.H. Role of gender in road accidents. Traffic Inj. Prev. 2003, 4, 64–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Cordellieri, P.; Baralla, F.; Ferlazzo, F.; Sgalla, R.; Piccardi, L.; Giannini, A.M. Gender effects in young road users on road safety attitudes, behaviors and risk perception. Front. Psychol. 2016, 7, 1–11. [Google Scholar] [CrossRef] [Green Version]
  71. Martin, J.L.; Lafont, S.; Chiron, M.; Gadegbeku, B.; Laumon, B. Différences entre les hommes et les femmes face au risque routier. Rev. Epidemiol. Sante Publique 2004, 52, 357–367. [Google Scholar] [CrossRef]
  72. Liu, C.; Ye, T.J. Run-Off-Road Crashes: An On-Scene Perspective; National Highway Traffic Safety Administration: Washington, DC, USA, 2011; p. 20590.
  73. Kebede, A.; Gebresenbet, G. Mapping out goods flow to Addis Ababa city, Ethiopia, and it impact on environment. Transp. Res. Procedia 2017, 25, 1008–1020. [Google Scholar] [CrossRef]
Applsci 11 08828 g001 550
Figure 1. Map of the study area in Ethiopia.
Figure 1. Map of the study area in Ethiopia.
Applsci 11 08828 g001
Applsci 11 08828 g002 550
Figure 2. Number of deaths recorded in Ethiopia from 2010–2017.
Figure 2. Number of deaths recorded in Ethiopia from 2010–2017.
Applsci 11 08828 g002
Applsci 11 08828 g003 550
Figure 3. (a) Result of spatial autocorrelation; (b) Attribute table result of spatial autocorrelation.
Figure 3. (a) Result of spatial autocorrelation; (b) Attribute table result of spatial autocorrelation.
Applsci 11 08828 g003
Applsci 11 08828 g004 550
Figure 4. RTC Hot spot map of Ethiopia regions and federal territories.
Figure 4. RTC Hot spot map of Ethiopia regions and federal territories.
Applsci 11 08828 g004
Applsci 11 08828 g005 550
Figure 5. (a) Total lives lost due to RTCs; (b) Percentage deaths (Female versus Male).
Figure 5. (a) Total lives lost due to RTCs; (b) Percentage deaths (Female versus Male).
Applsci 11 08828 g005
Applsci 11 08828 g006 550
Figure 6. (a) Crashes occurred in a day of the week: from 2014–2017; (b) Crashes occurred in the time of the day: from 2014–2017.
Figure 6. (a) Crashes occurred in a day of the week: from 2014–2017; (b) Crashes occurred in the time of the day: from 2014–2017.
Applsci 11 08828 g006
Applsci 11 08828 g007 550
Figure 7. Crashes analyzed by collision types: 2014–2017.
Figure 7. Crashes analyzed by collision types: 2014–2017.
Applsci 11 08828 g007
Applsci 11 08828 g008 550
Figure 8. The spatial autocorrelation of Oromia in (a) Graph and (b) Attribute table.
Figure 8. The spatial autocorrelation of Oromia in (a) Graph and (b) Attribute table.
Applsci 11 08828 g008
Applsci 11 08828 g009 550
Figure 9. Crash hot spot map of Getis Ord Gi* statistics of Oromia region.
Figure 9. Crash hot spot map of Getis Ord Gi* statistics of Oromia region.
Applsci 11 08828 g009
Table
Table 1. Crash Hot Spot Result of Getis Ord Gi*.
Table 1. Crash Hot Spot Result of Getis Ord Gi*.
OBJ_IDCrash SeverityShape_L. (m)Shape_A. (m2)GiZScoreGiPValueNNeighborsRegion_Names
1109,518106,098525,638,4012.872300.004073Addis Ababa
830,6616,403,850355,423,484,2082.558480.010513Oromia
911,1922,806,015117,263,152,8670.856760.391583SNNP
219,8312,689,742153,443,579,1030.494820.620735Amhara
1182011,580,72556,451,528,629−0.078040.937803Tigrai
414771,630,48048,889,173,519−0.548550.583313B_Gumuz
6845908,66525,649,364,273−0.862380.388483Gambella
316931,920,34195,242,894,663−0.947920.343175Afar
1021793,677,531278,073,581,426−1.409980.158553Somali R.
52780215,4101,507,085,646−1.772870.076254Dire-Dawa
7820193,895394,011,903−1.772870.076254Harari
Table
Table 2. A four-year number of traffic crashes reported in their respective severity type.
Table 2. A four-year number of traffic crashes reported in their respective severity type.
YearFatalSerious InjurySlight InjuryPDOTotal
20141310901110017225033
2015135662351514463940
2016118856851816263900
2017131972953218064386
Percentage (%)29.9716.3515.4438.24100
PDO: Property damage only.
Table
Table 3. The Attribute of Crash Hot Spot Result of Getis Ord Gi*.
Table 3. The Attribute of Crash Hot Spot Result of Getis Ord Gi*.
OBJ_IDCrash SeverityShape L. (m)Shape_A. (m2)GiZScoreGiPValueNNeighborsZone_Names
715201,269,9489,892,678,0542.905730.003668E.Shewa
2042960,22685,906,9782.800420.005107Burayu T.
61501913,42911,530,834,3482.131660.033037N.Shewa
514881,577,74314,806,415,0712.013160.044107W.Shewa
1571225,99329,858,2602.013160.044107Adama T.
131154681,3476,508,288,0322.005480.044918S-W.Shewa
19149756,6088,097,272,7561.596920.110285H.G.Wellega
178851,080,19911,776,723,8201.312450.189374W.Arsi
1011551,390,03818,239,926,9530.800060.423672E.Hararge
9724768,63116,523,003,2040.739690.459494W.Hararge
8798909,81620,696,957,1540.413900.678957Arsi
46341,009,02618,075,624,315−0.192530.847335Jimma_Zone
2597980,90113,830,420,375−0.193090.846896E.Wellega
14649985,71418,577,054,735−0.434890.663643Guji
1630338,11850,520,944−0.727800.466743Jimma_City
122463,121,70745,463,584,611−0.828270.407522Borena
18352642,8909,851,170,119−0.836080.403113Kelem_Wellega
1507892,52312,744,967,754−0.869700.384464W.Wellega
3632999,48116,516,931,736−1.249010.211666I.A.Bora
112181,434,62144,912,392,310−1.495460.134793Bale
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tola, A.M.; Demissie, T.A.; Saathoff, F.; Gebissa, A. Severity, Spatial Pattern and Statistical Analysis of Road Traffic Crash Hot Spots in Ethiopia. Appl. Sci. 2021, 11, 8828. https://doi.org/10.3390/app11198828

AMA Style

Tola AM, Demissie TA, Saathoff F, Gebissa A. Severity, Spatial Pattern and Statistical Analysis of Road Traffic Crash Hot Spots in Ethiopia. Applied Sciences. 2021; 11(19):8828. https://doi.org/10.3390/app11198828

Chicago/Turabian Style

Tola, Alamirew Mulugeta, Tamene Adugna Demissie, Fokke Saathoff, and Alemayehu Gebissa. 2021. "Severity, Spatial Pattern and Statistical Analysis of Road Traffic Crash Hot Spots in Ethiopia" Applied Sciences 11, no. 19: 8828. https://doi.org/10.3390/app11198828

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