DRAFT

1 Overview

Currently, there is no place where the public can readily find answers to general questions regarding the state of salmon in BC and the Yukon. It is difficult to find credible answers to questions such as:

  • How are salmon populations doing?
  • What are the major factors influencing the state of salmon?
  • What is being done - and what could be done - to mitigate threats and recover salmon populations and the habitats they depend on?

This lack of publicly accessible information about the state of Pacific salmon and opens the door to misinformation that can hinder and harm our ability to conserve salmon.

The Pacific Salmon Foundation (PSF) is developing an inaugural State of Salmon report, which will address these questions and be tailored to the questions of a salmon-informed public audience. The outputs will be updated regularly (e.g. annually or biennially) to reflect new data and changes to the state of Pacific salmon populations.

As part of the State of Salmon report, PSF’s Salmon Watersheds Program is developing simple metrics to report on the status and trends in regional abundance for each species. This document outlines the data sources and analytical approach used to estimate time series of spawner abundance and total returns (including catch) for each species and region, as well as draft results. This information is intended to be shared with external partners in order to garner feedback on preliminary results that will inform our approach moving forward. Our analyses are ongoing, and the information here is subject to change on a regular basis, but we will endeavor to keep it up-to-date.

2 Methods

2.1 General Approach

We report spawner abundance for each of nine regions that represent all major Pacific salmon-bearing watersheds in Canada: Yukon, Transboundary, Haida Gwaii, Nass, Skeena, Central Coast, Vancouver Island & Mainland Inlets, Fraser, and Columbia. These regions are also used to organize data in the Pacific Salmon Explorer. There are a relatively small number of Pacific salmon that spawn in the MacKenzie River basin in Arctic Canada that are currently not considered here.

We separate five species of Pacific salmon: Chinook, chum, coho, pink, and sockeye. We also report status of steelhead trout for regions where spawner data are are available and can reasonably be expanded to yield an estimate of regional abundance (Nass, Skeena, interior Fraser, and Columbia). When assessing biological status, pink salmon are often separated into even- and odd-year lineages due to their consistent 2-year life cycle. However, for the general overviews provided in State of Salmon we consider generational averages, which take the running average of even- and odd-year lineages for pink salmon. This approach of using generation running averages also smooths over dominant years for sockeye salmon, for which many populations display cyclic dominance. Shifts in dominance between even- and odd-year pink populations or declines in sub-dominant years of sockeye salmon are considered in a more nuanced way in our reporting when discussing how changes in abundance have been reflected in the diversity and distribution of each species within the region.

For each of these species, where data are available, we construct an index of spawner abundance at the regional scale. We focus on spawner abundance because these data are more readily accessible and easily summarised at different spatial scales. Spawners represent the abundance of salmon available to reproduce and contribute to future generations, as well as to meet cultural and ecological needs within watersheds, and thus provides a measure of status relevant to communities and ecosystems. We recognize that commercial catch has historically been a substantial portion of salmon that return to the coast, and that ignoring declines in catch will underestimate the declines in overall salmon abundance. Therefore, we include information on total return (i.e. spawners plus catch) for species and regions with reliable data. This assessment of the data landscape is ongoing, and one of the reasons why the results here should be considered preliminary.

We have strived to base our assessments on the best available data for each region and species. For example, stocks that are governed by international treaties may be monitored by the Pacific Salmon Commission (PSC), and tend to have reliable time series of abundance available at regional scales. These data sources are outlined in the Region-Specific Data section. For species and regions where aggregate abundance is not reported by the PSC or DFO at the scale needed, we adapted our approach to make the best use of available data. In most cases, this meant expanding spawner abundance from stream-level estimates to get a regional scale index of spawner abundance using two types of expansion factors (English et al. 2016). This expansion process is described in the Expansion Factors section.

R code and data for the analyses described here are available on GitHub at https://github.com/salmonwatersheds/state-of-salmon.

2.2 Region-Specific Data Sources

Here we describe specific data sources for the abundance of salmon and steelhead at the regional scale.

2.2.1 Yukon

The Canadian portion of the Yukon River is home to Chinook, chum, and coho salmon. Border escapement and total abundance (i.e., run size) of Canadian-origin Chinook and fall chum salmon in the mainstem Yukon River are counted at the Eagle Sonar station on the Yukon/Alaska border, and are available from the Yukon River Panel in their Joint Technical Committee (JTC) Reports. Specifically, we used Chinook RR Spawning escapement estimate and RR Canadian origin total run size estimate from Appendix B11 of JTC (2024). Chum salmon spawner data were from Appendix B16 and total return from Appendix B20 of JTC (2024).

There are also Chinook, chum, and coho salmon in the Canadian portion of the Porcupine River, which joins the mainstem Yukon River in Alaska. Data on escapement to the Porcupine River are more patchy and not currently included here.

2.2.2 Transboundary

For Chinook, coho, and sockeye, we assessed regional abundance in the Transboundary using estimates of border escapement provided in the PSC’s Joint Transboundary Technical Committee Reports. The regional spawner abundance and total return of Chinook was calculated as the sum of escapement and total return, respectively, to the Stikine, Alsek, and Taku watersheds, available from Table B2 of CTC (2023) and provided to PSF staff (E. Hertz) in Excel format following a data request to PSC’s Chinook Technical Committee. Sockeye spawner abundance and total return was the sum escapement and Terminal Run, respectively, from the Stikine River [Appendix B21; Transboundary Technical Committee (2022)] and the Taku River [Appendix D15; Transboundary Technical Committee (2022)]. Coho spawner abundance and total return are available for the Taku River only, taken from Appendix D20 of Transboundary Technical Committee (2022). The most recent TTC report (Transboundary Technical Committee 2022) includes data through 2021, and updated Appendices with data through 2023 were provided to PSF staff on request to Aaron Foos (DFO).

Pink and chum salmon are less extensively monitored in the Transboundary, with ongoing escapement available only from the Canyon Island fish wheel on the Taku River. We used the index of escapement from the Canyon Island fish wheel as an index of regional spawner abundance for pink and chum (Transboundary Technical Committee 2022). There are some historical data for pink salmon spawner abundance in the Nakina River of the Taku watershed, but this location has not reported data since 1998 and thus we did not include this stream survey in our index of regional abundance.

Steelhead trout have also been enumerated at the Canyon Island fish wheel since 1987, though recent estimates have been patchy and the timing of the fish wheel likely misses a substantial portion of the Taku summer run steelhead. Thus, we do not use the Canyon Island counts as an index of steelhead abundance.

We recognize that our approach in the Transboundary lacks information from many un-monitored watersheds, in particular the smaller watersheds of the Chilkat, Unuk, and Whiting Rivers. In the absence of better monitoring, we choose to report available data as a proxy for regional abundance, and note that the contribution of these smaller watersheds to total regional abundance for each species is likely small.

2.2.3 Haida Gwaii

For Haida Gwaii, we followed the Expansion Factors approach described below, with the exception of Chinook salmon for which the only data are from a single enumeration project on the Yakoun River. We show Yakoun River spawner abundance up to 2006, but no estimates were available for the most recent generation making the current status for Haida Gwaii Chinook “Unknown”. For other species, we compared our expansion approach to estimates of escapement to Area 1 reported in the TCNB (2023) (Appendix 30) but found the Area 1 estimates were much lower (except for Chinook, for which Appendix 30 matched the Yakoun River estimates exactly).

2.2.4 Nass

Chinook spawners and total return were taken from the PSC’s Chinook Technical Committee (CTC) data sets and are reported on in CTC (2024). Specifically, Nass River (Area 3) escapement and terminal run are provided in Table B3 (fields Esc and t.run respectively).

Chum and pink salmon were expanded from available spawner surveys, as described below in Expansion Factors. Note that the designation of indicator streams for the Nass region was based on English et al. (2018) rather than the indicator designation in NuSEDS.

Nass coho escapement and total return were reconstructed by English et al. (2023) for each of the three coho Conservation Units in the Nass region for 1992-2022. We summed the total escapement (TE) and total return (Total Run) across these three Conservation Units to yield the regional aggregate abundance.

Nass sockeye escapement and total return were derived from the Northern Boundary Sockeye Run Reconstruction (NBSRR) Model that estimates total escapement and total return for the Skeena and Nass Rivers (English et al. 2004; English et al. 2017). Specifically, we used fields TE (Total Escapement) and Total run for Area 3 sockeye from the 2022 update to the North and Central Coast (NCC) Salmon Database Version 2, maintained by LGL Ltd. with support from PSF (English et al. 2016; Challenger et al. 2018; English et al. 2018).

The index of Nass steelhead spawner abundance and total return are for the Nass Summer CU, developed in collaboration with the Nisga’a Fish and Wildlife Department and LGL Ltd. (English et al. 2023). There is another steelhead CU in the Nass region - Nass Winter - which is not well monitored and therefore not included in our index of spawner abundance.

2.2.5 Skeena

Chinook spawners were taken from the PSC’s Chinook Technical Committee (CTC) data sets and are reported on in CTC (2024). Specifically, Skeena River (Area 4) escapement is provided in Table B3 (field GSI esc). There is no total return for the Skeena River provided in the CTC data sets.

For Skeena chum salmon, we used escapement to Area 4 reported in Table 32 of TCNB (2023) rather than expanding escapement from spawner surveys because surveys of indicator stream have been increasingly patchy over the last 30 years, making expansions uncertain.

Skeena coho and pink salmon were expanded from available spawner surveys, as described below in Expansion Factors. Note that the designation of indicator streams for the Skeena region was based on English et al. (2018) rather than the indicator designation in NuSEDS.

Skeena sockeye escapement and total return were derived from the Northern Boundary Sockeye Run Reconstruction (NBSRR) Model that estimates total escapement and total return for the Skeena and Nass Rivers (English et al. 2004; English et al. 2017). Specifically, we used fields TE (Total Escapement) and Total run for Area 4 sockeye from the 2022 update to the North and Central Coast (NCC) Salmon Database Version 2, maintained by LGL Ltd. with support from PSF (English et al. 2016; Challenger et al. 2018; English et al. 2018).

The index of Skeena steelhead spawner abundance is derived from estimated escapement of Skeena Summer steelhead at the Tyee Test Fishery (1956 - present), provided by the Province via email. As for other regions, these estimates may not capture winter-run steelhead, for which data are not available.

2.2.6 Central Coast

Estimates of spawner abundance for all salmon species on the Central Coast were expanded from available spawner surveys, as described below in Expansion Factors. Note that the designation of indicator streams for the Central Coast region was based on English et al. (2018) rather than the indicator designation in NuSEDS.

There are two spawner surveys for Steelhead trout on the Central Coast that we report in the Pacific Salmon Explorer. However, these surveys have been patchy through time with the most recent estimates from 2016. Therefore, spawner abundance and total return for Central Coast steelhead was considered unknown.

2.2.7 Vancouver Island & Mainland Inlets

Estimates of spawner abundance for all salmon species and steelhead in Vancouver Island & Mainland Inlets were expanded from available spawner surveys, as described below in Expansion Factors. We began the time series of spawner abundance for this region in 1953 because earlier years had high expansion factor values indicating lack of monitoring on key indicator streams. We are exploring the potential use of well-monitored Chinook indicator stocks to inform an index of spawner abundance and total return for Chinook (CTC 2024).

2.2.8 Fraser

Escapement of Fraser Chinook is provided in the PSC’s Chinook Technical Committee (CTC) data sets and are reported on in CTC (2024). We took the sum of escapement to indicator stocks in Table B6 (all Spring/Summer, Harrison, Lower Shuswap, Nicola, Lower Chilcotin, and Chilko) as an index of spawner abundance. Catch of Fraser Chinook is provided in Table A14 but cannot be directly added to indicator escapement to yield total return. PSF is pursuing Chinook run reconstruction output for the Fraser to inform total return.

Fraser chum catch and escapement are reported by PSC’s Chum Technical Committee (TTCHUM 2023), however the most recently published report only contains estimates for 2010-2019 and the data were not otherwise made available to us. As such, we estimated spawner abundance from available spawner surveys in NuSEDS as described below in Expansion Factors. We note that, although most NuSEDS data was updated to include surveys in 2022, for Fraser chum the most recent data was from 2020.

For interior Fraser coho CUs, spawner abundance and total return were provided by DFO on data request (pers. comm. Marissa Glavas, Data Manager, Fraser River Stock Assessment). We report total return as the sum of “Total Pre-Fishery Abundance” to all CUs and spawner abundance as the sum of “Total Return” (i.e., final spawner estimate plus fish removed from the system by DFO Salmon Enhancement Program or First Nations) to all CUs (we note the term “total return” is used differently in this data set to represent spawners!). We recognize that using these estimates ignores coastal coho populations, for which monitoring has been patchy through time.

Data on total return and spawner abundance of pink salmon and sockeye salmon in the Fraser region are provided by the PSC and accessed through the Fraser Panel Annual Report: Data Application (Pacific Salmon Commission 2024). Total escapement for Fraser sockeye was not yet available for 2023, and 2023 total return is considered preliminary. We note that Fraser River pink salmon are only counted in the dominant, odd-year run.

Steelhead trout in the Fraser are monitored by the Province in at least 10 different streams, but these data are not readily available. Relatively reliable estimates of steelhead spawner abundance are available at the CU-level for interior Fraser steelhead from the Thompson Summer CU (monitored at the Thompson River) and Mid Fraser Summer CU (primarily monitored at the Chilcotin River; data in the Pacific Salmon Explorer). We used the sum of CU-level spawner abundance for these two CUs as an index of Fraser steelhead abundance. We note that this approach does not include more coastal populations, such as the Lower Fraser Summer steelhead monitored in the Coquihalla River or Boundary Bay Winter steelhead, which may not have declined to the same extent over the past decade. However, a lack of publicly accessible data on coastal Fraser steelhead has limited our ability to include these CUs in our index of abundance.

2.2.9 Columbia

We used CU-level estimates of spawner abundance (run reconstructions) sourced from DFO (Ogden, pers. comm.) for Chinook and Stockwell and Hyatt (2003) and subsequent updates for sockeye.

There is no monitoring of steelhead trout in the Canadian portion of the Columbia region, but the Okanagan Nation Alliance does enumerate steelhead in akskwəkwant (Inkaneep Creek) and estimate a Canadian portion of steelhead spawning abundance. These data can be found in associated report (e.g. OBMEP (2022)) and are available in the Pacific Salmon Explorer as the CU-level spawner abundance for the Mid Columbia Summer CU.

2.3 Expansion Factors

For species and regions that lacked reliable data on spawner abundance and total return at the appropriate scale, we estimated regional-scale abundance from stream-level surveys. We started with spawner survey data shown in the Pacific Salmon Explorer. Spawner surveys were each assigned to one of the nine regions we considered based on their geographic location. We note that this is slightly different from how spawner survey data are organized in the Pacific Salmon Explorer, where data are organized by Conservation Units (CUs) that may span regional boundaries (e.g. for pink salmon that have relatively geographically large CUs). In cases of trans-regional CUs, the spawner surveys appear in both regions in the Pacific Salmon Explorer, whereas here we assign spawner surveys to the region in which they fall geographically, regardless of the CU boundary.

Spawner survey data are largely derived from river-level estimates in DFO’s New Salmon Escapement Database System (NuSEDS), but are cleaned up to address issues of, for example, inconsistent naming of streams through time or duplicate data. The spawner survey abundance is equal to the MAX_ESTIMATE in NUSEDS for each year and river population, calculated as the maximum of all fields containing spawner abundance data (e.g. natural adult spawners, natural jack spawners, total broodstock removals). Each of these river populations has been designated as an indicator stream or non-indicator. Indicator streams are observed more consistently in recent decades, tend to have higher spawner abundance, and tend to be monitored using more intensive methods that provide greater accuracy [EnglishEtAl2016NorthCentralCoast]. For further information on the compilation of spawner survey data, see the Pacific Salmon Explorer Technical Report.

Expansion Factor 1, \(F_{1,y/d}\), expands the observed spawner abundances in indicator streams to account for indicator streams that are not monitored in a given year. It is calculated for each year \(y\) of the spawner time series, and relies on a decadal contribution of each indicator stream to the total escapement to all indicator streams, \(P_{d,i}\) in decade \(d\) [EnglishEtAl2016NorthCentralCoast]. The calculation of this decadal contribution requires at least one estimate from each indicator stream for the decade. If a decade does not contain sufficient information (i.e. one or more indicator streams are not monitored at all in a decade), then a reference decade is used to calculate \(P_{d,i}\). This reference decade is chosen to be: (1) the closest decade (historical or future) with sufficient information, or failing (1), (2) the 20-year period from 1980-1999 (Challenger et al. 2018).

For each decade (or reference decade if insufficient information) \(d\), the average number of spawners returning to indicator stream \(i\) is calculated as:

\[\bar{S}_{d,i} = \sum_{y = 1}^{Y_{d,i}} \frac{\hat{S}_{y/d, i}}{Y_{d,i}} \] where \(Y_{d,i}\) is the number of years for which spawner estimates are available within decade \(d\) for stream \(i\). From the average number of spawners for all indicator streams, the decadal proportional contribution of each indicator stream is calculated as:

\[P_{d,i} = \frac{\bar{S}_{d,i}}{\sum_{i=1}^{I} \bar{S}_{d,i}}\] where \(I\) is the total number of indicator streams.

Expansion Factor 1 is then calculated for each year within the decade \(y/d\) based on the decadal contributions and which streams were monitored or not in a given year:

\[F_{1,y/d}=\left( \sum_{i=1}^I P_{d,i} w_{y/d,i} \right)\] where \(w_{y/d,i}\) is 1 if stream \(i\) was monitored in year \(y\) and 0 if stream \(i\) was not monitored in year \(y\). Expansion Factor 1 is then multiplied by the sum of the observed spawners in all indicator stream to yield the expanded estimate of spawner abundances in all indicator streams in the region:

\[S'_{y} = F_{1,y/d} \sum_{i=1}^I \hat{S}_{y,i}\]

Expansion Factor 2 \(F_{2,d}\) expands the spawner abundance to all indicator streams, \(S'_{y}\), to account for non-indcator streams. Unlike Expansion Factor 1, this is calculated for each decade (rather than each year) and then applied to all years within a decade. Like Expansion Factor 1, there needs to be sufficient information within the given decade in order to calculate \(F_{2,d}\), or else a reference decade is chosen. See English et al. (2016) for detailed on how reference decades are chosen in that case.

Expansion Factor 2 is calculated as:

\[F_{2,d} = \frac{\sum_{i = 1}^I \bar{S}_{d,i} + \sum_{j = 1}^{J} \bar{S}_{d,j}}{\sum_{i = 1}^I \bar{S}_{d,i}}\] where \(\bar{S}_{d,i}\) and \(\bar{S}_{d,j}\) are the deacdal average number of spawners in indicator and non-indicator streams, respectively, calculated above. \(J\) is the total number of non-indicator streams. The adjusted total number of spawners in both indicator and non-indicator streams is then calculated as: \[ S''_{y} = F_{2,d} S'_{y} \]

Note that when expanding spawner abundance for spawner-recruit analysis, a third expansion factor is applied to account for streams that are never monitored and for observer (in)efficiency (Peacock et al. 2020). We did not apply this third expansion factor because it is highly undertain and we are interested in relative changes in abundance through time, so we do not require to expand to absolute abundance.

2.4 Smoothing

We smoothed time series of spawner abundance and total return using a right-aligned running geometric mean over the length of a generation. This reduces the influence of dominant years and produces an index that is less sensitive to stochastic interannual variability that is common in salmon population dynamics. The generation length was based on the dominant life-history type for each species in a particular region. The smoothed spawner abundance in year \(y\) given a generation length \(g\) was calculated as: \[ \bar{S_{y}} = \left( \prod_{t = y-g+1}^y {S''_t} \right)^{1/g} \]

We smoothed abundance from the first year to the most recent year of raw abundance data. If there were intermediate years with missing data, the smoothed abundance was still calculated using the available years (i.e. ignoring the missing data, with the exponent \(1/g\) adjusted so that \(g\) reflected the number of years with data in the generation). If all years in a generation were missing data, than the smoothed abundance was NA.

When plotting, we show the smoothed abundance relative to the long-term historical average, so that species that have vastly different abundances within a region can be plotted on the same y-axis for comparison.

2.5 Quantifying change

We summarize the time series of smoothed spawner abundance or total return at the regional scale using three different metrics:

  1. Current status relative to historical average;
  2. Long-term trend over the entire time series; and
  3. Short-term trend over the most recent three generations.

The currents status provides information on how the most recent spawner abundance compares to past values, while the two trend metrics provide information on the average direction of change and are less sensitive to the magnitude of current abundance.

The current status relative to historical average is calculated as: \[ (\bar{S_{y}} - \bf{S}) / \bf{S} \] where \(y\) is the most recent year for which the index of spawner abundance could be calculated and \(\bf{S}\) is the average smoothed spawner abundance over the entire time series.

The long-term trend is calculated following the recommendations of D’Eon-Eggertson et al. (2015), who found that the correct identification of declines in salmon population abundance may be most reliable when considering the entire time series, and applying regression-based estimates of change calculated from log-transformed and smoothed spawner abundances. As such, we fit a simple linear model to the time series of \(\log ( \bar{S}_{y} )\) over \(y\). The resulting slope, \(m\), gives an average estimate of annual change:

\[ m = \log \left( \frac{\bar{S}_{y}}{\bar{S}_{y-1}} \right)\].

With some rearranging, we can calculate the average annual percent change as: \[ \frac{\bar{S}_{y} - \bar{S}_{y-1}}{\bar{S}_{y-1}} = \frac{ e^m \bar{S}_{y-1} - \bar{S}_{y-1}}{\bar{S}_{y-1}} = e^m - 1\]. If the slope \(m\) is not significantly different from zero (\(p \geq 0.05\)), then the trend is classified as “stable” regardless of the magnitude of the estimated slope. If the slope is significantly different from zero (\(p < 0.05\)), then the trend is classified as increasing if \(m > 0\) or decreasing if \(m < 0\).

The short-term trend is calculated as described above, but only using the time series over the most recent three generations.

**Above:** Simulated time series of spawner abundance (grey line), smoothed with a running geometric mean over the generation length (4 years; black line). The three metrics are illustrated: **Current status** relative to historical average (black);  **Long-term trend** over the entire time series (green); and **Short-term trend** over the most recent three generations (clay). The shaded regions of the regression lines indicate the 95% prediction intervals. In this example, the short-term trend is not statistically significant (p > 0.05), and so the trend is characterized as 'stable' despite the negative slope.

Above: Simulated time series of spawner abundance (grey line), smoothed with a running geometric mean over the generation length (4 years; black line). The three metrics are illustrated: Current status relative to historical average (black); Long-term trend over the entire time series (green); and Short-term trend over the most recent three generations (clay). The shaded regions of the regression lines indicate the 95% prediction intervals. In this example, the short-term trend is not statistically significant (p > 0.05), and so the trend is characterized as ‘stable’ despite the negative slope.

2.5.1 Additional Considerations

The standardized approach to quantifying change described above was applied across regions and species. We then considered two additional attributes prior to reporting on the State of Salmon:

  1. If the species and region had no data in the most recent generation, the current status was “Unknown” even if historical abundance information was available. This was the case for Haida Gwaii Chinook, for which there were no reliable estimates of spawners since 2006.

  2. If the time series of abundance (spawners or total return) contained fewer than 20 years of data, we could not reliably establish a historical baseline and the current status was therefore “Unknown”. This applied to Columbia Chinook and Columbia steelhead, both of which have been monitored since 2006.

3 Results

NOTICE!

These are PRELIMINARY results and the following figures should not be taken as a statement on the status of salmon in these regions. At this time, we are working to verify datasets and refine our methods.

  • These trends in aggregated spawner abundance provide a broad overview of how salmon are doing, but must be interpreted with caution due to the diversity of populations with potentially different outlooks within each region and species. For example, Chinook in Vancouver Island & Mainland Inlets have an aggregate spawner abundance that is 95.7% above historical average, but some stocks (e.g., West Coast Vancouver Island) are of conservation concern. Preserving the population diversity within regions is key to the long-term resilience of salmon, particularly in the face of climate change.
  • The status of spawner abundance is, in most cases, more positive than the status in total return, highlighting that declines in fisheries have been pervasive.

3.1 All regions

Highlights

  • Across regions, pink and sockeye salmon seem to be doing the best (spawner abundance above historical average for over half of the regions assessed) and steelhead and chum are doing the worst (spawner abundance below historical average in almost all regions assessed).
  • The Fraser region has one of the most positive outlooks, with spawner abundances of coho and pink above historical average. Although the historically productive fisheries supported by this Region have been decimated, total return for coho and pink are now above below historical average, and short-term trends in both spawner abundance and total return for these species are stable, providing hope that declines can be stemmed and populations rebuilt.
  • Further north, in the Yukon and Transboundary, almost all species (with the exception of Transboundary sockeye salmon) are well below historical average with recent declines in either spawners, total return, or both. This finding is concerning, given that climate change may render southern habitats less suitable for salmon into the future and more northern populations must be nurtured as potential strongholds. Management of salmon harvest in these northern regions is complex, and involves international treaties that are not easily negotiated or adapted in response to environmental changes or emerging trends, challenging efforts to stem declines and recover these populations.

Below: Summary of current status of both spawner abundance and total return across regions and species. Regions and species below historical average are in red and above historical average are in green. Black cells indicate species that are not thought to be present in the given region, and grey cells indicate status was Unknown.

3.2 Yukon

Highlights

  • Chinook and chum salmon spawners have both experienced precipitous declines in recent years.
  • Trends in Chinook total return show statistically significant declines of 3% and 4% per year over long and short time frames, respectively, highlighting the significant, persistent decline that has led to closure of fisheries, including First Nations subsistence fishing, throughout in the Canadian portion of the Yukon.
  • Accompanying declines in abundance, there is also evidence that Yukon Chinook body size and reproductive capacity have been reduced (Ohlberger et al. 2020).
  • These shifts in quality of escapement, combined with warming rivers during spawning migrations and shifting ocean conditions with climate change, may continue to challenge recovery efforts for Yukon salmon (Schoen et al. 2023).
**Above:** Spawner abundance through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Spawner abundance through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

**Above:** Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Below: Summary of the three metrics for spawners and total return. For long-term trends and short-term trends, the percentage is the average annual change calculated over the respective time period. Grey values indicate the % change was not significantly different from zero (i.e., stable), red values indicate declining trends, and green values indicate increasing trends.

3.3 Transboundary

Highlights

  • Most species in the Transboundary are below historical average spawner abundance, except sockeye salmon which have shown a recent increase, driven by increases in river-type type sockeye in the Taku River.
  • As in many other regions, the most recent index of spawner abundance is lower than the index of total return for the three Transboundary species that have both types of data (Chinook: -49.5% compared to -56%, coho: -13.9% compared to -22.5%, and sockeye: 34.1% compared to -11.4%), indicating that declines in catch have also occurred in recent years.
**Above:** Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

**Above:** Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Below: Summary of the three metrics for spawners and total return. For long-term trends and short-term trends, the percentage is the average annual change calculated over the respective time period. Grey values indicate the % change was not significantly different from zero (i.e., stable), red values indicate declining trends, and green values indicate increasing trends.

3.4 Haida Gwaii

Highlights

  • All species for which we had reliable data are currently below their historical average, with long-term trends either stable or declining.
  • Unlike in other regions, where pink salmon have been doing relatively well, pink salmon have shown the most dramatic declines, and most recent spawner abundance was 77% lower than the historical average.
  • There is only one indicator stream for Chinook in Haida Gwaii,the Yakoun River, which has not been monitored since 2006 and although abundance was relatively high at that time, current status is unknown.
**Above:** Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

**Above:** Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Below: Summary of the three metrics for spawners and total return. For long-term trends and short-term trends, the percentage is the average annual change calculated over the respective time period. Grey values indicate the % change was not significantly different from zero (i.e., stable), red values indicate declining trends, and green values indicate increasing trends.

3.5 Nass

Highlights

  • While most species show declines over both long- and short-time frames, Nass pink and coho salmon have rebounded to above historical average abundance in the most recent generation with positive long-term trends.
  • Chum salmon are the species furthest below their historical average (-39.5%) with an average long-term decline of -0.7%.
  • Although sockeye total return is 40.4% below historical average with concerning negative short- and long-term trends, increases over the last two years have put spawner abundance above historical average with an increasing short-term trend over the most recent three generations.
**Above:** Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

**Above:** Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Below: Summary of the three metrics for spawners and total return. For long-term trends and short-term trends, the percentage is the average annual change calculated over the respective time period. Grey values indicate the % change was not significantly different from zero (i.e., stable), red values indicate declining trends, and green values indicate increasing trends.

3.6 Skeena

Highlights

  • Most species are below historical average spawner abundance, with chum salmon having declined by 89.2%. Exceptons include sockeye spawner abundance, which has recently increased to 30.5% above historical average in 2022 and pink salmon which were 79.7% in 2022.
  • Total return was only available for Skeena sockeye, and shows increases in the 1980s and 90s, with a decline from the peak in 1996. Slight increases over the last two years put current total return 11% below historical average, with stable long-term and short-term trends. Other research has found that the decline in sockeye abundance from the 1980s has been accompanied by erosion of diversity, with enhanced stocks increasingly dominating returns (Price et al. 2021).
  • Estimates of regional spawner abundance have become increasingly uncertain since 2000, especially for coho, due to a decline in monitoring of indicator streams.
**Above:** Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

**Above:** Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Below: Summary of the three metrics for spawners and total return. For long-term trends and short-term trends, the percentage is the average annual change calculated over the respective time period. Grey values indicate the % change was not significantly different from zero (i.e., stable), red values indicate declining trends, and green values indicate increasing trends.

3.7 Central Coast

Highlights

  • All species are below historical average abundance, with dramatic declines of culturally and economically important sockeye, Chinook, and chum salmon.
  • Pink salmon, which tend to fluctuate in abundance more than other species, have increased only recently from an all-time low in 2019, more than doubling between 2019 and 2021.
  • Expansion factors for regional-scale reconstructions have increased substantially since the mid-2000s, with a decline in monitoring creating uncertainties in our results for all species, but especially coho, sockeye, and Chinook.
**Above:** Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

**Above:** Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Below: Summary of the three metrics for spawners and total return. For long-term trends and short-term trends, the percentage is the average annual change calculated over the respective time period. Grey values indicate the % change was not significantly different from zero (i.e., stable), red values indicate declining trends, and green values indicate increasing trends.

3.8 Vancouver Island & Mainland Inlets

Highlights

  • Most species are currently below historical average, with the exception of Chinook and pink salmon.
  • Aggregate Chinook abundance is 95.7% above historical average, but some populations have concerning declines (e.g., West Coast Vancouver Island Chinook are at low enough abundance to have triggered a Recovery Potential Assessment in recent years) while others are increasing - including Cowichan River which has seen consistent increases in spawner abundance and total return since 2010 (Atlas et al. 2023).
  • Coho were below average and declining - poor outcomes under all 3 metrics.
  • Short-term trends are stable or increasing for Chinook, pink, and sockeye, suggestion possibility for recovery.
  • Dynamics through time are more stable for the regional aggregate because there are many component populations in this large area, but some local trends are concerning highlighting a need to understand and conserve the diversity represented by component populations (for more detailed information on populations within regions, see the Pacific Salmon Explorer).
**Above:** Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

**Above:** Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Below: Summary of the three metrics for spawners and total return. For long-term trends and short-term trends, the percentage is the average annual change calculated over the respective time period. Grey values indicate the % change was not significantly different from zero (i.e., stable), red values indicate declining trends, and green values indicate increasing trends.

3.9 Fraser

Highlights

  • Despite the uptick in Fraser sockeye spawner abundance in the early 2010s, the most recent spawner abundance and total return are below the devastating returns in 2009 that triggered a federal inquiry. The in-season estimate for Fraser sockeye total return in 2023 was 61.2% below historical average.

  • Chinook spawner abundance were below historical average with recent short-term declines. However some life-history types are doing better than others (Atlas et al. 2023), highlighting the need to conserve the within-species diversity that makes salmon resilient to environmental change.

  • Spawner abundance of coho and pink salmon are above historical average, with total return slightly above average.

  • Declines of Chinook, chum, pink, sockeye, and steelhead since the early 2000s may reflect the increasing pressures of climate change on these southern and interior populations.

**Above:** Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

**Above:** Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Below: Summary of the three metrics for spawners and total return. For long-term trends and short-term trends, the percentage is the average annual change calculated over the respective time period. Grey values indicate the % change was not significantly different from zero (i.e., stable), red values indicate declining trends, and green values indicate increasing trends.

3.10 Columbia

Highlights

  • The current status of Columbia Chinook and steelhead is categorized as Unknown because the are insufficient data to establish a meaningful baseline (fewer than 20 years; Chinook 2006-2022 and steelhead 2006-2021).
  • Columbia Chinook and sockeye show promise of rebuilding, with positive long-term trends.
  • The time series for Chinook and steelhead are short, making it difficult to establish a meaningful baseline. For this reason, we show the current status relative to the historical average as “Unknown” for Columbia Chinook and steelhead.
  • Although Chinook spawner abundance has increased from a low in 2011, the absolute abundance is still less than 1000 spawners and the Okanagan population of Chinook was listed as Endangered in 2017 by COSEWIC (COSEWIC 2017). Recent increases should be cautiously celebrated, knowing that this population is still at risk of extirpation.
  • Columbia steelhead have declined precipitously since a peak spawner abundance in 2012.
**Above:** Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Spawner abundance through time as a percentage of the long-term average. Spawner abundance is smoothed using a one-generation geometric running average (right-aligned), and plotted relative to the historical average for each species. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

**Above:** Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Above: Total return through time as a percentage of the long-term average. Closed points highlight the most recent index of spawner abundance. The percent shown is the percent change from the historical average to the most recent year.

Below: Summary of the three metrics for spawners and total return. For long-term trends and short-term trends, the percentage is the average annual change calculated over the respective time period. Grey values indicate the % change was not significantly different from zero (i.e., stable), red values indicate declining trends, and green values indicate increasing trends.

4 References

Atlas, W.I., Sloat, M.R., Satterthwaite, W.H., Buehrens, T., Parken, C.K., Moore, J.W., Mantua, N., Hart, J., and Potapova, A. 2023. Trends in Chinook salmon spawner abundance and total run size highlight linkages between life history, geography and decline. Fish and Fisheries 24(4): 595–617. doi:10.1111/faf.12750.
Challenger, W., Mochizuki, T., English, K., and Bychkov, Y. 2018. North and Central Coast Salmon Database and Analysis System User Manual. LGL Limited for Pacific Salmon Foundation, Sidney, BC.
COSEWIC. 2017. Chinook salmon (Oncorhynchus tshawytscha) Okanagan population: COSEWIC Assessment and Status Report 2017. Available from https://www.canada.ca/en/environment-climate-change/services/species-risk-public-registry/cosewic-assessments-status-reports/chinook-salmon-okanagan-population-2017.html [accessed 4 December 2023].
CTC. 2023. Annual Report of Catch and Escapement for 2022. Pacific Salmon Commission Joint Chinook Technical Committtee Report TCCHINOOK (23)-02: xiii + 202 p. Available from https://www.psc.org/publications/technical-reports/technical-committee-reports/chinook/.
CTC. 2024. Annual Report of Catch and Escapement for 2023. Pacific Salmon Commission Joint Chinook Technical Committtee Report TCCHINOOK (24)-01: xvii + 173 p. Available from https://www.psc.org/reports/tcchinook-24-01.
D’Eon-Eggertson, F., Dulvy, N.K., and Peterman, R.M. 2015. Reliable Identification of Declining Populations in an Uncertain World. Conservation Letters 8(2): 86–96. doi:10.1111/conl.12123.
English, K.K., Alexander, R.F., Ian A. Beveridge, Challenger, W., Percival, N., Hertz, E., and Atkinson, C. 2023. Preliminary Area 3 Salmon (2018–2022) and Nass River Summer-Run Steelhead (1994–2022) Escapement, Catch, Run Size, and Exploitation Rate Estimates. Prepared by LGL Ltd. and Nisga’a Lisims Government - Fisheries & Wildlife Department for the Pacific Salmon Foundation and Nisga’a-Canada-BC Treaty’s Joint Fisheries Management Committee. Available from https://salmonwatersheds.ca/document/lib_526/.
English, K.K., Gazey, W.J., Peacock, D., and Oliver, G. 2004. Assessment of the Canadian and Alaskan Sockeye Stocks Harvested in the Northern Boundary Fisheries using Run Reconstruction Techniques, 1982-2001. Available from https://salmonwatersheds.ca/document/lib_a_787_06/.
English, K.K., Noble, C., and Carr-Harris, C. 2017. Skeena Sockeye In-river Run Reconstruction Analysis Model and Analysis Results for 1982-2017. Report prepared for Pacific Salmon Foundation by LGL Ltd. and Fisheries and Oceans Canada. Available from https://salmonwatersheds.ca/document/lib_452/ [accessed 23 November 2023].
English, K.K., Peacock, D., Challenger, W., and Mochizuki, T. 2016. North and Central Coast Salmon Escapement, Catch, Run Size and Exploitation Rate Estimates for each Salmon Conservation Unit for 1954-2014. (November): 102.
English, K.K., Peacock, D., Challenger, W., Noble, C., Beveridge, I., Robichaud, D., Beach, K., Hertz, E., and Connors, K. 2018. North and Central Coast Salmon Escapement, Catch, Run Size and Exploitation Rate Estimates for each Salmon Conservation Unit for 1954-2017. Prepared by LGL Ltd. for the Pacific Salmon Foundation. 125p. Available from https://salmonwatersheds.ca/document/lib_451/.
JTC. 2024. Yukon River Salmon 2023 Season Summary and 2024 Season Outlook. Yukon River Joint Technical Committee Report YUKON JTC (24)-01: vi + 181 p. Available from https://www.yukonriverpanel.com/publications/yukon-river-joint-technical-committee-reports/.
OBMEP. 2022. 2021 Okanogan Subbasin Steelhead Spawning Abundance and Distribution. Okanogan Basin Monitoring and Evaluation Program. Colville Confederated Tribes Fish and Wildlife Department. Submitted to the Bonneville Power Administration., Nespelem, WA.
Ohlberger, J., Schindler, D.E., Brown, R.J., Harding, J.M.S., Adkison, M.D., Munro, A.R., Horstmann, L., and Spaeder, J. 2020. The reproductive value of large females: Consequences of shifts in demographic structure for population reproductive potential in Chinook Salmon. Canadian Journal of Fisheries and Aquatic Sciences 77(8): 1292–1301. doi:10.1139/cjfas-2020-0012.
Pacific Salmon Commission. 2024. PSC Annual Report Application. Release 2024.04 [Updated 2024 apr 03; accessed 2024 Apr 05]. Available from https://www.psc1.shinyapps.io/PSC_Annual_Fraser/.
Peacock, S.J., Hertz, E., Holt, C.A., Connors, B., Freshwater, C., and Connors, K. 2020. Evaluating the consequences of common assumptions in run reconstructions on Pacific salmon biological status assessments. Canadian Journal of Fisheries and Aquatic Sciences 77(12): 1904–1920. Canadian Science Publishing. doi:10.1139/CJFAS-2019-0432/ASSET/IMAGES/CJFAS-2019-0432IEQ43.GIF.
Price, M.H.H., Moore, J.W., Connors, B.M., Wilson, K.L., and Reynolds, J.D. 2021. Portfolio simplification arising from a century of change in salmon population diversity and artificial production. Journal of Applied Ecology 58(7): 1477–1486. doi:10.1111/1365-2664.13835.
Schoen, E.R., Howard, K.G., Murphy, J.M., Schindler, D.E., Westley, P.A.H., and von Biela, V.R. 2023. NOAA Arctic Report Card 2023: Divergent Responses of Western Alaska Salmon to a Changing Climate. [object Object]. doi:10.25923/F2HV-5581.
Stockwell, M.M., and Hyatt, K.D. 2003. A summary of Okanagan sockeye salmon (Oncorhynchus nerka) escapement surveys by date and river segment from 1947 to 2001. Can. Data. Rep. Fish. Aquat. Sci. 1106: 34 p. Available from https://salmonwatersheds.ca/library/lib_508/.
TCNB. 2023. U.S./Canada Northern Boundary Area 2022 Salmon Fisheries Management Report and 2023 Preliminary Expectations. Pacific Salmon Commission Joint Northern Boundary Technical Committee Report TCNB (23)-01: ii + 73 p. Available from https://www.psc.org/publications/technical-reports/technical-committee-reports/northern-boundary/.
Transboundary Technical Committee. 2022. Final Estimates of Transboundary River Salmon Production, Harvest and Escapement and a Review of Joint Enhancement Activities in 2020. Pacific Salmon Commission Transboundary Technical Committee Report TCTR (22)-01: xii + 162 p. Available from https://www.psc.org/publications/technical-reports/technical-committee-reports/transboundary/.
TTCHUM. 2023. 2019 Post-Season Summary Report. Pacific Salmon Commission Joint Chum Technical Committee 3-11: viii + 60. Available from https://www.psc.org/publications/technical-reports/technical-committee-reports/chum/.