Recent Crown Thinning in a Boreal Black Spruce Forest Does Not Reduce Spread Rate nor Total Fuel Consumption: Results from an Experimental Crown Fire in Alberta, Canada

A 3.6 ha experimental fire was conducted in a black spruce peatland forest that had undergone thinning the year prior. After 50 m of spread in a natural stand at 35–60 m min−1, the crown fire (43,000 kW m−1 intensity using Byram’s method) encountered the 50% stem removal treatment; spread rates in the treatment were 50–60 m min−1. Fuel consumption in the control (2.75 kg m−2) was comparable to the treatment (2.35 kg m−2). Proxy measurements of fire intensity using in-stand heat flux sensors as well as photogrammetric flame heights had detected intensity reductions to 30–40% of the control. Crown fuel load reductions (compensated by higher surface fuel load) appear to be the most significant contributor to the decline in intensity, despite drier surface fuels in the treatment. The burn depth of 5 cm in moss and organic soil did not differ between control and treatment. These observations point to the limited effectiveness (likely reductions in crown fire intensity but not spread rate) of stem removal in boreal black spruce fuel types with high stem density, low crown base height and high surface fuel load. The observed fire behaviour impacts differ from drier conifer forests across North America.


Introduction
Strategies for decreasing home and structure loss in the Wildland-Urban Interface (WUI) have been framed at a variety of scales including building materials [1], vegetation, and environmental characteristics immediately surrounding homes [2], as well as regional scales [3]. This increasing focus on firebrands (embers) as a primary mode of home ignition allows for a two-tiered approach to wildfire Fire 2020, 3, 28 2 of 18 risk mitigation: a simultaneous focus on homeowners tackling ignition sources on the home exterior and immediate surroundings [1], while wildfire management agencies and municipalities can reduce factors related to firebrand production and spread. Fuel treatments such as thinning and pruning are resource-intensive alternatives to surface prescribed burning in ecosystems with low crown base heights, high surface fuel loads, or dense understory ladder fuels [4]. A candidate ecosystem for thinning and pruning lies in the boreal forest, where the common spruce genus is not adapted to surface fires, and crown fire dominates [5]. In such stands, the goal of fuel treatments such as thinning and pruning is not the wholesale stopping of the fire, but rather an improved suppression environment characterised by lower fire intensity, lower rates of spread, and improved access for suppression resources such as ground crews of aircraft.
Increasing tree spacing alongside lower branch pruning can be expected to lower the likelihood of a surface fire accelerating into a crown fire [6]. The advantages gained from lower crown fuel loading and overall crown fuel reduction should outweigh the increased surface wind speeds, which act to marginally increase surface fire intensity [4]. Increasing tree spacing allows for increased ground-level evaporation and decreased surface fuel moisture via increases in solar radiation and wind [7]. At the same time, more rainfall reaches the surface, and root uptake of duff moisture by trees is decreased [8]. Specific to boreal conifer forests, shade-adapted common feather mosses can experience excess direct sun as a result of thinning operations that change the forest floor light regime. This photostress can lead to moss desiccation and reduced moisture content [9] with likely implications for increasing the depth or rate of duff consumption [10].
Thinning treatments are not expected to stop a wildfire, but rather reduce intensity and rate of spread so that suppression can be successful under conditions where a wildfire in a natural stand would be beyond suppression. Drier surface fuels are likely to occur in these more open thinned stands; as a result, a higher surface intensity has been observed compared to unthinned areas during moderate fire danger conditions [11]. Thinned stands are anticipated to be less flammable than natural stands during extreme wildfire danger when suppression success is unlikely. This reduction in crown fuel load directly leads to lower head fire intensity given the same rate of spread, resulting in increased initial attack success [12]. Lower fuel loads also assist in increasing fireline construction rates [13], allowing for faster mechanical line construction [14], and enhanced aerial suppression effectiveness [15]. Thinned stands are also more accessible for ground crew movement [16] and allow greater ground coverage if sprinklers are deployed [17].
Several case studies in boreal forests indicate that thinning can be effective. Two test burns in jack pine stands in northwestern Canada indicated that thinning and surface fuel cleanup could result in substantially slower and less intense fire behaviour compared to natural stands [18]. Similar tests have been attempted in black spruce stands and suggest similar outcomes as observed in thinned jack pine, albeit under less severe fire weather conditions [19,20]. A case study of a crowning wildfire burning into thinned black spruce indicated that thinning did enhance aerial and ground-based suppression efforts [21]. The impacts of thinning on fire behaviour in black spruce with its uniquely high surface fuel load and under extreme wildfire danger conditions remains unknown. In this study, an experimental fire in thinned black spruce is analysed to determine how the fire behaviour was affected by fuel treatment.

Site Description
The Pelican Mountain experimental fire research site is located in central Alberta, Canada (55.71 • N, 113.57 • W), in the southern boreal mixedwood forest of western Canada. The site itself is a 1.9 km 2 poorly drained peatland sitting at 640 m a.s.l ( Figure 1). The surrounding area is well-drained aspen forest with numerous headwater streams flowing north from the Pelican Mountain ridge which rises as high as 940 m a.s.l some 12 km to the south. The site is almost exclusively peat soil (organic soil greater than 40 cm thick) with a dominant black spruce (Picea mariana) cover and minor stands of larch (Larix laricina). Tree cover is generally dense, with canopy closure > 50%, and dominant tree heights between 5-10 m. Small trees and low branches (<1 m above the ground), even in dominant trees, provide continuous aerial fuels from the ground to tree top. Feather moss (Pleurozium schreberi (Brid.) Mitt.)) is the dominant ground cover, with minor Sphagnum spp. mosses and Cladonia spp. lichens in more open areas. The experimental fire that took place in May 2019 occurred in the northernmost portion of the Pelican Mountain site, where the peatland drains into natural beaver ponds with aspen and mixedwood forest beyond. The climate of the area is cool and continental, with a mean annual temperature of 1.7 • C and 460 mm of annual precipitation. The winters are dry (120 mm average snow water equivalent), with snowmelt typically occurring in early April. The month of April is typically dry (22 mm precipitation average, half as snow), with similarly limited rainfall in the first half of May. During this relatively dry period, deciduous vegetation (such as aspen forest, but also shrubs and grasses) are leafless, leading to high fire spread potential [22] and continuous crown fires in both upland forest as well as peatlands [23]. Fires in the region are 50% lightning-caused overall, with a median fire size (of fires over 50 ha) of 390 ha, and an overall annual burn rate of 0.5%, or a 200-year average fire return interval [24]. Human-caused spring fires in the region (defined as prior to deciduous forest leaf-out) are some of the largest and most impactful wildfires in Canada [25].
To test the effectiveness of tree removal thinning in black spruce, the Alberta Wildfire Management Branch selected the Pelican Mountain research site in 2015 as the location for multiple fuel and vegetation manipulation trials and experimental burns. The experimental block (Unit 5) burned in May of 2019 was thinned and pruned ( Figure 2) by hand crews over the winter of 2018/2019, and the perimeter was mulched by a mechanical drum mulcher in early 2018 ( Figure 1). Unlike mechanical thinning treatments using heavy machinery, the feather moss forest floor was not compacted or disturbed during the winter hand-thinning operation. All residue was piled and burned on-site, resulting in approximately 80 burned spots averaging 1.5 m in diameter.

Canopy Height and Fuel Loading
Measurements of fuel loading were conducted following the standard fixed radius plots in 2018 and spring 2019 in the treatment area before and after thinning. Pre-treatment measurements are assumed to be equal to those in the control.
Multispectral airborne LiDAR data for canopy height mapping and precise block layout was acquired over the Pelican Mountain Research site on 19 August 2018, by the University of Lethbridge ARTEMis lab. The westernmost portion of the hand-thinning treatment was not complete at the time, but most site preparation was complete at the time of the flight. The LiDAR survey was obtained using the Teledyne Optech Titan sensor mounted on a Piper Navajo aircraft. This sensor captures three wavelengths of 1550 nm, 1064 nm, and 532 nm. Preprocessing was also completed by the ARTEMis lab using LiDAR Mapping Suite (proprietary software from Teledyne Optech) to generate a calibrated point cloud for each channel. Accuracy for point cloud height data was completed using roof line analysis and had a delta height RMSE of 3.5 cm. For this research, the point clouds from each channel were merged into a single LAS dataset using US Forest Service FUSION/LDV 3.8 software [26]. FUISION/LDV was also used to create a 2.0 m resolution Digital Elevation Model (DEM) with the GridSurfaceCreate function. The CanopyModel tool was then used to create a 2.0 m resolution Canopy Height Model by using the highest return in each grid cell normalised with the DEM.
Monochromatic terrestrial LiDAR was collected using a GeoSLAM Zeb-Revo Handheld laser scanner. Scanning at 905 nm with an acquisition rate of 43,200 points/s and a maximum range of 30 m; before and after burn scans were collected from two control plots and seven test plots. The GeoSLAM collected point clouds with a relative accuracy of 2-3 cm. Data was collected in a multi-loop pattern to resurvey features as often as possible. This pattern resulted in minimised error and increased accuracy of the resulting point cloud. For this research, terrestrial LiDAR was implemented to assist with fuel consumption parameters. Point cloud data was preprocessed using GeoSLAM Hub software and further analysed in FUSION Cloudmetrics.
In the field, a nested plot model (see Supplementary Material for ground plot sampling methods) was used to measure all trees within a specified radius, depending on the diameter at breast height (DBH). Within the specified radius, all trees greater than 1.37 m in height were measured for height, live crown base height, and DBH. For trees with a DBH greater or equal to 9.0 cm, dead crown base height and canopy widths were also measured. For trees with a DBH less than 9.0 cm, a 3.57 m sampling radius was used unless there were fewer than 20 trees. In this case, a 5.64 m radius was used to provide a more representative sample of the trees in the area. For all trees greater than 1.37 m in height with a DBH greater or equal to 9.0 cm, a 11.28 m sampling radius was used. Canopy fuel load was calculated at the plot level using published allometric equations [27].
Ground cover was assessed visually using a 0.5 × 0.5 m vegetation quadrat placed on the ground along three 40 m transects within each treatment (Figure 1), with measurements every 2 m (for a total of 63 measurements per treatment). Moss cover was split into the two dominant moss types (genera) in black spruce dominated peatlands, Sphagnum and feather moss, with lichen (Cladonia/Cladina spp.) and litter also included as ground cover types. Moss cover assessments post-fire also accounted for moss burn severity with categories of lightly-and heavily-burned; whereas if moss capitula/moss structure is intact, but there is evidence of burning, the area is classified as lightly-burned (singed); if capitula/structure have been damaged considerably, the area is classified as heavily-burned (burned). Post-fire bare peat (where composition could not be determined) was also included. Depth of burn (DOB) was measured using 20 cm metal DOB pins placed in the centre of each vegetation survey plot. The top of the pin was made flush with the ground/moss surface, and so the height of the pin above the post-fire ground surface is equal to DOB. Surface volumetric water content (VWC; top 6 cm of moss and duff) was measured, pre-and post-fire (<10 h pre-fire and 1-2 days post-fire), using an ML3 ThetaProbe by Delta-T Devices (accurate to 1% VWC) calibrated using curves from [28]. VWC was then converted to gravimetric water content (GWC) using the bulk density of peat core samples taken in 0.1 m diameter PVC pipes of both Sphagnum and feather moss pre-fire. Depth to seasonal ground ice was measured at each vegetation plot (n = 63 each in treatment and control) using a thin metal probe, both 48 h prior to and within 48 h after the fire.

Weather and Fuel Moisture
Automatic weather stations have been located on site since 2016 to measure precipitation, air temperature, and humidity, though the forest opening is not sufficient to measure wind speeds. Overwinter precipitation was measured 32 km away at the Wabasca Airport (55.96 • N, 113.82 • W). On the day of the fire, an ATMOS-41 2D sonic anemometer and temperature/humidity probe was installed on a tripod 6.1 m above the ground (2 m above the surrounding trees) 40 m south of the ignition line ( Figure 1), with wind measurements taken every 10 s. Wind speeds were adjusted to 10 m standard heights using the logarithmic wind profile [29] with an apparent canopy height of 2.8 m, a zero-plane displacement height of 0.52 m, and a friction velocity of 0.93 m s −1 .
Windsond miniature weather balloons were launched from the southeast corner of the burn site at 14:50 (three hours prior to ignition) to measure the atmospheric conditions above the forest. Temperature, wind, and relative humidity were measured to determine turbulence and stability structures in the atmosphere. Additional measurements were made with a handheld air temperature, humidity, and pressure sensor in a helicopter recording airborne infrared measurements up to 700 m above ground level.
Destructive fuel moisture samples were collected between 13:00 and 15:00 on the day of ignition, including duff, feather moss tips (uppermost 2 cm), Sphagnum moss tips, black spruce needles, surface litter and surface mulch (to the north of the treatment). Each sample was weighed in the field and later dried in a convection oven for 24 h at 100 • C. Samples were collected in proximity to the ground-based sampling plots, and care was taken not to cause excessive compaction.

Ground Observations of Fire Spread and Intensity
In-fire video cameras, rate-of-spread data loggers, and heat flux sensors were placed around vegetation inventory plots or other observation points within the treated stand ( Figure 1). Video cameras were placed in passively cooled steel camera boxes [30] that could withstand the intensity of a crown fire. Ten video cameras located in the treated stand helped ascertain the speed, direction, intensity, and duration of the fire. The cameras were pointed towards the untreated area to capture the change in fire behaviour as the fire progressed from the untreated to the treated section of the plot.
A total of 31 temperature-triggered rate-of-spread data loggers (K-type thermocouple, sampling frequency 1 Hz) were laid out in a grid format in both untreated and treated parts of the stand and were used to obtain accurate time series fire progression data. The temperature probes from the rate-of-spread loggers were located approximately 30 cm above the surface. Five temperature-triggered aluminum heat flux sensors (radiometers) [31] (sampling frequency 1 Hz) were placed in both untreated and treated sections of the stand to capture the incident heat flux. The sensors were located approximately 50 cm above the surface.

Aerial Observations of Fire Spread and Intensity
Two helicopters (each nominally 500 m above the ground) hovered for the duration of the crown fire spread, with visible and infrared video taken to document fire behaviour. Infrared video was taken using a FLIR E40 handheld infrared camera with a dynamic temperature range of up to 0-600 • C. Infrared video was able to observe the location of the fire front through thick smoke. Both oblique visible and infrared images from selected video clips were georeferenced and orthorectified using the Monoplotting software [32] against the LiDAR digital elevation model and an orthophoto mosaic obtained from UAV flights in 2018. The fire front was located and digitised at 15 discrete timesteps (an average of once every 25 s) and digitised as a polyline with vertices every 10-20 m. A spatial database of the location of the fire front and associated time (in seconds) since the start of the ignition was used to build an interpolated surface of fire spread rate using a triangular interpolation technique in SAGA GIS to produce a 5 m resolution time of arrival grid. Half of the line vertices were removed, and the remainder converted to point locations prior to interpolation. GPS logger data from the helicopter ignition device was also used as point data in the interpolation. Additionally, flame height was measured using the same monoplotting software along the digitised fire front locations. Since the ground was not visible in the control, the flame base height in the control was sampled off the LiDAR canopy height raster at that location. For the treated thinned areas, the ground was visible, and flame heights were measured using the ground elevation as the height datum. Since the aircraft were observing the fire downwind, the flame angle and length were not quantifiable, but instead, a conservative estimate of flame length equal to flame height was made. Flame height was converted to a measure of fire intensity using the relation of [33] for crown fires:

Fuel Loading
The natural, unthinned (control) stand of 12,000 stems ha −1 of black spruce amounts to a crown fuel load of 1.3 kg m −2 of biomass, or a fuel load 60% greater than the standard C-2 boreal spruce fuel type crown fuel load of 0.8 kg m −2 . Litter and fine woody debris contributed an additional 0.39 kg m −2 to the surface fuel load, while coarse woody debris (> 1 cm diameter) averaged an additional 0.15 kg m −2 . Observed spatial variability in canopy structure varied widely, with the western edge of the treatment having taller canopy heights, greater fuel load, and lower canopy openness compared to the east ( Figure 1). Both the control and treatment ground cover were nearly identical, with between 86 and 88% cover of feather moss, and between 12 and 13% cover of Sphagnum mosses, in the control and treated areas, respectively.

Fuel Moisture
Litter and feather moss tips were similar in gravimetric moisture contents in the thinned stand (Figure 3a; 7 ± 2%), while feather moss tips were higher in moisture content in the control (17 ± 3%) compared to tips in the thinned or control block litter. However, volumetric moisture measurements of the entire upper 6 cm of feather moss in the control (5.8 ± 6.4%, n = 50) and treated areas (6.0 ± 4.5%, n = 45) were not different. Duff moisture measurements (below 6 cm) across all surface cover types were different, with values in the control (320 ± 100%) higher than those in the treatment (110 ± 50%).
One year or older black spruce foliage gravimetric moisture content was 79 ± 23%, showing considerable variability but on the lower range of the normal value used in the Canadian Forest Fire Behaviour Prediction System of 97% [34] (Figure 3).

Fire Weather and Atmosphere
Cumulative precipitation measured from 1 October 2018 until the ignition date of 11 May 2019 at the Wabasca weather station was 63 mm, or 37% of the 20-year normal of 169 mm of precipitation during the same period. Only 3 mm of precipitation was recorded at the onset of snowmelt around April 11 until May 11. Average air temperature during that same period was −6.3 • C, very similar to the 20-year normal of −6.1 • C. The Canadian Fire Weather Index System Drought Code (which correlates well with the water table in boreal peatlands, [35]) was 272 on 1 October 2018 at the nearby Sandy Tower, and reached 331 by May 11. A Drought Code value of 331 represents the 4th driest Drought Code on May 11 recorded for the 36-year record at Wabasca, corresponding to a 1 in 9-year dryness event using an empirical plotting position exceedance methodology [36]. Similarly, the Buildup Index (a harmonic mean of the Drought Code and Duff Moisture Code) was 71, the second highest on record and only 10 points lower than the May 11 2016 values, corresponding to the dry conditions causing extensive black spruce forested peatland burning during the very large 2016 Fort McMurray wildfire some 175 km to the northeast [37,38].
Synoptic conditions on May 11 were dominated by an upper air ridge that was present over the area with accompanied subsidence aloft. The subsidence allowed for clear conditions and, as a result, nighttime cooling reached a minimum air temperature of 2.8 • C at 06:00 local time on the morning of May 11, but relative humidity reached only a maximum of 76 % at the same time, suggesting limited nighttime recovery of fine fuel moisture. Clear skies and strong surface heating increased the air temperature to 20 • C by 13:00 local time with a relative humidity of 22 %, and light SE winds at 5 km with gusts to 22 km h −1 (Figure 4).
The 13:00 local time Canadian Fire Weather Index values were, daily, Fine Fuel Moisture Code 93.6, Buildup Index 69.6, and Fire Weather Index 27. Observed gustiness was driven by convective turbulence that was present during the latter half of the day, and during the time of the burn, due to the combination of a warm surface and warm air (Figure 4). Convective turbulence created a well-mixed boundary layer (the lowermost portion of the atmosphere), attributing to large variance in the mean wind speed at the surface. The top of the boundary layer was at roughly 2400 m above ground level, above which a weak inversion inhibited the rise of the convective thermals ( Figure 4). In the 30 minutes prior to ignition at 17:49, the 10-m wind speed was 10.8 km h −1 . During that 30 min period, only two 10 s gusts of 23 km h −1 were recorded. With the lighter observed winds of 5 km h −1 at 13:00 on May 11, the calculated Initial Spread Index of 8 increased to 18 if using gust wind speeds. This transient increase in winds increased the expected rate of spread in boreal spruce (C-2) fuels from 10 m min −1 to 28 m min −1 . By the time of ignition, air temperature was 22.2 • C and relative humidity 18.7%.   The southern polygon containing the ignition points is the control stand, and all areas past the east-west line 50 m into the control were treated. The pre-existing cutline running NW through the eastern treated area is also shown. Numbers correspond to points of known fire arrival time (seconds after aerial ignition commenced). Note that a fully continuous colour ramp from [39] is used for this and all figures, with mid-point colours of ranges shown in the legend for clarity.

Fuel Consumption and Fire Intensity
Canopy consumption was observed as nearly 100% in all areas of head and flank, but was not directly quantified. Surface depth of burn was measured immediately after the safe passage of the crown fire (within 30 minutes of ignition) when only patchy residual smouldering was occurring. Average depth of burn for the control was 5.9 ± 4.2 cm, and in the thinned treatment 5.4 ± 4.0 cm (Figure 6a). Nearly 100% of feather moss burned to some extent, while only 50% of Sphagnum moss areas had any combustion whatsoever. Feather moss-covered areas had the highest mean depth of burn of 6.8 ± 3.9 cm, followed by lichen-covered areas with an average of 4.9 ± 1.8 cm. Sphagnum-covered areas had a depth of burn of only 0.51 ± 0.94 cm, where Sphagnum in situ near-surface volumetric water content was 9.0 ± 6.8%, and feather moss was 5.9 ± 5.5%. Given an average bulk density for near-surface feather moss (the far dominant surface cover and largest depth of burn) of 18 kg m −3 [40], the mean total ground layer consumption was 1.06 kg m −2 of total biomass in the control and 0.97 kg m −2 in the treated area. Fifth percentile return heights from terrestrial LiDAR show trends similar to standard depth of burn measurements (Figure 6a,b). A qualitative relationship between depth of burn and change in return heights can be drawn as return height elevations range from 0 to 25 cm. In agreement with depth of burn metrics, post-burn scans show an overall decrease in most ground cover fuel elevation. Depth to the concrete layer of seasonal ground ice was variable between ground cover types, but not significantly different between the control and treatment (mean = 26 ± 8.2 and 26.5 ± 11 cm, respectively). Sphagnum groundcover was resistant to combustion without regard to depth to ice whereas depth of burn in feather moss tended to weakly increase as depth to ice increased (Figure 7).
Head fire intensity is measured typically using fuel consumption in both the canopy and all surface and ground fuels, including feather moss [41]. Total fuel consumption alongside the median gridded rate of spread from Figure 5 resulted in calculated head fire intensities (east and west portions only) of 43,000 kW m −1 in the control and 41,000 kW m −1 in the treatment, meaning that the treatment intensity was 94% of that in the control (Table 1). Incorporating only the canopy fuels, litter, and fine surface woody debris and using the same rate of spread values, the intensity in the control was calculated as 27,000 and 90% of the control value in the treated stand, or 24,000 kW m −1 . Median flank fire intensities calculated by Byram's method were between 6000 kW m −1 in the treatment and 14,000 m −1 in the control.  Median observed flame heights from the oblique aerial video (Figure 8) in the control area (n = 164) were 11.4 m tall, fully 6.8 m taller than the median tree height of 4.6 m in the control. Using the relation between flame length (conservatively estimating flame length equal to height) and intensity from [33], a median fire intensity of 16,500 kW m −1 was observed across the entire control area. Similarly, in the treatment, a median flame height of 5.5 m and a median intensity of 5500 kW m −1 was observed. Overall, the treated area was observed to have flame heights only 34% of those in the control (Table 1). Smaller flame heights in the mulch area showing head fire and flanking fire behaviour of 3.1 m and 1.9 m were observed, respectively. Alternate fully independent methods of computing the relative fire intensity between the control and treatment are given in Table 1. Though these methods do not directly yield the same units as Byram's head fire intensity, the relative change between control and treatment is comparable. In-fire totalizing heat flux sensor (recording both convection and radiation) recorded a total heat flux at 30 cm height of 3.9 MJ m −2 , and 29% of that value, or 1.1 MJ m −2 , in the treatment. Overall, a calculation of Byram's head fire intensity where only foliage is considered (with the treatment having 52% of its canopy fuels removed by thinning) showed the greatest similarity to the independently measured relative intensity values of 34% and 29% using the observed flame height and in-stand heat flux sensors, respectively.

Contrasts with Other Experimental Fires
Two test burns done in northwestern Canada in jack pine/black spruce stands that were also thinned resulted in a significant reduction of intensity and spread rate [18]. The natural stand in [18] and burn conditions were similar to those of the International Crown Fire Modelling Experiment (ICFME) test burns [27,41] and are as flammable as the Pelican Mountain fuels. Head fire intensities at ICFME exceeded 90,000 kW m −1 [41], with similar canopy fuel loads between the control plot and ICFME (1.3 vs. 1.6 kg m −2 ). The ICFME plots featured 5 cm of compact needle litter and a higher woody debris load in contrast to the 1 kg m −2 of consumed moss and duff and minor amounts of woody debris. The thinning work for the NWT burns was done during the summer, and full trees were skidded off site before being burned. As a result, the surface vegetation was heavily compacted, with sparse dead and down material either removed or lying flat on the surface. Comparatively, the surface fuels in this study were not treated or affected by tree removal since treatment was conducted manually and in winter. The importance of surface fuel treatment has also been shown in dry conifer work comparing fire behaviour between thinned stands with and without subsequent prescribed burns [42,43].
The comparable experimental fire in thinned black spruce in Alaska [19] is notable for the similar ecosystem, fuel loads, stem removal treatment, and control stand fire intensity over 30,000 kW m −1 . The large reduction in fire intensity down to a surface fire (160 kW m −1 ) observed in [19] may be in part due to the higher moisture content in the important upper live moss (36% vs. 7% in this study) and higher elevated fine fuel moisture (relative humidity 47% vs. 19% in this study). Moreover, greater fuel loads of live green vegetation (i.e., grasses and herbs) in the understory were likely present in [19] given the mid-June date of burn compared to mid-May in this study. Indeed, given the similarities in fuels and treatments, these two studies may represent both sides of a critical threshold of surface fuel moisture, below which this study was able to observe little impact of a fuel treatment compared to a control. Further examinations of fuel moisture and understory phenology effects need to be conducted in order to better understand the timing and duration of periods when thinning treatments in boreal spruce are expected to be of limited effectiveness.
The crown bulk density of the control stand (0.28 kg m −3 ; Table 1) is in the mid-range of observed experimental fire crown bulk densities in active crown fires spreading on the order of 50 m min −1 [44]. The treatment crown bulk density of 0.10 kg m −3 is lower than any other experimental fire value for continuous crown fires with such high spread rates. Video documentation of the fire in the thinned stand shows that despite the high crown fraction burned, the continuity of flaming combustion and the base of the flames was observed in the surface fuels, and thus should not be considered an active crown fire.

Fuel Moisture Considerations in Thinned Stands
In fuel beds dominated by dead forest floor fuels, the litter moisture content is modified by the density of the tree canopy, with more open stands seeing drier surface fuels [7], in this case limited to feather mosses. Process-based models of understory dead fuel moisture point to both increased solar radiation and surface wind speeds as drivers of fuel drying [45]. The dominant surface cover in this study was a mix of Sphagnum mosses and feather mosses, both of which are living moss layers rather than a thin layer of dead material. Feather moss is adapted to low light conditions in more closed canopy forests and maintains a lower moisture content compared to Sphagnum mosses [9]. This consistently lower moisture content is far more sensitive to drought in comparison to Sphagnum mosses, and leads to depths of burn that are over ten times greater in feather moss compared to Sphagnum, supporting findings from free-burning wildfires in similar stands [46,47]. This inability of the dominant surface cover of feather moss to resist the greater evaporation rates in the canopy likely lead to the lower surface moss moisture contents in the treated stands. While this lower moisture content at the surface did not ultimately lead to greater depth burn (Figure 6), this drier surface likely provided a modest contribution to the surface fire intensity, as suggested by both the Canadian empirical Fire Behaviour Prediction System (FBP) [48] as well as the lab-based Rothermal model [49]. The contrast in measured dead surface litter moisture content (not including any feather moss, Figure 3a) (10% vs. 7%) measured pre-fire in this study would ideally yield a smaller contrast of 700 kW m −1 vs. 1900 kW m −1 following the Canadian FBP, roughly consistent with observed intensity trends in this study. In the Canadian FBP System, the reduction of surface feather moss moisture content from 18% in the control to 6% as observed in the treatment (Figure 3b), if applied to the spread model of dead surface fuels in the understory of mature white or red pine (C-5) would see a contrasting fire intensity of a low 70 kW m −1 (control) to over 2400 kW m −1 (treatment). Similarly, a contrast in surface fire rate of spread would be expected from 0.24 m min −1 to 8 m min −1 , where none such occurred. Therefore, the outsized contrast in feather moss moisture content relative to the smaller differences in observed fire intensity suggests that increased feather moss consumption rates in the treatment did not translate to a faster rate of spread ( Figure 5) but merely a modest increase in fire intensity at the surface, which partially compensated for a lower canopy fuel load (Table 1).

Seasonally Frozen Organic Soils
The depth of the concrete layer of seasonally frozen organic soil (depth to ice) varied considerably throughout the stand as a result of differences in surface elevation of peatland microtopography, but was in line with other early season measurements of depth to ice in other non-permafrost boreal peatlands [50]. Although seasonally frozen soil can be an input of water into the upper layers of peat [50] and/or act to limit drainage [51], there is little field observation of the persistence of this added moisture in a sub-humid environment such as the Boreal Plains, especially in periods of high fire weather where evaporation rates are towards their upper limits. In this instance, there is the potential for the layer of frozen soil to act as a barrier to an upward movement of (unfrozen) water from the water table, leading to enhanced drying of the uppermost surface layers. Whilst Sphagnum mosses limited depth of burn under the full range of depth to ice, depth of burn in feather moss dominated plots generally increased as depth to ice increased, suggesting the ice layer may act to limit the propagation of vertical burning (Figure 7).

Evidence for Fire Intensity Changes
The two independent observation methods that showed a reduction in fire intensity (in-fire heat flux sensors and photogrammetric flame lengths, Figure 8) points to meaningful intensity reductions of 50-70% (30-50% of the untreated value) due to the tree thinning treatment applied the year prior. The photogrammetric flame lengths are not a systematic sample, however, and only represent snapshots in time during which aerial video was able to clearly observe flames during the helicopter's orbit around the fire. However, this intensity reduction still resulted in an Intensity Class 6 (>10,000 kW m −1 , as calculated using Byram's equation) active crown fire with nearly complete consumption of the thinned crown. Fires, even of this reduced intensity in natural spruce stands, are at the upper limit of consistently effective aerial suppression at the fire's head, while fires at or beyond the >17,000 kW m −1 range are beyond direct aerial suppression by water or retardant [52]. Conceptually, air tankers would deliver more water or retardant to the surface in thinned stands, analogous to enhanced rainfall, reaching the surface in more open forests [53].

Inferring Rate of Spread from Combined Aerial and Ground Data
Fire perimeter observations (infrared and oblique imagery) and point-based observations (thermocouples) of arrival time were combined into a discontinuous grid surface prior to interpolation. One complication with such inputs is that the line and point observations have inherently different levels of spatial and temporal accuracy. Historical fire behaviour sampling designs employed arrival times along a transect or sparse grid to estimate spread rate and acceleration on experimental fires [54,55]. In contrast, observations on free-burning wildfires often consisted of several points indicating head fire position at separate moments in time [56]. Grids permit more accurate and representative estimates of rate of spread than a transect [57], and increased sampling allows for a characterization of spread rate variability within the burn plot [58]. A study of flame front measurement observed arrival times on a fine spatio-temporal scale (10-cm pixels, 10 Hz) using an outdoor burn platform, and suggested that the various calculation methods produced significantly different estimates of the spread rate, in part due to differences in the estimated direction of spread vectors [59]. With the primary focus of this study being within plot variability of spread rate, the choice of calculation method is likely less important.

Conclusions
An experimental crown fire was ignited in a natural lowland boreal spruce forest in Alberta, Canada, in May 2019. The continuous crown fire burned into a thinned spruce stand with little alteration in spread rate or total fuel consumption. Spread rates in the treatment were similar to a natural boreal spruce fuel model (C-2) given the fuel moisture and wind speed including gusts. Despite drier surface fuels in the treated area due to increased sun and wind exposure, the final depth of consumption of the ground layer of mosses remained unchanged and was dominated by the consumption of feather moss. Any increase in surface fire intensity due to drier surface fuels did not compensate for the decrease in canopy fuel load (and consumption) from the thinning, leading to a net decrease in fire intensity as measured by flame heights and in-stand heat fluxes. Thinning work was done in winter and resulted in an uncompacted surface fuel bed that was vulnerable to drying and able to sustain the rate of spread. While fire spread rate changed little, the marked decrease in intensity and increased canopy openness from the thinning treatment suggests that direct suppression (i.e., sprinklers, mechanical fire guard, or aerial suppression with water, foam, or retardant) may be enhanced even under dry and windy conditions. Future studies should focus more on the role of active suppression measures and their interaction with fuel treatments, rather than fuel treatments as fully passive fire suppression strategies. New methodologies need to be developed to quantitatively integrate fuel consumption, fire energy, and rate of spread observations across the spectrum of flaming and smouldering. Additional work should be done on the potential for economically reducing the flammability of surface fuels in previously thinned black spruce stands.