Subpixel Analysis of Primary and Secondary Infrared Emitters with Nighttime VIIRS Data

: Biomass burning is a coupled exothermic/endothermic system that transfers carbon in several subtractive analysis is the key addition to VNF version 4, which will commence production later in 2021. In 2022, we will seek validation of the VNF v.4 from nighttime Landsat and other data sources.


Remote Sensing Background
Planck's Law [1] and its derivatives provide the theoretical basis for calculating temperatures, source areas and radiant heat for subpixel infrared (IR) emitters from remotely observed radiances. If an object completely fills a sensor's field of view, the temperature can be calculated from the radiance collected in a single spectral band, assuming blackbody behavior. This is the principal behind infrared brightness temperatures [2]. The situation is more complicated for kilometer-scale pixels collected by meteorological satellites, where fires and flares are nearly always subpixel. In this case, the vast majority of the pixel footprint is composed of a cooler background, such as the land surface, sea, and cloud. By assuming the subpixel IR emitters all have the same temperature, multispectral radiances can be used to model the Planck curves for a single IR emitter and background. The temperature of the IR emitter is calculated from the emitter's Planck curve using Wein's Displacement Law [3]. Subpixel IR emitters are "gray-bodies" because they cover only a small portion of the pixel footprint. The IR emitter source area can be estimated by multiplying the pixel's footprint area (Figure 1) multiplied by the "emission scaling factor" or esf, which is the ratio between the IR emitter's Planck curve amplitude and the full pixel Planck curve for an object at that temperature. The esf is a number between 0 and 1. Figure 2 shows the esf calculation graphically for a subpixel infrared emitter at 500 K. The subpixel emitter's Planck curve has the same shape as the full pixel Planck curve for a 500 K object. The emission scaling factor is the ratio between the observed Planck curve amplitude and the amplitude of the Planck curve for a full pixel having that temperature. With temperature and source area as inputs, the radiant heat is calculated with the Stefan-Boltzmann Law [4].
Fire 2021, 4, x FOR PEER REVIEW 2 of 32 situation is more complicated for kilometer-scale pixels collected by meteorological satellites, where fires and flares are nearly always subpixel. In this case, the vast majority of the pixel footprint is composed of a cooler background, such as the land surface, sea, and cloud. By assuming the subpixel IR emitters all have the same temperature, multispectral radiances can be used to model the Planck curves for a single IR emitter and background. The temperature of the IR emitter is calculated from the emitter's Planck curve using Wein's Displacement Law [3]. Subpixel IR emitters are "gray-bodies" because they cover only a small portion of the pixel footprint. The IR emitter source area can be estimated by multiplying the pixel's footprint area (Figure 1) multiplied by the "emission scaling factor" or esf, which is the ratio between the IR emitter's Planck curve amplitude and the full pixel Planck curve for an object at that temperature. The esf is a number between 0 and 1. Figure 2 shows the esf calculation graphically for a subpixel infrared emitter at 500 K. The subpixel emitter's Planck curve has the same shape as the full pixel Planck curve for a 500 K object. The emission scaling factor is the ratio between the observed Planck curve amplitude and the amplitude of the Planck curve for a full pixel having that temperature.
With temperature and source area as inputs, the radiant heat is calculated with the Stefan-Boltzmann Law [4].  In the early 1990s, this concept was applied to bi-spectral (MWIR and LWIR) data from the Advanced Very High-Resolution Radiometer (AVHRR) [5]. More recently, the VIIRS nightfire (VNF) algorithm expanded the suite of spectral bands used in calculating temperatures of combustion sources from space into the near infrared and short-wave infrared [6]. VNF uses radiances observed in up to nine spectral bands spanning the NIR to the LWIR. The NIR and SWIR band data VIIRS collects at night are particularly useful for fire detection. These are daytime channels that continue to collect at night. With sunlight eliminated, fire pixels stand out against the sensor's noise floor. The fire pixels' NIR and SWIR radiances can be fully attributed to the IR emitters present. VNF uses simultaneous dual Planck curve fitting to calculate temperatures and source areas for an IR emitter and a cool background. In the early 1990s, this concept was applied to bi-spectral (MWIR and LWIR) data from the Advanced Very High-Resolution Radiometer (AVHRR) [5]. More recently, the VIIRS nightfire (VNF) algorithm expanded the suite of spectral bands used in calculating temperatures of combustion sources from space into the near infrared and short-wave infrared [6]. VNF uses radiances observed in up to nine spectral bands spanning the NIR to the LWIR. The NIR and SWIR band data VIIRS collects at night are particularly useful for fire detection. These are daytime channels that continue to collect at night. With sunlight eliminated, fire pixels stand out against the sensor's noise floor. The fire pixels' NIR and SWIR radiances can be fully attributed to the IR emitters present. VNF uses simultaneous dual Planck curve fitting to calculate temperatures and source areas for an IR emitter and a cool background.
These considerations are moot for the satellite fire products that rely on a single midwave infrared (MWIR) channel for the detection of "hotspot" pixels containing IR emitters [7,8]. Here, the emitter's Planck curve is sampled along with the background emissions at a single MWIR wavelength, making it impossible to calculate the temperature of subpixel fires or flares. The VIIRS and MODIS operational fire products report fire radiative power (FRP), which has units of megawatts [9]. These are the same units as radiant heat, which is calculated based on fire size and temperature with the Stefan-Boltzmann Law. Because the MWIR fire detection algorithms are unable to calculate temperature, FRP typically assumes a temperature of 1000 K and relies on an empirical calibration derived from radiative output observations of fires where the mass lost to combustion has been measured [10]. The FRP temperature assumption results in under and overestimates of heat output in pixel where temperatures deviate from 1000 K. As a result, a high temperature version of FRP was developed for estimating the heat output from natural gas flares [11], which These considerations are moot for the satellite fire products that rely on a single mid-wave infrared (MWIR) channel for the detection of "hotspot" pixels containing IR emitters [7,8]. Here, the emitter's Planck curve is sampled along with the background emissions at a single MWIR wavelength, making it impossible to calculate the temperature of subpixel fires or flares. The VIIRS and MODIS operational fire products report fire radiative power (FRP), which has units of megawatts [9]. These are the same units as radiant heat, which is calculated based on fire size and temperature with the Stefan-Boltzmann Law. Because the MWIR fire detection algorithms are unable to calculate temperature, FRP typically assumes a temperature of 1000 K and relies on an empirical calibration derived from radiative output observations of fires where the mass lost to combustion has been measured [10]. The FRP temperature assumption results in under and overestimates of heat output in pixel where temperatures deviate from 1000 K. As a result, a high temperature version of FRP was developed for estimating the heat output from natural gas flares [11], which are nearly twice as hot as biomass burning. Again, an assumption is made regarding the temperature of the flare.
Is it possible to resolve more than a single temperature phase within individual pixel footprints? The answer to that question depends on the temperature separation of the phases, their source sizes and spectral bands available. If the temperatures of the phases are separated far enough, their Planck curves are offset from each other, with certain spectral ranges dominated by a particular phase. Planck curves for hotter sources are shifted towards shorter wavelengths and cooler sources are shifted to longer wavelengths ( Figure 3). Thus, a sensor with wide separation of spectral bands may produce data suitable for the analysis of multiple temperature phases present within individual pixels. phases, their source sizes and spectral bands available. If the temperatures of the phases are separated far enough, their Planck curves are offset from each other, with certain spectral ranges dominated by a particular phase. Planck curves for hotter sources are shifted towards shorter wavelengths and cooler sources are shifted to longer wavelengths ( Figure  3). Thus, a sensor with wide separation of spectral bands may produce data suitable for the analysis of multiple temperature phases present within individual pixels. For several years, we have been researching the subpixel analysis of two IR emitters based on the temperature differences between the phases. The initial development was conducted with nighttime Landsat data collected on smoldering peatland fires in Sumatra [12]. The initial indication that two fire temperature phases were present was that the SWIR fire features appeared visually different from the LWIR thermal anomalies. The SWIR anomalies tend to be brighter on the outer edges of the larger and more uniform LWIR features. In addition, patches of SWIR anomalies were found that lacked LWIR anomaly. This led to the "tip-and-tail" strategy for deriving flaming and non-flaming temperatures. The flaming-phase temperatures and source areas are estimated via Planck curve fitting with the NIR and SWIR radiances, with resulting temperature in the 780-980 K range. Conversely, the temperatures of heated soil surfaces with smoldering present were estimated as 320-455 K via Planck curve fitting with the two LWIR radiances. The LWIR channels are located on the trailing edge (tail) of the non-flaming-phase radiant emissions. For several years, we have been researching the subpixel analysis of two IR emitters based on the temperature differences between the phases. The initial development was conducted with nighttime Landsat data collected on smoldering peatland fires in Sumatra [12]. The initial indication that two fire temperature phases were present was that the SWIR fire features appeared visually different from the LWIR thermal anomalies. The SWIR anomalies tend to be brighter on the outer edges of the larger and more uniform LWIR features. In addition, patches of SWIR anomalies were found that lacked LWIR anomaly. This led to the "tip-and-tail" strategy for deriving flaming and non-flaming temperatures. The flaming-phase temperatures and source areas are estimated via Planck curve fitting with the NIR and SWIR radiances, with resulting temperature in the 780-980 K range. Conversely, the temperatures of heated soil surfaces with smoldering present were estimated as 320-455 K via Planck curve fitting with the two LWIR radiances. The LWIR channels are located on the trailing edge (tail) of the non-flaming-phase radiant emissions.
The Landsat tip-and-tail approach led to the development of the "flaming subtractive method" for subpixel analysis of up to three temperature phases present in VIIRS fire pixels from a peatland fire in Sumatra [13]. VIIRS has nearly the same SWIR and LWIR bands as Landsat, but also has two mid-wave infrared (MWIR) channels at 3.7 and 4.05 um ( Table 1). In the flaming subtractive method, the fire pixel's NIR and SWIR radiances are used to model the primary emitter Planck curves, yielding temperature, source sizes, plus the primary emitter radiances in all spectral bands. In this case, the primary emitter temperatures were in the 750-1400 K range. The primary emitter radiances are subtracted from the original radiances in all nine spectral bands. Then, dual Planck curve fitting is applied to the residual radiances to model Planck curves for a secondary emitter plus the background. The resulting temperatures and source sizes are then filtered to remove pixels yielding implausible results. The rejected misfit pixels revert back to the standard VIIRS nightfire analysis, with a single IR emitter and background. The Sumatra secondary emitters fell into two temperature ranges: one from 650 to 750 K and the other from 320 to 450 K.

Temperature Ranges of Fire Pixel Components
Biomass burning is a coupled exothermic/endothermic process that consumes biomass, leaving a mineral ash. The exothermic component is combustion, which produces the heat that drives the endothermic components, which prepares the fuels for combustion. The endothermic phases include preheating and pyrolysis. Preheating drives water from the fuel, preparing it for pyrolysis. In pyrolysis, the structural components of biomass, such as cellulose, hemicellulose and lignin are broken done to volatile gases (pyrolyzate) and char, leaving an ash residue. Char is the fuel of the smoldering phase and pyrolyzate is the fuel of the flaming phase. Flaming is the open air oxidation of volatile gases produced from pyrolysis. Smoldering is the oxidation of the char produced from pyrolysis. Smoldering has a characteristic orange glow present on the char surface. Smoke is the unoxidized fire components that escape to the air [14]. When all the carbon is gone, what remains is mineral ash-primarily oxides of calcium, sodium, and potassium [15]. The basic outlines of pyrolysis, smoldering and flaming are shown in Equations (1) Figure 4 show the temperature ranges for the basic components that may be present in a VIIRS fire pixel. Background, which normally occupies the largest part of the pixel footprint, includes clouds (162-300 K), sea (270-308 K) and land (248-318 K) [17][18][19]. The temperature ranges of the endothermic components of fire pixels include preheating from 323-453 K [14][15][16], slow pyrolysis from 453-673 K and fast pyrolysis from 773-1173 [14,20]. Smoldering temperatures range from 723-973 K [16]. Flaming biomass ranges from 923-1173 K [16,21]. Natural gas flaring is generally hotter than biomass flaming, with temperatures ranging from 1300-2000 K [22]. in a VIIRS fire pixel. Background, which normally occupies the largest part of the pixel footprint, includes clouds (162-300 K), sea (270-308 K) and land (248-318 K) [17][18][19]. The temperature ranges of the endothermic components of fire pixels include preheating from 323-453 K [14][15][16], slow pyrolysis from 453-673 K and fast pyrolysis from 773-1173 [14,20]. Smoldering temperatures range from 723-973 K [16]. Flaming biomass ranges from 923-1173 K [16,21]. Natural gas flaring is generally hotter than biomass flaming, with temperatures ranging from 1300-2000 K [22].

Study Objectives
To date, the flaming subtractive method has only been applied to a single orbit of VIIRS nighttime data, collected over peatland burning in Sumatra. The objective of this study is to evaluate its performance in temperate zone forest fires and natural gas flaring. In addition, because of the large number of misfit pixels found in the Sumatra testing [13], we decided it is important to test additional approaches to the unmixing of primary and secondary emitters from background. In this paper we report on an intercomparison of five methods for subpixel analysis of primary and smoldering combustion in the presence of background using nighttime VIIRS data collected on California mega-fires from 2018 and large natural gas flares present in southern Iraq. In addition, the temperature ranges of the emitters are used to identify the probable fire pixel components responsible for the radiant emissions based on Figure 4.

Study Objectives
To date, the flaming subtractive method has only been applied to a single orbit of VIIRS nighttime data, collected over peatland burning in Sumatra. The objective of this study is to evaluate its performance in temperate zone forest fires and natural gas flaring. In addition, because of the large number of misfit pixels found in the Sumatra testing [13], we decided it is important to test additional approaches to the unmixing of primary and secondary emitters from background. In this paper we report on an intercomparison of five methods for subpixel analysis of primary and smoldering combustion in the presence of background using nighttime VIIRS data collected on California mega-fires from 2018 and large natural gas flares present in southern Iraq. In addition, the temperature ranges of the emitters are used to identify the probable fire pixel components responsible for the radiant emissions based on Figure 4.

Materials and Methods
We tested and rated five different algorithms for separating primary and secondary emitters in three test sites listed in Table 2. The test sites include the Mendocino and Carr Fires in California plus a set of large natural gas flares in Iraq. Figures 5-7 show the set of spectral bands analyzed. Note that in all three test areas the fires and flares show up most clearly in the two SWIR bands. The Mendocino and Carr Fire detections gradually fade and partially disappear moving to either shorter or longer wavelengths. Both the Mendocino and Carr Fires show obvious fire fronts, indicating the fire is spreading outward into unburnt fuels. The Basra gas flares are still clearly present in the NIR bands, but progressively fade as in the LWIR channels.  Figure 5. The Mendocino fire study area. The fire is north and east of Clear Lake, which is marked with the letters "CL" on the M12 image.

VNF Detection Algorithms
VNF uses two types detection algorithms: one for the NIR and SWIR bands (M7, 8, 10, 11) and a second detector for the MWIR (M12-M13). The NIR and SWIR band are daytime channels, with detection limits and saturation levels set for recording reflected sunlight. At night, they primarily record the noise floor of the system, which is occasionally punctuated by high radiant emissions from fires or flares. Here, we set a detection threshold as the image mean plus four standard deviations. The MWIR signal is a mixture of radiant emissions from clouds, the earth's surface, and IR emitters. The VNF MWIR thermal anomaly detector relies on the fact that M12 and M13 radiances for background (land, sea, and clouds) are highly correlated, forming a dense diagonal on M12 versus M13 scattergrams ( Figure 8). The algorithm generates the MWIR scattergram for a granule or an aggregate, locates the diagonal, and detects pixels with thermal anomalies as outliers pulled away from the diagonal.

VNF Detection Algorithms
VNF uses two types detection algorithms: one for the NIR and SWIR bands (M7, 8, 10, 11) and a second detector for the MWIR (M12-M13). The NIR and SWIR band are daytime channels, with detection limits and saturation levels set for recording reflected sunlight. At night, they primarily record the noise floor of the system, which is occasionally punctuated by high radiant emissions from fires or flares. Here, we set a detection threshold as the image mean plus four standard deviations. The MWIR signal is a mixture of radiant emissions from clouds, the earth's surface, and IR emitters. The VNF MWIR thermal anomaly detector relies on the fact that M12 and M13 radiances for background (land, sea, and clouds) are highly correlated, forming a dense diagonal on M12 versus M13 scattergrams (Figure 8). The algorithm generates the MWIR scattergram for a granule or an aggregate, locates the diagonal, and detects pixels with thermal anomalies as outliers pulled away from the diagonal. . VIIRS has two mid-wave infrared (MWIR) channels that are closely spaced near 4 um. The scattergram shows M12 versus M13 radiances. At night, there is a prominent diagonal data cloud representing the temperature variations in the background. The presence of subpixel IR emitters pulls pixels away from the diagonal. The VNF MWIR detector locates the diagonal, draws an envelope around it, and labels the pixels outside the envelope as detections.

Types of VNF Detections
The VNF pixels are sorted into six types based on the spectral bands having detection.  . VIIRS has two mid-wave infrared (MWIR) channels that are closely spaced near 4 um. The scattergram shows M12 versus M13 radiances. At night, there is a prominent diagonal data cloud representing the temperature variations in the background. The presence of subpixel IR emitters pulls pixels away from the diagonal. The VNF MWIR detector locates the diagonal, draws an envelope around it, and labels the pixels outside the envelope as detections.

Types of VNF Detections
The VNF pixels are sorted into six types based on the spectral bands having detection.

Planck Curve Fitting
Planck curve fitting is accomplished using a simplex algorithm [23] to optimize the fit to the radiance data. The fitting starts with initial temperatures of 1000 K for primary emitter, 500 K for the secondary emitter and 300 K for background. The simplex algorithm iteratively modifies the temperature and emission scaling factor (esf) to minimize the residual radiance. The original VNF algorithm uses the fitting procedure to model Planck curves for a single IR emitter phase and a background. This is referred to as "dual-curve" fitting and is applied to the Type 1-3 detections. In the flaming subtractive processing, Type 4 detections are analyzed for three temperature phases: primary and secondary emitters plus the background. These pixels are then filtered to remove those having spurious results, the so called "misfits" (Figure 9). The misfit pixels are converted to Type 5 and are reprocessed with dual Planck curve fitting, solving for a single IR emitter and background.
Type 4-Has detection in two SWIR bands and MWIR. May have NIR detection as well. The primary and secondary emitter analysis is restricted to Type 4 detections. Type 5-Pixels that yield spurious results revert back to the original dual Planck curve fitting with a single IR emitter plus background.

Planck Curve Fitting
Planck curve fitting is accomplished using a simplex algorithm [23] to optimize the fit to the radiance data. The fitting starts with initial temperatures of 1000 K for primary emitter, 500 K for the secondary emitter and 300 K for background. The simplex algorithm iteratively modifies the temperature and emission scaling factor (esf) to minimize the residual radiance. The original VNF algorithm uses the fitting procedure to model Planck curves for a single IR emitter phase and a background. This is referred to as "dual-curve" fitting and is applied to the Type 1-3 detections. In the flaming subtractive processing, Type 4 detections are analyzed for three temperature phases: primary and secondary emitters plus the background. These pixels are then filtered to remove those having spurious results, the so called "misfits" (Figure 9). The misfit pixels are converted to Type 5 and are reprocessed with dual Planck curve fitting, solving for a single IR emitter and background. Figure 9. Scattergram of the unfiltered 2nd phase versus background temperatures from the three test sites. Secondary emitter temperatures below 350 K are associated with extremely large source areas fed by the incorporation of background radiance. The background size and temperature are forced down due to the loss of radiance to the secondary emitter. The misfit pixels are filtered out, relabeled as Type 5 and analyzed for a single IR emitter plus background. Figure 9. Scattergram of the unfiltered 2nd phase versus background temperatures from the three test sites. Secondary emitter temperatures below 350 K are associated with extremely large source areas fed by the incorporation of background radiance. The background size and temperature are forced down due to the loss of radiance to the secondary emitter. The misfit pixels are filtered out, relabeled as Type 5 and analyzed for a single IR emitter plus background.

Saturation
It is important to exclude radiances affected by saturation in the final Planck curve fitting. The only two spectral bands used in VNF where saturation can be found are M11 and M12 [24]. The general effect of saturation is to reduce the pixel radiance relative to the unaffected bands. Saturation can occur anywhere in the swath. In the nadir aggregation zone three detectors are averaged to make a pixel. Thus, if one of the three is saturated, the result is a radiance less than full-pixel saturation, corrupted by inclusion of a saturated detector radiance. This is also the case for the mid-swath aggregation zone where the signal from two detectors are averaged. Typically, only pixels from the outer aggregation zone reach the saturation radiances. Thus, there are two styles of saturation: "full-pixel saturation" and "subpixel saturation". Subpixel saturation does not occur in the outer aggregation zone. Full-pixel saturation is unambiguous because the recorded radiance matches the saturation radiances. To identify the presence of subpixel saturation, we look for anomalously low radiances in M11 and M13. For M11, we fit a Planck curve to the M7 -M11 radiances and compared the observed M11 radiance to that predicted from the Planck curve fit. Unsaturated M11 pixels fall on a diagonal in the M11 versus predicted M11 radiance scattergram. Pixels having subpixel M11 saturation are pulled away from the diagonal, with the predicted radiance higher than the observed radiance. For M12 VNF first runs Planck curve fitting using the radiances in all the M bands, reporting out the sum of the squared residuals (SSR). The presence of saturation or other radiance irregularity is indicated by SSR values greater than 2 watts/m2.sr.um. In this case, M11 and M12 are both tested for the presence of subpixel saturation by removing their radiances from Planck curve fitting to check if the SSR drops below 2 watts. The testing is conducted three ways, removal of the M11 radiance, removal of the M12 radiances, and removal of both M11 and M12 radiances. This testing makes it possible to discern whether the subpixel saturation is present in M11, M12, or both spectral bands. The final Planck curve fitting is performed without the spectral bands affected by saturation.

Atmospheric Correction
An atmospheric correction was developed for the three test areas using Libradtran [25] and applied to the M band radiances of the Type 4 detections. Libradtran was parameterized based on the terrain elevation, atmospheric path length, and atmospheric profiles of water vapor and temperature from simultaneously collected Advanced Technology Microwave Sounder (ATMS) data processed with the Microwave-Integrated Retrieval System (MIRS) [26].
The purpose of atmospheric correction is to derive the surface-emitted radiance Ls from the satellite observed Top-of-Atmosphere (TOA) radiance Lt. The Lt is a measurement of the combined surface-emitted radiance (Ls) and atmospheric up-(Lu) and downwelling (Ld) radiances transmitted through the atmosphere (τ): where ε is the target emissivity, and 1 − ε is the target albedo. To convert Lt to surface temperature, the Lu, Ld, τ, and ε need to be estimated. As suggested in [27], Lu and τ can be determined from Libradtran radiative transfer simulations for ε = 1 at two different surface temperatures Ls1 and Ls2 as the offset and slope of a linear fit of the surface and TOA signals: With the knowledge of Lu and τ, Ld can then be computed from the third Libradtran simulation with ε < 1 and Ls = 0: Ls can then be estimated from Equation (4).

Flaming Subtractive Method
The NIR and SWIR radiances are used to model the Planck curve of the primary emitter, yielding a temperature, source size and primary emitter radiances in all spectral bands. Residual radiances are calculated based on the subtraction of the primary emitter radiances. Then, dual Planck curve fitting is used to define the secondary emitter and background. The result is a set of temperatures and source areas for three phase: primary emitter, secondary emitter, and background.

Triple Curve Methods
Simultaneous Planck curve fitting for primary and secondary emitter plus background. Four varieties have been tested: A.
Constrained by the primary emitter temperature calculated from Planck curve fitting of the NIR and SWIR detection radiances. C.
Constrained by the local background temperature. D.
Constrained by the primary emitter temperature and local background temperature.

Rating
The five methods were rated based on the number of "misfit" detections. Misfits are recognized based on implausible results. Examples include unrealistically low secondary emitter or background temperatures ( Figure 9). Results from each of the tested methods were inspected and criteria were developed to filter out misfits. Since the same set of Type 4 pixels were processed in each method, the number of misfits found can be used to rate the methods. The rating is in the form of percent yield. Table 3 shows the scoring of the five methods, listing misfit tallies by type and the yield calculation. The yield is calculated as the percentage of the Type 4 detections that yield valid secondary emitter results. As with previous tests conducted on a Sumatra peatland fire, the flaming subtractive algorithm has the highest yield. The rest of this paper focuses on exploring the flaming subtractive results.   Figures 10 and 11 are from the Mendocino Fire and Figure 12 is from a Basra gas flare.

Examples of Good Fits
In Figure 10, VNF V.3 the IR emitter, calculated from the NIR, SWIR and MWIR radiances came out with a temperature of 894 K and area of 211 m 2 . The flaming subtractive method successfully split the fire radiant energy into primary and secondary IR emitters. Using the NIR and SWIR radiances to model the Planck curve of the primary emitter results in a slightly hotter and smaller source size, with a temperature of 1237 K and 25 m 2 . The secondary emitter, calculated from the residual radiances after subtracting the flaming phase, has a temperature of 490 K and source size of 4458 m 2 .    1 Figure 11. Spectral diagrams of a Type 4 pixel from the Mendocino Fire where a secondary emitter has been detected. Note that M12 has subpixel saturation and has been excluded from Planck curve fitting.  Figure 12 shows a gas flare detection from Basra. The V.3 algorithm calculated the IR emitter temperature as 1690 K and source size of 7.65 m 2 . By modeling the primary emitter with the NIR and SWIR radiances the temperature rises slightly to 1704 K and source size drops slightly to 7.26 m 2 . The dual Planck curve fitting on the residual radiances found a secondary emitter at 399 K and 1400 m 2 .

Temperature Versus Source Area Scattergrams and Temperature Histograms
The original VNF (v.3) performs dual Planck curve fitting on radiances spanning NIR to LWIR, calculating a single subpixel infrared emitter temperature and source size in the  the background temperature down to 223 K-an unrealistic value given the clear sky conditions. This style of misfit can be traced to the extremely high secondary emitter esf. This pulls the secondary emitter temperature down to the background level and the source size expands to nearly the full pixel footprint. Figure 14, from a Basra flare, shows the other common type of misfit, traced to a negative esf for the secondary emitter. From V.3 the IR emitter temperature 1604 K and source area is 129 m 2 . When the flaming subtractive method is applied the primary emitter comes out as 1566 K and 150 m 2 . The pixel is classified as a misfit because the secondary emitter's esf is negative, resulting in a source area of −530 m 2 .
Another example of a VNF v.4 misfit is shown in Figure 14. Here, the secondary emitter's esf is highly negative, resulting in the background and secondary emitter's Planck curves forming mirror images of each other-both with temperatures of 292 K.

Temperature versus Source Area Scattergrams and Temperature Histograms
The original VNF (v.3) performs dual Planck curve fitting on radiances spanning NIR to LWIR, calculating a single subpixel infrared emitter temperature and source size in the presence of a cool background (clouds and earth surfaces). In this case, the IR emitter temperatures and source areas are dominated by the flaming phase due to the T 4 term in the Stefan-Boltzmann Law. Figure 16 shows the VNF v.3 scattergram of temperatures and source areas, plus the corresponding temperature histogram. VNF v.3 produces a bimodal distribution-with a hot peak associated with natural gas flares near 1600 K and a peak centered near 900 K arising primarily from biomass burning. presence of a cool background (clouds and earth surfaces). In this case, the IR emitter temperatures and source areas are dominated by the flaming phase due to the T 4 term in the Stefan-Boltzmann Law. Figure 16 shows the VNF v.3 scattergram of temperatures and source areas, plus the corresponding temperature histogram. VNF v.3 produces a bimodal distribution-with a hot peak associated with natural gas flares near 1600 K and a peak centered near 900 K arising primarily from biomass burning. The flaming subtractive method splits a pixel's multispectral radiances into three phases-primary emitter, secondary emitter and background. This is accomplished by fitting a Planck curve to the NIR and SWIR radiances to represent the primary emitter. The primary emitter radiances are calculated in each spectral band from the primary emitter  The flaming subtractive method splits a pixel's multispectral radiances into three phases-primary emitter, secondary emitter and background. This is accomplished by fitting a Planck curve to the NIR and SWIR radiances to represent the primary emitter. The primary emitter radiances are calculated in each spectral band from the primary emitter Planck curve and are subtracted from the original radiances. Then, dual Planck curve fitting is applied to the residual radiances to derive a secondary emitter phase and background. Figure 17 shows the temperature versus area scattergram and temperature histogram from the flaming subtractive method. Only the pixels having "fits" are shown. The prominent v.3 peak near 900 K is split into a cooler peak near 550 K and warmer peak near 1000 K. The gas flaring peak near 1600 K is diminished by the low number of gas flaring detections found to have a valid secondary emitter fits.  Figure 18 shows a temperature montage of the VNF Type 4 detections from the three test sites. The VNF v.3 temperature images are on the left side. In the middle are the primary emitter temperature images and on the right are the secondary emitter temperature images. Note that flaming subtractive procedure raises the temperatures of the primary emitter slightly relative to v.3. In contrast, the secondary emitter temperatures are generally 400-500 K cooler than the v.3 temperatures for the same pixel.  Figure 18 shows a temperature montage of the VNF Type 4 detections from the three test sites. The VNF v.3 temperature images are on the left side. In the middle are the primary emitter temperature images and on the right are the secondary emitter temperature images. Note that flaming subtractive procedure raises the temperatures of the primary emitter slightly relative to v.3. In contrast, the secondary emitter temperatures are generally 400-500 K cooler than the v.3 temperatures for the same pixel.  Figure 19 is a scattergram of the flaming subtractive derived temperatures with and without atmospheric correction. The atmospheric correction had little influence on the primary emitter temperatures. We believe this is attributable to high atmospheric transmissivity in the VIIRS SWIR spectral bands, the dry atmospheres of the test sites, and near nadir location of fires, which minimizes variation in atmospheric path length. Significantly, the atmospheric correction doubled the number of secondary emitter fits and shifted the temperatures to form a histogram peak at 550 K.  Figure 19 is a scattergram of the flaming subtractive derived temperatures with and without atmospheric correction. The atmospheric correction had little influence on the primary emitter temperatures. We believe this is attributable to high atmospheric transmissivity in the VIIRS SWIR spectral bands, the dry atmospheres of the test sites, and near nadir location of fires, which minimizes variation in atmospheric path length. Significantly, the atmospheric correction doubled the number of secondary emitter fits and shifted the temperatures to form a histogram peak at 550 K.

Discussions
The VNF v.3 temperatures, source areas, and radiant heat are dominated by the high temperature exogenic combustion components. This is the case because the exogenic components are substantially hotter than the endothermic components. The exothermic component's dominance of the v.3 product can be attribute to the exponential rise in radiant emissions due to the T 4 term in the Stefan-Boltzmann Law. The high degree of similarity between the VNF v.3 data and the primary emitter from the flaming subtractive analysis can be seen in Figure 20, which shows temperatures and radiant heat for the Type 4 pixels from the study set. Using the short wavelength radiances to model the primary emitter Planck curve results in a slight increase in temperature and a slight decrease in source size as compared to the v.3 results.

Discussions
The VNF v.3 temperatures, source areas, and radiant heat are dominated by the high temperature exogenic combustion components. This is the case because the exogenic components are substantially hotter than the endothermic components. The exothermic component's dominance of the v.3 product can be attribute to the exponential rise in radiant emissions due to the T 4 term in the Stefan-Boltzmann Law. The high degree of similarity between the VNF v.3 data and the primary emitter from the flaming subtractive analysis can be seen in Figure 20, which shows temperatures and radiant heat for the Type 4 pixels from the study set. Using the short wavelength radiances to model the primary emitter Planck curve results in a slight increase in temperature and a slight decrease in source size as compared to the v.3 results.
Regarding the identity of the primary and secondary emitter components, let us examine the temperature ranges of the phases relative to the temperature ranges for fire pixel components shown in Figure 4. The biomass burning primary emitter peaks at 950 K, but spans 800-1300 K ( Figure 17). This corresponds to the flaming phase, with possible contribution from smoldering and fast pyrolysis. The secondary emitter peaks at 550 K and spans 350-750 K, corresponding to the preheating and slow pyrolysis components.
In the two biomass burn test cases (Mendocino and Carr), the secondary emitter is linked to the radiant emissions from the endothermic portion of the biomass burning system. The preheating and slow-pyrolysis radiant emissions are hidden under the higher radiant output of the exothermic components. Subtraction of the primary emitter radiances makes it possible uncover the presence of the secondary emitter, which then can be split from the background radiance through dual Planck curve fitting. The flaming subtractive processing result is a substantial increase in the total radiant heat ( Figure 20) and total source area ( Figure 21). Regarding the identity of the primary and secondary emitter components, let us examine the temperature ranges of the phases relative to the temperature ranges for fire pixel components shown in Figure 4. The biomass burning primary emitter peaks at 950 K, but spans 800-1300 K ( Figure 17). This corresponds to the flaming phase, with possible contribution from smoldering and fast pyrolysis. The secondary emitter peaks at 550 K and spans 350-750 K, corresponding to the preheating and slow pyrolysis components.
In the two biomass burn test cases (Mendocino and Carr), the secondary emitter is linked to the radiant emissions from the endothermic portion of the biomass burning system. The preheating and slow-pyrolysis radiant emissions are hidden under the higher radiant output of the exothermic components. Subtraction of the primary emitter radiances makes it possible uncover the presence of the secondary emitter, which then can be split from the background radiance through dual Planck curve fitting. The flaming subtractive processing result is a substantial increase in the total radiant heat ( Figure 20) and total source area ( Figure 21).  Looking at Figures 5-7, one might conclude that full pixels are on fire. However, amongst this test set, none of the VNF pixel emission scaling factor's (esf) rose above 0.1 or 10% of the pixel footprints ( Figure 22).  Looking at Figures 5-7, one might conclude that full pixels are on fire. However, amongst this test set, none of the VNF pixel emission scaling factor's (esf) rose above 0.1 or 10% of the pixel footprints ( Figure 22).

Conclusions
We have developed a flaming subtractive method for subpixel unmixing of radiant emissions from two IR emitter temperature phases present in fire pixels at night. The primary emitter temperature and source area are calculated using the Planck curve defined based on the NIR and SWIR radiances. The primary emitters Planck curve is used to calculate primary emitter radiances in all spectral bands, which are subtracted to create residual radiances. The residual radiances are analyzed for a secondary IR emitter and a cool background. Spurious secondary emitter results are filtered out based on unrealistic esf, temperature and source area values. The failure of the fitting indicates that a secondary emitter could not be found. These "misfit" pixels revert back to the dual-curve fitting for a single IR emitter and background-the standard processing from VNF v.3. The flaming subtractive analysis is the key addition to VNF version 4.
We have now applied the split IR emitter analysis to both nighttime Landsat and nighttime VIIRS data in a total of five test sites. The Landsat analysis used the "tip and

Conclusions
We have developed a flaming subtractive method for subpixel unmixing of radiant emissions from two IR emitter temperature phases present in fire pixels at night. The primary emitter temperature and source area are calculated using the Planck curve defined based on the NIR and SWIR radiances. The primary emitters Planck curve is used to calculate primary emitter radiances in all spectral bands, which are subtracted to create residual radiances. The residual radiances are analyzed for a secondary IR emitter and a cool background. Spurious secondary emitter results are filtered out based on unrealistic esf, temperature and source area values. The failure of the fitting indicates that a secondary emitter could not be found. These "misfit" pixels revert back to the dual-curve fitting for a single IR emitter and background-the standard processing from VNF v.3. The flaming subtractive analysis is the key addition to VNF version 4.
We have now applied the split IR emitter analysis to both nighttime Landsat and nighttime VIIRS data in a total of five test sites. The Landsat analysis used the "tip and tail" method to model the primary emitter Planck curve with NIR and SWIR radiances and the secondary emitter Planck curve is derived from the LWIR radiances. The "tip and tail" method makes no attempt to separate secondary emitter and background radiant emissions. The flaming subtractive method explicitly attempts to unmix temperatures and source areas for three temperature phases present in fire pixels: primary emitter, secondary emitter and background. The flaming subtractive analysis makes use of nighttime VIIRS data collected from nine spectral bands in four spectral ranges: NIR, SWIR, MWIR and LWIR.
The temperature of the secondary emitter has varied between the five test sites. In the Landsat 8 study of a smoldering peatland fires, the secondary emitter temperature was 320-450 K, indicating pre-heating. This is the result of smoldering in the soil column heating the soil surface. In the VIIRS analysis of a smoldering peatland fire in Sumatra, two secondary emitter temperature ranges emerged-one from preheating (320-450 K) and one from pyrolysis/smoldering from 650 to 750 K. In the 2018 California megafires, the secondary emitter had temperatures from 350 to 650 K, covering the temperature ranges typical of pre-heating and slow pyrolysis. Several of the gas flares in Basra had 2nd phases but the temperatures were highly variable, spanning 350-800 K. The identity of the secondary emitter found at gas flaring sites will require further investigation.
Our results indicate that an atmospheric correction is quite important for the flaming subtractive analysis. The atmospheric correction had minimal impact on the flaming-phase analysis in these three test cases. This may be due to high transmissivity in the SWIR spectral bands in the dry atmospheres present in Northern California and Southern Iraq. However, adding the atmospheric corrected resulted in a doubling in the number of good 2nd-phase fits. In addition, inclusion of an atmospheric correction results in a~150 K decline in the 2nd-phase temperature.
Our plan is to implement global near-real-time processing with the flaming subtractive algorithm as VNF version 4 for both VIIRS instruments (SNPP and NOAA-20). Production of the VNF v.3 product will continue. Later we plan to reprocess the VIIRS archive to extend VNF v.4 record back to the beginning of the VIIRS archive. The M11 collections extend back to the beginning of 2018. The algorithm can still work if there is detection in at least two of the three short wavelength channels (M7, M8 and M10) collected in the first five years of VIIRS data. This will cut the number of VNF detections that can be processed with the VNF v. 4 algorithm from 2012 through 2017.
By splitting the primary and secondary emitters apart, there is an approximate doubling of the VNF calculated source area and radiant heat. It will be important to examine the validity of VNF v.4 primary and secondary emitter results using independent data sources. This could include maps of flaming fire fronts and smoldering combustion made by firefighters, airborne multispectral data collected on active fires, and nighttime Landsat, which has many of the same set of spectral bands as VIIRS, though missing the MWIR bands.
The capability to distinguish multiple components of biomass burning from space has the potential to improve the identification of core zones having the largest mass of combustion biomass during fire suppression efforts, modeling of smoke production, the proportion of partially oxidized greenhouse gas emissions, and the severity of the burning. While there are several factors that affect the proportion of the components, fuel size is a major determining factor. Flaming predominates in grassland fires, where the fuels have large surface areas and are readily combusted [28]. Forest fires tend to start with flaming-phase combustion in fine dry fuels, such as beds of dry pine needles. Heat from the flaming combustion eventually induces pyrolysis and combustion of larger wood branches and trunks. The standard satellite hotspot is unable to distinguish multiple temperature phases present in biomass burning. However, this can be accomplished at night if there are multiple spectral bands that span a wide range of wavelengths.