Central British Columbia Research Articles

An inventory of rock glaciers in the central British Columbia Coast Mountains, Canada, from high resolution Google Earth imagery

ABSTRACT Little is known about the presence, distribution, age, or activity of rock glaciers in the British Columbia Coast Mountains of western Canada. Reflecting debris accumulation and mass wasting under a periglacial climate, these rock glaciers describe a geomorphic response to permafrost regimes that may or may not presently exist. An inventory of rock glacier landforms in the eastern front ranges of the Coast Mountains, using high-resolution Google Earth imagery, documented 165 rock glaciers between lat. 50°10ʹ and 52°08ʹ N. The majority of these rock glaciers occur at sites positioned between 1,900 and 2,300 m above sea level, where rain shadow effects and continental air masses result in persistent dry, cold conditions. Morphology and field observation suggest that these features contain intact ice. The rock glaciers occupy predominately northwest- to northeast-facing slopes, with talus-derived rock glaciers largely restricted to north-facing slopes. Glacier-derived features outnumber talus-derived features by a ratio of 5:1. Several of the inventoried rock glaciers were located up valley from presumed Younger Dryas terminal moraines, indicating that they formed after 9390 BP. Dendrogeomorphological investigations at one rock glacier record contemporary activity that resulted in 1.3 cm/yr of frontal advance since AD 1674. This inventory is the first to document the presence of rock glaciers in the Coast Mountains and supports preliminary understandings of permafrost distribution in the southwestern Canadian Cordillera.

ПРОРОК БИНИ И РЕЛИГИОЗНАЯ ЛИХОРАДКА У КЭРЬЕРОВ В 1834-1840-Х ГГ

The article deals with the phenomenon of religious fever among the indigenous people of central British Columbia in the second quarter of the 19th century. In a short time, the traditional belief system of the Carrier Indians adopted some of the Christian doctrines and symbolism. This took place in the period before the first missionaries came to the region. Paradoxically, at the initial stage, the enculturation of Christianity started without participation of Christians themselves. The transformation of the Indians’ religious views wasn’t accompanied by European cultural expansion. In the late 1860s, when the missionaries started their active agency, the Indians were ready to accept the new religion without resistance. This article is based on the analysis of folk tales, historical dates and researches devoted to this subject. The channels of cultural diffusion and the integration of elements of Christianity into traditional beliefs are revealed. The psychological reasons of acceptance of the new religion by the Indians are also in the focus of attention. The author determines the date of the beginning of the Prophet's movement among the Carrier. The key source of this dating was John McLean's journal, which has not previously been used in researches of the subjects. The author proposes accounting this cultural phenomenon as a spontaneous movement, not a cult of a religious leader.

Latest Neoproterozoic to Cambrian detrital zircon facies of western Laurentia

Abstract Late Neoproterozoic to Cambrian sandstone units are common in western Laurentia and record initial transgression of the craton after the formation of the western passive margin during the latest Neoproterozoic to earliest Cambrian. Detrital zircon measurements from 42 latest Neoproterozoic to Cambrian basal Sauk sequences and five older Neoproterozoic sandstone samples from a region extending from the Mexico–United States border to central British Columbia, Canada, are combined with previous results to characterize sediment source areas and dispersal systems. Detrital zircon populations in Neoproterozoic and Cambrian sedimentary rocks are divided into six facies based on a statistical comparison using multidimensional scaling. Detrital zircon facies are found in unique geographical regions reflecting proximity to the major tectonic provinces of Laurentia. Samples from northern regions are dominated by Archean and Paleoproterozoic zircons derived from Archean tectonic provinces and the orogenic belts that record the assembly of the Laurentian craton. More southerly sample locations show an increase in detrital zircons derived from younger Paleoproterozoic orogenic belts and early Mesoproterozoic intrusive suites. Detrital zircons from Grenville-aged sources are common in the south. The Transcontinental Arch, a feature interpreted to have controlled large-scale sediment dispersal patterns in the mid- to late Cambrian, likely played a major role in isolating the southern and northern signatures. Our data set can be used to test tectonic models for the Cordilleran orogen that invoke Jurassic or Cretaceous collision of a ribbon continent as the driving mechanism for orogenesis. Cambrian rocks of the Cassiar-Antler platform juxtaposed with North America during the hypothetical ribbon continent collision show the same geographic distribution of detrital zircon facies as similar-aged rocks from autochthonous and parautochthonous locations on the Laurentian margin. The concordance of detrital zircon facies across the proposed suture is a negative result for models that predict large dextral displacements, on the order of 2000 km, across the suture.

Site preparation severity influences lodgepole pine plant community composition, diversity, and succession over 25 years

Lodgepole pine (Pinus contorta var. latifolia Engelm.) ecosystems of central British Columbia face cumulative stresses, and management practices are increasingly scrutinized. We addressed trade-offs between “light-on-the-land” versus more aggressive silvicultural approaches by examining plant communities and indicator species (non-natives, berry producers, epiphytes, mycotrophs, pine rust alternate hosts) across a gradient of five or six site preparation treatments at the Bednesti trial (established 1987). We tested whether more severe site preparation (i) caused plant community composition to diverge from a 35- to 46-year-old reference forest, (ii) accelerated succession by hastening crown closure, or (iii) delayed succession by promoting seral species. Nonmetric multidimensional scaling (NMDS) ordination showed all treatments converging toward the reference forest composition. At 10 years, succession was incrementally delayed by more severe treatments; at 25 years, only burned windrows were still delayed. Mixed-effects models based on site preparation severity were better than crown closure models for 11 of 13 variables tested, suggesting that mostly belowground processes drive succession in these infertile ecosystems. Invasive hawkweeds persisted on all treatments at 25 years. Limited, contradictory data did not support using mechanical or fire treatments to reduce alternate hosts of pine stem rusts. Long-term trials such as Bednesti highlight the need for ecosystem-specific strategies and diverse approaches to accommodate conflicting benefits and risks of disturbance in forests.

Spatially and Temporally Associated Porphyry Deposits with Distinct Cu/Au/Mo Ratios, Woodjam District, Central British Columbia

Abstract The Woodjam district is a cluster of porphyry Cu-Au deposits of Early Jurassic age (~196 Ma) and is located in the Quesnel terrane in central British Columbia. Porphyry centers include the Southeast zone Cu-Mo porphyry, the Deerhorn and Megabuck Au-Cu porphyries, and the Takom and Three Firs Cu-Au porphyries. The Takomkane batholith, which intruded strata of the Nicola Group and is host to the Southeast zone, has characteristics of a calc-alkalic Cu-Mo porphyry. The Deerhorn, Megabuck, and Takom deposits are centered on narrow monzonite bodies with pencil-like geometries that intruded the Nicola Group volcanic sequence. These small volume intrusions have characteristics of high K calc-alkalic intrusions. The Southeast zone, Deerhorn, Megabuck, and Takom deposits have similar ages and their intrusive units can be divided into two groups with distinct geochemical characteristics. Differences between alteration and mineralization of the deposits are attributed to the magmatic evolution of the system as well as differences in the depth of emplacement and preservation. Based on stratigraphic relationships, the Deerhorn and Megabuck deposits are interpreted to be emplaced at the shallowest structural level in the district and have the highest Au/Cu ratios and the lowest temperature alteration assemblages. The Southeast zone is hosted by the felsic units in the Takomkane batholith and represents the deepest parts of the hydrothermal system. Cu-Au mineralization in the Takom deposit shows characteristics intermediate between the Southeast zone and the Deerhorn deposit. Based on stratigraphic interpretations the difference in depth of emplacement between the Southeast zone and Megabuck is about 1,400 m, with Takom being emplaced about 400 m above the Southeast zone. The Woodjam district illustrates the variety of styles of porphyry mineralization that can form over a restricted time interval (<1 m.y.) within a single district.

Quaternary geology of part of the McLeod Lake map area (NTS 093J), central British Columbia: lithostratigraphy, glacial history, and chronology

The coalescence and subsequent divergence of glaciers near the McLeod Lake map area during the last (Fraser) glaciation are recorded by landforms, ice-flow indicators, and the distribution of sediments. Ice initially flowed into the study area from the northwest, with at least one fluctuation in the ice-front position. Ice flow during glacial maximum was generally to the northeast and transitioned to the east during deglaciation. The Quaternary stratigraphic record spans the Fraser Glaciation and is represented by a sequence of advance glaciolacustrine sediments, multiple till units, retreat glaciolacustrine and glaciofluvial sediments, and associated postglacial aeolian material. The surficial geology is dominated by thick, streamlined till that thins where relief is high. Glaciofluvial outwash occurs mostly in northeast- and southeast-trending meltwater channels, while ice-contact glaciofluvial deposits and ablation till occur in depressions throughout the region. Extensive glaciolacustrine deposits blanket low-lying regions in the southern parts of the study area. Meltwater and ice-flow features suggest that deglaciation in the region was dominantly frontal retreat and that ice was largely active as it retreated through the study area. Postglacial aeolian activity was brief; optical dating on K-feldspar from aeolian landforms indicates that the landforms had stabilized by between 8.71–10.71 and 12.3–14.3 ka, and provides minimum ages for ice retreat. The data from the study area are evaluated with similar data from adjacent regions to develop an ice-flow history, and refine the conceptual model of deglaciation for the northern Interior Plateau.

Mixed-severity fire history at a forest-grassland ecotone in west central British Columbia, Canada.

This study examines spatially variable stand structure and fire-climate relationships at a low elevation forest-grassland ecotone in west central British Columbia, Canada. Fire history reconstructions were based on samples from 92 fire-scarred trees and stand demography from 27 plots collected over an area of about 7km2 . We documented historical chronologies of widespread fires and localized grassland fires between AD 1600 and 1900. Relationships between fire events, reconstructed values of the Palmer Drought Severity Index, and annual precipitation were examined using superposed epoch and bivariate event analyses. Widespread fires occurred during warm, dry years and were preceded by multiple anomalously dry, warm years. Localized fires that affected only grassland-proximal forests were more frequent than widespread fires. These localized fires showed a lagged, positive relationship with wetter conditions. The landscape pattern of forest structure provided further evidence of complex fire activity with multiple plots shown to have experienced low-, mixed-, and/or high-severity fires over the last four centuries. We concluded that this forest-grassland ecotone was characterized by fires of mixed severity, dominated by frequent, low-severity fires punctuated by widespread fires of moderate to high severity. This landscape-level variability in fire-climate relationships and patterns in forest structure has important implications for fire and grassland management in west central British Columbia and similar environments elsewhere. Forest restoration techniques such as prescribed fire and thinning are oftentimes applied at the forest-grassland ecotone on the basis that historically high frequency, low-severity fires defined the character of past fire activity. This study provides forest managers and policy makers with important information on mixed-severity fire activity at a low elevation forest-grassland ecotone, a crucial prerequisite for the effective management of these complex ecosystems.

An integrated modeling approach for assessing management objectives for mule deer in central British Columbia

ABSTRACTWe used an integrated Bayesian state‐space population model to assess whether management objectives were met before (1995–2003), during (2004–2010), and after (2011–2013) antlerless permits to harvest mule deer (Odocoileus hemionus) were increased in response to stakeholder concerns in central British Columbia, Canada. Data inputs included 19 years of harvest data, 7 years of autumn age–sex composition data, 17 years of spring age–sex composition data, and 15 years of a population index. Management objectives were to maintain a spring population of 7,000–9,000 deer and a posthunt adult sex ratio of 20–30 males:100 females. An 8.5‐fold increase in antlerless permits raised the antlerless harvest rates from 1.0% (1995–2003) to 4.8% (2004–2010). Antlerless harvest rates decreased to 2.8% following a 53% decrease in permits from 2011 to 2013. Population projections from 2014 to 2018 fell within the bounds of the management objectives, but 95% credibility intervals revealed great uncertainty in population size and composition. We recommend a structured, adaptive approach to mule deer management that includes annual adjustment of harvests, monitoring, and modeling, with an open‐ended stakeholder engagement process to ensure objectives remain relevant, measurable, and achievable. © 2017 The Wildlife Society.

Interannual climate variability drives regional fires in west central British Columbia, Canada

We investigated the influence of climate variability on forest fire occurrence at eight sites in west central British Columbia, Canada. Forty-six local fire years affecting a single site and 16 moderate fire years affecting two or more sites were identified (1600–1900 A.D.). Existing fire history data were incorporated to identify 17 regionally synchronous fire years (fires that affected ≥3 sites). Interannual and multidecadal relationships between fire occurrence and the Palmer Drought Severity Index (PDSI), El Nino–Southern Oscillation (ENSO), Pacific Decadal Oscillation (PDO), and the Pacific North American (PNA) pattern were examined, in addition to the effects of additive positive phases of ENSO and PDO. We examined multiple reconstructions of ENSO, PDO, and PNA and utilized three methodological approaches to characterize climate-fire relationships. We found that the influence of interannual climate, expressed as PDSI, increasingly synchronized the occurrence of fires when examined from local to regional scales. An association between local fires and positive antecedent moisture conditions suggests moisture-driven fine fuel development and the proximity of some sites to grasslands likely function as key determinants of local-scale fire activity. The relationships between regional fires and ENSO, PDO, and PNA suggest that large-scale patterns of climate variability exert a weak and/or inconsistent influence over fire activity in west central British Columbia between 1700 and 1900 A.D. Although inconsistent among reconstructions of climate patterns, we identified a significant relationship between regional fires and large-scale climate patterns when ENSO and PDO were both in positive phases.

Planning for Winter Road Maintenance in the Context of Climate Change

Abstract Winter weather creates mobility challenges for most northern jurisdictions, leading to significant expenditures on winter road maintenance (WRM) activities. While the science and practice of snow and ice control is continually evolving, climate change presents particular challenges for the strategic planning of WRM. The purpose of this study is 1) to develop a winter severity index (WSI) to better understand how winter weather translates into interannual variations in WRM activities and 2) to apply the WSI to future climate change projections to assist a northern community in preparing for climate change. A new method for creating a WSI model is explored, using readily available data from maintenance records and meteorological stations. The WSI is created by optimizing values for three levels of snowfall as well as potential icing events and is shown to have high predictive accuracy for WRM (coefficient of determination R2 of 0.93). The WSI is then applied to historic and future climate data in a municipality located in central British Columbia, Canada. Findings reveal that much of the variability in WRM can be attributed to weather. The results of the climate change analysis show that winter precipitation in the region is expected to increase by 5.2%–12.3%, and winter average temperatures are projected to increase by 1.5°–2.8°C in the 2050s, compared to the 1976–2000 baseline based on 65 GCMs. Based on the midrange (25th to 75th percentiles) of the 65 GCM projections, annual demand for WRM activities is estimated to decrease by 13.0%–22.0%.

Volumetric fault image enhancement — Some applications

The interpretation of faults on 3D seismic data is often aided by the use of geometric attributes such as coherence and curvature. Unfortunately, these same attributes also delineate stratigraphic boundaries (geologic signal) and apparent discontinuities due to cross-cutting seismic noise. Effective fault mapping thus requires enhancing piecewise continuous faults and suppressing stratabound edges and unconformities as well as seismic noise. To achieve this objective, we apply two passes of edge-preserving structure-oriented filtering followed by a recently developed fault enhancement algorithm based on a directional Laplacian of a Gaussian operator. We determine the effectiveness of this workflow on a 3D seismic volume from central British Columbia, Canada.

A glacial readvance during retreat of the Cordilleran Ice Sheet, British Columbia central coast

AbstractDescriptions of the Cordilleran Ice Sheet retreat after the last glacial maximum have included short-lived readvances occurring during the Older and Younger Dryas stadial periods and into the Holocene, but identification of these events has been largely limited to southwest and central British Columbia and northwest Washington State. We present evidence of a late Pleistocene readvance of Cordilleran ice occurring on the central coast of British Columbia on Calvert Island, between northern Vancouver Island and Haida Gwaii. Evidence is provided by sedimentological and paleoecological information contained in a sedimentary sequence combined with geomorphic mapping of glacial features in the region. Results indicate that a cold climate existed between 15.1 and 14.3 cal ka BP and that ice advanced to, and then retreated from, the western edge of the island between 14.2 and 13.8 cal ka BP. These data provide the first evidence of a major fluctuation in the retreating ice sheet margin in this region and suggest that a cold climate was a major factor in ice readvance. These data contribute to the understanding of past temperature, ice loading and crustal response, the nature of ice margin retreat, and the paleoenvironment of an understudied area of the Pacific Northwest.

Territorial black-capped chickadee males respond faster to high- than to low-frequency songs in experimentally elevated noise conditions.

Low-frequency urban noise can interfere with avian communication through masking. Some species are able to shift the frequency of their vocalizations upwards in noisy conditions, which may reduce the effects of masking. However, results from playback studies investigating whether or not such vocal changes improve audibility in noisy conditions are not clear; the responses of free-ranging individuals to shifted signals are potentially confounded by functional trade-offs between masking-related audibility and frequency-dependent signal quality. Black-capped chickadees (Poecile atricapillus) naturally sing their songs at several different frequencies as they pitch-shift to match conspecifics during song-matching contests. They are also known to switch to higher song frequencies in response to experimental noise exposure. Each male produces both high- and low-frequency songs and absolute frequency is not a signal of aggression or dominance, making this an interesting species in which to test whether higher-frequency songs are more audible than lower-frequency songs in noisy conditions. We conducted playback studies across southern and central British Columbia, Canada, using paired song stimuli (high- vs low-frequency songs, n = 24 pairs) embedded in synthetic background noise created to match typical urban sound profiles. Over the course of each playback, the signal-to-noise ratio of the song stimuli was gradually increased by raising the amplitude of the song stimuli while maintaining background noise at a constant amplitude. We evaluated variation in how quickly and aggressively territorial males reacted to each of the paired stimuli. We found that males responded more quickly to playbacks of high- than low-frequency songs when high-frequency songs were presented first, but not when low-frequency songs were first. This difference may be explained by high-frequency songs being more audible combined with a carry-over effect resulting in slower responses to the second stimulus due to habituation. We observed no difference in overall aggression between stimuli. These results suggest that high-frequency songs may be more audible under noisy conditions.

An investigation into the contrasting growth response of lodgepole pine and white spruce to harvest-related soil disturbance

Losses in forest productivity through poor soil management are typically evaluated by changes in crop production, but conflicting growth responses among co-occurring species can challenge criteria for sustainability. In this study, we evaluate species response to compaction and organic matter removal by contrasting the growth and foliar attributes at age 20 of lodgepole pine (Pinus contorta Dougl. ex Loud.) and hybrid white spruce (Picea glauca × engelmannii (Moench) Voss) in sub-boreal forests of central British Columbia. Bole volume increment was stable for lodgepole pine, averaging a 6 % difference across treatments. White spruce, in contrast, often had large growth increases on forest floor retained – compacted plots (average of 63 % gain) but reductions in productivity on the forest floor removed – compacted plots (40 % decline). Foliar nitrogen concentrations converged across treatments for both species, despite the removal of over 50 % of the site N capital. Some key differences in species traits were suggested by foliar δ15N, attributed to deeper rooting of pine, and foliar δ13C, indicating possibly higher CO2 assimilation potential for spruce. The implication is that standard metrics of sustainable forestry are contingent upon tree autecology and that comprehensive assessments of land management require a measure of collective species response.

Twenty years of ecosystem response after clearcutting and slashburning in conifer forests of central British Columbia, Canada.

Forests are being clearcut over extensive areas of western North America, but plant community response to harvesting and slashburning under varying climatic conditions in central British Columbia, Canada is still largely unknown. Evaluation of resilience is hampered by the short history of logging, lack of long-term experiments and methodological limitations. To test the effect of clearcut logging, prescribed burning and reforestation on forest resilience, we recorded vascular plant cover repeatedly after treatment between 1981 and 2008 in 16 permanent research installations in three biogeoclimatic zones: Engelmann Spruce-Subalpine Fir, Interior Cedar-Hemlock and Sub-Boreal Spruce. We created a plant-trait dataset for the 181 recorded species to define plant functional types representing groups of plants that behave in similar ways and/or produce similar ecological outcomes. These plant functional types, along with taxonomic analysis of diagnostic and indicator species, were then used to evaluate plant community response to disturbance. Twenty years post-treatment, species diversity increased in all zones and plant abundance was greatest in the Interior Cedar-Hemlock. Cover of understory plant functional types associated with mature conifer forests increased in all zones, constituting a significant proportion (> 40%) of the vegetation community by year 20. Response patterns varied by zone and with time. Understory species diagnostic of mature forests were present in all zones by year 20, but we identified indicator species sensitive to slashburning or requiring more time for recovery, including white-flowered rhododendron (Rhododendron albiflorum) and devil's club (Oplopanax horridus). Overall, loss of compositional or functional diversity following harvest and site remediation was not detected, suggesting that montane and subalpine forests in British Columbia are resilient to this treatment. However, because these forests can be slow to recover from disturbance, the post-disturbance assessment window of this study may not have been long enough to detect diminishment of ecosystem resilience.

Conifer Diet Choices Made by Mule Deer (Odocoileushemionus) of North Central British Columbia During a Cafeteria-Style Feeding Trial

The winter diets of northern mule deer are comprised predominantly of shrubs and trees. Included in winter diets are various species of conifers, some of which are important forestry and ornamental crop trees. Diet choices and feeding habits of mule deer utilizing conifer crops, however, have not been widely reported. Here, we document diet choices made by human-habituated mule deer from north central British Columbia during an exploratory cafeteria-style feeding trial using boughs of thirteen species of native North American and introduced European conifers. Mule deer consumed more biomass from exotic conifers (specifically pine) than from rarely encountered native conifers, and consumed the least biomass from those native conifers commonly encountered by deer in northern BC, a finding which appears to lend support to the Biotic Resistance Hypothesis. From most to least biomass consumed, deer selected Mugho pine, Scots pine, western hemlock, subalpine fir, Douglas fir, Norway spruce, Ponderosa pine, western white pine, lodgepole pine, common juniper, red cedar, black spruce, and hybrid white spruce. Analysis of Variance testing suggested that deer preferred to eat exotic Mugho and Scots pine and avoided eating commonly encountered, native black and hybrid white spruce. A relationship also existed between the biomass consumed and the number of bites taken by deer from each species, which generally indicated that deer took more bites from exotic and rarely encountered conifers, but not always. Our results suggest that northern mule deer, therefore, may be more likely to damage introduced pine and other exotics relative to native species planted in either ornamental or forest plantation settings when deer densities are high and/or browse resources are limited.

Differentiation of Alternate Harvesting Practices Using Annual Time Series of Landsat Data

Sustainable forest management practices allow for a range of harvest prescriptions, including clearcut, clearcut with residual, and partial or selective cutting, which are largely distinguished by the amount of canopy cover removed. The different prescriptions are aimed to emulate natural disturbance, encourage regeneration (seed trees), or offer other ecosystem services, such as the maintenance of local biodiversity or habitat features. Using remotely sensed data, stand-replacing disturbance associated with clearcutting is commonly accurately detected. Novel time series-based change detection products offer an opportunity to determine the capacity to detect and label a wider range of harvest practices. In this research, we demonstrate the capacity of time series imagery, spectral metrics, and related attributed change products, to distinguish between different harvesting practices over a study area in central British Columbia, Canada. Producer’s accuracy of harvest attribution was 79%, with 93% of harvest blocks >5 ha accurately identified. In relation to the amount of canopy cover removed, clearcut harvesting was the most accurately classified (84%), followed by clearcut with residual (79%), and partial cut (64%). Applying detailed spectral metrics derived from Landsat data revealed clearcut and partial cuts to be spectrally distinct. The annual nature of the Landsat time series also offers spatial harvest information within typical, often decadal, forest inventory update cycles. The statistically significant (p < 0.05) relationship between harvest practices and Landsat spectral information indicates a capacity to add increased attribution richness to remote sensing depictions of forest harvest.

Resource Management through Machine Learning

In the modern era of diminishing returns on fixed exploration budgets, challenging targets, and ever-increasing numbers of multi-parameter datasets, proper management and integration of available data is a crucial component of any resource exploration program. Machine learning algorithms have successfully been used for years by the technology sector to accomplish just this task on their databases, and recent developments aim at appropriating these successes to the field of natural resource exploration. Numerous algorithms have been attempted for resource prospectivity mapping in the past, and in this paper we apply a modified support-vector machine algorithm to a test dataset from the QUEST region in central British Columbia, Canada, to target undiscovered Cu-Au porphyry districts. The modified algorithm is designed to properly handle the highly variable uncertainty associated with both the training data (ie: geophysics, geochemistry, geological mapping) as well as the training labels (known Cu-Au porphyry targets in the region). Support vector machines are introduced, the challenges of working with geoscientific datasets are discussed, and finally results from applying the modified algorithm to the QUEST dataset are presented.

Population genetic isolation and limited connectivity in the purple finch (Haemorhous purpureus).

Using a combination of mitochondrial and z‐linked sequences, microsatellite data, and spatio‐geographic modeling, we examined historical and contemporary factors influencing the population genetic structure of the purple finch (Haemorhous purpureus). Mitochondrial DNA data show the presence of two distinct groups corresponding to the two subspecies, H. p. purpureus and H. p. californicus. The two subspecies likely survived in separate refugia during the last glacial maximum, one on the Pacific Coast and one east of the Rocky Mountains, and now remain distinct lineages with little evidence of gene flow between them. Southwestern British Columbia is a notable exception, as subspecies mixing between central British Columbia and Vancouver Island populations suggests a possible contact zone in this region. Z‐linked data support two mitochondrial groups; however, Coastal Oregon and central British Columbia sites show evidence of mixing. Contemporary population structure based on microsatellite data identified at least six genetic clusters: three H. p. purpureus clusters, two H. p. californicus clusters, and one mixed cluster, which likely resulted from high site fidelity and isolation by distance, combined with sexual selection on morphological characters reinforcing subspecies differences.

Are latrine sites an accurate predictor of seasonal habitat selection by river otters (Lontra canadensis) in freshwater systems?

Temporal variation in the availability of food resources is a likely driving factor influencing the distribution and habitat use of river otters (Lontra canadensis). Although latrine sites are commonly used to determine habitat selection, it is unclear if latrine sites are an accurate predictor or even a useful indicator of the seasonal habitat use and distribution of river otters. We apply resource selection functions (RSF) to both latrine and telemetry locations to investigate whether latrine sites identified along lake shorelines during the ice-free season are appropriate predictors of otter habitat selection along shorelines during the ice-free and ice-cover seasons in central British Columbia, Canada. We found that the top models describing otter latrine sites and telemetry locations during the ice-free season were similar. The top RSF models and associated coefficients for the ice-cover season differed, however, with otter presence being positively influenced by shallower water depths. For the spatial extrapolation of averaged RSF coefficients, we found that 21.4 and 69.3 % of predicted latrine habitat along lake shorelines overlapped with ice-cover and ice-free habitat generated from telemetry locations, respectively. The location and activity at latrine sites appear to be a useful method for monitoring otter distribution and habitat use during the ice-free, but not during the ice-cover season. The results of our RSF analyses as well as home range measurements of otters in our study area suggest that cold temperatures and ice cover may be a limiting factor for the distribution of otter populations at northern latitudes.