Next Article in Journal
Merging of Horizontally and Vertically Separated Small-Scale Buoyant Flames
Next Article in Special Issue
Assessing Potential Safety Zone Suitability Using a New Online Mapping Tool
Previous Article in Journal
A Historical Perspective of Landscape and Human Population Dynamics in Guimarães (Northern Portugal): Possible Implications of Rural Fire Risk in a Changing Environment
Previous Article in Special Issue
A Surrogate Model for Rapidly Assessing the Size of a Wildfire over Time
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Considerations for Categorizing and Visualizing Numerical Information: A Case Study of Fire Occurrence Prediction Models in the Province of Ontario, Canada

1
Aviation Forest Fire and Emergency Services, Ministry of Northern Development, Mines, Natural Resources and Forestry, Sault Sainte Marie, ON P6A 6V5, Canada
2
Aviation Forest Fire and Emergency Services, Ministry of Northern Development, Mines, Natural Resources and Forestry, Dryden, ON P8N 2Z5, Canada
3
Department of Statistical and Actuarial Sciences, University of Western Ontario, London, ON N6A 5B7, Canada
4
Great Lakes Forestry Centre, Canadian Forest Service, Natural Resources Canada, Sault Sainte Marie, ON P6A 2E5, Canada
5
Graduate Department of Forestry, University of Toronto, Toronto, ON M5S 3B3, Canada
6
Aviation Forest Fire and Emergency Services, Ministry of Northern Development, Mines, Natural Resources and Forestry, Peterborough, ON K9J 3C7, Canada
*
Author to whom correspondence should be addressed.
Submission received: 18 June 2021 / Revised: 6 August 2021 / Accepted: 12 August 2021 / Published: 18 August 2021
(This article belongs to the Special Issue Wildfire Hazard and Risk Assessment)

Abstract

:
Wildland fire management decision-makers need to quickly understand large amounts of quantitative information under stressful conditions. Categorization and visualization “schemes” have long been used to help, but how they are done affects the speed and accuracy of interpretation. Using traditional fire management schemes can unduly restrict the design of new products. Our design process for Ontario’s fine-scale, spatially explicit, daily fire occurrence prediction (FOP) models led us to develop guidance for designing new schemes. We show selected historical fire management schemes and describe our method. It includes specifying goals and requirements, exploring design options and making trade-offs. The design options include gradient continuity, hue selection, range completeness and scale linearity. We apply our method to a case study on designing the scheme for Ontario’s FOP models. We arrived at a smooth, nonlinear scale that accommodates data spanning many orders of magnitude. The colouring draws attention according to levels of concern, reveals meaningful spatial patterns and accommodates some colour vision deficiencies. Our method seems simple now but reconciles complex considerations and is useful for mapping many other datasets. Our method improved the clarity and ease of interpretation of several information products used by fire management decision-makers.

1. Introduction

Situational awareness and decision-making for operational wildland fire management is supported by a large amount of complex, numerical information, often covering large areas and sometimes spanning multi-day forecasts. Comprehending and interpreting that quantity of information under time-limited and stressful conditions is challenging. Among other ways, this task is commonly made faster and easier by categorizing and visualizing the numerical information. There are many ways to do so, but how it is done can help or hinder interpretation, highlight or obscure valuable information and accurately portray or distort the data.
The categorization of numeric indicators of potential fire activity has a long history in Canadian and other fire management agencies. Established categories are deeply integrated into fire operations and culture. Although categorization is useful, conforming to traditional schemes, such as a four-category blue–green–yellow–red sequence from low to extreme, may be less than ideal for new information products. This issue arose during our implementation of fine-scale, spatially explicit fire occurrence prediction (FOP) models. Our design process for a new scheme is the basis for this paper.
FOP is one of the pillars of situational awareness, a requirement for daily and multi-day preparedness planning [1] and a key component for modelling risk [2]. Ontario’s fire management agency has a lightning-caused FOP model [3] and, more recently, a human-caused FOP model [1] that was developed in collaboration between the agency’s science specialists and decision-makers and external researchers.
Our objective is to provide guidance for the design of the categorization and visualization of complex information used in fire management. We begin with an overview of selected historical categorization and colouring schemes used in fire management. We then describe our method in steps including (1) specifying design goals and requirements, (2) exploring design options and seeing how they interact and (3) making trade-offs. We present a case study on designing the scheme for displaying the output of Ontario’s FOP models. Although our method arose from this work, the considerations and general principles employed can be useful in many other situations. The Supplementary Material includes other applications of our method and extensions to some other considerations for the display of data. Note that we use the terms “category” and “categorization” synonymously with “class” and “classification”, the latter set being conventions in fire management.

Overview of Selected Historical Categorization and Colouring Schemes

Knowing the historical origin of the schemes aids us to understand their limitations and engender improvements. Prime categorization examples are those for the outputs of two Canadian Forest Fire Danger Rating System (CFFDRS) [4] subsystems: the Fire Weather Index (FWI) System [5] and the Fire Behaviour Prediction (FBP) System [6].
The FWI System accounts for effects of past and present weather on fuel ignitability and fire behaviour and has six main numeric outputs. Three track moisture in different depths or sizes of fuel: Fine Fuel Moisture Code (FFMC), Duff Moisture Code (DMC) and Drought Code (DC). The other three indicate potential fire behaviour: Buildup Index (BUI), Initial Spread Index (ISI) and Fire Weather Index (FWI). FWI indicates the potential fireline intensity in a standard pine (genus Pinus) stand [5]. FWI is also used in Ontario as a general indicator of fire hazard or danger and is mapped into ordinal adjective classes: Low–Moderate–High–Extreme. The classes accompanied the introduction of the FWI System [5] and were implemented early on by Ontario [7]. Fire danger schemes are intended to “sound an alarm” about potential extreme behaviour and difficulty of control [8]. The class boundaries were set so that Extreme covered the worst 2% of historical days, and the remaining lower boundaries were set by a geometric progression [5]. For a more detailed history of the FWI classification methods in Ontario, see [9]. Other fire agencies in Canada calculated different boundaries according to their data [10,11]. The FWI System and component classifications have also been applied outside North America. For examples, see [12,13]. The current danger classes and colours used in Ontario for the FWI System’s main outputs are outlined in Table 1. Examples of these are shown for roadside signs for public alerts about daily fire hazard (Figure 1a) and also for operational maps (Figure 1b). In mapping products such as Figure 1b, the FWI System values are calculated for weather station locations and are interpolated [14] across Ontario’s fire management area, which is almost 50% larger than France.
The other major CFFDRS subsystem that has outputs commonly communicated using classes is the FBP System, which provides quantitative estimates of fire behaviour outputs [6]. A primary output is fire intensity, the rate of energy or heat release per unit time per unit length of a spreading fire front [16], which ranges from 1 to ~100,000 kW/m in Canadian conditions. Fire intensity is also commonly categorized into fire intensity classes (ICs). Rather than adjective classes (i.e., low–extreme) they were given five numeric labels (IC 1–5 [17]), and later six (IC 1–6 [18]). Higher IC numbers correspond with the higher intensity values, but the boundaries are not evenly spaced (Table 2). The most commonly used IC boundaries in Canada delineate distinct differences in fire type characteristics (for example, surface, torching or crowning) in mature jack pine (Pinus banksiana Lamb.) stands and the corresponding general effectiveness of different types of fire suppression activities (for example, hand tools, pumps and hose, airtankers). These ICs are described in the Field Guide to the FBP System [19] and the ICs are also used to map fire intensity by many fire management agencies, for various purposes (Table 2).
Table 2 shows some of the different intensity values for higher-end class boundaries (i.e., additional thresholds beyond IC 6). The choice of colour when mapping can be operationally significant because colours covey information rapidly and have strong associations with levels of alarm—for example, red for danger [20] and green for calm [21]. Such psychological factors are not always considered in the visualization, however, which could lead to misinterpretation. For example, IC 4 in Table 2, which is associated with the upper limit of direct fire suppression effectiveness [17,22], is variously coloured a calming light green, a cautionary yellow or a warning orange. In Table 2 there are also cases where the same colour refers to different IC classes—for example, calming light green is used for IC 2 in Ontario, IC 3 in Alberta and IC 4 nationally.
An additional caveat in the classification is that simplifying numerical information by aggregation into few classes has a cost. For FWI and fire intensity, the wide range within some classes is operationally significant for some decisions—for example, FWIs of 11–22 become “High”. Furthermore, for interpolated maps, neighbouring points that are displayed as different classes will not have operationally meaningful differences. To compensate, maps often include the raw point values that were used for interpolation (Figure 1b). The numeric information cannot, however, be read and interpreted as quickly, thus reducing the benefit of categorization. The categorizing of data for spatial application such as FWI and FBP is common, and there are many recognized considerations (for example, [26]) and built-in solutions in geographic information systems. However, using a built-in classification option without a deep understanding may be unsuitable because potential distortions can lead to radically different interpretations, as others have noted [27]. Consequently, there is a strong need to use schemes that convey information accurately.

2. Methods

We propose five steps to categorize and visualize model outputs for use in fire situational awareness and decision-making:
  • Understanding and scoping the data
  • Understanding the decision-making uses of the information
  • Specifying the design goals and requirements
  • Designing the categorization and visualization scheme
  • Evaluating and revising the scheme
We first describe these steps in general, below, and then with further detail on their application, in our case study. Although the method is described in a linear sequence, the work is partly concurrent and highly iterative, especially within Step 4.

2.1. Step 1: Understanding and Scoping the Data

Our method applies to data with a continuous numerical scale of measure (real numbers). With minor modifications, it can also apply to data with a discrete numerical scale of measure (integers) and to ordinal categorical data (for example, Very Low, Low, Low–Moderate, …). The modification is that colouring with continuous gradients (described below) does not apply unless there are a great many discrete values or ordinal categories.
The technical details of the raw model output data may be straightforward, but unfamiliar units, scaling, storage or other conditions can lead to misinterpretation. The following need to be understood by the designers:
  • Units, including any scaling and transformations
  • Data storage type (for example, 64-bit floating point, signed long integer, string) and associated considerations (for example, storing a scaled real value as an integer, which truncates the precision)
  • The data’s range or anticipated range if using a static scale for all future maps
  • The frequency distribution of historical data
    It may or may not be useful to have more categories where there was a lot of data, and it may or may not be pointless to have multiple categories where there was little or no data (see examples in Section 2.2).
  • The data’s precision of measurement and storage and the data’s accuracy of measurement or estimation
    The stored precision may not correspond with the accuracy. Measured data such as weather observations have low precision (for example, to 0.1 °C), but calculated data such as FWI System values should be calculated and may be stored with full machine precision (~16 decimal digits), which is well beyond the accuracy of typical fire management data.
  • If the data are generated by a model, then the model’s meaning, structure, assumptions, limitations, precision and accuracy
Understanding data precision and accuracy is necessary for presenting information accurately and having decision-makers understand it easily and correctly. Regarding data storage and all subsequent calculations using data, full machine precision should be maintained to avoid accumulating rounding errors. Regarding the numbers displayed for decision-makers, the displayed precision should not exceed the data’s accuracy, because that could be misleading. The numbers displayed for decision-makers should ideally have the lowest precision that is operationally significant to minimize unnecessary mental processing. Further discussion and examples are given in the Supplementary Material.

2.2. Step 2: Understanding the Decision-Making Uses of the Information

The purpose of the information is to support decision-making, so it is necessary to know who is using the information and how it is used. Working directly with fire management staff to understand their needs is necessary for ensuring that new model outputs are effectively integrated into the decision-making process [28,29].
There are key questions to consider. What decisions are being supported? Are certain parts of the range more important, needing higher attention? Is a higher resolution (smaller class size) needed in some parts of the range rather than others? For example, consider the categorizing and mapping of the accumulated 24-h rainfall from a precipitation radar [30]. For differentiating the degree and duration of the reduced fire behaviour potential, a high resolution is useful at the low end but not at the high end. Conversely, for differentiating the degree and duration of flood potential, a high resolution is useful at the high end but not at the low end. Moreover, a higher top category is appropriate.

2.3. Step 3: Specifying the Design Goals and Requirements

As stated in the introduction, the high-level goals of categorization and its visualization are simply to show the information completely and accurately and to have the information be understood quickly and easily. These goals are elaborated into criteria as follows.
  • Regarding the complete and accurate display of information:
    Are the magnitudes shown with the original or reduced precision?
    Are the magnitudes undistorted or distorted by categorization, scale nonlinearity or truncation?
    Are the relative magnitudes evident by the colouring without or with referral to the legend?
  • Regarding the quick and easy understanding of information:
    Can the colouring be easily matched to the legend’s magnitude numbers?
    Does the colouring draw attention and convey a suitable psychological meaning for the degree of alarm?
    What is the overall ease of understanding?
Possible design requirements include the accommodation of colour vision deficiencies and other technical considerations such as the adequate appearance on low-quality displays, colour printers or standard photocopiers.

2.4. Step 4: Designing the Categorization and Visualization Scheme

There are four design options for categorizing and colouring the values in the scale (Figure 2):
  • Gradient continuity: whether to use the original values unaltered or categorized
  • Hue selection: the number and choice of colours and design of gradients
  • Range completeness: whether to show the full range or truncate the top or bottom of the range
  • Scale linearity: whether to have a linear or nonlinear scale or progression of category boundaries, and whether to have colour gradients that are linearly or nonlinearly proportional to the data magnitudes

2.4.1. Gradient Continuity

The alternatives for gradient continuity are either to use the original values or categorize them (Figure 2, parts 1a and 1b). For a continuous gradient, the percentage of colour saturation is proportional to the raw datum magnitude. Continuous gradients are therefore precise and accurate but harder to interpret using the legend compared to categorized gradients.

2.4.2. Hue Selection

Colours have strong psychological associations that affect the inferred meaning of information and its speed and ease of interpretation. Blue and green are associated with relaxation, calm and hope [21], and red is associated with danger [20] (Figure 2, part 2f). Historically, variations of a blue–green–yellow–orange–red sequence (which alludes to water, growing vegetation, dried vegetation and flame) have been used to represent escalating fire danger (see examples in Table 2). Accommodating colour vision deficiencies reduces the colour combination choices, particularly most of those in the traditional fire danger sequence [31]. Tools are available to assist with the testing of colour palettes for accessibility [31,32]. The design task is to choose colours appropriate for the implications of the data magnitudes, particularly the degree of attention or alarm.
Regarding the gradient design, there are a few distinct alternatives [33] (Figure 2, parts 2a–2e). A single sequential gradient is for data ranging over a meaning of zero or neutral to bad or good. A divergent sequential gradient is for data ranging over a meaning of good through neutral to bad. In Figure 2, parts 2b and 2c grade from a calming blue through white to an alarming red. The left and right ends each have a single hue. Part 2b is neutral in the middle, whereas part 2c has compressed and expanded ends. Part 2d, multiple sequential, is analogous to part 2a except that part 2d has multiple hues, which are in a rainbow spectrum in the example. Compared to a single hue, multiple hues provide more contrast over the range, making it easier to match the legend and signal levels of attention or alarm. Part 2e, multiple divergent, is analogous to part 2b except that part 2e has multiple hues for each side. The R software [34] package "inlmisc" [35,36] is useful for constructing continuous or categorized gradients.
A key concern is how these many design alternatives support or oppose the design goals. Only the single sequential and single divergent gradients (Figure 2, parts 2a and 2b) have the accuracy of continuous gradient continuity (Figure 2, part 1a), but they have a difficult interpretability. The remaining gradient alternatives require matching the legend to identify the magnitudes, but this can become quick and easy to interpret with familiarity and an effective use of colour psychology.

2.4.3. Range Completeness and Scale Linearity

These are described together because an incomplete range is an extreme form of nonlinearity. The alternatives for range completeness (Figure 2, part 3) are whether to show the full range of the data or truncate the top or bottom and group the truncated data. All resolution is lost beyond the truncation points. The alternatives for scale linearity (Figure 2, part 4) are for two components, independently:
  • Numeric scale: whether to have a linear (part 4a) or nonlinear numeric scale (parts 4b and 4c) or progression of category boundaries (if applicable)
    • For categorized gradient continuity, a nonlinear scale is used to vary the resolution over the range
  • Colour gradient: whether to have colour gradients that are linearly or nonlinearly proportional to the data magnitudes
    • For categorized gradient continuity, a nonlinear colour gradient is used to communicate the varying meaning or importance of the information over the range
For a continuous gradient, the same result can be achieved from nonlinearity in either of the above two components.
Nonlinear-systematic (part 4b) methods use a smooth function such as log or power to transform the output, while nonlinear-irregular methods use a non-smooth progression such as Jenks [37]. Nonlinear colouring requires the referral to the legend to understand the magnitudes. This is a trade-off between the goals of drawing attention to where it is needed and improving the speed and ease of understanding.

2.4.4. The Design Process

There are copious settings and combinations of alternatives for the four design options. Getting to a result is an iterative process, with analysts and subject matter experts trying alternatives and making trade-offs, hopefully avoiding the anchoring to tradition or early trials. We cannot recommend a path through the four design options other than saying it is iterative and concurrent. We do, however, recommend a starting point or baseline, which is the extreme of displaying all the data completely and accurately and ignoring the goals of quick and easy understanding. The baseline has continuous gradient continuity, an achromatic colour gradient of white through greys to black and a linear scale with no truncation (Table 3).
If categories are used, determining the number and their boundaries are fundamental design decisions [38]. The number of categories is a trade-off of accuracy (requiring more) and speed and ease of understanding (requiring fewer). Ideally, individual categories have no operationally significant physical differences, while adjacent categories do. In practice, all considerations require compromise. An example of determining categories for FWI System outputs based on physical differences is given by [9].

2.5. Step 5: Evaluating and Revising the Scheme

Once the design process is done and implemented, an essential further step is the ongoing work with decision-makers to evaluate the outputs and revise the design as necessary.

3. Case Study: Designing the Scheme for Ontario’s FOP Models

We now describe the application of the above method to categorizing and visualizing FOP data.

3.1. Case Study—Step 1: Understanding and Scoping the Data

The data are outputs from process and statistical FOP models for Ontario, so we begin with their description. The lightning- and the human-caused FOP models have been used operationally since the mid-2000s and 2015, respectively. The daily lightning-caused fire occurrence is modelled as two separate processes [3]: the probability of a lightning strike will lead to a holdover ignition and the probability that an existing ignition “arrives” (is reported). The ignition model is mainly driven by the forest floor organic layer’s moisture content, which determines the sustainability of smouldering and the survivability of the ignition. Additional factors are other moisture indicators, ecoregional modifiers and lightning strike polarity. The “arrival” model, which is conditional on a holdover ignition being present, is influenced by the surface litter moisture, organic layer moisture, wind speed and ecoregional differences. Ontario’s human-caused fire prediction system uses a set of logistic generalized additive models to model inherently nonlinear relationships with key drivers of human-caused fire occurrence, including seasonal and spatial patterns, fuel moisture and the characteristics of human land use [1]. The models are stratified regionally and by cause categories to account for different seasonal patterns in fire occurrence.
Both the lightning- and human-caused FOP models produce outputs for each of the 2574 cells in a grid that spans the province’s approximately 91.9 million ha wildland fire management area (Figure 3). Most of the cells are about 20 km × 20 km or 40,000 ha, with some fractional cells at the boundaries. The units of the FOP output data are interpreted as the expected number of fires per cell, regardless of the cell size. The model calculations and outputs have a 64-bit floating point precision, but the data are transferred to the mapping software via a text file holding up to 12 significant digits.
Regarding the range and frequency distribution, we analysed historical FOP data for each cell for each day from 15 May to 31 August; 2016–2018 for human-caused and 1992–2006 for lightning-caused fires. For this analysis, the start and end dates in each fire season were chosen to avoid the variability in spring and fall snow-free conditions, when the models are not making predictions for the entire province. The lower limit of the range is zero; there is no theoretical upper limit. Table 4 presents summary statistics for the data, with zeros excluded to characterize the important data more clearly. Most of the distribution statistics of the human- and lightning-caused data differ by about an order of magnitude.
Figure 4 shows the empirical probability distributions of the non-zero data. Those data were mostly clustered close to zero in both models, so we log-transformed the data for illustration (untransformed data are required for operational use). The magnitudes of data between the lower tails of the two distributions differ by orders of magnitude. This presents a challenge for categorizing and colouring the data on a common scale for mapping.

3.2. Case Study—Step 2: Understanding the Decision-Making Uses of the Information

We used a variety of methods to understand the FOP information needed and how it is used for daily decision-making: reviewing documentation, observing operational decision-making and (for some) working part-time in operational functions where FOP information is used. Most importantly, we held a series of engagements with fire management agency personnel. For example, we hosted a workshop in 2017 attended by agency personnel including regional and provincial Fire Intelligence Officers, external researchers and students. The workshop’s topics included the purposes and methods of subjective FOPs by experts.
To understand the agency’s use of FOP information, it is necessary to outline the agency’s hierarchical structure and the responsibilities of each level. Ontario’s fire management area (Figure 3) has two main parts, the Northeast and Northwest Regions, each of which is divided into six or seven Fire Response Sectors. There is also a Provincial level. Each level has a sole or shared responsibility for various decisions; most are made with consultation or coordination between adjacent levels. The Regions are primarily responsible for strategic fire response decisions and management, and the Sectors are primarily responsible for tactical fire response decisions and operations. The Province is primarily responsible for adjusting the near-term (1- to 21-day) capacity according to the demand via temporary commercial hiring and inter-provincial and international resource sharing.
Several decisions are directly dependent on the potential number of fires anticipated in various parts of the province. These decisions are made at the stated levels:
  • Prevention: for example, escalated and targeted messaging, temporary fire bans; made by the Province, Regions and Sectors
  • Preparedness: types, numbers, locations and readiness alert levels of firefighting resources (crews, helicopters, airtankers and engines); made by the Regions and Sectors
  • Detection: numbers, routes and times of aerial detection patrols described in [39]; made by the Regions, jointly
  • Dispatch: different resources may be sent now if more threatening fires are anticipated later; made by the Regions and Sectors
  • Inter-provincial and international resource sharing; made by the Province and Regions
The various FOP-dependent decisions have diverse needs in terms of the spatial extent and resolution of FOP information. For example, detection route designers can use relatively fine resolution information on the order of kilometres, while the Province needs only aspatial, numeric FOPs by region for resource-sharing decision-making.

3.3. Case Study—Step 3: Specifying the Design Goals and Requirements

The primary and conflicting goals are of course to display complete and accurate information and have the information quickly and easily understood. In twice-daily briefings, decision-makers have limited time (minutes) to view, interpret and absorb each of many information items regarding, for example, weather values, FWI and FBP System outputs, FOP, active fires, logistics and personnel. Completeness and accuracy are important because the information supports the many decisions described above, which are made under uncertainty and have potentially significant consequences.
Regarding specific information requirements:
  • There is a need for both maps and numeric subtotals and totals of fire occurrence by cause and location (i.e., Sectors, Regions, Province)
  • All the maps need to use the same categories and colours for FOP magnitudes
  • The FOP models’ output is the expected or average occurrence, but the actual occurrence varies around the average, so an indication of the variability is needed.
Decision-makers expressed strong preferences for the number of categories, ranging from three to many categories, and they desired a familiar colour sequence (blue–green–yellow–red). There was also a strong preference for integers for all numbers related to FOP. Finally, we wished to accommodate colour vision deficiencies.

3.4. Case Study—Step 4: Designing the Categorization and Visualization Scheme

Design has been described as a messy process with a tidy outcome. We do not detail our circuitous journey but show and describe some key alternatives, stages and considerations. Figure 5a illustrates the baseline alternative (Table 3) applied to the human-caused FOP and actual fire arrivals for a selected day. The range of the colour gradient is 0–3 fires/cell, the upper limit of which is between the maximum human- and lightning-caused FOP (Table 4). The map area looks mostly white, with three small, pale grey patches; there is little useful information, especially considering the six actual fires that day, which is a low-to-moderately busy day for this cause. Adding more hues alone would make no meaningful difference because the data are clustered very near zero. Categorizing at this stage would make it worse. Truncating the upper limit to a low value somewhat close to zero would add resolution and colour to the human-caused FOP here, but such truncation would lose all resolution of the rarer but critically important high lightning-caused FOP.
Our original solution was to use separate scales for the two fire occurrence causes (Figure 6). The lightning-caused FOP scale was linear. For the human-caused FOP scale, we led subject matter experts through a scenario for an area of Ontario that has a relatively high occurrence. When the FFMC (a strong indicator of sustainable ignition) was 90, that area was considered to have an elevated concern suitable for a High classification. We used the corresponding FOP magnitude (0.4 fires/cell) as the upper limit for the then 10-category scale. We determined Moderate similarly and interpolated with equal linear steps. In Figure 6, this scale is shown by the blue line, except that categories 1 to 10 in the original have been mapped to a 0 to 20 scale for comparability. Scaling the original total-fires map required a creative logic. The maps by individual cause were acceptable to decision-makers, but the inconsistency between the causes was ultimately unacceptable (and motivated the present work).
Several alternatives were considered for a unified scale that accommodated the conflicting needs of fine resolution at the low end and a high upper limit. Discussions with decision-makers indicated that concern increases relatively quickly as the likelihood of fire rises from zero. Providing a high resolution at the low end while retaining a high upper limit would require a great many colours (for example, like those of precipitation radar maps). That would be unfamiliar and confusing and would not correspond with psychological colour associations nor accommodate colour vision deficiencies. Piecewise linear and irregular scales were explored, but their abrupt changes made them difficult to interpret. We wanted a smooth, systematic progression of category boundaries and tested logarithmic and power functions. Those functions can be made to match fairly closely, but the power function had a more suitable shape at the low end. A power function takes the general form f ( x ) = a x b , with the shape controlled by parameters a and b . Our desired behaviour for the scale to increase quickly for low FOP but then increase progressively more slowly is provided when a > 0 and 0 < b < 1 , since this family of power functions is monotonically increasing and concave down. The FOP scale is x , and the category scale is f ( x ) . We developed a parametric scaling tool with three inputs to generate and plot boundaries: shape parameter, 1 b ; the upper limit of the FOP scale, F O P M a x ; and the number of categories, N u m C a t . The boundary for the top of category C a t is
B o u n d a r y C a t = F O P M a x · ( C a t N u m C a t ) 1 b   ,   C a t = 1 ,   2 ,   3 ,   ,   N u m C a t .
Any FOP > F O P M a x stays in the highest category. A category for true zero can be added if required. For convenience, we parameterized 1 b , for which we tested values in the range of 1.5 to 4.5; for example, 2 yields a square root shape. a works out to be N u m C a t F O P M a x b . The red curve in Figure 6 illustrates the boundaries for 1 b = 3 , F O P M a x = 3 and N u m C a t = 20 ; for example, the boundary for C a t = 11 (the top of category 11) is ~0.5. The tool lists the boundaries and graphs them as in Figure 6.
While working directly with subject matter experts, the tool facilitated the joint testing of alternatives for the number of categories, truncation and nonlinearity design options. These alternatives plus colouring needed to be adjusted simultaneously when trading off the goals because their effects interact. The FOP outputs for a set of representative days were mapped using R [35] for candidate sets of boundaries and colouring. Table 5 lists ways in which the alternatives for the number of categories, the amount of truncation and nonlinearity and number of hues generally interact in affecting several attributes related to completeness and accuracy of information or speed and ease of understanding. There are exceptions for some combinations and edge conditions. Every alternative for the design options improves some attributes and worsens others. The trade-off behaviour is more straightforward to work with than it may seem because the tool shows most of the trade-offs immediately. The difficulty lies in subjectively assessing the results and compromising on the attributes and goals.

3.5. Case Study—Step 4: Results of the Design Process

We state our current category and colouring design and give the rationale for the trade-offs made. The pressure to have a small number of categories and colours (~4) could not accommodate the need for fine resolution at the low end; we used 20 categories (plus a true zero if needed). Figure 5b–d show the same FOP data as Figure 5a but with 4, 10 and 20 categories, respectively. Only the largest number of categories reveals the meaningful network pattern of lines and nodes that corresponds roughly with roads and settlements.
In addition, a highly nonlinear scale was needed to show the patterns in the data. We used a nonlinear-systematic scale with boundaries obtained using Equation (1) with parameters 1 b = 3 , F O P M a x = 3 and N u m C a t = 20 . The boundaries are given in Table 6, stated in units of fires/cell and cells/fire. The final shape of the nonlinear scaling corresponds to parameter settings of a = 20 · 3 1 3 and b = 1 3 .
f ( x ) = { 20 · 3 1 3   x 1 3 , 0 x < 3 20 , x 3 ,
which is illustrated by the red curve in Figure 6.
Regarding truncation, we considered the rarity and operational importance of extreme FOP magnitudes. The 20th category ends at 3 fires/cell according to Equation (1), but that category is used for all higher FOP magnitudes. The maximum FOP in Table 4 is ≈3.9 fires/cell.
Regarding colouring, we used a nonlinear, divergent scale with one hue for the low end and multiple hues for the high end (Table 6). The hues transition from light blue to yellow through orange to red, which mostly follows the traditional blue-to-red progression. Avoiding green in that sequence accommodates some types of colour vision deficiency [31]. Even though there are 20 categories, having four main colours is easy to interpret and consistent with other CFFDRS outputs (for example, Table 1). The gradient within each main colour is difficult to match with the legend, but nonetheless reveals meaningful spatial patterns in the maps. Figure 7a–c show the final categorization and colouring scheme in example daily FOP maps of human- and lightning-caused fires and total fires, respectively. We intentionally show maps as formatted for operational use. They are intended for display on large monitors, but these reduced versions still show the colouring and categorization results adequately. Larger versions are provided in the Supplementary Material.
Note that a sequential scale is logical for FOP because any non-zero FOP is “bad” in this context. But the calming colours are assigned to very low magnitudes, and this provides a slightly greater distinction between the remaining colours. A much greater distinction could easily be achieved by adding more colours—for example, magenta–purple–black for the highest categories, which would also draw more attention to the critical extremes. This, however, would not accommodate colour vision deficiencies.
Broad adjective categories were added to emulate the familiar four- or five-category pattern and simplify interpretation. Those boundaries fit the general association between Ontario’s fire arrival density and fire situation severity. The three-significant-digit category boundaries were replaced in the legend by integers for selected category midpoints or boundaries (Table 6). Note that the units change from fires/cell to cells/fire to show integer magnitudes, which are far more meaningful than fractions.

3.6. Case Study—Step 5: Evaluating and Revising the Scheme

Several significant revisions were done in arriving at the current design in Figure 7. The first lightning-caused FOP maps had a linear scale with four equal-interval categories using traditional blue–green–yellow–red. When human-caused FOP was added, the scale was changed to 10 categories using a smooth green–yellow–orange–red transitional gradient and later to a special blue–yellow–orange–red to accommodate some colour vision deficiencies. As stated above, the inconsistent subjective scales were replaced completely in early 2020 per our method. The maps originally showed the expected number of fires/cell but were changed to show the density of the expected number of fires/unit area because of fractional cells. Note that the map legends intentionally omit this complication; density seems to be the automatic, intuitive interpretation. Additional revisions are planned, and further evaluation is always ongoing.

4. Discussion

While implementing the FOP models for operational evaluation, we arrived at and applied this method, which may seem well structured and straightforward now, but it was far from that during the design work. The method emerged as a by-product, having evolved during iterative deliberations. As illustrated in the introduction and identified in [9], the categorization of the FWI System outputs could benefit from this approach. We have since applied the method to categorize and colour outputs from other models for operational display [39,40,41], examples of which are in the Supplementary Material.
The many considerations discussed in this paper need to be addressed to ensure that the models are interpreted and used appropriately. We emphasize that any schemes including those built into software applications need this careful consideration. Incorporating the outputs from scientific models into operational decision-making is not straightforward. The key for a successful application is that researchers and practitioners work together closely throughout the process, from problem identification through to implementation and evaluation [1]. Through this collaborative approach, outcomes tend to have a higher acceptance and usefulness.
An additional factor not addressed by our method is that categories need to be meaningful for more than just map design, because there is a tendency to extend the use categories to guidelines and standard operating procedures and vice versa [9]. The classifications and their boundaries can also be used as unwritten mental shortcuts or heuristics in the place of a more deliberative consideration of complex information. Well-designed, science-based classifications can be consistent with less time-constrained situational analyses, while poor classifications may lead to suboptimal decision-making.
Design considerations for presenting quantitative data for decision support go beyond the categorization and colouring of numerical scales. Also important are options for spatial resolution and the use of simulated three-dimensional displays, examples of which are given in the Supplementary Material [42].
We presented a method of designing classification and colouring schemes that consider many complex factors, interactions and trade-offs. These satisfied the ultimate goals of showing decision-makers complete and accurate quantitative information that was understood quickly and easily.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/fire4030050/s1: other applications of the model (Figures S1 and S2), classification of numbers, additional considerations for the spatial display of data (Figures S3–S5) and larger versions of Figure 7a–c (Figures S6–S8).

Author Contributions

Conceptualization: D.B., J.E., C.C.H., C.B.M., D.G.W. and M.W. (Mike Wotton); methodology: D.B., C.B.M. and D.G.W.; software: J.E., A.S. and D.G.W.; formal analysis: D.B., C.C.H., C.B.M. and M.W. (Melanie Wheatley), D.G.W. and M.W. (Mike Wotton); investigation: D.B., J.E., C.C.H., C.B.M., A.S. and M.W. (Melanie Wheatley), D.G.W. and M.W. (Mike Wotton); writing—original draft preparation: D.B., C.C.H., C.B.M., M.W. (Melanie Wheatley) and D.G.W.; writing—review and editing: D.B., J.E., C.C.H., C.B.M., A.S., M.W. (Mike Wotton), D.G.W. and M.W. (Melanie Wheatley); visualization: D.B., J.E., C.C.H., C.B.M., M.W. (Melanie Wheatley), D.G.W., M.W. (Mike Wotton) and A.S.; project administration: D.B., C.B.M. and D.G.W. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge the support of the Natural Sciences and Engineering Research Council of Canada (NSERC) [RGPIN-2015–04221].

Acknowledgments

Thank you to the Ministry of Northern Development, Mines, Natural Resources and Forestry, Aviation Forest Fire and Emergency Services staff for important discussions and management for support. Thank you to Barry Graham, Dan Johnston, Dan Leonard and David Martell for insightful advice, Darren McLarty for the review, Jim Caputo, Natasha Jurko and Benito Russo for technical support and the Fire Intelligence Officer cadre for their regular engagement and feedback. A special thank you goes to Darryl Pajunen for supporting the model’s development. We thank the anonymous reviewers for their careful reading and helpful comments that improved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Woolford, D.G.; Martell, D.L.; McFayden, C.B.; Evens, J.; Stacey, A.; Wotton, B.M.; Boychuk, D. The Development and Implementation of a Human-Caused Wildland Fire Occurrence Prediction System for the Province of Ontario, Canada. Can. J. For. Res. 2021, 51, 303–325. [Google Scholar] [CrossRef]
  2. Johnston, L.M.; Wang, X.; Erni, S.; Taylor, S.W.; McFayden, C.B.; Oliver, J.A.; Stockdale, C.; Christianson, A.; Boulanger, Y.; Gauthier, S.; et al. Wildland Fire Risk Research in Canada. Environ. Rev. 2020, 28, 164–186. [Google Scholar] [CrossRef]
  3. Wotton, B.M.; Martell, D.L. A Lightning Fire Occurrence Model for Ontario. Can. J. For. Res. 2005, 35, 1389–1401. [Google Scholar] [CrossRef]
  4. Stocks, B.J.; Lynham, T.J.; Lawson, B.D.; Alexander, M.E.; Van Wagner, C.E.; McAlpine, R.S.; Dubé, D.E. The Canadian Forest Fire Danger Rating System: An Overview. For. Chron. 1989, 65, 450–457. [Google Scholar] [CrossRef]
  5. Van Wagner, C.E. Development and Structure of the Canadian Forest Fire Weather Index System; Forestry Technical Report; Canada Communication Group Publication: Ottawa, ON, Canada, 1987; ISBN 9780662151982. [Google Scholar]
  6. Forestry Canada Fire Danger Group. Development and Structure of the Canadian Forest Fire Behavior Prediction System; Information Report ST-X-3; Forestry Canada, Science and Sustainable Development Directorate: Ottawa, ON, Canada, 1992; ISBN 9780662198123.
  7. Stocks, B.J. Wildfires and the Fire Weather Index System in Ontario; Information Report O-X-213; Canadian Forestry Service, Great Lakes Forestry Centre: Sault Ste. Marie, ON, Canada, 1974.
  8. Taylor, S.W.; Alexander, M.E. Science, Technology, and Human Factors in Fire Danger Rating: The Canadian Experience. Int. J. Wildland Fire 2006, 15, 121. [Google Scholar] [CrossRef]
  9. Hanes, C.; Wotton, B.M.; McFayden, C.; Jurko, N. An Approach for Defining Physically Based Fire Weather Index System Classes for Ontario. Natural Resources Canada; Information Report GLC-X-30; Canadian Forest Service: Sault Ste. Marie, ON, Canada, 2021; in press.
  10. KiiI, A.D.; Lieskovsky, R.J.; Grigel, J.E. Fire Hazard Classification for Prince Albert National Park, Saskatchewan; Environ; Information Report NOR·X·S8; Environment Canada, Forestry Service, Northern Forest Research Centre: Edmonton, AB, Canada, 1973. [Google Scholar]
  11. Kiil, A.D.; Miyagawa, R.S.; Quintilio, D. Calibration and Performance of the CANADIAN Fire Weather Index in Alberta; Information Report NOR-X-173; Environment Canada, Canadian Forestry Service, Northern Forest Research Centre: Edmonton, AB, Canada, 1977. [Google Scholar]
  12. Alexander, M.E. Proposed Revision of Fire Danger Class Criteria for Forest and Rural Areas in New Zealand, 2nd ed.; National Rural Fire Authority: Wellington, New Zealand; Scion Rural Fire Research Group: Christchurch, New Zealand, 2008; 63p. [Google Scholar]
  13. De Groot, W.J.; Field, R.D.; Brady, M.A.; Roswintiarti, O.; Mohamad, M. Development of the Indonesian and Malaysian Fire Danger Rating Systems. Mitig. Adapt. Strat. Glob. Chang. 2007, 12, 165. [Google Scholar] [CrossRef] [Green Version]
  14. Flannigan, M.D.; Wotton, B.M. A Study of Interpolation Methods for Forest Fire Danger Rating in Canada. Can. J. For. Res. 1989, 19, 1059–1066. [Google Scholar] [CrossRef]
  15. Ontario Forest Fire Information Map. Available online: https://www.lioapplications.lrc.gov.on.ca/ForestFireInformationMap/index.html?viewer=FFIM.FFIM (accessed on 17 June 2021).
  16. Byram, G.M. Combustion of Forest Fuels. In Forest Fire: Control and Use; Davis, K.P., Ed.; McGraw-Hill: New York, NY, USA, 1959; pp. 61–89. [Google Scholar]
  17. Alexander, M.E.; De Groot, W.J. Fire Behavior in Jack Pine Stands as Related to the Canadian Forest Fire Weather Index (FWI) System; Poster with Text; Canadian Forest Service, Northern Forestry Centre: Edmonton, AB, Canada, 1988.
  18. Taylor, S.W.; Pike, R.G.; Alexander, M.E. Field Guide to the Canadian Forest Fire Behavior Prediction (FBP) System; Natural Resources Canada, Canadian Forest Service, Northern Forestry Centre and British Columbia Ministry of Forests, Research Branch: Victoria, BC, Canada, 1996.
  19. Taylor, S.W.; Alexander, M.E. A Field Guide to the Canadian Forest Fire Behavior Prediction (FBP) System, 3rd ed.; Natural Resources Canada, Canadian Forest Service, Northern Forestry Centre: Edmonton, AB, Canada, 2018.
  20. Pravossoudovitch, K.; Cury, F.; Young, S.G.; Elliot, A.J. Is Red the Colour of Danger? Testing an Implicit Red–Danger Association. Ergonomics 2014, 57, 503–510. [Google Scholar] [CrossRef]
  21. Kaya, N.; Epps, H.H. Relationship between Color and Emotion: A Study of College Students. Coll. Stud. J. 2004, 38, 396. [Google Scholar]
  22. Hirsch, K.; Martell, D. A Review of Initial Attack Fire Crew Productivity and Effectiveness. Int. J. Wildland Fire 1996, 6, 199. [Google Scholar] [CrossRef]
  23. Alberta Wildfire, Danger Forecast, Headfire Intensity. Available online: https://wildfire.alberta.ca/files/ahfi.gif (accessed on 17 June 2021).
  24. Canadian Wildland Fire Information System Fire Behaviour Maps. Available online: https://cwfis.cfs.nrcan.gc.ca/maps/fb (accessed on 17 June 2021).
  25. BC Wildfire PSTA Head Fire Intensity. Available online: https://catalogue.data.gov.bc.ca/dataset/bc-wildfire-psta-head-fire-intensity (accessed on 9 September 2020).
  26. Evans, I.S. The Selection of Class Intervals. Trans. Inst. Br. Geogr. 1977, 2, 98–124. [Google Scholar] [CrossRef]
  27. Beverly, J.L.; McLoughlin, N. Burn Probability Simulation and Subsequent Wildland Fire Activity in Alberta, Canada—Implications for Risk Assessment and Strategic Planning. For. Ecol. Manag. 2019, 451, 117490. [Google Scholar] [CrossRef]
  28. Martell, D. The Development and Implementation of Forest and Wildland fire Management Decision Support Systems: Reflections on Past Practices and Emerging Needs and Challenges. Math. Comput. For. Nat. Resour. Sci. 2011, 3, 18. [Google Scholar]
  29. Noble, P.; Paveglio, T.B. Exploring Adoption of the Wildland Fire Decision Support System: End User Perspectives. J. For. 2020, 118, 154–171. [Google Scholar] [CrossRef]
  30. Hanes, C.C.; Jain, P.; Flannigan, M.D.; Fortin, V.; Roy, G. Evaluation of the Canadian Precipitation Analysis (CaPA) to Improve Forest Fire Danger Rating. Int. J. Wildland Fire 2017, 26, 509. [Google Scholar] [CrossRef]
  31. Jenny, B.; Kelso, N.V. Color Design for the Color Vision Impaired. Cartogr. Perspect. 2007, 61–67. [Google Scholar] [CrossRef]
  32. Harrower, M.; Brewer, C.A. ColorBrewer.org: An Online Tool for Selecting Colour Schemes for Maps. Cartogr. J. 2003, 40, 27–37. [Google Scholar] [CrossRef]
  33. Zeileis, A.; Hornik, K.; Murrell, P. Escaping RGBland: Selecting Colors for Statistical Graphics. Comput. Stat. Data Anal. 2009, 53, 3259–3270. [Google Scholar] [CrossRef] [Green Version]
  34. The R Project for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 28 August 2020).
  35. inlmisc: Miscellaneous Functions for the USGS INL Project Office. Available online: https://cran.r-project.org/web/packages/inlmisc/index.html (accessed on 4 September 2020).
  36. TOL Color Schemes. Available online: https://waterdata.usgs.gov/blog/tolcolors/ (accessed on 4 September 2020).
  37. Jenks, G.F.; Caspall, F.C. Error on Choroplethic Maps: Definition, Measurement, ReductioN. Ann. Assoc. Am. Geogr. 1971, 61, 217–244. [Google Scholar] [CrossRef]
  38. Geospatial Analysis 6th Edition, 2020 Update. Available online: https://www.spatialanalysisonline.com/HTML/index.html (accessed on 4 June 2021).
  39. McFayden, C.B.; Woolford, D.G.; Stacey, A.; Boychuk, D.; Johnston, J.M.; Wheatley, M.J.; Martell, D. L Risk Assessment for Wildland Fire Aerial Detection Patrol Route Planning in Ontario, Canada. Int. J. Wildland Fire. 2020, 29, 28–41. [Google Scholar] [CrossRef]
  40. McFayden, C.B.; Boychuk, D.; Woolford, D.G.; Wheatley, M.J.; Johnston, L. Impacts of Wildland Fire Effects on Resources and Assets through Expert Elicitation to Support Fire Response Decisions. Int. J. Wildland Fire 2019, 28, 885–900. [Google Scholar] [CrossRef]
  41. McFayden, C.; Boychuk, D.; Stacey, A.; Woolford, D.G.; McLarty, D.; Leonard, D.; Evens, J.; Johnston, D.; McAlpine, R.S. Daily Spatial Response Objective and Speed Weight Indicators to Assist in Daily Wildland Fire Preparedness Planning for ONTARIO, Canada. (in prep)
  42. Albert-Green, A.; Braun, W.J.; Martell, D.L.; Woolford, D.G. Visualization Tools for Assessing the Markov Property: Sojourn Times in the Forest Fire Weather Index in Ontario. Environmetrics 2014, 25, 417–430. [Google Scholar] [CrossRef]
Figure 1. Examples of different displays of the Fire Weather Index (FWI) class used by the Ontario’s fire agency: (a) a public roadside sign; (b) an agency operational map showing interpolated, colour-coded FWI and the raw FWI values for each weather station.
Figure 1. Examples of different displays of the Fire Weather Index (FWI) class used by the Ontario’s fire agency: (a) a public roadside sign; (b) an agency operational map showing interpolated, colour-coded FWI and the raw FWI values for each weather station.
Fire 04 00050 g001
Figure 2. Illustration of the design options (rectangles), which are described in the text, and trade-offs (circle) that are made when categorizing and colouring data that are on a numerical scale of measure.
Figure 2. Illustration of the design options (rectangles), which are described in the text, and trade-offs (circle) that are made when categorizing and colouring data that are on a numerical scale of measure.
Fire 04 00050 g002
Figure 3. The approximately 91.9 million ha extent of Ontario’s fire management area and the 20 km × 20 km resolution of the fire occurrence prediction (FOP) grid for Ontario.
Figure 3. The approximately 91.9 million ha extent of Ontario’s fire management area and the 20 km × 20 km resolution of the fire occurrence prediction (FOP) grid for Ontario.
Fire 04 00050 g003
Figure 4. Empirical probability distributions of the non-zero, human- and lightning-caused fire occurrence predictions (FOPs) for Ontario, 15 May–31 August, for 2016–2018 (human-caused) and 1992–2006 (lightning-caused). FOPs for the two causes have distinctly different ranges and central tendencies because of the spatio-temporal processes of ignition. The data are transformed by log 10 for visualization here, but untransformed data are required for operational use.
Figure 4. Empirical probability distributions of the non-zero, human- and lightning-caused fire occurrence predictions (FOPs) for Ontario, 15 May–31 August, for 2016–2018 (human-caused) and 1992–2006 (lightning-caused). FOPs for the two causes have distinctly different ranges and central tendencies because of the spatio-temporal processes of ignition. The data are transformed by log 10 for visualization here, but untransformed data are required for operational use.
Fire 04 00050 g004
Figure 5. The human-caused fire occurrence prediction of 10 June 2018, mapped using the categorization and colouring schemes indicated in the legends. The black dots are the human-caused fires reported later that day. The scales range from zero to 3 fires/cell: (a) the baseline scheme (Table 3), which shows little useful information; (b) 4 + 1 categories; (c) 10 + 1 categories; (d) 20 + 1 categories. The additional category in (bd) are for a “no forecast model” or true zero. Using more categories and hues greatly increases the information portrayed but makes matching the colour to the legend more difficult. The spatial pattern in (d) corresponds with roads and settlements.
Figure 5. The human-caused fire occurrence prediction of 10 June 2018, mapped using the categorization and colouring schemes indicated in the legends. The black dots are the human-caused fires reported later that day. The scales range from zero to 3 fires/cell: (a) the baseline scheme (Table 3), which shows little useful information; (b) 4 + 1 categories; (c) 10 + 1 categories; (d) 20 + 1 categories. The additional category in (bd) are for a “no forecast model” or true zero. Using more categories and hues greatly increases the information portrayed but makes matching the colour to the legend more difficult. The spatial pattern in (d) corresponds with roads and settlements.
Fire 04 00050 g005
Figure 6. The original and unified fire occurrence prediction (FOP) classification scales. Because of order-of-magnitude differences in FOP, a separate scale was originally used for each cause: linear for lightning (black) and subjectively determined, irregular, nonlinear for human (blue). The unified scale (red) is a systematic, nonlinear one generated by a power function; a cube root in this example.
Figure 6. The original and unified fire occurrence prediction (FOP) classification scales. Because of order-of-magnitude differences in FOP, a separate scale was originally used for each cause: linear for lightning (black) and subjectively determined, irregular, nonlinear for human (blue). The unified scale (red) is a systematic, nonlinear one generated by a power function; a cube root in this example.
Fire 04 00050 g006
Figure 7. Examples of fire occurrence prediction maps that use the categorization and colouring in Table 6: (a) human-caused; (b) lightning-caused; (c) total. The legend uses units of fires/cell and cells/fire to make the magnitudes integer. Larger versions are in the Supplementary Material.
Figure 7. Examples of fire occurrence prediction maps that use the categorization and colouring in Table 6: (a) human-caused; (b) lightning-caused; (c) total. The legend uses units of fires/cell and cells/fire to make the magnitudes integer. Larger versions are in the Supplementary Material.
Fire 04 00050 g007
Table 1. Current danger classes and colour scheme [15] used in Ontario for each of the six Fire Weather Index System components: Fine Fuel Moisture Code (FFMC), Duff Moisture Code (DMC), Drought Code (DC), Buildup Index (BUI), Initial Spread Index (ISI) and Fire Weather Index (FWI).
Table 1. Current danger classes and colour scheme [15] used in Ontario for each of the six Fire Weather Index System components: Fine Fuel Moisture Code (FFMC), Duff Moisture Code (DMC), Drought Code (DC), Buildup Index (BUI), Initial Spread Index (ISI) and Fire Weather Index (FWI).
ClassColourFFMCDMCDCISIBUIFWI
Low 0–800–150–1400–2.20–200–3
Moderate 81–8616–30141–2402.3–5.012–364–10
High 87–9031–50241–3405.1–10.037–6011–22
Extreme ≥91≥51≥340≥10≥61≥23
Table 2. Examples of the diverse classification of fire intensity and colouring used in Canada. The Ontario, Alberta and national schemes are used in daily maps. The British Columbia (BC) scheme is used in a static map of the 90th percentile of historical fire intensity. The Field Guide to the Fire Behaviour Prediction System (Field Guide) scheme is used in printed tables. Intensity classes IC 1–IC 6 are as defined in [19]; the higher classes are informal. The colours are approximate.
Table 2. Examples of the diverse classification of fire intensity and colouring used in Canada. The Ontario, Alberta and national schemes are used in daily maps. The British Columbia (BC) scheme is used in a static map of the 90th percentile of historical fire intensity. The Field Guide to the Fire Behaviour Prediction System (Field Guide) scheme is used in printed tables. Intensity classes IC 1–IC 6 are as defined in [19]; the higher classes are informal. The colours are approximate.
Intensity Class (IC)Range (kW/m) 1Ontario MapAlberta Map [23]Field Guide 2 [19]National Map [24]BC 90th Percentile Map 3 [25]
IC 1<10 <1 k
IC 210–500 <1 k
IC 3500–2 k <1 k1 k–2 k
IC 42 k–4 k
IC 54 k–10 k>4 k 4 k–6 k6 k–10 k
IC 6>10 k 10 k–30 k10 k–18 k
7 >30 k18 k–30 k
8 30 k–60 k
9 60 k–100 k
10 >100 k
1 Range applies to the full row except where overridden by the text in the cell. 2 The colours are limited by the printing ink used. 3 BC’s IC 1–5 corresponds to BC’s colour progression rather than the row label.
Table 3. The baseline case alternatives for the design options. This is a starting point that presents complete and accurate information, while ignoring the goals of a quick and easy understanding of the information.
Table 3. The baseline case alternatives for the design options. This is a starting point that presents complete and accurate information, while ignoring the goals of a quick and easy understanding of the information.
Design OptionAlternativeDescription
1. Gradient continuityContinuousNo categorization; the percentage of grey saturation is proportional to the datum magnitude
2. Hue selectionSingle sequentialAchromatic gradient of white through greys to black; no chromatic psychological colour associations
3. Range completenessCompleteFull range of data magnitude is shown; no truncation or aggregation at the top or bottom
4. Scale linearityLinearNo distortion; the percentage of grey saturation equals the datum’s position within its range
Table 4. Summary statistics of historical fire occurrence prediction model data (expected number of fires/cell) generated from 15 May–31 August for 2016–2018 (human-caused) and 1992–2006 (lightning-caused), with zeros removed.
Table 4. Summary statistics of historical fire occurrence prediction model data (expected number of fires/cell) generated from 15 May–31 August for 2016–2018 (human-caused) and 1992–2006 (lightning-caused), with zeros removed.
Number of ObservationsFirst QuartileMedianMeanThird QuartileMaximum
Human607,3830.000030.000130.001070.000481.37168
Lightning1,026,8600.000900.002600.009550.007503.89129
Table 5. Tabulation of how alternatives for the number of categories, amount of truncation and nonlinearity and number of hues generally interact in affecting several attributes related to completeness and accuracy of information or speed and ease of understanding. There are exceptions for some combinations and boundary conditions. Every alternative improves some attributes (blue-grey shading) and worsens others (orange-tan shading).
Table 5. Tabulation of how alternatives for the number of categories, amount of truncation and nonlinearity and number of hues generally interact in affecting several attributes related to completeness and accuracy of information or speed and ease of understanding. There are exceptions for some combinations and boundary conditions. Every alternative improves some attributes (blue-grey shading) and worsens others (orange-tan shading).
GoalAttributeAlternative for Design Options
More CategoriesMore Truncation at High endMore NonlinearityMore Hues
Complete and accurate informationMore resolution at low endBetterBetterBetterNo effect
More resolution at high endBetterBetter in rangeWorseNo effect
NONE beyond truncation
Accurate, UndistortedBetterBetter in rangeWorseNo effect
ABSENT beyond truncation
Fast and easy to understandNeed to refer to legendWorseNo effect in rangeWorseWorse
Ease of matching legendWorseBetterWorseBetter
Attention drawn to important dataNo effectNo effect in rangeNo effectBetter
ABSENT beyond truncation
Has suitable psychological meaningNo effectNo effectNo effectBetter if using colour psychology
Table 6. Categories and colours used for mapping fire occurrence prediction model outputs. Note that the map legend has a highly simplified integer scale and broad adjective categories.
Table 6. Categories and colours used for mapping fire occurrence prediction model outputs. Note that the map legend has a highly simplified integer scale and broad adjective categories.
Category NumberCategory Upper Bound of Expected NumberColourAdjective CategoryMap Legend
(Fires/Cell)(Cells/Fire)
20≥3.00≤0.333 Extreme Fire 04 00050 i001
192.570.389
182.190.457 Very High
171.840.543
161.540.651
151.270.790
141.030.972 High
130.8241.21
120.6481.54
110.4992.00
100.3752.67
90.2733.66
80.1925.21 Moderate
70.1297.77
68.10 × 10−212.3
54.69 × 10−221.3
42.40 × 10−241.7 Low
31.01 × 10−298.8
23.00 × 10−3333 Very Low
13.75 × 10−42670
00 Nil
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Boychuk, D.; McFayden, C.B.; Woolford, D.G.; Wotton, M.; Stacey, A.; Evens, J.; Hanes, C.C.; Wheatley, M. Considerations for Categorizing and Visualizing Numerical Information: A Case Study of Fire Occurrence Prediction Models in the Province of Ontario, Canada. Fire 2021, 4, 50. https://doi.org/10.3390/fire4030050

AMA Style

Boychuk D, McFayden CB, Woolford DG, Wotton M, Stacey A, Evens J, Hanes CC, Wheatley M. Considerations for Categorizing and Visualizing Numerical Information: A Case Study of Fire Occurrence Prediction Models in the Province of Ontario, Canada. Fire. 2021; 4(3):50. https://doi.org/10.3390/fire4030050

Chicago/Turabian Style

Boychuk, Den, Colin B. McFayden, Douglas G. Woolford, Mike Wotton, Aaron Stacey, Jordan Evens, Chelene C. Hanes, and Melanie Wheatley. 2021. "Considerations for Categorizing and Visualizing Numerical Information: A Case Study of Fire Occurrence Prediction Models in the Province of Ontario, Canada" Fire 4, no. 3: 50. https://doi.org/10.3390/fire4030050

APA Style

Boychuk, D., McFayden, C. B., Woolford, D. G., Wotton, M., Stacey, A., Evens, J., Hanes, C. C., & Wheatley, M. (2021). Considerations for Categorizing and Visualizing Numerical Information: A Case Study of Fire Occurrence Prediction Models in the Province of Ontario, Canada. Fire, 4(3), 50. https://doi.org/10.3390/fire4030050

Article Metrics

Back to TopTop