“Pacific Rim National Park Reserve - Total Ecosystem Forest Carbon Density” is the annual carbon density (tonnes carbon per hectare) within Pacific Rim’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.

“Terra Nova 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 Terra Nova 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.

Exposure to wind-driven waves manifests an important physical gradient in the coastal zone that influences a variety of physical and biological processes (i.e., species distribution). Fetch, the unobstructed distance over which wind-driven waves can build, is a popular proxy for wave exposure at a given location commonly used for site-specific evaluations. Here, we provide two measures of fetch (unweighted fetch, effective fetch) and three fetch-derived indices of wave exposure (sum fetch, minimum fetch, and a relative exposure index) covering the coastal zones of two Canadian bioregions (Scotian Shelf-Bay of Fundy, Newfoundland-Labrador Shelves). For each region, we calculated fetch and exposure indices for input points in an evenly spaced fishnet grid (see linked records below for datasets by region). We provide unweighted fetch lengths (m) for 32 compass headings per input point (11.25° intervals), and effective fetch lengths for 8 headings per point (45° intervals). Effective fetch is a weighted average of multiple fetch measures around a given heading that reduces the influence of irregular coastline shape on exposure estimates. We also include calculations of the summed and minimum unweighted fetch lengths for each point that provide coarse proxies of exposure and distance to land, respectively. The relative exposure index (REI), provided as regional spatial layers in raster format, provides a more accurate metric of exposure by combining effective fetch with modelled wind speeds (m s-1) and frequency data. Users may also use fetch calculations to derive their own exposure layers using alternative sources of wind data, indices, or formulations. Detailed methodology on the calculations for fetch, effective fetch and REI are outlined in the Supplementary Information below. Citation information and differences in data pre-processing methods and spatial resolution of the regional analyses are described in their respective data records. The broad spatial coverage and high resolution offered by these indices are suitable to support regional-scale modelling and planning exercises. In particular, these indices will be of value to ongoing Marine Spatial Planning efforts, which includes regional conservation planning, that seek to evaluate the distribution of coastal species and overlap with human activities.

“Jasper National Park - Total Ecosystem Forest Carbon Density” is the annual carbon density (tonnes carbon per hectare) within Jasper’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.

“Gros Morne National Park - Total Ecosystem Forest Carbon Density” is the annual carbon density (tonnes carbon per hectare) within Gros Morne’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.

“Wood Buffalo National Park - Total Ecosystem Forest Carbon Density” is the annual carbon density (tonnes carbon per hectare) within Wood Buffalo’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.

“Forillon National Park - Total Ecosystem Forest Carbon Density” is the annual carbon density (tonnes carbon per hectare) within Forillon’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.

The Crop Stress Index is the ratio of actual evapotranspiration (AET) to potential evapotranspiration (PET) express as: CSI = 1-(AET/PET) AET and PET are calculated within the Versatile Soil Moisture Budget (VSMB) model using temperature and precipitation data and a crop-specific biometeorological time scale model to estimate growth stage (Robertson, 1968), with crop specific phenological and crop water extraction coefficients taken from Chipanshi et al 2013. The WDI ranges between 0 and 1, with a value closer to 1 indicating higher stress Crop Stress Index is modelled for each climate station using measured precipitation and temperature

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.