Since 2005, Fisheries and Oceans Canada has been collecting monitoring data for aquatic invasive species (e.g. https://open.canada.ca/data/en/dataset/8d87f574-0661-40a0-822f-e9eabc35780d, https://open.canada.ca/data/en/dataset/503a957e-7d6b-11e9-aef3-f48c505b2a29, https://open.canada.ca/data/en/dataset/8661edcf-f525-4758-a051-cb3fc8c74423). This monitoring data, as well additional occurrence information from online databases and the scientific literature, have been paired with high resolution environmental data and oceanographic models in species distribution models that predict the present-day and future potential distributions of 12 moderate to high risk invasive species on Canada’s east and west coasts. Future distributions were predicted for 2075, under Representative Concentration Pathway 8.5 from the Intergovernmental Panel on Climate Change’s fifth Assessment Report. Present-day and future richness of these species (i.e., hotspots) has also been estimated by summing their occurrence probabilities. This data set includes the occurrence locations of each species, the present-day and future species distribution modeling results for each species, and the estimated species richness. This research has been published in the scientific literature(Lyons et al. 2020).Lyons DA, Lowen JB, Therriault TW, Brickman D, Guo L, Moore AM, Peña MA, Wang Z, DiBacco C. (In Press) Identifying Marine Invasion Hotspots Using Stacked Species Distribution Models. Biological InvasionsCite this data as: Lyons DA., Lowen JB, Therriault TW., Brickman D., Guo L., Moore AM., Peña MA., Wang Z., DiBacco C. Data of: Species distribution models and occurrence data for marine invasive species hotspot identification. Published: November 2020. Coastal Ecosystems Science Division, Fisheries and Oceans Canada, Dartmouth, N.S. https://open.canada.ca/data/en/dataset/1bbd5131-8b34-4245-b999-3b4c4259d74f

Distribution of fulmar, shearwater and petrel bird species habitat in coastal British Columbia showing relative abundance (RA) by season and overall relative importance (RI). RI is based on project region and not on the province as a whole. CRIMS is a legacy dataset of BC coastal resource data that was acquired in a systematic and synoptic manner from 1979 and was intermittently updated throughout the years. Resource information was collected in nine study areas using a peer-reviewed provincial Resource Information Standards Committee consisting of DFO Fishery Officers, First Nations, and other subject matter experts. There are currently no plans to update this legacy data.

This dataset displays the Canadian geographic ranges of the priority species identified under the Pan-Canadian Approach for Transforming Species at Risk Conservation in Canada (“Pan-Canadian Approach”). These species include Barren-ground Caribou (including the Dolphin and Union population); Greater Sage-Grouse; Peary Caribou; Wood Bison; Caribou, Boreal population (“Boreal Caribou”); and Woodland Caribou, Southern Mountain population (“Southern Mountain Caribou”). The priority species were chosen following a number of criteria and considerations in collaboration with federal, provincial, and territorial partners. These include, but were not limited to, the species' ecological role on a regional or national scale, their conservation status and achievability of conservation outcomes, their social and cultural value (particularly to Indigenous peoples), and the leadership/partnership opportunities that they present. Delivering conservation outcomes for targeted priority species can have significant co-benefits for other species at risk, and wildlife in general. For more information on the Pan-Canadian Approach and the priority species, see https://www.canada.ca/en/services/environment/wildlife-plants-species/species-risk/pan-canadian-approach.html.This dataset includes: 1) the range for the Boreal Caribou (see https://species-registry.canada.ca/index-en.html#/consultations/2253); 2) the local populations for the Southern Mountain Caribou (see https://species-registry.canada.ca/index-en.html#/consultations/1309); 3) the range for the Greater Sage-Grouse (see https://species-registry.canada.ca/index-en.html#/consultations/1458); 4) local populations for the Peary Caribou (see https://species-registry.canada.ca/index-en.html#/consultations/3657); 5) range for the Barren-ground Caribou (see https://www.maps.geomatics.gov.nt.ca/Html5Viewer/index.html?viewer=NWT_SHV English only); 6) range for the Barren-ground Caribou, Dolphin and Union population (https://www.maps.geomatics.gov.nt.ca/Html5Viewer/index.html?viewer=NWT_SHV English only); 7) range for the Wood Bison (see https://species-registry.canada.ca/index-en.html#/consultations/2914).

Monitoring data from DFO invasive species monitoring programs, along with occurrence information from online databases and the scientific literature, have been paired with high resolution environmental data and oceanographic models in species distribution models that predict present-day and project future distributions of 24 non-indigenous species (NIS) on North America`s east coast, and 31 NIS on its west coast. Future distributions were predicted for 2100, under Representative Concentration Pathway 8.5 from the Intergovernmental Panel on Climate Change’s fifth Assessment Report. Present-day and future richness of these species (i.e., hotspots) have been estimated by summing the occurrence probabilities of NIS. This data set includes the present-day and year 2100 species distribution modeling results for each species, and the estimated species richness.Cite this data as: Lyons DA., Lowen JB, Therriault TW., Brickman D., Guo L., Moore AM., Peña MA., Wang Z., DiBacco C. Data of: Updated species distribution models for marine invasive species hotspot identification. Published: November 2023. Coastal Ecosystems Science Division, Fisheries and Oceans Canada, Dartmouth, N.S. https://open.canada.ca/data/en/dataset/1439dcb3-82a6-40fd-a9a4-8f045b20ff5b

The Aquatic Invasive Species Surveillance Database is a compilation of fish community and habitat data from DFO’s Aquatic Invasive Species and Invasive Carp Program early detection surveillance efforts in Canadian waters of the Great Lakes basin. Data includes: sampling site location, date, fish species and counts, and associated habitat information. Annual project-specific details including purpose/objectives and study methodology are often reported in the DFO Canadian manuscript report of fisheries and aquatic sciences series.

In this dataset, we share maps of annual dominant tree species (also known as leading tree species) from 1984-2022 covering the entirety of Canada's 650 Mha forested ecosystems using Landsat time-series imagery at a 30-m spatial resolution. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). Classifications are based on regionally representative Random Forests model using local training samples from Canada's National Forest Inventory (Hermosilla et al., 2024). Descriptive metrics provide information on spectral, geographic, climatic, and topographic characteristics. Initial annual tree species classifications were subjected to a time series post-classification process using the forward-backward Hidden Markov Model to improve the temporal consistency of tree species transitions within the time series. Assessment of the annual species maps using independent validation data resulted in an overall accuracy of 86.1% ± 0.14% (95%-confidence interval). These data allow consistent comparison of trends and rates of change in tree species composition nationally and across regions using a common time frame, spatial resolution, and analytical approach.Hermosilla, T., Wulder, M.A., White, J.C., Coops, N.C., Bater, C.W., Hobart, G.W., 2024. Characterizing long-term tree species dynamics in Canada's forested ecosystems using annual time series remote sensing data. Forest Ecology and Management, 122313. https://doi.org/10.1016/j.foreco.2024.122313 (Hermosilla et al. 2024)

Description:This dataset contains layers of predicted occurrence for 65 groundfish species as well as overall species richness (i.e., the total number of species present) in Canadian Pacific waters, and the median standard error per grid cell across all species. They cover all seafloor habitat depths between 10 and 1400 m that have a mean summer salinity above 28 PSU. Two layers are provided for each species: 1) predicted species occurrence (prob_occur) and 2) the probability that a grid cell is an occurrence hotspot for that species (hotspot_prob; defined as being in the lower of: 1) 0.8, or 2) the 80th percentile of the predicted probability of occurrence values across all grid cells that had a probability of occurrence greater than 0.05.). The first measure provides an overall prediction of the distribution of the species while the second metric identifies areas where that species is most likely to be found, accounting for uncertainty within our model. All layers are provided at a 1 km resolution.Methods:These layers were developed using a species distribution model described in Thompson et al. 2023. This model integrates data from three fisheries-independent surveys: the Fisheries and Oceans Canada (DFO) Groundfish Synoptic Bottom Trawl Surveys (Sinclair et al. 2003; Anderson et al. 2019), the DFO Groundfish Hard Bottom Longline Surveys (Lochead and Yamanaka 2006, 2007; Doherty et al. 2019), and the International Pacific Halibut Commission Fisheries Independent Setline Survey (IPHC 2021). Further details on the methods are found in the metadata PDF available with the dataset.Abstract from Thompson et al. 2023:Predictions of the distribution of groundfish species are needed to support ongoing marine spatial planning initiatives in Canadian Pacific waters. Data to inform species distribution models are available from several fisheries-independent surveys. However, no single survey covers the entire region and different gear types are required to survey the range of habitats that are occupied by groundfish. Bottom trawl gear is used to sample soft bottom habitat, predominantly on the continental shelf and slope, whereas longline gear often focuses on nearshore and hardbottom habitats where trawling is not possible. Because data from these two gear types are not directly comparable, previous species distribution models in this region have been limited to using data from one survey at a time, restricting their spatial extent and usefulness at a regional scale. Here we demonstrate a method for integrating presence-absence data across surveys and gear types that allows us to predict the coastwide distributions of 66 groundfish species in British Columbia. Our model leverages the use of available data from multiple surveys to estimate how species respond to environmental gradients while accounting for differences in catchability by the different surveys. Overall, we find that this integrated method has two main benefits: 1) it increases the accuracy of predictions in data-limited surveys and regions while having negligible impacts on the accuracy when data are already sufficient to make predictions, 2) it reduces uncertainty, resulting in tighter confidence intervals on predicted species occurrences. These benefits are particularly relevant in areas of our coast where our understanding of habitat suitability is limited due to a lack of spatially comprehensive long-term groundfish research surveys.Data Sources:Research data was provided by Pacific Science’s Groundfish Data Unit for research surveys from the GFBio database between 2003 and 2020 for all species which had at least 150 observations, across all gear type and survey datasets available.Uncertainties:These are modeled results based on species observations at sea and their related environmental covariate predictions that may not always accurately reflect real-world groundfish distributions though methods that integrate different data types/sources have been demonstrated to improve model inference by increasing the accuracy of the predictions and reducing uncertainty.

A derivative of DFO’s benthic species survey for the Strategic Program for Ecosystem-based Research and Advice (SPERA) (open data record ID: e736c0f0-b19e-4842-903d-28bfc756d48a), this benthic survey funded through the Canadian Healthy Oceans Network (CHONeII) looks at the presence/absence and abundance of two biogenic habitat-forming species that are listed as vulnerable to disturbance in a subset of 50 drift camera transects in the ‘Head Harbour/West Isles Archipelago/The Passages’ Ecologically and Biologically Significant Area (EBSA) in the Bay of Fundy, New Brunswick, Canada (~113km2). Presence/absence and abundance data of the stalked sea squirt (Boltenia ovifera) and horse mussel (Modiolus modiolus) were derived from the use of high-resolution Nikon D800 36.1 megapixel still images (n=2576, see link to parent record for more descriptive survey information) to be used in species distribution modelling. Image field of view (FOV) was estimated using a 10 cm-wide trigger weight for scale,and standardized across images using the average FOV estimate (0.75 x 0.5 m) across a subset of 200 images. Species counts were then converted to abundance estimates (number of individuals per square-meter) by dividing counts by 0.375m2. Boltenia ovifera was observed at densities reaching 456 ind./m2, while Modiolus modiolus density reached a maximum of 240 ind./m2.Cite this data as: Mireault C.A., Lawton P., Devillers R. and Teed L. Presence/absence and abundance of vulnerable marine ecosystem species Boltenia ovifera and Modiolus modiolus in the lower Bay of Fundy derived from high resolution still imagery. Published September 2023. Coastal Ecosystems Science Division, Fisheries and Oceans Canada, St. Andrews, N.B. https://open.canada.ca/data/en/dataset/152ae3f1-d2b9-43d9-a7b4-d769d9e9fc41

The Science Branch of Fisheries and Oceans Canada (DFO) in the Newfoundland and Labrador (NL) region has been conducting multispecies research vessel (RV) surveys using a stratified random survey design since the early 1970s. The DFO RV survey dataset represents the longest time series of species data in the NL region, making it ideal for mapping the average relative densities of species over time. Average relative density maps depict the interpolated densities (calculated from kg/tow) of fish species or functional groups. These densities are averaged over each time series (Engel and Campelen) and include data from all available seasons, so they represent persistent areas of relatively high and low densities for that species or functional group for the duration of the time series, independent of season. These maps are well suited as decision support tools related to conservation areas and marine spatial planning. These maps can also inform other processes that require information on areas important to marine fish, such as environmental assessments. Spring, fall, and winter data from the DFO RV survey between 1981 and 2017, inclusive, were used for the analysis. Due to a gear change from an Engel 145 Hi-Lift Otter Trawl to a Campelen 1800 Shrimp Trawl in 1995, the time series is treated as two separate datasets. NAFO Divisions 2J3KLNOP were sampled during the Engel time series and Division 2H was added for the Campelen time series. The data were filtered prior to use so that only core strata (areas consistently sampled across years) were included, resulting in most deep water and inshore sets being excluded in this analysis. Weight per tow (kg/tow; standardized for tow length for each gear type) data for fish, shrimp, and crab species were extracted from the database, and all successful sets from regular multispecies surveys were used for analyses. Eight fish functional groups (groups of species of similar size and diet) were identified based on the RV survey dataset: small benthivores, medium benthivores, large benthivores, piscivores, plank-piscivores, planktivores, shrimp, and forage fish. Data for each functional group were mapped three ways: all species, dominant species (i.e. top 90% biomass), and non-dominant species. In total, 40 dominant species and/or at-risk species (i.e. COSEWIC endangered, threatened, special concern; SARA; DFO/NAFO depleted) were mapped individually. To identify the average relative density, independent of seasonality, the spring, fall, and winter survey sets were compiled into a composite dataset using a log transformation on the biomass (kg/tow). For functional groups, these values were then standardized across each group. Absences (0 kg/tow catch values) were included. A continuous raster with a 4x4km resolution was generated through ordinary kriging. The raster was clipped to an 8-km buffer of the RV survey extent and the zero values were then removed. The results of this process are maps depicting the average relative density of fish functional groups and selected individual species during both the Engel (1981-1995) and Campelen (1995-2017) time series. Note that the original units (e.g. kg/tow) are no longer relevant due to data processing. Cell values are not comparable between groups or species; when mapping, all numeric values should be removed from the labels and legend, with relative qualifiers (“high” and “low”) used instead. More detailed information can be found in Wells et al. (2021). References: Wells, N.J., Pretty, C., Warren, M., Novaczek, E. and Koen-Alonso, M. 2021. Average Relative Density of Fish Species and Functional Groups in the Newfoundland and Labrador Shelves Bioregion from 1981-2017. Can. Tech. Rep. Fish. Aquat. Sci. 3427: viii + 76 p.

DescriptionConservation of marine biodiversity requires understanding the joint influence of ongoing environmental change and fishing pressure. Addressing this challenge requires robust biodiversity monitoring and analyses that jointly account for potential drivers of change. Here, we ask how demersal fish biodiversity in Canadian Pacific waters has changed since 2003 and assess the degree to which these changes can be explained by environmental change and commercial fishing. Using a spatiotemporal multispecies model based on fisheries independent data, we find that species density (number of species per area) and community biomass have increased during this period. Environmental changes during this period were associated with temporal fluctuations in the biomass of species and the community as a whole. However, environmental changes were less associated with changes in species’ occurrence. Thus, the estimated increases in species density are not likely to be due to environmental change. Instead, our results are consistent with an ongoing recovery of the demersal fish community from a reduction in commercial fishing intensity from historical levels. These findings provide key insight into the drivers of biodiversity change that can inform ecosystem-based management.The layers provided represent three community metrics: 1) species density (i.e., species richness), 2) Hill-Shannon diversity, and 3) community biomass. All layers are provided at a 3 km resolution across the study domain for the period of 2003 to 2019. For each metric, we provide layers for three summary statistics: 1) the mean value in each grid cell over the temporal range, 2) the probability that the grid cell is a hotspot for that metric, and 3) the temporal coefficient of variation (i.e., standard deviation/mean) across all years.Methods:The analysis that produced these layers is presented in Thompson et al. (2022). The analysis uses data from the Groundfish Synoptic Bottom Trawl Research surveys in Queen Charlotte Sound (QCS), Hecate Strait (HS), West Coast Vancouver Island (WCVI), and West Coast Haida Gwaii (WCHG) from 2003 to 2019. Cartilaginous and bony fish species caught in DFO groundfish surveys that were present in at least 15% of all trawls over the depth range in which they were caught were included. This depth range was defined as that which included 95% of all trawls in which that species was present. The final dataset used in our analysis consisted of 57 species (Table S1 in Thompson et al. 2022).The spatiotemporal dynamics of the demersal fish community were modeled using the Hierarchical Modeling of Species Communities (HMSC) framework and package (Tikhonov et al. 2021) in R. This framework uses Bayesian inference to fit a multivariate hierarchical generalized mixed model. We modeled community dynamics using a hurdle model, which consists of two sub models: a presence-absence model and a biomass model that is conditional on presence. Our list of environmental covariates included bottom depth, bathymetric position index (BPI), mean summer tidal speed, substrate muddiness, substrate rockiness, whether the trawl was inside or outside of the ecosystem-based trawling footprint, and survey region (QCS & HS vs. WCVI & WCHG)), mean summer near-bottom temperature deviation, mean summer near-bottom dissolved oxygen deviation, mean summer cross-shore and along-shore current velocities near the seafloor, mean summer depth-integrated primary production, and local-scale commercial fishing effort.Layers are provided for three community metrics. All metrics should be interpreted as the value that would be expected in the catch from an average tow in the Groundfish Synoptic Bottom Trawl Research Surveys taken in a given 3 km grid cell. Species density (sometimes called species richness) should be interpreted as the number of the 57 species that would be caught in a trawl. Hill-Shannon diversity is a measure of diversity that gives greater weight to communities where biomass is spread equally across species. Community biomass is the total biomass across all 57 species that would be expected to be caught per square km in an average tow. Data Sources:Research data was provided by Pacific Science's Groundfish Data Unit for research surveys from the GFBio database between 2003 and 2019 that occurred in four regions: Queen Charlotte Sound, Hecate Strait, West Coast Haida Gwaii, and West Coast Vancouver Island. Our analysis excludes species that are rarely caught in the research trawls and so our estimates would not include the occurrence or biomass of these rare species.Commercial fishing data was accessed through a DFO R script detailed here: https://github.com/pbsassess/gfdata. Local scale commercial fishing effort was calculated from this data. The substrate layers were obtained from a substrate model (Gregr et al. 2021). The oceanographic layers (bottom temperature, dissolved oxygen, tidal and circulation speeds, primary production) were obtained from a hindcast simulation of the British Columbia continental margin (BCCM) model (Peña et al. 2019).Uncertainties:Species that are not well sampled by the trawl surveys may not be accurately estimated by our model. The model did not include spatiotemporal random effects, which likely underestimates spatiotemporal variability in the region. It is also important to underline covariate uncertainty and model uncertainty. The hotspot estimates provide one measure of model uncertainty/certainty.