Report and assessment of prescribed fire effects on bird communities in Canada’s National Parks

This report is dynamically generated, meaning its results may evolve with the addition of new data or further analyses. For the most recent updates, refer to the publication date and feel free to reach out to the authors.
Abstract
Ecologists use prescribed fire as a management tool to maintain and restore ecosystem structure, composition, and function, particularly in fire-adapted landscapes. In Canada’s national parks, prescribed fire is increasingly used to reduce fuel loads, promote habitat heterogeneity, protect communities, and support biodiversity objectives. However, understanding how bird communities respond to prescribed fire remains essential for evaluating both ecological benefits and potential impacts. Over the past three years, Parks Canada has collected standardized bird survey data within prescribed fire areas across 9 national parks, encompassing 236 sampling locations and approximately 44449 species observation records. This report presents summary analyses of these data to characterize migratory bird communities before and after prescribed fire in terms of species composition and richness, and functional diversity. In addition, the analysis evaluates differences between pre- and post-fire conditions and assesses the effectiveness of the sampling design, including its ability to detect change and its statistical power. The results are intended to help guide future bird monitoring efforts within Parks Canada for adaptive fire management practices.
Land Acknowledgement
This work was conducted in national parks located on the traditional lands of Indigenous, First Nations and Métis Peoples who have stewarded these landscapes since time immemorial. We recognize and respect the enduring relationships Indigenous communities maintain with the land, water, and wildlife, and we acknowledge the value of Indigenous knowledge systems that continue to inform understanding, care, and management of these places today.
Kluane: We respectfully acknowledge that the lands of Kluane National Park Reserve where this study took place are the traditional territories of the Southern Tutchone people represented in the Kluane region by the Champagne and Aishihik First Nations and the Kluane First Nation. Champagne and Aishihik First Nations, Kluane First Nation and Parks Canada are jointly responsible for the management of Kluane’s natural and cultural resources.
Bannf, Lake Louise, Yoho, Kootenay: We respectfully acknowledge that Lake Louise, Banff, Yoho and Kootenay National Parks are located within the traditional lands of the Ktunaxa and Secwépemc peoples. We recognize their stewardship of the lands and waters in the area since time immemorial. Parks Canada is committed to reconciliation and renewed relationships with Indigenous peoples, based on a recognition of rights, respect, cooperation and partnership.
Riding Mountain: We respectfully acknowledge Riding Mountain National Park is located within Treaty 2 Territory and Parks Canada works with First Nations from Treaties 2, 4 and 1. We mention this to honor our relations and the contributions made to establish the park, the province, and Canada.
Elk Island: We acknowledge that Elk Island National Park is located within Amiskwaciy (Beaver Hills), a landscape that has been cared for and used by Indigenous peoples for thousands of years. We recognize the deep cultural, spiritual, and ecological connections Indigenous communities maintain with this land and affirm the importance of respecting Indigenous rights, knowledge, and stewardship in the ongoing care of this place.
Prince Albert: We respectfully acknowledge that Prince Albert National Park is located on Treaty 6 Territory, the traditional lands of the Cree, Dene, Dakota, Saulteaux, and Ojibwe peoples, and is within the homeland of the Métis.
La Mauricie: We respectfully acknowledge that La Mauricie National Park is part of the traditional territory of the Atikamekw Nation. For thousands of years, the Saint-Maurice River watershed and surrounding lakes were travelled and cared for by Indigenous Peoples, including the Atikamekw, who continue to maintain a deep and enduring relationship with this land.
Pukaskwa:
Grasslands: We respectfully honour and acknowledge that Grasslands National Park is located on the traditional territory of the Niitsitapi, Nakoda, Cree, and Saulteaux nations, as well as the Métis Homeland. We recognize the deep historical and cultural connections these Indigenous peoples have to this land, which they have stewarded for generations, and we are committed to moving forward in the spirit of reconciliation, respect, and stewardship.
Introduction
Ecologists use prescribed fire as a deliberate and carefully planned management practice to mimic natural disturbance regimes, reduce accumulated fuels, and promote ecological processes that sustain biodiversity. In many forested and grassland ecosystems, periodic fire plays a critical role in shaping vegetation structure, influencing habitat availability, and maintaining species assemblages. Prescribed fires are therefore widely used to support conservation objectives, particularly in protected areas where maintaining ecological integrity is a primary mandate. Within Parks Canada, prescribed fire is a key component of wildland fire management. The Parks Canada Fire Management Program coordinates fire management planning and operations across lands within the Parks Canada network, including the design and implementation of prescribed fires. While prescribed fire is expected to provide long-term ecological benefits, short- and medium-term responses of wildlife, particularly migratory birds, can vary depending on fire timing, intensity, extent, and pre-fire habitat conditions. As a result, systematic monitoring is required to evaluate both potential impacts and benefits of fire management actions.
Prescribed fire has a long history as a management tool in Canada, used to reduce fuel loads, restore fire-adapted ecosystems, and guide vegetation succession toward desired ecological states (Weber and Taylor (1992)). Research across forested and grassland systems shows that prescribed burning can substantially alter understory composition, soil seed banks, and habitat structure, with effects that vary by fire timing, intensity, and pre-disturbance conditions (Faivre et al. (2016); Travnicek, Lym, and Prosser (2005)). Beyond its biophysical effects, prescribed fire also operates within complex social and management contexts, shaping how landscapes are accessed, monitored, and interpreted by practitioners and researchers (Sutherland (2019)). These interacting ecological and operational dimensions demonstrate the importance of targeted, long-term monitoring to evaluate wildlife responses, particularly for mobile taxa such as migratory birds, whose responses to fire may be transient, indirect, or delayed.
Over the past years, Parks Canada staff have conducted bird surveys within prescribed fire areas to assess how migratory bird communities respond to fire. This dataset provides an opportunity to examine changes in bird species composition, functional communities and richness before and after prescribed fire, as well as to evaluate variation across parks and fire events. The primary objective of this report is to describe and compare migratory bird communities before and after prescribed fire and to assess whether detectable differences are evident in species-level and community-level metrics. A secondary objective is to evaluate the sampling strategy used, including its capacity to characterize bird communities and its statistical power to detect or model changes associated with prescribed fire. The results of this work will support Parks Canada in interpreting the ecological outcomes of prescribed fire and in refining future monitoring and fire management practices.
Methods
Data collection
Data were collected using passive autonomous recording units (ARUs; Shonfield and Bayne (2017)) at 236 locations (Figure 1).
Recordings were summarized by month and park to characterize the temporal distribution of sampling effort across each region and year (Figure 2). Sampling occurred primarily during the breeding season (April–July), with the majority of effort concentrated in June. Temporal coverage varied among parks with sampling in single years and others sampled repeatedly across multiple years (Table 1). Monthly sample sizes ranged from low replication in some early or late-season deployments to higher levels of effort in core sampling months, reflecting differences in deployment timing, park-level logistics, and local program objectives.
Code
datatable(locs_summary,
options = list(
searching = TRUE,
paging = TRUE,
pageLength = 10
)) |>
formatStyle(columns = colnames(locs_summary),
backgroundColor = styleEqual(c("NA"), "lightgray")) Data management and processing
The principal goal for data processing was to describe the acoustic community of species heard at locations while choosing a large enough subset of recordings for analyses. Data processing was completed in WildTrax with tags made using count-removal (see Farnsworth et al. (2002), Sólymos et al. (2018)) where tags are only made at the time of first detection of each individual heard on the recordings (MacPhail et al. 2026 In Review). Multiple recording lengths were used in analysis relative to each national park (Figure 3) We also verified that all tags that were created were checked by a second observer (n = 59.83) to ensure accuracy of detections. If amphibians were detected, abundance was estimated at the time of first detection using the North American Amphibian Monitoring Program with abundance of species being estimated on the scale of “calling intensity index” (CI) of 1 - 3. Mammals such as Red Squirrel, were also noted on the recordings. After the data are processed in WildTrax, the wildrtrax package is use to download the data into a standard format prepared for analysis. The wt_download_report() function downloads the data directly to a R framework for easy manipulation (see wildrtrax APIs).
Analysis
We compared species richness between pre-burn and post-burn surveys at locations of interest (Figure 4). For each location, year, and sample type, richness was calculated as the number of unique species observed. Mean richness per location was then computed across years. Paired differences between pre- and post-burn richness were tested using a Wilcoxon signed-rank test. Data were visualized with lines connecting pre- and post-burn values for each location, points for individual observations, and median values highlighted to illustrate paired changes across burn events.
We examined how species abundance varied with time since fire and site characteristics (Figure 5). Data were filtered to include relevant columns and distinct observations of species and abundance. For visualization, total and mean species abundances were summarized by years since fire. To test for statistical effects, a generalized linear model (GLM) with a Poisson distribution was fit, using individual abundance as the response variable and years since fire, vegetation type, and fire frequency as explanatory variables. This approach allowed us to assess the influence of fire history and site characteristics on species abundance.
We used redundancy analysis (RDA) to examine how species composition varied with burn status and vegetation type (Figure 6). Data were restricted to locations of interest and to pre- and post-burn samples. For each location, year, sample type, and vegetation type, species abundances were calculated as the number of occurrences per species. Data were then converted to a site-by-species matrix, with missing species values filled as zero. RDA was performed using sample type and vegetation type as explanatory variables. Site scores and species scores were extracted for visualization. The resulting RDA plot displays sites colored by vegetation type, ellipses representing burn status, and species vectors illustrating associations with axes. This approach highlights shifts in community composition associated with fire and vegetation characteristics.
We calculated the Shannon diversity index for each surveyed location and year to quantify species diversity. Only vertebrate species (excluding mammals, amphibians, insects, abiotic, and unknown detections) were retained. For each recording, the maximum observed abundance per species was recorded. Data were then aggregated by location and year to produce a species-by-location abundance matrix. Shannon diversity was calculated from these totals using the diversity function, and results were combined with vegetation type metadata for each site. Diversity patterns across years were visualized and grouped by Park (Figure 7) Shannon diversity index was compared between pre-burn and post-burn surveys at each location. Only sites with known vegetation type and sample type were included (Figure 8).
Results
Some of these analyses are still a work-in-progress. Check back soon for updates and additional details.
A summary of locations and covariates collected can be found in Table 2. Species richness overall did not change significantly following burns. Paired pre- and post-burn comparisons across locations showed no systematic increase or decrease in richness (Wilcoxon signed-rank test, V = 19.5, p = 0.44).
Wilcoxon signed rank test with continuity correction
data: paired_df$post_richness and paired_df$pre_richness
V = 19.5, p-value = 0.4443
alternative hypothesis: true location shift is not equal to 0
Discussion
Lorem ipsum.
Recommendations
Location metadata
Consistent site naming is essential for linking repeat surveys across years and ensuring data loads cleanly into WildTrax. Location names should follow a stable, documented convention that persists across seasons and years (e.g., EINP-MOS-3, AC-C2). This avoids fragmentation caused by ad-hoc renaming and allows records from different sampling periods to be reliably joined. Where possible, names should use a hierarchical structure, moving from broad to specific (e.g., project or region → sub-site → point → treatment). This structure supports automated grouping, summarization, and analysis without manual reconciliation. Do not include burn type, year, or post-burn state in the location name, as this information is managed separately in metadata. Including it directly in the name can create inconsistencies within WildTrax as different location names would be generated. However, each location can still store spatial coordinates, and locations can be merged as needed within the system.

Recording frequency and duration
Prescribed burns typically occur outside the migratory breeding bird season; therefore, within-season pre- and post-burn surveys are generally not required. However, if a burn occurs during the breeding season (approximately May 15–July 15 across most surveyed parks), additional surveys should be conducted during both the breeding period and the subsequent post-breeding period (July 15–August 31) to capture potential immediate and short-term responses.
Both resident and migratory species may respond to changes in burn state. Migratory species, in particular, may temporarily avoid burned areas or pass through them, although empirical evidence remains limited. Long-term acoustic monitoring using stationary ARUs placed within or adjacent to burn areas provides an effective approach for assessing site use and changes in community composition over time.
Recording effort should be consistent across the park system to support comparability among sites and years. For single-visit point counts, a minimum of 11 minutes of recording per session is recommended, which includes buffer time for the observer to initiate and complete a 10-minute acoustic survey. For surveys using pre-programmed ARUs, we recommend a staggered recording schedule distributed across seasons (winter, spring, summer, and post-breeding) and times of day (Table 3).
Recordings are organized into blocks representing short (3–10 minute) sampling periods spread across the season and diurnal cycle. This staggered design balances logistical constraints with ecological coverage, providing a robust dataset for characterizing bird communities and evaluating responses before, during, and after prescribed fire events.
| blocks | Julian_Days | Recording_min | Daily_Times |
|---|---|---|---|
| 1 | 90-104 | 60 | 00:00, 02:00 |
| 2 | 105-119 | 60 | 05:37, 07:07 |
| 3 | 120-139 | 60 | 12:00, 15:00 |
| 4 | 140-149 | 180 | 21:10, 23:10 |
| 5 | 150-159 | 180 | 00:00, 02:00 |
| 6 | 160-169 | 180 | 05:37, 07:07 |
| 7 | 170-179 | 180 | 12:00, 15:00 |
| 8 | 180-210 | 60 | 21:10, 23:10 |
| 9 | 90-139 | 60 | 00:00, 02:00 |
| 10 | 140-159 | 60 | 05:37, 07:07 |
| 11 | 160-179 | 60 | 12:00, 15:00 |
| 12 | 180-210 | 60 | 21:10, 23:10 |
Burn severity
For a field-usable and repeatable approach to assessing burn severity, which is likely a strong indicator of bird community composition (Knaggs et al. (2020)), the Composite Burn Index (CBI) provides a well-established framework. CBI is designed to characterize fire effects across vegetation and substrate layers in a way that is comparable among sites and directly compatible with remote-sensed burn severity products.
Assessments are conducted within a defined plot, typically on the order of 10-30 m in diameter to align with common satellite pixel sizes. Observations are made across five vertical strata: (1) substrate (litter, duff, soil), (2) herbaceous and low shrub layers, (3) tall shrubs and saplings, (4) understory trees, and (5) dominant and co-dominant trees. Within each stratum, multiple attributes, such as vegetation mortality, soil exposure, and char height, are scored on a standardized 0–3 scale, where 0 represents unburned conditions and 3 indicates the highest severity. Scores are averaged within strata and then combined to produce a composite severity value for the plot. To ensure consistency and transparency, field assessments should be recorded on standardized scoring sheets and supported by photographs taken from fixed viewpoints. In practice, field sheets should capture plot identifiers that align with the site naming convention, GPS coordinates, date and time, and observer names, followed by structured scoring for each vegetation and substrate layer. Additional observations, such as maximum char height, the presence of unburned refugia, or evidence of fire skips, provide valuable context for interpretation. For operational use, composite scores can be grouped into simplified burn severity classes (unburned, low, moderate, high). These classes retain the ecological meaning of the underlying CBI scores while providing a clear, interpretable summary that can be readily linked to management objectives and remote sensing products.