Land cover, 2019
Follow one 27 × 27 km window in the Churchill range. Every 30 m pixel carries a class from VLCE2, NRCan's national forest land-cover map.
Woodland caribou · Ontario · 2015–2025
Woodland caribou need large areas of undisturbed forest. We measured how much is left in seven Ontario ranges, every year from 2015 to 2025.
Caribou: Alain Caron, Wikimedia Commons, CC BY-SA 4.0 · Earth: NASA Blue Marble
Research questions
Scientific basis
Boreal woodland caribou are listed as Threatened in Canada. They persist where large areas of mature forest remain undisturbed. Logging, roads and fire create young forest that supports moose and deer; these sustain more wolves, and caribou calves are lost to the added predation (Environment Canada 2011).
Overview
Canada considers a caribou range with more than 35% of its area disturbed unlikely to keep a self-sustaining herd. Logging, roads and wildfire open up the forest and bring in moose, wolves and bears.
Managed forest, lakes included is the published study’s basis and our default. Whole range is how Environment and Climate Change Canada and Ontario report. Land only leaves out open water, where caribou do not live. See how our numbers compare.
The rings and range card at the top show this split for any year.
Caribou range
Discussion
Each source measures a different area. Switch between the managed forest, as in the published study, and the whole range, as in the federal and Ontario reports. Striped bars are figures we measured from that source's own maps or rebuilt with its rules.
References
Data sources and their providers are listed under Datasets.
Method and workflow
We use the method of Mackey et al. (2024). First we rebuilt their 2020 map from their data. Then we added newer logging (ARI) and fire (NBAC) records. The Data tab shows how we checked them.
Each circle is one type of data. The label on each line is the rule we use. Point at a circle to learn more.
Follow one 27 × 27 km window in the Churchill range. Every 30 m pixel carries a class from VLCE2, NRCan's national forest land-cover map.
Only conifer forest and treed wetland count as potential habitat. Other land-cover classes are excluded before disturbance is assessed.
Potential habitat: – of the window
A pixel is disturbed if any dataset maps disturbance there: logging since 1980 from FRI or NTEMS harvest, planted stands, roads, mines, towns, farms, corridors and OLCC's disturbance class, each with a 500 m buffer zone, plus wildfire. Overlaps count once.
Disturbed: – of the window
What remains is potential habitat that no disturbance reaches. In this window, – of the potential habitat is lost.
Undisturbed habitat: – of the window
A 30 m edge is removed and adjacent pixels are grouped into patches. Patches under 100 ha are excluded; the rest are grouped at 500, 1,000 and 10,000 ha.
Core habitat: – of the window · – of it in patches over 1,000 ha
The reproduction
Reproduction. Same inputs as the study: 17 Forest Resources Inventory packages, the national Landsat harvest layer 1985–2020, MNRF and national road networks, the Ontario Land Cover Compilation (classes 18, 25, 27, 28), and the National Fire Database 1985–2018. We checked our map pixel by pixel against the layers the authors shared. The remaining 1 to 3 points come mostly from roads, which are only published as current editions.
Overlap. On the 30 m grid inside each range's managed forest, every dataset's mapped area is split into the part also mapped by all the other datasets of the same type, by some of them, or by none. Another dataset counts as mapping a pixel when it maps something within 60 m, so that polygon edges, Landsat pixels and road centre lines drawn slightly apart still match. Comparisons that include ECCC 2020 are limited to the federal range, where ECCC maps.
Rules for the 2015–2025 maps. Fire counts from 1985 in every year (NFDB and NBAC). ARI logging from 2013 on is used only where at least half of a stand shows a clear drop in the forest signal on Landsat or Sentinel-2; 2,609 stands (40,946 ha) without a visible cut are excluded. OLCC's disturbance class (18) mixes cuts and burns, so where NFDB or NBAC maps a fire from 1985 on it is counted as fire (no 500 m buffer zone); everywhere else it stays anthropogenic, as in the published study.
Satellite check. Random spots inside ARI records: about 400 per record type on Landsat (logging only where FRI and NTEMS have none) and 190 to 350 per type on Sentinel-2, with cloud-free summer images from two years before to three years after the recorded year (Landsat 8/9 for 2013–2020, Sentinel-2 for 2018–2024). A spot shows forest loss if its Normalized Burn Ratio falls by more than 0.1 at some point in that six-year window, so cuts recorded a year or two early or late still count. Logging recorded by both FRI and NTEMS harvest is checked the same way as a reference: 94% of it shows loss, and the rest is mostly partial cuts, cuts smaller than a pixel, or stands already open. The full check then tested every ARI logging stand from 2013 on over its whole area, with Landsat and, from 2018, Sentinel-2. Blowdown and drought records are excluded because Canada's caribou method counts only human disturbance and fire.
Code and notebooks: github.com/RSL-SLU/caribou
Data sources
We start from the data a published study used, Mackey et al. (2024), and add two datasets that carry it to 2025: ARI for logging and NBAC for fire. Satellite images are used only to check the records. Below: which years each source covers, two checks of the data, and every dataset by name.
Check 1 · Agreement between datasets
For each disturbance type, datasets were compared pixel by pixel on the 30 m grid, with a 60 m tolerance. Each card gives the share of one dataset that the others also map.
Switch datasets on or off. Places mapped by two or more are green.
used goes into the habitat maps · check only is compared with them but not counted.
A high “mostly unique” share can reflect a more complete dataset, such as MNRF roads with its forest roads.
Check 2 · are the new records real?
NBAC fire is not tested here: it is itself mapped from satellite images, so the test would confirm it by construction. Check 1 shows it matches the NFDB and ARI fire records instead.
Sampled records: ARI records drawn at random for a point test, shown by range and record type; filled where loss is detected. All ARI logging stands: every logging stand recorded from 2013 to 2023, filled where the full-stand test retains it. Hover a shape for its forest signal; double-click to play its summers.
Datasets
The short name in the first column is the one used throughout this study.
| Name | Dataset and provider | Years | Format | Used for | Role |
|---|---|---|---|---|---|
| Ranges | Woodland caribou range boundaries Land Information Ontario | current | polygons | The seven ranges analysed | Published study |
| FMU | Forest Management Units Land Information Ontario | 2025 | polygons | The managed forest; every percentage is of the range inside it | Published study |
| VLCE2 | Canada forest land cover 2019 (Virtual Land Cover Engine) Natural Resources Canada | 2019 | 30 m raster | Potential caribou habitat: conifer forest and treed wetland | Published study |
| FRI | Forest Resources Inventory, Packaged Products v2 (17 management units) Ontario MNRF | logging 1980 to each inventory year | polygons | Logging and planted stands; unclassified corridors | Published study |
| NTEMS harvest | Landsat-derived forest harvest 1985–2020 (Hermosilla et al.) Natural Resources Canada | 1985–2020 | 30 m raster | Logging; its year is used where FRI and NTEMS harvest both record logging | Published study |
| MNRF roads | MNRF Road Segments Land Information Ontario | by year built, to 2025 | lines | Roads, 500 m buffer zone | Published study |
| NRN | National Road Network v18, Ontario Statistics Canada | current | lines | Roads, 500 m buffer zone | Published study |
| OLCC | Ontario Land Cover Compilation v2.0 (classes 18, 25, 27, 28) Ontario MNRF | imagery 2000–2011 | 15 m raster | Mines, towns, farms and land-cover disturbance, 500 m buffer zone. Its disturbance class (18) mixes cuts and burns: inside a fire from 1985 on (NFDB or NBAC) it counts as fire | Published study |
| NFDB | National Fire Database fire polygons Natural Resources Canada (CWFIS) | 1972–2024 | polygons | Fire from 1985 (1985–2018 for the reproduction) | Published study |
| ARI | Analysis Ready Inventory 2025 Ontario MNRF | annual reports to 2023 | polygons | Logging the other sources miss, and all logging after 2020 (from 2013 on, only stands seen from space); its fire and blowdown records are only compared | Added to reach 2025 |
| NBAC | National Burned Area Composite Natural Resources Canada (CWFIS) | 1980–2025 | polygons | Fire from 1985, including 2025; mapped from satellite images | Added to reach 2025 |
| NTEMS harvest 2022 | Landsat-derived forest harvest 1985–2022 Natural Resources Canada | 1985–2022 | 30 m raster | Checking ARI after 2020 and its year | Check only |
| ECCC 2020 | Anthropogenic disturbance footprint 2020, from 30 m Landsat-8 Environment and Climate Change Canada | 2020 | lines and polygons | Checking logging, roads, settlements and mines; compared with our results | Check only |
| Landsat | Collection 2 Level-2 surface reflectance NASA / USGS | 2011–2025 images | 30 m images | Checking ARI logging from 2013 on | Check only |
| Sentinel-2 | Level-2A surface reflectance ESA Copernicus | 2016–2025 images | 10–20 m images | Checking ARI logging from 2018 on; satellite base map (summer 2024) | Check only |