Geology topics

Research about British Columbia

Source-linked reports with geographic coverage including British Columbia.

At least 19 recordsLinked to original sources

USGS2024

Gene flow prevents genetic diversity loss despite small effective population size in fragmented grizzly bear (Ursus arctos) populations

Genetic monitoring is important in small, fragmented populations that rely on gene flow to maintain genetic diversity. The Selkirk, Yaak, and Cabinet grizzly bear ( Ursus arctos ) populations are among the smallest in North America and are near the southernmost extent of the species’ range. These populations received little to no effective migration for generations but have recently experienced increased gene flow through natural migration and a population augmentation program. A long-term dataset of grizzly bear microsatellite genotypes from 1973 to 2021 presented a unique opportunity to examine genetic trends in these populations over time. We used this dataset of 464 bears to evaluate if gene flow affected observed heterozygosity ( H O ), expected heterozygosity ( H E ), allelic richness ( A R ), and average pairwise relatedness ( r ) in each of these populations. We also estimated effective population size ( N e ) using the temporal and linkage disequilibrium (LD) methods. Post gene flow, A R increased in the Selkirk and Cabinet populations and r decreased in all three populations. We did not observe any significant changes in H E or H O , but H E values in our populations were significantly higher than those estimated using a model without gene flow. Our N e estimates were consistent between the temporal and LD methods and ranged from 15.2 to 15.8, 15.4–17.5, and 5.6–8.9 for the Selkirk, Yaak, and Cabinet populations, respectively. Overall, our findings indicate that gene flow is increasing or maintaining genetic diversity in these populations. However, N e remains low and additional connectivity or augmentation may be needed, particularly in the Cabinet population.

Idaho, Montana, Washington↗

USGS2024

Combined high rates of alternative breeding strategies unexpectedly found among populations of a solitary nesting raptor

Social monogamy is the prevalent mating system in birds, but alternative strategies of extra-pair paternity (EPP) and conspecific brood parasitism (CBP) occur in many species. Raptors are virtually absent in discussions of broad taxonomic reviews regarding EPP and CBP likely because these strategies are mostly absent or at low frequency; CBP is unreported in solitary nesting raptors. In contrast, we found high frequencies of EPP (16%–31%) and CBP (15%–26%) nests among three populations of Cooper's Hawks ( Accipiter cooperii ) across the northern breeding range of this solitary nesting, socially monogamous species. EPP and CBP combined occurred in 42%–46% of all nests among populations and hence unexpectedly were nearly equivalent to proportions of genetically monogamous nests. Select covariates failed to predict presence of EPP and CBP in part because virtually all extra-pair adults were uncaught and likely were floaters. We found no support for the hypothesis that territorial females traded copulations for food to maximize energy intake for increased production. Our unique discoveries enhance knowledge of the extent and diversity of alternative breeding strategies among groups of avian and other animal species.

North Dakota, Wisconsin↗

USGS2024

Effect of a dense inflow on the stratification of a steep-side lake

We detail the effect of a small stream of dense inflow that significantly altered the stratification and water quality in a constructed water body in northern British Columbia, Canada. As the dense inflow passed through the epilimnion of the steep-sided lake, it entrained relatively large quantities of water. The resulting mixture of dense inflow and entrained epilimnetic water sank to the bottom of the lake. The removal of water from the epilimnion due to entrainment reduced the epilimnetic thickness. This opposes the normal process of epilimnetic deepening due to wind and convective cooling. The flux of fluid entrained into the dense inflow was calculated to be between 4 and 14 times the inflow, depending primarily on the thickness of the epilimnion. The entrainment had four major effects: (1) it reduced the residence time of the epilimnion from half a year to less than a month; (2) it removed the freshwater cap that resulted from spring ice melt; (3) it enabled fall turnover, which further enhanced deep oxygen content and helped to prevent meromixis from developing in the lake; and (4) it produced a rapid decline in contaminant (i.e., zinc) concentrations in the epilimnion, which received dissolved metals inputs from oxidized sulfide minerals exposed in the subaerial walls of the lake. Given the wide variety of inflows to inland water bodies, some of which are at least seasonally dense, an understanding of the mechanisms detailed here can inform lake management in general, and more specifically, management of water quality in mine-impacted water bodies.

British Columbia↗

USGS2024

Climate and landform interact to control the source and transport of nitrate in Pacific Northwest rivers

The hydrological effects of climate change are documented in many regions; however, climate-driven impacts to the source and transport of river nutrients remain poorly understood. Understanding the factors controlling nutrient dynamics across river systems is critical to preserve ecosystem function yet challenging given the complexity of landscape and climate interactions. Here, we harness a large regional dataset of nitrate (NO 3 – ) yield, concentration, and isotopic composition (δ 15 N and δ 18 O) to evaluate the strength of hydroclimate and landscape variables in controlling the seasonal source and transport of NO 3 – . We show that hydroclimate strongly influenced the seasonality of river NO 3 – , producing distinct, source-dependent NO 3 – regimes across rivers from two mountain ranges. Riverine responses to hydroclimate were also constrained by watershed-scale topographic features, demonstrating that while regional climate strongly influences the timing of river NO 3 – transport, watershed topography plays a distinct role in mediating the sensitivity of river NO 3 – dynamics to future change.

Washington↗

USGS2023

High-precision ID-TIMS U-Pb geochronology of perovskite (CaTiO3) from the Ice River Complex, southeastern British Columbia

Uranium‑lead perovskite in situ geochronology has become a cornerstone technique for determining the emplacement timing of alkaline, ultrapotassic, and silica-undersaturated igneous rocks, kimberlites, and carbonatites. Accurate in situ dates are dependent on the availability of matrix matched mineral reference materials which themselves are chemically well characterized and dated accurately to the highest possible precision. When dating perovskite to high precision, such as by isotope dilution thermal ionization mass spectrometry (ID-TIMS), appropriately accounting and correcting for the quantity and isotopic composition of Pb incorporated into the crystal upon crystallization (Pb i ) is a large source of uncertainty and potential inaccuracy. Unfortunately, although ultra-high precision perovskite dates are attainable with modern mass spectrometry techniques, the accuracy of applied Pb i compositions, which can be a considerable percentage of total Pb in a crystal, has not kept pace, resulting in percent level inaccuracy on precisely measured isotopic ratios. In an effort to characterize the age and initial Pb isotopic composition of a readily available, relatively pure perovskite endmember (93–98%) which will be useful as a matrix matched age reference material for in situ U-Pb geochronology, we date crystals isolated from three samples from the Ice River Complex (samples 81IR6, I90.3, and I92.30) by ID-TIMS and evaluate their suitability as known-age reference materials. We directly determine the isotopic composition of Pb i in each sample by ID-TIMS measurement of cogenetic low 238 U/ 204 Pb (μ < 500) phlogopite, apatite, and/or clinopyroxene. Using these initial common Pb isotopic compositions, which are significantly more radiogenic than those predicted by age appropriate Pb evolution models, we obtain ultra-precise weighted-mean 206 Pb/ 238 U perovskite dates of 355.83 ± 0.14, 355.04 ± 0.15, and 357.34 ± 0.12 Ma on the three Ice River samples, respectively. These dates are consistent with the relative emplacement ages previously established for units in the Ice River layered ultramafic series based on field relationships. They also demonstrate the feasibility of attaining accurate perovskite dates in instances when maximum precision is critical by utilizing ID-TIMS analysis of cogenetic low uranium phases for direct measurement of Pb i .

British Columbia↗

USGS2022

Coupling validation effort with in situ bioacoustic data improves estimating relative activity and occupancy for multiple species with cross-species misclassifications

The increasing complexity and pace of ecological change requires natural resource managers to consider entire species assemblages. Acoustic recording units (ARUs) require minimal cost and effort to deploy and inform relative activity, or encounter rates, for multiple species simultaneously. ARU-based surveys require post-processing of the recordings via software algorithms that assign a species label to each recording. The automated classification process can result in cross-species misidentifications that should be accounted for when employing statistical modelling for conservation decision-making. Using simulation and ARU-based detection counts from 17 bat species in British Columbia, Canada, we investigate three strategies for adjusting statistical inference for species misclassification: (a) ‘coupling’ ambiguous and unambiguous detections by validating a subset of survey events post-hoc, (b) using a calibration dataset on the software algorithm's (in)accuracy for species identification or (c) specifying informative Bayesian priors on classification probabilities. We explore the impact of different Bayesian prior specifications for the classification probabilities on posterior estimation. We then consider how the quantity of data validated post-hoc impacts model convergence and resulting inferences for bat species relative activity as related to nightly conditions and yearly site occupancy after accounting for site-level environmental variables. Coupled methods resulted in less bias and uncertainty when estimating relative activity and species classification probabilities relative to calibration approaches. We found that species that were difficult-to-detect and those that were often inaccurately identified by the software required more validation effort than more easily detected and/or identified species. Our results suggest that, when possible, acoustic surveys should rely on coupled validated detection information to account for false-positive detections, rather than uncoupled calibration datasets. However, if the assemblage of interest contains a large number of rarely detected or less prevalent species, an intractable amount of effort may be required, suggesting there are benefits to curating a calibration dataset that is representative of the observation process. Our findings provide insights into the practical challenges associated with statistical analyses of ARU data and possible analytical solutions to support reliable and cost-effective decision-making for wildlife conservation/management in the face of known sources of observation errors.

British Columbia↗

USGS2022

Northern Cascadia Margin gas hydrates — Regional geophysical surveying, IODP drilling leg 311, and cabled observatory monitoring

This article reviews extensive geophysical survey data, ocean drilling results and long-term seafloor monitoring that constrain the distribution and concentration of gas hydrates within the accretionary prism of the northern Cascadia subduction margin, located offshore Vancouver Island in Canada. Seismic surveys and geologic studies conducted since the 1980s have mapped the bottom simulating reflector (BSR), detected gas hydrate occurrence and estimated gas hydrate and free gas concentrations. Additional constraints were obtained from seafloor-towed, controlled-source electromagnetic surveying. A component of these studies has been the examination of low-temperature seafloor vents and seeps that emit gas and fluids into the ocean. These features are identified seismically as chimney-like zones of reduced acoustic reflectivity within the sediment stratigraphy, functioning as conduits for gas and fluid migration from below the BSR to the seafloor. Gas hydrates have been recovered from the seafloor and from sediment cores at vent sites, mostly in massive (nodular) form and as a vein-like fracture filling. The Ocean Networks Canada cabled NEPTUNE observatory has gathered extensive continuous, long-term observations on gas hydrate dynamics at the seafloor and in boreholes at two nodes on the continental slope featuring high gas hydrate concentrations. Measurements taken at the observatory include a time-series of gas bubble emission rates, changes in the near-seafloor electromagnetic structure and seafloor compliance linked to gas hydrate formation and dissociation. Two Integrated Ocean Drilling Program (IODP) expeditions collected cores, measured downhole properties and deployed downhole instruments within the central accretionary prism. At IODP Site U1364, pore pressures are being monitored above and below the base of the gas hydrate stability zone at a slope setting using an “Advanced Circulation Obviation Retrofit Kit” (A-CORK). Downhole pore pressures, temperatures and electrical resistivities also are being monitored at IODP Site U1416 using the “Simple Cabled Instrument for Measuring Parameters In Situ” (SCIMPI) tool at a vent site from near-seafloor to just above the base of the gas hydrate stability zone.

British Columbia↗

USGS2021

Paleohydrological context for recent floods and droughts in the Fraser River Basin, British Columbia, Canada

The recent intensification of floods and droughts in the Fraser River Basin (FRB) of British Columbia has had profound cultural, ecological, and economic impacts that are expected to be exacerbated further by anthropogenic climate change. In part due to short instrumental runoff records, the long-term stationarity of hydroclimatic extremes in this major North American watershed remains poorly understood, highlighting the need to use high-resolution paleoenvironmental proxies to inform on past streamflow. Here we use a network of tree-ring proxy records to develop 11 subbasin-scale, complementary flood- and drought-season reconstructions, the first of their kind. The reconstructions explicitly target management-relevant flood and drought seasons within each basin, and are examined in tandem to provide an expanded assessment of extreme events across the FRB with immediate implications for water management. We find that past high flood-season flows have been of greater magnitude and occurred in more consecutive years than during the observational record alone. Early 20th century low flows in the drought season were especially severe in both duration and magnitude in some subbasins relative to recent dry periods. Our Fraser subbasin-scale reconstructions provide long-term benchmarks for the natural flood and drought variability prior to anthropogenic forcing. These reconstructions demonstrate that the instrumental streamflow records upon which current management is based likely underestimate the full natural magnitude, duration, and frequency of extreme seasonal flows in the FRB, as well as the potential severity of future anthropogenically forced events.

British Columbia↗

USGS2021

Earthquake magnitudes from dynamic strain

Dynamic strains have never played a role in determining local earthquake magnitudes, which are routinely set by displacement waveforms from seismic instrumentation (e.g., M L ⁠ ). We present a magnitude scale for local earthquakes based on broadband dynamic strain waveforms. This scale is derived from the peak root‐mean‐squared strains ( ⁠ A ⁠ ) in 4589 records of dynamic strain associated with 365 crustal earthquakes and 77 borehole strainmeters along the Pacific‐North American plate boundary on the west coast of the United States and Canada. In this data set, catalog moment magnitudes range from 3.5≤ M w ≤ 7.2 ⁠ , and hypocentral distances range from 6≤ R ≤500 km . The 1D representation of geometrical spreading and attenuation of A common to all strain data is log A 0 ( R )=−0.00072R−1.45log(R) . After correcting for instrument gain, site terms, and event terms, the magnitude scale, M DS =log A −log A 0 ( R )−log(3×10−9) ⁠ , scales as ≈0.92 M w with a residual standard deviation of 0.19. This close association with M w holds for events east of the −124° meridian; west of this boundary, however, a constant correction of 0.41 is needed to adjust for additional along‐path attenuation effects. As a check on the accuracy of this magnitude scale, we apply it to dynamic strain records from three strainmeters located in the near field of the 2019 M 6.4 and 7.1 Ridgecrest earthquakes. Results from these six records are in agreement to within 0.5 magnitude units, and five out of six records are in agreement to within 0.34 units.

California, Oregon, Washington↗

USGS2021

Channel-amphitheatre landforms resulting from liquefaction flowslides during rapid drawdown of glacial Lake Fraser, British Columbia, Canada

Unusual channel-amphitheatre landforms are present in Late Pleistocene–early Holocene, subaqueous fan and delta deposits in the glacial Lake Fraser basin, central British Columbia. The lake formed during the decay of the last Cordilleran Ice Sheet and drained ~11,500 years ago during a large outburst flood. The fronts of a delta and two subaqueous fans consisting of silt to fine sand are marked by branching networks of wide, nearly flat channels that terminate upstream in digitate, steep-walled amphitheatres. We propose that these channel-amphitheatre landforms formed by liquefaction flowslides that were induced by the rapid drawdown of glacial Lake Fraser during the outburst flood. Similar geomorphic forms, which we believe also to be associated with rapid drawdowns of large Late Pleistocene–early Holocene lakes, occur elsewhere in North America. A recent tailings dam failure and an intentional breaching of a 100-year-old hydroelectric dam provide insights into the processes responsible for the landforms. By using geomechanical analysis, we show how rapid lake drawdown can trigger liquefaction flowslides in deposits of silt to fine sand. The novelty of our approach lies in combining geomechanical reasoning with geomorphic analogues to understand histories of ancient glacier-dammed lakes and of the glacial lake outburst floods that are sourced from them.

British Columbia↗

USGS2021

Expansion of intertidal mussel beds following disease-driven reduction of a keystone predator

Disease shapes community composition by removing species with strong interactions. To test whether the absence of keystone predation due to disease produced changes to the species composition of rocky intertidal communities, we leverage a natural experiment involving mass mortality of the keystone predator Pisaster ochraceus from Sea Star Wasting Syndrome. Over four years, we measured dimensions of mussel beds, sizes of Mytilus californianus, mussel recruitment, and species composition on vertical rock walls at six rocky intertidal sites on the central California coast. We also assessed the relationship between changes in mussel cover and changes in sea star density across 33 sites along the North American Pacific coast using data from long-term monitoring. After four years, the lower boundary of the central California mussel beds shifted downward toward the water 18.7 ± 15.8 cm (SD) on the rock and 11.7 ± 11.0 cm in elevation, while the upper boundary remained unchanged. In central California, downward expansion and total area of the mussel bed were positively correlated with mussel recruitment but were not correlated with pre-disease sea star density or biomass. At a multi-region scale, changes in mussel percent cover were positively correlated with pre-disease sea star densities but not change in densities. Species composition of primary substrate holders and epibionts below the mussel bed remained similar across years. Extirpation of the community below the bed did not occur. Instead, this community became limited to a smaller spatial extent while the mussel bed expanded.

California, Oregon↗

USGS2021

It’s complicated…environmental DNA as a predictor of trout and char abundance in streams

The potential to provide inferences about fish abundance from environmental (e)DNA samples has generated great interest. However, the accuracy of these abundance estimates is often low and variable across species and space. A plausible refinement is the use of common aquatic habitat monitoring data to account for attributes that influence eDNA dynamics. We therefore evaluated the relationships between eDNA concentration and abundance of bull trout ( Salvelinus confluentus ), westslope cutthroat trout ( Oncorhynchus clarkii lewisi ) and rainbow trout ( Oncorhynchus mykiss ) at 42 stream sites in the Intermountain West (USA and Canada) and tested whether accounting for site-specific habitat attributes improved the accuracy of fish abundance estimates. eDNA concentrations were positively associated with fish abundance, but these relationships varied by species and site, and there was still considerable variation unaccounted for. Random site-level differences explained much of this variation, but specific habitat attributes of those sites explained relatively small amounts of this variation. Our results underscore that either eDNA sampling or environmental characterization will require further refinement before eDNA can be used reliably to estimate fish abundance in streams.

Idaho, Montana, Oregon, Washington↗

USGS2020

Observations on the May 2019 Joffre Peak landslides, British Columbia

Two catastrophic landslides occurred in quick succession on 13 and 16 May 2019, from the north face of Joffre Peak, Cerise Creek, southern Coast Mountains, British Columbia. With headscarps at 2560 m and 2690 m elevation, both began as rock avalanches, rapidly transforming into debris flows along middle Cerise Creek, and finally into debris floods affecting the fan. Beyond the fan margin, a flood surge on Cayoosh Creek reached bankfull and attenuated rapidly downstream; only fine sediment reached Duffey Lake. The toe of the main debris flow deposit reached 4 km from the headscarp, with a travel angle of 0.28; while the debris flood phase reached the fan margin 5.9 km downstream, with a travel angle of 0.22. Photogrammetry indicates the source volume of each event is 2-3 Mm3, with combined volume of 5 Mm3. Lidar differencing, used to assess deposit volume, yielded a similar total result; although error in the depth estimate introduced large error and masks expected increase due to dilation and entrainment. The average velocity of the rock avalanche-debris flow phases, from seismic analysis, was ~25-30 m/s, and the velocity of the 16 May debris flood on the upper fan, from super-elevation and boulder sizes, was 5-10 m/s. The volume of debris deposited on the fan was ~104 m3, 2-orders of magnitude less than the avalanche/debris flow phases. The 13 May landslide was apparently triggered by rapid snowmelt; with debuttressing triggering the 16 May event. While spring 2019 was warm, it wasn’t unusual. It is likely that progressive glacier retreat and permafrost degradation were the conditioning factors; precursor activity was noted at least 1 yr previous; thus, the mountain was primed to fail and average seasonal snowmelt tipped the balance.

British Columbia↗

USGS2020

Cryptic and extensive hybridization between ancient lineages of American crows

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow ( Corvus caurinus ) and American Crow ( Corvus brachyrhynchos ) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%‐divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late‐generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human‐modified habitats.

Washington↗

USGS2019

Looking for love under the ice: Using passive acoustics to detect burbot (Lota lota: Gadidae) spawning activity

Burbot ( Lota lota : Gadidae) is a difficult species to manage effectively due to its preference for deep‐water habitats and under‐ice spawning behaviour, resulting in a poor understanding of its reproductive activity. However, the use of acoustic signalling by burbot as part of their mating system has recently been described and this behaviour may provide a means of investigating questions regarding the spatial and temporal distribution of spawning aggregations using passive acoustic monitoring. We used audio and video recording to confirm that burbot vocalise and that these vocalisations can be detected under field conditions as well as to characterise the relationship between burbot acoustic signalling and spawning behaviour. We also evaluated the feasibility of locating and monitoring burbot spawning aggregations in real time using passive acoustics. Burbot vocalisations were difficult to identify with only about 6% of the recordings containing calls being successfully identified as such in the field. Burbot vocalised more often between sundown and sunrise than during daylight hours. Calls recorded at night tended to be lower frequency, longer duration, and have lower bandwidth than those made during the day. Burbot vocalisations could not be recorded in conjunction with video recordings of spawning activity, indicating that burbot may not call during active spawning, but may use acoustic communication to signal the onset of reproductive readiness and to form pre‐spawning aggregations. While burbot calls were readily identifiable, observers had a difficult time identifying burbot calls in real time under field conditions. Passive acoustic monitoring demonstrates considerable potential as a management tool to locate burbot spawning grounds and identify periods of activity, but may not be an appropriate technique for monitoring spawning activity in real time.

British Columbia↗

USGS2019

Two-dimensional seismic velocities and structural variations at three British Columbia Hydro and Power Authority (BC Hydro) dam sites, Vancouver Island, British Columbia, Canada

Summary In June, 2017, we acquired seismic data along five linear profiles at three British Columbia Hydro and Power Authority (BC Hydro, a Canadian provincial Crown Corporation) dam sites (John Hart, Ladore, and Strathcona Dams) on Vancouver Island, British Columbia, Canada. We also attempted to acquire linear seismic profiles at two additional BC Hydro dam sites (Ruskin Dam and Stave Falls Dam) east of the City of Vancouver, British Columbia, Canada; however, due to a seismograph programming error, little active-source data from Ruskin Dam and Stave Falls Dam were recorded. Thus, results from Ruskin Dam and Stave Falls Dam are not included in this report. At the three dam sites with successful data acquisition, we acquired both active- and passive-source data. Data acquisition details for each of the three dam sites varied in terms of seismic sources, the number of seismographs, and profile length and orientation. However, for active-source acquisition at each dam site, we acquired one or more linear seismic profiles ranging in length from about 150 to 400 meters (m), and along each profile, seismograph spacing was either 3 m or 5 m (see appendix 1). All data were recorded in three components (vertical and two horizontals). To greatly increase the resolution of the seismic velocity structure along these profiles, we co-located active sources at each seismograph.

British Columbia↗

USGS2018

Holocene surface rupture history of an active forearc fault redefines seismic hazard in Southwestern British Columbia, Canada

Characterizing the hazard associated with Quaternary‐active faults in the forearc crust of the northern Cascadia subduction zone has proven challenging due to historically low rates of seismicity, late Quaternary glacial scouring, and dense vegetation that often obscures fault‐related geomorphic features. We couple lidar topography with paleoseismic trenching across the Leech River Fault on southern Vancouver Island to produce the first detailed surface rupture history of an onland forearc fault in British Columbia, Canada. The results indicate that this fault produced three surface‐rupturing earthquakes in the last ∼9 kyr and is therefore capable of producing large (Mw>6) earthquakes in the future. We provide new constraints on the fault's length (∼130 km) and Holocene slip rate (≥0.2–0.3 mm/year) that, together with the earthquake ages, should be incorporated into new seismic hazard assessments and building code practices relevant to urban centers in southwestern British Columbia (Canada) and northwestern Washington State (United States).

British Columbia↗

USGS2018

Optimizing release strategies: A stepping-stone approach to reintroduction

Evaluation of alternative management strategies enables informed decisions to accelerate species recovery. For reintroductions, post-release survival to reproductive age is a key parameter influencing population growth. Here we trial a ‘stepping-stone’ method to maximize the success of captive-bred animals when the availability of more suitable wild-born release candidates is limited. Our approach makes use of relatively safe and established wild populations to prepare captive-bred individuals for eventual translocation to a final release destination, thus building resilience through establishment of multiple populations over time. We developed a novel multievent model integrating encounter history and biotelemetry data to evaluate reintroduction strategies for the critically endangered Vancouver Island marmot (Marmota vancouverensis). We compared post-release survival of 176 individuals (52 wild-born, 47 captive-bred marmots released directly to destinations, and 77 captive-bred marmots released with a stepping-stone approach). Survival estimates to prime breeding-age (PBA), were then used to quantify expected success of potential release strategies. Our analysis indicates that post-release survival varies by source population and release method, as well as age, season, year, and years since release. Conditional on an objective of maximizing survival to PBA, our results suggest that using wild-born marmots for translocations as often as possible, and stepping-stone captive-bred marmots prior to final release, will result in the best outcomes. There was a 0.86 probability that survival to PBA was greater for captive-bred marmots released as yearlings using a stepping-stone approach (survival to PBA mode = 0.13, 95% CRI = 0.05-0.30) than for captive-bred animals that were directly released to destination sites as one-year-olds (survival to PBA mode = 0.04, 95% CRI = 0.01-0.24). Consequently, the stepping-stone approach yields much higher population establishment or growth potential than previous release strategies that used captive-bred marmots. Optimizing the combination of release candidates, sites, and timing can thereby increase the effectiveness of reintroductions.

British Columbia↗