Assessment of the COVID-19 Lockdown Effects on Spectral Aerosol Scattering and Absorption Properties in Athens, Greece

COVID-19 is evolving into one of the worst pandemics in recent history, claiming a death toll of over 1.5 million as of December 2020. In an attempt to limit the expansion of the pandemic in its initial phase, nearly all countries imposed restriction measures, which resulted in an unprecedented reduction of air pollution. This study aims to assess the impact of the lockdown effects due to COVID-19 on in situ measured aerosol properties, namely spectral-scattering (bsca) and absorption (babs) coefficients, black carbon (BC) concentrations, single-scattering albedo (SSA), scattering and absorption Ångström exponents (SAE, AAE) in Athens, Greece. Moreover, a comparison is performed with the regional background site of Finokalia, Crete, for a better assessment of the urban impact on observed differences. The study examines pre-lockdown (1–22 March 2020), lockdown (23 March–3 May 2020) and post-lockdown (4–31 May 2020) periods, while the aerosol properties are also compared with a 3–4 year preceding period (2016/2017–2019). Comparison of meteorological parameters in Athens, between the lockdown period and respective days in previous years, showed only marginal variation, which is not deemed sufficient in order to justify the notable changes in aerosol concentrations and optical properties. The largest reduction during the lockdown period was observed for babs compared to the pre-lockdown (−39%) and to the same period in previous years (−36%). This was intensified during the morning traffic hours (−60%), reflecting the large decrease in vehicular emissions. Furthermore, AAE increased during the lockdown period due to reduced emissions from fossil-fuel combustion, while a smaller (−21%) decrease was observed for bsca along with slight increases (6%) in SAE and SSA values, indicating that scattering aerosol properties were less affected by the decrease in vehicular emissions, as they are more dependent on regional sources and atmospheric processing. Nighttime BC emissions related to residential wood-burning were slightly increased during the lockdown period, with respect to previous-year means. On the contrary, aerosol and pollution changes during the lockdown period at Finokalia were low and highly sensitive to natural sources and processes.


Introduction
The novel coronavirus SARS-COV-2 that was first evidenced in Wuhan, China, has evolved into the worst pandemic of the century so far due to widespread infection from late winter 2019 to this day [1,2]. On 11 March 2020, the World Health Organization (WHO) declared the coronavirus-related COVID-19 disease as a global pandemic [3] and,

Study Area, Instrumentation and Methodology
From 23 March 2020, a total lockdown was implemented in Greece, with restrictions in the mobility of private vehicles and public transportation that lasted to 3 May. Since 4 May, when the unrestricted mobility of citizens within the regional administrative boundaries was again allowed, to the end of May, some of the restriction measures have been progressively relaxed [40]. In the present study, we considered as lockdown period the time from 23 March to 3 May, which includes the period with the strictest restrictions in transportation and movement of private vehicles, given that one week later (11 May), the majority of commercial facilities re-opened and many citizens returned to work [40,57]. In order to assess the impact of the lockdown, the aerosol spectral-scattering and absorption properties were compared against pre-  March 2020) and post-(4-31 May 2020) lockdown periods. In addition, in order to eliminate seasonality effects, the aerosol properties in these three sub-periods of spring 2020 were compared with respective means during the previous years (2017-2019 or 2016-2019 for absorption).
Nephelometer and aethalometer measurements in Athens were performed at the Thissio Air Monitoring supersite (Figure 1), which is located atop a hill in the center of the Athens basin (37.97° N, 23.72° E, 105 m a.s.l.) and is considered representative of urban background conditions in central Athens [58]. The site is surrounded by the Acropolis and Pnyx hills and a low-density residential area in its south and west, while the nearby traffic emissions are rather limited since the city center and main avenues are at a distance of more than 500 m. Spectral-scattering aerosol coefficient (bsca(λ)) values (between 7° and 170°; in Mm −1 ) were measured at Thissio by means of a 3-wavelength (450, 550, 700 nm) TSI 3563 Integrating nephelometer (1 min resolution), covering the period 1 March-31 May 2020, as well as the same time frames in 2017-2019. The instrument operates at a relative humidity (RH) below 50% using a processor-controlled automatic dryer for preventing aerosol hy- Spectral-scattering aerosol coefficient (b sca(λ) ) values (between 7 • and 170 • ; in Mm −1 ) were measured at Thissio by means of a 3-wavelength (450, 550, 700 nm) TSI 3563 Integrating nephelometer (1 min resolution), covering the period 1 March-31 May 2020, as well as the same time frames in 2017-2019. The instrument operates at a relative humidity (RH) below 50% using a processor-controlled automatic dryer for preventing aerosol hygroscopicity effects that may increase particle scattering and modify the spectral-scattering dependence [59]. The nephelometer measurements were corrected for angular truncation errors following Anderson and Ogren [60], while the instrument was regularly calibrated using CO 2 as high span gas and air as low span gas [56]. The overall uncertainty in b sca(λ) is considered to be about 7% [61,62]. It should be noted that during all study periods, the b sca(λ) values were above 1 Mm −1 at all wavelengths, a threshold below which the instrumental noise becomes high [63]. Spectral absorption (b abs(λ) ) measurements at Thissio were performed using an AE-33 aethalometer at 7 wavelengths (370, 470, 520, 590, 660, 880 and 950 nm) and at 1 min resolution covering the spring seasons of 2016-2020. The b abs(λ) values were computed from the derived e-BC measurements using the mass absorption cross-section recommended by the manufacturer at each wavelength [64]. Furthermore, the source-specific components of BC, related to fossil fuel (BC ff ) and wood-burning (BC wb ), were computed using the biomass burning fraction (BB%), automatically provided by AE-33 that internally applies the aethalometer model [64,65].
Aethalometer (AE-33) measurements of BC and source-apportioned fractions (BC ff and BC wb ) were also measured at the regional background station in Finokalia, Crete ( Figure 1) during March-May 2020, as well as in the previous years (2016-2019), following the same procedure described above for Athens. Furthermore, during the same period, spectral b sca values were monitored at Finokalia with an integrating nephelometer (Aurora-3000 Ecotech, Melbourne, Australia) at three wavelengths (450, 525 and 635 nm). Calibrations were performed regularly using CO 2 as high span gas and zero particle air as low span gas. The measurements were also corrected for angular non-idealities [66]. During periods that Aurora-3000 was not operating, a similar single wavelength (525 nm) Integrating nephelometer (Aurora-1000 Ecotech, Melbourne, Australia) was used under the same operation protocol, however, providing only total scattering measurements. A few studies in the past have examined the near-surface aerosol scattering and absorption properties at this remote coastal site [67][68][69].
Using the spectral b sca and b abs measurements, we computed intensive aerosol optical properties related to particle size and shape of aerosols (scattering Ångström exponent-SAE), to aerosol chemical composition and sources (absorption Ångström exponent-AAE) and to the relative contribution of scattering and absorption to total aerosol extinction (single scattering albedo-SSA). SAE was computed in the spectral band of the nephelometer (450-700 nm), while AAE at 470-950 nm. The spectral-scattering at the 450-700 nm band was initially fitted by a 2nd order polynomial curve in log-log coordinates, i.e., lnb sca(λ) = A 2 *(lnλ) 2 + A 1 *(lnλ) + A o , since it was found to exhibit a slight curvature rather than a first-order linear wavelength dependence. The accuracy of the polynomial fitting in the 450-700 nm band was verified by the strong relationship between A 2 -A 1 and SAE (slope: 1.03; R 2 = 0.99; [70]). The second-order polynomial was used to extrapolate the nephelometer wavelengths to the aethalometer spectrum 370-950 nm [71]. Therefore, the spectral b sca values were computed at the seven AE-33 wavelengths, and by combining the spectral-scattering and absorption values, the SSA was also estimated at 7 wavelengths. Errors in these computations can be attributed to measurement uncertainties in b sca and b abs , as well as to uncertainties in the extrapolated b sca for wavelengths shorter than 450 and longer than 700 nm. These intensive aerosol properties were analyzed only in Athens, on an hourly basis.
In addition, hourly meteorological data (ambient temperature, RH, pressure, solar radiation, precipitation, wind speed and direction) were obtained during the examined spring periods (2016-2020) from the meteorological stations located nearby the sampling sites. Analysis of aerosol properties at both locations was performed for the same sub-periods, on an hourly basis, while in all cases, the statistical significance of the differences in aerosol properties between the examined periods was checked with t-tests at the 95% confidence level (C.L.).

Aerosol and Pollutant Changes at the Remote Background Site (Finokalia)
Previous works conducted in Athens have shown that a significant part of the aerosol load is due to regional sources [50][51][52][53]. In this respect, before focusing on the Athens urban environment, we examined the variability of b sca and BC components at the regional background site of Finokalia, Crete ( spring periods of 2016-2019, mostly due to strong influence from Saharan dust storms, which are very frequent in this time of the year and are often characterized by extreme intensity, as it was, for example, observed in March 2018 [72]. Therefore, the large b sca peaks on certain days due to natural causes could possibly mask the effects of changes in anthropogenic emissions. However, during the lockdown period, it can be seen in the graph that the b sca time-series was mostly within the shaded area corresponding to the standard deviation of the 2016-2019 mean. On the contrary, the BC levels did not significantly change during the lockdown period with respect to the 2016-2019 mean, showing an increase by only 5% (0.20 to 0.21 µg m −3 ; statistically non-significant), while the BC wb levels remained low (<0.15 µg m −3 ) throughout the examined period ( Figure 2d). However, on days with intense dust storms like on 22 March 2018, BC presented peak levels, which are likely attributed to dust absorption, along with stagnant conditions of enhanced pollution since both BC ff and BC wb exhibited peaks [73]. This is also seen for few days in 2020, despite the low BC levels (<0.6 µg m −3 ), as, for example, during 14-18 May 2020 when a Saharan dust episode affected Greece along with an unprecedented for that season heatwave [40], Mean BC concentrations exhibited a declining trend from pre-lockdown-to-lockdown (−11%) and from lockdown to post-lockdown (−2%). Additional analysis also showed low variability in gaseous pollutant concentrations during the lockdown period compared to pre-lockdown (−10% for NO 2 , −1% for CO and +4% for O 3 ). Overall, the remote site of Finokalia received negligible influence by the COVID-19 lockdown to aerosol scattering and BC absorption, which could indicate limited variability in the regional aerosol sources for the lockdown period compared to the 2016-2019 mean.

Meteorological Conditions in Athens during the Spring Season
Changes in meteorological parameters and boundary-layer dynamics may significantly affect the amount, types and properties of aerosols and mask the effect of changes

Meteorological Conditions in Athens during the Spring Season
Changes in meteorological parameters and boundary-layer dynamics may significantly affect the amount, types and properties of aerosols and mask the effect of changes in emissions [74,75]. Figure S1 shows the variability in the meteorological parameters at Thissio in spring 2020 (pre-, post-and lockdown periods) compared with those during the previous years.
The average air temperature during the lockdown period was slightly lower with respect to the 2016-2019 mean (−9.7%; statistically significant at 95% confidence level). RH showed larger variability with a slightly higher mean value in 2020 (1.5%), while the solar irradiance generally followed the 2016-2019 mean ( Figure S1). The average wind speed (WS) in the lockdown period was 3.0 ms −1 , slightly higher (7.2%; statistically significant) compared to the 2016-2019 mean (2.9 m s −1 ). More details about the variability and relative changes in meteorological parameters between the lockdown period and previous years can be found in Grivas et al. [40].
The wind-rose patterns for b sca and b abs were also used to compare the wind-vector dependence of aerosol properties during the lockdown and similar periods in previous years ( Figure S2). The results showed that the highest b sca,550 and b abs,520 levels during the lockdown were associated with low WS < 3 ms −1 , similarly to previous works [56,76]. A mostly similar pattern was observed for b abs,520 in the previous years, while enhanced values from southern directions and for WS of 4-6 ms −1 were observed for b sca,550 due to intense dust intrusions from the south, such as that on 25-27 March 2018 [77]. Therefore, the distributions of b sca and b abs values in the wind roses did not significantly change, apart from the lower levels due to lower emissions during the lockdown period.

Changes in Aerosol Scattering and Absorption Properties in March-May
A comparative analysis of the hourly time-series of aerosol scattering (b sca,550 ) and absorption (b abs,520 ) coefficients between March-May 2020 and the years 2016/2017-2019 in Athens is shown in Figure 3. The mean b sca,550 and b abs,520 values are given for each period, and the results show a notable decrease in b abs,520 during the lockdown period (15.3 ± 13.8 Mm −1 ) compared to pre-lockdown (25.2 ± 26.9 Mm −1 ), corresponding to an average reduction of −39% (statistically significant at 95% C.L.). In addition, the b abs,520 (15.3 Mm −1 ) during the lockdown period, was notably lower (29-54%) than the respective periods in all the previous years (22-28 Mm −1 ). On the contrary, the reduction in b sca,550 during the lockdown was much lower (−15%; statistically significant) compared to prelockdown and with respect to previous years (−12% to −25%). It should be noted that the decrease in b sca,550 during lockdown with respect to pre-lockdown is similar to the respective decrease (−18%) found for PM 2 . 5 at the same site [40]. On the other hand, significant variability is seen in b sca,550 and b abs,520 between the spring seasons, especially for the scattering coefficient for the period 1-22 March. This is likely attributed to higher RWB emissions for heating during the beginning of March [54,56] and to more unstable meteorological conditions. During May, the atmospheric conditions in Athens are more stable, with lower levels of b sca,550 and b abs,520 ( Figure 3) compared with the previous months. The progressive decrease in b abs,520 from March to May during 2016-2019, is attributed to the declining trend of the BC emissions ( Figure S3), as well as to better dilution processes due to progressively increasing MLH. A recent study has shown that the traffic-related pollutants in central Athens presented an important decrease during the lockdown period (32% for NO2, 35% for CO and 33% for BCff) compared with the pre-lockdown period [40]. These are comparable to the presently reported for babs,520 (−39%). The CO emissions during the non-heating period in Athens originate from the traffic sector and mainly from private vehicles [85,86]. On the other hand, NO2 is also a precursor of nitrate (NO3 − ) aerosols, which are considered as strongly scattering particles [55,87], and the reduction in their concentrations during the lockdown period is partly reflected in the lower bsca,550 values, along with primary organics and sulfates. Therefore, the reduction in bsca,550 during the lockdown period, should also be mainly attributed to lower primary emissions of aerosols related to the traffic sector.
The PM and gaseous pollutant concentrations were found to increase in the postlockdown period due to the re-opening of the economy and the progressive escalation in traffic and human activities [40]. On the contrary, bsca,550 and babs,520 presented only marginal differences (±4-5%) compared to lockdown, as they are significantly affected by the seasonality, regional background conditions and long-range transport of aerosols, mostly in the case of bsca,520. Therefore, in all the examined years, the bsca,550 on May 4-31 was lower than the preceding period (23 March-3 May), indicating a declining seasonal trend in scattering aerosols during the spring season in Athens. However, the mean decrease ( Figure  3a) compared to the previous period was less pronounced in 2020 (2.2 Mm −1 against 9.4 Mm −1 during 2017-2019 on average). A similar feature can be observed for babs,520 with significantly lower values by the end of spring during 2016-2019, as opposed to the slight increase in 2020 due to the notable increase in traffic emissions (Figure 3b). A photochemical pollution event and a concurrent Saharan dust transport episode during 14-18 May 2020, which increased the NO2 production and dust levels in Athens [40], slightly affected the bsca,550 values ( Figure 3a). The effect of this event was also detected at Finokalia ( Figure  2), indicating its regional character. The bsca,550 and babs,520 levels along with the intensive The lower reductions in b sca,550 than in b abs,520 during the lockdown period can be ascribed to the wider chemical composition range of scattering aerosols (including, e.g., organics, nitrate and sulfate) in Athens. Secondary and regionally transported aerosols, which are mostly of scattering nature [78,79], comprise large fractions of PM 2 . 5 concentrations, reaching up to 50% during the summer period [50,52,53]. Continuous monitoring of PM 2 . 5 samples at Thissio and analysis of chemical composition using ion chromatography and EC/OC thermal/optical analysis [53,76] revealed that during April 2020 (lockdown period), the mean concentrations of OC and SO 4 2were 3.4 µg m −3 and 2.7 µg m −3 , respectively. Although OC levels were very close to the long-term (2014-2019) monthly mean of 3.3 µg m −3 , sulfate was lower than the 2014-2019 mean of 3.7 µg m −3 (unpublished results). This decrease would account for~3 Mm −1 according to the parametrization of Kalivitis et al. [69], thus explaining a large part of the reduction in b sca . Although the inorganic species and a large fraction of organic aerosols in urban background Athens are highly water-soluble [53,80,81], the scattering measurements were regulated for RH < 50% levels and, therefore, the water uptake cannot inflate the b sca [59]. Similarly, at the regional background Finokalia site, the regional scattering aerosols did not present important changes, on average, during the lockdown period, and their presence was attributed to natural rather than anthropogenic sources. This justifies the larger reduction in absorption than in scattering in Athens, a fact also supported by a~50% decrease in mean elemental carbon (EC) concentrations during April 2020 (0.61 µg m −3 ) with respect to 2014-2019 mean of 1.2 µg m −3 (unpublished results). In Athens, the principal source of BC (~80%) in April is from fossil-fuel combustion emissions (mainly from the traffic sector) that significantly impact the aerosol absorption [54], and therefore, the drastic restrictions in the road transport sector had a large impact on lowering b abs . The vast majority (~90%) of registered private cars in the Greater Athens Area (around 2.7 million) are gasoline-powered [82]. There are also 0.3 million trucks and buses (mostly diesel-powered) and 0.7 million 2-wheelers [83]. Based on data for traffic volumes from the Greek Ministry of Infrastructure and Transport and Attiki Odos (the intra-city major tollway), the mean weekly vehicular traffic was reduced by about 40-70% during the lockdown period with respect to the same period in the previous two years [40,84].
A recent study has shown that the traffic-related pollutants in central Athens presented an important decrease during the lockdown period (32% for NO 2 , 35% for CO and 33% for BC ff ) compared with the pre-lockdown period [40]. These are comparable to the presently reported for b abs,520 (−39%). The CO emissions during the non-heating period in Athens originate from the traffic sector and mainly from private vehicles [85,86]. On the other hand, NO 2 is also a precursor of nitrate (NO 3 − ) aerosols, which are considered as strongly scattering particles [55,87], and the reduction in their concentrations during the lockdown period is partly reflected in the lower b sca,550 values, along with primary organics and sulfates. Therefore, the reduction in b sca,550 during the lockdown period, should also be mainly attributed to lower primary emissions of aerosols related to the traffic sector.
The PM and gaseous pollutant concentrations were found to increase in the postlockdown period due to the re-opening of the economy and the progressive escalation in traffic and human activities [40]. On the contrary, b sca,550 and b abs,520 presented only marginal differences (±4-5%) compared to lockdown, as they are significantly affected by the seasonality, regional background conditions and long-range transport of aerosols, mostly in the case of b sca,520 . Therefore, in all the examined years, the b sca,550 on May 4-31 was lower than the preceding period (23 March-3 May), indicating a declining seasonal trend in scattering aerosols during the spring season in Athens. However, the mean decrease ( Figure 3a) compared to the previous period was less pronounced in 2020 (2.2 Mm −1 against 9.4 Mm −1 during 2017-2019 on average). A similar feature can be observed for b abs,520 with significantly lower values by the end of spring during 2016-2019, as opposed to the slight increase in 2020 due to the notable increase in traffic emissions (Figure 3b). A photochemical pollution event and a concurrent Saharan dust transport episode during 14-18 May 2020, which increased the NO 2 production and dust levels in Athens [40], slightly affected the b sca,550 values ( Figure 3a). The effect of this event was also detected at Finokalia (Figure 2), indicating its regional character. The b sca,550 and b abs,520 levels along with the intensive aerosol properties (SAE, AAE, SSA) in the examined periods, the differences between lockdown and periods before and after it, as well as with respect to previous years, are summarized in Suppl. Tables S1 and S2. Figure 4 shows the correlations between b sca,550 and b abs,520 values for the pre-lockdown (a) and lockdown (b) periods in 2020, compared with the same periods in 2017-2019. In general, the graphs present considerable scatter (R 2 = 0.36-0.47). The correlations are significantly lower than those observed at other sites in the Mediterranean, where they are generally higher than R 2 = 0.7 [88][89][90][91]. The larger scatter in Athens during the spring season is likely attributed to the variety of aerosol sources, both natural and anthropogenic, and to the effect of transported dust that disproportionately increases the b sca compared to the b abs [62,92]. In addition, there were cases of enhanced b abs,520 for very low b sca,550 (<20 Mm −1 ) values during 1-22 March of 2017-2019 (Figure 4a). On the contrary, the appearance of scatter plots notably changed during the lockdown period, with a much stronger correlation (R 2 = 0.70) between b sca,550 and b abs,520 values, which reflects a more homogeneous atmospheric composition within the Athens basin and a considerable decrease in absorption peaks. Simultaneously, at Finokalia, b sca and b abs were highly correlated (R 2 = 0.80). However, in the case of the regional background site, similar correlations (R 2 = 0.69-0.80) were also observed for the pre-lockdown period and for previous years, indicating the homogeneity of regionally transported aerosol types during the spring season.
correlation (R 2 = 0.70) between bsca,550 and babs,520 values, which reflects a more homogeneous atmospheric composition within the Athens basin and a considerable decrease in absorption peaks. Simultaneously, at Finokalia, bsca and babs were highly correlated (R 2 = 0.80). However, in the case of the regional background site, similar correlations (R 2 = 0.69-0.80) were also observed for the pre-lockdown period and for previous years, indicating the homogeneity of regionally transported aerosol types during the spring season.   [93]. Recently, Kaskaoutis et al. [70] reported that the BC-dominated aerosol type that represents the traffic conditions in central Athens dominates in spring with a fraction of around 90%.   [93]. Recently, Kaskaoutis et al. [70] reported that the BC-dominated aerosol type that represents the traffic conditions in central Athens dominates in spring with a fraction of around 90%.

Changes in Intensive Aerosol Properties
The changes in SAE values (Figure 5a) between the different spring seasons (2017-2020) are generally larger than the slight increase observed during the lockdown period (+11.2% compared to pre-lockdown and +6.1% with respect to previous years; Tables S1 and S2). The slightly larger SAE (2.08 ± 0.26) during the lockdown is likely attributed to a smaller impact from transported dust (note the small SAE values during March 2018 due to several dust storms) and probably to less road dust resuspension due to restrictions in traffic. On the other hand, AAE presents a rather smooth, progressively declining trend from the beginning of March to the end of May (Figure 5b), which is attributed to the decreasing heating demand and the reduction in emissions from residential biomass burning [54]. In 2020, the large decrease in traffic-related absorbing particles, as a consequence of the drastic restrictions, led to higher AAE values (6.3%; statistically significant) during the lockdown period compared to the 2016-2019 mean (Table S1), taking in mind the low AAE values (1.0-1.1) of particles from vehicle exhaust emissions [94][95][96]. Furthermore, a slight intensification of wood-burning emissions was observed on certain nights during the lockdown period that may have also contributed to the increase in AAE, as will be analyzed in the next section. Finally, SSA 520 appeared relatively stable during the spring seasons (Figure 5c). During the lockdown period, SSA 520 slightly increased by 5.9% with respect to the 2017-2019 means and by 9.1% compared to pre-lockdown (Tables S1 and S2). However, in the latter case, the decrease in wood combustion between the pre-lockdown and lockdown periods may affect the results. Overall, the increase in SSA during the lockdown period reflects the larger reduction in levels of highly absorbing BC aerosols from vehicular emissions rather than of scattering-type aerosols [97,98]. The changes in SAE values (Figure 5a) between the different spring seasons (2017-2020) are generally larger than the slight increase observed during the lockdown period (+11.2% compared to pre-lockdown and +6.1% with respect to previous years; Tables S1

Changes in Diurnal Aerosol Patterns during Lockdown, Pre-and Post-Lockdown Periods
In order to evaluate the changes in extensive and intensive aerosol properties and their dependence on the lockdown restrictions, the mean diurnal patterns in the examined periods of 2020 are compared against those of the 3-4 previous years. Figure 6 shows the mean diurnal cycles of b sca,550 and b abs,520 in 2020 for pre-, post-and lockdown periods compared with the respective means in 2016/2017-2019, while similar graphs were produced for SAE 450-700 , AAE 470-950 and SSA 520 (Figure 7), as well as for the BC components ( Figure S3). their dependence on the lockdown restrictions, the mean diurnal patterns in the examined periods of 2020 are compared against those of the 3-4 previous years. Figure 6 shows the mean diurnal cycles of bsca,550 and babs,520 in 2020 for pre-, post-and lockdown periods compared with the respective means in 2016/2017-2019, while similar graphs were produced for SAE450-700, AAE470-950 and SSA520 (Figure 7), as well as for the BC components ( Figure  S3).  The most regional character of scattering aerosols in Athens is supported by the rather weak diurnal patterns of b sca,550 for all periods, since morning traffic peaks are nearly absent, while the changes between pre-and post-lockdown periods were small (mostly within ±10%) (Figure 6a,c). Lower vehicular emissions and road dust resuspension in morning hours (08:00-09:00 LST) [56] did not influence b sca,550 much (Figure 6b), indicating a weak impact of direct traffic emissions on aerosol scattering.
Although b abs,520 diurnal cycles in 2020 closely follow the 2016-2019 means in preand post-lockdown periods (Figure 6d,f), the traffic restrictions during lockdown lead to significant differences throughout the day, reaching −60% in the morning traffic hours (Figure 6e). In all examined periods, the morning traffic peak in absorption clearly prevails, while during the lockdown, the traffic peak appears largely reduced and an hour earlier (07:00 LST) (Figure 6e), probably due to the absence of go-to-work traffic congestion. Similarly, a flat diurnal cycle of BC was noticed in Suzhou, Yangtze River Delta, China during the lockdown period, indicating the remarkable decrease in traffic and the disappearance of the morning and evening traffic-related peaks [41]. During the lockdown night hours, b abs,520 displayed much lower decrease fractions (about 10% with respect to previous years), as the decrease in nighttime traffic emissions is somewhat counterbalanced by a slight intensification of wood combustion. This is also supported by the observed BC ff and BC wb diurnal patterns ( Figure S3). During the post-lockdown period, b sca,550 levels were mostly similar to slightly lower than during lockdown, while b abs,520 was again increased in the morning hours, signifying the free circulation of vehicles. The decrease in the b sca,550 and b abs,520 at noon and early afternoon hours, is attributed to boundary-layer dynamics related to a deeper MLH [99].
The average diurnal cycles of SAE 450-700 do not exhibit pronounced variability during the examined periods. In 2020, higher SAE values can be seen throughout the day (but not more than 15%) with respect to the 2017-2019 means (Figure 7a-c. Apart from a marginal SAE decrease during daytime in pre-lockdown, these mostly constant diurnal patterns imply small changes in aerosol size distribution or in the fine-to-coarse mode ratio, indicating a predominance of fine anthropogenic aerosols throughout the spring season [56,100,101]. Furthermore, the SAE changes during the lockdown period remain nearly constant (4-12%) throughout the day, indicating a limited effect from changes in the traffic sector. It is also characteristic that the highest (%) increase in AAE during the lockdown period is observed during the night hours due to the combined effects of the decrease in evening traffic and the increased RWB emissions from fireplaces on certain days of the lockdown period. On the other hand, secondary organic aerosol (SOA) may slightly contribute to the increase of AAE around midday in all examined periods, following the minimum AEs observed during the traffic rush hours (around 09:00 LST) (Figure 7d-f. The OC/EC ratio is commonly used to empirically discern between primary and secondary OC, estimating the primary and secondary organic carbon fractions (POC and SOC, respectively), using the EC tracer method in Athens [50,70,102]. The mean OC/EC ratio in PM2.5 was 6.2 in April 2020, compared to 2.9, averaged for April during 2014-2019 [53,103]. The large increase of the OC/EC ratio can be attributed to the smaller impact of primary combustion emissions and an enhanced presence of secondary organic aerosol [99,104]. Similarly, the increased organic aerosol fraction has been associated with larger values of AAE and SAE in Rome [105]. Chatterjee et al. [106] reported higher SOA formation during the lockdown in the eastern Indian Himalayas due to reduced NO and high concentrations of O3 and biogenic VOCs, implying It is also characteristic that the highest (%) increase in AAE during the lockdown period is observed during the night hours due to the combined effects of the decrease in evening traffic and the increased RWB emissions from fireplaces on certain days of the lockdown period. On the other hand, secondary organic aerosol (SOA) may slightly contribute to the increase of AAE around midday in all examined periods, following the minimum AEs observed during the traffic rush hours (around 09:00 LST) (Figure 7d-f. The OC/EC ratio is commonly used to empirically discern between primary and secondary OC, estimating the primary and secondary organic carbon fractions (POC and SOC, respectively), using the EC tracer method in Athens [50,70,102]. The mean OC/EC ratio in PM 2 . 5 was 6.2 in April 2020, compared to 2.9, averaged for April during 2014-2019 [53,103]. The large increase of the OC/EC ratio can be attributed to the smaller impact of primary combustion emissions and an enhanced presence of secondary organic aerosol [99,104]. Similarly, the increased organic aerosol fraction has been associated with larger values of AAE and SAE in Rome [105]. Chatterjee et al. [106] reported higher SOA formation during the lockdown in the eastern Indian Himalayas due to reduced NO and high concentrations of O 3 and biogenic VOCs, implying that the site characteristics and the relative fraction between urban/anthropogenic and biogenic VOCs are especially important for the SOA formation processes during the lockdown period.
The steep decrease in SSA values during the morning rush traffic hours in all the examined periods (Figure 7g-i reflects the strong impact of highly absorbing BC aerosols from vehicular combustion [90,107,108]. Following the morning minima, SSA displays an increasing tendency and a large plateau of higher values around noon to afternoon, which is likely ascribed to SOA formation that may increase SSA, as these particles are more efficient scatterers [56,98]. This large plateau continues to the evening/night hours in the post-lockdown period, while conversely, the enhanced RWB emissions at the beginning of March (pre-lockdown period) result in a decrease of SSA during the late evening hours due to absorbing OC from wood combustion [70,108,109]. Previous research in Athens has shown that SSA decreased with increasing BC/PM 1 ratio, while it remained practically unaffected by changes in the organic component [70]. This can explain the decrease in SSA during the morning traffic hours in all periods and the mostly neutral diurnal cycle in the rest of the day. During the lockdown period, SSA was always larger than the 2017-2019 mean, with the highest difference detected in the morning traffic hours (14-16%) due to the higher relative decrease in absorption than scattering (Figure 7h).
Supplemental Figure S3 presents the diurnal variations of the BC, BC ff and BC wb concentrations for the examined periods in 2020 against the mean diurnal patterns of the preceding years. The BC ff component showed a considerable decrease during the lockdown period (−25% to −40%), which maximized (−65%) during the morning traffic hours (08:00-9:00 LST) with respect to the 2016-2019 mean. On the contrary, the BC wb presented a continuous decreasing trend from March to May, while its levels in 2020 were above the 2016-2019 mean in the evening and night hours during the pre-lockdown period, while in lockdown and post-lockdown, BC wb levels were mostly within the range of the 2016-2019 means. The enhanced BC wb concentrations on certain days during the lockdown period (mainly at night hours) are attributed to escalated emissions from the residential sector due to higher demand for heating (also using fireplaces) from people staying or even working indoors [25].

Modifications in Spectral Aerosol Scattering and Absorption
This section examines the changes in spectral aerosol scattering, absorption coefficients and SSA in order to assess the lockdown effect through comparisons with the same period in previous years. The spectral aerosol properties are initially examined during the morning (06:00-10:00 LST) and evening/night (20:00-02:00 LST) hours (Figure 8), in order to quantify the effects of different emission sources, i.e., vehicular emissions during morning rush hours, traffic and biomass burning during the night [80].
The results show a very similar mean spectral dependence of scattering for morning and night hours during 2017-2019, with SAE means of 1.95 and 1.99, respectively. During the lockdown, the spectral b sca values were reduced, more during the morning hours, but at the same rate across the spectrum, as the Ångström power-law fit (b sca,λ = βλ -SAE ) estimated SAE values of 2.14 and 2.16 (Figure 8a). The slightly larger SAE during the lockdown implies a higher relative decrease in spectral b sca at longer wavelengths, which is likely attributed to less coarse material such as road dust due to restrictions in traffic [70]. More pronounced changes were observed in the spectral b abs between morning and night hours for the lockdown and 2016-2019 periods (Figure 8b). Spectral b abs maximizes during the morning hours, with a fitted AAE of 1.18 for the 2016-2019 mean corresponding to fossil-fuel combustion sources, while the mean spectral b abs during night exhibits a higher AAE = 1.33 due to enhanced emissions from RWB [99]. The modifications in the spectral b abs during the lockdown period are mostly emphasized by the higher decrease during the morning hours, but with a slightly higher spectral dependence (AAE = 1.25). For the night hours during the lockdown, the spectral b abs appears enhanced at shorter wavelengths (b abs,370 higher than the morning respective value), possibly indicating a stronger effect by BrC absorption (AAE = 1.52) from RWB emissions [110][111][112].  SSA is an important climatic factor, and its spectral dependence has also been widely used as an indicator of aerosol types [94,[113][114][115]. Decrease of SSA by wavelength is characteristic of urban-industrial aerosols with high BC content from fossil-fuel combustion [107,116], as these aerosols exhibit stronger wavelength dependence of scattering than absorption. On the contrary, external or internal mixing with other mostly scattering urban aerosols (organics, sulfate, nitrate) lowers the decreasing spectral SSA rate, while a high presence of BrC from biomass burning may result in dSSA/dλ ratios close to 0, or even positive [109,117,118]. The current analysis shows a rather subtle modification in the spectral SSA for day and night hours during the lockdown period, compared to the same periods in preceding years (Figure 8c). The spectral SSA is lower during lockdown morning compared to nighttime due to the strong effect that the traffic sector still has in spite of the restrictions ( Figure S3). It also exhibits a higher SSA Ångström exponent (SSAAE = 0.292) compared to night (SSAAE = 0.205). The comparison with respect to previous years reveals that the SSA during lockdown daytime is higher by ~9% than the mean due to reasons discussed above (Figures 5c and 7h), but it exhibits a similar wavelength dependence, indicating negligible changes in aerosol-type composition. The SSAAE values during night-time are slightly lower due to the higher presence of BrC, while the enhanced RWB during the lockdown nights (mean BB% = 36%) could explain the lower SSAAE of 0.205 [118][119][120] since light-absorbing BrC displays lower decreasing rates of SSA at longer wavelengths compared to BC [73,121,122]. In supporting this statement, Katsanos et al. SSA is an important climatic factor, and its spectral dependence has also been widely used as an indicator of aerosol types [94,[113][114][115]. Decrease of SSA by wavelength is characteristic of urban-industrial aerosols with high BC content from fossil-fuel combustion [107,116], as these aerosols exhibit stronger wavelength dependence of scattering than absorption. On the contrary, external or internal mixing with other mostly scattering urban aerosols (organics, sulfate, nitrate) lowers the decreasing spectral SSA rate, while a high presence of BrC from biomass burning may result in dSSA/dλ ratios close to 0, or even positive [109,117,118]. The current analysis shows a rather subtle modification in the spectral SSA for day and night hours during the lockdown period, compared to the same periods in preceding years (Figure 8c). The spectral SSA is lower during lockdown morning compared to nighttime due to the strong effect that the traffic sector still has in spite of the restrictions ( Figure S3). It also exhibits a higher SSA Ångström exponent (SSAAE = 0.292) compared to night (SSAAE = 0.205). The comparison with respect to previous years reveals that the SSA during lockdown daytime is higher by~9% than the mean due to reasons discussed above (Figures 5c and 7h), but it exhibits a similar wavelength dependence, indicating negligible changes in aerosol-type composition. The SSAAE values during night-time are slightly lower due to the higher presence of BrC, while the enhanced RWB during the lockdown nights (mean BB% = 36%) could explain the lower SSAAE of 0.205 [118][119][120] since light-absorbing BrC displays lower decreasing rates of SSA at longer wavelengths compared to BC [73,121,122]. In supporting this statement, Katsanos et al. [56] found lower SSAAE values during winter compared to the other seasons in Athens due to enhanced levels of biomass-burning aerosols and a pronounced increase in SSAAE values during the morning traffic hours. Furthermore, during the same time frame (20:00-02:00 LST) of 23 March-3 May in 2016-2019, the BC wb fraction was 22%, much lower than the respective one during the lockdown nights, resulting in a higher mean SSAAE of 0.234 (Figure 8c). Therefore, it is concluded that the SSAAE is rather sensitive to the relative changes between fossil fuel and wood-burning emissions [118,123].
Correlations of measured aerosol load, either as columnar AOD or as scattering and absorption coefficients, and its spectral dependence has been widely used for the identification of key aerosol types [89,91,[124][125][126]. Therefore, an investigation of the correlations b sca,550 vs. SAE 450-700 and b abs,520 vs. AAE 470-950 during the lockdown and the same period in previous years can reveal modifications in atmospheric composition, attributable to changes in emission sources and not so much to seasonality. This analysis ( Figure 9) shows considerable changes in the distribution of data points between scattering and absorption and between lockdown data (solid red circles) and data from previous years (open gray circles). Focusing on scattering (Figure 9a), the highest b sca,550 values (>80 Mm −1 ) are associated with very high SAE (>2) during the lockdown period, indicating the prevalence of fine-mode aerosols in the Athens atmosphere, while there is also a higher possibility of lower b sca,550 and higher SAE values compared to the same periods (23 March-3 May) in 2017-2019. In the previous years, several cases indicate an enhanced presence of coarsemode dust aerosols corresponding to very low SAE values (<0.5), mostly detected during the dust events in March 2018 [77]. Cases with SAE values in the range of 0.5-1.2, which mostly imply mixing of urban background with transported coarse particles, are also absent during the lockdown period but are fairly common during the spring season [70].
Atmosphere 2021, 12, x FOR PEER REVIEW 16 of 25 [56] found lower SSAAE values during winter compared to the other seasons in Athens due to enhanced levels of biomass-burning aerosols and a pronounced increase in SSAAE values during the morning traffic hours. Furthermore, during the same time frame (20:00-02:00 LST) of 23 March-3 May in 2016-2019, the BCwb fraction was 22%, much lower than the respective one during the lockdown nights, resulting in a higher mean SSAAE of 0.234 (Figure 8c). Therefore, it is concluded that the SSAAE is rather sensitive to the relative changes between fossil fuel and wood-burning emissions [118,123]. Correlations of measured aerosol load, either as columnar AOD or as scattering and absorption coefficients, and its spectral dependence has been widely used for the identification of key aerosol types [89,91,[124][125][126]. Therefore, an investigation of the correlations bsca,550 vs. SAE450-700 and babs,520 vs. AAE470-950 during the lockdown and the same period in previous years can reveal modifications in atmospheric composition, attributable to changes in emission sources and not so much to seasonality. This analysis (Figure 9) shows considerable changes in the distribution of data points between scattering and absorption and between lockdown data (solid red circles) and data from previous years (open gray circles). Focusing on scattering (Figure 9a), the highest bsca,550 values (>80 Mm −1 ) are associated with very high SAE (>2) during the lockdown period, indicating the prevalence of fine-mode aerosols in the Athens atmosphere, while there is also a higher possibility of lower bsca,550 and higher SAE values compared to the same periods (23 March-3 May) in 2017-2019. In the previous years, several cases indicate an enhanced presence of coarse-mode dust aerosols corresponding to very low SAE values (<0.5), mostly detected during the dust events in March 2018 [77]. Cases with SAE values in the range of 0.5-1.2, which mostly imply mixing of urban background with transported coarse particles, are also absent during the lockdown period but are fairly common during the spring season [70]. Regarding absorption (Figure 9b), a remarkable modification of the 2017-2019 scatter plot is observed during the lockdown period, with main characteristics the accumulation of data points at lower babs,520 values and an enhanced AAE (>1.6) for the highest babs,520 values (>50 Mm −1 ). The latter indicates cases with enhanced recreational RWB emissions (e.g., fireplaces), as citizens stayed more at their houses. Another major difference between the two examined periods is the near absence of data points characterized by high babs,520 (>60 Mm −1 ) and AAE around 1.1-1.2 during the lockdown. These data points typically correspond to enhanced emissions from fossil-fuel combustion along with weak influence from the secondary aerosol formation. In ambient conditions, the AAE related to fossilfuel (normally close to 1) may be modified depending on aerosol size, BC coating and Regarding absorption (Figure 9b), a remarkable modification of the 2017-2019 scatter plot is observed during the lockdown period, with main characteristics the accumulation of data points at lower b abs,520 values and an enhanced AAE (>1.6) for the highest b abs,520 values (>50 Mm −1 ). The latter indicates cases with enhanced recreational RWB emissions (e.g., fireplaces), as citizens stayed more at their houses. Another major difference between the two examined periods is the near absence of data points characterized by high b abs,520 (>60 Mm −1 ) and AAE around 1.1-1.2 during the lockdown. These data points typically correspond to enhanced emissions from fossil-fuel combustion along with weak influence from the secondary aerosol formation. In ambient conditions, the AAE related to fossil-fuel (normally close to 1) may be modified depending on aerosol size, BC coating and atmospheric processing [127,128] and increase up to 1.5 depending on the shell material and the lensing effect [129]. The modification of the b abs,520 vs. AAE 470-950 scatter plot indicates significant changes in the atmospheric chemical composition during the lockdown period in Athens.
The correlation between SAE and AAE values has been extensively used for the identification of key aerosol types [63,[130][131][132][133][134] and may reveal changes in atmospheric composition during the lockdown period, able to modulate the aerosol radiative effect [135]. In this respect, Figure 10a-c shows the correlations between SAE 450-700 and AAE 470-950 for the lockdown period and for the same periods in previous years (2017-2019), aiming to explore possible changes that signify a modification in the dominant aerosol types in Athens. The variety of aerosol types in Athens leads to a high scatter in both data cases. Focusing on the SAE vs. AAE scatter plots (Figure 10a-c, the results point to a higher homogeneity during the lockdown period with an absence of data points lying in the SAE < 1 area (mostly super-micron aerosols). In addition, there is a higher frequency of SAE > 1.5 and AAE > 1.5, thus increasing the fraction of the mixing of brown carbon (BrC) and BC aerosol [70,131] that corresponds to wood-burning ( Figure 10a). These data points are also associated with the highest b sca,550 values (Figure 10b), indicating turbid atmospheric conditions under RWB emissions on certain nights during the lockdown period ( Figure S4). During 2017-2019 (Figure 10c), high b sca,520 values were mostly related to low AAE (close to 1.1) and/or to low SAE (close to zero; dust events). atmospheric processing [127,128] and increase up to 1.5 depending on the shell material and the lensing effect [129]. The modification of the babs,520 vs. AAE470-950 scatter plot indicates significant changes in the atmospheric chemical composition during the lockdown period in Athens. The correlation between SAE and AAE values has been extensively used for the identification of key aerosol types [63,[130][131][132][133][134] and may reveal changes in atmospheric composition during the lockdown period, able to modulate the aerosol radiative effect [135]. In this respect, Figure 10a-c shows the correlations between SAE450-700 and AAE470-950 for the lockdown period and for the same periods in previous years (2017-2019), aiming to explore possible changes that signify a modification in the dominant aerosol types in Athens. The variety of aerosol types in Athens leads to a high scatter in both data cases. Focusing on the SAE vs. AAE scatter plots (Figure 10a-c, the results point to a higher homogeneity during the lockdown period with an absence of data points lying in the SAE < 1 area (mostly super-micron aerosols). In addition, there is a higher frequency of SAE > 1.5 and AAE > 1.5, thus increasing the fraction of the mixing of brown carbon (BrC) and BC aerosol [70,131] that corresponds to wood-burning ( Figure 10a). These data points are also associated with the highest bsca,550 values (Figure 10b), indicating turbid atmospheric conditions under RWB emissions on certain nights during the lockdown period ( Figure S4). During 2017-2019 (Figure 10c), high bsca,520 values were mostly related to low AAE (close to 1.1) and/or to low SAE (close to zero; dust events).

Conclusions
In this study, we examined the effects of the spring 2020 lockdown (23 March-3 May) on spectral aerosol scattering and absorption properties in Athens, Greece, with respect to the pre-(1-22 March) and post-(4-31 May) lockdown periods. The comparison was repeated for the same periods of the previous 3-4 years, allowing for an assessment of the effects of restriction measures on aerosol properties. Moreover, similar analyses were performed at the remote background site of Finokalia, Crete, to assess the lockdown impacts on a regional level. The analysis revealed a pronounced influence of reduced anthropogenic emissions in urban areas on aerosol properties, which in turn are closely related to radiative forcing. Remarkable changes were also observed for intensive properties related to aerosol size, shape and chemical composition, like scattering and absorption Ångström exponents (SAE, AAE) and single-scattering albedo (SSA).
At the remote site of Finokalia, the effect of lockdown did not significantly affect the b sca and BC components, while the variations in scattering and absorbing aerosols were mostly attributed to natural causes like dust events. This feature was different from that observed in Athens, where the locally emitted aerosols from the traffic sector mainly affected the light absorption rather than the scattering, with morning peaks in b abs being much more pronounced than those of b sca . Therefore, the drastic limitation of circulation during the lockdown period caused a much larger reduction in b abs (reaching 60% during the morning hours with respect to previous years) compared with b sca , which presented only a 25-30% decrease during the same time frame. Overall, the much lower mean reduction rate in scattering (−21%) compared to absorption (−36%) can be attributed to the more regional character of scattering aerosols around the greater Athens area. The absorbing properties seem to be highly related to vehicular emissions and mostly to BC from fossil fuel combustion, which also revealed a pronounced decrease reaching −65% during the morning traffic hours. On the contrary, the wood-burning BC fraction (BC wb ) did not present notable changes, while it was found to slightly increase during the night hours, with respect to 2016-2019 mean, as a consequence of increased emissions from residential wood-burning in fireplaces due to stay-at-home measures. In addition, the traffic restrictions affected the intensive aerosol properties in Athens, but to a lower degree, since an increase was observed for SAE (6.1%), AAE (6.3%) and SSA (5.9%) during the lockdown period with respect to previous years. This could be associated with the larger reductions in absorbing rather than scattering aerosols and lowering concentrations of resuspended road dust due to the limitation in road traffic. Consequently, the aerosol during the lockdown period in Athens was dominated by finer, less-absorbing particles and enhanced residential wood-burning emissions during nighttime.
Overall, the lower urban aerosol loading may also affect the aerosol radiative forcing (ARF), leading to smaller ARF values in the surface and top of the atmosphere due to the presence of less-absorbing aerosols during the lockdown period. These potential effects that have been considered of a transient character should be reevaluated given the renewed restriction in the second half of 2020 and in 2021. Since the start of November 2020, Greece has entered a second national lockdown, which is foreseen to extend well into the new year. Therefore, it will be challenging to extend the results of the present study, exploring possible modifications in atmospheric composition and aerosol properties and types during winter, when local pollution events are normally more frequent and more intense.  Table S1: Mean ± standard deviation values of scattering, absorption, SAE, AAE and SSA during the pre-lockdown (1-22 March 2020), lockdown (23 March-3 May 2020) and post-lockdown (4-31 May 2020) periods and associated changes at Thissio, Athens. The asterisk denotes statistically signifi-cant difference (0.95 confidence level), Table S2: Mean ± standard deviation values of scattering, absorption, SAE, AAE and SSA for 2020 and previous years, corresponding to the pre-lockdown , lockdown (23 March-3 May) and post-lockdown (4-31 May) periods. Displaying also the% differences in 2020 compared with the mean of previous years. The asterisk denotes statistically significant difference (0.95 confidence level), Figure  Funding: This research was funded by the project "PANhellenic infrastructure for Atmospheric Composition and climatE change" (MIS 5021516), which is implemented under the Action "Reinforcement of the Research and Innovation Infrastructure", funded by the Operational Program "Competitiveness, Entrepreneurship and Innovation" (NSRF 2014-2020) and co-financed by Greece and the European Union (European Regional Development Fund).

Informed Consent Statement: Not applicable.
Data Availability Statement: Data are available after request.