Paleoecological Investigation of Vegetation, Climate and Fire History in, and Adjacent to, Kootenay National Park, Southeastern British Columbia, Canada

Paleoecological investigation of two montane lakes in the Kootenay region of southeast British Columbia, Canada, reveal changes in vegetation in response to climate and fire throughout the Holocene. Pollen, charcoal, and lake sediment carbon accumulation rate analyses show seven distinct zones at Marion Lake, presently in the subalpine Engelmann Spruce-Subalpine Fir (ESSF) biogeoclimatic (BEC) zone of Kootenay Valley, British Columbia. Comparison of these records to nearby Dog Lake of Kootenay National Park of Canada in the Montane Spruce (MS) BEC zone of Kootenay Valley, British Columbia reveals unique responses of ecosystems in topographically complex regions. The two most dramatic shifts in vegetation at Marion Lake occur firstly in the early Holocene/late Pleistocene in ML Zone 3 (11,010–10,180 cal. yr. B.P.) possibly reflecting Younger Dryas Chronozone cooling followed by early Holocene xerothermic warming noted by the increased presence of the dry adapted conifer, Douglas-fir (Pseudotsuga menziesii) and increasing fire frequency. The second most prominent change occurred at the transition from ML Zone 5 through 6a (∼2,500 cal. yr. B.P.). This zone transitions from a warmer to a cooler/wetter climate as indicated by the increase in western hemlock (Tsuga heterophylla) and subsequent drop in fire frequency. The overall cooling trend and reduction in fire frequency appears to have occurred ∼700 years later than at Dog Lake (∼43 km to the south and 80 m lower in elevation), resulting in a closed montane spruce forest, whereas Marion Lake developed into a subalpine ecosystem. The temporal and ecological differences between the two study sites likely reflects the particular climate threshold needed to move these ecosystems from developed forests to subalpine conditions, as well as local site climate and fire conditions. These paleoecological records indicate future warming may result in the MS transitioning into an Interior Douglas Fir (IDF) dominated landscape, while the ESSF may become more forested, similar to the modern MS, or develop into a grassland-like landscape dependent on fire frequency. These results indicate that climate and disturbance over a regional area can dictate very different localized vegetative states. Local management implications of these dynamic landscapes will need to understand how ecosystems respond to climate and disturbance at the local or ecosystem/habitat scale.

Genomic evidence of an ancient Inland Temperate Rainforest

The disjunct temperate rainforests of the Pacific Northwest of North America (PNW) are characterized by late-successional dominant tree species western redcedar (Thuja plicata) and western hemlock (Tsuga heterophylla). The demographics of these species, along with the PNW rainforest ecosystem in its entirety, have been heavily impacted by the geological and climatic changes the PNW has experienced over the last 5 million years, including mountain orogeny and repeated Pleistocene glaciations. These environmental events have ultimately shaped the history of these species, with inland segments potentially being extirpated during the Pleistocene glaciation. Here, we collect genomic data for both species across their ranges in order to develop multiple demographic models, each reflecting a different hypothesis on how the ecosystem dominant species may have responded to dramatic climatic change. Results indicate that inland and coastal populations in both species diverged an estimated ~2.5 million years ago and experienced a decrease in population size during glaciation, with a subsequent population expansion. Importantly, we found evidence for gene-flow between coastal and inland populations during the mid-Holocene. It is likely that intermittent migration in these species has prevented allopatric speciation. In conclusion, the combination of genomic data and population demographic inference procedures involving machine learning establish that populations of the ecosystem dominants Thuja plicata and Tsuga heterophylla persisted in refugia located in both the coastal and inland regions, with populations expanding and contracting in response to glacial cycles with occasional gene-flow.

The Lari-Leuco Container: A Novel Collection Arena for Separating Insects Ascending or Descending From a Plant Foliage Sample

Efficient separation of insects from plant material for quantification and collection is an important component of entomological research. This paper reports on a novel, easily replicable container designed to efficiently collect two different biological control agents dispersing from hemlock (Tsuga spp.) foliage infested with the invasive hemlock woolly adelgid (HWA), Adelges tsugae Annand (Hemiptera: Adelgidae). The container utilizes a simplified Berlese-style funnel design to collect Laricobius spp. (Coleoptera: Derodontidae) larvae dropping from the foliage into a removable bottom jar, a central jar to house the foliage sample, and a removable top jar to collect adult silver flies (Leucopis spp., Diptera: Chamaemyiidae) emerging from puparia on the twigs. The efficacy of two designs (with and without a funnel leading to the top collection jar) was evaluated using western hemlock [Tsuga heterophylla (Raf.) Sarg.] foliage naturally colonized with HWA and the two predator genera. All Laricobius larvae were effectively collected in the bottom jar, and the addition of an inverted funnel leading to the top collection jar increased the proportion of Leucopis flies reaching the target jar from 60% to 94%. This ‘Lari-Leuco’ container is presented as a research and motoring tool to benefit the integrated pest management program for HWA in eastern North America and for potential use in simultaneously separating ascending and descending life stages in other insect-plant or predator-prey systems.

A species-specific, site-sensitive maximum stand density index model for Pacific Northwest conifer forests

Maximum stand density index (SDIMAX) models were developed for important Pacific Northwest conifers of western Oregon and Washington, USA, based on site and species influences and interactions. Inventory and monitoring data from numerous federal, state, and private forest management groups were obtained throughout the region to ensure a wide coverage of site characteristics. These observations include information on tree size, number, and species composition. The effects and influence on the self-thinning frontier of plot-specific factors such as climate, topography, soils, and geology, as well as species composition, were evaluated based on geographic location using a multistep approach to analysis involving linear quantile mixed models, random forest, and stochastic frontier functions. The self-thinning slope of forest stands dominated by Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) was found to be –1.517 and that of stands dominated by western hemlock (Tsuga heterophylla (Raf.) Sarg.) was found to be –1.461, leading to regionwide modelled SDIMAX values at the 95th percentile of 1728 and 1952 trees per hectare, respectively. The regional model of site-specific SDIMAX will support forest managers in decision-making regarding density management and species selection to more efficiently utilize site resources toward healthy, productive forests.

The importance of large-diameter trees to the creation of snag and deadwood biomass

Background Baseline levels of tree mortality can, over time, contribute to high snag densities and high levels of deadwood (down woody debris) if fire is infrequent and decomposition is slow. Deadwood can be important for tree recruitment, and it plays a major role in terrestrial carbon cycling, but deadwood is rarely examined in a spatially explicit context. Methods Between 2011 and 2019, we annually tracked all trees and snags ≥1 cm in diameter and mapped all pieces of deadwood ≥10 cm diameter and ≥1 m in length in 25.6 ha of Tsuga heterophylla / Pseudotsuga menziesii forest. We analyzed the amount, biomass, and spatial distribution of deadwood, and we assessed how various causes of mortality that contributed uniquely to deadwood creation. Results Compared to aboveground woody live biomass of 481 Mg ha−1 (from trees ≥10 cm diameter), snag biomass was 74 Mg ha−1 and deadwood biomass was 109 Mg ha−1 (from boles ≥10 cm diameter). Biomass from large-diameter trees (≥60 cm) accounted for 85%, 88%, and 58%, of trees, snags, and deadwood, respectively. Total aboveground woody live and dead biomass was 668 Mg ha−1. The annual production of downed wood (≥10 cm diameter) from tree boles averaged 4 Mg ha−1 yr−1. Woody debris was spatially heterogeneous, varying more than two orders of magnitude from 4 to 587 Mg ha−1 at the scale of 20 m × 20 m quadrats. Almost all causes of deadwood creation varied in importance between large-diameter trees and small-diameter trees. Biomass of standing stems and deadwood had weak inverse distributions, reflecting the long period of time required for trees to reach large diameters following antecedent tree mortalities and the centennial scale time required for deadwood decomposition. Conclusion Old-growth forests contain large stores of biomass in living trees, as well as in snag and deadwood biomass pools that are stable long after tree death. Ignoring biomass (or carbon) in deadwood pools can lead to substantial underestimations of sequestration and stability.

Temporal Asynchrony of Adult Emergence Between Leucopis argenticollis and Leucopis piniperda (Diptera: Chamaemyiidae), Predators of the Hemlock Woolly Adelgid (Hemiptera: Adelgidae), with Implications for Biological Control

Two species of silver fly, Leucopis argenticollis (Zetterstedt) and Leucopis piniperda (Malloch) (Diptera: Chamaemyiidae), from the Pacific Northwest region of North America have been identified as potential biological control agents of hemlock woolly adelgid (Hemiptera: Adelgidae: Adelges tsugae Annand) in eastern North America. The two predators are collectively synchronized with A. tsugae development. To determine whether adult emergence of the two species of silver fly are also synchronized with one another, we collected adult Leucopis which emerged from A. tsugae-infested western hemlock [Pinaceae: Tsuga heterophylla (Raf.) Sarg.] from four sites in the Pacific Northwest over a 29-d period. Specimens were collected twice daily in the laboratory and identified to species using DNA barcoding. The study found that more adult Leucopis were collected in the evening than the morning. Additionally, the daily emergences of adults over the 29-d sampling period exhibited sinusoidal-like fluctuations of peak abundance of each species, lending evidence to a pattern of temporal partitioning. This pattern could have logistical implications for their use as biological control agents in eastern North America, namely the need to release both species for maximum efficacy in decreasing A. tsugae populations.

Modelling the Crown Profile of Western Hemlock (Tsuga heterophylla) with a Combination of Component and Aggregate Measures of Crown Size

Research Highlights: We present statistical methods for using crown measurement data from multiple destructive sampling studies to model crown profiles in the Tree and Stand Simulator (TASS) and evaluate it using component (branch-level) and aggregate (tree-level) predictions. Combining data collected under different sampling protocols offered unique challenges. Background and Objectives: The approach to modelling crown profiles was based on Mitchell’s monograph on Douglas-fir growth and simulated dynamics. The functional form defines the potential crown size and shape and governs the rate of crown expansion. With the availability of additional data, we are able to update these functions as part of ongoing TASS development and demonstrate the formulation and fitting of new crown profile equations for stand-grown western hemlock (Tsuga heterophylla (Raf.) Sarg. Materials and Methods: Detailed measurements on 1616 branches from 153 trees were collected for TASS development over a 40-year period. Data were collected under two different sampling protocols and the methods were designed to allow the use of data from both protocols. Data collected on all branches were then introduced through the application of the ratio of length of each of the selected branches to the largest branch in the internode (RL). Results: A mixed-effects model with two random effects, which accounted for tree-level variation, provided the best fit. From that, a model that expressed one parameter as a function of another with one random effect was developed to complement the structure of the Tree and Stand Simulator (TASS). The models generally over-estimated crown size when compared to the projected crown area recorded from field measurements, and a scalar adjustment factor of 0.89 was applied that minimised mean-squared error of the differences. The new model is fit from direct measures of crown radius and predicts narrower crown shapes than previous functions used in TASS.

Deglacial landforms and Holocene vegetation trajectories in the northern interior cedar-hemlock forests of British Columbia

The northern Rocky Mountain Trench of eastern British Columbia is a broad valley mantled by glaciolacustrine terraces supporting a complex mix of mesic-temperate (“interior wetbelt”) forests that are strongly affected by terrain and substrate. Neither the geomorphic history during early-Holocene deglaciation nor the vegetation history of the origin of the Tsuga heterophylla and Thuja plicata populations in the interior wetbelt forest is well understood. Sediment cores were obtained from two lakes, 10 km apart and occupying different terraces (83 m elevational difference) and compared to existing fire-history and paleoclimate reconstructions. Radiocarbon dates and a mapped terrain classification indicate the upper terrace formed as a lacustrine and glaciofluvial kame terrace hundreds of years prior to a lower terrace formed by glaciolacustrine sediments of a proglacial lake. The minimum limiting ages of these terraces correlate with dated jökulhlaup deposits of the Fraser River. The upper site’s first detectable pollen at > 11.0 ka was dominated by light-seeded pioneer taxa (Poaceae, Artemisia, and Populus) followed by a peak in Pinus and finally dominance by Betula at 10.2 ka. Pollen data suggest an earlier invasion of T. heterophylla than previously understood. Wetlands on extensive poorly drained glaciolacustrine soils promoted the persistence of boreal taxa and open forests (e.g., Picea mariana) while the better-drained upper kame terrace promoted development of closed-canopy shade-tolerant taxa. Invasion and expansion of mesic cedar-hemlock taxa progressed since at least the middle Holocene but was highly constrained by edaphic controls.

Plant Water Stress and Soil Depletion in Variable-Density, Red Alder/Western Hemlock Coastal Oregon Plantations

Riparian ecosystems provide critical habitat and functions while being some of the most productive areas in forests. Both conifers and hardwoods contribute to maintenance of habitat and function. To determine the impact of water stress on growth of red alder (Alnus rubra Bong.) and western hemlock (Tsuga heterophylla [Raf.] Sarg.), we installed Nelder type 1a combined with replacement series plots on three Oregon Coast Range sites. Densities ranged from 988 to 85,400 trees/hectare, with ratios (hemlock:alder) of 100:0, 75:25, 50:50, 25:75, and 0:100. In the first 4 years after planting, alder used water in the growing season at greater depths earlier than western hemlock. Higher densities resulted in greater water stress later in the growing season in weeded areas (maintained by herbicide applications), but stress was similar across densities in unweeded areas. Water stress at early ages was correlated with decreased size 14 or 24 years after planting for both species, but these correlations were confounded with other effects of density. Increasing water availability in areas with low summer precipitation could enhance growth of red alder and western hemlock, even in highly productive riparian areas.