3D Spatial Analysis of Particulate Matter (PM 10 , PM 2.5 and PM 1.0 ) and Gaseous Pollutants (H 2 S, SO 2 and VOC) in Urban Areas Surrounding a Large Heat and Power Plant

: In many regions of the world, the winter period is a time of poor air quality, due primarily to the increased use of individual and district heating systems. As a consequence, the atmospheric air contains increased concentrations of both particulate matter and gaseous pollutants (as a result of “low” emissions at altitudes of up to 40 m and “high” emissions more than 40 m above ground level). In winter, the increased pollution is very often exacerbated by meteorological conditions, including air temperature, pressure, air speed, wind direction, and thermal inversion. Here, we analyze the concentrations of particulate matter (PM 10 , PM 2.5 , and PM 1.0 ) and gaseous pollutants (H 2 S, SO 2 , and VOC) in the immediate vicinity of a large solid fuel-ﬁred heat and power plant located in an urban agglomeration. Two locations were selected for analysis. The ﬁrst was close to an air quality measurement station in the center of a multi-family housing estate. The second was the intersection of two main communication routes. To determine the impact of “low” and “high” emissions on air quality, the selected pollutants were measured at heights of between 2 and 50 m using an unmanned aerial vehicle. The results were compared with permissible standards for the concentration of pollutants. Temperature inversion was found to have a strong inﬂuence on the level of pollutants at various heights, with higher concentrations of particulate matter registered at altitudes above 40 m. The source of PM, H 2 S, and SO 2 pollutants was conﬁrmed to be “low emission” from local transport, industrial plant areas, and the housing estate comprising detached houses located in the vicinity of the measuring points. “High emission” was found to be responsible for the high concentrations of VOC at altitudes of more than 40 m above the intersection and in the area of the housing estate.


Introduction
Air quality depends on the volume of pollutant emissions, the intensity and type of physico-chemical changes occurring in the atmosphere, and the large-scale movement of air pollutants. Pollutants occur in the air for natural reasons, independent of humans, and from anthropogenic sources. Anthropogenic pollutants occur at high concentrations in highly urbanized areas, where there is often high population density, and, thus, have a particularly negative impact on human health and quality of life. Increased levels of pollutants are mainly caused by "low emissions" [1], mainly from road transport and household and municipal waste, including from heating individual homes [2,3].
Among the most dangerous pollutants emitted into the atmosphere are particulate matter pollutants (PM 10 , PM 2.5 , PM 1.0 ) and gaseous pollutants: (sulphur dioxide, carbon dioxide, carbon monoxide, hydrogen sulphide, or odors) [4][5][6][7]. Particulate matter pollutants have a negative impact on human health, both directly by penetrating the body (causing allergies and lung diseases), and indirectly by acting as carriers for heavy metals, movement of pollutants [26,27]. The transport of pollutants is influenced by the layout and coefficient of the shape of streets, the geometry of the buildings, the source of pollution (height and location), and the direction and speed of the wind [25]. Due to the complexity of the phenomena and the number of parameters that should be taken into account, field studies are a very important source of information [28][29][30][31]. The purpose of this study was to determine spatial changes in the state of air pollution in a city in relation to sources of "low" and "high" emissions.
Currently, unmanned aerial vehicles are increasingly used for various purposes, such as mapping, transport, monitoring. The development and miniaturization of air quality measuring devices allowed the use of UAVs for their transport and measurement of pollutant concentrations in the air table [32][33][34]. Currently, research on this subject is carried out in many research centers. Most of them concern the stage of UAV and measuring equipment construction [35,36]. In the article, the authors present the previously unpublished use of UAV with measuring equipment, as well as the method of analyzing and interpreting the obtained results concerning 3D spatial air quality analysis.

Methodology Measurement Area
Measurements were made in the city of Lodz, which is the third largest settlement unit in Poland (Central-Eastern Europe) in terms of the number of inhabitants (population 677,286, population density 2309.6 people/km 2 ). The area of analysis covers the area of a municipal communication junction, which is located next to a large housing estate in the immediate vicinity of a solid fuel-fired heat and power plant. The critical points in this area are a housing estate with an air quality measurement station owned and run by the Voivodship Inspectors of Environmental Protection (VIEP) (Figure 1 point 1), the intersection of Puszkina street and Przybyszewskiego street ( Figure 1, point 2), and the EC4 heat and power plant ( Figure 1, point 3). The VIEP station (station ID: PL0096A) is a container station that automatically measures the parameters for PM 10 , PM 2.5 , and O 3 . Around the air quality measurement station is an area of multi-family housing, characterized by a predominance of five-story residential buildings (15-20 m high) and a smaller number of 10-storey residential buildings (30-40 m high) (133,855 inhabitants, 1.45 km from the heat and power plant). The selected area allowed for the assessment of the impact of "low" and "high" emissions on residential areas. Moreover, the nearby air quality measurement station made it possible to verify the measurement results. The second critical point is an important north-south and east-west communication junction for residential and industrial areas, and an access point to the A1 motorway. It is a single-level intersection of multi-lane roads with two-lane circular traffic, located between the high emitter (EC4) and the air quality measurement point (VIEP station). The total measurement area is about 2 ha, surrounded from the north by collective housing estates and from the south by a cemetery, as well as industrial and commercial areas (in the north-south direction, 0.8 km from the main emitter of the heat and power plant). The second selected area, which is a communication node, made it possible to refer to linear emissions from road transport, and as the area of CHP plants located closer to it, it allows for the assessment of the direct impact of EC4 on pollutant emissions. The gross development intensity index for the area ranges from 0.05 to 1.0, which defines "the intensity of development of concentrated single-family housing, areas of large housing estates, areas of large industrial plants and areas of old suburban buildings". Thus, both the VIEP station and the intersection are in the zone of influence of the pollutant emitter, i.e., the EC4 heat and power plant. of light fuel. This translates into average emissions of 9.86 kg/h of PM10, 4.23 kg/h of PM2.5, and 250 kg/h of SO2 (data from January to March 2019). The maximum average annual emissions from the plant of PM10, PM2.5 and SO2 at 1.5 m above ground level are 0.013 μg/m 3 , 0.005 μg/m 3 , and 0.827 μg/m 3 , respectively. Measurements were made at three levels (20,35, and 50 m) above the intersection, using an DJI Matrice 600 Pro unmanned aerial vehicle (UAV) (MOT < 25 kg, diameter approximately 2 m) equipped with instruments for measuring particulate matter (PM10, PM2.5, PM1) and gas pollutants (H2S, SO2, and VOC). The minimum height of 20 m was chosen to ensure safe flight conditions over the traffic route, street lighting, and tram lines. The maximum height of 50 m was prescribed by the controlled traffic region (CTR) around Wladyslaw Reymont airport in Lodz, located 9 km away. A grid of points was established in the area of the intersection, where the UAV made its measurements ( Figure 2B). The EC4 power plant is one of two coal-fired CHP plants operating in the city of Lodz, producing electricity in a cogeneration system (thermal power 820 MW, electric power 198 MW). The main fuel used in EC4 is hard coal. The plant uses approximately 530,000 Mg/year of hard coal, in addition to 395,000 Mg/year of biomass and 400 Mg/year of light fuel. This translates into average emissions of 9.86 kg/h of PM 10 , 4.23 kg/h of PM 2.5 , and 250 kg/h of SO 2 (data from January to March 2019). The maximum average annual emissions from the plant of PM 10 , PM 2.5 and SO 2 at 1.5 m above ground level are 0.013 µg/m 3 , 0.005 µg/m 3 , and 0.827 µg/m 3 , respectively.
Measurements were made at three levels (20,35, and 50 m) above the intersection, using an DJI Matrice 600 Pro unmanned aerial vehicle (UAV) (MOT < 25 kg, diameter approximately 2 m) equipped with instruments for measuring particulate matter (PM 10 , PM 2.5 , PM 1 ) and gas pollutants (H 2 S, SO 2 , and VOC). The minimum height of 20 m was chosen to ensure safe flight conditions over the traffic route, street lighting, and tram lines. The maximum height of 50 m was prescribed by the controlled traffic region (CTR) around Wladyslaw Reymont airport in Lodz, located 9 km away. A grid of points was established in the area of the intersection, where the UAV made its measurements ( Figure 2B).
The measurements at a height of 20 m reveal the state of air quality in the area of so-called "Low emission" (up to about 40 m relative height). The measurements taken at a height of 50 m represent the area of "high emissions". In the vicinity of the measuring station ( Figure 1, point 1), measurements were made at a height of 2 m, and then from 5 m to 50 m with a gradient of 5 m (Figure 2A). The measurements at a height of 20 m reveal the state of air quality in the area of socalled "Low emission" (up to about 40 m relative height). The measurements taken at a height of 50 m represent the area of "high emissions". In the vicinity of the measuring station ( Figure 1, point 1), measurements were made at a height of 2 m, and then from 5 m to 50 m with a gradient of 5 m (Figure 2A).
The measurement apparatus was an eight-sensor module, which enabled the quantification of particulate matter and gaseous pollutants in the atmospheric air. Using a laser scattered (LS) sensor, the module measured PM10, PM2.5, and PM1.0 (10,000 particles per second). Metal oxide semiconductor (MOS) type sensors measured the concentration of organic solvents (Ethanol, Iso-Butane, H2, 0-500 ppm) and odors (0.5-1000 ppm isobutanol). ElectroChemical (EC) type sensors measured H2S (3 ppb-1 ppm), O2 (0.20-100%) and SO2 (0.5-2000 ppm). The measuring apparatus was equipped with a probe, 1.5 m long, which eliminated the influence of the UAV on air turbulence and, therefore, any possible impact on the results (Figure 3). The measurement apparatus was an eight-sensor module, which enabled the quantification of particulate matter and gaseous pollutants in the atmospheric air. Using a laser scattered (LS) sensor, the module measured PM 10 , PM 2.5 , and PM 1.0 (10,000 particles per second). Metal oxide semiconductor (MOS) type sensors measured the concentration of organic solvents (Ethanol, Iso-Butane, H 2 , 0-500 ppm) and odors (0.5-1000 ppm isobutanol). ElectroChemical (EC) type sensors measured H 2 S (3 ppb-1 ppm), O 2 (0.20-100%) and SO 2 (0.5-2000 ppm). The measuring apparatus was equipped with a probe, 1.5 m long, which eliminated the influence of the UAV on air turbulence and, therefore, any possible impact on the results (Figure 3).
Validation of the measurement data of particulate matter was performed on the basis of data from an accredited measuring station VIEP (the method equivalent to the reference method), while the gaseous pollutants were validated in relation to the VEGA-GC microchromatograph (equipped with a thermal conductivity detector TCD, minimum concentration of 500 ppb (0.005 ppm)). The analysis was based on the results of the parameters of particulate matter PM 10 , PM 2.5 , PM 1.0 , as well as H 2 S, SO 2 and VOC. Measurement data from station (VIEP) was obtained from a publicly accessible internet database [37]. Validation of the measurement data of particulate matter was performed on the basis of data from an accredited measuring station VIEP (the method equivalent to the reference method), while the gaseous pollutants were validated in relation to the VEGA-GC microchromatograph (equipped with a thermal conductivity detector TCD, minimum concentration of 500 ppb (0.005 ppm)). The analysis was based on the results of the parameters of particulate matter PM10, PM2.5, PM1.0, as well as H2S, SO2 and VOC.
Measurement data from station (VIEP) was obtained from a publicly accessible internet database [37].

Analysis of Results
Data analysis was carried out using the ArcGIS Pro 2.6 program (Advanced license). Data containing the GPS coordinates of the measurement points from the measuring equipment were entered in the form of a vector layer in the EPSG: 2180 system. The frame of reference was changed to WGS 1984 Web Mercator for spatial analysis. The attributes of the points were the concentration values of the measured pollutant parameters. Using ArcGIS Pro software, the spatial distribution of pollutants was determined based on data interpolation. Interpolation was carried out using the Empirical Bayesian Kriging 3D method. The Kriging method has the advantage over other interpolation methods (Inverse Distance Weighting 'IDW', Triangulated Irregular Network 'TIN') that it treats the observed variable as a random variable [38,39]. The weighting factors are estimated by minimizing the sum of squared deviations for regression and using spatial autocorrelation (semi-variogram) of the examined feature. This improves the quality of spatial prediction.

Meteorological Conditions
Instantaneous measurements were made during the "winter period" (from 30 September 2020 to 1 April 2021). This period was chosen because it is commonly associated with increased air pollution, especially due to the combustion of solid fuels, especially hard coal. Poor air quality is one of the most significant problems in Polish cities (of the 50 most polluted cities in the EU, 36 are in Poland). The level of pollution in the "winter period" is strongly related to meteorological conditions, in particular air temperature,

Analysis of Results
Data analysis was carried out using the ArcGIS Pro 2.6 program (Advanced license). Data containing the GPS coordinates of the measurement points from the measuring equipment were entered in the form of a vector layer in the EPSG: 2180 system. The frame of reference was changed to WGS 1984 Web Mercator for spatial analysis. The attributes of the points were the concentration values of the measured pollutant parameters. Using ArcGIS Pro software, the spatial distribution of pollutants was determined based on data interpolation. Interpolation was carried out using the Empirical Bayesian Kriging 3D method. The Kriging method has the advantage over other interpolation methods (Inverse Distance Weighting 'IDW', Triangulated Irregular Network 'TIN') that it treats the observed variable as a random variable [38,39]. The weighting factors are estimated by minimizing the sum of squared deviations for regression and using spatial autocorrelation (semi-variogram) of the examined feature. This improves the quality of spatial prediction.

Meteorological Conditions
Instantaneous measurements were made during the "winter period" (from 30 September 2020 to 1 April 2021). This period was chosen because it is commonly associated with increased air pollution, especially due to the combustion of solid fuels, especially hard coal. Poor air quality is one of the most significant problems in Polish cities (of the 50 most polluted cities in the EU, 36 are in Poland). The level of pollution in the "winter period" is strongly related to meteorological conditions, in particular air temperature, wind direction and wind strength. The analysis presented here focuses on three representative measurement series, taken on 11 December 2020 (series A), 18 December 2020 (series B), and 8 January 2021 (series C). The three series of measurements reflect typical conditions for the period. The measurements were all taken between 9:30 am and 12:00 am. During measurement series A and B, the weather conditions were similarly characterized by fog, Energies 2021, 14, 4070 7 of 21 low cloud cover, high relative humidity (above 80%), and average temperatures above 0. They differed in the direction of the wind. During series A, the wind was blowing from the south-east-i.e., from the EC4 pollutant emitter. During series B, a south-westerly wind predominated-i.e., from the side of the cemetery, industrial areas and single-family houses. Series C was selected for analysis because it showed a similar wind direction to series B, and a similar wind speed to series A. However, the weather conditions differed from those in both series A and series B. On the day of the measurements and on preceding days, there were freezing temperatures, frost, light snowfall, and a low and high base of cloud cover. Figure 4 shows photos of the weather conditions taken by the on-board UAV camera. The meteorological data are summarized in Table 1. representative measurement series, taken on 11 December 2020 (series A), 18 December 2020 (series B), and 8 January 2021 (series C). The three series of measurements reflect typical conditions for the period. The measurements were all taken between 9:30 am and 12:00 am. During measurement series A and B, the weather conditions were similarly characterized by fog, low cloud cover, high relative humidity (above 80%), and average temperatures above 0. They differed in the direction of the wind. During series A, the wind was blowing from the south-east-i.e., from the EC4 pollutant emitter. During series B, a south-westerly wind predominated-i.e., from the side of the cemetery, industrial areas and single-family houses. Series C was selected for analysis because it showed a similar wind direction to series B, and a similar wind speed to series A. However, the weather conditions differed from those in both series A and series B. On the day of the measurements and on preceding days, there were freezing temperatures, frost, light snowfall, and a low and high base of cloud cover. Figure 4 shows photos of the weather conditions taken by the on-board UAV camera. The meteorological data are summarized in Table 1.

Results
Each series of measurements began at the voivodeship measuring station "1" (VIEP station). This was so that the values detected by its measuring equipment for PM 10 and PM 2.5 could be used as a reference ( Table 2). The percentage difference between the data obtained from our measurements and those of the VIEP station did not exceed the permissible level of ±10%. The average for all series was approximately ±2%.
The measurements performed by the station also allowed us to observe changes in the air quality around the area of dense multi-family housing located about 1.45 km south of the EC4 emitter and industrial facilities. According to the Chief Inspectorate for Environmental Protection [37], a PM 10 level below 20 µg/m 3 is very good; up to 50 µg/m 3 is good; up to 80 µg/m 3 is moderate; up to 110 µg/m 3 is poor; up to 150 µg/m 3 is bad; above 150 µg/m 3 is very bad. The distribution of PM 10 pollution observed in series C was different from those observed in series A and B. The highest concentrations, reaching 38.4 µg/m 3 , were recorded at heights from 25 to 35 m above ground level. Above 35 m and below 25 m, the PM 10 concentrations were about 30% lower and ranged from 21 to 33 µg/m 3 .
The concentration of PM 10 pollutants in series A and B were found to exceed the permissible level of 50 µg/m 3 at heights above 25 m ( Figure 5). The concentrations of PM 10 were on average 8% higher in series A than in series B. The highest concentration was 66.5 µg/m 3 , which is 133% of the PM 10 standard. Closer to the ground (<25 m), in both series A and B the dust content in the air was also above the permissible standard, with 5-perc values as high as 40.7 µg/m 3 (series A) and 37.7 µg/m 3 (series B).

Results
Each series of measurements began at the voivodeship measuring station "1" (VIEP station). This was so that the values detected by its measuring equipment for PM10 and PM2.5 could be used as a reference ( Table 2). The percentage difference between the data obtained from our measurements and those of the VIEP station did not exceed the permissible level of ±10%. The average for all series was approximately ±2%. The measurements performed by the station also allowed us to observe changes in the air quality around the area of dense multi-family housing located about 1.45 km south of the EC4 emitter and industrial facilities. According to the Chief Inspectorate for Environmental Protection [37], a PM10 level below 20 μg/m 3 is very good; up to 50 μg/m 3 is good; up to 80 μg/m 3 is moderate; up to 110 μg/m 3 is poor; up to 150 μg/m 3 is bad; above 150 μg/m 3 is very bad. The distribution of PM10 pollution observed in series C was different from those observed in series A and B. The highest concentrations, reaching 38.4 μg/m 3 , were recorded at heights from 25 to 35 m above ground level. Above 35 m and below 25 m, the PM10 concentrations were about 30% lower and ranged from 21 to 33 μg/m 3 .
The concentration of PM10 pollutants in series A and B were found to exceed the permissible level of 50 μg/m 3 at heights above 25 m ( Figure 5). The concentrations of PM10 were on average 8% higher in series A than in series B. The highest concentration was 66.5 μg/m 3 , which is 133% of the PM10 standard. Closer to the ground (<25 m), in both series A and B the dust content in the air was also above the permissible standard, with 5-perc values as high as 40.7 μg/m 3 (series A) and 37.7 μg/m 3 (series B).   The concentration of pollutants above 25 m was between 120% and 205% of the permissible PM 2.5 level of 25 µg/m 3 . Better air quality in terms of these pollutants was observed in series C. The highest concentrations of PM 2.5 and PM 1.0 , at 28 and 27 µg/m 3 , respectively, were measured 20-30 m above ground level. Interestingly, the lowest concentrations of PM 2.5 were recorded above 45 m (21 µg/m 3 on average). Below 10 m, the average was 25 µg/m 3 .    [40]. In series A and B, the influence of the EC4 emitter on air quality can be excluded, because the chimney outlet is located at a height of 260 m, which is significantly above the temperature inversion layer. However, in the case of series C, there was no clear inversion layer. This translated into the presence of a noticeable air layer in the range of 20 to 30 m, characterized by lower concentrations of particulate matter pollutants. This was probably due to the formation of turbulence and the local transport of pollutants at    [40]. In series A and B, the influence of the EC4 emitter on air quality can be excluded, because the chimney outlet is located at a height of 260 m, which is significantly above the temperature inversion layer. However, in the case of series C, there was no clear inversion layer. This translated into the presence of a noticeable air layer in the range of 20 to 30 m, characterized by lower concentrations of particulate matter pollutants. This was probably due to the formation of turbulence and the local transport of pollutants at  [40]. In series A and B, the influence of the EC4 emitter on air quality can be excluded, because the chimney outlet is located at a height of 260 m, which is significantly above the temperature inversion layer. However, in the case of series C, there was no clear inversion layer. This translated into the presence of a noticeable air layer in the range of 20 to 30 m, characterized by lower concentrations of particulate matter pollutants. This was probably due to the formation of turbulence and the local transport of pollutants at the height of the roofs of the surrounding buildings. At higher altitudes, the pollutants dispersed, as evidenced by the fact the lowest concentrations in series C were recorded at heights above 40 m.  The concentrations of gaseous pollutants were also measured in the area of VIEP stations. No clear differences in concentration versus height were observed in any of the measurement series (Table 3). This was probably due to the distance of the measurement site from the main sources of gaseous pollutants, allowing them to mix with air. The only reproducible pattern was the increased concentrations of both hydrogen sulphide and sulphur dioxide at heights below 10 m, by up to 30% compared to the average value. This probably resulted from air stagnation caused by the presence of buildings in the measurement area. In all measurement series, the concentration of pollutants at many of the measurement points was below the minimum measuring threshold of the equipment. Only the average concentration of hydrogen sulphide significantly exceeded acceptable levels. In series A and B, the average concentration of hydrogen sulphide was about 0.036 mg/m 3 , which is 77% higher than the permissible value of 0.02 mg/m 3 . In series C, the concentration was as much as 750% of the permissible level, at 0.170 mg/m 3 , although it should be emphasized that this is significantly below lethal levels (1400 mg/m 3 ). Overall, in series C The concentrations of gaseous pollutants were also measured in the area of VIEP stations. No clear differences in concentration versus height were observed in any of the measurement series (Table 3). This was probably due to the distance of the measurement site from the main sources of gaseous pollutants, allowing them to mix with air. The only reproducible pattern was the increased concentrations of both hydrogen sulphide and sulphur dioxide at heights below 10 m, by up to 30% compared to the average value. This probably resulted from air stagnation caused by the presence of buildings in the measurement area. In all measurement series, the concentration of pollutants at many of the measurement points was below the minimum measuring threshold of the equipment.
Only the average concentration of hydrogen sulphide significantly exceeded acceptable levels. In series A and B, the average concentration of hydrogen sulphide was about 0.036 mg/m 3 , which is 77% higher than the permissible value of 0.02 mg/m 3 . In series C, the concentration was as much as 750% of the permissible level, at 0.170 mg/m 3 , although it should be emphasized that this is significantly below lethal levels (1400 mg/m 3 ). Overall, in series C the mean concentration of H 2 S was five times higher than in series A and B and the concentration of SO 2 was twice as high. It can therefore be assumed that the results were influenced by the emission of pollutants from the EC4 CHP plant, especially since in series C there was no temperature inversion. The high concentrations of H 2 S were likely emitted from the numerous vents from sewage systems in the densely packed residential buildings that surround the vicinity of the measuring station. This supports the results of previous research [11]. After the tests within the VIEP station, measurements were taken at the intersection of Puszkina street and Przybyszewskiego street. Clear changes were observed in the spatial distribution of the pollutants on the plane and with changes in height. For example, in series A the highest concentration of PM 10 (83 µg/m 3 ) was measured at 50 m (average 63 µg/m 3 ), whereas the lowest average (52 µg/m 3 ) was measured at an altitude 35 m (Figure 9). Concentrations higher than 50 µg/m 3 occurred in series A from the windward direction-i.e., from the open space side. On the opposite side of the intersection, where there are multi-family buildings and a larger number of trees, the concentration was up to 50% lower. The lowest concentration of PM 10 was recorded in the upper right corner of the area of analysis, where there are large numbers of trees and a green area. The influence of this area was visible even at a height of 50 m. An identical relationship was observed in the other series of measurements (B and C). In both series B and C the wind was blowing from the direction of the cemetery and industrial facilities. In both series (B and C), the average concentration of PM pollutants was 35% lower than in series A (Table 4). It can be assumed that this was due to the roughness of the terrain in the cemetery, which is covered with tall trees that trapped or accumulated the pollutants. In series A, the wind was blowing from the EC4 side, and probably transported particulate matter pollutants directly from the hard coal storage area. In series B, no significant differences (±2%) were found between the mean concentration of PM 10 at different heights. This was probably due to the wind speed, which was below 6 km/h and did not cause either violent air turbulence or PM 10 transport. The highest concentration of PM 10 , about 40 µg/m 3 , was recorded in series B towards the boundary of the cemetery and the surrounding buildings. As noted in [41],  In series C, the wind was almost three times stronger than in series B, and the influence of the buildings on the PM10 concentration was even more visible (Figure 9). At a height of 20 m, the mean concentration of 40-50 μg/m 3 was predominantly in the area sheltered by buildings, while average concentrations below 40 μg/m 3 were measured only in the unprotected area. At a height of 35 m, practically the entire area of analysis was  In series C, the wind was almost three times stronger than in series B, and the influence of the buildings on the PM 10 concentration was even more visible (Figure 9). At a height of 20 m, the mean concentration of 40-50 µg/m 3 was predominantly in the area sheltered by buildings, while average concentrations below 40 µg/m 3 were measured only in the unprotected area. At a height of 35 m, practically the entire area of analysis was characterized by an average concentration of 40-50 µg/m 3 , probably due to the formation of vortices. At an altitude of 50 m, the PM 10 concentration dropped to below 40-30 µg/m 3 . This can be explained by the lack of thermal inversion in series C. In series A and B, when thermal inversion occurred, there was an accumulation of pollutants at higher altitudes.
The results for PM 2.5 and PM 1 were similar to those for PM 10 (Figures 10 and 11). A strong correlation (>0.50) was found between the concentrations of the particulate matter pollutants. In series A, the highest concentration of dust pollutants at 75 µg/m 3 was found windward from the heat and power plant and industrial areas, especially at a height of 50 m. Closer to the buildings, the PM content in the air fell by up to 40%, which may have been caused by the formation of vortices on the leeward side. As the measurement altitude decreased, the concentration of PM 2.5 and PM 1 fell by about 20% in comparison to the concentration at 50 m. This may indicate the accumulation of dust pollutants due to temperature inversion in the upper air layers. Peng [32] came to similar conclusions. In series B and C, the concentrations of PM 2.5 and PM 1 were about 30% lower than in series A. The content of pollutants in the air was below the permissible level of 25 µg/m 3 only in the wooded area. In other places, the concentration of PM 2.5 exceeded the permissible level, reaching 175% of the permissible maximum. Gao [43] reported concentrations of PM 2.5 during the winter period reaching up to 500% of the limit defined in WHO guidelines. This can be explained by the lack of thermal inversion in series C. In series A and B, when thermal inversion occurred, there was an accumulation of pollutants at higher altitudes.
The results for PM2.5 and PM1 were similar to those for PM10 (Figures 10 and 11). A strong correlation (>0.50) was found between the concentrations of the particulate matter pollutants. In series A, the highest concentration of dust pollutants at 75 μg/m 3 was found windward from the heat and power plant and industrial areas, especially at a height of 50 m. Closer to the buildings, the PM content in the air fell by up to 40%, which may have been caused by the formation of vortices on the leeward side. As the measurement altitude decreased, the concentration of PM2.5 and PM1 fell by about 20% in comparison to the concentration at 50 m. This may indicate the accumulation of dust pollutants due to temperature inversion in the upper air layers. Peng [32] came to similar conclusions. In series B and C, the concentrations of PM2.5 and PM1 were about 30% lower than in series A. The content of pollutants in the air was below the permissible level of 25 μg/m 3 only in the wooded area. In other places, the concentration of PM2.5 exceeded the permissible level, reaching 175% of the permissible maximum. Gao [43] reported concentrations of PM2.5 during the winter period reaching up to 500% of the limit defined in WHO guidelines.   Sources of gaseous pollutants include technological processes, fuel combustion and heavy transport. One of the contaminants analyzed was H2S (Figure 12). High concentrations of H2S were observed in series A, reaching 0.20 ppm (average 0.042 ppm) at a height of 50 m. At lower altitudes, the concentration decreased by 35% to an average value of 0.027 ppm at a height of 20 m. This can be explained by the wind direction and pollution accumulation due to thermal inversion.
In contrast to series A, in series B and C the highest concentrations of H2S (over 0.06 ppm) were recorded at 20 m ( Figure 12). This is consistent with the fact that H2S is heavier than air and, therefore, accumulates near the ground. High concentrations of H2S have also been observed near the ground by Cichowicz and Dobrzański [44]. An interesting phenomenon was the increased concentrations of H2S at the exits of the intersection. This can be seen by comparing the maps of the H2S distribution ( Figure 12) in series A and B with the layout of the roads at the intersection. Given the temperature inversion in series A and B, it can be assumed that the H2S in the air came from low emissions from the traffic rand industrial plants. Kourtidis et al. [45] suggest that a source of H2S may be emissions from car catalysts. It is also possible to exclude the CHP as a source of H2S, as in series C (no inversion) there was no increase in the concentration compared to series A and B. It should be emphasized that in all series the average concentration of H2S was alarmingly high, exceeding the permissible level of 0.02 mg/m 3 . In series A, the average concentration was 0.058 mg/m 3 (Table 4), which is almost 300% of the permissible level [14]. In series B, the average concentration was double the permissible level. In series C, it was about 40% higher than 0.02 mg/m 3 . Sources of gaseous pollutants include technological processes, fuel combustion and heavy transport. One of the contaminants analyzed was H 2 S (Figure 12). High concentrations of H 2 S were observed in series A, reaching 0.20 ppm (average 0.042 ppm) at a height of 50 m. At lower altitudes, the concentration decreased by 35% to an average value of 0.027 ppm at a height of 20 m. This can be explained by the wind direction and pollution accumulation due to thermal inversion.
In contrast to series A, in series B and C the highest concentrations of H 2 S (over 0.06 ppm) were recorded at 20 m ( Figure 12). This is consistent with the fact that H 2 S is heavier than air and, therefore, accumulates near the ground. High concentrations of H 2 S have also been observed near the ground by Cichowicz and Dobrzański [44]. An interesting phenomenon was the increased concentrations of H 2 S at the exits of the intersection. This can be seen by comparing the maps of the H 2 S distribution ( Figure 12) in series A and B with the layout of the roads at the intersection. Given the temperature inversion in series A and B, it can be assumed that the H 2 S in the air came from low emissions from the traffic rand industrial plants. Kourtidis et al. [45] suggest that a source of H 2 S may be emissions from car catalysts. It is also possible to exclude the CHP as a source of H 2 S, as in series C (no inversion) there was no increase in the concentration compared to series A and B. It should be emphasized that in all series the average concentration of H 2 S was alarmingly high, exceeding the permissible level of 0.02 mg/m 3 . In series A, the average concentration was 0.058 mg/m 3 (Table 4), which is almost 300% of the permissible level [14]. In series B, the average concentration was double the permissible level. In series C, it was about 40% higher than 0.02 mg/m 3 . Another of the analyzed gas pollutants was sulphur dioxide, which is also heavier than air ( Figure 13). In series B and C, high concentrations of SO2 (more than 0.20 ppm, which is approximately 0.57 mg/m 3 ) were measured at a height of 20 m. The average concentration fell at higher altitudes, by 17% and 25% in series B and C, respectively. However, in series A, where there was temperature inversion and a strong wind, the highest concentration (above 0.40 ppm) occurred mainly at a height of 50 m. Analyzing the area of the intersection, it can be seen that in series A and B there was a space with an increased SO2 concentration compared to the surroundings. In series A, this is the exit from the intersection towards the city center, whereas in series B it is the exit from the intersection towards the residential part of the city. Note that these areas are in line with the direction air travel from the intersection. Interestingly, in the absence of thermal inversion (series C) no such relationship was observed. In most of the crossing area, the permissible SO2 concentration of 0.35 mg/m 3 was not exceeded. Only at the exits from the intersection was the permissible concentration of SO2 exceeded, reaching about 130%. However, increased concentrations of SO2 are usually observed in the winter [46], probably due to the increased production of thermal energy. Another of the analyzed gas pollutants was sulphur dioxide, which is also heavier than air ( Figure 13). In series B and C, high concentrations of SO 2 (more than 0.20 ppm, which is approximately 0.57 mg/m 3 ) were measured at a height of 20 m. The average concentration fell at higher altitudes, by 17% and 25% in series B and C, respectively. However, in series A, where there was temperature inversion and a strong wind, the highest concentration (above 0.40 ppm) occurred mainly at a height of 50 m. Analyzing the area of the intersection, it can be seen that in series A and B there was a space with an increased SO 2 concentration compared to the surroundings. In series A, this is the exit from the intersection towards the city center, whereas in series B it is the exit from the intersection towards the residential part of the city. Note that these areas are in line with the direction air travel from the intersection. Interestingly, in the absence of thermal inversion (series C) no such relationship was observed. In most of the crossing area, the permissible SO 2 concentration of 0.35 mg/m 3 was not exceeded. Only at the exits from the intersection was the permissible concentration of SO 2 exceeded, reaching about 130%. However, increased concentrations of SO 2 are usually observed in the winter [46], probably due to the increased production of thermal energy. The final analyzed gas compound was VOC ( Figure 14). There are no specific permissible levels for VOC. The distributions of VOC observed in series A and B were completely different to those in series C. Under the meteorological conditions in series A and B, with the visible influence of thermal inversion, a very low VOC concentration of less than 0.004 ppm (0.014 mg/m 3 ) was recorded, but in most of the studied area the value was even lower, at less than 0.002 ppm. In series A and B, up to 50% higher concentrations were recorded at a height of 20 m, averaging 0.003 ppm at 35 and 50 m. The increased concentrations of VOC in series A and B were mainly associated with the communication routes, so the source of pollution was "local" low emission. The situation was completely different in series C (Figure 14), where the highest concentrations of up to 0.016 ppm (0.056 mg/m 3 ) were found at an altitude of 50 m. However, at 35 m and 20 m the concentration of VOC decreased up to 7-fold. This may be evidence of the migration of VOC pollutants as a result of high emissions, although it would only be possible to test such a hypothesis under stable temperature gradient conditions. The final analyzed gas compound was VOC ( Figure 14). There are no specific permissible levels for VOC. The distributions of VOC observed in series A and B were completely different to those in series C. Under the meteorological conditions in series A and B, with the visible influence of thermal inversion, a very low VOC concentration of less than 0.004 ppm (0.014 mg/m 3 ) was recorded, but in most of the studied area the value was even lower, at less than 0.002 ppm. In series A and B, up to 50% higher concentrations were recorded at a height of 20 m, averaging 0.003 ppm at 35 and 50 m. The increased concentrations of VOC in series A and B were mainly associated with the communication routes, so the source of pollution was "local" low emission. The situation was completely different in series C (Figure 14

Conclusions
This study set out to perform a spatial analysis of the distribution pollutants caused by "low" and "high" emissions. In series A and B, "low" emissions were associated with high concentrations of PM, SO2 and H2S at heights of 30-50 m. This was probably due mainly to the occurrence of temperature inversion. In series C there was no temperature inversion, which translated into lower concentrations of pollutants and a different spatial distribution to those in series A and B. The concentrations of pollutants were higher when temperature inversion occurred (Series A and B) than when there was no temperature inversion (Series C). This indicates that the emitters/chimneys of the coal-fired heat and power plant (EC4) had little effect on air quality in the analyzed area. It can be assumed that the EC4 CHP plant only contributed to high concentrations of VOC at altitudes above 40 m.
In series C, increased levels of suspended dust in the air were recorded at heights of 25 and 35 m above ground level. This may indicate an air corridor transporting PM pollution at this altitude. In all the measurement series, a strong correlation (above 0.5) was found between the PM10, PM2.5, and PM1 dust fractions. A particularly high correlation coefficient (close to 1) was calculated for PM2.5 and PM1 (Table 5). Similar observations were made by Li et al. (2015), obtaining a correlation of 0.71-0.77 between PM. In all three series, concentrations of PM were measured that exceeded permissible values by up to 175% (PM2.5). Presumably, this was related to season, with increased production of individual and industrial thermal energy during the winter period. In the housing estate, the concentration of gaseous pollutants remained the same in the height range from 2 to 50 m. This may be explained by turbulence from buildings, which mixed the gaseous pollutants with air, and by the distance from the source.

Conclusions
This study set out to perform a spatial analysis of the distribution pollutants caused by "low" and "high" emissions. In series A and B, "low" emissions were associated with high concentrations of PM, SO 2 and H 2 S at heights of 30-50 m. This was probably due mainly to the occurrence of temperature inversion. In series C there was no temperature inversion, which translated into lower concentrations of pollutants and a different spatial distribution to those in series A and B. The concentrations of pollutants were higher when temperature inversion occurred (Series A and B) than when there was no temperature inversion (Series C). This indicates that the emitters/chimneys of the coal-fired heat and power plant (EC4) had little effect on air quality in the analyzed area. It can be assumed that the EC4 CHP plant only contributed to high concentrations of VOC at altitudes above 40 m.
In series C, increased levels of suspended dust in the air were recorded at heights of 25 and 35 m above ground level. This may indicate an air corridor transporting PM pollution at this altitude. In all the measurement series, a strong correlation (above 0.5) was found between the PM 10 , PM 2.5 , and PM 1 dust fractions. A particularly high correlation coefficient (close to 1) was calculated for PM 2.5 and PM 1 ( Table 5). Similar observations were made by Li et al. (2015), obtaining a correlation of 0.71-0.77 between PM. In all three series, concentrations of PM were measured that exceeded permissible values by up to 175% (PM 2.5 ). Presumably, this was related to season, with increased production of individual and industrial thermal energy during the winter period. In the housing estate, the concentration of gaseous pollutants remained the same in the height range from 2 to 50 m. This may be explained by turbulence from buildings, which mixed the gaseous pollutants with air, and by the distance from the source. The measurements taken at the intersection of communication routes, located closer to the industrial areas, showed clear differences in terms of the distribution of gaseous pollutants at different heights. In the case of H 2 S and SO 2 , the distribution of concentrations was consistent with the road system and the orientation towards industrial areas. This indicates that these pollutants were caused by local low emissions. High VOC concentrations were only recorded at altitudes of more than 40 m within the intersection and housing estate areas in Series C (no temperature inversion), indicating that they originated from high emissions. The permissible level for H 2 S was exceeded by up to 300% both in the area of the housing estate and at the intersection. The concentrations of H 2 S were higher in the housing estate than in the area around the intersection. This was probably due to the numerous drainage vents located on the roofs of buildings on the housing estate. There was no clear correlation between the levels of gaseous pollutants in the three series. A weak correlation of 0.44 was found between H 2 S and SO 2 levels in the three series. There was also no correlation between dust and gaseous pollutants (correlation coefficients below 0.5).
Summing up, the spatial analysis of pollutants with the use of UAVs gives the opportunity to get acquainted with the distribution of pollutants concentration in a selected area in detail. It should be emphasized that it is not possible to perform such measurements using traditional, stationary measuring stations located at a constant height from the ground, with a low spacing density. The method of measurement and analysis proposed by the authors can be implemented all over the world, therefore, it seems to be an important supplement to the existing methods of air quality measurement. The research methods presented in the article are now successfully used by us in air quality analyzes in various areas of Poland. The most important advantage of the obtained results is their high resolution in a small area, directly answering the question of air quality. Funding: This study was conducted as part of the research project entitled "Spatial analysis of air pollution changes in the Lodz agglomeration (in Polish: Analiza przestrzenna zmian stanu zanieczyszczenia powietrza w aglomeracji łódzkiej)", which was co-financed by approximately 80% by the Provincial Fund for Environmental Protection and Water Management in Lodz (in Polish: Wojewódzki Fundusz OchronyŚrodowiska i Gospodarki Wodnej w Łodzi).
Institutional Review Board Statement: Not applicable.