Monitoring of Emission of Particulate Matter and Air Pollution using Lidar in Belgorod , Russia

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.


INTRODUCTION
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 PM 2.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 .PM 2.5 composition include water-soluble ions like SO 4 2-, NO 3 -, and NH 4 + (Yang et al., 2017).In EU countries, 3% of mortality from cardiovascular pathology and 5% of mortality from lung cancer are caused by PM 2.5 , and generally 3.1 million deaths per annum are caused by air pollution from PM 2.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.
In 2010 in Russia hygienic standards (НS 2.1.6.2604-10) were introduced to establish TLV for РМ 10 and PM 2.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 PM 10 and PM 2.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 PM 10 and PM 2.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 (PM 10 , PM 2.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 PM 10 and PM 2.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 PM 10 and PM 2.5 , with recommendations for mandatory use of the lidar measurements in case of detection exceedances of reference surface concentrations of TSP at monitoring stations.

Methodology
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 PM 10 and PM 2.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), V 1 is the dust and air mixture consumption (m 3 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 are the aerosol and molecular backscattering coefficients, and σ a (λ k , z) and σ m (λ k , z) 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 z g is the anchor point in which β a (λ k , z) 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 (t, z).
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 semiautomated dust meter and a laser particle size analyzer of ground-level concentrations of TSP, PM 10 and PM 2.5 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 РМ 2.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 PM 10 and PM 2.5 along the residential area of the city, and the identification of priority sources for controlling emissions of PM 10 and PM 2.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 РМ 2.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.

Methods
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 (98 th percentiles) of shortterm 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 onetime 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 focusedbeam 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, mathand-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.
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: 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 at wavelength 1064 nm 4 nm at wavelength 532 nm 2.5 nm at wavelength 355 nm 3.5 nm Signal processing mode analogue the aerosol were estimated with an upper limit on the particle size, V 2.5 and V 10 .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).

RESULTS AND DISCUSSION
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 PM 2.5 /PM 10 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 TLV ms and annual mean concentrations of TLV (TLV am ) 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.8times and 1.5-2.2times were obtained for TLV ms РМ 10 and TLV ms РМ 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.
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.3TLV ms , while the maximum annual average concentrations occur at 0.2-0.9TLV am .The highest level of local pollution of TSP at 1.7-2.3TLV ms 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.3TLV ms and РМ 2.5 at 0.4-2.8TLV ms and an average annual concentration of РМ 10 at 0.3-2.0TLV am and РМ 2.5 at 0.4-2.2TLV am .By results of GISanalysis, 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 TLV ms were caused, mainly, by tall sources of cement emissions (Fig. 2(f)).Their contribution to the values in excess of the TLV ms was over 50%.
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 TLV ms and 0.7 TLV am for РМ 10 and to 0.9 TLV ms and 0.8 TLV am 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 PM 2.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 TLV am 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.
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 P B (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, P S (R).
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 P B and P S conditionally.The signals for the lines АВ and CD are: where P B (R) and P S (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  In the case of a small fluctuations, Eq. ( 5) provides the ratio for determination of τ S as: 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 r 32 (µm), area of particles per unit volume of air S (µm 2 cm -3 ), volume of particles per unit volume of air V (µm 3 cm -3 ), 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.005iand 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 r 32 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 r 32 , 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: • РМ 2.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 РМ 2.5 through a crosssection 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.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 PM 10 and PM 2.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.

CONCLUSIONS
In this paper uses an example of air pollution assessment for an industrial urban center of Russia (Belgorod) to demonstrate the importance of monitoring PM 10 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 PM 10 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 TLV am , the maximum TSP, caused by industrial and vehicular emissions, exceeded the TLV ms by 1.7-2.3times during the summer.The stable PM 10 /TSP (0.59) and РМ 2.5 /TSP (0.39) ratios simplified the procedure for assessing the PM 10 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 PM 2.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 threewavelength (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.

Fig. 1 .
Fig.1.Zones on the territory of the city of Belgorod, where there are excesses of the maximum single TLV ms of PM 10 and PM 2.5 : 1: 1.0 TLV ms of PM 10 ; 2: 1.0 TLV ms of PM 2.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 TLV ms ).

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

Fig. 3 .
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 P S , probing the "clean atmosphere" along the line LB gives the signal P B ; (b) the determination of the optical thickness (shown by arrows) from the comparison of signals P S and P B .

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

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.

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

Table 4 .
Reference surface concentrations of TSP, PM 10 , PM 2.5 at the stationary observation sites at environmental control of PM 10 and PM 2.5 emissions from CJSC "Belgorod cement"