Eskers are complex geological formations shaped with a linear accumulation of sand and gravel under the glaciers during the last ice age and that provide crucial resources such as drinking water, sand/gravel, outdoor recreational sites, and productive forests. Surrounding sand and gravel and connection with the groundwater influences the physicochemical properties of lakes on esker which can benefit different biotic communities in the food web. The sustainable management of resources provided by eskers requires baseline ecological knowledge of these ecosystems. However, very little information exists about the ecology of freshwater ecosystems on eskers. This study uses a food web approach to identify the environmental variables, biological diversity, and indicator species associated with esker lakes to better understand their ecological functioning and biodiversity patterns to benefit their sustainable management and conservation. Fifty lakes were sampled in the Abitibi-Témiscamingue region (Canada), half on eskers and half on the surrounding boreal clay belt to include the most abundant lake ecosystems of the region. Physicochemical, environmental, and anthropogenic variables measured in the two lake types showed that esker lakes differed markedly from clay lakes. Nutrient concentrations, conductivity, and macrophyte cover were significantly lower in esker lakes than in clay lakes, whereas dissolved oxygen saturation and concentration showed the opposite trend. Three interconnected trophic levels of the esker lake food webs—waterbird, fish, and macroinvertebrate communities—were characterized for biological diversity and the associated species. We found a significantly lower Shannon diversity index for waterbirds (mean ± standard deviation; 0.7 ± 0.2), fish (0.4 ± 0.3), and a tendency for a lower value for macroinvertebrates (0.9 ± 0.3) in esker lakes than the clay lakes (1.1 ± 0.4, 0.9 ± 0.3, and 1.3 ± 0.5, respectively). Common goldeneye (Bucephala clangula) and Canada goose (Bucephala clangula) were associated significantly with esker lakes and identified as indicator species for esker lakes. In contrast, ring-necked duck (Aythya collaris) and hooded merganser (Lophodytes cucullatus) were associated significantly with clay lakes. Perlidae was similarly associated with esker lakes as an indicator for macroinvertebrates. Anthropogenic activities such as forest harvesting have altered the waterbird community, and recreational activities around the lakes have modified the fish and macroinvertebrate communities. We conclude that esker lakes differ from other regional lakes and are associated with specific environmental and biological variables and indicator species. The biological diversity in esker lakes is lower than that of clay lakes for all studied trophic levels of the food web, but these waterbodies provide preferential habitats for some species. This research provides the first baseline ecological information necessary to establish sustainable management and conservation strategies for this vulnerable ecosystem.

Wetlands and predation in boreal ecosystems play essential roles throughout the breeding season for bird assemblages. We found a positive association of beaver activity and a negative influence of American red squirrels (Tamiasciurus hudsonicus) on bird assemblages. We used a multispecies hierarchical model to investigate whether bird communities differ between two major wetland habitats in boreal Canada: beaver ponds and peatland ponds. In addition to including variables such as forest cover and latitude, we adopted a structural equation model approach to estimate the occupancy of American red squirrels and its potential influence on bird communities. Using automated recording stations deployed at 50 ponds, we detected 96 bird species in 2018 and 2019. Bird species were grouped into four taxonomic guilds according to their habitat successional requirements: early successional species, late successional species, generalists, and wetland species. Beaver ponds harbored higher species richness, a pattern driven primarily by early successional species. The occupancy of almost a quarter of the species was lower in the presence of red squirrels. Late successional species responded positively to the cover of forest surrounding the pond. Our results highlight the value of considering acoustic data of red squirrels to quantify habitat quality in boreal forests. We conclude that beaver activity shapes bird assemblages through modification of their habitat, and that some bird guilds are associated negatively with the presence of American red squirrels.

Les milieux humides et la prédation dans les écosystèmes boréaux jouent des rôles essentiels tout au long de la saison de reproduction pour les assemblages d’oiseaux. Nous avons trouvé une association positive entre les assemblages d’oiseaux et l’activité des castors, mais une association négative entre les assemblages d’oiseaux et la présence des écureuils roux d’Amérique (Tamiasciurus hudsonicus). Nous avons utilisé un modèle hiérarchique multi-espèces pour étudier si les communautés d’oiseaux diffèrent entre deux importants types de milieux humides du Canada boréal: les étangs de castors et les étangs de tourbières. Au-delà de l’inclusion de variables telles que le couvert forestier et la latitude, nous avons adopté une approche de modèle d’équation structurelle pour estimer l’occupation de l’écureuil roux d’Amérique et son influence potentielle sur les communautés d’oiseaux. En utilisant des stations d’enregistrement automatisées déployées sur 50 étangs, nous avons détecté 96 espèces d’oiseaux en 2018 et 2019. Les espèces des oiseaux ont été regroupées en quatre guildes taxonomiques en fonction de leurs exigences en matière de succession de l’habitat: espèces de début de succession, espèces de fin de succession, généralistes et espèces des milieux humides. Les étangs de castors abritaient une plus grande richesse d’espèces et une composition d’espèces différente de celle des étangs de tourbières, une tendance principalement liée aux espèces de début de succession. L’occupation de près d’un quart des espèces était plus faible en présence de l’écureuil roux. Les espèces de fin de succession ont répondu positivement à la couverture de la forêt entourant l’étang. Nos résultats soulignent l’intérêt de considérer les données acoustiques de l’écureuil roux pour quantifier la qualité de l’habitat dans les forêts boréales. Nous concluons que l’activité des castors façonne les assemblages d’oiseaux en modifiant leur habitat et que certaines guildes d’oiseaux sont associées négativement à la présence de l’écureuil roux d’Amérique.

Linear feature networks are the roads, trails, pipelines, and seismic lines developed throughout many commercial boreal forests. These linear features, while providing access for industrial, recreational, silvicultural, and fire management operations, also have environmental implications which involve both the active and non-active portions of the network. Management of the existing linear feature networks across boreal forests would lead to the optimization of maintenance and construction costs as well as the minimization of the cumulative environmental effects of the anthropogenic linear footprint. Remote sensing data and predictive modelling are valuable support tools for the multi-level management of this network by providing accurate and detailed quantitative information aiming to assess linear feature conditions (e.g., deterioration and vegetation characteristic dynamics). However, the potential of remote sensing datasets to improve knowledge of fine-scale vegetation characteristic dynamics within forest roads has not been fully explored. This study investigated the use of high-spatial resolution (1 m), airborne LiDAR, terrain, climatic, and field survey data, aiming to provide information on vegetation characteristic dynamics within forest roads by (i) developing a predictive model for the characterization of the LiDAR-CHM vegetation cover dynamic (response metric) and (ii) investigating causal factors driving the vegetation cover dynamic using LiDAR (topography: slope, TWI, hillshade, and orientation), Sentinel-2 optical imagery (NDVI), climate databases (sunlight and wind speed), and field inventory (clearing width and years post-clearing). For these purposes, we evaluated and compared the performance of ordinary least squares (OLS) and machine learning (ML) regression approaches commonly used in ecological modelling—multiple linear regression (mlr), multivariate adaptive regression splines (mars), generalized additive model (gam), k-nearest neighbors (knn), gradient boosting machines (gbm), and random forests (rf). We validated our models’ results using an error metric—root mean square error (RMSE)—and a goodness-of-fit metric—coefficient of determination (R2). The predictions were tested using stratified cross-validation and were validated against an independent dataset. Our findings revealed that the rf model showed the most accurate results (cross-validation: R2 = 0.69, RMSE = 18.69%, validation against an independent dataset: R2 = 0.62, RMSE = 20.29%). The most informative factors were clearing width, which had the strongest negative effect, suggesting the underlying influence of disturbance legacies, and years post-clearing, which had a positive effect on the vegetation cover dynamic. Our long-term predictions suggest that a timeframe of no less than 20 years is expected for both wide- and narrow-width roads to exhibit ~50% and ~80% vegetation cover, respectively. This study has improved our understanding of fine-scale vegetation dynamics around forest roads, both qualitatively and quantitatively. The information from the predictive model is useful for both the short- and long-term management of the existing network. Furthermore, the study demonstrates that spatially explicit models using LiDAR data are reliable tools for assessing vegetation dynamics around forest roads. It provides avenues for further research and the potential to integrate this quantitative approach with other linear feature studies. An improved knowledge of vegetation dynamic patterns on linear features can help support sustainable forest management.

Structural complexity generated by forest development processes and tree species compositional changes provide key habitat features for vertebrate communities that rely upon tree size and decay processes for foraging, denning or nesting. Complexity of forest structure in old stands could not only be key for harboring increased taxonomic species diversity but also greater functional diversity through more complexity in networks of tree cavity dependent species. Using a nest web approach that hierarchically links cavity-bearing trees with cavity formation agents (natural decay processes and avian excavators) and cavity users (non-excavator species), we compared network characteristics of nest webs along a time since fire gradient in a naturally disturbed boreal mixedwood forest landscape in eastern North America. Since 2003, twelve 24 to 40 ha plots ranging from 61 to more than 245 years after fire were surveyed at the Lake Duparquet Research and Teaching Forest in Abitibi, Quebec, Canada to detect active nesting, and denning cavities. We found that network complexity both in terms of number of vertebrate species and number of interactions among species, increased along the age gradient and was significantly higher in the older stands than predicted by chance. Whereas cavity-nesting communities in old forests used a higher diversity of tree species over a wide range of decay stages, trembling aspen remained a key cavity-bearing tree throughout the age gradient. Woodpeckers were the main cavity formation agents whereas less than 1% of cavities originated from natural decay. The structural development of older forests is thus a driver for functional diversity in cavity-using vertebrate communities through higher interaction richness in nest webs, among cavity-bearing trees, excavators and non-excavating users. The pivotal contribution of the entire gradient of old forest cover types to the overall complexity of nest webs in the boreal mixedwood zone is also a key for the resilience of the cavity-using vertebrate community to natural disturbances. We discuss how such resilience may be compromised by even-aged industrial timber harvesting with short rotations that shifts the age structure of boreal landscapes toward regenerating and young pole forests whereas old forest cover types become below their historical range of variability.

Interspecific interactions can mediate site occupancy of sympatric species and can be a key factor in habitat use patterns. American martens (Martes americana) and Fishers (Pekania pennanti) are two sympatric mesocarnivores in eastern North American forests. Due to their larger size, fishers have a competitive advantage over martens. We investigated site occupancy of martens and fishers in temperate deciduous forests of Québec, an environment modified by forest management and climate change. We formulated hypotheses on the spatial distribution of the studied species based on the knowledge of local trappers and on the scientific literature regarding forest cover composition, habitat fragmentation, and competitive relationships. We used a network of 49 camera traps monitored over two fall seasons to document site occupancy by both species. We used two-species site occupancy models to assess habitat use and the influence of fishers on martens at spatial grains of different sizes. None of the habitat variables that we considered explained site occupancy by fishers. Availability of dense old coniferous stands explained the spatial distribution of martens both at the home range grain size and at the landscape grain size. We identified the characteristics of habitat hotspots based on the knowledge of trappers, which highlighted the importance of stand composition, height, age, and canopy closure. The characteristics of habitat hotspots for martens in temperate deciduous forests refine the habitat suitability model for American martens that was originally developed for boreal forests of Québec.