Characterisation of Beaver Habitat Parameters That Promote the Use of Culverts as Dam Construction Sites: Can We Limit the Damage to Forest Roads?

: The use of forest roads as foundations for dam construction by beavers is a recurrent problem in the management of forest road networks. In order to limit the damage to forest roads, our goal was to calculate the probability of beaver dam installation on culverts, according to surrounding habitat parameters, which could allow for improvement in the spatial design of new roads that minimise conﬂicts with beavers. Comparisons of culverts with ( n = 77) and without ( n = 51) dams in northwestern Quebec showed that catchment surface, cumulate length of all local streams within a 2-km radius, and road embankment height had a negative effect on the probability of dam construction on culverts, while ﬂow level and culvert diameter ratio had a positive effect. Nevertheless, predicted probabilities of dam construction on culverts generally exceeded 50%, even on sites that were less favourable to beavers. We suggest that it would be more reasonable to take their probable subsequent presence into account at the earliest steps of road conception. Installing mitigation measures such as pre-dams during road construction would probably reduce the occurrence of conﬂicts with beavers and thus reduce the maintenance costs of forest roads.


Introduction
In Canada, transportation of wood between forests and mills is mainly overland, which requires the construction and maintenance of an important road network.To prevent road damage, several forest industries and associations throughout North America need to decrease the densities of beavers, which represent a significant issue wherever they are naturally established [1][2][3][4].Indeed, beavers frequently use roads as the foundations for dam construction, which is a recurrent problem in the management of forest road networks.Dams are usually located on the upstream side of roads, usually a few metres away or even within culverts, i.e., structures that are supposed to allow the free movement of water between the two sides of roads.Road embankments also facilitate dam construction, because they prevent water flow and may seem attractive to beavers.Embankments have the effect of raising the water table, while lowering the water table downslope, and concentrating the flow at the culvert level [5].When the culvert is blocked and water accumulates upstream, significant infiltration will occur, leading to destabilisation of roads, frequent flooding and even complete road erosion at the time of severe flooding [6].This damage obviously represents significant road network maintenance costs, often at the expense of forest companies, in several Canadian provinces.
Following a major population decline that was related to trapping from the beginning of European settlement of northeastern regions of North America, several beaver populations have apparently increased at a time when the forest road network in the boreal zone is still rapidly growing.For example, according to the report of [7] and data obtained from the Ministry of Natural Resources and Wildlife of Quebec [8], the network length of permanent forest roads throughout Quebec has tripled over the last 30 years, while timber harvest increased only about 60%.Indeed, the road construction rate increased from 150 km/year in the late 1970s to 300 km/year in the late 1990s.As a result, an estimated 10,000 culverts are installed on new roads or repaired within the existing road network every year throughout the province of Quebec [7].Within this context, potential conflicts between beavers and users of forest roads are likely to increase in the coming years.
In order to limit the damage that is caused to roads by beavers, the effectiveness of different control systems has been the subject of several recent studies [2,4,9].Control measures have included systems that are aimed at preventing culvert blockage by beavers or for maintaining water level control.These systems require rigorous maintenance and other measures, such as trapping.These approaches are effective occasionally and, therefore, short-lasting [6].Thus, long-term management approaches that limit the risk of damage to forest roads by beavers are greatly needed.
One of the first steps toward longer-term management to limit risks to road infrastructure is to increase our understanding and approaches to modelling beaver habitat.Generally, geomorphological variables are more useful than variables reflecting food availability in explaining habitat use by beavers [10,11].We consider the following variables to be appropriate for dam-site selection: stream gradient [12,13], dam upstream catchment area [13,14], Strahler order [15], channel geomorphic unit (stream width, depth or cross-sectional area [12,16], floodplain size [13,17], riverbank slopes [18], and substrate type [19].Yet some studies also suggest that the presence, abundance, size and distribution of deciduous stems that are most frequently used by beavers, either as food or as building materials, are also crucial variables that should be considered within habitat models [20,21].In order of preference, the most important species in the beaver's diet are trembling aspen (Populus tremuloides Michx.), willow (Salix spp.), balsam poplar (Populus balsamifera L.) and alder (Alnus spp.) [22].According to [23], some woody plants, such as alder, are harvested by beavers mainly for the construction of dams rather than being utilised as a food source.The use of woody stems for shelter or food appears to depend upon the diameter of the stems rather than the species of plant [21].
Some studies have considered the effect of road or culvert presence on the distribution of dams within the landscape and most have found a positive association [13,17,24,25].Two of these studies were conducted in boreal Alberta.In the first study, dams were observed on all inventoried roads, while in the second, the attractive effect of culverts was verified, though only at a scale of 300 m [24,25].Furthermore, research in New York State was undertaken specifically on habitat use by beavers near culverts; however, these studies dealt with highways rather than forest roads [6,26].Nevertheless, the research has shown that oversized culverts allowing greater water flow were avoided by beavers [6].Moreover, it demonstrated that the probability of site occupation by beavers was inversely related to the proportion of the roadside area that was devoid of woody vegetation and an increasing (>3%) stream gradient [26].
In Quebec, the Abitibi-Témiscamingue region exhibits the highest beaver densities, with around 5.5 colonies/10 km 2 [27], together with the highest dam densities, with 4.01 to 6.22 dams/10 km 2 [10].Also, forest roads are widely established across the territory, harvesting activities are widely distributed due to their economic significance, and road density is constantly increasing [7], which has led to a growing number of conflicts with beavers.A better understanding of factors affecting the use of culverts by beavers on forest roads should allow managers to improve the localisation of new roads as well as culverts and, if necessary, target the most profitable locations for installing mitigation control systems.The objectives of this study were thus to identify the contribution of geomorphological, hydrological, anthropic and food factors, together with surrounding dam densities, to the risk of dam construction by beavers on new culverts.To do this, we compared with a model selection approach the relative weight that was assigned to each of the hypotheses that were identified a priori to potentially explain the presence of dams, over 128 culverts where the presence or absence of beaver dams had been confirmed on site.Our main prediction was that geomorphological factors would predominate explanations regarding the presence of beaver dams.

Study Area
The study area covers 8000 km 2 and is located in northwestern Quebec, within the heart of the boreal forest.Study sites were about 50 km away from the city of Rouyn-Noranda (79.06 • W; 48.12 • N).Located in the Abitibi region and James Bay lowlands, topographic relief of the study territory is relatively uniform and flat (235 m to 542 m above sea level) and mostly covered by glaciolacustrine clays left by the withdrawal of proglacial Lake Barlow-Ojibway [28].Among other features, the region is characterised by a well-developed river system [27].The extent of the study area is limited on its western edge by the Province of Ontario, by the Clay Belt at the southern and eastern boundaries, and by agricultural areas around Lake Abitibi at the northern side (Figure 1).In addition, the sites are mostly located in the western balsam fir-white birch bioclimatic subdomain, as described by [29].In boreal mixedwoods, black spruce (Picea mariana [Mill.]BSP), jack pine (Pinus banksiana Lamb.), white or paper birch (Betula papyrifera Marsh.), and trembling aspen generally dominate the landscape.Balsam fir (Abies balsamea [L.] Mill.), white spruce (Picea glauca (Moench) Voss) and eastern white cedar (Thuja occidentalis L.) are also present in the territory.Throughout the study area, there are registered traplines.Yet beaver trapping is generally minimal due to the low market value of the furs.
As a first step, all practicable forest roads in the study area were visited to select sites.Each location corresponds to a culvert, which was characterised by the presence or absence of a beaver dam.When blocked culverts cause damage to roads, managers usually undertake dam removal and then engage professional trappers to eliminate the associated colony.This combined action often prevents fast rebuilding of the dam.However, such dam removals leave evidence (drained ponds, cluster swathes and eaten vegetation) in the area and we considered culverts showing such signs of beaver occupancy.Similarly, culverts exhibiting control of beaver flooding were also considered among the sites that beavers had occupied.Of 128 sampled culverts, 77 were used as dam building sites, while 51 were not.

Geomorphological and Hydrological Variables
Hydrographic data include catchment area (SBD), the ratio between flow level and culvert diameter (Ddia), catchment area gradient (inclBD), stream gradient (inclR), Strahler stream order (Order), cumulative length of all local (2-km radius) streams (kmR2km), and an index of riparian slopes (CoteP).All variables except CoteP and culvert diameter (field data) were extracted from a digital elevation model (DEM) available from the topographic database of Quebec (BDTQ), with 10-m resolution.Using the BASIN (version EX_02, ESRI, Redlands, CA, USA) (BASIN script allows catchment delineation from a point on a stream and calculation of the catchment's characteristics, together with generation of a stream network, display of a raindrop path traced from a point and display of elevations extracted from a grid (cell value) or from a surface, for a point defined by the user.)script in ARCGIS (http://desktop.arcgis.com/fr/arcmap/10.3/tools/spatial-analysttoolbox/basin.htm), a drainage pattern for the study area was produced.Drainage basins are created by locating the pour points at the edges of the analysis window (where water would pour out of the raster), together with sinks, then identifying the contributing area above each pour point.When creating the linear layer representing the streams, the script automatically assigns an order to each segment using Strahler's method [30].This information was retrieved directly for stream segments that contained a sampled culvert.Culverts that were located more than 150 m from a stream, depending upon the related upstream alignment, were not considered to be associated with a mapped stream and were thus assigned a value of zero in Strahler's order.By selecting a point on the map, BASIN also allows the estimation of areal surface that is being drained to that point.The drainage basins are identified within the analysis window by identifying ridgelines between basins.The theoretical surface of the catchment area was calculated for each culvert using such methods.The flow rate was calculated using the GSFdebit tool (version 10.3.002,Groupe Système Forêt, Québec, QC, Canada) for ArcGIS, which was developed by Groupe Système Forêt, Lévis, Québec (http://www.gsf.ca/en-ca/applications/gsf-debit-pour-arcgis.aspx).The tool also calculates the average gradient of the catchment area, together with the main stream gradient.Since we knew that flow rate variation would be small due to regional topography and the streams that were sampled, we used the ratio between the estimated flow level and culvert diameter as an indicator of the amount of water available versus bottlenecks to be blocked.
Forests 2017, 8, 494 4 of 13 assigned a value of zero in Strahler's order.By selecting a point on the map, BASIN also allows the estimation of areal surface that is being drained to that point.The drainage basins are identified within the analysis window by identifying ridgelines between basins.The theoretical surface of the catchment area was calculated for each culvert using such methods.The flow rate was calculated using the GSFdebit tool (version 10.3.002,Groupe Système Forêt, Québec, Canada) for ArcGIS, which was developed by Groupe Système Forêt, Lévis, Québec (http://www.gsf.ca/en-ca/applications/gsfdebit-pour-arcgis.aspx).The tool also calculates the average gradient of the catchment area, together with the main stream gradient.Since we knew that flow rate variation would be small due to regional topography and the streams that were sampled, we used the ratio between the estimated flow level and culvert diameter as an indicator of the amount of water available versus bottlenecks to be blocked.

Figure 1.
Locations in Abitibi, Québec, of 128 culverts that were sampled for the presence or absence of beaver dams.
Disparities were observed when overlap was calculated between readily available "base de données topographiques du Québec" (BDTQ) water systems and that created with BASIN.Thus, several catchment limit differences were observed, together with those for other resulting variables: flow level, catchment area gradient and stream gradient.To take into account the reality that we felt best represented by the BDTQ water system (derived through photo interpretation), we retained coarse variables from the DEM and compared them to those that were obtained using water systems that were available in the BDTQ database to validate a systematic method of extraction.The aim was to develop a practical tool that could be used by land managers, considering the extremely useful potential applications of current software developments.In order to see whether the density of Disparities were observed when overlap was calculated between readily available "base de données topographiques du Québec" (BDTQ) water systems and that created with BASIN.Thus, several catchment limit differences were observed, together with those for other resulting variables: flow level, catchment area gradient and stream gradient.To take into account the reality that we felt best represented by the BDTQ water system (derived through photo interpretation), we retained coarse variables from the DEM and compared them to those that were obtained using water systems that were available in the BDTQ database to validate a systematic method of extraction.
The aim was to develop a practical tool that could be used by land managers, considering the extremely useful potential applications of current software developments.In order to see whether the density of available streams had an effect on the density of beavers that were likely to colonise the culverts, the cumulative length of all local (2-km radius) streams of each sampled culvert was also calculated using BASIN script.Lastly, the CoteP riparian slope index was measured using a clinometer for both riverbanks that were upstream of roads.Gradient degrees were divided into four classes: Class 0 (null slope), Class 1 (0.1 to 3%), Class 2 (3.1 to 6%) and Class 3 (>6%).Classes of each riverbank were then summed to provide a rating of riparian slope for each culvert.We expected that higher values would be associated with lower dam-site selection by beavers, because of higher risks of floods.

Variables Associated with the Availability of Food and Dam Building Materials
Dietary variables and building material availability were estimated according to hardwood dominance (Hardw) and alder (Alnus sp.) presence (Al).This information was derived from current provincial forest inventory ecoforestry maps.Areas (ha) that were covered by deciduous or mixed type stands-dominated by trembling aspen or white birch in our study area-were added to a sampling area measuring 200 m in radius around each culvert.Sampling area radius was established based upon consultations with trappers regarding their experience with beaver foraging behaviour in our study area.The presence of alders over at least 1 ha was subsequently recorded and compiled as a categorical variable (presence or absence) for each culvert.

Anthropogenic Variables and Characteristics of Culverts
These data refer to the height of the road (embankments) and road age category.Road embankment height was measured directly on sites.Road embankment height was measured upstream, from contact with the water stream to the highest point that was perpendicular to the road bed.To determine the construction year of roads, we used the harvest year and sylvicultural activities close to the roads that were registered in annual forest operation reports.A spatial chronosequence (1985, 1995, 2000 and 2005) from classified Landsat images [31] was also used to validate and complete the information for more recent roads.Given the narrow range that was obtained for the age distribution of roads, with recent roads much less represented in our sample, we decided to use two age groups in subsequent analyses: group 0 for recent roads (construction year ≥1985) and group 1 for older roads (<1985).

Local Abundance of Beaver Dams
Beaver density is generally assessed by counting the number of lodges, which was information not readily available in our study area.However, precise locations of beaver dams (at least 10 m in length) were available on the BDTQ map of streams.Although we are aware that dam abundance data might not provide a good surrogate for beaver colonies, we hypothesised that the local abundance of existing beaver dams might influence the likelihood of using a culvert as a new dam site.The local abundance of beaver dams was thus estimated according to the number of surrounding dams that were within two kilometres of each culvert (NbB).

Statistical Analyses
Logistic regression was used to predict the probability of dams near culverts.Each model represents a hypothesis that could explain the presence of beavers near culverts, and which contains one or more variables based upon biological knowledge of beaver habitat.Models refer to geomorphology and hydrology, food and dam-building materials, human factors, and surrounding beaver populations (Table 1).Two models were compared in the case of geomorphology and hydrology: a first model was obtained from DEM data (HYDRO1) and a second one was obtained from BDTQ (HYDRO2), which also involves two global models (Global1 and Global2), totalling eight models.Lastly, to take into account the effect of time since road construction, the time that had elapsed between road construction and sampling was added to each model as an offset.
All models subsequently were compared with one another, using the Akaike information criterion (AICc) [32].Weighted estimates and unconditional standard errors were calculated for each variable in the best models.The effects of variables were determined according to a 95% confidence interval.
To illustrate the magnitude of the effect of each variable, multi-model inference was determined for the predicted probabilities from each model, for the range of values that were covered by each variable of interest, by keeping the other variables fixed.All statistical analyses were performed in R (v. 2.7.2; http://probability.ca/cran/).

Results
Regardless of whether dams were incorporated into culverts or not, the range of estimated values for each variable (HYDRO1, HYDRO2) was similar in both cases (Table 2).Further, when comparing the average value with the range of values for each variable, few sampled culverts had values close to the upper limit (Table 2).Indeed, most sampled culverts exhibited values closer to the lower limit.Following the comparison of eight candidate models using AICc, two plausible models emerged from among those that were set a priori to explain the presence of dams near culverts.The first model contains only geomorphological and hydrological BDTQ variables (HYDRO2), with a delta AICc weight of 0.46 (Table 3).The second model includes the height of the embankment (ANTHROPO), with an AICc delta of 0.34 and an AICc weight of 0.38 (Table 3).According to evidence ratio statistics, the HYDRO2 model is approximately equivalent to the ANTHROPO model and roughly seven times better than the next model, i.e., the global model with BDTQ data (GLOBAL2).All other models can be dismissed as being less plausible, since AICc deltas are close to or above 4, and AICc weights are near or below 0.05 (Burnham and Anderson, 2002).Multi-model inference allowed us to identify catchment area surfaces (SBDbdtq), cumulative length of all local (2-km radius) streams (kmR2km) around study culverts and road embankment heights (Embank) as all having a negative effect on the probability of encountering beaver dams (Table 4; Figure 2).BDTQ level flow and culvert diameter ratio (Ddiacor) had, in turn, a positive effect (Table 4; Figure 2).All other variables in the best two models had no effect, since they all included 0 in their respective confidence intervals (Table 4).Finally, the adjustment test of Le Cessie and Van Houwelingen [33] over the HYDRO2 model with BDTQ data confirms a good fit (p = 0.987).Moreover, according to Cox and Snell's R test [34] for the same model, about 18.7% of the probability estimates for culvert dams may be explained by the study variables (R 2 = 0.187).

Variables Predictive Power and Support of Submitted Hypotheses
As suggested by reports in the literature [12,13,15,21] and our hypotheses, the effects of geomorphological and hydrological parameters were the most significant in predicting the presence of beaver dams near culverts.In addition, we observed that coarse variables from the DEM

Variables Predictive Power and Support of Submitted Hypotheses
As suggested by reports in the literature [12,13,15,21] and our hypotheses, the effects of geomorphological and hydrological parameters were the most significant in predicting the presence of beaver dams near culverts.In addition, we observed that coarse variables from the DEM (HYDRO1) were not sufficient to make valid predictions.Indeed, the model (HYDRO2) with the BDTQ water system data was the most suitable, whereas models with the DEM alone were the least plausible.Similarly, the model that was based upon anthropogenic factors also has explanatory power, yet assumptions on dietary variables and the local abundances of beaver dams were not supported by the model selection process.

Effects of Geomorphic and Hydrologic Variables
Some geomorphology and hydrology variables that were included in the model (HYDRO2) appear to exert a greater effect than others when predicting the probability of dams in culverts.This is particularly true regarding the values of catchment area surface (SBDbdtq), level flow and culvert diameter ratio (Ddiabdtq), together with the cumulative length of all local (2-km radius) streams (kmR2km).
A high proportion of beaver dams was located within medium-sized catchment areas between 500 and 5000 ha, as was already observed in a previous study that was conducted in the upper Coastal Plain of South Carolina [13].At the lower end of this catchment area gradient, the authors concluded that small watersheds are inadequate for the establishment of beavers due to the intermittent presence of surface water.On the higher end of this gradient, excessive water velocity and other physical problems that are related to the amplitude of discharge reduced the suitability of the habitat for beavers in large river basins.Since the observed values in this study ranged from 0 to 3000 hectares, we expected that the relationship between the probability of occurrence and the catchment area surface would be positive.The negative relationship that we obtained could be due to the fact that beavers may create ponds, even on culverts that are located in a low catchment area, provided that there is a minimal water flow, which likely exists given the boreal context of our study area.
Since level flow variation range is relatively small in relation to the topography and study streams, we evaluated the relationship between level flow and culvert diameter.As expected, this ratio results in a positive relationship with the likelihood of beaver dam construction.On one hand, a strong ratio suggests a significant amount of water (high flow level) or a small culvert diameter; on the other hand, these characteristics were both likely to make sites attractive to beavers.Indeed, high flow levels ensure considerable water input throughout the year, and small culverts may limit the efforts to build dams.Finally, we had hypothesised that the number of streams near culverts could have increased the chances of settlement of one or more colonies and that of colonisation of culverts.However, the negative effect of this variable suggests otherwise: where streams are naturally rare, fewer potential spots for dam settlement exist, which thus increases the likelihood of dam construction on culverts.

Effects of Anthropogenic Variables
A positive selection effect has previously been shown for dam sites that allow the creation of larger ponds [35].Culvert sites with high road embankments would thus be more likely to be colonised by beavers since they allow greater water blockage.Indeed, there is generally a strong correlation between embankment height and culvert diameter.In this sense, we suggest that the negative relationship that was obtained with road embankments may not be linked to the avoidance of sites allowing the creation of larger ponds, but rather to sites with greater diameter culverts [6].

Effects of Dietary Related Variables
Like the data that were obtained by [10], our results suggest that the availability of hardwood species neither presents a factor limiting food in the study area, nor do they represent a constraint on the availability of construction materials.It is even likely that beavers sometimes use coniferous stems for the construction of dams or huts [25].We also have no information regarding the purpose of the ponds that were created by roadside dams.Indeed, we know that beavers build a primary pond (one containing a hut) and, sometimes, secondary ponds that provide access to more food or which maintain more stable water levels upstream [36].In our case, we noted that most ponds that were observed in the field did not contain a hut, while other habitat studies have generally related to colonies and, therefore, primary ponds.It is thus possible that some variables that proved significant in these other studies do not provide information when it comes to secondary ponds, which could be the case for hardwoods that are used primarily as food.It is also suggested in the literature that dietary parameters are difficult to consider in the case of beaver habitat studies, since beavers greatly shape their habitat themselves [15].Consequently, riparian vegetation can be greatly modified.

Study Limitations
In this study, adjusted estimates were very low (<1; Table 3), even for all four of the most important variables, which suggests weak relationships.We believe that within the range of observed values for each variable, there is no threshold at which it would be more difficult for beavers to create dams, which reduces the discrimination capacity of the study.This can be explained by the fact that sampling sites are limited to culverts, that is to say, small-scale streams, as compared to sites where priority is given to the construction of short bridges.
It should also be mentioned that all geomorphological and hydrological data for this study came from a digital elevation model that was designed for level curves, with ten-metre intervals, from Quebec's topographic database (BDTQ).Since the topography is flat, it may be necessary to travel long horizontal distances to attain a ten-metre rise in elevation.In such cases, basins or topographical variations may exist, but they are much more subtle than level curve detail and, thus, not adequately represented.This could explain the presence of catchment areas with an area close to zero, together with the resulting variables (flow, catchment area and stream slope).Given the progress and opportunities created by new technology, the promising use of LiDAR-based remote sensing, for example, could contribute to improving our models [37].
Furthermore, classification errors may have occurred in terms of data regarding the presence of dams in culverts.Indeed, there is no record of culvert maintenance.Certain site indices allowed us to know whether culverts had been recently used by beavers or not, but false absences may occur in some cases.In addition, culverts are sometimes replaced.Site modifications and machinery may lead us to believe that dams were absent, whereas culverts are often replaced precisely to eliminate the harmful presence of a beaver dam.
Finally, we may wonder whether a road itself might represent a very influential factor in explaining the presence of dams in the culverts.If so, the presence of the road alone in many cases may account for a larger portion of the phenomenon better than all other variables together, thereby justifying in part the low predictive power of the models (Cox and Snell's R equivalent to 18.7%).

Management Implications
One of our main conclusions is that the predictive power of our models proved to be insufficient to make them reliable for implementation in the context of road network management.First, regarding the extent of the observed values for the study variables, the probability of dam presence is always relatively high (>50%, see Figure 2), suggesting that virtually all territory exhibits a sufficiently high potential for dam construction and the risk is significant.Also, for all four variables that accounted for the most important effects, the confidence interval on the predicted probability of dam presence varies from nearly 50% to 100% over the whole range of observed values (Figure 2).In that sense, it would be appropriate to take into account the potential presence of a dam when planning roads, rather than seeking places that would allow for minimisation of the problem.
We suggest that the incorporation of mitigation measures, such as pre-dams, should be planned at the time of construction of roads.Pre-dams are short semicircular, circular, or rectangular fences that divert the beavers' attention away from a culvert to the fence [38].These devices are also called diversion dams because, unlike beaver deceivers which function by preventing damming, pre-dams are intended to be dammed.We believe that in many cases, this would reduce both the cost of culvert maintenance, and that related to the installation of the pre-dam itself.Indeed, the installation costs are often lower when performed at the time of construction, since the machinery is already present on sites.

Conclusions
In short, among the factors that were considered in this study, i.e., geomorphology and hydrology, food availability, human factors and the density of surrounding dams, the study demonstrated that geomorphology and hydrology variables are the most important in explaining the presence of dams in culverts.However, the predictive power of our models was insufficient and unreliable in the context of road network management.These variables were considered, since the catchment area upstream of the road and the cumulate length of all local (2-km radius) streams have a negative effect on the probability of the presence of a dam, whereas the ratio between flow and culvert diameter has a positive effect.In addition, the height of the embankment, which most often refers to the diameter of the culvert, also has a negative effect on the probability of dam presence.
In most cases, the predicted probability of dam construction in culverts was equal to or greater than 50%, which obviously means that the entire study area is highly susceptible to beaver establishment.The promise offered by LiDAR remote sensing technology could contribute to improving our models.Until then, it would be preferable to install mitigation measures such as pre-dams during road construction.This would likely reduce the occurrence of conflicts with beavers and, thus, reduce the maintenance costs of forest roads.In addition, we know that the pre-dam system is less expensive when done at the time of construction of roads compared to a breach in the road caused by flooding (washout), because the machinery required for the operation is already in place (St-Amant, pers.comm.).
Research Perspective -Conducting the same kind of research in a hilly region could establish thresholds and confirm whether study variables are truly potential explanatory variables.-Using more specific topographic data (though more expensive), such as those acquired with LIDAR, to determine whether a more detailed and more realistic geomorphological and hydrological characterisation of the landscape would provide a better explanation of observed phenomena.More accurate data would also allow us to calculate the potential catchment flooded area according to the height of embankments.Indeed, we know that flood area has a positive effect on the presence of beaver dams, as has been determined in other habitat studies for this species.-Assessing the contribution of the presence of the road as such to explain the presence of dams at culverts.Indeed, it is believed that in some cases, beavers preferentially set their dams at the side of roads, although other aspects of habitat would be less favourable.

Figure 1 .
Figure 1.Locations in Abitibi, Québec, of 128 culverts that were sampled for the presence or absence of beaver dams.

ForestsFigure 2 .
Figure 2. Probability of dam presence near culverts according to (a) BDTQ catchment area upstream of the road, (b) the number of km of streams within 2-km radius around the culvert, (c) embankment height and (d) BDTQ flow rate and culvert diameter ratio.The dotted lines represent 95% confidence bands and dots represent sample culverts (n = 128).

Figure 2 .
Figure 2. Probability of dam presence near culverts according to (a) BDTQ catchment area upstream of the road, (b) the number of km of streams within 2-km radius around the culvert, (c) embankment height and (d) BDTQ flow rate and culvert diameter ratio.The dotted lines represent 95% confidence bands and dots represent sample culverts (n = 128).

Table 1 .
Variables included in candidate models to test different hypotheses explaining presence of dams.

Table 2 .
Descriptive data for 12 study variables that were used to explain the presence or absence of dams.
* See Table 1 for description of variables.** The median is given for CoteP and AGE, and proportion of AI variable.

Table 3 .
Results of model selection using the Akaike information criterion (AIC).

Table 4 .
Results of multi-model inference for the variables included in the top two models (HYDRO2, ANTROPHO).
* See Table1for description of variables.