Roadkill-Data-Based Identiﬁcation and Ranking of Mammal Habitats

: Wildlife–vehicle collisions, as well as environmental factors that affect collisions and mitigation measures, are usually modelled and analysed in the vicinity of or within roads, while habitat attractiveness to wildlife along with risk to drivers remain mostly underestimated. The main goal of this study was the identiﬁcation, characterisation, and ranking of mammalian habitats in Lithuania in relation to 2002–2017 roadkill data. We identiﬁed habitat patches as areas (varying from 1 to 1488 square kilometres) isolated by neighbouring roads characterised by at least one wildlife– vehicle collision hotspot. We ranked all identiﬁed habitats on the basis of land cover, the presence of an ecological corridor, a mammalian pathway, and roadkill hotspot data. A ranking scenario describing both habitat attractiveness to wildlife and the risk to drivers was deﬁned and applied. Ranks for each habitat were calculated using multiple criteria spatial decision support techniques. Multiple regression analyses were used to identify the relationship between habitat ranks, species richness, and land cover classes. Strong relationships were identiﬁed and are discussed between the habitat patch ranks in ﬁve (out of 28) land cover classes and in eight (out of 28) species (97% of all mammal road kills). We conclude that, along with conventional roadkill hotspot identiﬁcation, roadkill-based habitat identiﬁcation and characterisation as well as species richness analysis should be used in road safety infrastructure planning.


Introduction
Wildlife-vehicle collisions (WVCs) pose a threat to human life and biological diversity and result in damage to property [1][2][3][4][5][6]. Over the last two decades in Lithuania, while the overall number of road traffic accidents has continuously decreased, road accidents involving wildlife have increased [7].
To mitigate mammal-vehicle collisions (MVCs), fencing, underpasses, gates, and jump-out ramps are used as the most common mitigation measures in the country [8]. Additional road safety infrastructure elements such as repellents, reflectors, noise, and natural predators can also be used; these focus on a single and/or multiple wildlife species. They repel, attract, or redirect wildlife with different ecological and financial efficiencies [9][10][11][12][13][14][15][16][17][18]. The selection of tangible multi-scale [19], multi-objective, and multifunctional WVC mitigation measures is the focus of a considerable research challenge [20].
The identification of roadkill hotspots (road sections where collisions occur more frequently than expected) is the first step of the highway safety management process. However, erroneous hotspot identification [21] as well as gaps in roadkill data [1] may result in inefficient use of resources for safety improvements [22]. There are many generalised linear models [23] that can be used to identify hotspots, such as ecological niche modelling [24], kernel density estimation [25][26][27], distance-based approaches [28], and methods based on modelling the number of collisions in a road section assuming a Poisson distribution [21,[29][30][31][32]. These methods use roadkill data to detect collision hotspots as well as their risk to drivers. In order to assess habitat attractiveness to wildlife and the associated habitat risk to drivers, it is important to understand where mammals cross roads more frequently.
Field research usually brings disparate results [47] of differential scale and quality [1]. Consequently, the results are frequently not fit for deriving habitat patch characteristics and assessing habitat attractiveness to wildlife. This would require standardised habitat data that are usually lacking.
Habitats can be characterised using behavioural and spatiotemporal events, landscape connectivity and fragmentation, species richness, animal abundance, and other field research data. Large scale, long-term, and accurate data that can characterise habitats usually require methodologically robust and expensive research. Employing the available roadkill data from police reports would decrease (not replace) the amount of field research required in cases when there is insufficient habitat data available. Multiple, long-term, and standardised habitat characteristics (criteria) can help us to identify MVC mitigation measures focused on single or multiple species.
Decision-makers often deal with problems that involve multiple criteria [48][49][50][51]. Identification of the primary sources of MVCs, namely the habitats that are highly attractive to wildlife and simultaneously of high risk to drivers, is also a multiple criteria analysis problem. Therefore, we selected Simple Additive Weighting (SAW) and Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS) [49,52] multi-criteria spatial decision support techniques for the ranking of habitats. The habitat ranking outcomes can be considered reliable if both methods generate similar ranking results [8].
Wildlife-vehicle collisions, as well as environmental factors that affect collisions and mitigation measures, are usually modelled and analysed at the level of the roads themselves [7,8,[53][54][55][56], while the wider issues of adjacent habitat attractiveness to wildlife and its risk to drivers remain underestimated. MVCs with wild species accounted for about 91% of all WVCs in Lithuania in 2002-2017 [57]. There is a need for a framework that helps us to unify the disparate results emerging from different data sources and field studies on MVC occurrence, which allows for roadkill-based identification, characterisation, and multi-objective ranking of mammalian habitats by their attractiveness to wildlife and their risk to drivers. Here, we understand "risk to the driver" as a derivative of the cluster strength in KDE+ [26,27].
In this study, habitat identification, characterisation, and ranking are based on the definition of habitats as "areas isolated by neighbouring roads that have at least one hotspot (a road section where MVCs occur more frequently than expected)" and the assumptions that (1) highly attractive habitats for wildlife increase the risk of MVCs on adjacent roads; (2) habitats that are surrounded by roads with an abundance of MVC hotspots are of high attractiveness to wildlife movement; and (3) road kills in the hotspots can help us to identify species richness within adjacent habitats. However, the accuracy of such estimations depends on the completeness of MVC data [1].
The overall purpose of this study, therefore, is to: • Identify habitat patches that are surrounded by roads with kernel density estimation (KDE+)-based [27] MVC hotspots; • Characterise habitat patches using the properties of adjacent habitats, hypothetical corridors and wildlife pathways, hotspots, and land cover data; • Define ranking scenarios (criteria utility functions and criteria weights) to detect habitat patches that are highly attractive to wildlife [37] and pose a risk to drivers [27]; • Rank habitat patches using two different multiple criteria spatial decision support techniques: SAW and TOPSIS; and • Find relationships between habitat ranks, species richness, and land cover classes for use in the planning of multispecies MVC mitigation measures using multiple linear regressions.

Study Area
Our study area covers the entire territory of Lithuania (Figure 1), which can be characterised as mostly a plain. It represents a surface area of 65,286 square kilometres. In 2012, 33% of the surface was covered by arable land and permanent crops, 27% by semi-natural vegetation, 33% by forested land, 3% by developed (artificial) areas, and 4% by water bodies and other land. The land cover change (a 0.48% change rate per year) in the country is slowing, mainly due to a rapid decrease in the intensity of forest conversion [58].
Land 2021, 10, x FOR PEER REVIEW 3 of 33 • Rank habitat patches using two different multiple criteria spatial decision support techniques: SAW and TOPSIS; and • Find relationships between habitat ranks, species richness, and land cover classes for use in the planning of multispecies MVC mitigation measures using multiple linear regressions.

Study Area
Our study area covers the entire territory of Lithuania (Figure 1), which can be characterised as mostly a plain. It represents a surface area of 65,286 square kilometres. In 2012, 33% of the surface was covered by arable land and permanent crops, 27% by semi-natural vegetation, 33% by forested land, 3% by developed (artificial) areas, and 4% by water bodies and other land. The land cover change (a 0.48% change rate per year) in the country is slowing, mainly due to a rapid decrease in the intensity of forest conversion [58]. In 2017, there were 21,244 km of State-owned roads of national significance (excluding roads in cities): 1751 km of main roads; 4925 km of national roads; and 14,568 km of regional roads [59]. In this study, we analysed 1784 roads (21 main/highway, 13 national, and 1631 regional) as shown in Figure 1.
In the period 2002-2017, the annual average daily traffic (AADT) increased from 5600 to 11,000 vehicles a day on main roads, from 2200 to 2900 vehicles a day on national roads, and remained at up to 500 vehicles a day on regional roads [60]. In 2017, there were 21,244 km of State-owned roads of national significance (excluding roads in cities): 1751 km of main roads; 4925 km of national roads; and 14,568 km of regional roads [59]. In this study, we analysed 1784 roads (21 main/highway, 13 national, and 1631 regional) as shown in Figure 1.

Mammal-Vehicle Collision Data
In the period 2002-2017, the annual average daily traffic (AADT) increased from 5600 to 11,000 vehicles a day on main roads, from 2200 to 2900 vehicles a day on national roads, and remained at up to 500 vehicles a day on regional roads [60].

Mammal-Vehicle Collision Data
According to the data from the Lithuanian Police Traffic Supervision Service and the Nature Research Centre, a total of 24,083 WVCs were recorded over the period 2002-2017 in Lithuania [57]. These numbers may, however, have a bias regarding taxonomic groups and not account for all accidents as reporting to the authorities is not mandatory in Lithuania. The Traffic Supervision Service registers only road kills from those accidents that were reported by drivers; therefore, their data are biased to larger species. Small mammals are represented exclusively in the data from the Nature Research Centre, which registered all road kills.
Out of all WVCs, we selected 21,911 WVC reports that involved mammals. A total of 19,622 reports included accurate information relating to 32 wild mammal species (Table 1). Of these reports, we mapped the 18,218 reports that included precise information on location ( Figure 1, Table 1).

Clustering of Collision Data
Using a clustering method, habitats were identified according to the location of hotspots. The literature contains many different spatial techniques for identifying short, significant road segments where collisions occur more frequently than usual [21,27,32,[61][62][63][64][65][66]. We utilized the KDE+, which analyses MVCs that are represented as point features and are located along the roads represented as line features (Figure 2a). The KDE+ algorithm finds locations (clusters) with statistically significant concentrations of collisions and assigns strength values (measured from 0 to 1) showing the risk severity to drivers [27,36] (Figure 2b). We performed MVC clustering analysis and created MVC clusters Using a clustering method, habitats were identified according to the location of hotspots. The literature contains many different spatial techniques for identifying short, significant road segments where collisions occur more frequently than usual [21,27,32,[61][62][63][64][65][66]. We utilized the KDE+, which analyses MVCs that are represented as point features and are located along the roads represented as line features (Figure 2a). The KDE+ algorithm finds locations (clusters) with statistically significant concentrations of collisions and assigns strength values (measured from 0 to 1) showing the risk severity to drivers [27,36] (Figure 2b). We performed MVC clustering analysis and created MVC clusters using the KDE+ parameters derived from the road network properties (KDE+ bandwidth-150 metres, Monte Carlo simulations-800, and minimal cluster strength-0.2).

Figure 2.
Roadkill-data-based identification, characterisation, and ranking of mammalian habitats: (a) MVC reports (small dots) with different species (dots marked as X and Y) placed within a road network (double-arrowed and labelled lines); (b) KDE+ clusters (short thick lines) labelled with underlined integer numbers show the length and non-integer numbers the strength of a cluster, small grey and white dots represent MVCs that did not form a cluster, and dashed lines represent the roads without clusters that did not form habitat patches; (c) Numbers show areas of habitat patches; (d) Numbers represent the length of hypothetical wildlife pathways (single-arrowed lines); (e) Numbers represent the length of hypothetical mammal corridors (double-arrowed lines); (f) Larger dots (habitat patches), darker colours of habitat patches, and thicker lines (corridors) represent a higher risk to drivers and higher attractiveness to wildlife, red lines (roads) highlight the highly ranked adjacent habitat patches, white-and black-coloured bars illustrate the share of species richness (for species X and Y), and labels show the number of mammals involved in MVCs (within the clusters); (b-d) Red dots represent the centre points of clusters; (c-f) Habitat patches (large green dots and polygons) labelled as ABCDE are represented by centre points.

Figure 2.
Roadkill-data-based identification, characterisation, and ranking of mammalian habitats: (a) MVC reports (small dots) with different species (dots marked as X and Y) placed within a road network (double-arrowed and labelled lines); (b) KDE+ clusters (short thick lines) labelled with underlined integer numbers show the length and non-integer numbers the strength of a cluster, small grey and white dots represent MVCs that did not form a cluster, and dashed lines represent the roads without clusters that did not form habitat patches; (c) Numbers show areas of habitat patches; (d) Numbers represent the length of hypothetical wildlife pathways (single-arrowed lines); (e) Numbers represent the length of hypothetical mammal corridors (double-arrowed lines); (f) Larger dots (habitat patches), darker colours of habitat patches, and thicker lines (corridors) represent a higher risk to drivers and higher attractiveness to wildlife, red lines (roads) highlight the highly ranked adjacent habitat patches, white-and black-coloured bars illustrate the share of species richness (for species X and Y), and labels show the number of mammals involved in MVCs (within the clusters); (b-d) Red dots represent the centre points of clusters; (c-f) Habitat patches (large green dots and polygons) labelled as ABCDE are represented by centre points.

Definition of Mammalian Habitats and Movement Patterns
Our conceptual model for the identification of mammalian habitats is shown in Figure 2a-c, characterization in Figure 2d,e, and ranking of habitat patches in Figure 2e. MVC reports with different species were mapped on the road network. Road sections where MVCs occurred more frequently were identified using the KDE+ clustering method [27].
We assumed that roadkill clusters are important indicators not only of risk to drivers [27], but also indicate locations where important mammalian pathways and roads intersect. We identified habitat patches as areas that are bounded (surrounded) by neighbouring road sections characterised by at least one cluster. We merged habitats having no clusters with neighbouring habitats iteratively until a merged habitat patch had a road with at least one neighbouring cluster. In our study, urban areas and urban clusters were excluded and not used for the identification of habitats. Identified habitats were used for their characterisation and, later, for ranking.
Hypothetical wildlife pathways were created by connecting the Clementini [67] centroids of habitat patches and cluster centroids using spider lines illustrating the shortest (Euclidean) distances. Hypothetical mammal corridors were created using the triangulated irregular network (TIN) between the Clementini [67] habitat patch centroids as peaks [37,42].

Characterisation of Mammalian Habitats
The habitat patches ( Figure 2c) were characterised using topological relationships between habitat patches, hypothetical pathways, and corridors. Each cluster centroid illustrates a "gateway" that mammals use to traverse from one habitat patch to another.
Following this conceptual framework (Figure 2), we identified and collected the necessary network-based criteria ( Table 2) for each habitat patch. Later, the habitat patches were ranked according to their attractiveness to wildlife and risk severity to drivers. Table 2. Criteria used for ranking the habitat patches described in Figure 2. They are considered as attractive to wildlife. ii Habitat patches with a higher number of shorter mammalian pathways and longer and stronger KDE+ clusters are characterised by higher numbers of collisions. They are considered as being a more severe risk to drivers. iii The number of species is an important indicator for both ( i , ii ) modelling assumptions, since a higher number of species within a certain habitat patch (species richness) simultaneously indicates a higher attractiveness to wildlife and a higher risk to drivers.

Objective Functions and Criteria Importance
The objective criterion importance (weights) for all criteria was calculated based on criteria utility (minimisation/maximisation) functions using SortViz for the ESRI inc. ArcGIS desktop software add-in [37,68].
Using the same ArcGIS desktop software add-in, we ranked habitat patches based on criteria derived from the individual ( Figure 2) and spatial connectivity properties (Table 2) of the habitat patch. In order to find habitat patches that were simultaneously the most attractive to wildlife and of most severe risk to drivers, modelling assumptions (see Table 2's footnote) and objective (utility) functions (Table 2) were set.

Ranking of Habitats and Ecological Corridors
Criterion importance values, defined as weights (Table 2), were then used as an input for ranking the habitat patches using the SAW and TOPSIS [49,52] methods. Both ranking approaches use the same input habitat data ( Table 2). The final SAW and TOPSIS values ranged from 0 (worst) to 1 (best) and altogether built the so-called 'composite indicator' of habitat attractiveness to wildlife (mammals) and risk to drivers. A higher rank value means higher attractiveness to wildlife and a higher risk to drivers. The SAW and TOPSIS values for each habitat patch were separately calculated and compared with each other ( Table 2).
Average SAW and average TOPSIS rank values, derived from the habitat patches connecting the two ends of the corridor, were allocated to each ecological corridor to determine the relative importance of the corridor (Figure 2f, Table 3).  Table 3) illustrate higher and more intense mammalian locomotion patterns [69] and risk to drivers.
We analysed land cover classes and species that had a strong (r > 0.50) relationship to habitat ranks (SAW and TOPSIS values). Habitat ranks were used as intercept and land cover classes and species as independent regressors.

Habitats and Habitat Characteristics
We identified 281 state-owned roads with at least one KDE+ cluster ( Figure 3): 18 main roads/highways; 107 national roads; and 156 regional roads (85.7%, 81.1%, and 9.6% of all roads in their respective category). The rest of the roads (thin grey lines in Figure 4) were not taken into account. tered in the road kills (Table 1) had no impact on the location and number of identified clusters.
We identified 3171 hypothetical mammalian pathways (thin grey lines in Figure 4), 672 corridors (dashed lines in Figure 4), and 243 habitat patches ( Figure 4). The hypothetical mammalian pathways (Figures 2d and 4) and corridors (Figures 2e and 4) were used for the characterisation and collection of criteria (Tables A1 and A2) (Tables 1 and 2).
Using the KDE+ method, we found 1642 mammalian clusters ( Figure 3), of which 22 (1.3%) were located in urban areas and therefore were excluded from further analyses. A total of 28 out of the 32 road-killed mammal species were identified within the clusters. Four small-sized mammals (M. glareolus, S. araneus, R. rattus, N. fodiens) were only registered as road kills outside the clusters. However, small numbers of these species registered in the road kills (Table 1)

Criteria Weights, Habitat Ranks, Ecological Corridors, and Movement Patterns
In order to rank the identified habitat patches (Table A1), the criteria weights (Table  4) were calculated using objective functions ( Table 2). The most important criterion for assessment was the shortest length of adjacent pathways, while the least important was the number of adjacent corridors. Table 4. Criteria (Table A1) and criteria weights used for ranking (following the same objective functions as in Table 2) the habitat patches ( Figure 3 Table 2, footnote identifier provided in the superscript. Figure 2, part identifier provided in the brackets.

Criteria Weights, Habitat Ranks, Ecological Corridors, and Movement Patterns
In order to rank the identified habitat patches (Table 1), the criteria weights (Table 4) were calculated using objective functions ( Table 2). The most important criterion for assessment was the shortest length of adjacent pathways, while the least important was the number of adjacent corridors. Table 4. Criteria (Table 1) and criteria weights used for ranking (following the same objective functions as in Table 2) the habitat patches ( Figure 3) in Lithuania.  Following objective functions, the SAW ( Figure 5) and TOPSIS ( Figure 6) ranks (Table 2) were assigned to each habitat patch. The average rank values were calculated for the corridors as well. The labels (Figures 5 and 6) identify the main roads.

Criteria Name * Weight (Index)
The highest SAW and TOPSIS rank values assigned to the habitat patches were 0.7 and 0.6, respectively. Furthermore, the SAW and TOPSIS ranks of habitats had a very strong correlation (r = 0.86), which means that the ranking results are similar and can be trusted.
The habitat patches contained from 3 to 477 MVCs and from 1 to 20 road-killed mammal species (Table A2). The corridor links (Figures 5 and 6) indicate the most probable movement patterns. The highly ranked corridors that intersect main roads highlight the highest potential risk to drivers and wildlife. Consequently, the MCV clusters that are on the roads with such intersections are of the highest importance for MVC mitigation actions.  The highest SAW and TOPSIS rank values assigned to the habitat patches were 0.7 and 0.6, respectively. Furthermore, the SAW and TOPSIS ranks of habitats had a very strong correlation (r = 0.86), which means that the ranking results are similar and can be trusted.
The habitat patches contained from 3 to 477 MVCs and from 1 to 20 road-killed mammal species (Table 2). The corridor links ( Figures 5 and 6) indicate the most probable movement patterns. The highly ranked corridors that intersect main roads highlight the highest potential risk to drivers and wildlife. Consequently, the MCV clusters that are on the roads with such intersections are of the highest importance for MVC mitigation actions.

Relationship between Habitat Ranks, Species Richness, and Land Cover Classes
Inside the clusters, MVCs with C. capreolus, S. scrofa, V. vulpes, L. europaeus, E. concolor, N. procyonoides, A. alces, M. putorius, and Martes sp. had strong relationships (r > ~0.5) with habitat patch ranks, showing the high severity risk to drivers and wildlife ( Figure 7). All other species had a weak or no relationship with habitat patch ranks. Five of these species, B. bonasus, L. lynx, M. erminea, L. lutra, and L. timidus, are rare in nature (Table 1), while others are small in size and their road kills were most probably under-registered.
Land cover classes such as road and rail networks, transitional woodland-shrub areas, mixed forest, broad-leaved forest, pastures, complex cultivation patterns, and discontinuous urban fabrics showed strong relationships (r > 0.5) with habitat patch ranks (Figure 8). All other land cover classes had a weak or no relationship with habitat patch ranks, indicating that these land cover classes do not pose a severe risk to drivers and wildlife.

Relationship between Habitat Ranks, Species Richness, and Land Cover Classes
Inside the clusters, MVCs with C. capreolus, S. scrofa, V. vulpes, L. europaeus, E. concolor, N. procyonoides, A. alces, M. putorius, and Martes sp. had strong relationships (r >~0.5) with habitat patch ranks, showing the high severity risk to drivers and wildlife ( Figure 7). All other species had a weak or no relationship with habitat patch ranks. Five of these species, B. bonasus, L. lynx, M. erminea, L. lutra, and L. timidus, are rare in nature (Table 1), while others are small in size and their road kills were most probably under-registered.  Land cover classes such as road and rail networks, transitional woodland-shrub areas, mixed forest, broad-leaved forest, pastures, complex cultivation patterns, and discontinuous urban fabrics showed strong relationships (r > 0.5) with habitat patch ranks ( Figure 8). All other land cover classes had a weak or no relationship with habitat patch ranks, indicating that these land cover classes do not pose a severe risk to drivers and wildlife. The results of multiple linear regression analyses (Table 5) indicate that broad-leaved forests and transitional woodland-shrub areas bordered by road and rail networks are characterised by the highest risk to drivers and wildlife. In the vicinity of such habitats, MVCs mostly occur with C. capreolus and S. scrofa. MVCs with other species such as A. alces, V. vulpes, Martes sp., M. putorius, L. europaeus, and E. concolor are also likely.  The results of multiple linear regression analyses (Table 5) indicate that broad-leaved forests and transitional woodland-shrub areas bordered by road and rail networks are characterised by the highest risk to drivers and wildlife. In the vicinity of such habitats, MVCs mostly occur with C. capreolus and S. scrofa. MVCs with other species such as A. alces, V. vulpes, Martes sp., M. putorius, L. europaeus, and E. concolor are also likely. 0.000803 ± 0.000 **** 0.0004290 ± 0.000 **** F (16,226) 98.02606 ± 0.000 **** 132.64396 ± 0.000 **** R 2 0.874 0.904 *-p < 0.10, **-p < 0.05, ***-p < 0.01, ****-p < 0.001. NS -not significant.
Using SAW ( Figure 9) and TOPSIS ( Figure 10) values, we created heat maps that show the potential risk severity to drivers and wildlife (urban areas were used as a reference). For better visual representation, the maps were created using the inverse distance weighed (IDW) interpolation method. The IDW method is used to interpolate spatial data and is based on the concept of distance weighting [74,75].   The SAW-based habitat patch heat map (Figure 9) shows more severe risk habitat patches than the TOPSIS-based heat map ( Figure 10) due to the differences in the ranking methods. However, both heat maps identified the same highly severe locations for drivers and wildlife.
Following the results from both ranking methods (Figures 9 and 10), we identified that the habitat patches with the unique identification numbers 577 and 2248 (Figure 3, Tables A1 and A2) posed the most severe risk to drivers and wildlife. For instance, around the top-ranked habitat patch (id: 577), which is bordered by the A14 main road and national roads 114, 111, 102, and 108, we found MVC clusters including 20 different mammal species (Figure 11). Most of the MVC clusters were found on A14. Clusters on the roads at the edge and within the habitat patch were also present. Due to the low traffic intensity there, we did not find any cluster on the national road 173, which is within the habitat patch. The SAW-based habitat patch heat map (Figure 9) shows more severe risk habitat patches than the TOPSIS-based heat map ( Figure 10) due to the differences in the ranking methods. However, both heat maps identified the same highly severe locations for drivers and wildlife.
Following the results from both ranking methods (Figures 9 and 10), we identified that the habitat patches with the unique identification numbers 577 and 2248 (Figure 3, Tables 1 and 2) posed the most severe risk to drivers and wildlife. For instance, around the top-ranked habitat patch (id: 577), which is bordered by the A14 main road and national roads 114, 111, 102, and 108, we found MVC clusters including 20 different mammal species ( Figure 11). Most of the MVC clusters were found on A14. Clusters on the roads at the edge and within the habitat patch were also present. Due to the low traffic intensity there, we did not find any cluster on the national road 173, which is within the habitat patch.

Habitat Risk Severity to Drivers
In order to plan MVC mitigation measures, spatial habitat characteristics together with MVCs and MVC cluster data are needed [76]. Habitat characteristics and factors that allow us to predict MVCs are important, but are usually the missing component. This can be explained by the disparate character of field research data [47]. Thus, the framework we propose may help to identify and characterise the missing components in a unified form.
Our results on habitat risk severity to drivers and wildlife at the local level are based on a long-term mammal roadkill dataset [77,78]. The main A14 road, delimiting the topranked habitat patch (id: 577, see Figure 11), is one of the most frequently checked for roadkill [1]. Because of ongoing long-term reconstruction of the A14 road (until 2030), short-term redirection of traffic onto national road 173 might be foreseen, thereby increasing the traffic intensity on that road, thus also increasing the likelihood of more MVCs than before and a higher risk to drivers.

Habitat Attractiveness to Mammals
The rates of annual land cover change in Lithuania are decreasing, dropping from 0.48% in 1990 to 0.18% in 2012 [58]. This indicates that the habitat composition has remained stable over time. A growing rate of forest land (woodland) and a rapid decline in active farming [58] has improved habitat attractiveness to different wildlife species, especially for forest dwellers. The increasing MVC numbers in all categories of roads and the

Habitat Risk Severity to Drivers
In order to plan MVC mitigation measures, spatial habitat characteristics together with MVCs and MVC cluster data are needed [76]. Habitat characteristics and factors that allow us to predict MVCs are important, but are usually the missing component. This can be explained by the disparate character of field research data [47]. Thus, the framework we propose may help to identify and characterise the missing components in a unified form.
Our results on habitat risk severity to drivers and wildlife at the local level are based on a long-term mammal roadkill dataset [77,78]. The main A14 road, delimiting the topranked habitat patch (id: 577, see Figure 11), is one of the most frequently checked for roadkill [1]. Because of ongoing long-term reconstruction of the A14 road (until 2030), short-term redirection of traffic onto national road 173 might be foreseen, thereby increasing the traffic intensity on that road, thus also increasing the likelihood of more MVCs than before and a higher risk to drivers.

Habitat Attractiveness to Mammals
The rates of annual land cover change in Lithuania are decreasing, dropping from 0.48% in 1990 to 0.18% in 2012 [58]. This indicates that the habitat composition has remained stable over time. A growing rate of forest land (woodland) and a rapid decline in active farming [58] has improved habitat attractiveness to different wildlife species, especially for forest dwellers. The increasing MVC numbers in all categories of roads and the increase in annual average daily traffic [79] have coincided with an enlargement of wildlife populations in the country. Species richness (the number of species in Table 2) has a strong relation-ship (r = 0.72) with the number of MVCs (the number of MVCs within adjacent clusters in Table 1).
We assumed that larger values of species richness indicate higher habitat diversity [80], suitability [81], and attractiveness to wildlife. We found 20 different species within MVC clusters that are adjacent to habitat patch id: 577, which means that road 137 ( Figure 11) is more dangerous than the roads adjacent to habitats with a smaller number of species. However, species richness does not take into account the abundances of the species or their relative abundance distributions. The proposed framework allows for the accurate identification of species richness in relation to MVCs that are in the vicinity of the particular habitat. This information is especially useful when wildlife observation data (ground-truth) are not available at all. However, the accuracy of the result is very much dependent on the quality of the available police registered reports and professional field research data [1].
Habitats were defined and characterised across all territory of Lithuania; therefore, the validation of our model is possible: (i) using data from a similar territory, such as a neighbouring country; or (ii) using data from Lithuania from a different time period, e.g., 2018-2021 (our model covered 2002-2017). At the moment, however, such a dataset is not available.
Species richness may be validated by intensive roadkill counts or using wildlife cameras to check for animal movement across roads.

Multi-Objective Mitigation Measures
The only effective mitigation of road kills in a multi-species animal community is a complex of wildlife fencing (with a sufficient number of wildlife underpasses and overpasses according to the length of the fence) and active driver warning systems on roads without wildlife fences. We did not manage to find tangible research on other effective multi-species and multi-objective mitigation measures for large and small mammal species.
Mitigating MVCs on road 173 ( Figure 11) may be challenging, as the MVC-targeting measures are likely to focus de minimis on ungulates, namely C. capreolus, S. scrofa, and A. alces (Tables 2 and 5), rather than on the other 17 large and small body size mammal species recorded ( Table 2). Numbers of carnivore road kills also grow in areas with a higher abundance of small mammal species [55]. MVC clusters found in different locations can help us to select species-specific mitigation measures. However, due to the high cost of the abovementioned complex of measures and the low traffic intensity on roads other than A14 and 102 ( Figure 11), implementation of such measures is not possible in the near future. Therefore, our method currently may serve as part of the toolbox to identify the most dangerous roads and the most important habitat patches.
The observation of near misses (road 173 in Figure 11) might provide further input for the task. Field studies should incorporate long-term data collection, before the mitigation measure is applied [18]. Last, but not least, clearing vegetation along roads can also help to lower the MVC risk [54,82]. The mitigation measures for managing the risks to drivers and wildlife may be more challenging when many species are present. This may result in higher road reconstruction costs. The lack of data on the effectiveness of road mitigation measures [18,20] is a further obstacle to decision-making. The most common MVC mitigation measure in Lithuania is fencing. Short wildlife fences may not sufficiently reduce the risk of MVCs, but they are economically more affordable. Long fences are less efficient economically, but may perform better [9][10][11]17] on the roads with the highest traffic intensity. Therefore, we conclude that, even when involving all habitat data, the selection of multi-objective MVC mitigation measures in a dynamic environment still remains a considerable research challenge.

Conclusions
This study developed models that allow for the identification, characterisation, and ranking of habitats based on mammal roadkill data. The main conclusions are:

1.
Habitats were characterised by connectivity, land cover, roadkill, roadkill cluster, and mammal species and ranked using multiple criteria for the identification of habitat risk severity to drivers and attractiveness to wildlife; 2.
Despite the potential limitations of the scope of the roadkill data, our habitat ranking suggests that this procedure can provide information on habitats, habitat locations, species richness, habitat risk severity to drivers, and attractiveness to wildlife; 3.
Strong relationships were identified and discussed between the habitat patch ranks, five (out of 28) land cover classes, and eight (out of 28) species (97% of all mammal road kills); 4.
This methodology facilitates decision-making on the habitats that must be prioritized to preserve wildlife in the vicinity of roads that are prone to MVCs. It is also suitable for the planning of multi-objective mitigation measures to improve road security in a dynamic environment.

Acknowledgments:
The authors thank Jos Stratford for linguistic revision of the manuscript as well as the anonymous reviewers for their suggestions and contributions.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A