Columnar and surface urban aerosol in the Moscow megacity according to measurements and simulations with the COSMO-ART model

. Urban aerosol pollution was analyzed over the Moscow megacity region using the COSMO-ART (COSMO – COnsortium for Small-scale MOdelling, ART – Aerosols and Reactive Trace gases) online coupled mesoscale model system and intensive measurement campaigns at the Moscow State University Meteorological Observatory (MSU MO, 55.707 ◦ N, 37.522 ◦ E) during the April–May period in 2018 and 2019. We analyzed mass concentrations of particulate matter with diameters smaller than 10 µm (PM 10 ), black carbon (BC) and aerosol gas precursors (NO x , SO 2 , CH x ) as well as columnar aerosol parameters for ﬁne and coarse modes together with different meteorological parameters, including an index characterizing the intensity of particle dispersion (IPD). Both model and experimental datasets have shown a statistically signiﬁcant linear correlation of BC with NO 2 and PM 10 mass


Introduction
Anthropogenic aerosol pollution has a complex impact on the atmosphere, significantly affecting solar radiation, air temperature and humidity and resulting in noticeable climatic effects (IPCC, 2022(IPCC, , 2013Jacobson, 2004;Bond et al., 2013). Aerosol particles at the surface level also have a harmful effect on human health (Manisalidis et al., 2020;Lu et al., 2015). Radiative effects of the anthropogenic aerosol are negative and exceed 1 W m −2 , partially compensating for the increase in the air temperature in the troposphere due to the rise in the concentration of greenhouse gases. However, the uncertainty of aerosol climate impact estimates remains quite high (IPCC, 2021;Myhre et al., 2013). These uncertainties are associated with a wide variety of optical and microphysical characteristics of aerosol (Seinfeld and Pandis, 2016) and its significant temporal and spatial variations. Anthropogenic aerosol is considered to be smaller in size and is more absorbing than natural aerosol (Myhre, 2009;Su et al., 2013;Kinne et al., 2013).
The intensive aerosol studies concern the optical properties of urban aerosol, its relation with meteorological characteristics as well as the relationship between surface concentration and columnar aerosol content and the emission sources (Segura et al., 2017;Zhuang et al., 2018;Wang et al., , 2019Zhdanova et al., 2020). However, in most publications the authors consider natural and urban columnar aerosol in polluted areas without highlighting its urban component (Kumar et al., 2019;Chou et al., 2006;Zhuang et al., 2018;Segura et al., 2017). In only a few papers were the urban component of columnar aerosol optical depth (AOD) and other aerosol properties evaluated (Zawadska et al., 2013;Chubarova et al., 2011b;Zhdanova et al., 2020). At the same time, the detection of the urban aerosol component and its relationship with anthropogenic emissions of gas precursors are critical for assessing aerosol radiation forcing and its climate effect (IPCC, 2022).
When analyzing urban aerosol pollution, it is important to understand the relationship between the content of groundlevel aerosol and the aerosol content in the atmospheric column. However, such an analysis has been made only in a small number of publications (Segura et al., 2017;Gubanova et al., 2018;Wang et al., 2019).
Black carbon (BC) is a particularly important urban aerosol component which absorbs visible radiation and contributes to the heating of the atmosphere, in contrast to most other aerosol species (Jacobson, 2004(Jacobson, , 2006Ramanathan and Carmichael, 2008;Bond et al., 2013). The urban environment is the main source of BC emissions due to the use of diesel fuel (Weingartner et al., 1997). The contribution of these emissions from heavy vehicles (trucks, buses, etc.) with diesel engines can reach 42 % of the total mass of black carbon emissions into the atmosphere (Reddy and Venkataraman, 2002). To date, measurements of black carbon are irregular and quite rare, especially in cities, where the emissions are particularly high. This leads to a lack of understanding of the processes occurring in the atmosphere, and, in particular, to a possible error in estimating the contribution of BC to the balance between heating and cooling rates (Bond et al., 2013), which, in turn, increases the uncertainty of its climatic effect. The possible influence of BC on the absorbing properties of the atmosphere has been analyzed in only a few publications (Kozlov et al., 2016;Markowicz et al., 2017;Rajesh and Ramachandran, 2018). Thus, measurements and modeling of black carbon and the evaluation of its influence on the atmospheric properties are important tasks, especially in urban conditions (Gilardoni et al., 2011;Kozlov et al., 2008;Lugon et al., 2021;Tang et al., 2021). The quantification of urban aerosol and its properties is also important for improving the accuracy of the meteorological forecast of temperature and humidity of the atmosphere, which noticeably depends on aerosol amount and its characteristics (Toll et al., 2016;Wang et al., 2020;Huang and Ding, 2021;Poliukhov and Blinov, 2021).
Chemical transport models (CTMs) coupled with weather prediction models or with the data of re-analysis are used for studying aerosol pollution and its influence on meteorological parameters (Evans et al., 2003;WMO-COST, 2008;Vogel et al., 2010;Baklanov et al., 2017). However, for reliable estimations of aerosol parameters, the simulations should be carefully tested against measurements. For example, in Ukhov et al. (2020), the application of the WRF-Chem model over the Middle East and its examination against observations provided reliable assessment of pollution by mineral and sulfate aerosol over the urban area in this region.
According to AEROCOM modeling data and CMIP5 assessments, the anthropogenic component of AOD at a wavelength of 550 nm is 0.03 ± 0.01, which is 24 ± 6 % of the total AOD (IPCC, 2013). This is smaller compared with the satellite retrievals, which provide the estimates of about 0.06 (or 20 %-40 % of the total AOD) over land (Loeb and Su, 2010;Bellouin et al., 2013). According to the latest estimates, there is still an underestimation in simulated AOD of about 21 % (Gliß et al., 2021). Hence, the measurements of the urban aerosol component and the evaluation of its ratio in total AOD may provide a helpful test bed for aerosol urban modeling.
The Moscow megacity with its population of about 13 million people and with about 7 million vehicles is one of the largest urban agglomerations in the world. As the capital of the Russian Federation, it is a large financial and administrative center. Air pollutant levels in Moscow are consistent with those in other megacities in Europe and North America (Elansky, 2014;Elansky et al., 2018). BC measurements in the center of Moscow have revealed substantially lower levels of air pollution than in Beijing (Golitsyn et al., 2015). The levels of BC concentration are similar to those in European cities, which indicates the comparable impacts of major urban sources such as traffic with intensive implementation of modern environmental requirements, heating power plants and manufacturing industries (Popovicheva et al., 2020(Popovicheva et al., , 2022.
The main tasks of this paper concern the analysis of the aerosol properties at the surface and in the atmospheric column, their relationship with meteorological parameters in Moscow, and evaluation of aerosol pollution as the difference between the aerosol properties in the Moscow megacity and in pristine conditions using the results of the chemical transport model and the data obtained during the intensive measurement campaigns over the Moscow region in the spring periods of 2018 and 2019.

Measurements
The MSU MO is located in the territory of the MSU Botanical Garden in the park area at a distance of several kilometers from the local sources of emissions (power stations).
The nearest highways are about 300-450 m away from the site. During the intensive measurement campaigns, the mass concentrations of various gas aerosol species at the Earth's surface and aerosol characteristics in the total atmospheric column were studied together with meteorological observations. The gas and aerosol measurements at the surface consisted of mass concentrations of particulate matter with diameters smaller than 10 µm (PM 10 ) and BC, different aerosol gas precursors and other gas species (NO, NO 2 , SO 2 , volatile organic compounds -VOCs -marked as CH x , CO) with a 20 min time resolution. These observations were operated by the Mosecomonitoring State Environmental Protection Agency (2022). For PM 10 measurements, the TEOM 1400a (Thermo Environmental Instruments Inc., USA) was used. Internationally certified OPTEC Russian instruments (http://www.optec.ru, last access: 13 January 2022; OPTEC, 2022) were applied to measure gas species mass concentrations of NO, NO 2 , SO 2 and CO. The Gamma-ET instrument (http://etek-ltd.ru, last access: 26 January 2022; Gamma-ET, 2022) was used for the CH x measurements. The description of the quality assurance (QA) procedures is given at http://mosecom.mos.ru (last access: 26 January 2022).
Aerosol equivalent BC (eBC) mass concentrations were measured with 1 min resolution using a custom-made portable aethalometer (Popovicheva et al., 2017). With this instrument the light attenuation caused by the particles depositing on a quartz fiber was analyzed at three wavelengths (450, 550 and 650 nm). The eBC concentrations were determined by converting the time-resolved light attenuation to eBC mass at 650 nm and were characterized by a specific mean mass attenuation coefficient, as described in Popovicheva et al. (2017). The calibration parameter for quantification of eBC mass was derived during parallel longterm measurements against an AE33 aethalometer (Magee Scientific) that operates at the same three wavelengths. More details can be found in Popovicheva et al. (2020).
Aerosol measurements in the atmospheric column were carried out using the CIMEL sun/sky photometer, which has been operated at Moscow State University in the framework of the AERONET program since 2001 (Holben et al., 1998;Chubarova et al., 2011a). We analyzed AOD in the spectral range from 340 to 1020 nm, fine and coarse AOD modes at a wavelength of 500 nm (O'Neill et al., 2001), the Ångström extinction exponent (AEE) and the Ångström absorption exponent (AAE) in the spectral range of 440-870 nm, singlescattering albedo (SSA) at 675 nm and an asymmetry factor at 675 nm for various aerosol modes calculated in accordance with the AERONET algorithms (Dubovik and King, 2000). The latest version 3.0 level-2 AERONET dataset was used in the analysis (Giles et al., 2019). Its detailed testing revealed that the new algorithm of automatic cloud filtration, applied in this version, worked much better than the old one, with the exception of the winter months (Chubarova, 2020). As a result, we did not apply additional cloud filtering, as has been previously done (Chubarova et al., 2016). During daytime the dataset has 15 min resolution of the aerosol characteristics obtained from the direct sun measurements and 1 h resolution from sky measurements. For comparisons with model AOD at 550 nm (AOD550), we recalculated measured AOD at 500 nm (AOD500) to AOD550 using the AEE parameter.
We also used meteorological observations (air temperature, atmospheric pressure, wind speed, wind direction) with 1 min resolution from the Vaisala MAWS-301 automatic weather station as well as the standard meteorological MSU MO measurements with 3 h resolution. In addition, the index, characterizing the intensity of particle dispersion (IPD) proposed in Kuznetsova et al. (2014), has been applied in the analysis. This index is calculated using a set of meteorological parameters, including atmospheric pressure conditions, the type of atmospheric circulation, the stratification of the atmosphere, wind speed up to 850 hPa and level of precipitation. The IPD index varies from 1 to 3. The conditions with IPD = 1 are characterized by a stable stratification of the atmosphere, a low-gradient baric field, low wind speed conditions and the absence of precipitation. At IPD = 3 the opposite picture is observed with intensive air mixing conditions, high wind speed, precipitation, unstable stratification and frontal zones. In our study the 1 h resolution IPD index is evaluated using the 24 h COSMO mesoscale model forecast. Finally, all the data were combined in the 1 h resolution dataset.

COSMO-ART mesoscale model system and numerical experiments
The urban aerosol characteristics were calculated using the COSMO-Ru-ART model system (Vogel et al., 2010;Vil'fand et al., 2017;Rivin et al., 2019) with a horizontal grid step of 7 km and a time resolution of 40 s over the 1000 × 1000 km area. In this model system meteorological simulations were performed by the COSMO mesoscale model (http://www. cosmo-model.org/, last access: 27 January 2022; COSMO, 2022), an operational weather prediction model at the Russian Hydrometeorological Centre (Rivin et al., 2019).
Since the main objective of our study was to evaluate aerosol pollution in Moscow at the MSU MO and in upwind background conditions, we did not focus on detailed variations of aerosol inside the city and used a 7 km grid step for model simulations.
The gas aerosol concentrations are simulated using the ART chemical transport model, which is coupled with the COSMO model. The COSMO-Ru-ART mesoscale model system reproduces chemical transformations of substances in the gas-phase and heterogeneous reactions, photolysis, nucleation, coagulation, condensation, emissions of various types of aerosols and dry and wet aerosol deposition (Vogel et al., 2010;Vil'fand et al., 2017). About 172 chemical reactions are used in the ART model to describe chemical processes in the troposphere. One of the most important fea-tures of COSMO-Ru-ART is the parallel calculation of meteorological parameters and chemical transformations at each time step, which allows a user to take into account the reverse effect of aerosols on radiation and meteorological characteristics of the atmosphere.
The simulations of reactive gaseous and particulate matter are based on the enhanced KAMM/DRAIS/MADEsoot/dust model (Riemer et al., 2003;Vogel et al., 2006Vogel et al., , 2010. In MADEsoot (Modal Aerosol Dynamics Model for Europe extended by soot), all aerosol modes are represented by lognormal distributions. Five modes for Aitken and sub-micron aerosol particles include one pure soot mode, secondary particles consisting of sulfate, ammonium, nitrate, organic compounds (SOA), as well as the modes representing aged soot particles, consisting of sulfate, ammonium, nitrate, organic compounds, water, and soot. It also includes the coarse particle mode, which contains additional anthropogenic emitted particles. All aerosol fractions are subject to coagulation and condensation following Whitby et al. (1991), Kerminen and Wexler (1994), Binkowski and Shankar (1995), Odum et al. (1996) and Schell et al. (2001). The soot particles are directly emitted into the atmosphere. Coagulation and condensation are accounted for in transfer of soot from external to internal mixtures. The Aitken and sub-micron particles are formed due to the aging process. For each mode prognostics equations for the number density and the mass concentration are solved numerically. The standard deviations are kept constant. Since the number densities of the coarse mode are small, the inter-modal coagulation between the coarse mode and the other modes and the intra-modal coagulation of the coarse-mode particles are both neglected. Additionally, the aerosol distributions are modified by the sedimentation, advection and turbulent diffusion processes. More details can be found in Vogel et al. (2010). The resuspension of urban dust with stronger winds currently is not taken into account. The chemical reactions of the gaseous species are calculated using the chemical mechanism RADMKA (Regional Acid Deposition Model Version Karlsruhe) based on RADM2 (Regional Acid Deposition Model; Stockwell et al., 1990), with the important updates described in Vogel et al. (2010). The photolysis frequencies were simulated according to Vogel et al. (2009). For the evaluation of the aerosol optical properties (the extinction coefficient, the SSA and the asymmetry factor), a special parameterization scheme is used based on the a priori calculations with the application of the approach described in Bohren and Huffmann (1983) and precalculated aerosol distributions. This procedure is based on typical size distributions and chemical compositions, which are simulated in the model domain.
As a result, using the COSMO-Ru-ART model system, we can quantify the rate of formation of new aerosol particles and aerosol gas precursors in the polluted urban atmosphere in real atmospheric conditions, which in turn are modified by the updated chemical composition. The data of the forecast of the COSMO-Ru system (Rivin et al., 2019) and the global ICON model were used as initial and boundary conditions.
In addition, for the operation of the ART model, the data from the Global Land Cover 2000 project on land use and inventory data from TNO2010 (Netherlands Organization for Applied Scientific Research, https://www.tno.nl/en/, last access: 28 January 2022; TNO, 2022) were applied to determine anthropogenic emissions of pollutants. The 1 h resolution TNO2010 emission inventory has been developed using official reported emissions data by source category and combining them with other estimates where needed (Kuenen et al., 2014). Urban aerosol sources according to TNO2010 include direct emission of particulate matter of undefined composition (as a dust fraction) and soot as well as the gas aerosol precursors of sulfate, organic, and nitrate secondary aerosols. The model also accounts for relatively small natural biogenic emissions of non-methane volatile organic compounds from the Global Land Cover 2000 project, which are the gas precursors of organic aerosol. Testing the model estimates with the TNO2010 and TNO2003-2007 inventory datasets against observations provided much better agreement for urban aerosol when TNO2010 was used (Chubarova, 2020). This enables us to apply the TNO2010 inventory in this study. The preliminary comparisons with the modern CAMS inventory dataset for 2019 also showed agreement of the urban aerosol estimates. The spatial distribution of monthly mean aerosol gas precursor urban emissions, including SO 2 , NO x , non-methane volatile organic compounds (NMVOCs) as well as direct emissions of BC and PM 10 over the Moscow area for April and May, is shown in Fig. 2. One can see that most of the urban emissions are observed over the center of the Moscow megacity due to the influence of traffic. In April the emissions of SO 2 and NO x are larger due to the heating season.
Aerosol concentrations at the borders of the simulated area were assumed to be close to zero to exclude the influence of regional background aerosol and aerosol gas precursor effects. So, the simulated gas and aerosol concentrations are associated mainly with the urban Moscow emissions, except for small biogenic ones. The time set for aerosol generation was equal to or higher than 31 h, in accordance with the recommendations of the model developers (Bernhard Vogel, private communication, 2020). Thus, we mainly consider the simulation of the anthropogenic components of the surface mass concentrations of PM 10 , BC and aerosol gas precursors as well as the columnar urban component of aerosol optical depth at 550 nm and SSA.  (see Fig. 1). Due to prevailing westerlies and the location of the ZSS site far from local anthropogenic emissions (see Fig. 2), it can be characterized as a background pristine site. This kind of diagnostic provides us with reliable estimates of the urban aerosol effect over the Moscow megacity.

Evaluation of the urban aerosol component
The urban component of columnar aerosol optical depth from both measurements and modeling was estimated as the difference between the data at the MSU MO (marked as Moscow) and at the background ZSS site (marked as Zven): In a similar manner, we estimated the urban components of some other columnar aerosol parameters, such as AEE and the fine and coarse modes of AOD at 500 nm. At ZSS the AERONET measurements have been in operation since 2006; therefore, the parallel measurements between Moscow and Zvenigorod were analyzed for the 2006-2020 period. The time difference between the two instant measurements in these sites is only 3 min.
A similar approach was used for evaluating the urban component of PM 10 mass concentration: In Zvenigorod the PM 10 mass concentration was also measured with the help of the TEOM 1400a instrument by the Mosecomonitoring Agency. Since the data were available only for 2018, the comparisons of the PM 10 urban component were made only for this year. We consider our BC measurements in Moscow to provide the pure BC urb component, whereas black carbon is mainly formed and emitted in the urban environment (see Fig. 2).
The joint use of aerosol measurements and modeling in the atmosphere of the Moscow region provides the more reliable assessment of urban aerosol pollution.
In order to have more accurate evaluation of the urban aerosol component, it was necessary to remove the cases with the influence of smoke aerosol, which has different optical properties (Dubovik et al., 2002;Liu et al., 2018). For the initial information on the location of the biomass burning event, we used the FIRMS (FIRe Monitoring Service) dataset (https://firms.modaps.eosdis.nasa.gov/, last access: 31 January 2022). After identification of biomass burning spots, we applied the backward-trajectory analysis using the READY system (Rolph et al., 2017) with the help of the HYSPLIT visualization model at a height of 0.5-3 km (Stein et al., 2015). We consider the air mass to be affected by biomass burning aerosol if the cases are detected within 50 km of the line of particle motion. If the number of biomass burning spots was smaller than 5, in addition, we analyzed AAE measurements from AERONET at the MSU MO and used the threshold of AAE < 1 to reveal typical aerosol for the Moscow area. During the low-temperature biomass burning process, the AAE values should be higher than 1 because of much intensive absorption at shorter wavelengths (at 440 nm in our case) by organic carbon (Kirchstetter et al., 2004;Sun et al., 2017). More details of this method can be found in Chubarova et al. (2021).
It should be noted that this procedure has been applied only for the comparisons between modeled and measured aerosol parameters, since no fire emissions have been accounted for in the model simulations.

Aerosol characteristics in Moscow according to long-term AERONET measurements
In order to understand whether or not the aerosol features during the intensive experiment were representative for the whole warm period, we analyzed the results of longterm aerosol measurements using the MSU MO AERONET dataset from 2001 to 2020. Figure 3 presents the seasonal variability of the AOD at 500 nm (AOD500), its fine and coarse modes and the AEE according to long-term observations and, in particular, for April and May in 2018 and 2019. One can see a noticeable AOD500 increase during the warm period. The spring maximum of AOD500 is associated with the descent of snow cover and the effects of seasonal agricultural biomass burning in conditions with low precipitation typical for this period . The elevated spring AOD500 values are accompanied by a lower AEE, which also is in agreement with the slightly smaller fraction of the fine-mode AOD500. The summer AOD500 maximum is associated with the active formation of submicron aerosol with a fine-mode AOD500 fraction higher than 80 %. The April-May period of 2018-2019 is characterized by slightly lower AOD500, which is in agreement with a negative AOD500 trend in Moscow in recent years (Chubarova et al., 2016;Zhdanova et al., 2020). The lower fraction of fine-mode aerosol (64 % compared with 71 %-73 %) may also indicate the decrease in the formation of secondary aerosol due to the effective reduction of urban gas precursor emissions in Moscow . However, in general, the aerosol conditions in April and May of 2018-2019 correspond to those during the warm period, with a slightly reduced AOD500 and its fine-mode fraction, which corresponds to the observed trends of purification of the Moscow atmosphere in recent years.
3.2 The main characteristics of aerosol, aerosol gas precursors and their relationship according to the intensive measurement campaigns of 2018-2019 Table 1 shows the statistics of aerosol and gas parameters of the atmosphere during the spring intensive measurement campaigns in 2018 and 2019. The median value of AOD at 500 nm is small (0.12), corresponding to its level in central and northern Europe (Chubarova, 2009;Filonchuk et al., 2019) with the predominance of the fine aerosol mode. The median PM 10 value of 25 µg m −3 is also relatively small and significantly lower than the PM 10 concentrations in Chinese megacities, where average concentrations exceed 100 µg m −3 (Kislov, 2017). However, for some days (16 April 2018, 22, 25 and 27 April 2019), we observed elevated PM 10 levels exceeding the threshold of the daily maximum allowable concentration of 60 µg m −3 adopted as the Russian standard. The median value of SSA = 0.94 is typical for slightly absorbing aerosol, which is in agreement with the rather low BC/PM 10 ratio (4.3 %) and relatively low mean concentrations of BC (1.03 µg m −3 ). However, during some conditions in Moscow we observe an increase in hourly BC up to 8.9 µg m −3 . This corresponds to high BC concentrations of about 5.5-9 µg m −3 , which are observed in Dhanbad, India (Singh et al., 2015), and Guangzhou, China (Wu et al., 2013). Over remote unpolluted areas, however, BC concentrations are about 0.4-0.5 µg m −3 (Bhugwant and Brémaud, 2001;Herich et al., 2011). Due to the predominance of the fine aerosol fraction, the asymmetry factor of the aerosolphase function is relatively small (about 0.63 when considering both the fine and coarse aerosol modes), which also corresponds to relatively high AEE values (Dubovik et al., 2002). The analysis of the aerosol gas precursors has revealed very low concentrations of sulfur dioxide in Moscow, while nitrogen oxides are traditionally high due to heavy traffic in the city and emissions from power plants (Kulbachevsky, 2020). Figure 4 presents the time series of daily mean AOD at 500 nm, PM 10 , BC/PM 10 as well as the concentrations of the main aerosol gas precursors during the intensive campaigns. To characterize meteorological conditions, we also show daily variability of water vapor content W and the IPD in- Table 1. Statistics of hourly mean aerosol characteristics in the total column of the atmosphere and at the surface, including aerosol optical depth (AOD at 500 nm), Ångström extinction exponent, fine-mode fraction at 500 nm, single-scattering albedo (SSA) at 675 nm, factor of asymmetry (ASY) for fine, coarse and total aerosol at 675 nm, water vapor content (W , cm) and surface mass concentrations of PM 10 , BC, different aerosol gas precursors and the BC/PM 10 ratio. dices. There are large variations in both surface and columnar aerosol characteristics of the atmosphere during these periods. In the stable atmosphere with a daily mean IPD of about 2, elevated columnar and surface aerosol loadings are observed (for example, on 12-16 April 2018, 14-16 May 2018 and 20-22 April 2019). However, during the days affected by the advection of biomass burning aerosol (for example, 1 May 2018 or 27 April 2019), there is high aerosol loading even in good air mixing conditions at IPD = 3. These days are also characterized by the elevated NO x concentrations due to the active chemical transformation affected by forest fires (Jin et al., 2021). Note that the high NO x level is observed in spite of low traffic due to weekends (27 April 2019) or holidays (1 May 2018). A correlation matrix has been estimated for evaluating the relationship between different columnar and surface aerosol characteristics, aerosol gas precursors and meteorological parameters (Table 2). There is a high correlation between AOD500 and fine AOD500 mode, which is dominant in central and eastern Europe (Logothetis et al., 2020). The prevailing fine-aerosol-mode fraction is also observed in PM 10 for urban conditions over central and northern Europe (see, for example, Fig. 10 in Wu and Boor, 2021). Relatively high correlation is detected between surface measurements of PM 10 with BC and aerosol gas precursors, except SO 2 , which indicates the importance of these substances for aerosol formation. We also obtained a statistically significant but not very high correlation of columnar AOD500 with surface PM 10 and BC. Fine AOD500 mode has a slightly higher correlation with BC, which could be explained by the finemode BC composition (Bond et al., 2013). The importance of secondary urban aerosol in columnar fine-mode AOD500 (Dubovik et al., 2002) has also been proven by a statistically significant correlation between fine AOD500 mode and aerosol gas precursors (NO x , SO 2 , CH x ); however, the correlation coefficients are not high due to the complexity of the chemical and meteorological processes.
A positive correlation between water vapor content W in the atmospheric column and the aerosol parameters has revealed more favorable processes of aerosol formation in relatively warmer and wetter air masses. In addition, the advection of cold air masses with small W and aerosol loading from northern regions may also be an important cause of this correlation (Szkop et al., 2016). There is also a statistically significant but not very high negative correlation of AOD with surface wind speed due to the ventilation effect in the urban environment, which occurs due to blowing of the urban aerosol out of Moscow. A decrease in AEE and, correspondingly, decrease in the fine AOD fraction with the increase in wind speed may also be associated with less effective fine-mode aerosol generation due to better ventilation conditions and with possibly more effective urban coarse aerosolmode dust resuspension in stronger winds (Hosiokangas et al., 2004;Amato et al., 2009). This is also in accordance with statistically significant correlation between wind speed and aerosol gas precursors. The exception is sulfur dioxide, whose concentrations are extremely small in Moscow (see Table 1), and therefore large errors can be observed when detecting these relations. The pronounced negative correlations with wind speed were found for surface aerosol species, such as PM 10 and, especially, BC. Negative correlation between BC and wind speed was also shown in Chen et al. (2014) and Popovicheva et al. (2020). Note that the observed negative correlation of the BC/PM 10 ratio with wind speed may lead to the decrease in the absorbing properties of the atmosphere in case of high wind speed.
A statistically significant but not high correlation between surface aerosol gas precursors and IPD indices confirms more favorable conditions for intense air mixing, which, as a result, provide a decrease in aerosol gas precursor mass concentration. However, the correlation of columnar AOD and PM 10 with the IPD index, in contrast to wind speed, is not statistically significant, probably due to the considerable contribution of natural aerosol in AOD and PM 10 in Moscow. The closer relationship of wind speed and IPD with BC compared with PM 10 indicates a more important role of local meteorological situations for black carbon, since urban emissions of pollutants are the main source of BC, while PM 10 has a regional aerosol source, which undergoes significant variations (Air quality in Europe, 2020). The smaller PM 10 negative correlation with wind speed could also be explained by the effects of dust resuspension in stronger winds (Hosiokangas et al., 2004;Amato et al., 2009).  A more detailed analysis of the relationship between AOD500 and PM 10 surface mass concentrations, shown in Fig. 5a, demonstrates that, along with the existence of generally not very high correlation (see Table 2), there is a split into two types of dependences at a point of bifurcation of PM 10 ∼ 0.05 mg m −3 . A weaker AOD500 dependence versus PM 10 characterizes the accumulation of PM 10 only in the low layer (due to the local emission sources near the surface) in the absence of the considerable AOD500 increase, with many cases at IPD = 1 relating to the low intensity of particle dispersion. A more pronounced dependence between AOD500 and PM 10 is associated with the influence of air mass advection, when the concentration of surface particles increases simultaneously with AOD500. In this case only a few cases at IPD = 1 are observed (Fig. 5a). The increase in PM 10 is also connected with the significant increase in the fine-mode AOD500 fraction and the total absence of its low values at high PM 10 levels (Fig. 5b). The existence of these two dependences may explain the not very high correlation between AOD500 and PM 10 for the whole dataset.
There are also noticeable variations in the BC/PM 10 ratio, depending on PM 10 and IPD (Fig. 5c). In well-mixed air conditions (IPD = 3), much lower values of the BC/PM 10 ratio are observed: in most cases, they are smaller than 0.01 and decrease with the growth of PM 10 . This corresponds to the situation when there is advection of air outside of Moscow, with a high natural aerosol content but with a relatively low BC content. On average, at IPD = 3, the BC/PM 10 ratio is equal to 0.7 %. At the same time, with the IPD decrease, the BC/PM 10 ratio becomes higher with mean values of 5.5 % and 5.9 %, respectively, at IPD = 2 and IPD = 1. Thus, the use of IPD data may significantly refine the BC/PM 10 level and, as a result, the absorbing properties of the atmosphere. Figure 6 presents the scattering diagrams of BC mass concentration as a function of PM 10 , NO 2 and SO 2 for different IPD regimes obtained according to both measurements and COSMO-ART simulations. Model simulations confirm close relationships of BC with PM 10 and NO 2 . At the same time, the correlation of BC with sulfur dioxide has been revealed only by modeling at relatively high concentrations of SO 2 , which are not observed in Moscow (Kislov, 2017;Kulbachevsky, 2019). This indicates that the data on SO 2 emissions in Moscow according to the TNO inventory were overestimated. The main source of SO 2 emissions is usually the coal fuel at power plants, which is not used in the Moscow region, except for situations of extremely cold winters (Kislov, 2017). In addition, the Euro-5 motor fuel standard, which has been used in Moscow since 2016, provided low-SO 2 emissions in the atmosphere.
Since the measurements of black carbon are very sparse, in some cases it may be necessary to evaluate its concentration according to the available measurements of the gas composition at environmental monitoring stations. According to our measurements, hourly values of BC mass concentration (µg m −3 ) can be evaluated from PM 10 (µg m −3 ) or NO 2 (µg m −3 ) using the following regression equations: where R is the Pearson correlation coefficient. These regression dependences can be used as a first approximation for the estimates of BC concentrations during the warm period with relatively high temperatures and high solar radiation providing favorable conditions for photochemistry, which is important for the NO 2 production. The similar results, obtained by modeling, confirmed the possibility of using these regression dependences. Note that the account of IPD can additionally increase the correlation between BC and PM 10 (R = 0.94 for IPD = 1, R = 0.81 for IPD = 3, compared with R = 0.64 for the whole dataset). A similar but smaller increase in correlation is observed after the IPD account in the relationship between BC and NO 2 (R = 0.74 for IPD = 1, R = 0.85 for IPD = 3 compared with R = 0.7 for the whole dataset).

Relationships between the aerosol single-scattering albedo and BC/PM 10 ratio
We noted earlier that the BC/PM 10 ratio may characterize the absorbing properties of aerosol, since BC is almost the only source of solar radiation absorption in the visible spectral range, and its high concentrations can lead to a decrease in aerosol SSA (Kozlov et al., 2016) and to significant radiative heating of the atmosphere. As a result, we proposed using the BC/PM 10 ratio in a first approximation for estimating the SSA values. This approach might also be useful in different atmospheric tasks, since in the standard AERONET algorithm there is a strong limitation on cloud-free conditions and relatively high AOD for SSA retrievals (Dubovik and King, 2000). However, it is necessary to understand how accurately the BC/PM 10 ratio at the surface captures the conditions of the entire column of the atmosphere. According to our observations, the restriction of air mixing conditions at IPD = 3 was not enough to obtain the relationship between these characteristics. In addition, we applied a limitation to the daytime period (±3 h around the solar noon), when a significant increase in air convection is observed during the warm period. The application of this additional restriction provides the dependence between SSA and BC/PM 10 (Fig. 7), which is close to the results obtained in the previous experiment in Moscow (Chubarova et al., 2013). The dependence was not strong, possibly due to the large uncertainty (about 0.03) of the SSA AERONET retrievals (Dubovik and King, 2000) and the relatively small statistics. Model estimates of SSA dependence on the BC/PM 10 ratio provide a much more significant relationship with the correlation coefficient R = 0.87, but the values themselves are lower and the SSA sensitivity to the BC/PM 10 value is higher. Thus, to better attribute this dependence, further analysis with more statistics is required.

Aerosol urban pollution based on comparisons
between Moscow and background conditions at the ZSS As described in Sect. 2, we estimated the urban aerosol pollution in the Moscow megacity as the difference of aerosol characteristics between MSU MO and the Zvenigorod site (see Eqs. 1 and 2). Figure 8 shows the annual mean total urban component of AOD, the fine and coarse modes of AOD urb , AEE urb , and the AOD urb / AOD ratio for the entire period of parallel AERONET observations in Moscow and Zvenigorod from 2006 to 2020. On average, total AOD urb at 500 nm is about 0.025, with a predominant fine-mode AOD urb = 0.021, which is in agreement with the positive sign of AEE urb . The AOD urb / AOD ratio at 500 nm is about 19 %. No statistically significant difference in the coarse AOD mode is found between Moscow and the clean unpolluted site. The inset in Fig. 8 shows the AOD urb spectral dependence, which is characterized by larger values at  shorter wavelengths corresponding to the fine aerosol mode (Chubarova et al., 2011b). The inverse dependence of AOD in the UV region with smaller AOD urb at 340 nm may be observed due to a slight underestimation of nitrogen dioxide content in the atmospheric column in Moscow, which is used in AOD evaluation. Note that this underestimation was much larger in the previous version 2.0 AERONET dataset (see the discussion in Chubarova et al., 2011b). For the period of the experiment (April-May 2018 and 2019), we compared the observed AOD urb , which was not affected by the smoke air advection from the areas of forest and agricultural fires, with the parallel model AOD urb values. In addition, we removed the cases with cloud amount N > 5. The latter filtering is necessary, since in cloudy conditions there is too active aerosol generation in the COSMO-ART model system (Chubarova, 2020). Figure 9a shows the time series of the measured and modeled AOD urb at a wavelength of 550 nm and aerosol optical thickness AOD550 from the observations in Moscow. The model AOD urb varies mainly in the range of 0.05, reaching in some cases 0.1-0.17. The measured AOD urb varies in the larger range from −0.12 to +0.14. Negative AOD urb may be associated with the influence of the advection of polluted Note: we used a filter of total cloud amount N < 5. We also consider only the cases with the AOD data during daytime for evaluating the PM 10 statistics and no effects from biomass burning aerosol.
air from Moscow, which will be analyzed later. On average, model and measured AOD urb values are in good agreement, comprising 0.015-0.016 (see Table 3) or 15 %-16 % of the total AOD550. Note that these estimates include the cases with urban air advection from Moscow to the Zvenigorod area. Figure 9b presents the time series of the PM 10 mass concentration in Moscow and its modeled and measured urban components. We see significant variations in the measured PM 10urb , which can be negative during the night or early in the morning. On average, the model PM 10urb is lower than the measured PM 10urb values (6.0 and 15.9 µg m −3 , respectively; see Table 3). The higher values of the measured PM 10urb provide a larger PM 10urb / PM 10 ratio of about 70 %, while according to the model estimates, this ratio is much smaller (about 27 %). This may happen due to some underestimation of urban aerosol and gas emissions in Moscow megacity conditions, which should be studied further.
Since BC is the almost purely urban aerosol component in the absence of smoke aerosol advection, the model BC urb / BC ratio comprises more than 93 % of the total BC. We also see higher model BC concentrations compared with the measurements (1.6 and 0.95 µg m −3 , respectively), which may also result in the too low model single-scattering albedo in urban conditions shown in Fig. 7. This may happen due to the overestimation of BC emissions in the TNO2010 inventory dataset.
Urban aerosol may have a relationship with natural aerosol, since they are both determined by the chemical composition of the atmosphere and meteorological conditions. To evaluate their relations, we analyzed the dependences between the urban aerosol component and its total amount. Figure 10a and b present modeled and measured AOD urb and the AOD urb / AOD550 ratio as a function of AOD550 according to the MSU MO measurements. There are large variations in AOD urb obtained from measurements and modeling. According to the model results, there is a slight positive AOD urb in- Figure 8. Annual mean urban components of different aerosol parameters -AOD urb at 550 and 500 nm, fine-mode AOD urb , coarse-mode AOD urb , urban component of AEE urb and AO-Durb / AOD500 ratio with confidence intervals at the 0.05 significance level. The inset shows the mean AOD urb spectral dependence. Comment: we show the AOD at two wavelengths to provide more convenient comparisons with the CIMEL sun-photometer observations (AOD at 500 nm) and model results (AOD at 550 nm). Moscow. 2006Moscow. -2020 crease at AOD550 > 0.2. The absence of the dependence for measured AOD urb versus AOD at 550 nm (see Fig. 10a) can be observed due to a significant contribution of the advection of natural aerosol with high AOD. Figure 10c and d present similar dependencies of PM 10urb and the PM 10urb / PM 10 ratio as a function of the observed PM 10 in Moscow. A positive correlation dependence between PM 10urb and PM 10 is more pronounced, especially that obtained from observations. This can be explained by higher concentrations of aerosol gas precursors at the surface and, as a result, more favorable conditions for the secondary aerosol generation. Note that there is a significant correlation between the concentrations of aerosol gas precursors and PM 10 (see Table 2).
The analysis of the AOD urb / AOD550 and PM 10urb / PM 10 ratios has revealed a tendency to decrease at high aerosol content. This may mean that large aerosol content in Moscow is observed due to advection, while the largest contribution of urban aerosol (up to 50 %-100 %) is detected at a relatively low aerosol level of about AOD550 = 0.1-0.2 or PM 10 < 0.04 mg m −3 .
As mentioned earlier, the increase in aerosol loading over the Moscow suburbs can occur due to the advection of polluted air from Moscow. Therefore, the calculated and measured aerosol urban components were compared in addition for the cases when the influence of the air advection from the Moscow megacity was not observed. To remove such cases from the sample, we applied the HYSPLIT model ensembles of the 24 h forward trajectories (Stein et al., 2015) at 500-1000 m for noon conditions. We consider that, since the Zvenigorod site is located directly to the west of Moscow's city center, the air quality was not affected by Moscow pollution if the trajectories were in the zone from 0 to 180 • . Figure 11 shows the dependence between modeled and measured AOD urb and PM 10urb for all the cases ( Fig. 11a and  b) and for the cases without air advection from Moscow ( Fig. 11c and d). The analysis was made for the same AOD and PM 10 statistics obtained during daytime and in sunny conditions. We see that, after the removal of the Moscow influence, far fewer cases with the measured negative AOD urb values are observed (compare Fig. 11a and c), and the remaining negative AOD urb values do not exceed 0.01, which is the uncertainty of AOD measurements in AERONET (Holben et al., 1998). As a result, after removing the Moscowaffected cases, the mean AOD urb value is equal to 0.019, which is slightly higher compared with AOD urb = 0.016 obtained for all the cases (see Table 3). Similarly, Fig. 11b and d present the relationships between the calculated and measured PM 10urb for all the cases and for the cases without the Moscow influence. Interestingly, during daytime there were no negative PM 10urb values when PM 10 zven > PM 10 Moscow , as shown in Fig. 9b. The observed mean PM 10urb is equal to 15.9 µg m −3 , which is higher than the model PM 10urb estimates (6 µg m −3 ), while for BC we have, by contrast, model overestimation (1.5 against 0.95 µg m −3 ).
Since the emissions and intensity of particle dispersion have a pronounced daily course, the urban aerosol component may also have significant differences. Figure 12 shows the composite diurnal cycles of AOD, PM 10 and BC at different IPDs as well as the primary emissions of black carbon and PM 10 according to the TNO2010 inventory. In general, there are noticeable diurnal changes in modeled and experimental data at the surface layer, which has some specific features depending on IPD. One can see the accumulation of PM 10 and BC at night below the inversion layer in the stable atmosphere, which is characterized by IPD = 1. Note that, during daytime (from 10:00 to 17:00), the conditions with IPD = 1 were never recorded because of warming of the surface and the amplification of convection.
As for the columnar AOD characteristic, there is no evident diurnal cycle of measured AOD urb during daylight hours; however, model AOD urb values demonstrate a small increase at night, especially under conditions with IPD = 1. Figure 12b and c show a noticeable dependence of BC urb and PM 10urb on the IPD index, especially for night and early morning conditions. The elevated values of the surface urban aerosol at night in conditions with IPD = 1 reach 30-40 µg m −3 for PM 10urb and 3-3.5 µg m −3 for BC urb .
The PM 10 and especially BC diurnal cycles are mainly determined by variations of the boundary layer of the atmosphere. In warm periods there is an increase in its height during daytime, which contributes to the processes of dilution and strengthening of convective processes due to the additional heating by solar radiation, which leads to a decrease in the concentration of BC at the surface (Ramachandran and Rajesh, 2007;Kozlov et al., 2011;Chen et al., 2014). There is the absence of a morning BC maximum in Moscow during rush hours, which was observed in many other cities, for example, in Tomsk (Kozlov et al., 2011), in Ahmedabad (Ramachandran andRajesh, 2007) and in Athens (Diapouli et al., 2017). It can be explained by the specific regulation of diesel-heavy trucks, which have permission of entry only at night in Moscow (Popovicheva et al., 2020). In addition, it is necessary to account for a specific location of the MSU MO at a distance from the direct sources of urban emissions (highways).

Discussion
The analysis of the urban aerosol pollution was done for a large agglomeration of the Moscow megacity using the COSMO-Ru-ART mesoscale model system and intensive measurement campaigns in April-May of 2018 and 2019 for a wide range of meteorological and atmospheric conditions. We showed that, on average, the columnar aerosol characteristics during the intensive spring campaigns were similar to those observed during the warm period. However, AODs were slightly smaller compared with corresponding monthly average values over the 2001-2020 period. This is in agreement with the observed negative AOD trend in the Moscow megacity (Chubarova et al., 2016;Zhdanova et al., 2020). A reduction in the fine-mode AOD fraction may be associated with a decrease in the emissions of urban aerosol gas precursors in recent years . A weak aerosol absorption in Moscow with relatively high values of SSA = 0.94 corresponds to relatively small concentrations of black carbon (for urban areas) and its low BC/PM 10 ratio. The BC mass concentration is consistent with the estimated BC values of 1.1 µg m −3 in the GADS database for the Moscow region during the warm period (Koepke et al., 1997). This is twice as high compared with the BC concentrations in clean unpolluted regions (Herich et al., 2011) and more than 5 times smaller than in the polluted areas of India and China (Wu et al., 2013;Singh et al., 2015). We show that Moscow is characterized by reduced concentrations of sulfur dioxide and by increased concentrations of nitrogen oxides due to the emissions from transport and power plants (Kulbachevsky, 2019(Kulbachevsky, , 2020. Mean PM 10 concentrations in Moscow correspond to those in large European cities (about   15-30 µg m −3 ) and are significantly smaller than those in Asian industrial centers (Kislov, 2017).
The analysis of the experimental data during the intensive spring campaigns shows a noticeable day-to-day variability of different gases, aerosol characteristics and IPD indices. In some cases, especially high concentrations are observed at IPD = 1, but in some cases smoke advection, even in conditions with IPD = 3, provides elevated levels of aerosol pollution.
We obtained high correlations between AOD500 and the fine AOD500 mode, which were typical for central and eastern Europe (Logothetis et al., 2020). Fine aerosol-mode fraction also dominates in PM 10 in urban regions of central and northern Europe according to Wu and Boor (2021). We found a relatively high correlation between surface measurements of PM 10 and BC with aerosol gas precursors, except for SO 2 . This may indicate the importance of secondary aerosol generation in the urban atmosphere of Moscow. According to the model simulations, secondary aerosols are about 55 % of the total urban aerosol content, which is generally in accordance with Huang et al. (2014), with the maximum contribution of the organic component (30 %).
The atmospheric water vapor content, which can be used as an indicator of the air mass (Myachkova, 1983), has a relationship with the aerosol parameters in the atmospheric column, which indicates more favorable processes of aerosol formation in relatively warmer and wetter air masses (Chubarova, 2009). Wind speed has a statistically significant correlation with almost all surface and columnar aerosol characteristics as well as aerosol gas precursor species. This could be explained by the effects of ventilation of the urban environment, its better dilution and movement of the urban polluted air to the suburbs. The closer relationship between wind speed and BC compared with PM 10 indicates a more important role of local meteorological conditions for black carbon. This happens since BC is characterized by urban origin, while for PM 10 the regional background aerosol variations may also play an important role. The obtained negative relationship of BC and the BC/PM 10 ratio with wind speed is especially important, since it can serve as an indicator of changes in the absorbing properties of the atmosphere. We also found the influence of IPD on the BC/PM 10 mean ratio, which is 0.7 % at IPD = 3, 5.5 % at IPD = 2 and 5.9 % at IPD = 1. The small BC / PM 10 ratio at IPD = 3 can be explained by the intensive air advection from the clean areas outside Moscow with low BC concentrations and by strong vertical mixing with the upper layers of the atmosphere with smaller BC concentrations. By contrast, at IPD = 1, in conditions of a stable stratified atmosphere and the absence of ventilation, we observe a strong increase in BC/PM 10 . Thus, the use of the IPD index may significantly specify the BC/PM 10 ratio and, as a result, the absorbing properties of the atmosphere.
The relationship between columnar AOD and surface PM 10 concentrations has a point of bifurcation around 0.05 mg m −3 , revealing two types of dependences, which has a physical explanation. The lower dependence characterizes the growth of PM 10 only in the close to surface atmospheric layer with predominant IPD = 1 in the absence of the pronounced AOD increase. The upper dependence is associated with the influence of air mass advection, when the concentration of surface particles increases simultaneously with AOD. In Gubanova et al. (2018), close links between the AOD and PM 2.5 concentrations on monthly scale have also been obtained, especially during the warm period. In addition, we show a significant increase in the fine-mode AOD fraction with PM 10 , and in conditions with PM 10 > 80 µg m −3 , the fine-mode AOD fraction is always higher 70 %.
In the analysis, more attention is paid to BC, since this aerosol component significantly absorbs visible solar irradiance (Jacobson, 2004(Jacobson, , 2006Ramanathan and Carmichael, 2008). Based on the measured data, we have obtained simple regression equations to quantify BC concentration using the observed PM 10 or NO 2 mass concentrations for warm period conditions. These relationships have also been confirmed by the model simulations. In addition, the account of IPD indices in these dependences has led to a stronger connection between these characteristics.
We show that in the well-mixed atmosphere there is an inverse dependence of columnar aerosol SSA on the BC/PM 10 ratio according to both model and measurements. Model estimates provide much higher correlation between these parameters and stronger sensitivity of SSA to the BC / PM 10 ratio. These dependences should be studied further using larger statistics, since they can be important for indirect evaluation of absorbing aerosol properties in cloudy conditions, when the observations of SSA are not available from AERONET.
According to the parallel simultaneous AERONET measurements in Moscow and at the background Zvenigorod site over the 2006-2020 period, the annual mean AOD urb at 500 nm is 0.025, corresponding to 19 % of the total AOD500. We also determined that AOD urb consists of 85 % of the fine-mode fraction. These AOD urb estimates are in agreement with the results of the first years of the urban aerosol studies of Moscow (Chubarova et al., 2011b). Note that this is the only pair of long-term simultaneous AERONET measurements in the world, which are located inside and outside the large megacity and operated by the same type of sun photometer. The estimates of the urban AOD of about 0.02 were also made using a handheld sun photometer in Warsaw (Zavadzka et al., 2013). Warsaw with a population of less than 2 million is much smaller than the Moscow megacity. However, the active use of coal with large SO 2 emissions may provide a significant increase in sulfate aerosol generation and, hence, may lead to the increase in AOD urb . For Moscow conditions the results of MODIS satellite estimates using the MAIAC multi-angle algorithm with accounting for surface reflectivity (Lyapustin et al., 2018) showed the urban AOD550 effect of about 0.01 . In contrast, in Li et al. (2018), over Berlin according to the MODIS MYD04_3 data, a much higher urban AOD (about 0.08) was obtained. These AOD urb values are likely to be overestimated, since Berlin with a population of 3.6 million people and restrictions on emissions should have lower aerosol pollution compared with the Moscow megacity.
We found that the AOD urb / AOD550 and PM 10urb / PM 10 ratios had a tendency to decrease at high aerosol content. This may mean that the large aerosol content in Moscow was observed mainly due to advection, while the urban aerosol contribution of up to 50 %-100 % was important at a relatively low aerosol level of about AOD = 0.1-0.2 or PM 10 < 40 µg m −3 .
The increase in aerosol loading in the background clean conditions in Zvenigorod can also occur due to the advection of polluted air from Moscow. After the removal of the cases affected by Moscow pollution, the amount of negative AOD urb became much smaller, and it did not exceed 0.01, which is the uncertainty of sun-photometer measurements. As a result, the average aerosol urban pollution, according to the observations, was estimated as AOD urb = 0.019, PM 10urb = 15.9 µg m −3 and BC urb = 0.95 µg m −3 . The urban AOD fraction (AOD urb / AOD meas ) was equal to 19 %, which was in agreement with the annual mean over the 2006-2020 period as well as with the estimates in Chubarova et al. (2011). The model PM 10urb provided some underestimation compared with measurements (6 and 15.9 µg m −3 , respectively). By contrast, there was a model overestimation of BC concentration (1.6 µg m −3 compared with 0.95 µg m −3 ), which may be the cause of the too low values of the model single-scattering albedo. This difference of model estimates with the observations occurred likely due to underestimation of primary aerosol emissions and overestimation of BC emissions in the TNO2010 inventory for the Moscow megacity, which should be studied further.
The analysis of the composite diurnal cycles for PM 10urb and BC urb at different IPD at the surface layer have detected their noticeable changes in both model and experimental data. We revealed the significant accumulation of PM 10 and BC at night below the inversion layer at IPD = 1. In these conditions at night the increase in concentration reached 4 times for PM 10urb (30-40 µg m −3 ), and 3 times for BC (up to 3-3.5 µg m −3 ) compared with the conditions at IPD = 3. During the warm period there was a noticeable increase in the height of boundary layer during daytime, which contributed to the processes of dilution and amplification of convection, when the Earth's surface was heated by solar radiation (Ramachandran and Rajesh, 2007;Kozlov et al., 2011;Chen et al., 2014) providing a decrease in surface concentrations of different aerosol and gas species. The observed AOD urb values did not have a clear diurnal cycle; however, according to model calculations, they were slightly higher at night, especially in conditions of the stratified atmosphere at IPD = 1.

Conclusions
We have presented a detailed analysis of surface and columnar aerosol measurements and model simulations in urban and clean background conditions, which allows us to obtain reliable quantitative estimates of the urban component of the aerosol pollution at the surface and in the atmospheric column and to identify the relationships between them in different meteorological conditions.
We found the predominance of the fine AOD500 mode in AOD500 and a statistically significant though not very high correlation between columnar AOD500 and surface PM 10 mass concentrations with the effect of splitting the depen-dence into two different ones. Relatively high correlations between surface measurements of PM 10 and BC were observed with aerosol gas precursors, except for SO 2 .
Both model and experimental datasets have shown a statistically significant linear correlation of BC with NO 2 and PM 10 mass concentrations, which indicated mostly common sources of emissions of these substances.
The annual mean urban component of AOD (AOD urb ) at 500 nm in the Moscow megacity over the 2006-2020 period was about 0.025 (or 19 % of total AOD), with more than 85 % of the fine-mode fraction. According to the measurements, we have obtained the mean estimates of AOD urb = 0.019, PM 10urb = 15.9 µg m −3 and BC = 0.95 µg m −3 after removing the cases affected by smoke aerosol advection. A similar level of AOD urb = 0.015 has been evaluated using model simulations.
There was a pronounced increase in the BC/PM 10 ratio from 0.7 % to 5.9 % with the decrease in the IPD index related to the amplification of the atmospheric stratification. We also found an inverse dependence between the BC/PM 10 ratio and columnar single-scattering albedo (SSA) for the intense air mixing conditions. This dependence together with the obtained negative correlation between wind speed and the BC/PM 10 ratio may serve as an indicator of changes in the absorbing properties of the atmosphere due to meteorological factors.
A pronounced diurnal cycle of PM 10urb and BC urb has been obtained with a strong dependence on IPD indices. At night a significant accumulation of PM 10urb and BC urb below the inversion layer was observed in conditions with IPD = 1, reaching 4 times for PM 10urb and 3 times for BC compared with the conditions at IPD = 3. The observed AOD urb did not have a clear diurnal cycle; however, according to model calculations, it was slightly higher at night, especially in poorly mixed conditions at IPD = 1.
In future work, we plan to use the obtained results for evaluating the radiative effects of the urban aerosol pollution and for identifying its influence on meteorological parameters and weather forecast.
Author contributions. The conceptualization, data analysis, and final text writing were done by NEC. AAK, GSR, BV and HV designed the model experiments and performed the simulations. OBP contributed with the BC dataset, and EEA contributed with data analysis and the design of the figures. NEC prepared the manuscript with contributions from all the co-authors.
Competing interests. The contact author has declared that none of the authors has any competing interests.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Special issue statement. This article is part of the special issue "Pan-Eurasian Experiment (PEEX) -Part II". It is not associated with a conference.