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.

“Thousand Islands National Park - Total GHG Emissions” datasets consist of estimates of GHG emissions (carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O)) in carbon dioxide equivalents (CO2e) from forested ecosystems in Thousand Islands National Park from 1990 to 2020 (tonnes carbon dioxide equivalent per hectare). Total GHG emissions for 31 national parks were estimated using the Generic Carbon Budget Model (GCBM), a spatially explicit carbon budget model developed by Canadian Forest Service which uses forest inventory, disturbance, and mean annual temperature data along with yield data to estimate growth and merchantable volume for dominant tree species. Species- and Ecozone-specific equations are then used to convert merchantable volume to aboveground and belowground biomass carbon. The GCBM simulates carbon dynamics to produce spatially explicit estimations of carbon stocks and fluxes. The model simulates and tracks carbon stocks, transfers between Intergovernmental Panel on Climate Change (IPCC)-defined pools, and other metrics including net ecosystem production, net biome production, and emissions of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) in annual time steps. The stocks and fluxes are also tracked by disturbance event (e.g., forest fires). Total GHG emissions include those from natural processes like respiration and decomposition and those due to natural and anthropogenic disturbances, including wildfires, prescribed burns, and insect outbreaks. These were calculated as the sum of CO2, CH4, and N2O emission estimates in tonnes carbon (tonnes C) generated by the GCBM. Emissions estimates were then converted to carbon dioxide equivalents (CO2e) using the 100-year Global Warming Potential (IPCC Fourth Assessment Report) factors for CH4 (25) and N2O (298). These products have a spatial resolution of 30m. This information is part of the Parks Canada Carbon Atlas Series. To obtain a copy of this report, please contact changementclimatique-climatechange@pc.gc.ca. When using this data, please cite as follows: Sharma, T., Kurz, W.A., Fellows, M., MacDonald, A.L., Richards, J., Chisholm, C., Seutin, G., Richardson, K., Keenleyside, K. (2023). Parks Canada Carbon Atlas Series: Carbon Dynamics in the Forests of Canada’s National Parks. Scientific Report. Parks Canada Agency, Gatineau, QC, Canada, 104 p.

“Gwaii Haanas National Park Reserve – Total GHG Emissions” datasets consist of estimates of GHG emissions (carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O)) in carbon dioxide equivalents (CO2e) from forested ecosystems in Gwaii Haanas National Park Reserve from 1990 to 2020 (tonnes carbon dioxide equivalent per hectare). Total GHG emissions for 31 national parks were estimated using the Generic Carbon Budget Model (GCBM), a spatially explicit carbon budget model developed by Canadian Forest Service which uses forest inventory, disturbance, and mean annual temperature data along with yield data to estimate growth and merchantable volume for dominant tree species. Species- and Ecozone-specific equations are then used to convert merchantable volume to aboveground and belowground biomass carbon. The GCBM simulates carbon dynamics to produce spatially explicit estimations of carbon stocks and fluxes. The model simulates and tracks carbon stocks, transfers between Intergovernmental Panel on Climate Change (IPCC)-defined pools, and other metrics including net ecosystem production, net biome production, and emissions of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) in annual time steps. The stocks and fluxes are also tracked by disturbance event (e.g., forest fires). Total GHG emissions include those from natural processes like respiration and decomposition and those due to natural and anthropogenic disturbances, including wildfires, prescribed burns, and insect outbreaks. These were calculated as the sum of CO2, CH4, and N2O emission estimates in tonnes carbon (tonnes C) generated by the GCBM. Emissions estimates were then converted to carbon dioxide equivalents (CO2e) using the 100-year Global Warming Potential (IPCC Fourth Assessment Report) factors for CH4 (25) and N2O (298). These products have a spatial resolution of 30m. This information is part of the Parks Canada Carbon Atlas Series. To obtain a copy of this report, please contact changementclimatique-climatechange@pc.gc.ca. When using this data, please cite as follows: Sharma, T., Kurz, W.A., Fellows, M., MacDonald, A.L., Richards, J., Chisholm, C., Seutin, G., Richardson, K., Keenleyside, K. (2023). Parks Canada Carbon Atlas Series: Carbon Dynamics in the Forests of Canada’s National Parks. Scientific Report. Parks Canada Agency, Gatineau, QC, Canada, 104 p.

“Yoho National Park - Total GHG Emissions” datasets consist of estimates of GHG emissions (carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O)) in carbon dioxide equivalents (CO2e) from forested ecosystems in Yoho National Park from 1990 to 2020 (tonnes carbon dioxide equivalent per hectare). Total GHG emissions for 31 national parks were estimated using the Generic Carbon Budget Model (GCBM), a spatially explicit carbon budget model developed by Canadian Forest Service which uses forest inventory, disturbance, and mean annual temperature data along with yield data to estimate growth and merchantable volume for dominant tree species. Species- and Ecozone-specific equations are then used to convert merchantable volume to aboveground and belowground biomass carbon. The GCBM simulates carbon dynamics to produce spatially explicit estimations of carbon stocks and fluxes. The model simulates and tracks carbon stocks, transfers between Intergovernmental Panel on Climate Change (IPCC)-defined pools, and other metrics including net ecosystem production, net biome production, and emissions of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) in annual time steps. The stocks and fluxes are also tracked by disturbance event (e.g., forest fires). Total GHG emissions include those from natural processes like respiration and decomposition and those due to natural and anthropogenic disturbances, including wildfires, prescribed burns, and insect outbreaks. These were calculated as the sum of CO2, CH4, and N2O emission estimates in tonnes carbon (tonnes C) generated by the GCBM. Emissions estimates were then converted to carbon dioxide equivalents (CO2e) using the 100-year Global Warming Potential (IPCC Fourth Assessment Report) factors for CH4 (25) and N2O (298). These products have a spatial resolution of 30m. This information is part of the Parks Canada Carbon Atlas Series. To obtain a copy of this report, please contact changementclimatique-climatechange@pc.gc.ca. When using this data, please cite as follows: Sharma, T., Kurz, W.A., Fellows, M., MacDonald, A.L., Richards, J., Chisholm, C., Seutin, G., Richardson, K., Keenleyside, K. (2023). Parks Canada Carbon Atlas Series: Carbon Dynamics in the Forests of Canada’s National Parks. Scientific Report. Parks Canada Agency, Gatineau, QC, Canada, 104 p.

“Prince Albert National Park - Total Ecosystem Forest Carbon Density” is the annual carbon density (tonnes carbon per hectare) within Prince Albert’s forested ecosystems over a 31-year period from 1990 to 2020. Total Ecosystem Forest Carbon Density includes aboveground and belowground biomass, soil carbon, and dead organic matter. Total Ecosystem Forest Carbon Density was estimated for 31 national parks using the Generic Carbon Budget Model (GCBM), a spatially explicit carbon budget model developed by Canadian Forest Service which uses forest inventory, disturbance, and mean annual temperature data along with yield data to estimate growth and merchantable volume for dominant tree species. Species- and Ecozone-specific equations are then used to convert merchantable volume to aboveground and belowground biomass carbon. Ecozones were classified according to Canada Ecological Land Classification Level 1. The GCBM simulates carbon dynamics to produce spatially explicit estimations of carbon stocks and fluxes. The model simulates and tracks carbon stocks, transfers between Intergovernmental Panel on Climate Change-defined pools, and other metrics including net ecosystem production, net biome production, and emissions of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) in annual time steps. The stocks and fluxes are also tracked by disturbance event (e.g., forest fires, insect outbreaks). Total Ecosystem Forest Carbon Density accounts for the effects of natural and anthropogenic disturbances, including wildfires, prescribed burns, and insect outbreaks. These products have a spatial resolution of 30m. This information is part of the Parks Canada Carbon Atlas Series. To obtain a copy of this report, please contact changementclimatique-climatechange@pc.gc.ca. When using this data, please cite as follows: Sharma, T., Kurz, W.A., Fellows, M., MacDonald, A.L., Richards, J., Chisholm, C., Seutin, G., Richardson, K., Keenleyside, K. (2023). Parks Canada Carbon Atlas Series: Carbon Dynamics in the Forests of Canada’s National Parks. Scientific Report. Parks Canada Agency, Gatineau, QC, Canada, 104 p.

“Kejimkujik National Park and National Historic Site - Total GHG Emissions” ddatasets consist of estimates of GHG emissions (carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O)) in carbon dioxide equivalents (CO2e) from forested ecosystems in Kejimkujik National Park and National Historic Site from 1990 to 2020 (tonnes carbon dioxide equivalent per hectare). Total GHG emissions for 31 national parks were estimated using the Generic Carbon Budget Model (GCBM), a spatially explicit carbon budget model developed by Canadian Forest Service which uses forest inventory, disturbance, and mean annual temperature data along with yield data to estimate growth and merchantable volume for dominant tree species. Species- and Ecozone-specific equations are then used to convert merchantable volume to aboveground and belowground biomass carbon. The GCBM simulates carbon dynamics to produce spatially explicit estimations of carbon stocks and fluxes. The model simulates and tracks carbon stocks, transfers between Intergovernmental Panel on Climate Change (IPCC)-defined pools, and other metrics including net ecosystem production, net biome production, and emissions of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) in annual time steps. The stocks and fluxes are also tracked by disturbance event (e.g., forest fires). Total GHG emissions include those from natural processes like respiration and decomposition and those due to natural and anthropogenic disturbances, including wildfires, prescribed burns, and insect outbreaks. These were calculated as the sum of CO2, CH4, and N2O emission estimates in tonnes carbon (tonnes C) generated by the GCBM. Emissions estimates were then converted to carbon dioxide equivalents (CO2e) using the 100-year Global Warming Potential (IPCC Fourth Assessment Report) factors for CH4 (25) and N2O (298). These products have a spatial resolution of 30m. This information is part of the Parks Canada Carbon Atlas Series. To obtain a copy of this report, please contact changementclimatique-climatechange@pc.gc.ca. When using this data, please cite as follows: Sharma, T., Kurz, W.A., Fellows, M., MacDonald, A.L., Richards, J., Chisholm, C., Seutin, G., Richardson, K., Keenleyside, K. (2023). Parks Canada Carbon Atlas Series: Carbon Dynamics in the Forests of Canada’s National Parks. Scientific Report. Parks Canada Agency, Gatineau, QC, Canada, 104 p.

Greenland Sharks (Somniosus microcephalus) are estimated to have the highest longevity of any invertebrate (392 ± 120 years), making bycatch a significant concern (Nielsen et al. 2016). However, in the Newfoundland and Labrador (NL) region, accurate estimates of bycatch are not often available for the species (Simpson et al. 2021). To address this, species distribution models (SDMs) were generated to delineate habitat suitability for the species throughout the NL region in order to identify areas where a higher rate of bycatch is expected to occur. Observations of Greenland shark bycatch recorded by At-Sea Observers (ASOs) in NL (1983 – 2019), Spain (1999 – 2017), and by the Northwest Atlantic Fisheries Organization (NAFO) Secretariat (2014 – 2019) were compiled to generate a presence-only dataset. Multiple environmental variables were assessed for collinearity, and non-collinear variables (Bathymetry and mean monthly bottom temperature for March and November (1990 – 2015)) were retained for use in the SDM. MaxEnt (maximum entropy) software was used to model habitat suitability because it is a presence-only modelling program that is able to account for a lack of absence data by comparing the environmental conditions at occurrence locations to those at randomly selected background points. Overall, the results indicated that habitat suitability for Greenland Shark was highest in deeper waters along the shelf edge in NAFO Divisions 3OP, and the Laurentian Channel, but also extended along the edge of the Labrador shelf, the Grand Banks, and deeper areas along the continental shelf such as the Hawke Channel, Funk Island Deep, and the slopes of Saglek, Nain, and Hamilton Banks. Beyond the Economic Exclusive Zone (EEZ) and within the NAFO regulatory area (NRA), suitable habitat was also present within the Flemish Pass, and along the slope of the Flemish Cap and shelf edge in NAFO Divisions 3NO (Simpson et al. 2021). More detailed information can be found in Simpson et al. 2021. References: Nielsen, J., R. B. Hedeholm, J. Heinemeier, P. G. Bushnell, J. S. Christiansen, J. Olsen, C. B. Ramsey, R. W. Brill, M. Simon, K. F. Steffensen, J. F. Steffensen. 2016. Eye lens radiocarbon reveals centuries of longevity in the Greenland shark (Somniosus microcephalus). Science 353 (6300):702-704 Simpson, M. R., Gullage, L., Konecny, C., Ollerhead, N., Treble, M.A., Nogueira, A., González-Costas, F. 2021. Spatial-temporal variation in Greenland shark (Somniosus microcephalus) bycatch in the NAFO Regulatory Area. NAFO SCR Doc. 21/028

“Elk Island National Park – Total GHG Emissions” datasets consist of estimates of GHG emissions (carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O)) in carbon dioxide equivalents (CO2e) from forested ecosystems in Banff National Park from 1990 to 2020 (tonnes carbon dioxide equivalent per hectare). Total GHG emissions for 31 national parks were estimated using the Generic Carbon Budget Model (GCBM), a spatially explicit carbon budget model developed by Canadian Forest Service which uses forest inventory, disturbance, and mean annual temperature data along with yield data to estimate growth and merchantable volume for dominant tree species. Species- and Ecozone-specific equations are then used to convert merchantable volume to aboveground and belowground biomass carbon. The GCBM simulates carbon dynamics to produce spatially explicit estimations of carbon stocks and fluxes. The model simulates and tracks carbon stocks, transfers between Intergovernmental Panel on Climate Change (IPCC)-defined pools, and other metrics including net ecosystem production, net biome production, and emissions of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) in annual time steps. The stocks and fluxes are also tracked by disturbance event (e.g., forest fires). Total GHG emissions include those from natural processes like respiration and decomposition and those due to natural and anthropogenic disturbances, including wildfires, prescribed burns, and insect outbreaks. These were calculated as the sum of CO2, CH4, and N2O emission estimates in tonnes carbon (tonnes C) generated by the GCBM. Emissions estimates were then converted to carbon dioxide equivalents (CO2e) using the 100-year Global Warming Potential (IPCC Fourth Assessment Report) factors for CH4 (25) and N2O (298). These products have a spatial resolution of 30m. This information is part of the Parks Canada Carbon Atlas Series. To obtain a copy of this report, please contact changementclimatique-climatechange@pc.gc.ca. When using this data, please cite as follows: Sharma, T., Kurz, W.A., Fellows, M., MacDonald, A.L., Richards, J., Chisholm, C., Seutin, G., Richardson, K., Keenleyside, K. (2023). Parks Canada Carbon Atlas Series: Carbon Dynamics in the Forests of Canada’s National Parks. Scientific Report. Parks Canada Agency, Gatineau, QC, Canada, 104 p.

“Kejimkujik National Park and National Historic Site - Total Ecosystem Forest Carbon Density” is the annual carbon density (tonnes carbon per hectare) within Kejimkujik’s forested ecosystems over a 31-year period from 1990 to 2020 (tonnes carbon per hectare). Total Ecosystem Forest Carbon Density includes aboveground and belowground biomass, soil carbon, and dead organic matter. Total Ecosystem Forest Carbon Density was estimated for 31 national parks using the Generic Carbon Budget Model (GCBM), a spatially explicit carbon budget model developed by Canadian Forest Service which uses forest inventory, disturbance, and mean annual temperature data along with yield data to estimate growth and merchantable volume for dominant tree species. Species- and Ecozone-specific equations are then used to convert merchantable volume to aboveground and belowground biomass carbon. Ecozones were classified according to Canada Ecological Land Classification Level 1. The GCBM simulates carbon dynamics to produce spatially explicit estimations of carbon stocks and fluxes. The model simulates and tracks carbon stocks, transfers between Intergovernmental Panel on Climate Change-defined pools, and other metrics including net ecosystem production, net biome production, and emissions of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) in annual time steps. The stocks and fluxes are also tracked by disturbance event (e.g., forest fires, insect outbreaks). Total Ecosystem Forest Carbon Density accounts for the effects of natural and anthropogenic disturbances, including wildfires, prescribed burns, and insect outbreaks. These products have a spatial resolution of 30m. This information is part of the Parks Canada Carbon Atlas Series. To obtain a copy of this report, please contact changementclimatique-climatechange@pc.gc.ca. When using this data, please cite as follows: Sharma, T., Kurz, W.A., Fellows, M., MacDonald, A.L., Richards, J., Chisholm, C., Seutin, G., Richardson, K., Keenleyside, K. (2023). Parks Canada Carbon Atlas Series: Carbon Dynamics in the Forests of Canada’s National Parks. Scientific Report. Parks Canada Agency, Gatineau, QC, Canada, 104 p.

“Pukaskwa National Park - Total GHG Emissions” datasets consist of estimates of GHG emissions (carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O)) in carbon dioxide equivalents (CO2e) from forested ecosystems in Pukaskwa National Park from 1990 to 2020 (tonnes carbon dioxide equivalent per hectare). Total GHG emissions for 31 national parks were estimated using the Generic Carbon Budget Model (GCBM), a spatially explicit carbon budget model developed by Canadian Forest Service which uses forest inventory, disturbance, and mean annual temperature data along with yield data to estimate growth and merchantable volume for dominant tree species. Species- and Ecozone-specific equations are then used to convert merchantable volume to aboveground and belowground biomass carbon. The GCBM simulates carbon dynamics to produce spatially explicit estimations of carbon stocks and fluxes. The model simulates and tracks carbon stocks, transfers between Intergovernmental Panel on Climate Change (IPCC)-defined pools, and other metrics including net ecosystem production, net biome production, and emissions of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) in annual time steps. The stocks and fluxes are also tracked by disturbance event (e.g., forest fires). Total GHG emissions include those from natural processes like respiration and decomposition and those due to natural and anthropogenic disturbances, including wildfires, prescribed burns, and insect outbreaks. These were calculated as the sum of CO2, CH4, and N2O emission estimates in tonnes carbon (tonnes C) generated by the GCBM. Emissions estimates were then converted to carbon dioxide equivalents (CO2e) using the 100-year Global Warming Potential (IPCC Fourth Assessment Report) factors for CH4 (25) and N2O (298). These products have a spatial resolution of 30m. This information is part of the Parks Canada Carbon Atlas Series. To obtain a copy of this report, please contact changementclimatique-climatechange@pc.gc.ca. When using this data, please cite as follows: Sharma, T., Kurz, W.A., Fellows, M., MacDonald, A.L., Richards, J., Chisholm, C., Seutin, G., Richardson, K., Keenleyside, K. (2023). Parks Canada Carbon Atlas Series: Carbon Dynamics in the Forests of Canada’s National Parks. Scientific Report. Parks Canada Agency, Gatineau, QC, Canada, 104 p.