<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>Novosibirsk, Russia
" alexdergunov@icm.krasn.ru (A. V. Dergunov); oleg@icm.krasn.ru (O. E. Yakubailik)</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>Analysis of the influence of temperature inversions on the ecological situation in Krasnoyarsk</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Alexander V. Dergunov</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Oleg E. Yakubailik</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Federal Research Center Krasnoyarsk Science Center SB RAS</institution>
          ,
          <addr-line>Krasnoyarsk</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Institute of Computational Modelling SB RAS</institution>
          ,
          <addr-line>Krasnoyarsk</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2021</year>
      </pub-date>
      <volume>000</volume>
      <fpage>0</fpage>
      <lpage>0002</lpage>
      <abstract>
        <p>The paper analyzes the meteorological conditions in the city of Krasnoyarsk in the period from January 1, 2019, to December 31, 2020. The relationship between temperature inversions in the surface layer of the atmosphere and air pollution by suspended solid particles PM25 is investigated. The paper uses a set of meteorological data of the NCEP GFS weather forecast model on the air temperature on three isobaric surfaces of 1000, 925, and 850 Mb; on wind gusts and the height of the atmospheric boundary layer. Data on PM25 solid particle concentrations and wind speed are provided by the air monitoring system of the KSC SB RAS. The relationship between the presence of temperature inversions in the lower layer of the atmosphere and periods of significant air pollution is shown, as well as the dependence of changes in wind speed and the height of the boundary layer of the atmosphere with changes in the average daily PM25 concentration. The results of the data analysis allow us to conclude that there is a high degree of correlation between these parameters. The possibility of using the meteorological data of the NCEP GFS model to study the surface layer of the atmosphere and the periods of its pollution, predicting the deterioration of the environmental situation in Krasnoyarsk, is demonstrated.</p>
      </abstract>
      <kwd-group>
        <kwd>eol&gt;Temperature inversion</kwd>
        <kwd>unfavorable meteorological conditions</kwd>
        <kwd>PM25</kwd>
        <kwd>GFS</kwd>
        <kwd>air pollution</kwd>
        <kwd>Krasnoyarsk</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>
        Atmospheric pollution is one of the main problems of large cities in the world. As a result of
anthropogenic factors, many harmful substances that pollute the air enter the atmosphere [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ].
Studies show an association between increased concentrations of particulate matter (PM) in the
air and the deterioration of human health [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. During periods of unfavorable meteorological
conditions (UMC), there is a sharp increase in the concentration of pollutants in the lower layer
of the atmosphere. We are talking about industrial and automobile emissions, furnace heating,
and so on [
        <xref ref-type="bibr" rid="ref3 ref4">3, 4</xref>
        ]. One of the determining factors of UMC is the temperature inversion in the
atmosphere.
      </p>
      <p>
        Temperature inversions occur because layers of greater or lesser thickness are located at
diferent heights in the surface layer of the atmosphere. The temperature decrease slows
down very much, stops, or, conversely, instead of the temperature decreasing with altitude, it
increases. Their properties are the height of the lower and upper boundaries, vertical power
(layer thickness), and intensity (inversion value). Temperature inversions are divided into
3 types: surface (the lower boundary is located at ground level), elevated (the lower boundary
is located at some height from the surface), and inversions of the free atmosphere (the height
can vary greatly) [
        <xref ref-type="bibr" rid="ref5 ref6 ref7">5, 6, 7</xref>
        ].
      </p>
      <p>There is a need to study atmospheric processes with the help of meteorological data sets to
analyze periods leading to environmental degradation.</p>
      <p>Sets of meteorological data can be obtained, for example, from ground-based environmental
monitoring stations. The advantages of such weather information include a fairly high accuracy
of measurements on the spot. However, the disadvantages of such data are that basic parameters
limit such data sets, and the ground stations themselves can be located at a great distance from
each other, which imposes certain restrictions on the spatial coverage of the problem being
solved.</p>
      <p>
        An alternative option is the weather data of various weather forecast models, for example, the
Global Forecast System model. These sources provide many diferent layers of meteorological
information on dozens of vertical levels, plotted on regular rectangular grids covering the entire
Globe. The spatial resolution of such data can vary from 2.5 to 0.25 degrees and higher [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ].
Developers try to improve this parameter regularly.
      </p>
      <p>The purpose of this work is to analyze the influence of temperature inversions in the surface
layer of the atmosphere on the formation of an unfavorable environmental situation in the city
of Krasnoyarsk in 2019–2020 based on high-spatial-resolution meteorological data of the NCEP
Global Forecast System weather forecast model. Identification of the relationship between the
periods with the average daily concentrations of PM25 particulate matter suspended in the air,
exceeding the average daily maximum permissible concentrations (MPC), the speed, wind gusts,
and the height of the atmosphere’s boundary layer.</p>
      <p>The research area in this work is the city of Krasnoyarsk and part of the adjacent territories
(Figure 1).</p>
      <p>
        No. Dates
12 periods of UMC were established in the city of Krasnoyarsk from 2019 to 2020 according
to the Ministry of Ecology and Rational Nature Management data of the Krasnoyarsk
Territory [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ]. The periods of UMC are characterized by low wind speed and high particulate matter
concentrations, significantly exceeding the average daily MPC equal to 0.035 mg/m 3. The list of
UMC periods for the studied time period is presented in Table 1.
      </p>
    </sec>
    <sec id="sec-2">
      <title>2. Materials and methods</title>
      <p>
        The work uses meteorological information from the dataset of the Global Forecast System (GFS)
model. This is a numerical weather forecasting system containing a global computer model
and variational analysis performed by the US National Weather Service. This model combines
four separate models: atmosphere, ocean, land/soil, and sea ice. Dozens of atmospheric and
ground-soil variables are available in this data set, from temperature, wind, and precipitation, to
soil moisture indicators and atmospheric ozone concentration. This is one of the most famous
meteorological models in the world. Global data analysis and forecasting are carried out 4 times
a day. The weather forecast is available up to 16 days in advance [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ].
      </p>
      <p>The accuracy of the GFS model is constantly improving. In particular, data with a horizontal
resolution of 1 degree has been available since March 2004, with a resolution of 0.5 degrees
since January 2007. Since January 2015, the horizontal resolution has been 0.25 degrees (about
25 km at the latitude of the city of Krasnoyarsk).</p>
      <p>
        The GFS model data is presented in the grib2 format. This data format is standardized by the
World Meteorological Organization (WMO) and is intended for storing historical and forecast
weather data [
        <xref ref-type="bibr" rid="ref10">10</xref>
        ]. Each file contains more than 500 layers of various meteorological information
on more than a hundred vertical levels.
      </p>
      <p>This work used the actual analysis data on the air temperature at three vertical levels
corresponding to three isobaric surfaces: 1000, 925, and 850 Mb. Data on the height of the atmospheric
boundary layer and information on wind gusts from January 1, 2019, to December 31, 2020,
obtained from the GFS model data set, were also analyzed.</p>
      <p>During the preliminary processing of the meteorological information of the GFS model, the
data was loaded and cropped according to the specified coordinates corresponding to the city
of Krasnoyarsk, and the necessary layers with information about the air temperature at the
selected vertical levels, as well as about the height of the boundary layer of the atmosphere and
wind gusts. Since the spatial resolution of the data is 0.25 degrees, two cells of the regular grid
of the GFS model covered the entire city of Krasnoyarsk and part of its surroundings (Figure 1).</p>
      <p>
        Then the received data was converted to the tabular CSV format. These procedures were
performed using special scripts written in Python using the wgrib2 program. The wgrib2
program is specially developed by the manufacturer of the GFS model data for reading and
converting files of the grib2 format [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ].
      </p>
      <p>Inversion layers were detected by obtaining the diference between temperature data at
diferent vertical levels when the diference was negative. For example, with a negative diference
between temperatures on isobaric surfaces of 1000 and 925 Mb, the inversion was considered
ground-level or elevated; if between 925 and 850 Mb, it was elevated or high-altitude (inversion
of the free atmosphere). If the surface elevated and high-altitude inversions were recorded
simultaneously, this was considered a powerful inversion.</p>
      <p>We also used data on the concentrations of PM25 in the air and wind speed according to
ground-based monitoring stations obtained at the ICM SB RAS geoportal.</p>
      <p>Data on PM25 concentrations were averaged per day for all available observation posts in the
city. Thus, the average daily values of the PM25 concentration for the entire city of Krasnoyarsk
were obtained.</p>
    </sec>
    <sec id="sec-3">
      <title>3. Results and discussion</title>
      <p>As a result, an archive of meteorological information was created according to the GFS model for
the territory of the city of Krasnoyarsk for the period from 2019 to 2020. The archive includes
data on the air temperature on three isobaric surfaces of 1000, 925, and 850 Mb, the height
of the atmosphere’s boundary layer, and wind gusts. The archive also contains information
from ground-based monitoring stations on the average daily concentrations of PM25 particulate
matter suspended in the air and wind speed.</p>
      <sec id="sec-3-1">
        <title>3.1. Analysis of temperature inversions and periods of high atmospheric pollution in Krasnoyarsk</title>
        <p>Based on the analysis of air temperature data at various vertical levels, the number of days was
compared with temperature inversions by month for 2019 and 2020 (Figure 2, a). You can pay
attention to the fact that in April, May, and November, the number of days with inversions
is greater in 2020, in the remaining months in 2019, there are more days with temperature
inversions, except for March when both in 2019 and in 2020 their number was the same and
equal to 10 days.</p>
        <p>A comparison of the number of days during which the average daily concentration of PM25
exceeded the average daily MPC by month for 2019 and 2020 is shown in Figure 2, b. In 2019,
there were more days with pollution of the atmosphere than in 2020 by 15 days.</p>
        <p>Table 2 provided information on the total number of days with temperature inversions and
periods when the average daily concentration of PM25 exceeded the average daily MPC in 2019
and 2020.</p>
        <p>The number of days of the study period is 731 days. Of these, 29% are days with a temperature
inversion, and 14% are days with a high concentration of PM25 in the city’s atmosphere.</p>
        <p>The total duration of days in the UMC periods indicated in Table 1 was 47 days. Of these, the
days when PM25 pollution exceeded the average daily MPC is 38.</p>
        <p>Monthly comparison of the number of days characterized by an average daily concentration
of PM25 exceeding the average daily MPC and the number of days with the same characteristic,
but with the existence of a temperature inversion for 2019 and 2020, was carried out (Figure 3).
The months when increased air pollution was not recorded are not displayed.</p>
        <p>In 2019, increased atmospheric pollution was observed for 7 months; in 2020, during 5 months.
Data analysis showed that in 2019, the diference between the number of days characterized by
the average daily concentration of PM25, exceeding the average daily MPC, and the number of
days with the same characteristic, but with the existence of a temperature inversion, was 19 days.
A big contribution to this diference was made by July when 12 days with increased atmospheric
pollution were registered, while there were only 2 days with a temperature inversion in such a
period; this is due to non-regular reasons — smoke from remote forest fires that covered the
city at this time. In 2020, the diference was 10 days.</p>
      </sec>
      <sec id="sec-3-2">
        <title>3.2. The relationship of meteorological parameters with periods of high air pollution</title>
        <p>Data analysis showed a stable relationship between changes in the average daily concentration
of PM25 particulate matter suspended in the air and wind gusts according to the GFS model. As
an example, this relationship is presented for February 2020 (Figure 4).</p>
        <p>In Figure 4, the orange stripe indicates the period when the concentration of PM25 exceeded
the average daily MPC. The blue stripe indicates the oficial periods of UMC. Temperature
inversions were observed on all days in the selected periods, except for February 5 and 26, as
well as PM25 high concentrations. It is shown that when the wind speed increases, the level of
pollution decreases and vice versa.</p>
        <p>The level of correlation between the wind speed obtained from the data of the ground
monitoring station and wind gusts according to the GFS model is 0.82.</p>
        <p>
          The variation of the height of the planetary boundary layer of the atmosphere is compared
with the wind speed obtained from the data of the ground monitoring station. The correlation
coeficient was 0.79. This is because, at high wind speed, the atmosphere’s stability decreases,
turbulence increases, which leads to an increase in the height of the boundary layer [
          <xref ref-type="bibr" rid="ref12">12</xref>
          ].
        </p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>4. Conclusions</title>
      <p>As a result of the work done, the relationship between temperature inversion in the surface
layer of the atmosphere and the periods of significant air pollution by PM 25 solid particles in
the city of Krasnoyarsk in 2019 and 2020 was shown.</p>
      <p>The largest number of days with a temperature inversion is observed in the cold period
(December, January, February), the smallest — in the warm season (May, June, July). In 2019,
inversions were registered in all 12 months of the year, and in 2020, there were no inversions
only in June and July.</p>
      <p>The total number of days with an inversion of air temperature in 2019 was 110 and in 2020 —
103. The total number of days with inversion was 29% of the number of days of the entire study
period. The most common type of inversion is surface or elevated.</p>
      <p>The greatest number of days with the average daily concentrations of PM25 particulate
matter suspended in the air, exceeding the average daily MPC, were recorded in the cold
season (December, January, February). In April, May, June, September, and October, periods of
significant air pollution were not recorded in 2019 and 2020.</p>
      <p>The total duration of the days of the UMC periods was 47 days. Of these, the days when the
PM25 pollution exceeded the average daily MPC is 38.</p>
      <p>The total number of days with polluted air in Krasnoyarsk in 2019 is 58. Of these, 39 days
during the temperature inversion. In 2020 — 43 days, of which 33 days during the temperature
inversion, the total number of days with average daily concentrations of PM25 exceeding the
average daily MPC was 14% of the number of days of the entire study period.</p>
      <p>The relationship between changes in wind speed and the height of the atmospheric boundary
layer with changes in the average daily PM25 concentration in Krasnoyarsk in February 2020 is
shown. As the wind speed increases, the level of pollution decreases, which is also true in the
opposite direction.</p>
      <p>The correlation between the wind speed obtained from the data of the ground monitoring
station and the wind gusts obtained from the data of the GFS weather forecast model is 0.82.
The diference between the wind speed and the height of the boundary layer is 0.79. Data on the
height of the boundary layer can serve as an additional indicator when studying UMC periods
in Krasnoyarsk.</p>
      <p>The analysis of meteorological data of the GFS model contributes to solving problems related
to the study of the lower layer of the atmosphere, its pollution and can play an important role in
more accurately identifying periods of unfavorable meteorological conditions and forecasting
them.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>G.</given-names>
            <surname>Zarubin</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Novikov</surname>
          </string-name>
          ,
          <article-title>Hygiene of the city</article-title>
          , Medicine, Moscow,
          <year>1986</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>Y.</given-names>
            <surname>Kaufman</surname>
          </string-name>
          , D.
          <string-name>
            <surname>Tanré</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          <string-name>
            <surname>Boucher</surname>
          </string-name>
          ,
          <article-title>A satellite view of aerosols in the climate system</article-title>
          ,
          <source>Nature</source>
          <volume>419</volume>
          (
          <year>2002</year>
          )
          <fpage>215</fpage>
          -
          <lpage>223</lpage>
          . doi:
          <volume>10</volume>
          .1038/nature01091.
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>A.</given-names>
            <surname>Becker</surname>
          </string-name>
          , T. Agaev,
          <article-title>Protection and control of environmental pollution</article-title>
          , Gidrometeoizdat, Leningrad,
          <year>1989</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>A.</given-names>
            <surname>Eremkin</surname>
          </string-name>
          , I. Kvashnin,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Junkerov</surname>
          </string-name>
          ,
          <article-title>Standardization of emissions of pollutants into the atmosphere</article-title>
          ,
          <source>Publishing house of Ass. builds. universities</source>
          , Moscow,
          <year>2001</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>L.</given-names>
            <surname>Matveev</surname>
          </string-name>
          ,
          <article-title>General meteorology course</article-title>
          .
          <source>Atmospheric Physics</source>
          , Gidrometeoizdat, Leningrad,
          <year>1984</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>E.</given-names>
            <surname>Bezuglaya</surname>
          </string-name>
          , G. Rastorgueva,
          <string-name>
            <surname>I. Smirnova</surname>
          </string-name>
          ,
          <article-title>What does an industrial city breathe?</article-title>
          , Gidrometeoizdat, Leningrad,
          <year>1991</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>V.</given-names>
            <surname>Zuev</surname>
          </string-name>
          ,
          <string-name>
            <given-names>N.</given-names>
            <surname>Krasnenko</surname>
          </string-name>
          ,
          <string-name>
            <given-names>V.</given-names>
            <surname>Fedorov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Fursov</surname>
          </string-name>
          ,
          <article-title>Acoustic sounding of the atmospheric boundary layer</article-title>
          ,
          <source>Dokl. Akad. Nauk SSSR</source>
          <volume>257</volume>
          (
          <year>1981</year>
          )
          <fpage>1092</fpage>
          -
          <lpage>1096</lpage>
          . URL: http://mi.mathnet.ru/ eng/dan44384.
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>The</given-names>
            <surname>Global Forecast System (GFS) Documentation</surname>
          </string-name>
          ,
          <year>2021</year>
          . URL: https://is.gd/KJHKbK.
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          <article-title>[9] Unfavorable meteorological conditions in Krasnoyarsk. regional departmental information and analytical data system on the state of the environment of the Krasnoyarsk Territory</article-title>
          ,
          <year>2021</year>
          . URL: http://www.krasecology.ru/Nmu.
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          <source>[10] NCEP WMO GRIB2 Documentation, Version 23.0.0</source>
          ,
          <year>2019</year>
          . URL: https://www.nco.ncep. noaa.gov/pmb/docs/grib2/grib2_doc/.
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11] wgrib2:
          <article-title>Utility to read and write grib2 files</article-title>
          ,
          <year>2021</year>
          . URL: https://www.cpc.ncep.noaa.gov/ products/wesley/wgrib2/.
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <given-names>Y.</given-names>
            <surname>Wang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>H.</given-names>
            <surname>Sayit</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Mamtimin</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Zhu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>C.</given-names>
            <surname>Zhou</surname>
          </string-name>
          ,
          <string-name>
            <given-names>W.</given-names>
            <surname>Huo</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Yang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>X.</given-names>
            <surname>Yang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Gao</surname>
          </string-name>
          ,
          <string-name>
            <given-names>X.</given-names>
            <surname>Zhao</surname>
          </string-name>
          ,
          <article-title>Evaluation of five planetary boundary layer schemes in WRF over China's largest semi-fixed desert</article-title>
          ,
          <source>Atmospheric Research</source>
          <volume>256</volume>
          (
          <year>2021</year>
          ). doi:
          <volume>10</volume>
          .1016/j.atmosres.
          <year>2021</year>
          .
          <volume>105567</volume>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>