Fedor Lisetskii This email address is being protected from spambots. You need JavaScript enabled to view it., Аndrei Borovlev


Federal-Regional Center of Aerospace and Surface Monitoring of the Objects and Natural Resources, Belgorod State National Research University, Belgorod 308015, Russia


Received: December 25, 2017
Revised: November 16, 2018
Accepted: January 29, 2019

Download Citation: ||https://doi.org/10.4209/aaqr.2017.12.0593  

  • Download: PDF

Cite this article:

Lisetskii, F. and Borovlev, А. (2019). Monitoring of Emission of Particulate Matter and Air Pollution using Lidar in Belgorod, Russia. Aerosol Air Qual. Res. 19: 504-515. https://doi.org/10.4209/aaqr.2017.12.0593


  • We studied emissions from all sources of air pollution in the industrial center.
  • Particles PM2.5 were detected beyond 3.5 km from high emission sources.
  • The ratios of PM10/TSP and PM2.5/TSP were obtained.
  • Monitoring emissions from high sources is promising when using multiwave lidars.
  • A new software for estimating PM10 and PM2.5 share in the industrial emissions flare.


Fine suspended particulate matter with an aerodynamic diameter smaller than 10 (РМ10) or 2.5 µm (РМ2.5) can be a dangerous air pollutant necessitating operational monitoring. Of the 1113 major Russian cities, however, only a few monitor industrial emissions of PM10 and PM2.5. Here, we develop an approach to using mobile multi-wave (1064, 532, and 355 nm) lidar to estimate the concentration of PM10 and PM2.5. This approach was implemented for Belgorod, where 1378 sources of air pollution with anthropogenic dust, primarily of carbonate composition, were registered. We have developed algorithms with seven stages of assessing the spatial distribution and monitoring of РМ10 and РМ2.5, which made it possible to establish that fine-mode particles from tall sources of cement and construction material production (pipes with a height of ≥ 50 m) contributed 39% of the total particulate matter emissions. Using GIS to map the fields of the total suspended particulate matter (TSP) and determining the ratios of РМ10/TSP and РМ2.5/TSP, excesses in РМ10 and РМ2.5 up to 2.5 and 2.8 times greater, respectively, than the maximum threshold limit were observed. Tall sources’ contribution to emissions increased in proportion to the distance from the source, resulting in 40–85% of the РМ10 and 43–91% of the PM2.5. We demonstrate how lidar can be applied to optimize a particulate matter emissions monitoring network for environmental policy making.

Keywords: Particulate matters; Urban air pollution; Tall sources; Lidar measurements.


A number of studies have shown that neither cloud cover nor precipitation can fully account for the observed temperature trends and that aerosols need to be considered (Portmann et al., 2009; Tosca et al., 2017). Fine suspended particulate matter (PM) with an aerodynamic diameter less than 10 µm (РМ10) and 2.5 µm (РМ2.5) are considered, with good reason, as significant pollutants, requiring systematic monitoring. In the USA and EU, threshold limit values (TLV) for РМ10 and PM2.5 have been in existence for many years. Economic growth in Asian countries has stimulated the active development of control methods for РМ for this region (Karimian et al., 2016; Lu et al., 2016; Munir et al., 2017). A special concern arises from enhanced mortality rates caused by the inhalation of РМ10 and РМ2.5 (Samoli et al., 2008). Atmospheric PM may adsorb a great number of pollutants from the air, which together may be ingested causing a toxic action, leading to development of cardiovascular, respiratory and cancerous diseases (Reche et al., 2012). РМ2.5 poses a more serious threat to health than РМ10. PM2.5 composition include water-soluble ions like SO42–, NO3, and NH4+ (Yang et al., 2017). In EU countries, 3% of mortality from cardiovascular pathology and 5% of mortality from lung cancer are caused by PM2.5, and generally 3.1 million deaths per annum are caused by air pollution from PM2.5 (Cohen, 2004). Moreover, there is now no certain threshold below which РМ2.5 would not pose some threat (European Union, 2008). A comparison of the current regulatory requirements for the content of РМ10 and РМ2.5 and also for the content of the total suspended PM (TSP) in the air is presented in Table 1.

Table 1. The Threshold Limit Values (TLV) and other criteria for the content of suspended particles in the air for the USA, Russia and on the recommendation of the WHO.

In 2010 in Russia hygienic standards (НS were introduced to establish TLV for РМ10 and PM2.5. However, Russian systems monitor a maximum one-time concentrations of pollutants in the air (i.e., the particulate content over a 20-minute averaging time). Presently, these methods for monitoring emissions into the air for PM10 and PM2.5 are insufficient. The subsequent accumulation of air pollutants caused by urban soils and other natural constituents is the basis for monitoring of urbanized ecosystems (Bulygin and Lisetskii, 1996). Therefore, new measurement approaches to assess compliance within these TLV for the protection of human health are necessary.

For many Russian industrial centers, including Belgorod, where monitoring of air pollution with PM10 and PM2.5 does not exist, monitoring of the spatial distribution of PM with a consideration to the degree of its dispersion in air and an assessment of total emissions is urgent (Sadovnikov et al., 2013). To solve the problem, we propose to integrate monitoring and the lidar method of probing, which gives a number of apparent advantages for operational non-contact collection of necessary data for aerosol concentration and emission sources in contrast to standard gravimetric methods.

Monitoring using dispersion models (calculated monitoring) implies regular efforts to determine the spatial and temporal characteristics of air pollution based on transfer and diffusion, while taking into account emissions sources and meteorological conditions. Mathematical models are used for monitoring primarily due to the fact that they are cheaper to deploy and use than instrument-based monitoring. For example, for a city the size of Belgorod, calculated monitoring will be 11.2 times cheaper (compared to servicing four monitoring stations) (Borovlev et al., 2013b).

Both single- and multi-wavelength lidars may be also used for these tasks. Probing at several wavelengths provides additional information of PM size and enables the discriminating of gas pollutions. In some instances, both single- and multi-wavelength lidars may be used (e.g., Behrendt et al., 2004; Weitkamp, 2005; Matvienko et al., 2014; Mammez et al., 2015). For example, such applications can be as atmospheric wind speed measurements (Lolli et al., 2013), aerosol detection and transport (Campbell et al., 2016), and gas concentration (Queißer et al., 2016; Dandocsi et al., 2017). The potential for adapting lidar measurements for the assessment of the spatial distribution of suspended PM is hypothesized here. This study aims to develop a methodological approach to using a mobile multi-wavelength lidar and mathematical simulation of PM distribution to control emissions, and the subsequent application of the results to an existing system of ecological control of the atmospheric boundary layer at monitoring stations in Belgorod. This allows us to obtain results of the pollutant distribution pattern over different functional zones of the industrial center, and a solution for improved efficiency in operational control over PM emissions (PM10, PM2.5) into the air of a city where cement and construction are primary particulate sources.

Achieving this goal ensures a consistent solution to the following main tasks: (1) creation of a fully functional database of all sources of air pollution emissions in the city of Belgorod for calculating monitoring; (2) identification of priority sources, which provide the greatest contribution to urban air pollution by particles of PM10 and PM2.5, based on the GIS analysis of the mathematical simulation of the geographical distribution of suspended PM; (3) developing proposals for the organization of lidar control of source emissions of PM10 and PM2.5, with recommendations for mandatory use of the lidar measurements in case of detection exceedances of reference surface concentrations of TSP at monitoring stations.



Most monitoring stations administered by the Russian Federal Service for Hydrometeorology and Environmental Monitoring (Roshydromet) in the territory of Belgorod, as well as in many Russian cities, are located near highways, which in most cases result in overestimated pollution indexes. Because of this, the effectiveness of the Belgorod monitoring network cannot be considered satisfactory because of its relatively low informative value, both qualitatively and spatially. Therefore, to analyze urban industrial aerosols, including PM10 and PM2.5, it is advisable to use monitoring methods that allow for characterizing sources of industrial emissions, which make the greatest local contribution (for the functional zones of the city). We refer to this as priority sources for environmental control. Thus, if active monitoring is linked to a centralized monitoring system that allows for quickly mitigating priority sources of air pollution, the practical use of monitoring data is significantly expanded, which will improve the efficiency of the existing network.

This research is focused on the development of such a methodological approach, which could support an integrative system of lidar measurements and results of PM modeling within the existing monitoring network. The mathematical model, which is based on the solution of the atmospheric diffusion equation (Berlyand, 1991), is the one approved procedure for Russian monitoring operation (OND-86, 1987). Maximum one-time pollutant concentrations in OND-86 are calculated under relatively unfavorable meteorological conditions, meaning the absence of precipitation that could induce particle fallout, and thus the characteristics of vertical and horizontal turbulent mixing are combined unfavorably (Berlyand, 1991).

In accordance with this model, maximum one-time pollutant concentrations (см, mg m–3) at the discharge of dust and air mixture from a single-point source is reached at a threshold wind speed uм at a distance xм from the source, determined by:   


where А is a coefficient that depends on the temperature stratification of the atmosphere, which determines the conditions for horizontal and vertical dispersion of pollutants in the air, М is the mass of pollutants that is released into the air in a unit of time (injection capacity; g s–1), F is a dimensionless coefficient that takes into account the sedimentation rate of pollutants (gaseous and aerosols, including particulate matter) in the air, m and n are dimensionless coefficients that take into account the conditions of emission at the source, η is a dimensionless coefficient that takes into account the influence of the earth surface (geomorphology), H is the height of the emission source (m), V1 is the dust and air mixture consumption (m3 s–1), and ΔT is the difference between the temperature of the gas-air mixture and the air temperature (°С).

Lidar profiling provides valuable information in addition to in situ measurements, since it allows for continuous long-term measurements with relatively high spatial and temporal resolution. Multi-wavelength lidar measurements contain information not only about the spatial distribution of the aerosol, but also about its microphysical parameters, such as the size distribution of aerosol particles and the refractive index of the aerosol substance (Weitkamp, 2005). The problem of determining the aerosol parameters using the backscatter coefficients at several wavelengths, represents a class of inverse problems that are solved using algorithms (Burton et al., 2014).

The optical parameters of the aerosol are found from the solution of the lidar equation (Measures, 1992), which connects the measured signals at the wavelength λk (k = 1, …, Nλ) with the optical parameters of the atmosphere, as:

where F(λkz) = A(λkPL(λkzz2A(λk) is the hardware constant, z is distance along the sounding path, PL(λkz) is the lidar signal, βa(λkz) and βm(λkz) are the aerosol and molecular backscattering coefficients, and σa(λkz) and σm(λkz) are the attenuation coefficients from aerosol and molecular scattering.

The solution of Eq. (2) for a certain wavelength, λk, is obtained using transformations such as Fernald-Klett (Fernald, 1984; Klett, 1985), which can be written as follows (to reduce the notation for Eq. (3), we omit the variable λk):

Here t is an arbitrary point on the sounding path, θa(z) = βa(z)/σa(z) and θm = βm(z)/σm(z) are values that are called “lidar ratios” for aerosol and constant molecular backscatter to extinction coefficient (a synonym is the inverse ratios) (Hinkley, 1976; Measures, 1992). Point zg is the anchor point in which βa(λkz) is assumed to be known. For t ≤ z, the value of T(t, z) is interpreted as the transparency of the atmosphere in the segment of the path (tz).

For assessment of spatial distribution of suspended PM (taking into consideration their degree of dispersion), and an arrangement of monitoring of emissions of РМ10 and РМ2.5, an algorithm was suggested, which includes seven stages of study:

  1. Identification of priority sources, that is, those that have a determining contribution to total emissions of solid pollutants into the air, using source parameters in the created database.
  2. Determination of size distribution of anthropogenic dust for these sources (using a focused-beam ion/electron microscope).
  3. Measurements with the use of a combined semi-automated dust meter and a laser particle size analyzer of ground-level concentrations of TSP, PM10 and PM5 particles in the functional zones of the city for determination of the relations between РМ10/TSP and РМ2.5/TSP (conversion factors).
  4. Numerical modeling of surface concentrations of TSP for the warm (April–November) and cold (January–March, December) seasons.
  5. GIS-mapping of TSP concentrations in the atmospheric boundary layer to obtain РМ10 and РМ5 concentration fields, taking into account the coefficients РМ10/TSP and РМ2.5/TSP for recalculating the TSP concentration for fine particles.
  6. Analysis of the modeled results and GIS-mapping of surface concentrations of PM10 and PM5 along the residential area of the city, and the identification of priority sources for controlling emissions of PM10 and PM2.5, as well as the benchmark of surface concentrations of TSP at the pollution monitoring stations of the urbanized territory.
  7. Lidar monitoring of emissions of РМ10 and РМ5 from local tall sources (pipes with a height of ≥ 50 m) is conducted when benchmark pollution levels of surface TSP at the pollution monitoring stations are exceeded.


Software for evaluating air pollution, widely used in Russia and other Commonwealth of Independent States members, was developed on the basis of OND-86. These, specifically, include the program “Ecologist” (Firm “Integral”, St. Petersburg) that was used to derive estimates in our study. The regulatory dispersion model in OND-86 simulates the field of annual upper limits (98th percentiles) of short-term concentrations. On 1 January 2018 in Russia, a new regulatory document, “Method for calculation of diffusion of emissions of harmful substances (pollutants) into the air” (MCD), was introduced replacing OND-86. MNRE (2017) contain formulae for calculating the maximum one-time mean and annual mean PM concentrations, which allows them to be used to assess long-term effects of air pollution on the environment and to assess and minimize health risks. However, the model for the calculation of the maximum one-time ground-level concentrations, including the emissions of a single-point source of the fixed height, is the same as in OND-86.

The size distribution of dust in emissions from tall sources has been determined by experiment using a 3D focused-beam ion/electron microscope (https://www.fei.com/products/
sem/quanta-sem). Measurements of TSP mass concentration were carried out by gravimetric method with the use of a combined semi-automated dust meter, and a laser particle size analyzer was used for measurements of concentration of РМ10 and РМ2.5. Furthermore, the following methods have been applied in the study: cartographic analysis, math-and-stats methods, mathematical simulation and systemic functional analysis. The analysis of spatial distribution of ground-level concentrations of РМ10 and РМ2.5 regionally were carried out by methods of overlay, interpolation and spatial statistics and with the use of GIS-technologies.

Study Area and Data Sources

The subjects of the study were sources of air pollution located within the Belgorod agglomeration, having the center town Belgorod with an area of 15,600 ha and a population of 385,000. The Belgorod agglomeration includes 600,000 people total, and is a capital of the region with intensive development of iron ore (i.e., Kursk Magnetic Anomaly) where environmental problems were previously recorded (Lisetskii et al., 2011), rational land use has been studied (Goleusov and Lisetsky, 2014) and impacts have been reported on water use (Lisetskii et al., 2014; Yermolaev et al., 2015). However, the high level of air pollution resulting from emissions of combustible products (caused by explosions in open pits) and metallurgy is also a problem.

In Belgorod alone, there exists mainly construction materials production factories, like for cement and chalkstone (Borovlev et al., 2013a). The annual emission of pollutants in Belgorod is 63,600 t or 3.7 kg per person. The share of stationary sources in the volume of emissions is 11.5%, and the share of motor vehicles is 88.5%. About 98.0% of the volume of emissions from stationary sources comes from 56 large industrial enterprises (construction materials production factories, machine-building and metal-working factories, food industry factories and others). Gaseous and liquid fractions dominate emissions (97%). As for solid pollutants, the main contribution to the contamination of the air basin is by emissions from stationary sources (98%). Of these, 75.6% comes from cement (“Belgorod cement”, Closed Joint-Stock Company [CJSC]).

We developed a database of pollution sources in Belgorod, which was used for the analysis of the distribution of volumes of emissions of PM. This database includes 2726 stationary sources, 1378 of which are sources of solid pollutants (including emissions from 11 large construction materials production factories), and 256 motorways with traffic intensity of more than 300 cars per hour (for peak traffic).

Three stationary observation sites (in system of Roshydromet) are located at distances of 50–60 m from highways and 3500–4000 m from the building materials and cement production, facilities which in most cases cause local overestimation of pollution by gaseous pollutants and underestimation by the amount of solid particles.

Table 2. Technical characteristics of the lidar МWL-60-МОB

A mobile lidar system (МWL-60-МОB) mounted to a motor vehicle was applied. Lidar can operate in two modes of measurements: multi-wavelength modes (i.e., the measurements are made at wavelengths 355, 532 and 1064 nm) and polarization mode (at 532 nm). The main technical characteristics of the lidar, МWL-60-МОB, are given in Table 2. The source transmission unit of the lidar is a 1064 nm laser, which is doubled and tripled into separate channels. The 532 nm channel has a Glan-laser prism installed to calibrate the polarized receiving channels. The receiving telescope and transmitter are installed on an azimuth rotating unit, allowing for scanning at angles from 10 to 90° and in azimuth from 0 to 360°. The goal of monitoring the rate of dust emissions is to determine the integral of the mass (volume) concentrations of aerosol, η, across an arbitrary cross-section of the plume and then multiply the result by the wind speed, U. The resulting value found in this manner is the pollutant flow, which for a stationary case will be equal to the rate or intensity of the emission M, as:

 = Uη(xz)dz,                                                      (4)

where η is the linear integral of the concentrations of particles along the segment of the probing path z that passes through the plume along the x-axis and perpendicular to the middle direction of the wind.

The aerosol concentration is estimated based on the results of lidar measurements of its integral parameters. Multi-wavelength sensing, which makes it possible to determine the attenuation ratios for several wavelengths, was used to account for the variation in the micro-physical characteristics of the aerosol. This type of sensing was complemented by polarized profiling to identify the shape of the particles (Lidar, 2005; Matvienko et al., 2014; Chew et al., 2016). Eq. (3) in the case of multi-wavelength sensing is solved for all wavelengths jointly by an iterative method (Beard, 1974), and lidar relations are determined in the solution process. In order to determine the emission rate for finely dispersed particles, the volume concentrations of

the aerosol were estimated with an upper limit on the particle size, V2.5 and V10. The wind speed was determined based on measured data from a meteorological station, using the known links between wind speeds at the ground level and at an altitude (Nieuwstadt and Van Dop, 1982).


The prevailing influence on pollution in the air above the city Belgorod was exerted by seven sources over 50 m high from the construction materials industry (“tall” sources according to OND-86). They included: six stacks from cement production factories, with heights of 65–96 m, and one stack from a construction materials factory, with a height of 65 m. The contribution to solid pollutant emissions from stationary sources was about 41.0%. Of them, emissions from “tall” sources representing cement production factories reach 39.0% relative to all stationary sources. The size distributions of dust in the emissions from these tall sources, using the focused-beam ion/electron microscope showed that the share of РМ2.5 in the air emissions from the cement production factories was 42.0%, and the share of РМ10 was 51.0%. For emissions from the construction materials production factories, the share of РМ2.5 was 45.0%, and the share of РМ10 was 50.0%. Dust composition showed that the fraction of fine particles from construction materials tall sources was 39% relative to all PM.

Ground-level concentrations of TSP, РМ10 and РМ2.5 from twelve monitoring points in different functional zones around the city (residential, recreational, industrial, transport infrastructure) were available. During a calendar year (2015), 87 samples were taken for measurements of mass concentration of TSP by gravimetric method and by a semi-automated dust meter, and for measurements of РМ10 and РМ2.5 using a DustTrak Aerosol Monitor. The use of gravimetric method is necessary for calibration of the optical device, inspection of accuracy of its operation and determination of TSP. Of these, 15 samples characterized dust emissions when the dust collectors are turned off and the CO concentration reaches very high, dangerously explosive levels. The maximum emission of solid particles from tall sources located at altitude increases from 5 to 335 g s–1.

It has been suggested (Kim et al., 2006; Li et al., 2017; Munir, 2017; Xu et al., 2017) to use PM2.5/PM10 ratios as a relative indicator of fine particles. We analyzed 72 samples that characterized the standard emissions situation at the production facilities. The following ratios were obtained: РМ10/TSP = 0.59 and РМ2.5/TSP = 0.39. In this case, excess of TLVms and annual mean concentrations of TLV (TLVam) for РМ10 and РМ2.5 were not detected. These ratios were characterized by certain stability and varied depending on TSP concentration, which gives grounds for simplification of the procedure for assessment of concentrations of РМ10 and РМ2.5 by means of their conversion from the determined TSP concentration. It should be noted that the ratios РМ10/TSP and РМ2.5/TSP were derived from measurements featuring cement emissions. For such cases, excesses of 1.2–1.8 times and 1.5–2.2 times were obtained for TLVms РМ10 and TLVms РМ2.5, respectively. For this situation, the majority of particles of size > 10 µm occurred on the underlying surface at a distance of up to 0.5 km from the emission sources, and the proportion of particles of size < 10 µm increased with distance (Table 3). The typical duration of salvo cement emissions at the factories is 40–90 min at a frequency of 20–40 times a year, which negates the differences in the fractional composition of the air in the situations we looked at. This means we can use the coefficients 0.59 and 0.39 we obtained in the simplified procedure for estimating concentrations of РМ10 and РМ2.5 from the concentration of TSP we determined.

Table 3. Distribution of suspended particles in the atmospheric air of Belgorod, 2015.

The results of estimation-based monitoring showed that the maximum one-off concentrations of TSP in residential Belgorod occur at levels of 0.3–2.3 TLVms, while the maximum annual average concentrations occur at 0.2–0.9 TLVam. The highest level of local pollution of TSP at 1.7–2.3 TLVms for residential areas was recorded in the summer (Fig. 1). GIS-mapping of the ground-level TSP concentration fields, taking into account the ratio of РМ10/TSP and РМ2.5/TSP, established maximum one-off concentrations of РМ10 at 0.39–2.3 TLVms and РМ2.5 at 0.4–2.8 TLVms and an average annual concentration of РМ10 at 0.3–2.0 TLVam and РМ2.5 at 0.4–2.2 TLVam. By results of GIS-analysis, a sector rigid 73.5 ha (2.2% of the residential area) in the central district of the city was evaluated (Fig. 2). Excesses over TLVms were caused, mainly, by tall sources of cement emissions (Fig. 2(f)). Their contribution to the values in excess of the TLVms was over 50%.

Fig. 1. Zones on the territory of the city of Belgorod, where there are excesses of the maximum single TLVms of PM10 and PM2.5: 1: 1.0 TLVms of PM10; 2: 1.0 TLVms of PM2.5. F: CJSC “Belgorod cement”: A: industrial zone; B: residential zone; C: woodland; F: cement production factory (CJSC “Belgorod cement”); — isolines of air pollution of TSP (the share of TLVms).Fig. 1. Zones on the territory of the city of Belgorod, where there are excesses of the maximum single TLVms of PM10 and PM2.5: 1: 1.0 TLVms of PM10; 2: 1.0 TLVms of PM2.5. F: CJSC “Belgorod cement”: A: industrial zone; B: residential zone; C: woodland; F: cement production factory (CJSC “Belgorod cement”); — isolines of air pollution of TSP (the share of TLVms).

Fig. 2. Maximum surface concentrations of TSP in the city of Belgorod (the warm period of the year), the share of TLVms: A: zone of elevation of TLVms; B: zone of non-elevation of TLVms; C: woodland; D: recreational area; — isolines of air pollution of TSP (the share of TLVms).Fig. 2. Maximum surface concentrations of TSP in the city of Belgorod (the warm period of the year), the share of TLVms: A: zone of elevation of TLVms; B: zone of non-elevation of TLVms; C: woodland; D: recreational area; — isolines of air pollution of TSP (the share of TLVms).

To solve the environmental problems caused by cement emissions, new standards were developed in 2017 that take into account targets for reducing dust emissions from tall sources (15–30 m). Modeling estimates that we carried out indicate that the measures put in place to reduce the emissions of cement production have succeeded at reducing concentrations of pollutants in the residential areas of the city to 0.8 TLVms and 0.7 TLVam for РМ10 and to 0.9 TLVms and 0.8 TLVam for РМ2.5. At the same time tall source emissions of fine PM should be closely watched. They make a sizable contribution to air pollution of fine PM, increasing in proportion to the distance from the source, resulting in 40–85% of the РМ10 and 43–91% of the PM2.5. The relevance of monitoring of emissions of fine PM with tall sources is also borne out by the fact that the target emissions for these sources can be exceeded due to dust collector malfunctions.

These air pollution estimates can be significantly changed for the worse, given the trend towards harmonization of environmental standards in Russia with the EU. So, we will experience exceedances of average annual mean concentration for РМ10 by 1.4 times and РМ2.5 by 1.7 times if we use those TLVam that are accepted in the EU.

Information on compliance with PM emission standards is evaluated by comparison of TSP concentrations, which are measured maximum surface concentrations (reference concentrations) at stationary observation sites. Measured reference concentrations for the warm and cold seasons are shown in Table 4. If the TSP concentration exceeds the reference one, the violation of the established emission standards is counted. It is therefore necessary to carry out complementary monitoring using lidar of tall source emissions for solid pollutants into the air, which contribute greatly to the value of the reference measurements.

Table 4. Reference surface concentrations of TSP, PM10, PM2.5 at the stationary observation sites at environmental control of PM10 and PM2.5 emissions from CJSC “Belgorod cement”

In Belgorod city, such monitoring was conducted for tall source emissions of industrial cement using lidar, which was located on high ground at a distance about 1.5 km from sources at the aforementioned production facility. The average temperature of the dust and gas mix that was coming out of the stacks into the atmosphere, as per the production process, was 100°С. As a result of mixing with cold air, water vapor would condense resulting in the emergence of a visible plume. Profiling was performed in the turbulent diffusion zone of the mixture outside the observable plume. The distribution of the smoke stacks and residential structures in the area determined the minimal angle of the place (height) for the laser beam. Depending on the azimuth of the observations, this height was 30–50 m. When the visible plume was sensed with a laser beam, a significant contribution was apparent in the back scattering by water vapor particles. For this reason, in order to estimate the integral parameters of the solid aerosol emissions, sensing was carried out under the lower edge of the visible plume in the turbulent mixture diffusion area.

Fig. 3(a) presents a schematic diagram of the lidar measurements. The lidar is located at the point L, and sensing directions LA and LB differ in azimuth at the fixed position angle. Thus, the route passes LB in the area that is free of emissions from pipes and gives a background signal of the inverse dispersion PB(R) from a relatively “clean atmosphere,” depending on the distance R along the path. The sensing route passes LA under the plume and yields a signal, PS(R).

Fig. 3. The scheme for conducting lidar measurements to determine the optical thickness of the flare: (a) measurement scheme: sounding under the flare along the LA line gives the signal PS, probing the “clean atmosphere” along the line LB gives the signal PB; (b) the determination of the optical thickness (shown by arrows) from the comparison of signals PS and PB.
Fig. 3. The scheme for conducting lidar measurements to determine the optical thickness of the flare: (a) measurement scheme: sounding under the flare along the LA line gives the signal PS, probing the “clean atmosphere” along the line LB gives the signal PB; (b) the determination of the optical thickness (shown by arrows) from the comparison of signals PS and PB.

Fig. 3(b) demonstrates the method for determining the under-plume aerosol optical thickness and its mean aerosol extinction coefficient, and shows the adjusted square distances of the logarithms of the signals PB and PS conditionally. The signals for the lines АВ and CD are:

(AB)  lnPs(R) = lnPB(R) + K + δ(lnPB(R)) and
(CD)  lnPs(R) = 2τS+ lnPB(R) + K + δ(lnPB(R))              (5)

where PB(R) and PS(R) are signals of the inverse dispersion from “clean atmosphere” and under the plume respectively with range RτS is the optical thickness of the aerosol “trace” under the plume, K is independent from R constant, and δ(ln PB(R)) is signal fluctuations caused by error of signals record and small-scale fluctuations of optical characteristics of background aerosol environment. In the case of a small fluctuations, Eq. (5) provides the ratio for determination of τS as:

 = 0, 5[〈lnPS(R) – lnPB(R)⟩AB] – [〈lnPS(R) – lnPB(R)⟩CD],                         (6)

where the symbol 〈 ⟩ denotes averaging over segments AB or CD.

The procedure for estimating the percentage of liquid and solid aerosols in the industrial emissions involves the use of the polarized lidar signal. Polarized sensing was done for three typical cases: clean air, a visually observable aerosol plume coming out of the smoke stack, and an aerosol trail under the plume. The results indicate that the observable part of the plume is a dense aerosol comprised primarily of spherical water vapor droplets. Outside the visible plume, there is aerosol comprised of non-spherical particles that are spreading as a result of turbulent diffusion. The signal depolarization from the aerosol was 39–40%.

Measurements of the optical thicknesses on three wavelengths (355, 532, and 1064 nm) allowed us to construct a particle size distribution function and then estimate the integral parameters of the aerosol, including particle radius r32 (µm), area of particles per unit volume of air S (µm2 cm3), volume of particles per unit volume of air V (µm3 cm3), and aerosol extinction coefficient at a wavelength of 532 nm (km–1). The integral parameters were determined for dry and liberally moisturized states of the dust (refractive indices of n = 1.55–0.005i and n = 1.41–0.002i, respectively).

Analysis of the experimental data showed that the average mass concentration varied by no more than 15%, which means that variations in the refractive index are not critical relative to the aerosol mass concentration estimates for r32 in the area 0.6–1.2 µm. The difference obtained for the specific surface area of the particles (S) was even less. Microscopic analysis of the dust samples showed that most of the particles had a size of 2–4 µm and elongated shapes (or radii not > 2 µm when replaced with equivalent spheres), which agrees satisfactorily with the results of the lidar measurements.

On the whole, the results of the study we conducted suggest that with three wavelength sensing measurements we can estimate the integral parameters of aerosol, Sσ (532), V and r32, on condition that the refractive index n and its error ± 0.05 are known in advance. In this case, the aforementioned parameters can be estimated with an error of no more than 25%. Analysis of our data, obtained for five recent years of lidar observations, makes it possible to come up with the following general conclusions:

  • РМ5 are detected at a distance of over 3.5 km from tall emission sources. If the wind speed exceeds 2 m s–1, the width of the airborne plumes may increase from 50 to 700 m;
  • The main air pollutant in the residential zone of Belgorod is industrial cement production artifact, and the concentration of suspended PM at the height of 100–300 m reaches 0.003 µg m–3;
  • The size distribution of dominant PM from tall sources of the cement production factories is within the range of 2–4 µm.

The idea of controlling emissions of industrial by products from tall sources (detection of excess aerosol emissions) is based on taking lidar measurements when benchmark concentrations of TSP are exceeded at monitoring stations for specific wind directions (from tall emission sources to the monitoring stations).

For processing our lidar measurements of fine suspended PM in the plume of industrial emission, a special software package (Borovlev and Kungurtsev, 2008; Borovlev et al., 2013b) was used, which included:

  • Using a series of horizontal contours and a utility for said application for determining the distribution of concentrations of suspended PM in the industrial plume (an example is shown mapping the distribution of РМ10 concentrations on a map of Belgorod in Fig. 4);
  • Determining the integral of the mass (volume) concentration of РМ10 and РМ5 through a cross-section of the plume and then multiplying it by the wind speed, as well as getting estimates for emissions rates of РМ10 and РМ2.5 and comparing them with the benchmark emission rate.

Fig. 4. Distribution of concentrations of РМ10 (C·103, µg m–3) at a height of 52 m using the results of lidar measurements. Sub-base: vector map of the territory of Belgorod.Fig. 4. Distribution of concentrations of РМ10 (C·103, µg m–3) at a height of 52 m using the results of lidar measurements. Sub-base: vector map of the territory of Belgorod.

Thus, an integrative approach has been suggested for using mobile multi-wave (1064, 532, and 355 nm) lidar (МWL-60) and specialized software in estimating the emission rates of PM10 and PM2.5, and for adapting lidar measurements to the monitoring emission of TSP by gravimetric method common in Russia. As a result, standard and actual values of emissions of РМ10 and РМ2.5 from tall sources, obtained using laser based monitoring and mathematical simulation, may be compared. For the task of improving the promptness of monitoring PM emissions, a methodological approach is suggested, which is based on comparison of the results of lidar measurements and mathematical simulation of distribution of suspended PM concentrations.


In this paper uses an example of air pollution assessment for an industrial urban center of Russia (Belgorod) to demonstrate the importance of monitoring PM10 and РМ2.5 emissions relative to established standards for air quality. When these emissions remain unmonitored (which is typical for many industrial centers in Russia), capturing the spatial distribution of suspended PM, taking into consideration its degree of aerial dispersion and the available monitoring within the city, including lidar technologies, becomes urgent. Hence, an algorithm for estimating the spatial distribution of the РМ10 and РМ2.5 based on numerical modeling was developed.

The results indicate that fine suspended particles from tall sources (≥ 50 m) contributed 39% of the total solid pollutants emitted from stationary sources in Belgorod. Factories producing cement and construction materials contributed 50–51% of the PM10 and 42–45% of the РМ2.5. Although the mean annual concentrations for РМ10 and РМ2.5 in the residential zone did not exceed the TLVam, the maximum TSP, caused by industrial and vehicular emissions, exceeded the TLVms by 1.7–2.3 times during the summer. The stable PM10/TSP (0.59) and РМ2.5/TSP (0.39) ratios simplified the procedure for assessing the PM10 and РМ2.5 levels via the TSP concentration, which was obtained by computational monitoring.

Lidar observations revealed that when emissions from high-production cement factories did not exceed the permitted levels, the dominant solid particles were 2–4 µm in diameter. Furthermore, a series of measurements established that in cases of salvo emissions and exceedances, most particles that were > 10 µm fell on the ground approximately 0.5 km from the factory, and the fraction composed by these particles increased with the distance from the source.

Excessive РМ10 and PM2.5 was detected by comparing TSP concentrations observed at three stationary observation sites with the maximum measured values (reference concentrations). When the measured concentration exceeded the reference value, it was deemed necessary to conduct operational monitoring, first to measure emissions from the main sources of solid pollutants and then to assess the risk of pollutant deposition in the residential area of the city. Lidar is a promising technology for monitoring emissions from tall sources. In our study, we used mobile multi-wave lidar (МWL-60) and special software. The use of three-wavelength (1064, 532, and 355 nm) and polarized (532 nm) lidar sensing enabled us to evaluate the size distribution of PM in industrial plumes and to discern complex particle shapes, which is a good addition to the numerical simulations of pollutant distribution patterns over different functional zones of the industrial center.

It has been suggested that the effectiveness of the existing monitoring system could be boosted by comparing the TSP concentrations measured at stationary monitoring stations with benchmark concentrations calculated for warm and cold periods. This study was conducted to integrate the data from lidar observation of РМ10 and РМ2.5 with TSP monitoring systems, thus increasing the effectiveness of the latter.


This work was supported by the Ministry of Science and Education of the Russian Federation (the 2017–2019 state assignment for the Belgorod State National Research University, Project No. 5.4711.2017/6.7). The authors are very grateful to Dr. James Campbell and anonymous reviewers for their valuable comments and suggestions on our paper.

Aerosol Air Qual. Res. 19 :504 -515 . https://doi.org/10.4209/aaqr.2017.12.0593  

Don't forget to share this article 


Subscribe to our Newsletter 

Aerosol and Air Quality Research has published over 2,000 peer-reviewed articles. Enter your email address to receive latest updates and research articles to your inbox every second week.

Latest coronavirus research from Aerosol and Air Quality Research

2018 Impact Factor: 2.735

5-Year Impact Factor: 2.827

SCImago Journal & Country Rank

Aerosol and Air Quality Research (AAQR) is an independently-run non-profit journal, promotes submissions of high-quality research, and strives to be one of the leading aerosol and air quality open-access journals in the world.