Fire Diurnal Cycle Derived from a Combination of the Himawari-8 and VIIRS Satellites to Improve Fire Emission Assessments in Southeast Australia

: The violent and persistent wildﬁres that broke out along the southeast coast of Australia in 2019 caused a large number of pollutant emissions, which seriously affected air quality and the global climate. The existing two methods for estimating combustion emissions based on burned area and ﬁre radiative power mainly use a medium resolution imaging spectrometer (MODIS) on the Aqua and Terra satellites. However, the low temporal resolution of MODIS and insensitivity to small ﬁres lead to deviation in the estimation of ﬁre emissions. In order to solve this problem, the Visible Infrared Imaging Radiometer Suite (VIIRS) with better performance is adopted in this paper, combined with the ﬁre diurnal cycle information obtained by geostationary satellite Himawari-8, to explore the spatio-temporal model of biomass combustion emissions. Using this, a high-spatial-and -temporal-resolution ﬁre emission inventory was generated for southeastern Australia from November 2019 to January 2020, which aims to fully consider the highly dynamic nature of ﬁres and small ﬁres (low FRP) that are much lower than the MODIS burned area or active ﬁre detection limit, with emphasis on dry matter burned (DMB). We found that during the study period, the ﬁre gradually moved from north to south, and the diurnal cycle of the ﬁre in the study area changed greatly. The peak time of the ﬁre gradually delayed as the ﬁre moved south. Our inventory shows that the DMB in southeast Australia during the study period was about 146 Tg, with major burned regions distributed along the Great Dividing Range, with December 2019 being the main burning period. The total DMB we calculated is 0.5–3.1 times that reported by the GFAS (Global Fire Assimilation System) and 1.5 to 4 times lower than that obtained using the traditional “Burned Area Based Method (FINN)”. We believe that the GFAS may underestimate the results by ignoring a large number of small ﬁres, and that the excessive combustion rate used in the FINN may be a source of overestimation. Therefore, we conclude that the combination of high-temporal-resolution and high-spatial-resolution satellites can improve FRE estimation and may also allow further veriﬁcation of biomass combustion estimates from different inventories, which are far better approaches for ﬁre emission estimation.


Introduction
Biomass combustion is an important global source of atmospheric emissions [1,2], which has a significant impact on air quality, climate change and human health [2][3][4][5][6]. Estimating the pollutants emitted into the atmosphere by biomass combustion is very important for studying its impact on air quality and climate change. Greenhouse gases produced by biomass combustion include carbon dioxide, carbon monoxide, methane, nitric oxide and atmospheric particulate matter, etc. [7]. Modeling studies indicate that these emissions can lead to severe regional air pollution events. Huang et al. [8] suggest This study focuses on a better understanding of the diurnal variation and spatiotemporal dynamic characteristics of fire, so as to improve the estimation results of fire combustion emissions. We use the new fire diurnal cycle information obtained from the geostationary satellite Himawari-8, to reduce the error caused by the low-time-resolution monitoring results of the medium resolution imaging spectrometer. The spatial variation of this fire diurnal cycle and its drivers is discussed in detail and is combined with the twice-daily FRP information provided by the VIIRS product. We aimed to reduce the deviation and uncertainty caused by using only one FRP data type, to take into account the burning even in the short term and to reduce the omission of active fire detection. We used land use data to mask the study area and chose the most suitable smoke emission factor in order to calculate the final total fire emissions based on the estimated value of DMB obtained by FRE.
The paper is organized as follows: Section 2 introduces the study area and data characteristics used in this study. Section 3 introduces the methods of estimating biomass combustion emissions. Section 4 shows the relevant experimental results and comparative analysis. Section 5 concludes with the major findings and problems of this study.

Study Area and Data
Section 2.1 is a brief description of the study area. Section 2.2 introduces VIIRS data, one of the data sources for estimating emissions, and makes a comparative analysis with MODIS data that is widely used at present. Himawari-8 satellite with higher time resolution for fitting Gaussian Curves is presented in Section 2.3. Land-use-type data used to distinguish different land types is described in Section 2.4. Section 2.5 describes other datasets for validation, such as GFAS and FINN data.

Study Area
Southeast Australia includes New South Wales and Victoria, with a population of about 3/5 of Australia's total population, with undulating terrain in the southeast. A vast desert is in the west and the Great Dividing Range is in the east, which runs through southeastern Australia from north to south. Forests are distributed narrowly along the Great Dividing Range, covering about 28.4 million hectares of woodland, accounting for about 18% of the national forest area (Figure 1). The main tree species are eucalypts, high in oil and flammable. Australia has an annual fire season, but since September 2019 wildfires have broken out in many parts of Australia, lasting for a long time and proving difficult to extinguish, covering a total area of more than 6 million hectares. According to the latest annual climate report of 2019 released by the Australian Meteorological Office, the average annual temperature in 2019 was 1.52 • C higher than the record set in 2013 (1.33 • C), and the average maximum temperature was 2.09 • C higher, the highest on record. At the same time, the precipitation was also the lowest on record. In 2019, the entire Australian rainfall was 62% below average, making it the driest spring on record (1967 was the lowest level). Low rainfall has exacerbated pre-existing insufficient rainfall and low soil moisture, and exacerbated meteorological and hydrological drought conditions, meaning that forest fuels are still dry (http://www.bom.gov.au/climate/current/statements/scs72.pdf (accessed June 28,2021).The research period of this article is from 7 November 2019 to 15 January 2020. During this period, the monthly average surface temperature in most parts of the southeast was above 50 • C, the monthly average wind speed reached 7 m/s and the coastal temperature was slightly lower (GES DISC (accessed March 15,2021)). High temperature, little rain and strong wind are factors that cause large-scale landscape fires to continue [38].

Polar-Orbiting VIIRS-IM Fire Radiative Power
The Visible Infrared Imaging Radiometer Suite (VIIRS) instrument onboard the polarorbiting Suomi NPP (since 2011) and NOAA-20 (since 2017) satellites has transitioned much of the capability of the AVHRR and Moderate Resolution Imaging Spectroradiometer (MODIS) instrument for environmental monitoring [40] (Table 1). VIIRS provides a total of 21 bands (plus a day and night band, DNB [41]), spanning from visible light to longwave infrared, providing better spatial resolution. The nadir pixel resolution is between 375 m and 750 m [42]. The width is 3000 km. Compared with MODIS, VIIRS has a wider bandwidth and provides gapless global coverage [43]. In addition, VIIRS pixels are generated using a pixel aggregation scheme. The samples include: the sample at the lowest point along the scanning direction (scan angle <31.72 • ), the sample in the middle area (scan angle: 31.7244.86 • ) and the sample at the edge of the scan (scan angle >44.86 • ). Therefore, the pixel size of VIIRS decreases significantly as the viewing angle increases, and the pixel size is limited to less than four times as the scanning angle increases. When the scanning angle of MODIS changes from the lowest point to the scanning edge, the pixel size along the scanning direction increased by three to nine times [44]. In addition, the VIIRS airborne bow-tie removal algorithm removes two and four over-sampled scan lines in the central gathering area and the scan edge, respectively [45], thereby reducing repeated fire detection between adjacent scans near the lowest point. MODIS does not handle it. VIIRS I4 band (3.55 to 3.93 µm) spans the wavelengths of peak spectral radiance for blackbodies emitted at temperatures between 737 and 817 K, in combination with the improved spatial resolution, the I4 band is therefore well suited for distinguishing pixels of sub-resolution combustion components [46]. Schroeder et al. [46] showed that compared with MODIS AF products, VIIRS AF products have more sensitive "small fire" active detection performance. Due to the relatively low saturation temperature of VIIRS I band, combined with improved spatial resolution, it is easy to cause frequent fire pixel saturation, so it needs to be used in combination with the M band.

Geostationary Himawari-8 Fire Radiative Power
The new-generation geostationary meteorological satellite Himawari-8 was launched on 7 October 2014. It has 16 observation belts with a spatial resolution of 0.5 or 1 km in visible and near infrared bands, a spatial resolution of 2 km in infrared bands and a temporal resolution of 10 min. It provides data covering the entire global hemisphere of East Asia and Australia, and has provided remote sensing observation data to the world since 2015, including real-time and archived formats. They are available through the official JAXA website. In order to convert a twice-daily VIIRS-IM FRP product into a daily integrated FRE, information on the fire diurnal cycle is required [37,47,48]. Apparently, compared with the observation frequency of the MODIS satellite, four times a day, the high time resolution provided by Himawari-8 can obtain the above information more accurately. However, the spatial resolution of geostationary satellites is relatively coarse, and the presence of mixed pixels increases the deviation of FRP values. Therefore, we combine VIIRS data with higher spatial resolution to reduce the deviation and uncertainty caused by using FRP data alone.

Land Cover Data
We use the Global 30 m Fine Surface Coverage Product in 2020 [39] to classify land types for the study area in southeast Australia. This dataset is based on the 2015 Global Fine Surface Coverage Product, combined with the time-series surface reflectance data of the Landsat satellite from 2019 to 2020, Sentinel-1 SAR data, DEM terrain elevation data, global thematic auxiliary dataset and prior knowledge datasets, etc., which produced the Global 30 m Fine Surface Coverage Product in 2020. Figure 1 shows the spatial distribution of the land cover in southeast Australia.

FINN and GFAS Emissions Inventory Data
We used two state-of-the-art global fire emission databases (Fire Inventory from NCAR and Global Fire Assimilation System) to compare the results of this study. FINNv2.2 infers the combustion zone based on MODIS and VIIRS active fire detection data. Since the dataset does not have the combined data of MODIS and VIIRS in 2020, we only use FINNv2.2 based on MODIS data. The FINN assumes that each active fire detection has a pixel area of one square kilometer, and the burning area is further scaled according to the exposed coverage percentage of the VCF (MODIS vegetation continuous field) product at the fire point [49]. FINNv2.2 provides estimates of global fire emissions with a resolution of 1 km per day, including trace gases and aerosols. The emissions are calculated as follows: where E i is the emission of species i, A X,t the area burned at location x and moment t, B X the biomass loading, FB the fraction of that biomass that is burned in the fire (burning ratio) and EF i the emission factor. The biomass loading is assumed to be based on a specific land cover type, using the estimation method of Hoelzemann et al. [20], and a small number of updates are described by Wiedinmyer et al. [49]. According to the theory found by Ito and Penner [19], the fraction of biomass burned (FB) is related to vegetation cover, which is further described by Wiedinmyer et al. [49,50]., and the emission factor comes from the literature investigation by Akagi et al. [51].
In contrast to the FINN, the GFAS calculates biomass combustion emissions by assimilating fire radiative power observations from the MODIS instruments onboard the Terra and Aqua satellites [52]. The GFAS forms a fire database based on active fire (AF) detection and integrates into the Copernicus Atmosphere Monitoring Service (CAMS) system for near-real-time atmospheric composition monitoring and forecasting. Using the Kalman filtering algorithm developed by Kaiser [25], the FRP data observed by MODIS four times a day are integrated into daily FRE. Then, DMB can be calculated by multiplying FRE by a conversion factor. Finally, DMB is linearly transformed into emission fluxes of 40 smoke constituents (0.1 • grid cell) using emission factors derived from field measurements [15].

Methodology
The assessment process of pollutants emitted from biomass combustion can be divided into two parts. The first part is the daily FRE generated by the Gaussian curve of the fire diurnal cycle based on the observation results of Himawari-8 and the observation results of VIIRS. The second part is to combine FRE with corresponding emission factors to generate various pollutant emissions.

Fire Diurnal Cycle and Daily FRE Generation
It has been found from previous studies [27,47] that the diurnal variation of fire has very obvious Gaussian characteristics in a day (0:00-23:00, local time).At present, most of the studies [27][28][29]53] on pollutant emission estimation are based on MODIS FRP fire diurnal cycle fitting. Vermote et al. [27] found that the parameters related to the fire diurnal cycle have a high correlation with the ratio of FRP Terra to FRP Aqua (T/A). However, for some highly dynamic fires, the low time resolution of MODIS obviously has a high error uncertainty, and the correlation between T/A and h (the peak hour of the fire diurnal cycle) is only 0.2. Although Vermote suggested that h has little influence on the final estimated DMB, it is obvious that h is also a source of error.
Here, we represent the fire diurnal cycle as a Gaussian function, parameterized using the Himawari-8 FRP diurnal cycle superimposed on a fixed baseline. Chathura et al. [54] found that the active fire detection ability of Himawari-8 under different FRP is relatively consistent with that of VIIRS. Based on the method of Andela et al. [55], the fire diurnal information obtained by Himawari-8 was mixed with VIIRS data to obtain the parameters related to the fire diurnal cycle. For a given grid cell j at instantaneous time t, the mixed FRP is calculated by: where FRP j,t is the instantaneous mixed FRP (MW) for grid cell j at time t, FRP j,VIIRS-night is the night-time (~01:00 LST) VIIRS-IM FRP value for grid cell j, FRP j,peak-Him is the peak value of the fire diurnal cycle fitted based on the hourly data of Himawari-8 FRP for grid cell j, t j is the local solar time at time step t, t j,peak is the local time at which the Himawari-8 FRP diurnal cycle peaks for grid cell j and σ is the standard deviation (dependent on fire duration) of the main peak of the Himawari-8 FRP diurnal cycle calculated by fitting a Gaussian function (using non-linear least squares) to the Himawari-8 FRP diurnal cycles. For each grid cell j, daily FRE was calculated by:

Conversion to Smoke Emissions
Wooster et al. found that the relationship between FRE and dry matter burned (DMB) was linear and significant (R 2 = 0.98), as shown by the following equation [23]: DMB j is total dry matter burned (kg) for grid cell j, FRE j represents the fire radiation energy for grid cell j and CR represents the burning conversion ratio (0.368 kg MJ −1 (±0.015 kg MJ −1 )) derived by Wooster et al. [23] from a series of outdoor experimental fires. To convert the resultant DMB into smoke emissions, we used the emission factors from Kaiser et al. [25], as is used in the GFAS V1.0 ( Table 2). DMB is then multiplied by appropriate particulate and gaseous emission factors in order to estimate the emissions from fire burning. However, the DMB is obtained according to the Himawari-8 FRP product. The spatial resolution of the land-type cover product used in this study is much higher than that of DMB, so it is necessary to calculate the mixed emission factor of each grid cell. The total area of land use type within each grid cell j was calculated by: where A j is the total area of all land-use types in a grid cell j and A i is the area of the i land-use type in grid cell j.
Mixed emission factors are calculated by multiplying the emission factors of each landuse type by the corresponding area weight in the grid cell and then cumulative calculation: where EF j,mixed is the mixed emission factor of grid cell j and EF i is the emission factor of the i type of land use in grid cell j.

Differences in Fire Monitoring
Wooster et al. [56] found that it is difficult for MODIS to detect FRP when below 10 MW, while Schroeder et al. [46] show that VIIRS is more sensitive for small fires. In order to quantify this difference, we compared the spatial distribution of FRP below 10 MW of VIIRS and MODIS Aqua during the study period. For better visualization, FRP is divided into two categories: less than 5 MW and between 5-10 MW. It can be clearly found from Figure 2 that the monitoring result FRP of MODIS Aqua below 5 MW is far worse than that of VIIRS. During the study period, a total number of 269,753 fires with FRP lower than 5 MW were detected by VIIRS, and 2624 fires with FRP lower than 5 MW were detected by MODIS Aqua. 83,900 fire points of FRP between 5-10 MW were detected by VIIRS, and 6387 fire points of FRP between 5-10 MW were detected by MODIS Aqua. Although the spatial resolution of VIIRS is higher than that of MODIS, which means that multiple VIIRS AF pixels may be represented by one MODIS AF pixel, the band feature that VIIRS is sensitive to small fires is the main reason for the higher number of fire pixels than MODIS.
Wooster [55] and Ichoku [56] also suggested that low FRP fires are the most common type of fire events at present. If this part is ignored, it may have a significant impact on the estimated results of fire emission pollutants. Obviously, the estimation of fire emission pollutants based on VIIRS has more advantages than on MODIS.

The Fire Diurnal Cycle Characteristic
The combustion process of the fire is highly dynamic; the study of the variation in characteristics of fire behavior has great value for the accurate estimation of pollutant emissions and future fire extinguishing management. As mentioned in Section 3.1, the diurnal variation of the fire has an obvious Gaussian law [27]; however, there are a few issues in its accurate estimation [28,29]. Therefore, we try to describe the variation in characteristics of fire behavior more accurately by using geostationary satellite Himawari-8.
The parameters of the fire diurnal cycle have been present in the report. The duration and peak time of fire were estimated based on the hourly Himawari-8 FRP time series. Figure 3 shows an example of a grid cell (2 km) in northeastern New South Wales from 7 November 2019 to 11 November 2019, where the last several days of the fire can be seen. The results shown in Figure 3 show that no matter the intensity of the fire, the Gaussian function can describe fire activities well. At the same time, it is also obvious that observations from MODIS sampling intervals did not observe highly dynamic changes in fire well. Based on the method of Vermote et al. [27] described in Section 3.1, we use MODIS data to estimate the duration and peak time of the fire diurnal cycle within five days of the selected grid. Compared with the results obtained from the Gaussian curve fitted with the Himawari-8 data, the method based on MODIS overestimates σ, that is, overestimates the fire duration. In this example, σ is overestimated by about 1-2 h. Moreover, the peak time is also greatly different from the results of Himawari-8, in this case, there is a difference of about 0.5-2 h across each of the five days. The time resolution of Himawari-8 is once every 10 min, so the observation of Himawari-8 is more consistent with the actual fire combustion process, which will make the error smaller in the process of integrating FRP into FRE. In addition to an observed variability in the fire diurnal cycle seen on different days, we found an interesting spatial distribution in the duration and peak time of fire ( Figure 4). We found that the peak time of fire is usually after 14:00 local time. From November 2019 to January 2020, with the southward movement of the fire, the peak time of the fire is gradually later. In November 2019, the peak time of fire was mostly 15:00. In December 2019, the peak time of fire was mostly 16:00. In January 2020, the peak time of fire mainly occurs after 18:00, and the peak time in the early morning was a higher ratio than that in November and December 2019. According to the research of Hély et al. [57], fuel moisture, wind and other ambient atmospheric conditions are far less conducive to intense fire activity by night than by day, therefore, the fire activity at night will generally be significantly reduced. The continuous burning at night in Victoria indirectly reflects the fierce fire in Victoria, which is hard to control. In addition, most σ values are between 1-3.5 h. As mentioned in Section 3.1, the value of σ is determined by the duration of the fire. Compared to December 2019, the fire in northeast New South Wales in November 2019 lasted for a shorter period of time. In December 2019, most of the fires occurred in the central-eastern and central-southern areas of New South Wales, with a σ value of 2-3 h. As the fire moves to the south, the 0-2 h σ value of January 2020 is much higher than that of November 2019 and December 2019, indicating that there are more short-duration fires in January 2020.

DMB Comparisons to the FINN and GFAS
The emission estimation results of FRP generated by VIIRS and Himawari-8 are compared with those of the GFAS and FINN. DMB is used as a common comparison indicator because it eliminates the difference in emission factors between different datasets. Overall, compared with the results of the GFAS and FINN, the estimated results of this study are of a very similar consistency with the former two in terms of time variation ( Figure 5).

DMB Comparisons to the FINN and GFAS
The emission estimation results of FRP generated by VIIRS and Himawari-8 are compared with those of the GFAS and FINN. DMB is used as a common comparison indicator because it eliminates the difference in emission factors between different datasets. Overall, compared with the results of the GFAS and FINN, the estimated results of this study are of a very similar consistency with the former two in terms of time variation ( Figure 5). The estimated DMB based on the research method in this paper has four centralized production periods, namely, from 7 November to 13 November 2019, from 3 December to 7 December 2019, from 16 December to 21 December 2019, and from 30 December 2019 to 4 January 2020 (Table 3). In the last concentrated production period, 37.97 Tg DMB was produced, accounting for about 26% of the total DMB. The first three concentrated output periods produced 25 Tg, 23 Tg and 18 Tg DMB, respectively. The DMB obtained through this paper during the study period is about 1.13 times higher than that reported by the GFAS. As described in Section 4.1, the MODIS data used by the GFAS have been confirmed to have a large number of small fire omissions, which makes the GFAS underestimate the emission results to some extent. The estimation result of pollutant dataset FINN based on burned area is higher than that of this study and the GFAS, which is about 1.4 times the estimated results of this study. Zhang et al. [58] pointed out that FB is usually overestimated when used for pollutant estimation. As described in Equation (1), FB is the key parameter in the estimation of pollutant emissions [59]. The combustion efficiency of biomass always changes during the combustion process, and is affected by the surrounding environment, such as wind speed, combustibles, humidity, etc., which is difficult to accurately measure [27], so it may be the source of the reason why the FINN result is higher. Furthermore, we focused on comparing the spatial distribution of emissions from the three datasets. The spatial resolution of our DMB is 2 km × 2 km, the spatial resolution of GFAS is 0.1 • × 0.1 • , the spatial resolution of FINN is 1 km × 1 km. The third one seems to be sparser than the former two, the reason being that the spatial distribution of FINN obtained by resampling the combustion biomass of each fire to a grid of 1 km × 1 km from the center of the fire, and so the actual range will be wider than Figure 6g-i. It can be seen that our results are consistent with the spatial distribution of GFAS and FINN ( Figure 6). The regions with higher DMB in November 2019 were mainly located in the northeast of NSW, Australia, and the major regions with higher DMB include Clarence Valley, Greater Taree, Hawkerbury, Kempsey, Lithgow, Port Macquarie-Hastings Richmond Valley, Tenterfield and Walcha, which produced about 45 Tg of DMB. GFAS is slightly overvalued compared with the results of this study in Shoalhaven. In December 2019, the regions with higher DMB were mainly located in the east, middle and southeast of New South Wales, including Bega Valley, Eurobodalla, Lithgow, Muswellbrook, Shoalhaven, Singleton and Wollondilly, which produced a total of 74 Tg of DMB. In January 2020, the areas with higher DMB were mainly located in the southeast of New South Wales and the northeast of Victoria. The major cities with higher DMB included Bega Valley, Eurobodalla, Tumbarumba, Tumut, Alpine, East Gippsland and Towong, which produced about 28 Tg of DMB. Finally, based on the emission factors and methods of each land-use type described in Section 3.2, we converted DMB into the estimated value of each pollutant emitted by fire ( Figure 7). A total of 229064 Gg CO2, 14640 Gg CO, 1807 Gg PM2.5, 78 Gg BC was emitted during the study period. It can be seen from Figure 7 that forest combustion emissions are the main source of pollution (specifically, 80.5% of CO2, 85% of CO, 90% of PM2.5 and 85.4% of BC). Figure 7 also shows that November 2019 was the month with the highest fire emissions during the study period, and the total 115966 Gg CO2 and 7437 Gg CO, 920 Gg PM2.5 and 40 Gg BC emissions in January 2020 were halved compared with December 2019. The reason for this is that the fire in January 2020 gradually extinguished after 15 January, and the fire was almost always burning in December 2019, so the amounts of pollutants emitted was high. Finally, based on the emission factors and methods of each land-use type described in Section 3.2, we converted DMB into the estimated value of each pollutant emitted by fire ( Figure 7). A total of 229064 Gg CO 2 , 14640 Gg CO, 1807 Gg PM 2.5 , 78 Gg BC was emitted during the study period. It can be seen from Figure 7 that forest combustion emissions are the main source of pollution (specifically, 80.5% of CO 2 , 85% of CO, 90% of PM 2.5 and 85.4% of BC). Figure 7 also shows that November 2019 was the month with the highest fire emissions during the study period, and the total 115966 Gg CO 2 and 7437 Gg CO, 920 Gg PM 2.5 and 40 Gg BC emissions in January 2020 were halved compared with December 2019. The reason for this is that the fire in January 2020 gradually extinguished after 15 January, and the fire was almost always burning in December 2019, so the amounts of pollutants emitted was high.

Conclusions
In this paper, we used a new approach to improve the estimation of fire pollutan emissions by combining VIIRS and Himawari-8 FRP products and produced a more ac curate fire emission inventory for the period of November 2019 to 15 January 2020 in southeast Australia. We have adopted the medium resolution satellite VIIRS with more superior performance and improved monitoring capacity for small fires. Compared with the results obtained by MODIS, we found that the results based on Himawari-8 are more accurate. We also found that there are significant spatial changes in the fire diurnal cycle The peak time of fire usually occurs after 15:00, but as the fire moves southward, the peak time is delayed. For example, for the fire in northeastern Victoria in January 2020, the peak time is usually after 19:00. σ is mostly 1-3.5 h, and the fire duration varies from region to region. The fire diurnal cycle characteristics determined by using the Gaussian law with VIIRS and Himawari-8 datasets can constrain the observations of fire behavior more con sistent with the actual burning process, which makes the calculation error of FRE smaller The detailed comparison with the existing inventory shows that in terms of spatia distribution, our emission inventory is in good agreement with the existing inventory, the main burning cities distributed along the Great Dividing Range. In terms of time change of emission results, our emission inventory is 0.5-3.1 times that reported by the GFAS The MODIS data used by the GFAS have been confirmed to have a large number of smal fire omissions, which makes the GFAS underestimate the emission results to some extent Zhang et al. discussed the shortcomings of the GFAS in more detail [58]. The total emis sions shown in the emission inventory of this study are 1.5-4 times lower than those ob tained using the burned-area-based approach (FINN). Studies have shown that the FINN often overestimates the burn rate of a fire or uses an outdated or inappropriate burn rate which may be the source of the reason why the FINN result is higher [58]. This shows the need for great care in estimating the inventory based on the burned area, and the uncer tainty of this method is relatively large. The satellite remote sensing method based on FRP

Conclusions
In this paper, we used a new approach to improve the estimation of fire pollutant emissions by combining VIIRS and Himawari-8 FRP products and produced a more accurate fire emission inventory for the period of November 2019 to 15 January 2020 in southeast Australia. We have adopted the medium resolution satellite VIIRS with more superior performance and improved monitoring capacity for small fires. Compared with the results obtained by MODIS, we found that the results based on Himawari-8 are more accurate. We also found that there are significant spatial changes in the fire diurnal cycle. The peak time of fire usually occurs after 15:00, but as the fire moves southward, the peak time is delayed. For example, for the fire in northeastern Victoria in January 2020, the peak time is usually after 19:00. σ is mostly 1-3.5 h, and the fire duration varies from region to region. The fire diurnal cycle characteristics determined by using the Gaussian law with VIIRS and Himawari-8 datasets can constrain the observations of fire behavior more consistent with the actual burning process, which makes the calculation error of FRE smaller.
The detailed comparison with the existing inventory shows that in terms of spatial distribution, our emission inventory is in good agreement with the existing inventory, the main burning cities distributed along the Great Dividing Range. In terms of time change of emission results, our emission inventory is 0.5-3.1 times that reported by the GFAS. The MODIS data used by the GFAS have been confirmed to have a large number of small fire omissions, which makes the GFAS underestimate the emission results to some extent. Zhang et al. discussed the shortcomings of the GFAS in more detail [58]. The total emissions shown in the emission inventory of this study are 1.5-4 times lower than those obtained using the burned-area-based approach (FINN). Studies have shown that the FINN often overestimates the burn rate of a fire or uses an outdated or inappropriate burn rate, which may be the source of the reason why the FINN result is higher [58]. This shows the need for great care in estimating the inventory based on the burned area, and the uncertainty of this method is relatively large. The satellite remote sensing method based on FRP can clearly express the intensity of fire, and combined with the remote sensing observation of a geostationary satellite, can observe the dynamic nature of fire more accurately, which is a recommended fire emission estimation method.
Author Contributions: F.Y. and Y.Z. designed the experiment. J.L. provided data and supported the data analysis. Y.Z. and J.L. carried out experiments and verification. Y.Z. drafted the paper. F.Y., H.J. and X.F. helped with editing and finalizing the paper. All authors have read and agreed to the published version of the manuscript.