=Paper=
{{Paper
|id=Vol-3006/41_regular_paper
|storemode=property
|title=Satellite monitoring of the environment in the area of the Iskitim coal mines in 2013-2020
|pdfUrl=https://ceur-ws.org/Vol-3006/41_regular_paper.pdf
|volume=Vol-3006
|authors=Natalia V. Rodionova
}}
==Satellite monitoring of the environment in the area of the Iskitim coal mines in 2013-2020 ==
Satellite monitoring of the environment in the area of the Iskitim coal mines in 2013–2020 Natalia V. Rodionova1 1 Kotel’nikov FIRE RAS, Fryazino, Russia Abstract The paper considers the use of multispectral data from the Landsat-8, Sentinel-2, Aqua and Terra satellites for monitoring pollution in the areas of open-pit coal mines in the Iskitim district of the Novosibirsk region for the period 2013–2020. The change in the values of the reflection coefficient (RC) from the surface and water bodies, the snow index NDSI during the snowmelt period, the change of NDVI in the summer, in the area of Kolyvansky and Vostochny coal mines and in the area of the Linevo village are considered. The dynamics of the aerosol optical thickness (AOT) changes, CO and CH4 concentrations in the atmosphere of the Iskitim district using the Giovanni data analysis and visualization system are shown. Keywords Remote sensing, multispectral images, coal mine, surface and air pollution. 1. Introduction Open-pit coal mining, its transportation have a negative impact on the environment, polluting the surface and air with emissions of solid substances in the form of coal and inorganic dust, coal ash and soot. In addition, more active coal mining is the cause of man-made earthquakes due to the movement of rocks [1]. Over the past 10 years, the frequency of earthquakes in the Novosibirsk region has increased-this is directly related to active coal mining. One of the last earthquakes occurred in September 2020 in the vicinity of the Kolyvansky coal mine (CM) near the village of Ust-Chem. Satellite images make it possible to assess the area and degree of pollution of the territories adjacent to coal mining, mainly during the snowmelt period. In this paper, satellite monitoring of snow cover pollution in the Iskitimsky district of the Novosibirsk region near Kolyvansky and Vostochny CM and in the vicinity of the Novosibirsk Electrode Plant in the village of Linevo was carried out using multispectral data from the Landsat-8 (L8) and Sentinel-2 (S2) satellites to determine the dynamics of changes in the values of the reflection coefficient from the surface, the NDSI snow index during snowmelt and the NDVI vegetation index in summer for the period 2013–2020. The work also makes an assessment of atmospheric pollution based on the use of data from the Terra/Aqua satellites. The estimated parameters include the aerosol optical thickness, CO and CH4 concentrations in the atmosphere of Iskitim district using the Giovanni data analysis and visualization system [2]. SDM-2021: All-Russian conference, August 24–27, 2021, Novosibirsk, Russia " rnv1948123@yandex.ru (N. V. Rodionova) © 2021 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings http://ceur-ws.org ISSN 1613-0073 CEUR Workshop Proceedings (CEUR-WS.org) 353 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 2. Materials and methods 2.1. Input data and study area The work uses Landsat-8 OLI C2 L2 images with radiometric calibration and atmospheric correction, which can be accessed for free, for example, through the EarthExplorer service [3]. The channels used (B2–B7) of each image were converted from the brightness values DN to the reflectivity values of the underlying surface. The spatial resolution for multispectral channels is 30 m. Cloudless images of early April (snowmelt) and summer (July–August) 2013–2020 were used. The work also uses Sentinel-2 multispectral data with high temporal, spatial, and spectral resolution. The multispectral camera has 13 channels with different spatial resolutions from 10 to 60 meters. We used data from the S2 survey system with the L2A processing level in the form of an albedo at the lower boundary of the atmosphere (with atmospheric correction). S2 data processing was carried out by the program SNAP [4]. S2 images were used in cases when Landsat 8 images were missing or covered in clouds. For the analysis of atmospheric pollution, we used the NASA Goddard online visualization and analysis tool Giovanni [2] which allows usage, comparison, and analysis voluminous remote sensing (RS) data (NASA MODIS (Terra and Aqua), ESA MERIS (ENVISAT)) without downloading. The study area is located in the Iskitimsky district of the Novosibirsk region (Eastern Siberia). Iskitim district is a territory with a high technogenic load. There are several coal basins with open-pit coal mining, a cement plant (Iskitim), the Novosibirsk Electrode Plant (Linevo), the raw material base for which is anthracite deposits. All this, as well as all the processes associated with the operation of coal mines, transportation of coal, lead to negative consequences for the natural environment. Figure 1 shows the study area in the image S2, channel B4, the date of the survey 02.04.2020, with a characteristic view of snow cover contamination in the area of Iskitim (cement plant) (1), Figure 1: Image S2 of Iskitim region, channel B4, 02.04.2020. 354 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 the village of Linevo (Novosibirsk Electrode Plant) (2), the areas of the Kolyvansky and Vostochny coal mines (3), Urgunsky, Gorlovsky coal mines and highways along which coal is transported. The upper left part of the image shows a fragment of the Ob reservoir. 2.2. Assessment of anthropogenic impact on surface pollution It is shown [5] that the minimum time interval for the study of man-made impacts from mining on the environment is a period of 5 years. In this paper, this period is equal to 8 years from 2013 to 2020. 2.2.1. Snow cover during snowmelt To assess the level of surface pollution in the winter-spring period, the assessment of snow cover pollution during the snowmelt period is used. Snow cover, having a cumulative effect, allows to get the real total amount of precipitation of pollutants. The evaluation criteria are the reflective characteristics of the snow cover in various spectral channels and snow indices, one of which is the snow index NDSI (normalized difference snow index) [6], defined as GREEN − SWIR1 NDSI = . GREEN + SWIR1 The threshold for identifying snow is the value NDSI ≥ 0.4 [6]. Let’s compare the reflection coefficient (RC) from the snow cover and NDSI near Linevo (Figure 1, (2)) and Kolyvansky coal mine (Figure 1, (3)) in 2013 and 2020. The number of coal mines in the Kolyvansky district increased from 2013 to 2020 (Figure 2, L8 images). Figure 3 shows graphs of changes in the reflection coefficient and NDSI along the profiles in the Kolyvansky and Vostochny CM region during the snowmelt for 09.04.2013 (L8) (top row) and 02.04.2020 (S2) (bottom row). Profile consists of 14 points. The distance in km from the “zero point” is plotted along the abscissa axis (the value of 0 km corresponds to the point with coordinates 54∘ 39’47"N, 83∘ 37’30"E, near the road to CM. The NDSI value for all points of the 14.07.2013 11.06.2016 02.08.2020 Figure 2: Dynamics of changes of the Kolyvansky and Vostochny coal mine territory from 2013 to 2020 according to L8 data. 355 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 Figure 3: Graphs of the RC and NDSI change for the profile in the area of Kolyvansky and Vostochny coal mines for April 2013 and 2020. Figure 4: RC and NDSI change along the profile for 09.04.2013 (L8) and 02.04.2020 (S2). 2013 profile is greater than 0.5, i.e. snow is identified at these points. For 2020, snow has not been identified for profile points up to 2.5 km in the direction of Ust-Chem (from 0 point) and about 1.2 km (from 0 point) in the western direction. 356 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 Snow pollution in the area of Linevo is associated with the work of the Novosibirsk electrode plant (NEP). Dust removal from high-rise pipes of the plant, coal processing, quarry operations, and technological route in the vicinity of the plant are the main sources of pollution in this area. Figure 4 shows the graphs of RC and NDSI changes for April 9, 2013 (L8) (top row) and April 2, 2020 (S2) (bottom row) (snowmelt) along the 8-point profile passing through Linevo. The distance in km from the “zero kilometer” is plotted along the abscissa axis (the value of 0 km corresponds to the point with coordinates 54∘ 31’15"N, 83∘ 24’05"E, near the NEP). The snow-free zone in 2020 compared to 2013 has practically not changed in the direction to the west from the 0 point, at the same time it has expanded in the eastern direction by 4–5 km. 2.2.2. Vegetation To assess the state of vegetation, integral characteristics are used, expressed in the form of vegetation indices (VI). The basis for the choice of VI is the features of the reflectivity curves of vegetation, namely, in the red region of the spectrum (RED), the maximum absorption of solar radiation by chlorophyll is noted, in the near infrared (NIR), the maximum reflection of energy by the cellular structure of the leaf. As a result, the normalized difference index NIR − RED NDVI = NIR + RED is the most common method for quantifying the phytomass density. Technogenic factors associated with the open-pit coal mining method have a negative impact on the vegetation cover. Pollution of the atmosphere, soil, water environment, etc. are factors of impaired growth and development of vegetation. In this paper, the values of RC and NDVI for the profile in the village of Linevo are compared according to L8 data for the summer time period for the dates 30.07.2013 and 02.08.2020 (Figure 5, the profile is the same as in Figure 4). Note in Figure 5 the local minimum of the values of the RC at the point of “zero” km (near the NEP). The NDVI values for all points of the profile, except for the 9.2 km point, for 2013 exceed the NDVI values for 2020. Figure 6 shows RGB images of the studied territory for the summer period, demonstrating the dynamics of changes in the values of NDVI = (B5 − B4)/(B5 + B4) (L8 channels) for Figure 5: RC and NDVI changes for 30.07.2013 and 02.08.2020. 357 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 30.07.2013 04.08.2018 02.08.2020 Figure 6: RGB images of NDVI value changes. dates 30.07.2013, 04.08.2018 and 02.08.2020 with the following encoding of NDVI values: red — NDVI < 0.4, green — NDVI > 0.7, blue — NDVI ≥ 0.4 and NDVI ≤ 0.7. There is a tendency from 2013 to 2020 to increase the territories with NDVI < 0.7. 2.2.3. Water bodies In Figure 7, the numbers in the image Landsat-8 (B5, 2.8.2020) indicate the monitoring water objects. These are: 1 — the Ob reservoir; 2 — the Berd River near the bridge; 3 — the Berd River; 4 — the Linevo Lake; 5 — the Koinikha Lake. The spectral curves for the wavelengths of Figure 7: Location of water bodies monitoring points. 358 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 L8 B2–B7 are shown in Figure 8 for four dates 30.07.2013, 08.07.2014, 31.07.2019 and 02.08.2020. The spectral curves for 2014 and 2019 are presented here to demonstrate a significant difference between the spectral curves for 2013 and for 2014–2020 in terms of the RC value in the B5 channel (NIR). The data in 2013 are characterized by a significant increase in RC in three and more times in NIR channel for the considered water bodies except for the Ob reservoir. In the work [7], it is shown that the two spectral zones NIR and SWIR are very sensitive to the presence of impurities in a water medium. In closed reservoirs (Lake Linevo), such impurities settle to the bottom, while for water bodies with an active current of water the particles are mixed, which causes a jump in the amount of RC. Moreover, the presence of suspended particles in the water medium increases the RC in NIR spectrum, and organic substances in SWIR1. For 2013, there are two peaks in the RC values: in the green region of the B3 spectrum and in NIR (B5). The first maximum is present on all the charts presented. The increase RC in B3 channel is associated with an increase in the volume of phytoplankton, while there is a decrease in the red zone of the B4 spectrum, which is associated with the absorption of chlorophyll by plankton and so a decrease RC in the B4 channel. For all graphs, the smallest changes in channels B3 and B4 were found for object 4 — Lake Linevo (except 2013). The second maximum (channel B5) is associated with the presence of suspended impurities in the water medium. In 2013, the value of the RC for B5 is in three and more times higher than the RC for 2014–2020. The exception is the value of the RC for the Ob reservoir, for which the values of the RC in the channel B5 are the smallest relative to the other considered water bodies. The second maximum, pronounced in 2013, gradually comes to naught since 2014, and the RC values gradually decrease during the transition from the B5 channel to the B7 channel. Figure 8: Spectral curves for water bodies: 1 — Ob reservoir; 2 — the Berd River near the bridge; 3 — the Berd River; 4 — the Linevo Lake; 5 — the Koinikha Lake. 359 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 Based on the values of the RC for the B5 channel, the Ob reservoir is the cleanest. The present contamination is suspended particles, inorganic impurities. 2.3. Assessment of air pollution Iskitim is one of the cities with a high level of air pollution (AP). The substances that determine the level of AP in Iskitim are benz(a)pyrene, suspended substances (dust), CO — carbon monoxide, NO2 /NO — dioxide/nitrogen oxide, SO2 — sulfur dioxide, hydrogen sulfide and soot. Iskitim is located in a low-lying area, which contributes to the accumulation of harmful substances in the air above the city. This paper examines the dynamics of changes in the aerosol optical thickness (AOT), CO and CH4 values in the atmosphere of Iskitim district in 2013–2020 using the Giovanni data analysis and visualization system [2]. 2.3.1. Seasonal and interannual variations of AOT Among the radiation-significant components of the atmosphere, along with greenhouse gases and clouds, an important role plays atmospheric aerosol, the main optical characteristic of which is the AOT of the atmosphere. Aerosol particles have a significant effect in the temperature balance of the atmosphere, as they are able to absorb, disperse, reflect solar radiation and initiate cloud formation processes. The main natural factors that affect the long-term variability of AOT are volcanic eruptions [8] and forest fires [9]. In open-pit coal mining, the main pollutants are solid emissions: inorganic dust, coal dust, coal ash, and soot. For the analysis of spatiotemporal variations of AOT, data from satellite measurements, in particular the MODIS radiometer of the Terra/Aqua satellites, are increasingly used. To estimate atmospheric aerosol variations, we used data from the level 3 (average daily) collection 6.1 of the MODIS radiometer (Terra satellite) — the product Terra MOD08_D3 v6. 1 (Deep Blue, Land only) at a wavelength of 550 nm. The AOT data obtained with MODIS is a two-dimensional field (image) integrated in height. The dimensions of the individual elements (pixels) of the map are 1 × 1∘ (latitude × longitude). The average monthly AOT values for the period April–October 2013–2020 were calculated, averaged for a 1 × 1∘ domain with a center coordinates 54.5∘ N, 83.5∘ E, using the Giovanni data analysis and visualization system [2]. Figure 9 shows the average long-term AOT values for the months of April–October 2013–2020 and SD AOT. We note the seasonal course of atmospheric turbidity with a maximum in the spring-summer period (April and July) and a minimum in the autumn (there are no data for the winter period). Similar trends have been observed for the Baikal region [10]. An increased AOT values in April is observed for the entire study period. In July the increased AOT values most pronounced for 2019, and slightly weaker for 2016. We should note also the strong spread of AOT values in July 2019 (see SD AOT values). One of the main reasons for the April AOT maximum in Iskitim region is the release by snowmelt of aerosol accumulated during the winter period in the snow cover. The local minimum in May can be associated with the formation of fresh vegetation cover and an increase in rain level that washes out the aerosol. The AOT maximum in July is most likely associated with an increase in temperature, a decrease in humidity and drought. 360 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 Figure 9: Average AOT and SD AOT values for April–October 2013–2020 for the study area according to MODIS (Terra) data. Figure 10: Interannual variations of AOT for the period 2013–2020. Figure 10 shows a graph of the interannual variations of AOT for the period 2013–2020. No annual trends were found. Large variations of the interannual AOT values have been detected. Only for 2018, the average AOT value is close to the average annual global value of 0.16, obtained from satellite data MODIS, TOMS, AVHRR [11]. For the rest of the dates for the Iskitim region, the average long-term value of AOT is significantly higher than average annual global value and the value for the Baikal region (data from MODIS and ground-based solar photometer measurements [10]). 2.3.2. Seasonal and interannual variations of CO The main sources of CO are the burning of fossil fuels, biomass, and the atmospheric oxidation of methane and other hydrocarbons. Carbon monoxide is a colorless and odorless gas. CO is formed by the combustion of organic fuels (wood, coal, paper, oils, gasoline, etc.). The dominant part (∼90%) of atmospheric CO is formed as a result of natural processes (volcanic and swamp gases, forest and steppe fires, the vital activity of terrestrial and oceanic flora and fauna). At the same time, about 70% of the CO in the atmosphere is the product of methane oxidation in the troposphere with the participation of the hydroxyl radical (OH). On a global scale, atmospheric methane (which is oxidized by OH) is the main source of CO in the atmosphere. In large cities and industrial regions, the preponderance may be on the side of CO, directly emitted from 361 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 anthropogenic sources. Hundreds of millions of tons of CO enter the atmosphere annually as a result of human activities: motor transport, rail and sea transport; malfunction of gas pipelines and gas equipment; metallurgy, chemical, coal mining, tobacco, bread production, light copying; waste processing; fuel combustion in everyday life. In industry, CO is produced by incomplete oxidation of natural gas or by gasification of coal and coke. CO enters the atmospheric air during any type of combustion. In cities CO gets into the air , mainly in the composition of exhaust gases from internal combustion engines. CO has a long life time — about a month [13]. Daily data of the AIRS/Aqua instrument product AIRS3STM v006, averaged over an area of 54–55∘ N, 83–84∘ E were used to study variations in the atmospheric CO content. Figure 11, a shows the average monthly values of the CO concentration and Figure 11, b shows the SD values for the period 2013–2020 (common for all the 8 years from 2013 to 2020 are 1) the highest concentrations of CO occur in the winter months, as well as in the spring months of March and April. The summer months are characterized by the lowest concentrations of CO. The CO concentration values in the autumn months are higher than in the summer months, but lower than in the winter months. Figure 12 shows the average annual CO concentrations for 2014–2020. The averaging was carried out for 11 months (the month of July was not included due to the lack of data for 2014–2018). The chart does not include data for 2013, because, except for July, there is no data a b Figure 11: Average monthly changes in carbon monoxide concentration, averaged over the site 54–55∘ N, 83–84∘ E. Figure 12: Average annual CO concentrations for 2014–2020. 362 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 for May for this year. The trend line with the coefficient of determination 𝑅2 = 0.95 indicates a stable decrease in the CO concentration for 7 years from 2014 to 2020. In 2014, the average annual CO concentration is 0.145 ppm = 0.168 mg/m3 . The average annual value for 2020 is 0.136 ppm = 0.158 mg/m3 . The MAC (Maximum Allowable Concentration) for CO is the average daily MACad = 0.05 mg/m3 and the MAC single MACs = 0.15 mg/m3 [13]. 2.3.3. Seasonal and interannual variations of CH4 Methane enters the atmosphere as a result of both natural and anthropogenic processes. The main natural source of methane emissions into the atmosphere is wetlands. Methane is formed primarily as a result of the activity of bacteria involved in the decomposition of organic matter. Therefore, the main places of methane formation are swamps, garbage dumps, etc. In addition, methane is released through cracks in the earth’s crust at the bottom of water bodies, during mining operations, and during the burning of forests. Anthropogenic sources of methane in the atmosphere are energy production from coal and natural gas, landfills, sewage treatment plants, agricultural (manure processing), rice farming, and biomass burning. Methane is released during the production and transportation of coal, natural gas, and oil. After the release, methane remains in the atmosphere for approximately 8.4 years, after which it is removed mainly by chemical oxidation in the troposphere. The main mechanism for removing methane from the atmosphere is its oxidation in the upper atmosphere by the hydroxyl radical OH, which is formed by sunlight from ozone and water vapor. Figure 13 shows a histogram of the average monthly values of CH4 concentrations and SD value, based on the daily average data of the AIRS (Atmospheric Infrared Sounder) infrared spectrometer installed on the EOS AQUA satellite, AIRS3STD v.7.0 product. These are maps of daily methane distributions with a resolution of 1 × 1∘ at a baric altitude of 1000 hPa, available at the link [2]. There are no values for the summer months in the product. According to Figure 13, an increase in the methane concentration was observed for the autumn months (September–November). So the maximum value for September is 1.99 ppm, for October 1.97 ppm, for November 1.966 ppm. The decline in values was noted in April–May (there is no data for the summer months). The minimum value for May is 1.87 ppm. The maximum values of the SD do not exceed 3% of the average monthly values of methane concentration. The average concentration of CH4 in the atmosphere is 1.86 ppm [14]. For the autumn months, Figure 13: Monthly methane concentrations and SD values for 2013–2020. 363 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 Figure 14: Average annual values of CH4 concentration (no summer months). the methane concentration for the study area exceeds the average value for the atmosphere. Figure 14 shows a graph of the average annual values of the CH4 concentration averaged over 8 months (January–April, September–December), where the trend line with the coefficient of determination 𝑅2 = 0.79 shows a steady increase in the methane concentration over the period 2013–2020. For 2013, the average methane concentration is 1.92 ppm = 1.3 mg/m3 , and for 2020, 1.957 ppm = 1.32 mg/m3 . 3. Conclusions According to satellite monitoring data of the Iskitimsky district with a high technogenic load, the following changes have been outlined in the ecology of the region for the period 2013–2020. The dirty snow zone expanded eastward by 3–5 km over the period 2013–2020 in the area of Linevo and Kolyvansky and Vostochny coal mines. A decrease in the NDVI values was noted in these areas. It is shown the dynamics of changes in the aerosol optical thickness (AOT), CO and CH4 values in the Iskitim district. The seasonal change of AOT with the maximum values in April and July is shown. Only for 2018, the average AOT value is close to the average annual global value of 0.16. For the rest of the dates for the Iskitim region, the average long-term value of AOT is significantly higher than average annual global value. The highest concentrations of CO occur in the winter months, as well as in the spring months of March and April. The summer months are characterized by the lowest concentrations of CO. The CO concentration values in the autumn months are higher than in the summer months, but lower than in the winter months. The trend line with the coefficient of determination 𝑅2 = 0.95 indicates a stable decrease in the CO concentration for 7 years from 2014 to 2020, and increase in the methane concentration with the coefficient of determination 𝑅2 = 0.79. Acknowledgments The work was carried out within the framework of the state assignment of the Kotel’nikov FIRE RAS on the topic 0030-2019-0008 “Space”. 364 Natalia V. Rodionova CEUR Workshop Proceedings 353–365 References [1] https://ngs.ru/text/gorod/2020/09/28/69479455. [2] http://giovanni.gsfc.nasa.gov. [3] https://earthexplorer.usgs.gov. [4] https://sentinel.esa.int/web/sentinel/toolboxes/sentinel-1. [5] Purkis S., Klemas V. Remote sensing and global environmental change. UK, USA: Wiley- Blackwell, 2011. [6] Hall D.K., Riggs G.A., Salomonson V.V. Development of methods for mapping global snow cover using Moderate Resolution Imaging Spectroradiometer (MODIS) data // Remote Sensing of Environment. 1995. Vol. 54(2). P. 127–140. [7] Oparin V.N., Potapov V.P., Giniyatullina O.L. O kompleksnoj ocenke sostoyaniya okruzhayushchej sredy po dannym distancionnogo zondirovaniya Zemli v regionah s vysokoj tekhnogennoj nagruzkoj // Fiziko-Tekhnicheskie Problemy Razrabotki Poleznyh Iskopaemyh. 2014. No. 6. P. 199–209. [8] Bryson R.A., Goodman B.M. Volcanic activity and climatic changes // Science. 1980. Vol. 207. P. 1041–1044. [9] Sakerin S.M., Gorbarenko E.V., Kabanov D.M., Osobennosti mnogoletnei izmenchivosti aerozol’noi opticheskoi tolshchiny atmosfery i otsenki vliyaniya razlichnykh faktorov // Optika Atmosfery i Okeana. 2008. Vol. 21. No. 7. P. 625–631. [10] Tashchilin M.A., Yakovleva I.P., Sakerin S.M. Prostranstvenno-vremennye variacii aero- zol’noj opticheskoj tolshchi v Bajkal’skom regione // Sovremennye Problemy Distan- cionnogo Zondirovaniya Zemli iz Kosmosa. 2021. Vol. 18. No. 1. P. 219–226. [11] Chubarova N., Nezval Ye., Sviridenkov M., Smirnov A., Slutsker I. Smoke aerosol and its radiative effects during extreme fire event over Central Russia in summer 2010 // Atmospheric Measurement Techniques. 2012. Vol. 5. P. 557–568. [12] Deryugina A.B. Analiz lokal’nyh izmerenij koncentracii CO v g. Petergof. 2017. Available at: http://elib.rshu.ru/files_books/pdf/rid_8cfe04d7c0304f5587627e49f20d451a.PDF. [13] https://vozdyx.ru/article/pdk-vrednyx-veshhestv-v-atmosfernom-vozduxe. [14] Repina I. Remote studies of methane concentrations and fluxes in the atmosphere. Available at: http://d33.infospace.ru/d33_conf/tarusa2015/07.pdf. 365