1. Introduction
Forest fires are considered common and even positive phenomena; fires free forests from old trees and makes it possible for young trees to develop and absorb CO2. The greenhouse gases released by fires, primarily carbon dioxide, are therefore part of the natural global carbon cycle. The above, however, is only true in a stable climate. Under the influence of human activity, the climate is changing, for example, temperatures are rising, which leads to an increase in the intensity and number of droughts that facilitate fires [
1]. The contribution of wildfires to the carbon balance may be increasing [
2,
3]. Carbon dioxide (CO2) and methane (CH4) add to anthropogenic emission. Rising CO
2+CH
4 concentrations increase temperature trend further, as a positive feedback [
4].
In recent years, fires in Russian forests have become known as catastrophic: they span over vast areas, creating thick smoke that spreads to regions adjacent to the fires [
5]. Moreover, longer living gaseous emissions from Siberian fires fill almost completely the extratropical latitudes of the Northern Hemisphere [
6,
7]. The 2021 fires in Yakutia (Sakha) and the Krasnoyarsk region of Russia attracted particular attention [
8,
9].
By a coal fire we mean the uncontrolled burning or smoldering of a coal seam [
10]. Coal fires are nothing new. Since the 60s of the last century, coal seams have been constantly burning near the city of Centralia in eastern Pennsylvania, USA [
11]. Attempts to extinguish this fire were unsuccessful, and residents were forced to leave the city. The most famous fire in the World is the so-called "Burning Mountain" in south-eastern Australia [
12]. The underground fire there is believed to be smoldering at a depth of about 30 m and has been going on for the last 6,000 years. Coal fires are reported from China, India, USA, Australia, Indonesia, Venezuela, South Africa, Russia, Croatia and other countries ( [
10] and references therein). The fires of this type can start from spontaneous combustion or from human exposure. The reaction of carbon in coal with oxygen in the air is an exothermic process that releases heat as carbon and oxygen react to form CO2, CO and a variety of other products [
12]. Coal fires have been known to burn for hundreds of years and date back to the Pleistocene [
13].
Remote sensing research on coal fires mainly focuses on detecting thermal anomalies associated with coal combustion, detecting land cover changes under mining conditions, detecting land surface subsidence due to volume loss underground, and developing methods to identify risk areas for coal fires [
11]. In addition, in-depth thermal analysis such as subpixel thermal mapping, emissivity-based geological mapping, and coal fire quantification have been conducted [
14].
Remote measurements of CO2 from satellites require extremely high precision. For CO these requirements are dramatically softer. Let’s compare satellite measurements of CO and CO2 in wildfires. The percentage change in concentration caused by the fire is mostly important. For example, the catastrophic fires in Indonesia in 2015 resulted in an increase in average whole-atmosphere CO2 concentration of 2 ppm (2000 ppb) from a background concentration (outside the fire zone) of ~400 ppm [
15]. Thus, the impact of fire is only 0.5%. This value is close to accuracy of the most advanced satellite CO2 instrument OCO-2 (0.3% or ~1 ppm) [
15]. Fires of similar intensity in Yakutia in 2021, according to the AIRS satellite sounder (see below), led to an increase in the monthly Xco to 160–170 ppb. For the same fire the more accurate TROPOMI ( see Section 3.5 for abbreviations) sounder detected an increase up to 200 ppb with a background value for CO of just ~80 ppb. Therefore, fire signal for CO amounts to 100-200% whilst CO2 increase is only 0.5%.
There is a large number of studies of CO emissions from wildfires based on measurements of CO from satellites based on data from TIR instruments (4.6 μm of wavelength) [
16,
17], The weak point of TIR technique is a low sensitivity to the surface layers of the atmosphere. Significant progress in this area was achieved in 2017 with the launch into orbit of the Sentinel-5p satellite equipped with the TROPOMI SWIR sounder. A set of publications devoted to analysis of TROPOMI data for wildfires appeared on last four years [
18,
19,
20,
21,
22]. Combination of TROPOMI CO measurements with an approach that uses satellite-observed fire radiative power (FRP) as proxy for fuel availability was examined by [
23].
The main goals of this article is to demonstrate the capabilities of TROPOMI for studying pollution in the low troposphere and relationships between coal and forest fires. Also estimates for CO fire emission from the extra-tropical latitudes of the NH, updated to present time, are reported.
2. Materials and Methods
TROPOMI is an imaging spectrometer on the ESA Sentinel-5 Precursor platform in orbit at 824 km altitude. Its 2600 km swath provides global daily coverage at very high spatial resolution, of 7.2 × 7.2 km
2 at nadir [
24] I It measures radiation in the ultraviolet, visible and shortwave infrared spectra. Total CO column (TC) values, expressed in mol/m
2 and converted into Xco, column-averaged volumetric mixing ratio, were obtained from measurements of reflected solar SWIR radiation in the first overtone of the CO fundamental band in the spectral range of 2.3 μm [
25]. Over land Xco were retrieved from solar radiance for clear sky and cloudy conditions. Measurements of CO over water are possible only when clouds are in the field of view [
25]. Otherwise, due to the low reflectivity of open water in the SWIR range, the spectra of reflected light are too noisy. A special procedure is applied to estimate the fractional CO column below the cloud top [
25]. A priori profiles are generated using the global chemical transport model TM5 [
26]. The profiles vary by place, month and year.
AIRS is a diffraction grating spectrometer launched into a sun-synchronous polar orbit in May 2002 on board the Aqua satellite [
27]. The instrument has a mirror scanning ±48.3° from nadir across the orbit. It provides nearly complete global daily coverage. The spectral resolution is 1.79 cm
-1 in the fundamental absorption band of CO in the wavelength region of 4.6 μm. The instrument has a spatial resolution of 13.5 km at nadir. Currently (December 2023) AIRS is still operational. Monthly and daily average "AIRS-only" VMR profiles (version 7, level 3 [
28]) for the period October 2002 to November 2023 for ascending orbits are available online in a 1° × 1° latitude-longitude grid at
https://disc.gsfc.nasa.gov/datasets/ (accessed December 10, 2023), AIRS3STM_7.0 (monthly), and AIRS3STD_7.0. (daily). Xco were calculated as weighted averages of CO VMR values as a function of pressure for 24 layers of different thicknesses, used as weighting factors.
This article is based on two remote sensing experimental data sets. A comprehensive validation of the TROPOMI/Centinnel-5p TC-averaged mixing ratio of carbon monoxide and methane (Xco and Xch4) has been published [
29]. 28 ground-based TCCON and 24 NDACC stations, included in that study, were located mostly in remote areas and characterize background concentrations of greenhouse gases. The stations were equipped with Fourier Transform interferometers of high spectral resolution. A system of mirrors tracked the Sun. The data were available only for clear sky conditions. Excellent agreement between TC-averaged mixing ratios of gases measured from the satellite and from the ground was documented.
AIRS Xco were validated by comparison with data from 9 TCCON sites [
7]. These comparisons reveal a significant scatter of points and underestimation of AIRS’ Xco.
Appendix A). The mean slope of linear regression line was found to be 0.73, and the standard deviation for the slopes was estimated as ± 19%. The measured HNH anomalies of CO burden were corrected accordingly. The low Xco for AIRS are explained by reduced sensitivity of the retrieval algorithm to lower atmospheric altitudes. The AIRS sensitivity to the middle layers, namely 3-5 km of altitude, was verified by comparisons with aircraft CO in situ data in this range of altitude from 32 sites in USA and Canada (6905 profiles) [
30] (
Figure A2). A significant scatter for measurement days > 10 ppb was found. On average the binned values of AIRS Xco differed from the in situ aircraft data by less than 10 ppb within an Xco range of 90 to 140 ppb.
3. Results
3.1. Two Decades of Wildfire Emissions
Monthly and yearly CO emissions from wildfires in the HNH were calculated by [
7] using a two-box (0 - 30 N and 30 - 90 N) mass balance model for the period between January 2003 and December 2021 (for model description see [
16,
17]). In 2021 a record high annual emission of 149 Tg CO yr
-1 was found. This calculation was updated here using the same algorithm until October 2023 (
Figure 1). The Global Fire Emission Database, version 4c (GFED4c, described by [
31]) specified even a higher magnitude for 2021: 192 Tg CO yr
-1 [
7]. The plot of
Figure 1 represents cumulative CO mass emitted in the HNH from the beginning of the given year. 2003, 2012, 2021 and 2023 were the years of "mega-fires".
Figure 1 clearly demonstrates also a range of year-to-year variability of CO annual emissions for the "normal" fires. The time series of CO annual emissions is shown in
Appendix B (
Figure A3 ). A (5.1 ± 2.8) % yr
-1 up-going trend for "normal" years was observed
The design of our box model does not allow to derive spatial distribution of emission inside the selected box. On the contrary, GFED4c approach provides global maps of CO fire emissions since 1997 with time resolution of one month. We used these maps to illustrate temporal variations of CO discharge from wildfires for various areas in the Russian Federation (
Figure 2). Emission from Central Siberian forest/tundra domain (for boundaries see maps of
Figure 3d and
Figure A4) was compared with that from the rest of boreal Russia after 1997. The most striking feature of this plot is growing CO emission at the Central Siberia, especially after 2015 and up to maximum in 2021. In that year a contribution of the Central Siberia reached ~ 90% of total Russian pyrogenic CO output. Before 2015 emissions from that area for the most of previous period (12 years of total 19 ) were insignificant. Anthropogenic CO production in Russian Federation is in the range between 25 and 30 Tg CO yr
-1 [
32], i.e. it is comparable with those from normal fires or lower than CO discharge from mega-fires.
3.2. Mega-Fire in Central Siberia in 2021 as Seen Remotely by TROPOMI and AIRS.
In this and the next sections CO measurements provided by two satellite-borne sounders for fire regions will be analyzed. Shortcomings of TIR and SWIR techniques will be discussed below in the
Section 4.1. Geometry and spectroscopy of these two instruments are different. The TROPOMI registers Solar SWIR radiation penetrating the total atmospheric column and reflected by the surface. CO spectral features allow to calculate total number of CO molecules along this track. Dividing this number by the number of air molecules (excluding water vapor molecules) on the same way gives us total-column-average Volume Mixing Ratio (VMR), measured in parts per billion (ppb), usually denoted as Xco. The optical path can be shorter or longer depending on cloudiness and aerosol multiple scattering. The problem of accurate determination of the optical path is a central one for any SWIR or visible solar technique.
The AIRS grating spectrometer measures spectrally resolved Thermal IR (TIR) outgoing terrestrial radiation. The initial source of this radiation is the Earth surface, on the way to the satellite the radiation is absorbed and re-emitted by CO. The shape of absorption/emission spectral lines depends on air pressure and air temperature. They are different for different altitudes. A vertical profile of CO VMR for any location on the Globe and any time may be retrieved. A special complicated mathematical algorithm for solving this problem (called inverse problem) is applied. A significant drawback of the TIR approach is a dependence of the sensitivity on the atmospheric temperature profile: as a result, the sensitivity to CO in the lowest atmospheric layer diminishes. On the contrary, the sensitivity of SWIR techniques to all altitudes is almost the same for clear sky. This difference in sensitivity opens a door to separate TC-averaged (Xco) and low tropospheric (ideally, boundary layer) VMR. The latter is highly important for investigation of wildfires, namely to distinguish between flaring and smoldering types of burning. The shortcomings of the specific datasets, however, will be discussed below.
Emission of CO from the Siberian 2021 fire maximized in July and August. Monthly August Xco measured by TROPOMI and AIRS (
Figure 3 a and b, respectively) have similarities and differences. Both sounders registered increased Xco over Central Siberia domain bounded by the rivers Yenisey (on the West) and Lena (on the East). TROPOMI’s Xco for this area exceeds that measured by AIRS. Both instruments measured similar levels of Xco over the Pacific ocean and areas beyond East and/or Central Siberia. Pyrogenic CO from Siberia was transported northward and polluted the Arctic ocean area and the Bering sea. TROPOMI clearly "see" the polluted area in China.
Figure 3c displays the differential between these two data sets (TROPOMI minus AIRS). The ΔXco varies between ~0 and 100 ppb. If ΔXco = 0, then it means lowest pollution cases. In other words, real CO profile does not differ from a priori one. The a priori is obtained from a cluster of available aircraft and surface in situ measurements for unperturbed atmosphere. Therefore, ΔXco can be treated as a measure for deviation of lower tropospheric VMR from a priori. Mostly important is a character of emission mechanism. Urban emissions are connected with low-temperature processes. Flaring wildfires, on the contrary, generate very hot burning products that are mixed throughout the atmospheric column effectively and easily reach the altitudes of several kilometers. Intermediate examples are smoldering peat and coal fires. Using TIR and SWIR remote sensing is promising for elucidation of the fire’s nature. Shortcomings of the current versions of retrieval techniques will be discussed in the
Section 4.
The map of
Figure 3d is a zoom-in of
Figure 3a. It illustrates a general similarity in the shape of enhanced Xco area and locations of two coal basins, Tunguska (TCB) and Lena (LCB). Maximal monthly average Xco (up to 240 ppb) were recorded by TROPOMI to the North of Yakutsk (the capital of the Sakha Republic) near the small settlement Sangar. An underground coal fire there began in May, 2000, it has not been extinguished by 2011 [
33]. Citizens of Sangar suffer pollution from this fire year-round still now [
34].
3.3. Mega-Fire in Australia in 2019/2020.
The summer season of 2019/2020 was characterized by record severe wildfires in Australia [
35].
Figure 4, similarly to
Figure 3, displays Xco for January, 2020, provided by two sensors, together with a differential between them. The color scale for
Figure 4, a and b, is from 40 ppb to 160 ppb, slightly different from that for
Figure 3 (CO background concentration in the Southern hemisphere is lower than in Northern). The color scale for
Figure 4c is between -10 and 100 ppb, i.e., the same as for Siberia. A plume of fire products, including CO, stretched across the Pacific Ocean until South America’s shore. Origin of the plume was in South-Eastern Australia near Sydney city. These three maps (
Figure 4, a, b, and c) evidence in favor of CO eastward transfer through the lowermost atmospheric layers.
Figure 4d is a map of individual TROPOMI retrievals (pixels) for December 30, 2019. The burning area coincided with the shape of the Sydney coal Basin.
3.4. Abbreviations
TROPOMI: . TROPOspheric Monitoring Instrument
AIRS Atmospheric InfraRed Sounder
TCCON: Total Carbon Column Observing Network
NDACC: Network for the Detection of Atmospheric Composition Change
HNH: High Northern Hemispheres
TCB: Tunguska Coal Basin
LCB: Lena Coal Basin
NOAA: National Oceanic and Atmospheric Administration
SWIR: Short-Wave InfraRed
TIR: Thermal InfraRed
OCO-2: Orbiting Carbon Observatory-2
ESA: European Space Agency
VMR: mixing ratio by volume