1. Introduction
Climate and weather impact the frequency and behavior of wildfires [
1,
2]. Wildfires also provide substantial short and long-term feedbacks to the atmosphere and land-surface [
3]. The burning of biomass increases the output of emission products, heat and water vapor fluxes into the atmosphere; the short-term impact on the atmosphere is an immediate small-scale fire caused by local heat and moisture fluxes such as fire winds and pyro-cumulus clouds. Long-term atmospheric cooling can be caused by indirect and direct atmospheric aerosol radiative scattering or warming caused by released greenhouse gases such as carbon monoxide, methane and nitrous oxide [
4,
5,
6,
7]. Every year, global fires directly emit 200−300 Gkg(
g) of carbon into the atmosphere [
7]. Emissions from Canadian boreal biomass burning have regularly been observed in the upper troposphere and even in the lower stratosphere over Europe and the USA [
8]. Our capacity to understand and manage wildfires will remain imperfect as long as our knowledge is confined by our limited observation and restricted understanding of wildfire mechanisms. One of the areas that requires investigation is the transition of a wildfire to a large fire with a high intensity, spread and burned area. The current study offers a numerical simulation case-study and discussion of the dynamics and physics of the transition of wildfire to pyroconvection. The research also suggests that a coupled atmosphere–fire prediction system could provide additional guidance for decision-makers.
Certain atmospheric conditions can result in excessive fuel drying and extreme pyroconvection activities [
9]. Moisture in the mid-atmosphere tends to be resilient to rapid changes in the absence of dry/moist large-scale advection [
10,
11], while wildfire can change the atmosphere moisture content rapidly [
12]. A lack of wildfire–atmosphere heat and moisture interaction can throw off any non-coupled weather or wildfire forecasts.
The wildfires associated with pyroconvection behavior cannot be predicted by a simple two-dimensional physical or semi-physical spread model. Although the simplified fire spread models are a useful measure of fire front spread on the ground in the absence of the forceful and rapid changes in the wind and moisture of the atmosphere boundary layer, they cannot foresee the time and location of occurrence of pyroconvection and extreme forces.
Primary knowledge [
13] of mesoscale convection systems (MCSs) suggests three requirements for the development of deep convections: (1) instability, (2) an abundance of boundary layer moisture and (3) a lifting mechanism (trigger) to produce vertical motion, allowing a parcel to attain free convection. The requirements are fundamental to the development of all MCSs (or pyroconvection mesoscale convection systems, PMCSs) from short-lived single cells to destructive supercells. What separates the environments of varying storm types from each other is the vertical wind shear within the pre-storm environment. Buoyancy primarily determines the updraft strength, and wind shear determines the storm evolution [
13]. Wildfire surface heating can act as a trigger, allowing free convection and adding moisture (via fuel combustion) to the air in a short period of time, which will increase buoyancy. Wildfires can modify the pre-existing environment toward a favorable deep convection condition by warming and moistening the lowest level of the stratosphere.
The pyroconvection grows with very strong winds which change directions and drive with lateral vorticity and vertical wind shear. Badlan et al. [
14] performed an idealized numerical simulation to study the role of deep flaming in pyroconvection by using WRF-Fire version 3.9. They repeated the simulation with variations of fire intensities, spatial extensions and horizontal resolutions and three different atmospheric profiles to explore the dynamics of pyroconvection. These simplified simulations gave an informative insight into the dynamics of pyroconvection without considering the moisture fluxes. Lareau and Clements [
12] showed the atmospheric vertical thermodynamics and dynamics profile, and the daily variation had a significant impact on the onset and development of pyroconvection, as the wildfire moisture and heat fluxes from below were shown to be able to expedite the process. They established direct observations using a scanning Doppler light detection and ranging (LiDAR) and mobile radiosonde system during two large wildfires in northern California. The results suggested that the plume condensation level was substantially higher than the ambient lifting condensation level (LCL), which implies that the lifting air in the fire plume reflects the actual properties of the lower plume. The plume must be cooled largely near its base due to entrainment as it decays towards adiabatic ascent higher up the convective condensation level (CCL).
To examine the transition process of a wildfire to a pyroconvection, we chose to study the early stages of the Horse River wildfire which eventually burned an area of 589,522 hectares. The emissions from the Horse River fire had a pronounced impact on the measured annual concentrations and health monitoring sites (NH
3 was up to 14 times higher than non-fire impacted periods for some pollutants) [
6].
We are interested in the first days of the event because that is when the fire made the first transition to pyroconvection. We simulated the fire in hindcast mode using the numerical weather and fire prediction model. The choice of model depends on the questions at hand and the availability of the model and resources. After some investigations, we used WRF-Fire [
15].
The objective of this paper is to investigate the transition of a wildfire to a large wildfire and the initiation of pyroconvection. We also investigate the practicality of predicting a large fire. We study the transition mechanism by using WRF 3.9.1. We present a Horse River wildfire (Fort McMurray, Alberta Canada) hindcast simulation, as well as our results, with a focus on atmosphere profiles and storm onset. The current research is not an observational or numerical advancement to wildfire models; it is rather an application and a study on the practicality and efficacy of wildfire forecasts. Lastly, we discuss the existing challenges to a fully functional operational coupled atmospheric–wildfire forecast system.
2. Horse River Wildfire 2016
The 2015–16 El Niño (NOAA El Niño Advisories) caused a warmer winter over Canada, at almost 5 C warmer than the climatological average across all provinces. The 2016 spring remained warmer and drier than the normal climate in western Canada and the Prairies, while an inflow of cold Arctic air brought colder than normal temperatures over eastern and northeastern Canada. The average precipitation for the period 1981–2010 for Fort McMurray in July was 80.7 mm and the maximum temperature was 24 C. Fort McMurray international airport weather station recorded winter temperatures which were about 4 C warmer than normal in 2015–2016, while precipitation was around half that of the normal records. This led to an extremely dry biomass condition combined with gusty winds at the beginning of May, creating critical conditions for wildfire with extreme fire danger indices in the area on May 1, 2016. The subject area is located in the northeast of Alberta, south-west of Fort McMurray. The city is in the middle of the Athabasca oil sands, surrounded by boreal forest of white spruce, jack pine, trembling aspen and peatlands. The overall topography of the area is relatively flat, while the four rivers that pass through the city of Fort McMurray create slopes with riparian vegetation of white spruce and balsam fir.
A severe wildfire in May 2016 led to the largest evacuation in the history of Canada and caused widespread damage. There were five active fires between April 29 and May 1, 2016 in the area. Four of them were within close distances from the populated area of the city of Fort McMurray; those fires (MMD001-004) were suppressed and controlled by Alberta Fire crews less than a day after their ignitions. The Horse River fire (identified by MWF009-16) started in the south of the city and swept through the city in less than 72 h. It destroyed over 2400 homes and buildings. It continued to spread across Alberta and into Saskatchewan. The Horse River fire started on May 1, 2016, and a mandatory evacuation was issued on May 3. The fire spread reached 30 to 40 m per minute and created fire-driven weather patterns and winds at its peak. The satellite image from the second day of the fire displayed the pyro-cumulus clouds. The burned area stretched west and east towards the Saskatchewan border, covering 589,552 hectares. The wildfire was declared to be under control on 4 July 2016 after burning for 65 days. During the fire, the lack of knowledge regarding the wind direction (due to the coarse observation and weather forecast biases) made it difficult to predict the fire behavior. The forecast uncertainty was large, which was partially due to fire wind modification and associated heat fluxes. The scientific questions regarding the fire are the following: What was the main trigger that unleashed the Horse River fire (MWF009-16)? What was the main difference between the Horse River fire and the other four fires in the area that had been suppressed? and could we have predicted the Horse River wildfire in a timely manner? To answer these questions, we examine the weather and fire as one system.
4. Results
The Fort McMurray event was the costliest natural disaster in Canada and a fitting case to investigate the viability and serviceability of using a simulation to save lives and costs associated with a wildfire. However, this was also a difficult case to study due to the limited data, which is a common issue in extreme events. During this study, the only available data included the satellite images, hotspots and total daily burned area images. The WRF-Fire simulation was set up for five domains (
Figure 1) for a period (lead time) of 24 h. The 24 h hindcast simulation took 6 h to complete using 144 computational processing units (CPUs). The reasonably fast performance shows the functionality of the model for real-time forecasts. The simulation started a few hours before the fire, early on May 3. The final stage of simulation at the time of pyroconvection onset is presented here.
Figure 2 is an illustration of hotspots, reported fires and vertical profile locations. The vertical profile location is in the vicinity of the fire front on May 3. The model wind speed output is illustrated in
Figure 3. Wind speed increased in the burned area and reached its maximum where the heat flux was at maximum, as illustrated in
Figure 3b. The convergence and divergence became visible and more pronounced in the finer-resolution inner domains.
Figure 3a demonstrates the 1 km resolution wind at 10 m above-ground at 00 UTC on May 3, 2016.
Figure 3b simultaneously shows the 200 m resolution wind speed and direction. The comparison of the two images shows that the wind speed at the burned area was up to 3 times faster than the surrounding environmental wind speed.
Figure 4 is a 3D illustration of the vertical wind speed (updraft). The north–south view of Fort McMurray city demonstrates many up and down draft cells in the burned area. The strongest updraft is presented in red, and downdrafts are in the blue spectrum. The vertical velocities between −1 to 1 m/s were eliminated by the reduction of color transparency in this figure, allowing a better presentation of severe pyroconvection updrafts. Rising air occurs in the region of maximum vorticity, and sinking air is associated with negative vorticities. The next figures shows variational views at the location of the strongest model output for updraft.
Figure 5 is a thermo-diagram (Skew-T diagram) of model outputs at the location of 56.69
N and 111.39
W (a few kilometers to the north of the burned area). The thermo-diagram shows a well-established environmental convective available potential energy (CAPE) by 00 UTC on May 3.
Figure 6 shows two thermo-diagrams at the same location but at a different time to show the atmosphere profile changes due to daily warming and synoptic forces. The simulation demonstrated a rapid change in environmental CAPE starting at 18 UTC on May 2. The CAPE reached 887 J/kg in less than 5 h (
Figure 5 and
Figure 6b). The unstable atmosphere was caused by the diurnal heating and moistening of the lower planetary boundary layer (PBL) While the nocturnal thermal inversion disappeared at around 18 UTC (
Figure 6a), the dew point temperature increased at the lower PBL via synoptic moisture advection. The hourly observation from YMM (Fort McMurray airport) demonstrated the same trend of temperature and moisture increasing with slightly lower dew-point temperatures as model outputs.
The wind direction changed around 18 UTC from south westerlies to south easterlies in the burned area. The high bulk Richardson number of 93 (BRN) at 23 UTC (
Figure 5) was an indicator of the small vertical wind shear and high CAPE; it was also an indicator of strong buoyant convection. The mid-level lifted index (LI, which is an index used to assess low level parcel stability) changed from positive (
) to negative (
) during the same period of time. A negative LI indicates that the PBL is unstable with respect to the middle troposphere. An LI magnitude of
to
is classified as a large, unstable BPL. This is an environment in which convection/pyroconvection can occur. The more negative the LI, the more unstable the atmosphere will be for a rising parcel of air from the PBL. At the same time, the lifting condensation level (LCL) and level of free convection (LFC) merge at the 800 mb pressure level (2–2.5 km above the ground) with an equilibrium level (EL) about 270 mb or 11 km above ground level. Our simulation indicated that the pyro-cumulus cell onset was as late as 21 UTC on May 2, when the CAPE reached 650 J kg
, and the convection inhibition (CIN, which is the energy that is needed to lift an air parcel upward adiabatically to the LCL) was 0 J kg
. The dew-point temperature of the atmosphere profile above 50 kPa increased to a maximum of 10
C around 8 km above-ground in less than 13 h (
Figure 6a,b, green line). Moistening the entire atmospheric profile increased the atmospheric buoyancy.
Figure 7 gives a demonstration of three cross-sections at the grid points [170, 186], where the highest heat flux and wind speed occurred. The figure includes three snapshots at 00 UTC on May 3, 2016 in three different vertical planes at 163, 0 and 90
from north at the the grid point with the highest wind speed.
Figure 7a is identical to a classical MCS cross-section with two downdraft areas behind and in front of the main updraft cell. The extension of the updraft 5 km above ground level shows the extension of pyro-cumulus in the back and front of the main cell. The two other cross-sections (
Figure 7b,c) are the meridional and zonal views of the cell at 00 UTC on May 3, 2016. The pyroconvection activity slowed above the freezing level, but updrafts stretched to the top of the troposphere (
Figure 7b,c). The vertical updraft reached the tropopause as early as 6 h after the fire ignition in the model (not shown here).
Figure 8 and
Figure 9 are two close-ups of the horizontal planes at model levels 1 (second level of 60 model levels) and 10 (boundary layer).
Figure 8 displays a detailed perspective of horizontal wind, temperature and humidity slightly above-ground within the center of the pyroconvection updrafts and downdrafts.
Figure 9 is similar to
Figure 8 at model level 10, located around 2.5–3 km above the ground. The comparison of the vertical and horizontal views indicates the developed cells, associated downdrafts, dry air, wind convergence and divergence locations. The thermal gradient aloft reached 2
C at 3.5 km above the ground (
Figure 9b). The thermal gradient created strong convergence and divergence areas, strong surface winds, updrafts and downdrafts (
Figure 8a and
Figure 9a). The updraft reached a magnitude above 11 m/s at around 4 km of altitude (
Figure 7c), which is well above the LFC (
Figure 5). The unstable PBL with the fire wind created large cells of pyrocumulus. The downdrafts of the main cell created smaller cells in the direction of the lower-level (0–6 km above the ground) vertical mean wind. On the other hand, the heat fluxes from the fire front created a horizontal thermal gradient as large as 10
C in a distance of less than 200 m (
Figure 8b) at the surface.
Figure 10 demonstrates the burned area on May 1 and 2.
Figure 10a presents the WRF-Fire burned area outputs at 00 UTC on May 3. The overlay of the May 2, hotspots and WRF-Fire burned areas is acceptably matched. WRF-Fire outputs gave a demonstration of the rate of spread, the fire perimeter as wind direction changed and pyroconvection activities. They captured the change of wind direction and timing when the eastward fire front slowed down and the fire progressed westwards towards the river. The simulated fire front shifted almost 2 km west of the actual fire perimeter in this simulation on the west and east sides. The bias is mainly due to the ignition locations, time and biomass map resolution, which is discussed in
Appendix A.
The NASA image on May 3, 2016 (
Figure 11) demonstrates the large pyrocumulus cloud over the area with a distinguished rotation and overshooting area.
Figure 12 shows two different views of the Horse River wildfire at early development time (23:51 UTC May 1). The updraft and downdraft can be identified at early stages. There is no report on lightning during this fire (a strong downdraft can produce more fires via lightning).
5. Conclusions
Despite the shortcomings of initialization (hotspot acquisitions) and high-resolution vegetation and biomass fuel data, we were successful in simulating the pyroconvection activity with good accuracy in a timely manner. The 24 h simulation was completed in 6 h using 144 computational processing units (CPUs) from the resources of Compute Canada. We captured the increased horizontal wind speed and updraft, moisture fluxes, strong thermal gradient, updraft and onset of pyroconvection. None of the above could be captured with an atmosphere or wildfire model. A coupled atmosphere–wildfire model is required.
The use of the vertical atmosphere profile is crucial to foresee the onset of pyroconvection, unstable atmosphere and buoyancy. Model outputs can be a suitable substitute for the atmospheric profile around a burning area where there is no sounding. The simulation demonstrated that the combination of the day-time elimination of nocturnal capping inversion, Horse River wildfire pyroconvection (increasing heat and moisture at the surface) and moisture advection increased the CAPE. The 887 J/kg total increases of CAPE in less than 5 h and the high bulk Richardson number of 93 (BRN) were indicators of the growing pyro-cumulus cell. The simulation showed the adequacy of the current WRF-Fire approach for real-time applications.
We have made an effort to overcome all compatibility, technicality and mismatch issues between different data sources and the WRF-Fire model, such as adapting the Canadian biomass fuel classification for WRF-Fire. However, there is a larger stumbling block to overcome on the road to an operational coupled model. Wildfires and associated problems are classified as ground-based issues, which are managed by provincial, territorial or state levels of governments. Unlike fire, the atmosphere is a concern for national and international organizations and management levels. The separation of authorities and their abilities inti the two key fields is an obstacle on the road to an operational coupled atmosphere–wildfire model that requires actions and collaborations from outside of scientific communities. While it sounds simple to add a fire model to a weather model (and it has been done in the past), in a complicated and well established organized NWP developing centre, no national or international weather forecasting center has crossed into the area of wildfire forecasting. As a result, natural resource organizations with limited meteorological staff or developers focus on 2D spread models, and weather prediction organizations use a dispersion model approach for wildfire emission forecasts with biases due to a lack of wildfire heat, moisture and emissions coupling. Similar to weather and climate, wildfires know no national boundaries; thus, national and international cooperation at a large scale is essential for the development of an advanced operational coupled atmosphere–wildfire prediction system equipped with data assimilation and numerical weather forecast experts. Wildfire management is done mainly by provincial authorities in Canada, and they receive national help in times of emergency at their request. However, it unclear how they can foresee the requirement help without an accurate prediction method. No province in Canada has the potential to adopt an operational coupled weather–wildfire model simply because of the lack of an operational numerical weather prediction (NWP) center with data assimilation capacity and sufficient expertise. The temporary solutions in the past, such as short-term contracts, graduate research grants and IT developers, did not result in an evolving operational solution that can be updated with advancements in observation and fire science.
Regarding the areas of improvement within the current paper, we recognize that the verification of a model is not scientifically possible at this time due to limited data and accessibility. We also realize that satellite fire detection has limitations including time lapses, false alarms and undetected fires or misrepresentations of the actual burned area and fire size. Additionally, we found that the fuel classification can be improved by using a dynamic approach with high-resolution vegetation satellite images; this will be another research topic in future. The uncertain initialization and initial conditions lead to biased results, and the bias increases by the lead time in simulation. Advancements in initialization improve the overall outcomes.