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.

Jeux de données disponibles au téléchargement

Informations supplémentaires

Champ Valeur
Dernière mise à jour octobre 20, 2025, 02:40 (TU)
Créé octobre 20, 2025, 02:40 (TU)
Domaine / Sujet
Domaine ou sujet du jeu de données catalogué.
Oceans
Titre
Titre du jeu de données.
Groundfish biodiversity change in northeastern Pacific waters under projected warming and deoxygenation
Description
Une description du jeu de données.

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.
Étiquettes / Mots-clés
Mots-clés/étiquettes catégorisant le jeu de données.
Format (CSV, XLS, TXT, PDF, etc.)
Format de fichier du jeu de données.
Taille du jeu de données
Taille du jeu de données en mégaoctets.
Identifiant des métadonnées
Identifiant des métadonnées — peut être utilisé comme identifiant unique pour l’entrée du catalogue.
Date de publication
Date de publication du jeu de données.
2023-11-24
Période temporelle des données (date de début)
Date de début des données dans le jeu de données.
Période temporelle des données (date de fin)
Date de fin des données temporelles du jeu de données.
Zone géospatiale couverte
Région spatiale ou lieu nommé couvert par le jeu de données.
Champ Valeur
Catégorie d’accès
Type d’accès accordé pour le jeu de données (ouvert, restreint, service, etc.).
Licence
Licence utilisée pour accéder au jeu de données.
Open Government Licence - Canada
Restrictions d’utilisation
Restrictions d’utilisation des données.
Emplacement
Emplacement du jeu de données.
https://open.canada.ca/data/en/dataset/661f58cc-e1c3-4249-ba62-08de5b464fe1
Service de données
Service de données pour accéder à un jeu de données.
Propriétaire
Propriétaire du jeu de données.
Fisheries and Oceans Canada | Pêches et Océans Canada
Point de contact
Qui contacter concernant l’accès ?
Government of Canada; Fisheries and Oceans Canada; Pacific Science/Ecosystem Science Division/Marine Spatial Ecology & Analysis Section, 604-999-3490, [email protected]
Courriel du point de contact
Courriel à utiliser pour l’accès ?
Éditeur
Éditeur du jeu de données.
Courriel de l’éditeur
Adresse courriel de l’éditeur.
[email protected]
Auteur
Auteur du jeu de données.
Courriel de l'auteur
Adresse courriel de l’auteur.
Date d’accès
Date à laquelle les données et métadonnées ont été consultées.
Champ Valeur
Identifiant
Identifiant unique du jeu de données.
Langue
Langue(s) du jeu de données
Lien vers la description du jeu de données
Une URL vers un document externe décrivant le jeu de données.
Identifiant persistant
Les données sont identifiées par un identifiant persistant.
Identifiant globalement unique
Les données sont identifiées par un identifiant persistant et globalement unique.
Contient des données sur des individus
Le jeu de données contient-il des données sur des individus ?
Contient des données identifiables sur des individus
Le jeu de données contient-il des données identifiables sur des individus ?
Contient des données autochtones
Le jeu de données contient-il des données sur des communautés autochtones ?
Type de portail
Type de plateforme du portail source.
Champ Valeur
Version
Version du jeu de données
None
Source
Source du jeu de données.
None
Notes de version
Notes de version à propos du jeu de données.
Est une version d’un autre jeu de données
Lien vers le jeu de données dont celui-ci est une version.
Autres versions
Lien vers les jeux de données qui en sont des versions.
Texte de provenance
Texte de provenance des données.
URL de provenance
URL de provenance des données.
Résolution temporelle
Décrit la granularité des données temporelles dans le jeu de données.
Résolution géospatiale en mètres
Décrit la granularité (en mètres) des données géospatiales dans le jeu de données.
Résolution géospatiale (par régions)
Décrit la granularité (par régions) des données géospatiales dans le jeu de données.
Champ Valeur
Permission de la communauté autochtone
Qui détient la permission de la communauté autochtone. Qui contacter pour l’accès à un jeu de données concernant des communautés autochtones.
Permission communautaire
Permission communautaire (qui a accordé la permission).
Les communautés autochtones concernées par le jeu de données
Communautés autochtones dont sont issues les données.
Champ Valeur
Nombre de lignes de données
Pour un jeu de données tabulaire, nombre total de lignes.
Nombre de colonnes de données
Pour un jeu de données tabulaire, nombre total de colonnes uniques.
Nombre de cellules de données
Pour un jeu de données tabulaire, nombre total de cellules contenant des données.
Nombre de relations de données
Pour un jeu de données RDF, nombre total de triplets.
Nombre d’entités
Pour un jeu de données RDF, nombre total d’entités.
Nombre de propriétés de données
Pour un jeu de données RDF, nombre total de propriétés uniques utilisées par les triplets.
Qualité des données
Décrit la qualité des données du jeu de données.
Métrique de qualité des données
Une métrique utilisée pour mesurer la qualité des données, telle que les valeurs manquantes ou les formats invalides.

0 Commentaires

Please login or register to comment.