Department of geoinformatics, CDV – Transport Research Centre, Olomouc branch

Papers in international peer-reviewed journals

Articles in press

  • Bíl, M., Andrášik, R., Duľa, M., Sedoník, J., in press. On reliable identification of factors influencing wildlife-vehicle collisions along roads. Journal of Environmental Management.


  • Bíl, M., Andrášik, R., Sedoník, J., 2019. Which curves are dangerous? A network-wide analysis of traffic crash and infrastructure data. Transportation Research Part A 120C, 252–260.
    Show abstract |

    We conducted spatial analyses of traffic crashes, which took place in Czechia over 2010–2016, with respect to the road geometry data. The aim of the work was to identify hazardous road sub-segments where higher than expected numbers of traffic crashes occur.
    The entire Czech road network (58,200 km) was segmented at intersections into 39,074 between-intersection segments of varying lengths. Each road segment was further automatically sectioned, according to its horizontal alignment, into geometry-homogenous units (horizontal curves and tangents). Overall, 257,101 curves, defined as curved sections with radii below 2100 m, and 136,388 tangents, were identified. Subsequently, traffic crashes were joined to the respective geometrical units to determine their hazardousness. The degree of hazardousness was determined relatively, on a segment-by-segment basis, in order to eliminate the lack of precise traffic exposure data. In addition, the exact binomial test and Bayesian inference were used to identify the most hazardous horizontal curves.
    It was found that, in general, the curves with a higher crash risk have lower radii than the other curves. We identified the geographical locations of all curves with a high crash hazard. We also ranked the curves according to the crash hazard. Approximately ten percent of road segments contained at least one hazardous horizontal curve. 6943 significantly hazardous curves were identified by the use of the exact binomial test. The Bayesian inference reduced this number to 1395 (0.31% of the entire road network) and ranked them according to the Bayes factor. The most hazardous curve was 45 m long and contained 8.7 traffic crashes per year. Its hazard rate accounted for 37.4. This state-wide analysis of primary data was conducted over an extremely short time (up to 3 days) as the result of an application of an efficient algorithm for automatic road curvature determination.


  • Bíl M, Andrášik R, Sedoník J, Cícha V (2018) ROCA – An ArcGIS toolbox for road alignment identification and horizontal curve radii computation. PLoS ONE 13(12): e0208407.
    Show abstract | | Download

    We present the ROCA (ROad Curvature Analyst) software, in the form of an ESRI ArcGIS Toolbox, intended for vector line data processing. The software segments road network data into tangents and horizontal curves. Horizontal curve radii and azimuth of tangents are then automatically computed. Simultaneously, additional frequently used road section characteristics are calculated, such as the sinuosity of a road section (detour ratio), the number of turns along an individual road section and the average cumulative angle for a road section. The identification of curves is based on the naïve Bayes classifier and users are allowed to prepare their own training data files. We applied ROCA software to secondary roads within the Czech road network (9,980 km). The data processing took less than ten minutes. Approximately 43% of the road network in question consists of 42,752 horizontal curves. The ROCA software outperforms other existing automatic methods by 26% with respect to the percentage of correctly identified curves. The segmented secondary roads within the Czech road network can be viewed on the web-map application. We combined data on road geometry with road crashes database to develop the crash modification factors for horizontal curves with various radii. We determined that horizontal curves with radii of 50 m are approximately 3.7 times more hazardous than horizontal curves with radii accounting for 1000 m. ROCA software can be freely downloaded for noncommercial use from website.

  • Zahradníček, P, Münster, P, Bíl, M., Skalák, P., Štěpánek, P., Farda, A., Panský, M., Brzezina, J., Bílová, M., Kubeček, J., 2018. The December 2014 glaze event in the Czech Republic: predictability and impacts. Weather 73, 375–382.
  • Vodák, R., Bíl, M., Křivánková, Z., 2018. A modified ant colony optimization algorithm to increase the speed of the road network recovery process. International Journal of Disaster Risk Reduction 31C, 1092–1106.
    Show abstract | | Download

    When a disaster strikes many roads are blocked and the affected network may break up into a number of isolated parts. The reconnection of the network is therefore necessary for both relief distribution and planning of construction work. Shortening the time during which the road network is separated into isolated parts helps decrease indirect losses from disasters. The obstacles usually faced during the process of reconstruction include both the large number of blocked links and extensive affected areas (road networks).
    A reduction of the network into a much smaller complete graph and metaheuristic based on an ant colony optimization has been introduced to overcome this issue. We demonstrate that, for small networks, the metaheuristic produces the same results as other deterministic algorithms. We further show that the method is still a viable approach for large networks (723 nodes and 974 links, where we artificially blocked 46 links) when the NP-hard nature of this problem began to affect the computational time of the deterministic algorithms.
    We demonstrate how the various scenarios can be included into the algorithm. We finally introduce a new ranking of feasible solutions which enables the algorithm to minimize the time of reconstructions for all repair units. Reasonable results were obtained after five minutes of computation. There is nevertheless an up-to-38% improvement of the initial solution. The algorithm can also be used for both relief distribution, when no roads were damaged, and for planning of construction work when damaged roads occur.

  • Favilli, F., Bíl, M., Sedoník, J., Andrášik, R., Kasal, P., Agreiter, A., Streifeneder, T., 2018. Application of KDE+ software to identify collective risk hotspots of ungulate-vehicle collisions in South Tyrol, Northern Italy. European Journal of Wildlife Research 64:59.
    Show abstract |

    This paper is the first dealing with animal-vehicle collisions (AVC) with red and roe deer in South Tyrol, Northern Italy. The Autonomous Province of Bolzano (South Tyrol) has been collecting AVC data since 2012 on the entire provincial road network. Each year, AVC data accounted for more than 700 cases per year, with several socioeconomic and ecological implications. The aim of this research is to identify the locations where AVC occur more frequently than expected (hotspots) and better outline subsequent implementation of mitigation measures. For an effective identification of AVC hotspots, we applied a combined methodology of temporal and spatial analysis on AVC data collected on the South Tyrol road network in the years 2012–2014. AVC data enabled the identification of the temporal patterns, which showed different behaviors of the two target species in close proximity of the road network and throughout the 12 months. The KDE+ software applied to the 2012–2014 AVC database allowed for spatial analysis and the identification of hotspots, i.e., the road sections having the highest risk for drivers. The integration of the results, coming from the abovementioned methodologies, contributes to a detailed assessment of roads that would allow the identification of the local contributing factors and a base-line of potential problematic areas that will highlight the need for further investigation to assess whether the risk-rank is accurate and allocate effectively limited resources to a feasible number of identified hotspots and reduce the current degree of AVC in the South Tyrolean road network.

  • Bartonička, T., Andrášik, R., Duľa, M., Sedoník, J., Bíl, M., 2018. Identification of Local Factors Causing Clustering of Animal-Vehicle Collisions. Journal of Wildlife Management 82, pp. 940–947
    Show abstract |

    Effective measures reducing risk of animal-vehicle collisions (AVC) require defining high-risk locations on roads where AVCs occur. Previous studies examined factors explaining locations of individua l AVCs; however, some AVCs can form hotspots (i.e., clusters of AVCs) that can be explained by local factors. We therefore applied a novel kernel density estimation (KDE) method to AVCs for the Czech Republic from October 2006 to December 2011 to identify AVCs hotspots along roads. Our main goal was to identify local factors and their effect on the non-random (clustered) occurrenc e of AVCs. The remaining solitary AVCs occurred randomly and are likely induced by other human factors on the global scale. The hotspot identification method followed by the selected data mining methods (KDE þ methods) identified factors causing local clustering of AVCs. Distance from forest ( < 350 m) or linear vegetation were important factors for estimating presence of clusters of AVCs; in open areas, AVC clusters were absent. Further research on effectiveness of measures reducing risk of AVC should focus on clusters of AVCs, not on the individua l AVC. We recommend that state transportati on agencies focus mitigation actions in forested areas.

  • Bíl, M., Dobiáš, M., Andrášik, R., Bílová, M., Hejna, P., 2018. Cycling Fatalities: When A Helmet is Useless and when it Might Save Your Life. Safety Science 105C, 71–76.
    Show abstract | | Download

    Autopsy reports of 119 cyclists who died in two Czech regions between 1995 and 2013 as a result of traffic crashes were studied. In all the study cases, pathologists analyzed whether a helmet could have helped the cyclists survive the crash or not. The crash circumstances from the police reports were then evaluated.
    The results indicate that helmets could have helped the most in cases of single-vehicle crashes when cyclists fell off their bicycles or hit obstacles and in certain cases when an intracranial injury was the primary cause of death. Altogether 44 cyclists (37%) from this study could have survived if they had been wearing helmets during the crashes.
    Helmets would not have helped cyclists in most high-energetic crashes, especially when motor-vehicles or trains were involved. Some rear-end crashes outside urban areas also resulted in injuries when a helmet would not have helped.
    This study concludes that cyclists should wear helmets, but they should also be aware that it cannot protect them in particular situations. These facts should be incorporated into safety campaigns to prevent cyclists from feeling protected in such situations when helmets cannot help. Our results also support the building of cycling paths separate from traffic, particularly outside of urban areas.

  • Bíl, M., Andrášik, R., Bartonička, T., Křivánková, Z., Sedoník, J., 2018. An Evaluation of Odor Repellent Effectiveness in Prevention of Wildlife-Vehicle Collisions. Journal of Environmental Management 205C pp. 209–214.
    Show abstract |

    Wildlife-vehicle collisions (WVC) amount to 11 % of all registered traffic crashes in the Czech Republic causing, apart from numerous deaths and serious injuries to animals, property damage and injuries to car passengers. Odor repellents have the potential to lower the overall number of WVC and allow animals to cross roads at the same time.
    We tested the effectiveness of odor repellent preparation in prevention of WVC. 18 places were selected on the Czech road network where WVC were concentrated on the basis of traffic crash data. Control sections on the same road segments were also delimited in order to keep the traffic intensities constant.
    We applied a Before-After-Control-Impact (BACI) study design to control not only the effect of the measures but also the expected natural variations in wildlife populations over time. Data were compared before and after odor repellent installations. Wildlife carcass gathering was carried out during the spring and autumn. We also used the police crash database to supplement carcass data when no field works were carried out. 201 killed mammals (roe deer and wild boars) were identified in total over 47 months.
    We applied a Bayesian approach as only a limited numbers of WVC were available. A WVC decrease between 26 – 43 % can be expected on the treated road sections. These numbers are, however, up to three-times lower than those claimed by producers of odor preparations.


  • Bíl, M., Andrášik. R., Nezval, V., Bílová, M., 2017. Identifying Locations along Railway Networks with the Highest Tree Fall Hazard. Applied Geography 87, 45–53.
    Show abstract |

    Disruptions of railway traffic have many reasons. Tree falls onto railway tracks or overhead lines rank among the most common causes of disruptions of a natural origin. 2039 tree-fall events, containing up to 70 individual trees per event, were registered on the Czech railway network between 2012 and 2015. 32% of them were directly caused by 14 weather extremes during which more than 20 concurrent tree-fall events were registered. Moreover, 12 train derailments due to fallen trees were registered on Czech railways within the same period.
    We combined land use data along railway tracks and data on tree falls. Land use and railway tracks data were obtained from a freely available Open Street Map database. The tree fall hazard was then computed using empirical data, data on land use and a generalized rule of succession. The clustering approach was also applied to focus on localities where tree falls were concentrated regardless of the resulting segment hazard. There were 59 rail track segments (out of 2960) with the highest tree fall hazard and 267 clusters were finally identified. The clusters and the most hazardous railway segments will be among the first in the process of line side vegetation monitoring in order to minimize potential losses from tree fall. The presented method can be widely applicable elsewhere.

  • Bíl, M., Kubeček, J., Sedoník, J., Andrášik. R., 2017. A system for evidence of animal-vehicle collisions along transportation networks. Biological Conservation 213PA, pp. 167–174.
    Show abstract |

    We present a system for state-wide evidence of animal-vehicle collisions (AVC). The primary part of this system is a geographic database which is connected to a web-map application. AVC data come from the Police via an online system of traffic incidents (JSDI) and from volunteers through a web or mobile interface. Data are processed using automatic scripts which identify data errors and perform spatial analyses.
    The application automatically computes AVC hotspots every midnight and crash densities along road sections. Hunter area administrators consequently have an overview of their areas. More than 40,000 records are currently included in this database. 50% of them were added over the last two years when it was launched. The majority of data (90%) came from JSDI. The species is known for 44% of JSDI records. The majority of the identified species were roe deer (75%), followed by wild boar (15%).
    Roe deer crashes occur most frequently in May within 2 h after sunset. 32.5% (56.4%) of these crashes occur within 1 h (2 h) before or after sunset or sunrise. For wild boar, the values are less distinctive (19.4% and 37.7%). Approximately 1800 AVC hotspots, which cover 0.5% of the Czech road network, were detected and visualized on a map.


  • Andrášik, R.; Bíl, M., 2016. Efficient Road Geometry Identification from Digital Vector Data. Journal of Geographical Systems 18(3), 249–264. DOI 10.1007/s10109-016-0230-1
    Show abstract |

    A new method for the automatic identification of road geometry from digital vector data is presented. The method is capable of efficiently identifying circular curves with their radii and tangents (straight sections). The average error of identification ranged from 0.01 to 1.30 % for precisely drawn data and 4.81 % in the case of actual road data with noise in the location of vertices. The results demonstrate that the proposed method is faster and more precise than commonly used techniques. This approach can be used by road administrators to complete their databases with information concerning the geometry of roads. It can also be utilized by transport engineers or traffic safety analysts to investigate the possible dependence of traffic accidents on road geometries. The method presented is applicable as well to railroads and rivers or other line features.

  • Bíl, M., Andrášik, R., Zahradníček, P., Kubeček, J., Sedoník, J., Štěpánek, P., 2016. Total water content thresholds for shallow landslides, Outer Western Carpathians. Landslides 13, 337–347.
    Show abstract |

    Shallow landslides are fairly frequent natural processes which emerge as a result of both rainfall and rapid snowmelt in the Flysch Belt of the Outer Western Carpathians. We estimated the total water content thresholds for the previously defined seven phases of increased landsliding which took place between 1939 and 2010 around the Napajedla meteorological station. The time series were reconstructed on the basis of data from surrounding stations. Rainfalls with the highest intensities (>1 mm/min) were removed from the set. Rainfall of such an intensity primarily causes overland flow and soil erosion and does not contribute to landslide threshold. The snow water equivalent was computed on the basis of the snow height, and possible errors were evaluated as interval estimations. An interval of 10 days before a landslide phase was selected for the total water content threshold. The resulting lower boundary (67.0 mm/10 days) and upper boundary (163.3 mm/10 days) thresholds of water infiltrated into soil during an event shall be part of the prepared online warning system in this area.

  • Bíl, M., Bílová, M., Dobiáš, M., Andrášik, R., 2016. Circumstances and Causes of Fatal Cycling Accidents in the Czech Republic. Traffic Injury Prevention 17 (4), 394–399. doi:10.1080/15389588.2015.1094183.
    Show abstract |

    Objectives: The circumstances and causes of death of 129 cyclists registered in the Olomouc and the Zlín regions, the Czech Republic, between 2005 and 2013 were the subject of this study.
    Methods: We analyzed the autopsy reports, where the principal cause of death was stated, and obtained a detailed description of the circumstances recorded by the police officers.
    Results: Eighty-three cases (64.3% of the set) were collisions involving a motor vehicle. The driver was the guilty party in 57 cases (68.7%) and the cyclist in the remaining 26 cases (31.3%). The most frequent cause of the crash was connected with right of way (29 cases). Cars were involved in 52 cases; heavy vehicles, including buses, in 26 cases; and motorcycles in 5 cases. Single-vehicle crashes consisted of 43 (33.3%) cases. We divided this group into 3 subgroups based on whether the particular case could be attributed to a cyclist having lost control of the bicycle (31 cases) or to other particular causes. Sixty-eight cases (52.7%) of fatal outcomes were directly linked to intracranial injuries. Multiple injuries were the principal cause of death in 19 cases (14.7%), followed by hemorrhagic traumatic shock (12 cases, 9.3%). Seventy-two (55.8%) cyclists died immediately after the crash and 23 (17.8%) cyclists died within a day of the accident.
    Conclusions: Trucks were more dangerous to cyclists than cars at intersections, whereas cars were more dangerous on straight sections. The most important pattern was identified as a motor vehicle hitting a cyclist from behind on a straight road section. We identified a strong underestimation of natural death as a cause of cycling fatalities in the official police reports.

  • Bíl, M., Andrášik, R., Svoboda, T., Sedoník, J., 2016. The KDE+ software: a tool for effective identification and ranking of animal-vehicle collision hotspots along networks. Landscape Ecology 31, 231–237.
    Show abstract |

    Context: Objective identification of locations on transportation networks, where animal-vehicle collisions (AVC) occur more frequently than expected (hotspots), is an important step for the effective application of mitigation measures.
    Objectives: We introduce the KDE+ software which is a programmed version of the KDE+ method for effective identification of traffic accident hotspots. The software can be used in order to analyze animal-vehicle collision data.
    Methods: The KDE+ method is based on principles of Kernel Density Estimation (KDE). The symbol ‘+’ indicates that the method allows for the objective selection of significant clusters and for the ranking of the hotspots. It is also simultaneously applicable to an unlimited number of road segments.
    Results: We applied the KDE+ method to the entire Czech road network. The hotspots were ranked according to their significance. The resulting hotspots represent a short overall road length which should require a more detailed assessment in the field. The 100 most important clusters of AVC represent, for example, only 19.7 km of the entire road network (37,469 km).
    Conclusions: We present an objective method for hotspots identification which can be used for AVC data. This method is unique because it determines the significance level of hotspots in an objective way. The prioritization of hotspots allows a transportation manager to effectively allocate resources to a feasible number of identified hotspots. We describe the software, data preparation and present the KDE+ application to AVC data.


  • Bíl, M., Andrášik, R., Kubeček, J., 2015. How comfortable are your cycling tracks? A new method for objective bicycle vibration measurement. Transportation Research Part C: Emerging Technologies 56, 415–425.
    Show abstract |

    Cycling comfort consists of several factors. Their relevant values are important in the process of bicycle facility planning. Poor surface pavement quality manifests itself in terms of vibrations of a bicycle. This strongly influences the perception of a cycle track, general cycling comfort and the route choice as well. We introduce dynamic comfort index (DCI) which is capable of objectively describing the vibration properties of surface pavement on a track. The DCI is derived from data gathered when riding a bicycle equipped with a GPS device and an accelerometer. The most common types of devices were selected to make the DCI widely applicable. We tested DCI values on various bicycles and surface pavements. DCI values on individual cycling tracks were compared with the subjective feelings of 43 cyclists via questionnaires. A strong correlation (−0.94) was obtained between the objectively measured DCI values and the subjectively assessed evaluations. This makes the DCI approach transferable to any other environment. This method has been applied to an entire road network within the historical center of the city of Olomouc (Czech Republic). It can further be used by bicycle track administrators to monitor surface quality, by planners to obtain relevant surface pavement values, and by individual cyclists for optimal route choice.

  • Bíl, M., Vodák, R., Kubeček, J., Bílová, M., Sedoník, J., 2015. Evaluating Road Network Damage Caused by Natural Disasters in the Czech Republic between 1997 and 2010. Transportation Research Part A: Policy and Practice 80, 90–103.
    Show abstract |

    Road networks play a vital role in maintaining a functioning modern society. Many events perceptibly affect the transport supply along these networks, especially natural disasters such as floods, landslides, and earthquakes. Contrary to more common disruptions of traffic from accidents, or maintenance closures, natural disasters are capable of destroying large numbers of roads and usually cover vast areas. When evaluating network damage no single measure alone is able to describe the full extent of network destruction. In this study, we investigated six highly damaging natural disasters, which occurred in the Czech Republic between 1997 and 2010. They were all induced by extreme rainfall or by rapid snowmelt and resulted in floods and landslides. Their impacts are evaluated with respect to the damage to road networks and decreased serviceability. For mutual comparison of the impacts and their analysis we used several criteria, described in the paper, related to economic impacts, physical harm to individuals and infrastructures, and the effects on connectivity and serviceability. We also introduced a new measure based on the network efficiency index which takes into account the importance of nodes based on their population. Moreover, we provide a detailed analysis of one such event in July 1997 that significantly affected the road network of the Zlín region.

  • Vodák, R., Bíl, M., Sedoník, J., 2015. Network robustness and random processes. Physica A 428, 368–382.
    Show abstract |

    We introduce two new measures of network robustness and apply them to four different strategies. The measures are independent from the number of nodes in the network and have the strong potential to cover a large portfolio of applications. Using the Monte-Carlo methods, we demonstrate how to approximate the measures. The methods are based on random interruption of links with suitable constraints which represent the above-mentioned strategies. We introduce two networks with obvious varying robustness to demonstrate the measures. We also demonstrate how to employ the measures in order to improve the robustness of the networks by adding one new link. We further indicate that the measures are able to identify the infrequently connected parts of the network and suggest the most appropriate improvement. We also discuss the consequences of the obtained results and the possible applications of the measures.


  • Bíl, M., Krejčí, O., Bílová, M., Kubeček, J., Sedoník, J., Krejčí V., 2014. A Chronology of Landsliding and its Impacts on the Village of Halenkovice, Outer Western Carpathians, Geography 119 (4), 342–363.
    Show abstract |

    The article is based on an investigation of landsliding chronology coducted at the village of Halenkovice. (Outer Western Carpathians, Czechia). On the basis of historical data, (chronicles and other archive sources, air photos, old maps), field mapping and interviews of eyewitnesses, we determined six major and seven minor phases of landsliding in the village and its immediate neighborhood for the period 1915–2010. Inactive and active landslides currently cover 20% of the Halenkovice cadastral area. Landslides have strongly affected the development of the village itself. Since 1941, at least 9 houses have been destroyed by landslides, with many other local buildings and roads suffering significant damage. We also documented two examples of periodic reactivation of landslides.

  • Baroň, I, Bíl, M., Bábek, O., Smolková, V., Pánek, T., Macur, L., 2014. Effect of slope failures on river-network pattern: A river piracy case study from the Flysch Belt of the Outer Western Carpathians. Geomorphology 214, 356–365.
    Show abstract |

    Landslides are important geomorphic agents in various mountainous settings. We document here a case of river piracy from the upper part of the Malá Brodská Valley in the Vsetínské Mts., Czech Republic (Rača Unit of the flysch Magura Group of Nappes, flysch belt of the Outer Western Carpathians) controlled by mass movement processes. Based on the field geological, geomorphological and geophysical data, we found out that the landslide accumulations pushed the more active river of out of two subparallel river channels with different erosion activity westwards and forced intensive lateral erosion towards the recently abandoned valley. Apart from the landslide processes, the presence of the N-striking fault, accentuated by higher flow rates of the eastern channel as a result of its larger catchment area, were the most critical factors of the river piracy. As a consequence of the river piracy, intensive retrograde erosion in the elbow of capture and also within the upper portion of the western catchment occurred. Deposits of two landslide dams document recent minimum erosion rates to be 18.8− 1 in the western (captured) catchment, and 3.6− 1 in the eastern catchment respectively. The maximum age of the river piracy is estimated to be of the late Glacial and/or the early Holocene.

  • Bíl, M., Kubeček, J., Andrášik, R., 2014. An Epidemiological Approach to Determining the Risk of Road Damage due to Landslides. Nat Haz. 73 (4), 1323–1335.
    Show abstract |

    Disruption of segments of roads can have a significant impact on the vulnerability of the entire network. Natural disasters are frequent causes of disruptions of this kind. This article focuses on determining the risk of road disruptions due to landslides. Our approach is based on methodology widely used in the field of epidemiology. We had available data on the location of the landslides, the road network and a list of the disrupted road segments. With the use of a 2 × 2 table, we determined the relationship between landslide data and road segment disruptions and derived the risk coefficient based on the number of landslides in the vicinity of the road and its length. The result is a disruption risk map with risk coefficients ranging from 0 to 47.94. In order to distinguish the most risky segments, we calculated a threshold of 12.40 with the use of a risk breakdown in a group of segments without damage. Nineteen percentage (402 km) of the road network in the Zlín region (Czech Republic), where the methodology was applied, is located beyond this threshold. The benefits of this approach stem from its speed and potential to define the most risky areas on which a detailed geomorphologic analysis can be focused.


  • Bíl, M., Andrášik, R., Janoška, Z., 2013. Identification of Hazardous Road Locations of Traffic Accidents by means of Kernel Density Estimation and Cluster Significance Evaluation. Accident Analysis and Prevention 55, 265–273.
    Show abstract |

    This paper proposes a procedure which evaluates clusters of traffic accident and organizes them according to their significance. The standard kernel density estimation was extended by statistical significance testing of the resulting clusters of the traffic accidents. This allowed us to identify the most important clusters within each section. They represent places where the kernel density function exceeds the significance level corresponding to the 95th percentile level, which is estimated using the Monte Carlo simulations. To show only the most important clusters within a set of sections, we introduced the cluster strength and cluster stability evaluation procedures. The method was applied in the Southern Moravia Region of the Czech Republic.


  • Bíl, M., Kubeček, J., 2012.  Piping in loess-like and loess-derived soils: case study of Halenkovice site, Czech Republic. Annales Societatis Geologorum Poloniae 82, 45–50.
    Show abstract |

    The soil piping that occurs on luvisols in the vicinity of the village of Halenkovice was studied for 5 years. These piping phenomena can only be found where arable land meets the forest or a belt of shrubbery. If there is a scarp in the locality, which usually changes from 6° in the field to approximately 30° in the forest, soil pipes are more likely to occur. Before the scarp, the slope flattens out and it is almost horizontal. This factor makes it possible for the overland flow to seep into the slope. This seepage results in soil piping, which is formed in loess loam and colluvial deposits. There are about 15 sites in the vicinity of the village of Halenkovice, where soil piping occurs. In one of them, Halenkovice 1 (an area of 900 m2) we closely studied 47 partial cavities. Their internal volume is 3.8 m3. The volume of the sink holes is 23 m3. There are two types of soil pipes – vertical, which on average tend to be shorter (40 cm) and lead the water under the surface, and soil pipes parallel with the slope, which are on average 81 cm long. Water flows through the pipes during a thaw or precipitation, which often takes away the top soil. The intensity of this process depends on the intensity of precipitation, which occurs outside the growing season, when there are no crops in the fields.

  • Bíl, M., Bílová, M., Kubeček, J., 2012.  Unified GIS database on cycle tourism infrastructure. Tourism Management 33, 1554–1561.
    Show abstract |

    Bicycles are used in the Czech Republic for commuting to work and for leisure time activities. This is reflected in the cycle trail administrators' offer to make the existing network denser, design new routes, mark their courses and install a complementary cycle infrastructure.
    However, extensive growth of these activities in the last years has led to the loss of overview information on the overall cycle trail network. That is the reason for producing a methodology of capturing and representing the information. A unified GIS database on the cycle infrastructure (UDCI) was created and includes a data collection system with the use of GPS, the coding of descriptive information on cycle trail segments and the administration of GIS layers in a topologic data model. The methodology of the UDCI application is demonstrated with a specific example of a cycle trail network in the South Bohemian region.


  • Bíl, M., Bílová, M., Müller, I., 2010. Critical factors in fatal collisions of adult cyclists with automobiles. Accident Analysis and Prevention 42, 1632–1636.
    Show abstract |

    This article evaluates, by means of multivariate regression, critical factors influencing the collisions of motor vehicles with adult (over 17 years) cyclists that result in fatal injury of cyclists. The analysis is based on the database of the Traffic Police of Czech Republic from the time period 1995–2007. The results suggest that the most consequential categories of factors under study are: inappropriate driving speed of automobile; the head-on crash; and night-time traffic in places without streetlights. The cyclists’ faults are of most serious consequence on crossroads when cyclists deny the right of way. Males are more likely to suffer a fatal injury due to a collision with a car than females. The most vulnerable age group are cyclists above 65 years. A fatal injury of a cyclist is more often driver's fault than cyclist's (598 vs. 370).
    In order to reduce the fatal risk, it is recommended to separate the road traffic of motor vehicles from bicyclists in critical road-sections; or, at least, to reduce speed limits there.

Papers at international peer-reviewed conferences

  • Bíl, M., Andrášik, R., Kubeček, J., Křivánková, Z., Vodák, R., 2017. RUPOK: An Online Landslide Risk Tool for Road Networks. In: Mikoš M., Vilímek V., Yin Y., Sassa K. (eds) Advancing Culture of Living with Landslides. WLF 2017. Springer, Cham, pp 19–26.
  • Andrášik, R., Bíl, M., Slovák, R., 2016. How (not) to work with small probabilities: Evaluating the individual risk of railway transport. Risk, Reliability and Safety: Innovating Theory and Practice – Proceedings of the 26th European Safety and Reliability Conference, ESREL 2016. 672–676
  • Vodák, R., Bíl, M., Křivánková, Z., Andrášik, R., Kubeček, R., 2016. Hidden risks in electric grids due to dependency on transportation networks. Risk, Reliability and Safety: Innovating Theory and Practice – Proceedings of the 26th European Safety and Reliability Conference, ESREL 2016. 1439–1442.
  • Andrášik, R., Bíl, M., 2016. Traffic accidents hotspots: Identifying the boundary between the signal and the noise. Risk, Reliability and Safety: Innovating Theory and Practice – Proceedings of the 26th European Safety and Reliability Conference, ESREL 2016. 1634–1637.
  • Andrášik, R., Bíl, M., 2015. Traffic accidents: Random or pattern occurrence? Safety and Reliability of Complex Engineered Systems. Podofiliny et al. (Eds.), ISBN 978-1-138-02879-1
  • Vodák, R., Andrášik, R., Bíl, M., Sedoník, J., 2015. Network robustness analysis based on current road incident data. Safety and Reliability of Complex Engineered Systems – Proceedings of the 25th European Safety and Reliability Conference, ESREL 2015.
  • Bíl, M., Vodák, R., 2014. The stochastic approach in road network vulnerability analysis. Safety and Reliability: Methodology and Applications – Nowakowski et al. (Eds), 929–932.