Groundfish biodiversity change in northeastern Pacific waters under projected warming and deoxygenation

Description:

In the coming decades, warming and deoxygenation of marine waters are anticipated to result in shifts in the distribution and abundance of fishes, with consequences for the diversity and composition of fish communities. Here, we combine fisheries-independent trawl survey data spanning the west coast of the USA and Canada with high-resolution regional ocean models to make projections of how 34 groundfish species will be impacted by changes in temperature and oxygen in British Columbia (BC) and Washington. In this region, species that are projected to decrease in occurrence are roughly balanced by those that are projected to increase, resulting in considerable compositional turnover. Many, but not all, species are projected to shift to deeper depths as conditions warm, but low oxygen will limit how deep they can go. Thus, biodiversity will likely decrease in the shallowest waters (less than 100 m), where warming will be greatest, increase at mid-depths (100–600 m) as shallow species shift deeper, and decrease at depths where oxygen is limited (greater than 600 m). These results highlight the critical importance of accounting for the joint role of temperature, oxygen and depth when projecting the impacts of climate change on marine biodiversity.

The rasters available in this dataset project the occurrence of each of the 34 groundfish species in a 3 km^2 grid cell for the historical baseline, as well as for two emissions scenarios, from each of the two regional ocean models (BCCM and NEP36). Each projection layer is provided as the mean projected occurrence as well as the lower and upper 95% confidence interval of projected occurrence.

Methods:

Estimated species response curves:

We estimated how the observed distribution of groundfish species is determined by temperature, dissolved oxygen and seafloor depth using data from fisheries-independent scientific research trawls spanning the entire American and Canadian west coast. We included data from 4 surveys (NOAA West Coast, NOAA Alaska, NOAA Bering or DFO Pacific) from 2000 to 2019. For each species, we modelled occurrences in the coastwide trawl dataset using a generalized linear model (GLM) using the sdmTMB package in R v. 4.0.2. The predictors were temperature, log dissolved oxygen, log depth and survey. We included quadratic terms for temperature and log depth to allow species occurrences to peak at intermediate values. We fitted a breakpoint function for log dissolved oxygen to reflect the fact oxygen is a limiting factor. We assessed the forecasting accuracy of the SDM by comparing how well a model fitted to only data from 2000 to 2010 could forecast species’ occurrences in trawls within our focal region for the period of 2011–2019. We assessed all 77 groundfish species that were present in the overall trawl dataset, however the final analysis included only the 34 species for which the models had adequate forecasting ability.

Projecting groundfish biodiversity changes:

We based our groundfish biodiversity change projections on two regional models that downscale climate projections: the British Columbia Continental Margin model (BCCM) and the North-Eastern Pacific Canadian Ocean Ecosystem model (NEP36-CanOE). We used a historical baseline of 1986–2005 and future projected values for 2046–2065 based on RCP 4.5 and 8.5 emissions scenarios. Using the models that we validated in our forecasting accuracy assessment, we projected the occurrence of each species in each 3 km^2 grid cell for the historical baseline, as well as for two emissions scenarios, from each of the two regional ocean models.

Uncertainties:

Source survey data was collected by consistent methods with survey-grade GPS for all years included. Data quality is expected to be high. Modeled data are at 3 km resolution. Outputs are as accurate as source input models and are deemed to be of high quality and accurate based upon the precision of model inputs.

Projecting biodiversity responses to climate change involves considerable uncertainty and our approach allows us to quantify some aspects of this. Of the uncertainty that we could quantify, roughly half was due to uncertainty in our SDMs and the remainder was due to regional ocean model uncertainty or scenario uncertainty. This amount of uncertainty in the SDMs is typical, stemming from the fact that contemporary species distributions are also influenced by other factors that we have not included in our model. In addition, although oxygen demand is understood to vary with temperature, limitations in the implementation of breakpoint models prevented us from estimating a temperature-dependent oxygen breakpoint. However, although somewhat unrealistic, this limitation is unlikely to have greatly increased the uncertainty in our SDMs because low oxygen concentrations occurred almost exclusively at depths where temperature variation and projected change was small.

To reduce uncertainty due to year-to-year variation in climate, our model projections are based on 20-year climatologies with a future period that is far enough ahead to ensure that changes are unambiguously due to greenhouse gases. We have made projections based on two different emissions scenarios, and two different regional ocean models that are both downscaled from the same global model, the second generation Canadian Earth System Model (CanESM2), using different downscaling techniques. While the BCCM model was run inter-annually and then averaged to produce the climatologies, the NEP36 model used atmospheric climatologies with augmented winds to force the ocean model and produce representative climatologies. Comparing these regional projections provides an estimate of the uncertainty across different regional downscaling models and methods. We find that the projected impacts of climate change on the groundfish community are more sensitive to the differences in the regional ocean models than they are to the emissions scenarios used. However, these differences are in magnitude (changes tend to be larger based on NEP36 compared with the BCCM) rather than in direction, with both models resulting in similar overall patterns of biodiversity change and turnover for the groundfish community. Over the 60-year time period (1986–2005 versus 2046–2065) used in our study, our projections suggest that groundfish community changes are similar regardless of the scenario used.

Datasets available for download

Additional Info

Field Value
Last Updated October 20, 2025, 02:40 (UTC)
Created October 20, 2025, 02:40 (UTC)
Domain / Topic
Domain or topic of the dataset being cataloged.
Oceans
Title
Title for the Dataset.
Groundfish biodiversity change in northeastern Pacific waters under projected warming and deoxygenation
Description
A description of the dataset.

Description:

In the coming decades, warming and deoxygenation of marine waters are anticipated to result in shifts in the distribution and abundance of fishes, with consequences for the diversity and composition of fish communities. Here, we combine fisheries-independent trawl survey data spanning the west coast of the USA and Canada with high-resolution regional ocean models to make projections of how 34 groundfish species will be impacted by changes in temperature and oxygen in British Columbia (BC) and Washington. In this region, species that are projected to decrease in occurrence are roughly balanced by those that are projected to increase, resulting in considerable compositional turnover. Many, but not all, species are projected to shift to deeper depths as conditions warm, but low oxygen will limit how deep they can go. Thus, biodiversity will likely decrease in the shallowest waters (less than 100 m), where warming will be greatest, increase at mid-depths (100–600 m) as shallow species shift deeper, and decrease at depths where oxygen is limited (greater than 600 m). These results highlight the critical importance of accounting for the joint role of temperature, oxygen and depth when projecting the impacts of climate change on marine biodiversity.

The rasters available in this dataset project the occurrence of each of the 34 groundfish species in a 3 km^2 grid cell for the historical baseline, as well as for two emissions scenarios, from each of the two regional ocean models (BCCM and NEP36). Each projection layer is provided as the mean projected occurrence as well as the lower and upper 95% confidence interval of projected occurrence.

Methods:

Estimated species response curves:

We estimated how the observed distribution of groundfish species is determined by temperature, dissolved oxygen and seafloor depth using data from fisheries-independent scientific research trawls spanning the entire American and Canadian west coast. We included data from 4 surveys (NOAA West Coast, NOAA Alaska, NOAA Bering or DFO Pacific) from 2000 to 2019. For each species, we modelled occurrences in the coastwide trawl dataset using a generalized linear model (GLM) using the sdmTMB package in R v. 4.0.2. The predictors were temperature, log dissolved oxygen, log depth and survey. We included quadratic terms for temperature and log depth to allow species occurrences to peak at intermediate values. We fitted a breakpoint function for log dissolved oxygen to reflect the fact oxygen is a limiting factor. We assessed the forecasting accuracy of the SDM by comparing how well a model fitted to only data from 2000 to 2010 could forecast species’ occurrences in trawls within our focal region for the period of 2011–2019. We assessed all 77 groundfish species that were present in the overall trawl dataset, however the final analysis included only the 34 species for which the models had adequate forecasting ability.

Projecting groundfish biodiversity changes:

We based our groundfish biodiversity change projections on two regional models that downscale climate projections: the British Columbia Continental Margin model (BCCM) and the North-Eastern Pacific Canadian Ocean Ecosystem model (NEP36-CanOE). We used a historical baseline of 1986–2005 and future projected values for 2046–2065 based on RCP 4.5 and 8.5 emissions scenarios. Using the models that we validated in our forecasting accuracy assessment, we projected the occurrence of each species in each 3 km^2 grid cell for the historical baseline, as well as for two emissions scenarios, from each of the two regional ocean models.

Uncertainties:

Source survey data was collected by consistent methods with survey-grade GPS for all years included. Data quality is expected to be high. Modeled data are at 3 km resolution. Outputs are as accurate as source input models and are deemed to be of high quality and accurate based upon the precision of model inputs.

Projecting biodiversity responses to climate change involves considerable uncertainty and our approach allows us to quantify some aspects of this. Of the uncertainty that we could quantify, roughly half was due to uncertainty in our SDMs and the remainder was due to regional ocean model uncertainty or scenario uncertainty. This amount of uncertainty in the SDMs is typical, stemming from the fact that contemporary species distributions are also influenced by other factors that we have not included in our model. In addition, although oxygen demand is understood to vary with temperature, limitations in the implementation of breakpoint models prevented us from estimating a temperature-dependent oxygen breakpoint. However, although somewhat unrealistic, this limitation is unlikely to have greatly increased the uncertainty in our SDMs because low oxygen concentrations occurred almost exclusively at depths where temperature variation and projected change was small.

To reduce uncertainty due to year-to-year variation in climate, our model projections are based on 20-year climatologies with a future period that is far enough ahead to ensure that changes are unambiguously due to greenhouse gases. We have made projections based on two different emissions scenarios, and two different regional ocean models that are both downscaled from the same global model, the second generation Canadian Earth System Model (CanESM2), using different downscaling techniques. While the BCCM model was run inter-annually and then averaged to produce the climatologies, the NEP36 model used atmospheric climatologies with augmented winds to force the ocean model and produce representative climatologies. Comparing these regional projections provides an estimate of the uncertainty across different regional downscaling models and methods. We find that the projected impacts of climate change on the groundfish community are more sensitive to the differences in the regional ocean models than they are to the emissions scenarios used. However, these differences are in magnitude (changes tend to be larger based on NEP36 compared with the BCCM) rather than in direction, with both models resulting in similar overall patterns of biodiversity change and turnover for the groundfish community. Over the 60-year time period (1986–2005 versus 2046–2065) used in our study, our projections suggest that groundfish community changes are similar regardless of the scenario used.
Tags / Keywords
Keywords/tags categorizing the dataset.
Format (CSV, XLS, TXT, PDF, etc)
File format of the dataset.
Dataset Size
Dataset size in megabytes.
Metadata Identifier
Metadata identifier – can be used as the unique identifier for catalogue entry
Published Date
Published date of the dataset.
2023-11-24
Time Period Data Span (start date)
Start date of the data in the dataset.
Time Period Data Span (end date)
End date of time data in the dataset.
GeoSpatial Area Data Span
A spatial region or named place the dataset covers.
Field Value
Access category
Type of access granted for the dataset (open, closed, service, etc).
License
License used to access the dataset.
Open Government Licence - Canada
Limits on use
Limits on use of data.
Location
Location of the dataset.
https://open.canada.ca/data/en/dataset/661f58cc-e1c3-4249-ba62-08de5b464fe1
Data Service
Data service for accessing a dataset.
Owner
Owner of the dataset.
Fisheries and Oceans Canada | Pêches et Océans Canada
Contact Point
Who to contact regarding access?
Government of Canada; Fisheries and Oceans Canada; Pacific Science/Ecosystem Science Division/Marine Spatial Ecology & Analysis Section, 604-999-3490, [email protected]
Contact Point Email
The email to contact regarding access?
Publisher
Publisher of the dataset.
Publisher Email
Email of the publisher.
[email protected]
Author
Author of the dataset.
Author Email
Email of the author.
Accessed At
Date the data and metadata was accessed.
Field Value
Identifier
Unique identifier for the dataset.
Language
Language(s) of the dataset
Link to dataset description
A URL to an external document describing the dataset.
Persistent Identifier
Data is identified by a persistent identifier.
Globally Unique Identifier
Data is identified by a persistent and globally unique identifier.
Contains data about individuals
Does the data hold data about individuals?
Contains data about identifiable individuals
Does the data hold identifiable data about individual?
Contains Indigenous Data
Does the data hold data about Indigenous communities?
Portal Type
Platform type of the source portal.
Field Value
Version
Version of the datatset
None
Source
Source of the dataset.
None
Version notes
Version notes about the dataset.
Is version of another dataset
Link to dataset that it is a version of.
Other versions
Link to datasets that are versions of it.
Provenance Text
Provenance Text of the data.
Provenance URL
Provenance URL of the data.
Temporal resolution
Describes how granular the date/time data in the dataset is.
GeoSpatial resolution in meters
Describes how granular (in meters) geospatial data is in the dataset.
GeoSpatial resolution (in regions)
Describes how granular (in regions) geospatial data is in the dataset.
Field Value
Indigenous Community Permission
Who holds the Indigenous Community Permission. Who to contact regarding access to a dataset that has data about Indigenous communities.
Community Permission
Community permission (who gave permission).
The Indigenous communities the dataset is about
Indigenous communities from which data is derived.
Field Value
Number of data rows
If tabular dataset, total number of rows.
Number of data columns
If tabular dataset, total number of unique columns.
Number of data cells
If tabular dataset, total number of cells with data.
Number of data relations
If RDF dataset, total number of triples.
Number of entities
If RDF dataset, total number of entities.
Number of data properties
If RDF dataset, total number of unique properties used by the triples.
Data quality
Describes the quality of the data in the dataset.
Metric for data quality
A metric used to measure the quality of the data, such as missing values or invalid formats.

0 Comments

Please login or register to comment.