Aerosol optical properties calculated from size distributions, ﬁlter samples and absorption photometer data at Dome C, Antarctica, and their relationships with seasonal cycles of sources

. Optical properties of surface aerosols at Dome C, Antarctica, in 2007–2013 and their potential source areas are presented. Scattering coefﬁcients ( σ sp ) were calculated from measured particle number size distributions with a Mie code and from ﬁlter samples using mass scattering efﬁciencies. Absorption coefﬁcients ( σ ap ) were determined with a three-wavelength Particle Soot Absorption Photometer (PSAP) and corrected for scattering by using two different algorithms. The scattering coefﬁcients were also compared with σ sp measured with a nephelometer at the South Pole Station (SPO). The minimum σ ap was observed in the austral autumn and the maximum in the austral spring, similar to other Antarctic sites. The darkest aerosol, i.e., the lowest single-scattering albedo ω o ≈ 0 . 91, was observed in September and October and the highest ω o ≈ 0 . 99 in February and March. The uncertainty of the absorption Ångström exponent α ap is high. The lowest α ap monthly medians were observed in March and the highest in August–October. The equivalent black carbon (eBC) mass concentrations were compared with eBC measured at three other Antarctic sites: the SPO and two coastal sites, Neumayer and Syowa. The maximum monthly median eBC concentrations are almost the same ( ∼ 3 ± 1 ng m − 3 ) at all these sites in October–November. This suggests that there is no signiﬁcant difference in eBC concentrations between the coastal and plateau sites. The seasonal cycle of the eBC mass fraction exhibits a minimum f (eBC) ≈ 0 . 1 % in February–March and a maximum ∼ 4 %–5 % in August–October. Source areas were calculated using 50 d FLEXPART footprints. The highest eBC concentrations and the lowest ω o were associated with air masses coming from South America, Australia and Africa. Vertical simulations that take BC particle removal processes into account show that there would be essentially no BC particles arriving at Dome C from north of latitude 10 ◦ S at altitudes < 1600 m. The main biomass-burning regions Africa, Australia and Brazil are more to the south, and their smoke plumes have been observed at higher altitudes than that, so they can get transported to Antarctica. The seasonal cycle of BC emissions from wildﬁres and agricultural burning and other ﬁres in South America, Africa and Australia was calculated from data downloaded from the Global Fire Emissions Database (GFED). The maximum total emissions were in August–September, but the peak of monthly average eBC concentrations is observed 2–3 months later in November, not only at Dome C, but also at the SPO and the coastal stations. The air-mass residence-time-weighted BC emissions from South America are approximately an order of magnitude larger than from Africa and Oceania, suggesting that South American BC emissions are the largest contributors to eBC at Dome C. At Dome C the maximum and minimum scattering coefﬁcients were observed in austral summer and winter, respectively. At the SPO σ sp was similar to that observed at Dome C in the austral summer, but there was a large difference in winter, suggesting that in winter the SPO is more inﬂuenced by sea-spray emissions than Dome C. The seasonal cycles of σ sp at Dome C and at the SPO were compared with the seasonal cycles of secondary and primary marine aerosol emissions. The σ sp measured at the SPO correlated much better with the sea-spray aerosol emission ﬂuxes in the Southern Ocean than σ sp at Dome C. The seasonal cycles of biogenic secondary aerosols were estimated from monthly average phytoplankton biomass concentrations obtained from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) satellite sensor data. The analysis suggests that a large fraction of the biogenic scattering aerosol observed at Dome C has been formed in the polar zone, but it may take a month for the aerosol to be formed, be grown and get transported from the sea level to Dome C.


Introduction
The Antarctic interior region has scarce observations of atmospheric constituents, and many aspects of the atmospheric properties are underdetermined. The Antarctic dome or the polar vortex, which is much stronger than its northern counterpart and present throughout the year (Karpetchko et al., 2005), at most times efficiently prevents transport into the Antarctic troposphere from lower latitudes. However, wildfires and agricultural burning emissions from Africa, South America and Australia do affect vast regions of the Southern Hemisphere, including Antarctica. For instance, Hara et al. (2010) found that haze episodes at Syowa Station, during which visibility can drop to 10 km for periods of ∼ 30 h, were caused by biomass-burning aerosol from South America transported to the Antarctic coast via the eastward approach of cyclones. At the Neumayer station large-scale meridional transport of biomass-burning-derived black carbon, preferentially from South America, seems to determine the BC burden and causes a distinct and consistent spring/early summer concentration maximum (Weller et al., 2013).
Concordia Station lies on Dome C (75 • 06 S, 123 • 23 E) at 3233 m a.s.l. (above sea level) on the East Antarctic Plateau, about 1100 km from the nearest coastline, the Ross Sea. The base is French and Italian operated, with research fields within astronomy and glaciology as well as atmospheric sciences. The atmospheric instrumentation is located in a small cabin southwest of the main base (at the site described by Udisti et al., 2012), where it is upwind of the base at the prevailing wind directions. Concordia is one of only three permanent year-round stations operated on the Antarctic Plateau, the others being the American Amundsen-Scott Observatory (South Pole Station (SPO), 2835 m a.s.l., about 1300 km from the nearest open sea, 1600 km away from Dome C) and the Russian Vostok station (78 • 28 S, 106 • 51 E, 3488 m.a.s.l., 600 km away). Thus, there are large spatial distances between the continuous atmospheric observations. However, properties of the Antarctic atmosphere tend to extend over both longer temporal and spatial scales than elsewhere (Fiebig et al., 2014), suggesting that the scarce observations that exist can be assumed to be representative of larger areas than is typical in other climate regions. This would imply that Dome C is an important indicator for the entire Antarctic inland. Though measurement conditions are harsh, the continuous long-term monitoring provided here can be a baseline for the aerosol optical properties of the Antarctic inland and may provide indications of changes in atmospheric constituents and aerosol levels.
There are several studies on the aerosol chemical composition at Dome C (e.g., Jourdain et al., 2008;Becagli et al., 2012Becagli et al., , 2021Udisti et al., 2004Udisti et al., , 2012Legrand et al., 2016Legrand et al., , 2017a, and the aerosol optical depth (AOD) has also been measured there (Tomasi et al., 2007). However, in situ surface aerosol scattering and absorption coefficients at Dome C have not been presented. The light absorption coefficient and particle number size distributions (PNSDs) have been measured continuously with a three-wavelength Particle Soot Absorption Photometer (PSAP) and a differential mobility particle sizer (DMPS) since 2007. The PNSD data have already been used in several papers. Järvinen et al. (2013) analyzed the seasonal cycle and modal structure of PNSDs measured with the DMPS, Chen et al. (2017) analyzed the num-ber size distribution of air ions measured with an air ion spectrometer (AIS) and the PNSD measured with the DMPS, and Lachlan-Cope et al. (2020) used the Dome C DMPS data for comparing with the PNSD measured at the coastal site Halley. The PSAP data, however, have not been presented in detail. Caiazzo et al. (2021) used some of the PSAP data mainly for evaluating elemental carbon (EC) sample contamination. Grythe (2017) used the data from 2007 to 2013 as part of his PhD thesis, but in the present paper we will analyze that period in more detail. Here we will describe the methods for measuring absorption and calculating scattering from the size distributions and filter samples.
The goals of the paper are to present descriptive statistics of extensive and intensive aerosol optical properties at Dome C in 2007Dome C in -2013, their seasonal cycles, and the relationships between the seasonal cycles of major sources of absorbing and scattering aerosols. The aerosol optical properties (AOPs) will be compared with other observations from other Antarctic sites, in most detail the scattering coefficients measured at the South Pole.

Sampling site
Concordia Station is a permanently operated French and Italian Antarctic research base on the East Antarctic Plateau. The observations are performed at isolated sites around the main base. The Dome C sampling site is the same as used by Udisti et al. (2012), Becagli et al. (2012), and Järvinen et al. (2013). It is located about 1 km southwest of the station's main buildings, upwind in the direction of the prevailing wind. The northeastern direction (10-90 • ) has been declared the contaminated sector. Below the validity of the contaminated sector will be analyzed by using the absorption photometer data. For in situ aerosol instrumentation the sample air was taken at the flow rate of 5 L min −1 from the roof of the cabin with a straight 2 m-long 25 mm diameter stainless steel tube inlet. It was covered with a protective cap to protect against snowfall and ice buildup.

Aerosol measurements
Light absorption by particles was measured with a Radiance Research 3λ PSAP at three wavelengths, λ = 467, 530, and 660 nm. There was no nephelometer measuring scattering coefficient, so it was calculated from particle size distributions and filter sample data as described below. Particle number size distributions were measured at 10 min time resolution in the size range 10-620 nm with a custom-built differential mobility particle sizer (DMPS) as described by Järvinen et al. (2013) and in the size range 0.3-20 µm with a Grimm model 1.108 optical particle counter (OPC) in 2007-2009. RH was not measured in the Dome C sample air, but it can be safely claimed that it was dry. The absolute humidity in the air on the upper plateau is very low, and temperature varies from colder than about −20 • C in the austral summer down to about 80 • C in the austral winter. When air is sampled to the instruments in the measurement containers where temperature is > +10 • C, RH decreases to very low values. In addition to the in situ instruments, PM 1 and PM 10 filter samples were collected for chemical analyses by ion chromatography. The length of the sampling period of the PM 1 and PM 10 samples was 3 or 4 and 1 d, respectively.
The data coverage for the PSAP, the DMPS, the OPC, and the PM 1 and PM 10 filter sample data are presented in Fig. 1. The number of hours of accepted data and the number of samples are shown in parentheses for the continuous instruments and the filter samplers, respectively. The filtering criteria will be presented below (Sect. 2.4).

Meteorological measurements
Ambient air temperature (t), relative humidity (RH), wind speed (WS) and wind direction (WD) data were from the routine meteorological observation at the Concordia Station as part of the IPEV/PNRA project -a collaborative project between Programma Nazionale di Ricerche in Antartide (PNRA) and Institut Polaire Français Paul-Emile Victor (IPEV) (https://www.climantartide.it/, last access: 10 April 2022).

Data processing 2.3.1 Mass concentrations from size distributions
Sixty-minute average size distributions n(D p ) were first calculated from the original 10 min data and corrected for standard temperature and pressure (STP) (p = 1013 hPa, T = 273.15 K). The DMPS n(D p ) data were corrected for diffusion losses during the inversion (Järvinen et al., 2013;Chen et al., 2017). Mass concentrations were calculated from the number size distributions measured with the DMPS from m(DMPS) = ρ p V (DMPS) = ρ p 620 nm 10 nm π 6 D 3 p n(D p )dD p , (1) where the density ρ p = 1.7 g cm −3 was used. For a particle density of 1.7 g cm −3 , the particle diameter 620 nm corresponds to the aerodynamic diameter D a = ρ p /ρ 0 D p = √ 1.7 × 620 nm ≈ 808 nm, where ρ 0 = 1 g cm −3 . To be consistent with the definitions of filter-sample size ranges that typically show the upper aerodynamic diameter of a sampler inlet, the mass concentration calculated from Eq. (1) will be referred to as m(DMPS,PM 0.8 ) and the volume concentration as V (DMPS,PM 0.8 ).
In December 2007-July 2009 particles were also measured with the Grimm 1.108 OPC that measures number concentrations of particles in the D p range of 0.3-20 µm. The

5036
A. Virkkula et al.: Optical properties of aerosols at Dome C Figure 1. The periods of the PSAP, the DMPS, the Grimm OPC and the PM 1 and PM 10 filter sample data. The number of hours of accepted data and the number of samples are shown in parentheses for the continuous instruments and the filter samplers, respectively. particle number concentrations in the size range D p > 1 µm were first corrected for WS-dependent and particle-diameterdependent inlet and sampling tube losses by dividing the raw, non-corrected number concentrations n(D p , OPC, noncorrected) with the combined inlet and tube transmittance f inlet,tubing (WS, D p ), as described in the Supplement. The number concentrations were very small in the size ranges where the transmittance losses were significant. In a large fraction of data n(D p , OPC, non-corrected) was zero in the particle size range where f inlet,tubing is small. If the true concentration was larger than zero but the raw concentration in the OPC data was zero due to the instrument sensitivity and sampling losses, then the corrected concentration would also be zero even if the raw concentration was multiplied by a very large number 1/f inlet,tubing . Consequently, the number concentrations and the derived mass concentrations and scattering coefficients in the large-particle size range would be underestimated. The underestimation could in principle be estimated by using a collocated more sensitive instrument sampling air through a well-defined inlet with minimal particle losses. These were not available, so a detailed analysis of the underestimations of the derived quantities was omitted from the paper.
The three largest channels of the OPC measure the number concentrations in the D p range of 7.5-20 µm. For an assumed density ρ p = 1.7 g cm −3 the diameter D p = 7.5 µm corresponds to the aerodynamic diameter D a = 9.8 µm. Assuming that ρ p is constant over the whole size range, the mass concentration of particles smaller than D a = 10 µm is calculated from the number size distributions by excluding the three largest particle OPC channels as m n(D p ), PM 10 = ρ p V n(D p ), PM 10 = ρ p V (DMPS,PM 0.8 ) + V (OPC,PM 0.8−10 ). (2) The fraction of volume concentration measured by the DMPS equals This fraction was calculated from data collected during the simultaneous operation of the DMPS and the OPC. The monthly average f V (DMPS) values presented in Table 1 were used for the period 2008-2013 to calculate mass concentrations in the size range D a < 10 µm from In other words, the variable names m(DMPS,PM 0.8 ) and m(DMPS,PM 10 ) will be used below to emphasize that these mass concentrations were calculated from DMPS data. The mass concentrations m(DMPS,PM 0.8 ) and m(DMPS,PM 10 ) can be considered to be the lower and upper estimates of m.

Scattering coefficients from the size distributions
Scattering coefficients were calculated using the 60 min average size distributions from where Q s is the scattering efficiency calculated using the Mie code by Barber and Hill (1990), m is the refractive index, λ is the wavelength, and n(D p ) is the particle number size distribution. Analogously to the mass concentrations, the scattering coefficients were determined from the simultaneous DMPS and OPC measurements in December 2007-July 2009 from σ sp n(D p ), PM 10 = σ sp (DMPS,PM 0.8 ) where σ sp (OPC,PM 0.8 ) and σ sp (OPC,PM 0.8−10 ) are the scattering coefficient calculated from the particle number size distributions in the size ranges measured by the DMPS and the OPC, respectively. As explained above, the number size distributions for D p > 1 µm were corrected for the inlet and sampling tube losses. For σ sp (DMPS,PM 0.8 ) the refractive index of sulfuric acid (SA, H 2 SO 4 , m r = 1.426 + 0i, Seinfeld and Pandis, 2006) was used. This refractive index is slightly lower than that estimated for submicron aerosols at two low-altitude Antarctic stations, Aboa and Neumayer in Queen Maud Land. Virkkula et al. (2006) measured particle number size distributions in the size range D p < 800 nm With an iteration procedure matching nephelometer-measured and sizedistribution-derived scattering coefficients, the real refractive indices were 1.43 ± 0.07 and 1.45 ± 0.04 at λ = 550 nm for all data and excluding new particle formation, respectively. Jurányi and Weller (2019) measured size distributions with an SMPS and a laser aerosol spectrometer (LAS) for a full year at the coastal site Neumayer and by fitting data of the two instruments in the overlapping range of 120-340 nm obtained m r = 1.44 ± 0.08. Considering that both Aboa and Neumayer are closer to sources of ammonia, which neutralizes aerosol and increases the refractive index above that of pure sulfuric acid (1.426), it was assumed here that the use of 1.426 for the calculation of σ sp from the size range measured with the DMPS is reasonable. For the larger particle size range, σ sp (OPC,PM 0.8−10 ), the refractive index of NaCl (m r = 1.544, Seinfeld and Pandis, 2006) was used. This value is in line with the average refractive index of 1.54 with a range from 1.50 to 1.58 in the particle size range 0.3-12 µm in impactor samples taken at the South Pole (Hogan et al., 1979) and with the supermicron particle refractive index of 1.53 ± 0.02 calculated from the chemical composition of 12-stage impactor samples taken at the coastal site Aboa . The fraction of the scattering coefficient measured by the DMPS was calculated from f σ sp (DMPS, λ) = σ sp (DMPS,PM 0.8 , λ) σ sp n(D p ), PM 10 , λ .
The wavelengths of λ = 467, 530, and 660 nm were used to match the PSAP data. Similar to f V (DMPS), f σ sp (DMPS,λ) was calculated from data collected during the simultaneous operation of the DMPS and the OPC, the seasonal monthly statistics were calculated (Table 1), and the respective monthly averages were applied to the period 2008-2013 to calculate σ sp in the size range D a < 10 µm from σ sp (DMPS,PM 10 , λ) = σ sp (DMPS,PM 0.8 , λ) f σ sp (DMPS, λ) .
The wavelength symbol λ will be used below only when necessary. The variable names σ sp (DMPS,PM 0.8 ) and σ sp (DMPS,PM 10 ) will be used to emphasize that these scattering coefficients were calculated from DMPS data in the aerodynamic particle size ranges D a < 0.8 and D a < 10 µm. The scattering coefficients σ sp (DMPS,PM 0.8 ) and σ sp (DMPS,PM 10 ) can also be considered to be the lower and upper estimates of σ sp at the given wavelength. Figure 2 shows the average particle number, volume, and scattering size distributions at λ = 530 nm in the size range 10 nm-10 µm and the respective normalized cumulative size distributions in the size range of 10 nm-7.5 µm during the period from 14 December 2007 to 14 July 2009 in summer and in winter. Figure 2a and b show that for the number concentrations the OPC size range plays an insignificant role, whereas the larger particles contribute significantly to both total particle volume concentration ( Fig. 2c and d) and scattering coefficients ( Fig. 2e and f), and that this contribution varies seasonally. The contributions of f V (DMPS) and f σ sp (DMPS,λ) were calculated for hourly averaged size distributions from Eqs. (3) and (7), and the monthly seasonal statistics were calculated and presented in Table 1 and as a box plot in Fig. 3. Both the table and the box plot show that both f V (DMPS) and f σ sp (DMPS,λ) have max-  (b, d, f): cumulative fractions of the respective parameters in the size range D p < 7.5 µm, which corresponds to the aerodynamic particle size range D a < 9.8 µm.
ima in summer and minima in winter. They also show that the ranges are large. Consequently, the use of the monthly averages presented in Table 1 for calculating m(DMPS,PM 10 ) and σ sp (DMPS,PM 10 ), Eqs. (4) and (8), creates an additional uncertainty in the results. Another important result is that the wavelength dependency of f σ sp (DMPS,λ) is clear, and it also has a seasonal cycle.

Absorption coefficients and equivalent black carbon concentrations
The PSAP data were first corrected for flow and spot size.
where A is the filter spot area, Q the flow rate, the transmittance, f (Tr) the loading correction function, and t the count integration time. The PSAP reports σ ap,nsc with a 0.1 Mm −1 resolution at a 1 s time resolution. Averaging the 1 s data is not good enough since at Dome C absorption coefficients are most of the time clearly lower than 0.1 Mm −1 . Therefore, the signal and reference counts SIG and REF were used in Eq. (10) with t = 60 min. Manufacturer-cut spots of the standard filter material Pallflex E70-2075W were used in the PSAP. The spot diameter was measured to be 4.9 ± 0.1 mm, so the spot area A was 18.9 ± 0.6 mm 2 . The uncertainty of A is ∼ 3 %. Transmittance is reduced mainly by light absorption but also by scattering aerosol, which results in the so-called apparent absorption and has to be taken into account in the data processing. There are different algorithms for processing PSAP data, e.g., by Bond et al. (1999), Virkkula et al. (2005), Müller et al. (2014), and Li et al. (2020). Here we will use both the algorithms presented by Bond et al. (1999)  B1999) with the adjustment presented by Ogren (2010), and the algorithm presented by Virkkula et al. (2005) with the constants updated by Virkkula (2010) (here V2010): where is the single-scattering albedo and k 0 , k 1 , h 0 , h 1 , and s are wavelength-dependent constants. In the rest of the paper the symbol σ ap,nsc will be used to present the non-scatteringcorrected absorption coefficient, corrected with the constants and formula in Eq. (11) excluding the subtraction of σ sp .
Since there are the above-explained size-dependent uncertainties of the scattering coefficient, additional absorption coefficient estimates were calculated by using both algorithms. The upper estimates of absorption coefficients σ ap (σ sp (DMPS,PM 0.8 )) were calculated by using the lower estimate of the scattering coefficient σ sp = σ sp (DMPS,PM 0.8 ) in the scattering corrections in Eqs. (11) and (12), and the lower estimates of the absorption coefficient σ ap (σ sp (DMPS,PM 10 )) were calculated by using the upper estimate of the scattering coefficient σ sp = σ sp (DMPS,PM 10 ) in the scattering corrections. Consequently, the lower and upper estimates of ω o are denoted ω o (σ sp (DMPS,PM 0.8 )) and ω o (σ sp (DMPS,PM 10 )), respectively. They were calculated by using both Eqs. (11) and (12) for calculating σ ap .
Considering that the period with the simultaneous measurements with the DMPS and the OPC showed that the DMPS size range always leads to an underestimation of both aerosol mass and scattering coefficient, it is likely that σ ap corrected for scattering with σ sp (DMPS, PM 10 ) is closer to the true σ ap than that corrected with σ sp (DMPS, PM 0.8 ). In the results σ ap,nsc , σ ap (σ sp (DMPS,PM 0.8 )), and σ ap (σ sp (DMPS,PM 10 )) will be presented to evaluate the effect of using only the size range measured with the DMPS for the scattering correction.
Similarly to σ sp , the wavelength dependency of light absorption by particles can roughly be described by the absorption Ångström exponent α ap = − ln σ ap (λ 1 )/σ ap (λ 2 ) ln (λ 1 /λ 2 ) that was calculated by using λ = 467 and 660 nm for σ ap,nsc , σ ap (σ sp (DMPS,PM 10 )), and both Eqs. (11) and (12). The variable names α ap (σ ap,nsc ), α ap (σ sp (DMPS,PM 10 ),B1999), and α ap (σ sp (DMPS,PM 10 ),V2010, respectively, will be used to denote the α ap calculated in different ways. These calculations were conducted to study the uncertainty of α ap due to scattering corrections. The absorption coefficient was used to estimate the concentration of equivalent black carbon, eBC , from where MAC is the mass absorption coefficient. For freshly emitted BC the MAC value is approximately 7.5 m 2 g −1 at λ = 550 nm (Bond et al., 2013). By assuming a wavelength dependency of λ −1 , this corresponds to MAC ≈ 7.8 m 2 g −1 at λ = 530 nm. This can be considered to yield an upper estimate for eBC concentrations since for coated BC particles MAC is larger (Bond et al., 2013). eBC was calculated by using σ ap,nsc and σ ap (σ sp (DMPS,PM 10 )) calculated with both algorithms, Eqs.
(1) and (2). The corresponding variable names eBC(σ ap,nsc ) and eBC(σ sp (DMPS,PM 10 )) will be used below for them. The scattering-corrected eBC(σ sp (DMPS,PM 10 )) can be considered to be closer to the true eBC concentration. The reason for also presenting eBC(σ ap,nsc ) is that an estimate of BC concentrations is often needed even if it is known that it is an upper estimate (Caiazzo et al., 2021). It is also comparable with the eBC often presented from Aethalometer measurements. Presenting both yields a quantitative estimate of the bias due to not correcting the data for scattering. The eBC mass fractions in the two size ranges D a < 0.8 and D a < 10 µm were calculated from f eBC (m(DMPS,PM 0.8 )) = 100 % eBC m(DMPS,PM 0.8 ) , f eBC (m(DMPS,PM 10 )) = 100 % eBC m(DMPS,PM 10 ) , where the mass concentrations m(DMPS,PM 0.8 ) and m(DMPS,PM 10 ) were defined in Eq. (4) and eBC calculated from Eq. (15). Mass fractions were calculated for eBC(σ ap,nsc ) and eBC(σ sp (DMPS,PM 10 )).

Noise of scattering and absorption coefficients and eBC
The uncertainty of scattering coefficients should in principle be calculated from the error propagation formula δσ sp 2 = ∂σ sp /∂x i 2 (δx i ) 2 , where δx i is the uncertainty of variable x i in calculating σ sp (e.g., Sherman et al., 2015). That would require taking into account all uncertainties of the size distribution measurements and Mie modeling. However, a simplified approach was used here. The σ sp calculated from the size distribution data and the uncertainty of the size distribution range were used for calculating lower and upper estimates of σ sp as explained above. In addition to that, the noise of σ sp was estimated from the average of the absolute differences of both consecutive hourly averaged scattering coefficients δσ sp (average,1 h) = average(| σ sp (1 h)|) = average(|σ sp (t i+1 )σ sp (t i )|). The average noise of 24 h averages was calculated from δσ sp (24 h) = δσ sp (average,1 h) / √ 24. The noises were calculated for both σ sp (DMPS,PM 0.8 )) and σ sp (DMPS,PM 10 ). The noises are presented in Table 2. Note that the difference |σ sp (t i+1 ) − σ sp (t i )| is not only due to random noise, so higher | σ sp | values are observed when σ sp is in reality increasing or decreasing, so the true random noise is slightly lower. When σ sp is used in calculating the scattering correction of σ ap in B1999 (Eq. 11), σ sp is multiplied by 0.0164. Consequently, the σ sp noise for the 24 h averages results in a 0.0164σ sp noise for σ ap . These noises are also presented in Table 2.
The uncertainty of the absorption coefficient should also be calculated from the error propagation formula, similarly to Sherman et al. (2015). However, here only the uncertainties of the spot size (∼ 3 %) and the statistical noise are taken into account. The noise of the non-scattering-corrected hourly σ ap,nsc was estimated from the average of the absolute differences of both consecutive absorption measurements δσ ap,nsc (average) = average(| σ ap,nsc |) = average(|σ ap,nsc (t i )σ ap,nsc (t i )|), similarly to the noise estimate of σ sp . The noise of 24 h averages was estimated from δσ ap,nsc (24 h) = δσ ap,nsc (average,1 h) / √ 24. The noises in the scattering-corrected absorption coefficients were calculated from δσ ap = δσ ap,nsc + 0.0164δσ sp for both σ ap (σ sp (DMPS,PM 0.8 )) and σ ap (σ sp (DMPS,PM 10 )) and for 1 and 24 h averages ( Table 2). The noise determined this way is formally correct only for σ ap calculated with the B1999 formula, Eq. (11), not for V2010. However, calculated directly from the absolute differences, the average | σ ap (B1999)| ≈ average | σ ap (V2010)|, but the contribution of scattering to the noise was only determined for B1999, as explained above. For V2010, Eq. (12), a formal error propagation calculation is more complicated due to the iterative form of the procedure, and it is beyond the scope of the present paper. The noise of eBC was calculated from δ(eBC(σ ap )) = δσ ap / MAC for both non-scatteringcorrected and scattering-corrected eBC. The detection limits were defined as 2 × δ(eBC(σ ap )). The results are presented in Table 3.
The largest uncertainty factor for σ ap , ω o , α ap , and eBC is not related to noise. It is due to the uncertainty of the refractive index and size distributions used for calculating σ sp and the algorithm. This was evaluated by calculating σ ap by using the lower and upper estimates of σ sp in both scattering correction algorithms. These four values were used then for calculating ω o , α ap , and eBC, and they are presented below in relevant tables and figures.

Filtering and preprocessing the in situ data
Both PSAP absorption and DMPS-derived scattering coefficient data were filtered manually by removing rapidly changing values since they can be assumed to result from contamination from the station or from some technical problem. The PSAP transmittance data were used to filter out data measured at Tr < 0.7 following recommendations in WMO/GAW Report No. 227 (2016) and the PSAP handbook (Springston, 2018). During most of 2010 the PSAP flow was extremely unstable, so practically the whole year was removed.
All major sources of light-absorbing aerosol other than the Dome C base are so far away that rapid variations in σ ap,nsc are due to either instrument malfunction or influence from the base, for instance, emissions from vehicles. Further filter- Table 2. Noise of scattering and absorption coefficients calculated from the particle number size distributions and the PSAP data for two averaging times ( t = 1 and 24 h). Noise was estimated for the scattering coefficients in the two size ranges (σ sp (DMPS,PM 0.8 ) and σ sp (DMPS,PM 10 )) and for absorption coefficients without scattering corrections (σ ap,nsc ) and with scattering corrections (σ ap (σ sp (DMPS,PM 10.8 )) and σ ap (σ sp (DMPS,PM 10 ))) as explained in Sect. 2.5.
Average noise of absorption coefficients without and with scattering corrections, Mm −1 δσ ap,nsc δσ ap (σ sp (PM 0.8 )) δσ ap (σ sp (PM 10 )) ing of the data was done by removing data in which 10 min averages of σ ap,nsc were more than 10 times larger than the hourly σ ap,nsc . This was done to remove short events that are local but that do not appear to come directly from the base, based on wind direction. In all, roughly 13 % of the data were deemed contaminated. Additionally, wind data were used to remove clear contamination from the station. The sampling site is located upwind of the base itself by the prevailing wind directions. The base has a year-round diesel generator, and vehicles operated within the base area move around the base from November to February. Figure 4 shows the distribution of σ ap,nsc in 5 • wind direction (WD) sectors at wind speed WS > 2 m s −1 . The generator at the base is clearly observed as a pronounced peak of σ ap,nsc at WD 60 • . If the 75th percentile of the σ ap,nsc cumulative distribution were used as the criterion for the contaminated sector when the wind direction was 30 • < WD < 90 • , then 6 % of the data would be filtered. If the 99th percentile of σ ap,nsc is used, the contamination sector is wider, 20 • < WD < 110 • , and 10 % of the data would be filtered. Here the latter, i.e., the stricter criterion, was used. The distribution of σ ap,nsc in the same WD sectors at several wind speed intervals is shown in Supplement Fig. S1. It is obvious that at low wind speeds contaminated air can come from all directions. Therefore, when WS < 2 m s −1 , all data were filtered out, regardless of WD.
Since the size distribution and absorption measurements are done in the same cabin, the DMPS and OPC data were also removed when the PSAP observations indicated contamination. Figure 1 shows the instruments' operational time in hours. The DMPS measurements had more gaps than the PSAP. The three instruments required for a valid measure-ment were not always operational at the same time. After filtering, altogether 15 815 h of data remained for the statistical analyses. No filtering was applied to the PM 1 and PM 10 filter samples. The contamination is mainly BC, so it was to be assumed that the effect on ion concentrations was not significant.
The calculations were done using hourly averaged data. These data were filtered to remove contaminated data as explained above. The filtered data were then averaged over 24 h to reduce noise and improve detection limits. In the discussions below, the running 24 h averages were used, centered at each hour, i.e., σ ap (t, 24 H) = average(σ ap (t − 12, 1 H), . . ., σ ap (t + 11, 1 H)), which means, for instance, that at noon σ ap (t = 12, 24 H) = average(σ ap (t = 0, 1 H), . . . , σ ap (t = 23, 1 H)), so the noon average represents all absorption coefficients measured during that day. If, during any period to be averaged, there were less than 12 h of noncontaminated data, then that 24 h average was excluded from further analysis.

Filter sample analyses and data processing
There were two samplers in the immediate vicinity of the cabin where the other in situ measurements were made. There was a PM 10 sampling head operating following the CSN EN 12341 European Standard. The PM 1 samples were collected on the backup filter of a Dekati PM 10 impactor. In both of these, particles were sampled on Teflon filters (Pall-Gelman, 47 mm diameter, 2 µm nominal porosity). PM 10 and PM 1 load is obtained by summing the mass of the ions determined on Teflon filters. Note that this can be considered to  Table 3. Noise and detection limits of eBC concentrations calculated from the noise of the absorption coefficients presented in Table 2. eBC noise, ng m −3 Averaging time δ(eBC(σ ap,nsc )) δ(eBC(σ ap (σ sp (DMPS,PM 0.8 )))) δ(eBC(σ ap (σ sp (DMPS,PM 10 )))) Averaging time 2δ(eBC(σ ap,nsc )) 2δ(eBC(σ ap (σ sp (DMPS,PM 0.8 )))) 2δ(eBC(σ ap (σ sp (DMPS,PM 10 )))) be the lower estimate since there could be unidentified compounds, such as organic carbon on the filters. Just before the analysis, half of each filter was extracted with 10 L of ultrapure water (18 M Milli-Q) in an ultrasonic bath for 20 min. Every filter manipulation was carried out under a class-100 laminar-flow hood to minimize contamination risks. Inorganic anions and cations, as well as selected organic anions, were simultaneously measured by using a three Thermo Scientific Dionex ion-chromatograph (IC) system, equipped with electrochemically suppressed conductivity detectors. The sample handling during the IC injection was minimized by using a specifically designed Flow-Injection Analysis (IC-FIA) device (Morganti et al., 2007). Cations (Na + , NH + 4 , K + , Mg 2+ , and Ca 2+ ) were determined by using a Thermo Scientific Dionex CS12A-4 mm analytical column with 20 mM H 2 SO 4 eluent. Inor-ganic anions (Cl − , NO − 3 , SO 2− 4 , and C 2 O 2− 4 ) were measured by a Thermo Scientific Dionex AS4A-4 mm analytical column with a 1.8 mM Na 2 CO 3 /1.7 mM NaHCO 3 eluent. F − and some organic anions (acetate, glycolate, formate, and methanesulfonate) were determined by a Thermo Scientific Dionex AS11 separation column by a gradient elution (0.075-2.5 mM Na 2 B 4 O 7 eluent). Further details on the ion chromatographic measurements are reported in Udisti et al. (2004) and Becagli et al. (2011Becagli et al. ( , 2021. All concentrations were corrected to STP (1013.25 hPa and 273.15 K). The ion data used in the present work are a subset of the data from 2005 to 2013 that Becagli et al. (2021) used for an analysis of the relationships between non-sea-salt sulfate, MSA, biogenic sources, and environmental constraints.
In addition to calculating scattering coefficients from the DMPS data, PM 1 and PM 10 mass concentrations were also used for calculating scattering coefficients. The scattering coefficients were calculated by multiplying the mass concentrations by mass scattering efficiencies (MSEs) presented by Hand and Malm (2007). The PM 10 mass concentrations were multiplied by the mass scattering efficiency of 1.9 m 2 g −1 , and the PM 1 concentrations were multiplied by 3.6 m 2 g −1 . These are the MSE for "total mixed" aerosol and "fine mixed" aerosol in Table 5 in Hand and Malm (2007), respectively. It has to be kept in mind that the MSE values in the above-mentioned paper were derived from measurements in the continental USA, so they most likely have a high uncertainty when applied to the Dome C aerosol. The MSE values presented by Quinn et al. (2002) were used for calculating the scattering coefficient of nss sulfate in PM 1 filters.

Scattering data from the South Pole
At the SPO the light-scattering coefficient has been measured for more than 40 years. An integrating nephelometer was installed in 1979 and used to measure σ sp at four wavelengths (450, 550, 700, 850 nm). This nephelometer (Meteorology Research Inc. (MRI), Altadena, CA) was used until its failure in 2002, and a TSI Model 3563 three-wavelength nephelometer (λ = 450, 550, and 700 nm) replaced it in November 2002 (Sheridan et al., 2016). Running 24 h averages of σ sp (550 nm) were calculated for the years 2007-2013 the same way as was done for the Dome C data. The data were used for comparisons with σ sp calculated from the Dome C data.

Source area analyses
The air-mass history and transport of aerosols to Dome C were calculated with the Lagrangian dispersion model FLEXPART (Stohl et al., 2005;Pisso et al., 2019). ECMWF reanalysis meteorology was used to run 60 000 trajectories every 6 h 50 d backwards from Dome C to make a statistical sampling of the air measured there. The FLEXPART trajectories follow the mean flow of the atmosphere plus random perturbations to account for turbulence.
In backward mode, the FLEXPART output is emission sensitivity S that is proportional to residence time within a grid cell (Stohl et al., 2005;Hirdman et al., 2010;Pisso et al., 2019). Depending on the settings, the output unit of FLEX-PART in the backward runs can be s, s m 3 kg −1 , or s kg m −3 . In the present work the unit of S is seconds (s). When coupled with emissions, FLEXPART emission sensitivity creates a concentration at the release point that is equivalent to forward simulations from emissions, except for some small numerical differences (Seibert and Frank, 2004). One advantage of using a backward simulation in a case like this is that the emission sensitivity fields can be used not only to simulate concentrations, but also directly to quantitatively describe exactly where the air that reaches Dome C originates, and, thus, potential emissions influences. Emission sensitiv-ity close to the surface -here at levels < 1000 m a.g.l. -is often called the footprint (e.g., Hirdman et al., 2010). If a footprint were multiplied by emission mass flux in kg m −3 s −1 at some grid cell, the result would be a concentration due to that emission at the receptor site (Stohl et al., 2005). In the present work, this step was not done.
To investigate the role of removal processes during transport, for all model runs, two different tracers were used, one atmospheric tracer with no removal and simulated BC particles with a lognormal size distribution (geometric mean diameter = 150 nm, geometric standard deviation 1.5) experiencing both dry and wet deposition. All tracers were run backwards for 50 d, in most cases sufficient for the aerosol tracer to have less than 1 × 10 −12 of the emission sensitivity of the inert air tracer, meaning any emission prior to this would have been removed by the time of arrival at Dome C. The wet removal differentiates removal within and below clouds, also considering the water phase of the clouds and the precipitation type. The FLEXPART removal parameters are the efficiency of aerosols in serving as cloud condensation nuclei (CCN eff ) and ice nuclei (IN eff ). The values used for them were CCN eff = 0.9 and IN eff = 0.1 as in Table 4 of Grythe et al. (2017). The FLEXPART below-cloud scavenging is a scheme based on Laakso et al. (2003) and Kyrö et al. (2009), both described in Grythe et al. (2017). The model includes a realistic distribution of clouds by incorporating three-dimensional cloud information from ECMWF. For a detailed description, see Grythe et al. (2017).

Footprint difference calculations
A statistical analysis was applied to differentiate types of air pathways using a method derived from Hirdman et al. (2010). With the main aim of investigating the different pathways to Dome C, each 6 h interval was given a rank with regards to eBC concentration and single-scattering albedo. The emission sensitivity of the 50 d transport for an aerosol tracer was sorted according to its relative type. The emission sensitivities of the highest (S H ) and lowest (S L ) 10 % of eBC concentration and ω o were calculated by averaging their emission sensitivities for a given grid cell I , j , and m by where M is the number of measurements and S * can be any of the sorting criteria. The relative difference between two emission sensitivities S 1 and S 2 in percentage is then calculated as In the calculation the emission sensitivities close to the surface at < 1000 m a.g.l. were used, and so Eq. (19) can be called the relative difference of footprints. This analysis of Table 4. Statistical summary of mass concentrations estimated from particle number size distributions, sums of ion concentrations of PM 1 and PM 10 filter samples, and the PSAP data at Dome C in 2008-2013. The statistical values of the PM 1 and PM 10 are those of all individual filters, and the statistical values calculated from the DMPS and PSAP data are those of running 24 h-averaged data; see details in the text. m(DMPS,PM 0.8 ): mass concentration calculated from the particle number size distributions measured with the DMPS assuming a particle density of 1.7 g cm −3 ; eBC: equivalent black carbon concentration calculated from the absorption coefficients at λ = 530 nm calculated by using the B1999 algorithm without any scattering corrections and with B1999 and V2010 algorithms using σ sp = σ sp (DMPS,PM 10 ) for the scattering corrections and assuming MAC = 7.78 m 2 g −1 . f PM 0.8 : scattering-corrected eBC mass fraction calculated from (eBC/m(DMPS,PM 0.8 )) × 100 %; f PM 10 : scattering-corrected eBC mass fraction calculated from (eBC/m(DMPS,PM 10 )) × 100 %.
Percentiles N Avg ± SD 5th 25th 50th 75th 95th the footprints can be used to differentiate between different influencing factors on the air mass. This can be either the influence of transport, removal or combination of these (transport efficiency), or the emission strength.

Emissions used for interpreting the footprint statistics and observed seasonal cycles
The Global Fire Emissions Database (GFED) is a satellite information-based fire activity map. Monthly gridded burned area and emissions from fires are included in the product (http://www.globalfiredata.org, last access: 10 April 2022). Emitted BC is calculated based on emission factors, which depend on the type of vegetation that is burning. Satellites give snapshots collected to give pseudo global coverage and not continuous coverage. GFED v3.1 is based on the area burned, which is derived by coupling Moderate Resolution Imaging Spectroradiometer (MODIS) fire pixel counts with surface reflectance images (Giglio et al., 2006(Giglio et al., , 2009(Giglio et al., , 2010. This widely used emission inventory has uncertainties that arrive both from the emission factors and also from the amount of burnt material. A comparison of this bottom-up inventory with top-down inventories found large regional differences, and top-down estimates were about 30 % higher (Bond et al., 2013).
For the scattering aerosol two sources were considered. An offline tool (FLEX-SSA) developed by Grythe et al. (2014) and Grythe (2017) to simulate sea-spray aerosol (SSA) with FLEXPART was used. It uses inputs from the ECMWF model. These inputs are the wind speed at 10 m above the surface (U10) and the sea surface temperature (SST). The tool takes into account the sea ice fraction which is important to the Southern Ocean SSA emissions. The other major marine scattering aerosols discussed below are biogenic secondary aerosols. Behrenfeld et al. (2017) Tables 4  and 5. Several observations can be made from the time series in Fig. 5. First, the scattering coefficients calculated from the size distributions and the filter samples follow each other relatively well. There is a clear seasonal cycle of both σ sp and σ ap,nsc . It is clearly seen that σ ap,nsc follows the temporal variations of σ sp (DMPS): the high and low values occur mainly simultaneously, which is good considering that these two AOPs were measured with independent instruments. Since the PSAP and other filter-based absorption photometers are sensitive not only to absorbing but also scattering aerosol, and since Dome C is far from BC sources, it is possible that the good correlation is due to the apparent absorption only. Below, this will be studied simply by using Eqs. (11) and (12) to account for the scattering artifact in the absorption measurement.
The σ sp at Dome C and the SPO agrees better in austral summer than in winter. However, many high-concentration episodes are also observed in winter almost simultaneously at Dome C and the SPO. As an example, a 4-month pe-  Figure 6a shows the upper estimate σ sp = σ sp (DMPS,PM 10 ), the corresponding lower estimate of σ ap = σ ap (σ sp (DMPS,PM 10 )) (corrected according to B1999, Eq. 11), and the upper estimate of ω o . Figure 6b presents the lower estimate of σ sp = σ sp (DMPS,PM 0.8 ), the corresponding upper estimate of σ ap = σ ap (σ sp (DMPS,PM 0.8 )), and the lower estimate of ω o . In both Fig. 6a and b the non-scattering-corrected absorption coefficient σ ap,nsc is also shown.
In Fig. 6a, the period denoted by (1) shows an episode in which ω o decreases significantly for several days, being an example of long-range-transported eBC. Episodes (2) and (3) are examples of periods when σ sp is approximately an order of magnitude higher at the SPO than at Dome C. There are also events such as episode (3) when σ sp is approximately the same at both sites. The peaks often seem to appear slightly earlier at the SPO than at Dome C, suggesting transport from the SPO to Dome C rather than the other way around. An example of this is shown in the footprint (Fig. 7) calculated for the episode denoted by (3) in Fig. 6. The footprint shows that the air masses came from the direction of the Antarctic peninsula via the SPO to Dome C. Air flow from the direction of the Weddell Sea to the SPO and then to Dome C is consistent with a very-long-known winter-time circulation pattern (Alt et al., 1959) as reviewed by Shaw (1979). During the event denoted by episode (3), σ ap,nsc was also high. However, when the scattering correction (Eq. 11) was applied, the resulting σ ap was not especially high and ω o was in the range of 0.98-1.00 for both the upper and lower estimates of σ sp , which indicates that non-scattering-corrected absorption coefficients may be considerably overestimated when σ sp is high.
The scattering coefficients calculated from the size distributions, averaged over the filter sampling periods, correlate positively with the scattering coefficients calculated from the PM 1 and PM 10 filters (Fig. 8). Ordinary least squares regression was used here. The main purpose of the regression was to study whether there is a statistically significant correlation between the scattering coefficients calculated from the size distributions and the filter samples. According to the slopes 0.78 ± 0.02 and 0.76 ± 0.04 of the regression lines in Fig. 8a and c, σ sp (DMPS, PM 0.8 ) also seems to be the lower estimate of σ sp when it is compared with the filter-samplederived σ sp . According to the slope of 1.29 ± 0.04 in Fig. 8b, σ sp (DMPS, PM 10 ) is an upper estimate of σ sp compared with σ sp (PM 10 ), but when σ sp (DMPS, PM 10 ), it is compared with σ sp (PM 1 ) and the slope is 1.01 ± 0.06, which appears to be somewhat controversial. There are also other peculiarities in the scatter plots. The scatter plots of σ sp (DMPS, PM 0.8 ) vs. σ sp (PM 10 ) (Fig. 8a) have data points where σ sp (DMPS) is low, in the range of ∼ 0.02-0.03 Mm −1 , but σ sp (PM 10 ) varies in a much larger range from ∼ 0.02 to ∼ 0.9 Mm −1 . This also occurs when σ sp (DMPS, PM 10 ) is compared with σ sp (PM 10 ) (Fig. 8b). The pattern could be explained by too low values of both σ sp (DMPS, PM 0.8 ) and σ sp (DMPS, PM 10 ) or by too high values of σ sp (PM 10 ). A similarly suspicious pattern is not observed in the comparison with the PM 1 filters ( Fig. 8c and d), suggesting the problem may be with σ sp (PM 10 ). It is clear that this is not a calibration of either the size-distribution-derived or filter-sample-derived σ sp , but the main message of the regressions is that the values are of the same order of magnitude and that there is a statistically significant positive correlation between them which increases confidence in the results. When the regressions are compared with each other, it has to be kept in mind that the sampling periods and the number of samples of the PM 1 and PM 10 data were not the same.
Other reasons for the wide scatter of the data points are the mass scattering efficiencies (MSEs) used for calculating scattering coefficients from the filter samples (see Sect. 2.5), uncertainties in ion analyses from the filters, and uncertainties in calculating scattering coefficients from the size distributions, especially the estimation of σ sp (DMPS,PM 10 ) from size distributions measured with the DMPS only. In spite of all these uncertainties, the statistical values (averages and percentiles of the cumulative distributions) of the scat-tering coefficients are reasonably similar. For instance, the medians of σ sp (PM 10 , λ = 550 nm), σ sp (PM 1 , λ = 550 nm), σ sp (DMPS, PM 0.8 , λ = 530 nm), and σ sp (DMPS, PM 10 , λ = 530 nm) were 0.24, 0.24, 0.15, and 0.23 Mm −1 , respectively (Table 4). The fact that the medians of σ sp (PM 10 ) and σ sp (PM 1 ) are the same is somewhat suspicious: it would be expected that σ sp (PM 1 ) < σ sp (PM 10 ). At this point it is worth paying attention to the statistics of the mass concentrations calculated from the size distributions and from the sum of ions in the filter samples (Table 5). The median mass concentrations of the PM 1 and PM 10 filters were 66 and 126 ng m −3 , respectively, in the expected order. These mass concentrations are also in reasonably good agreement with the median m(DMPS, PM 0.8 ) of 70 ng m −3 and median m(DMPS, PM 10 ) of 110 ng m −3 (Table 5). This suggests that the MSE values used for calculating scattering coefficients from the filter masses were not correct. As was written in Sect. 2.5, the MSE values were taken from Hand and Malm (2007), who derived them from measurements conducted mainly in US national parks. Considering this, the agreement of the filtersample-derived with size-distribution-derived σ sp is reasonable.

Seasonal cycles of scattering and absorption coefficients
The seasonal cycles of scattering and absorption coefficients are presented in Fig. 9. The SPO scattering coefficients presented in Fig. 9a-d were measured using the TSI nephelometer, and the Dome C scattering coefficients were calculated using the PM 1 (Fig. 9a) and PM 10 (Fig. 9b) filter sample data as explained in Sect. 2.5 and from the number size distributions ( Fig. 9c-d). The maximum and minimum monthly average and median scattering coefficients were observed in austral summer and winter, respectively. At the SPO the scattering coefficient was similar to that at Dome C in austral summer, but there was a large difference in austral winter. At the SPO the maximum monthly average scattering coefficients were observed in austral winter but at Dome C in austral summer. This suggests that in austral winter SPO is more influenced by sea-spray emissions than Dome C. However, even though the averages and medians are lower at Dome C, high scattering coefficients are also occasionally observed there in austral winter, as is shown by the 95th percentiles in Fig. 9c and d and above in the time series of winter 2011 (Fig. 6). The data do not explain the reasons for the difference between Dome C and the SPO in austral winter. It may either be due to different geographical locations, different size ranges measured by the instruments, or both. A hypothetical explanation for the difference between the scattering coefficients at the SPO and Dome C could be that in the very dry conditions the particles are not spherical. It is true that the shape of particles affects light scattering. However, it mainly affects the polarization of scattered light: spherical particles do not change the state of the polarization of scattered light, but nonspherical particles do. This is used, for example, in polarization lidars to discriminate ice crystals, dust particles, and droplets. However, integral photometric characteristics, such as extinction, scattering, and absorption cross sections and single-scattering albedo, do not depend significantly on particle shape, as is shown in chapter 10 of the textbook by Mischenko et al. (2002). Therefore, nonsphericity is not a likely explanation for the difference.
The minimum monthly means and medians of σ ap at Dome C were observed in austral autumn (MAM) and the maximum monthly means and medians in austral spring (SON), which is different than the seasonal cycle of σ sp (Fig. 9e and f, Tables S2 and S4). As a result, the seasonal cycle of the single-scattering albedo ω o is such that the darkest aerosol, i.e., the lowest ω o , is observed in September and October and the highest ω o in February and March ( Fig. 9g and h, Table S5). When the lower estimate for σ sp (i.e., σ sp (DMPS, PM 0.8 )) is used for the scattering correction (Eqs. 11 and 12), the October monthly medians of ω o are 0.862 and 0.868 when using the B1999 and V2010 algorithms, respectively, and when the upper estimate σ sp (DMPS, PM 10 ) is used for the scattering corrections, the October monthly medians of ω o are 0.911 and 0.916 when using the B1999 and V2010 algorithms, respectively (Table S5). The highest monthly median single-scattering albedos are ∼ 0.98 and > 0.99 with both algorithms when using the σ sp lower and upper estimates for the scattering corrections, respectively. These results show that when σ sp is not measured but calculated from the size distributions, the σ ap and ω o are clearly less sensitive to the selection of the algorithm (B1999 or V2010) than to the scattering coefficient used for the scattering correction. However, as was noted in Sect. 2.3.4, it is likely that σ ap (σ sp (DMPS,PM 10 )) is closer to the true absorption coefficient than σ ap (σ sp (DMPS,PM 0.8 )), so we can also consider the seasonal cycles presented in Fig. 9d, f, and h to be the closest to the true ones.

Seasonal cycles of scattering and absorption Ångström exponents
The wavelength dependency of both scattering and absorption have clear seasonal cycles. The average scattering Ångström exponent of particles in the DMPS size range, α sp (DMPS,PM 0.8 ), varies from ∼ 2.6 in austral summer (DJF) to ∼ 2.1 in austral winter (JJA), indicating that in austral summer the size distributions are dominated by smaller particles than in winter (Fig. 10a, Table S3). This cycle is much clearer when α sp is calculated from the upper estimate of scattering: average α sp (DMPS,PM 10 ) varies from ∼ 1.9 in austral summer to ∼ 0.8 in winter. The seasonal cycle of α sp (DMPS,PM 10 ) is actually strikingly similar to the seasonal cycle of α sp of σ sp measured at the SPO. This supports the use of the wavelength-dependent formula (Eq. 8) for calculating σ sp (DMPS,PM 10 , λ) from σ sp (DMPS,PM 0.8 , λ). The range of α sp is much larger at the SPO than at Dome C, however. The main reason is probably that when σ sp (DMPS,PM 10 , λ) was calculated with Eq. (8), only the monthly averages of f σ sp (DMPS,λ) (Eq. 7) were used, but the fσ sp (DMPS,λ) range is actually quite large (Fig. 3). The SPO values were calculated from direct PM 10 scattering measurements from a nephelometer. the absorption Ångström exponent α ap was calculated for the σ ap without scattering correction and by using the B1999 and V2010 algorithms with scattering corrected using σ sp = σ sp (DMPS,PM 10 ). In all of them the data with σ ap > 3δσ ap were used.
The absorption Ångström exponent α ap was calculated for the non-scattering-corrected absorption coefficient σ ap,nsc and for the scattering-corrected σ ap (σ sp (DMPS,PM 10 )) with the two algorithms. Close to the σ ap detection limit the ratios of σ ap at two wavelengths are very noisy, so Fig. 10b and Tables 4 and S6 present α ap statistics of absorption coefficients for σ ap > 3 × δσ ap , where δσ ap is the wavelength-dependent 24 h average noise at λ = 467 and λ = 660 nm ( Table 2). Note that the number of accepted data points is lower for the scattering-corrected than non-scattering-corrected α ap (Table 4). The reason is that the scattering correction often decreases σ ap below 3 × δσ ap .
The first observation that can be made from looking at the statistics (Fig. 10b, Tables 4 and S6) is that α ap (σ sp (DMPS,PM 10 ),V2010) is always larger than α ap (σ ap,nsc ) and α ap (σ sp (DMPS,PM 10 ),B1999). The main explanation for this is that the constants in the V2010 algorithm (Eq. 12) depend on wavelength, but the B1999 algo-rithm (Eq. 11) uses the same constants for all wavelengths. The differences between the α ap obtained from different algorithms were also discussed by Backman et al. (2014) and Luoma et al. (2021).
The seasonal cycles of α ap (σ ap,nsc ) and α ap (σ sp (DMPS,PM 10 ),B1999) are qualitatively similar: the lowest medians are observed in March and the maxima in August-October. This cycle is approximately anticorrelated with the ω o seasonal cycle: in March the median ω o is the highest, and the lowest is in August-October. In March the median α ap (σ ap,nsc ) and α ap (σ sp (DMPS,PM 10 ),B1999) were ∼ 0.6 and 0.37 and in August-September 0.96 and ∼ 0.92-0.96, respectively (Table S6), essentially the value generally used for pure BC. The seasonal cycle of α ap (σ sp (DMPS,PM 10 ),V2010) is a little bit different: the minimum median of ∼ 1.2 is in February, and the maximum of ∼ 1.7 occurs in June (Table S6).
The interpretation of α ap is complicated. The α ap is related to the dominant absorbing aerosol type, but physical properties of the particles also affect it. For externally mixed BC particles it is generally assumed to be around 1 (Hegg et al., 2002;Bond and Bergstrom, 2006;Bond et al., 2013) and higher for some organic aerosol from biomass smoke and mineral dust (Kirchstetter et al., 2004;Russell et al., 2010;Devi et al., 2016). However, α ap also depends on the size of BC cores and coating thickness. It is easy to show with Mie models that for single non-coated BC particles with D p <∼ 20 nm, α ap is indeed close to 1, but when D p ≈ 100 nm α ap ≈ 1.3 depending on the wavelength pair used for the calculation and < 1 when D p >∼ 150 nm. For BC particle size distributions the width and the dominant particle size affect α ap . Coating of BC cores affects α ap even more: when BC particles are coated either with a light-absorbing shell or even with a light-scattering shell, α ap can be clearly larger than 1 (e.g., Gyawali et al., 2009;Lack and Cappa, 2010;Virkkula, 2021). Core-shell simulations of size distributions of BC particles coated with a light-scattering shell show that α ap >∼ 1.4 for the wavelength pair of 470/950 nm could be obtained for BC particle size distributions when the shell volume fraction is >∼ 89 %-90 %, and the geometric mean diameter of the BC particles is in the range of ∼ 70-100 nm (Virkkula, 2021). Higher α ap would also be obtained by coating with a light-absorbing shell such as brown carbon. In the present work such α ap values were obtained for α ap (σ sp (DMPS,PM 10 ),V2010) for the wavelength pair 467/660 nm. So, if these values are closer to the truth, it seems that the BC particles that are observed at Dome C are thickly coated and their dominant particle size is <∼ 100 nm. On the other hand, if the average α ap ≈ 0.8 obtained for α ap (σ ap,nsc ) and α ap (σ sp (DMPS,PM 10 ),B1999) is closer to the truth, the core-shell simulation of Virkkula (2021) suggests that BC particle size distributions would be dominated by thinly coated particles in the size range > 100 nm.

Seasonal cycles of mass concentrations, eBC mass concentrations, and mass fractions
The seasonal cycles of the mass concentrations m(DMPS,PM 0.8 ) and m(DMPS,PM 10 ) (see Sect. 2.3.2), the mass concentrations of the PM 1 and PM 10 filter samples, the mass fraction of the sum of secondary sulfur ions, the eBC mass concentrations, and the eBC mass fractions f eBC(m(DMPS,PM 0.8 )) and f eBC(m(DMPS,PM 10 )) (Eqs. 16 and 17) are presented in Fig. 11 and in Tables S1, S7, S8, and S9. Some corresponding published Antarctic data are also plotted in Fig. 11 for comparison. The m(DMPS,PM 0.8 ) and the m(PM 1 ) are consistent with each other in that the minimum median mass concentrations are observed in May and June and the maximum medians in February. This cycle is very similar to that observed at the Norwegian station Troll in 2007-2011 (Fiebig et al., 2014). The monthly average volume concentrations of particles in the size range 33-830 nm in Fig. 9 of Fiebig et al. (2014) were digitized and multiplied by the same particle density ρ = 1.7 g cm −3 that was used for the Dome C data and plotted in Fig. 11a. The average (± standard deviation) of the ratio m(DMPS,PM 0.8 ,Dome C)/m(DMPS,Troll) of the monthly averages is ∼ 0.6 ± 0.2, i.e., about 40 % lower at Dome C. Fiebig et al. (2014) reasoned that the seasonal cycle of particles in the size range measured by the DMPS, i.e., m(DMPS), is controlled by photo-oxidation-limited aerosol formation. This is obviously true for Dome C also. In February, when the maximum monthly average PM 1 and PM 10 concentrations were observed, the contribution of the sum of secondary sulfur ions (nssSO 2− 4 + MSA − ) was also the highest (Fig. 11b): the average (± standard deviation) contributions to the sum of ions in the PM 1 and PM 10 filters was then ∼ 81 ± 12 % and ∼ 61 ± 23 %, respectively. The concentrations and the contributions of nssSO 2− 4 + MSA − were lowest in July, ∼ 9 ± 5 % and ∼ 5 ± 5 % for PM 1 and PM 10 , respectively.
The seasonal cycle of larger particles (m(PM 10 ) and m(DMPS,PM 10 )) is much weaker (Fig. 11a) than the m(PM 1 ) and m(DMPS,PM 0.8 ) cycle. The explanation is that the contribution of sea salt to aerosol mass is the highest in winter (Fig. 11b) and that a large fraction of sea-salt particles is in the supermicron size range, in line with other studies of aerosols at Dome C (e.g., Jourdain et al., 2008;Udisti et al., 2012;Legrand et al., 2017a, b). Note that the seasonal cycle of the mass fraction of secondary sulfur ions is qualitatively similar to the seasonal cycle of the scattering Ångström exponent α sp (Fig. 10a): both have the highest values in the austral summer and the lowest values in the austral winter. This is especially clear for the PM 10 filters. The small insert in Fig. 11b shows the scatter plot of the monthly average α sp (DMPS,PM 10 ) vs. (nssSO 2− 4 + MSA − )/PM 10 . The relationship is essentially linear, and the correlation coefficient is high, r 2 = 0.93. Since the usual interpretation of the size dependence of α sp is that it is inversely proportional to the dominating particle size, it indicates that when the mass fraction of secondary aerosol is the highest, the dominating particle size is the smallest. As such this is not a surprising observation, but it is an additional piece of information that links the chemical composition and aerosol optical properties.
The estimated m(PM 10 ) values are consistent with the concentrations measured gravimetrically by Annibaldi et al. (2011) in December 2005-January 2006. The average PM 10 mass concentration they obtained was 134 ± 12 ng m −3 at p = 1013 hPa and T = 298 K, which equals 146 ± 12 ng m −3 at p = 1013 hPa and T = 273 K used in the present paper. The average (and median) PM 10 mass concentrations in the present work were 167 ng m −3 (140 ng m −3 ) and 167 ng m −3 (143 ng m −3 ) in December and January, respectively (Table S1), in good agreement with the gravimetric measurement of Annibaldi et al. (2011) even though their measurements were not conducted in the same period as ours. In the austral summer the mass concentration calculated from the size distributions (m(DMPS,PM 0.8 ) and m(DMPS,PM 10 )) were ∼ 100 ng m −3 higher than the sums of ions in the PM 1 and PM 10 filters (Table S1). Part of the explanation could in principle be that the density 1.7 g cm −3 used for calculating mass concentrations from the size distributions was too high, but it cannot explain all of it. Another possible explanation is that there were organic compounds not observed with ion chromatography. Caiazzo et al. (2021) took filter samples at Dome C in a different period, December 2016-January 2018, and analyzed them for organic and elemental carbon with an OC/EC analyzer. The average OC concentration was 86 ± 29 ng m −3 , approximately the concentration difference between the size-distribution-derived ions and the sums of ions in the filter samples.
The eBC concentrations eBC(σ ap,nsc ) and eBC(σ sp (DMPS,PM 10 )) were calculated from Eq. (15). For the scattering-corrected σ ap , the two algorithms,  11) and (12), yielded essentially the same absorption coefficients at λ = 530 nm. Therefore, only one of them is shown in the seasonal cycle plot in Fig. 11c, but both are presented in Supplement Table S7. On the other hand, eBC(σ ap,nsc ) is also plotted to show how much the scattering correction affects the calculated eBC concentrations in different seasons. For comparison, published monthly median eBC seasonal cycles at three other Antarctic sites are plotted in Fig. 11c: at Neumayer, a coastal site in Queen Maud Land, using two methods, an Aethalometer, and a MAAP (Weller et al., 2013), at Syowa, another Queen Maud Land coastal site using an Aethalometer (Hara et al., 2019), and at the SPO using an Aethalometer (Sheridan et al., 2016). The maximum median eBC concentrations are observed in October-November at all sites. The maximum eBC in October-November is ∼ 3 ± 1 ng m −3 , quite similar at all sites. For eBC it appears that there is no significant difference between the coastal and plateau sites. The highest monthly median eBC concentrations are those measured with the MAAP at Neumayer in October but, for the same month, the median Aethalometer-derived eBC at Neumayer is the lowest. The lowest monthly median eBC concentrations are observed in April-May at Neumayer, the SPO, and Dome C and 3 months earlier in February at Syowa. The lowest monthly medians, ∼ 0.2 and ∼ 0.3 ng m −3 , were observed at Dome C and the SPO in May, respectively. The minima were higher at the coastal sites. Note, however, that the eBC concentrations measured with the Aethalometer in Fig. 11c were not corrected for scattering. This correction was done only for the PSAP data from Dome C and automatically for the MAAP data from Neumayer. After the corrections the Dome C monthly median eBC(σ sp (DMPS,PM 10 )) ranged from ∼ 0.2 in May to ∼ 3 ng m −3 in October-November, i.e., approximately by an order of magnitude and approximately the same as at the SPO. The range is smaller at the coastal sites. This might be due to not correcting for the scattering artifact even though the range of MAAP-derived eBC concentrations at Neumayer is also smaller than on the plateau sites.
The seasonal cycle of eBC is somewhat different from that of the mass concentration. Consequently, the minimum eBC mass fractions in both size ranges (f eBC(m(DMPS,PM 0.8 ) and f eBC(m(DMPS,PM 10 )), Eqs. (16) and (17), were in February-March and the maxima in August-October (Fig. 11d, Tables S8 and S9). The eBC mass fractions during this peak were actually quite high. In particular, if it is assumed that all eBC is in the size range measured with the DMPS, even for the scattering-corrected eBC monthly medians and averages of f eBC varied around 4 %-5 % and the 75th percentiles around 6 %-7 % by using both algorithms (Table S8). These are BC mass fractions typically observed in urban locations (e.g., Liu et al., 2014;Shen et al., 2018), in airborne measurements over Europe (McMeeking et al., 2010), and in biomass-burning plumes (Pratt et al., 2011), suggesting that in these periods a large fraction of aerosol was long-range-transported aerosol from other continents or highly processed air with larger, more scattering aerosol preferentially removed. The highest eBC monthly average and median mass concentrations were observed in November, but then f eBC was lower than its maximum. This can be explained by the increase in the amount of new, non-absorbing natural secondary particles and condensational growth of BC cores by compounds originating from the sea austral during spring and summer. Järvinen et al. (2013) classified new particle formation (NPF) events observed at Dome C, and the highest fraction of new particle formation events was in November, while in austral spring the particle growth rate was also the highest. The minimum f eBC monthly averages were <∼ 0.5 % and medians <∼ 0.3 % in February-March (Tables S8 and S9). This minimum also occurs simultaneously with the minimum eBC concentrations. This suggests that during this time of the year the amount of long-range transported aerosol from other continents is at a minimum at the same time when the biogenic aerosol production from the oceans is still high.
The seasonal cycles of single-scattering albedo (Fig. 9) and eBC mass fraction (Fig. 11d) are anticorrelated with each other. It is logical: the lower the f eBC is, the higher the fraction of scattering aerosol and ω o is. Their relationships can be used for assessing whether their observed seasonal cycles could be explained by internal mixing of BC particles and scattering components. Linear regressions of monthly average and median ω o vs. f eBC yield high correlation coefficients, but the regression lines would yield negative values at f eBC = 100 %. So, an exponential function of the form of ω o (f eBC) = ω o (f eBC=0)exp(−k· f eBC) was fitted with the data (Fig. 12). The correlation coefficients were slightly worse, ∼ 0.89 ± 0.01 for ω o vs. f eBC(m(DMPS,PM 0.8 ) (Fig. 12a), than the ∼ 0.98 ± 0.01 for ω o vs. f eBC(m(DMPS,PM 10 ) (Fig. 12b). If the fitted exponential functions were valid up to f eBC = 100 %, the ω o (f eBC(m(DMPS,PM 0.8 ) would predict that the average ω o ≈ 0.2 and the ω o (f eBC(m(DMPS,PM 10 ) would predict that ω o (f eBC = 100 %) ≈ 0.06. These are reasonable values for pure BC: it has been measured that for fresh pure BC ω 0 is approximately 0.2 ± 0.1 (e.g., Bond and Bergstrom, 2006;Mikhailov et al., 2006;Bond et al., 2013).
To evaluate whether these relationships could be explained by coated spherical BC particles, lognormal size distributions n(D p , GMD, GSD) were generated where GMD is the geometric mean diameter, GSD is the geometric standard deviation, and the D p range is 3 nm-10 µm. In the simulations all particles had a BC core and a scattering shell, and the core volume fraction was the same for all particles in the size distribution. The core-shell model N-Mie (Voshchinnikov and Mathis, 1999) that is based on a recursive algorithm of Wu and Wang (1991) was used for calculating scattering and absorption efficiencies Q s and Q a and scattering and absorption coefficients. See Virkkula (2021) for a detailed description of the simulations. The BC core volume fraction varied from 0.25 % to 4 % to cover the range of f eBC(m(DMPS,PM 10 )) shown in Fig. 12b. The simulated ω o vs. f eBC(PM 10 ) agrees best with the observation-derived values when GSD = 1.8 and GMD = 200 nm. Smaller GMD and GSD yield lower ω o vs. f eBC(PM 10 ). A rigorous error-minimizing procedure was not done since the main goal of this simulation was to show that core-shell geometry can explain the observed ω o vs. f eBC(PM 10 ) relationship. However, a deeper analysis of the modal structure of the size distributions was considered to be beyond the scope of the present paper. Note that it was not assumed above that all particles observed at Dome C had a BC core. The simulations did not use any of the measured number concentrations or absorption coefficients. They were pure simulations with generated size distributions and variable BC core and scattering shell volume fractions. They were conducted to study whether the observed ω o vs. f eBC(PM 10 ) relationship could be explained with size distributions of coated BC particles. However, considering that the BC particles observed at Dome C most probably have been transported from the surrounding continents for thousands of kilometers, it is very likely that most of the observed BC particles have been coated during aging processes.

Aerosol sources and transport and their seasonal cycles
Source areas were modeled by using FLEXPART as described in Sect. 2.7. The 60 000 50 d back trajectories were used for estimating vertical transport routes and the source areas of the observed aerosols. Sources of individual plumes are not analyzed: the main goal is to find explanations for the observed seasonal cycles of aerosol optical properties.

Vertical and horizontal transport
The average residence time of the trajectories in each grid cell in the altitude range 50-14 400 m a.g.l. south of the Equator is shown for the inert tracer and for the simulated BC aerosol tracer (Fig. 13). The average altitude of both the inert tracer and the BC tracer decreases sharply as the latitude approaches 75 • S, which shows that on average air masses descend over the high plateau. This is in line with several studies that show that stratospheric air is brought down into the lower troposphere by descending air masses above the Antarctic continent (Ito, 1989;James, 1989;Parish and Bromwich, 1991;King and Turner, 1997;Krinner and Genthon, 2003;Stohl and Sodemann, 2010). Even though the average altitude of the trajectories shows descending air masses above to Dome C for the inert tracer, the low-level residence times are high also at low levels (< 1000 m a.g.l.) at latitudes south of about 40 • S, indicating that it is not uncommon that air masses flow near the surface up to Dome C. When transport was modeled with FLEXPART using incloud and below-cloud scavenging of aged BC particles as described by Grythe et al. (2017), for the BC tracer the residence times are lower than for the inert tracer except near Dome C. Additionally, the average trajectory altitude is clearly different than for the inert tracer. Far from Antarctica the BC tracer residence times are clearly smaller than for the inert tracer. This is due to the removal of BC particles by in-cloud and below-cloud scavenging and dry deposition. As a result there would be essentially no BC particles arriving at Dome C from north of latitude 10 • S at altitudes < 1600 m, which is indicated by the respective area in Fig. 13b as just white. Consequently, those BC particles that do get lifted above the clouds increase the average altitude of the BC tracer back trajectories compared with the inert tracer average altitude, as is shown by the respective black and white lines in Fig. 13b. At latitudes ∼ 20-50 • S the BC tracer average altitude is approximately 1000 m higher than that of the inert tracer.
Between Fig. 13a and b there are grey text boxes that show the approximate latitude range of some major BC-emitting regions. Indonesia is north of 10 • S, so BC particles should rise to about 2 km altitude in order to get transported to Dome C. The Indonesian smoke plumes are generally lower than that (Tosca et al., 2011), so their contribution to eBC observed at Dome C is probably negligible. The other three main biomass-burning regions Africa, Australia, and Brazil are more to the south, and the smoke plume heights have been observed at clearly higher altitudes (Pereira et al., 2016;Rémy et al., 2017;Gonzalez-Alonso et al., 2019). Smoke plumes were recently observed at a height of 4 km over Brazil and even higher, at 6 km over Australia (Shikwambana and Kganyago, 2021). Comparison with Fig. 13b shows that from these altitudes it is likely that BC particles also reach Dome C. All this is in agreement with Fiebig et al. (2009), who established a source-receptor relationship between biomassburning events in central Brazil and the aerosol observed at the Norwegian station Troll.
Near Dome C and over all continental Antarctica the BC tracer average altitude is lower than that of the inert tracer. When BC has a lower average altitude than the tracer, it means that the air at that point has undergone more removal above than below the average altitude. This suggests that near Dome C most removal takes place in the clouds, higher up than the release point. This implies that there is more incloud removal than precipitation scavenging or dry deposition over the plateau.

Source areas of eBC and high and low ω o using footprint differences
Source areas were next estimated by using FLEXPART and calculating the footprint differences as described in Sect. 2.7.1. The relative difference (RD) (Eq. 19) of the emission sensitivities in the two lowest layers (< 1000 m a.g.l.) was calculated for the highest and lowest 10 % of eBC concentration and ω o . As was shown in Sect. 3.3, the highest and lowest ω o percentiles correspond to the lowest and highest eBC mass fraction percentiles.
The highest eBC concentrations have an increased signal from the surrounding continents Australia, Africa, and South America (Fig. 14a). However, high RD values are also associated with large areas over the Pacific Ocean, the Indian Ocean, and the Atlantic Ocean. These high RD values on the oceans may in principle be due to shipping emissions, but it is more likely that they are due to BC emissions from the continents where the emissions are considerably higher than from shipping. Smoke from South America is predominantly transported southeast into the South Atlantic (Freitas et al., 2005;Reid et al., 2009;Ulke et al., 2011) and further to the east. African smoke first gets transported westward to the Atlantic Ocean by easterly winds (Reid et al., 2009). Most of this smoke wraps around the South Atlantic Subtropical High, gets mixed with the South American smoke plume, and flows east in the subtropical jet region towards Australia and even further to the remote South Pacific (Singh et al., 2000;Staudt et al., 2002;Chatfield et al., 2002). During the transport to higher latitudes, the smoke typically rises to about 5-7 km (Singh et al., 2000;Staudt et al., 2002;Reid et al., 2009).
In addition to the eastward transport in the subtropical jet smoke outflow associated with the trade winds, the dominating easterlies from South America towards the west to the Pacific Ocean have been observed between the Equator and 20 • S (Freitas et al., 2005;Ulke et al., 2011;Bourgeois et al., 2015). Bourgeois et al. (2015) studied aerosol transport over the Andes from the Amazon Basin to the remote Pacific Ocean by using 6 years (2007-2012) of data from the CALIOP satellite instrument. They found that mainly biomass-burning particles emitted during the dry season in the Amazon Basin are lifted in significant amounts over the Andes. Aerosols reaching the top of the Andes, at altitudes typically between 4 and 5 km, are entrained into the free troposphere (FT) over the southeastern Pacific Ocean, and they can be traced on average over 4000 km away from the continent, indicating an aerosol residence time of 8-9 d in the FT over the Pacific Ocean.    RD,Eq. (19), of transport for the two lowest layers (> 1000 m a.g.l.) between the highest and lowest 10 % of (a) eBC concentration and (b) single-scattering albedo. Positive values indicate a relative increase in transport and negative values a decreased transport from a given area by the highest 10 % of measurements.
to negative values in the South Atlantic immediately to the east of the continent. Such a decrease is not in the GISS sensitivity simulations. So, the RD values over the western half of the South Atlantic are the most suspicious ones in Fig. 14a. They suggest that the lowest eBC concentrations observed at Dome C were associated with air masses from there. This is not reasonable since the studies cited above show that South American smoke plumes generally flow to the South Atlantic. This discrepancy remains unexplained.
The highest ω o are obviously related to SSA and biogenic secondary aerosol emissions from the Southern Ocean (Fig. 14b). The differences between the oceanic sectors are not big, but there are some: in the Atlantic and Pacific sectors the RD of ω o is clearly positive south of latitude 30 • S, but in the Indian Ocean sector RD is positive south of latitude ∼ 40 • S. This can be compared with the chlorophyll-a (Chl-a) concentrations in five oceanic sectors presented by Becagli et al. (2021). They determined the seasonal cycle of Chl-a concentrations in five sectors: Weddell Sea (WS, 60 • W-20 • E), Indian Ocean (IO, 20-90 • E), western Pacific Ocean (WP, 90-160 • E), Ross Sea (RS, 160-130 • W), and the combined Bellingshausen and Amundsen seas (BA, 130-60 • W). They found that the highest Chl-a concentrations were in the WS, RS, and BA sectors and the lowest ones in the IO and WP sectors. This is in line with the above-described geographical distribution of the RD of ω o in Fig. 14b, suggesting that the differences in the RD distributions are due to biogenic aerosols.

Seasonal cycles of air-mass transport from different regions
Seasonal cycles of meteorological processes and conditions also affect properties of aerosols observed at Dome C. The most evident is solar radiation that affects new particle formation (e.g., Ito, 1989;Fiebig et al., 2014;Järvinen et al., 2013), but important are also the seasonal cycles of transport from the surrounding continents and from the marine boundary layer of the Southern Ocean . The seasonal contributions of air masses from different regions were calculated as the fraction of time the FLEXPART trajectories spent over Antarctica and the surrounding oceans and continents (Fig. 15).
In the austral summer (DJF) the contribution of continental air other than Antarctica is low, and its minimum is in April. This is also true for marine air masses other than the Southern Ocean. In summer the contribution of Antarctic air is at its highest. Then in May the air-mass transport from other continents increases clearly, and it reaches a maximum in July, a smaller fraction in August, and a second maximum in September. The seasonal cycle of transport from all the surrounding continental areas is qualitatively similar. It can be explained by the seasonal cycle of the tropospheric polar vortex. Waugh et al. (2017) described the structures, seasonality, and dynamics of the stratospheric and tropospheric polar vortices and their connections to extreme events at Earth's surface. Estimating from their Fig. 2, the approximate edge of the tropospheric polar vortex of the Southern Hemisphere is at the latitudes of 50 • S in January and 30 • S in July. This means that in the austral winter the southern parts of South America, Africa, and Australia are within the tropospheric polar vortex but not in the austral summer. This is in line with Fig. 15. The next question is how well the seasonal cycles of the aerosol optical properties observed at Dome C follow the emissions of major absorbing and scattering aerosols in different regions.

Relationships of seasonal cycles of BC emissions
in the surrounding continents and eBC at Dome C Monthly BC emissions from wildfires and agricultural burning and other fires in South America, Africa, and Oceania (Australia, New Zealand, Melanesia, Micronesia, and Polynesia) in 2006-2012 were downloaded from GFED v3.1. The seasonal cycles are presented in Fig. 16a together with the monthly average eBC concentrations -corrected for scattering -observed at Dome C. It is obvious in Fig. 16a that the maximum of the sum of the emissions occurs 3 months earlier than the maximum eBC concentration at Dome C. The seasonal cycles of air mass residence time over each of the major source areas play a crucial role in explaining the relationship between eBC at Dome C and BC emissions. The BC emissions were multiplied by the fraction of the time (f (t)) the FLEXPART trajectories spent over each of the three major source areas shown in Fig. 15. These quantities (f (t) × BC emissions) will be called residence-timeweighted BC emissions, RTW BC emissions. They are plotted in Fig. 16b together with the normalized eBC concentrations. Linear regressions between the eBC concentration and RTW BC emissions were calculated, but only the respective squared correlation coefficients R 2 are shown. The correlation between eBC concentrations and RTW BC emissions is weak, R 2 < 0.3 (Fig. 16b). However, the main reason is that the residence times shown in Fig. 15 present the fraction of time that trajectories arriving at Dome C at any given month have spent over each of the regions before arrival at Dome C. So, next it was assumed that the transport time from each of the source areas to Dome C is 1 month, so f (t) of each month was multiplied by the BC emission 1 month earlier. This shifts the RTW BC emission peaks of all the source areas so that all R 2 increase somewhat (Fig. 16c). When it is assumed that the transport time is 2 months, R 2 increase clearly, 0.688 for Oceania and 0.665 for South America, and for Africa correlation still remains low (Fig. 16d). In October-November, when the eBC concentrations at Dome C are the largest, the 2-month-shifted RTW BC emissions from South America are approximately an order of magnitude larger than from Africa and Oceania, suggesting that South American BC emissions are the largest contributor to eBC at Dome C.
The highest correlation between the eBC concentrations and RTW BC emissions with the 2-month time shift suggests further that transport time is considerably longer than the 30 d Stohl and Sodemann (2010) obtained in a 5.5-year climatology of atmospheric transport into the Antarctic troposphere. The 2-month transport time naturally implies that BC aerosol residence time is also so long. Williams et al. (2002) estimated by modeling that in the upper atmosphere aerosols of 0.065 µm in size have residence times of approximately 1 month and can be transported on a hemispheric scale. Papastefanou (2006) determined residence times of tropospheric aerosols with different methods and obtained a much shorter mean value of 8 d. Kristiansen et al. (2012) estimated accumulation-mode aerosol removal times from the atmosphere using a global measurement data set collected over several months after the accident at the Fukushima Dai-ichi nuclear power plant in March 2011. After correction for radioactive decay, the 137 Cs/ 133 Xe ratios can be assumed to reflect the removal of aerosols by wet and dry deposition. Using this method, Kristiansen et al. (2012) obtained removal times for 137 Cs of 10.0-13.9 d. The simulations of Croft et al. (2014) with instantaneous injections into layers at 5 and 7 km yielded mean lifetimes of 14.7 and 21.1 d, respectively. At higher altitudes removal mechanisms are weaker, so lifetimes are longer. These are altitudes where biomass-burning smoke plumes have been observed, as discussed above. Note that these lifetimes are the e-folding lifetimes τ from C(t) = C(t 0 ) exp(−t/τ ). That means, for example, that if biomass-burning smoke where BC concentration is 1000 ng m −3 -similar to concentrations observed in a Brazilian biomass-burning measurement campaign (Brito et al., 2014) -rises to 5 km and gets transported for 60 d, the BC concentration decreases to 16.9 ng m −3 . With τ = 10 d the concentration would decrease to 2.5 ng m −3 . At Dome C the 25th to 75th percentile range of eBC concentration in November is ∼ 2 − 6 ng m −3 (Table S7). As the simple esti-mation shows, this range is consistent with BC emissions of about 1 µg m −3 and transport for 2 months.
The above discussion does not mean that pollution episodes from the surrounding continents could not be transported more quickly than in 2 months. For instance, Fiebig et al. (2009) calculated the source areas and transport routes of an elevated aerosol concentration event observed at the Norwegian Troll station in Queen Maud Land. The analysis showed that the origin of the aerosol was biomass burning in Brazil 11-12 d before the arrival at Troll. Hara et al. (2010) analyzed several haze events observed at Syowa, the Japanese Antarctic coastal station. They found that the traveling time from southern Africa to Syowa varied in the range of 4 to 8 d, even more quickly than to Troll in the episode analyzed by Fiebig et al. (2009). The explanation for the considerably longer transport time to Dome C, 2 months, is that it was obtained from the correlation of seasonal cycles of BC emissions and observed concentrations, so it is more like a statistical estimate of seasonal transport times. That suggests that on average the BC emitted in the biomass burning in the surrounding continents circulates the Southern Hemisphere and migrates slowly towards the Antarctic upper plateau. Figure 16 also shows that the observed concentration varies by an order of magnitude only, although the emissions vary by 2 orders of magnitude. This suggests that at the observed seasonal minimum eBC in March-April the sources are something else than wildfires and agricultural burning BC emissions in the GFED data. Possible sources can be other anthropogenic emissions, such as traffic including shipping, industrial emissions, and heating. Stohl and Sodemann (2010) estimated that ship emissions south of 60 • S account for half of the total BC concentrations in the lowest 1000 m of the atmosphere south of 70 • S in December. However, due to the particle removal processes as described in Sect. 3.4.1, the contributions of ship emissions to BC at Dome C are much lower. No quantitative estimations were calculated. Xu et al. (2021) recently presented a new global BC emission inventory for all countries in the world for the years 1960-2017. In an Excel file supplement, Xu et al. (2021) present the emissions classified into energy sector, industrial, residential, on-road motor vehicles, other anthropogenic, and wildfires. The sums of South American, Southern Hemisphere African, and Oceanian emissions in the period 2008-2013 and the respective contributions of wildfire emissions were calculated. They were ∼ 40 ± 8 %, ∼ 49 ± 2 %, and ∼ 59 ± 9 % for South American, Southern Hemisphere African, and Oceanian emissions, respectively. However, only wildfire-emitted BC-particles have the potential to rise high enough to avoid in-cloud or below-cloud scavenging during the transport, as was discussed above. SSA emissions calculated with the FLEX-SSA offline tool (Sect. 2.7.2) show that SSA emissions from the Southern Ocean peak in the austral winter (Fig. 17a). Sea ice and the glacier surfaces, especially in the coastal regions, are also important sources of sea-salt aerosols (e.g., Frey et al., 2020;Hara et al., 2020). However, the aim was to find a qualitative explanation for the high scattering coefficients observed in winter, especially at the SPO. The FLEX-SSA offline tool was used for that, and it yields a seasonal cycle that at least partially explains the observations. There was no tool for estimating seasonal cycles of sea-salt emissions from sea ice.
Secondary marine aerosols are formed from the oxidation of emission products of phytoplankton. Becagli et al. (2021) discussed the relationships between biogenic aerosols (nssSO 2− 4 and MSA − ), chlorophyll-a concentration, which is a proxy for phytoplankton, sea ice extent and area, source areas, and the transport routes to Dome C and the Southern Annular Mode (SAM). They showed that these relationships are different in early, middle, and late austral summer. Here a subset of the data is analyzed in a method that is complementary to that of Becagli et al. (2021). The phytoplankton concentration maximum is in summer and the minimum in winter in the satellite data analysis of Behrenfeld et al. (2016) (Fig. 17a). However, phytoplankton concentration alone does not explain the seasonal cycle of secondary aerosol scattering coefficients. The DMS oxidation is a photochemical process, so a simplified method was used for estimating solar radiation intensity. Global radiation intensity I T (W m −2 ) was calculated by using a clear-sky model for direct and diffuse insolation (Bird and Hulstrom, 1981). I T was calculated at the surface (p = 1000 mbar) for each hour of the year at the prime meridian (0 • ) at three latitudes, 50, 60, and 70 • S, corresponding to the central latitudes of Zone 1 (45-55 • S), Zone 2 (55-65 • S), and the polar zone, Zone 3 (65-75 • S). The monthly averages of I T at 10:00-14:00 on each day of a month were calculated. The reasoning for using this time is that new particle formation is typically, although not exclusively, a daytime phe- Figure 18. Normalized monthly average scattering coefficients calculated from the nss sulfate concentrations of the PM 1 filter samples σ sp (nssSO 2− 4 ,PM 1 ) at Dome C vs. normalized monthly averages of C phyto PAR in the polar latitude zone with and without a time shift. nomenon at numerous locations around the world, also in Antarctica (e.g., Weller et al., 2015;Kerminen et al., 2018;Kim et al., 2019;Brean et al., 2021). The monthly product, C phyto I T , was then used as a proxy of biogenic secondary aerosol emissions. Behrenfeld et al. (2016) also presented the photosynthetically active radiation (PAR) in the polar zone obtained from the MODIS Aqua ocean color sensor acquired from http://oceancolor.gsfc.nasa.gov/cms (last access: 10 April 2022). In addition to the product C phyto I T , the monthly product C phyto PAR was also calculated for Zone 3 and used as a proxy for biogenic secondary aerosol emissions. In Zone 3 the two emission proxies C phyto I T and C phyto PAR agree well (Fig. 17a).
The normalized seasonal cycles of the SSA flux and the secondary aerosol emission proxy C phyto PAR in the polar zone are plotted together with the normalized seasonal cycles of the scattering coefficients at Dome C (Fig. 17b) and at the SPO (Fig. 17c). At Dome C both σ sp (DMPS,PM 0.8 ) and σ sp (DMPS,PM 10 ) approximately follow the seasonal cycle of the secondary aerosol emissions in October-April. In the austral winter the effect of SSA is visible, especially in the high monthly σ sp averages. The medians are clearly lower, and they follow the secondary aerosol emission proxy even in May. Also, the normalized seasonal cycle of σ sp (nssSO 2− 4 ,PM 1 ), the scattering coefficient calculated from non-sea-salt sulfate concentrations in PM 1 filter samples, is presented in Fig. 17b. It also follows approximately the seasonal cycle of C phyto PAR for the whole year. However, there appears to be a time lag. C phyto PAR grows quickly in October, but σ sp (nssSO 2− 4 ,PM 1 ) in November. When C phyto PAR is shifted by 1 month, they agree clearly better. The linear regression of σ sp (nssSO 2− 4 ,PM 1 ) vs. C phyto PAR yields the correlation coefficient R 2 = 0.65 when there is no C phyto PAR time shift and R 2 = 0.92 when C phyto PAR is shifted by 1 month (Fig. 18). This suggests that on average it takes ap-proximately a month for the secondary aerosol to be formed, grown, get mixed in the upper atmospheric layers, and be transported to the upper plateau. As was discussed above, Williams et al. (2002) estimated by modeling that in the upper atmosphere aerosol residence time can be ∼ 1 month, and Croft et al. (2014) obtained e-folding lifetimes of ∼ 15 and ∼ 21 d at 5 and 7 km altitudes, respectively, so it is reasonable.
There is one thing that the time shift cannot explain. The maximum σ sp -and mass concentrations as discussed in Sect. 3.3 -was in February. This cannot be explained by the seasonal cycle of the biogenic secondary aerosol emission proxy: the maximum C phyto PAR was in December, so even a 1-month transport time does not make February the maximum. Above it was assumed that the proxy is simply the product C phyto PAR. However, the underlying hypothesis of this proxy is that the emissions of the precursor gases of new particle formation, mainly DMS, are linearly related to the phytoplankton mass concentrations. This is probably not the case, and there are other factors. Another explanation could be meteorological. Becagli et al. (2021) showed that the MSA and nssSO 2− 4 maxima in February could possibly be explained by faster transport from the surrounding oceans with air masses traveling at lower elevation than in December-January. Further analyses of that relationship are out of the scope of the present paper, however. For more details, see Becagli et al. (2021). Note that σ sp (nssSO 2− 4 , PM 1 ) correlates well also with C phyto I T of zones 1 and 2 shown in Fig. 17a, which means the secondary marine aerosol may have been formed and grown also further to the north than in the polar zone.
At the SPO, σ sp is more strongly affected by the SSA emissions than at Dome C, especially in the austral winter and spring (Fig. 17c). Note that at the SPO the seasonal cycle of the normalized median σ sp follows roughly the normalized SSA emission flux in July-December. The seasonal cycle of the normalized average σ sp is different because averages are affected more by individual high transport events and medians represent better the prevailing background aerosol of the season. As explained above, the difference between Dome C and the SPO may either be due to geographical locations, different size ranges measured by the instruments, or both.

Summary and conclusions
Aerosol optical properties have been measured at several Antarctic sites, but scattering and absorption data measured at Dome C were not examined in detail earlier. This work fills that gap using light absorption from a 3λ PSAP and lightscattering coefficients (σ sp ) calculated from particle number size distributions measured with a differential mobility particle sizer (DMPS) and an optical particle counter (OPC). Additionally, single-scattering albedo (ω o ), absorption Ångström exponent (α ap ), scattering Ångström expo-nent (α sp ), and equivalent black carbon (eBC) concentrations were calculated. The sources of the aerosol were estimated by calculating footprints with FLEXPART and by calculating seasonal cycles of transport of both scattering and absorbing particles from different source areas.
Aerosol light-scattering coefficients were calculated from the DMPS size distributions using two different refractive indices and the contribution of scattering due to particles larger than those measured with the DMPS to provide a lower and an upper estimate of σ sp . Light scattering was also estimated using mass scattering efficiencies in conjunction with mass concentrations obtained from PM 1 and PM 10 filter samples. The two most frequently used algorithms to calculate aerosol absorption coefficients (σ ap ) from PSAP measurements require scattering coefficients. Both algorithms were applied and both the upper and lower estimates of σ sp were used in order to provide understanding of the differences between the algorithms and the impact of the scattering adjustment. The absorption coefficient calculated using the σ sp upper estimate was considered to be the best σ ap estimate.
There were clear seasonal cycles of σ sp and σ ap at Dome C. The maximum and minimum of σ sp were observed in austral summer and winter, respectively. The Dome C scattering coefficients were also compared with σ sp measured with a nephelometer at the South Pole Station (SPO). At the SPO the scattering coefficient was similar to that measured at Dome C in austral summer, but there was a large difference in the austral winter. At the SPO the maximum monthly averages were observed in austral winter. This suggests that, in winter, the SPO is more influenced by sea-spray emissions than Dome C. At Dome C the σ ap exhibited a different seasonal cycle of than σ sp -the minimum σ ap was observed in the austral autumn and the maximum in spring. As a result, the seasonal cycle of the single-scattering albedo ω o is such that the darkest aerosol, i.e., the lowest ω o ≈ 0.91, is observed in October and the highest ω o > 0.99 in February and March.
The scattering Ångström exponent α sp calculated from the σ sp lower estimate varied from ∼ 2.6 in austral summer to ∼ 2.1 in austral winter, indicating that in austral summer the size distributions are dominated by smaller particles than in winter. For the σ sp upper estimate α sp varied from ∼ 1.9 in austral summer to ∼ 0.8 in winter. This seasonal cycle is quite similar to the seasonal cycle of α sp of σ sp measured at the SPO. The uncertainty of the absorption Ångström exponent α ap is high, particularly in the clean conditions existing in Antarctica. However, despite the high uncertainties the seasonal cycles of α ap , with and without the scattering correction, are qualitatively similar: the lowest monthly medians are observed in March, while the maxima occur in August-October. This cycle -even that of the non-scatteringcorrected α ap -is anticorrelated with the ω o seasonal cycle.
The eBC mass concentrations were compared with eBC measured at three other sites: the South Pole and two coastal sites Neumayer and Syowa. The maximum monthly median eBC concentrations are almost the same (∼ 3 ± 1 ng m −3 ) at all these sites in October-November. This suggests that, as far as eBC is concerned, there is no significant difference between the coastal and plateau sites. The seasonal cycle of eBC is slightly different from the mass concentration calculated from the number size distributions measured with the DMPS. Consequently, the seasonal cycle of the eBC mass fraction (f (eBC)) is such that f (eBC) minimum is in February-March and the maximum is in August-October, anticorrelating with ω o . The eBC mass fractions in this peak are actually quite high and vary around 4 %-5 % with the 75th percentiles > 6 %-7 %, only slightly depending on PSAP correction algorithm. These levels of eBC mass fractions are typically observed in polluted air, suggesting that in these periods a large fraction of aerosol is long-range transported aerosol from other continents.
Source areas were calculated with 50 d FLEXPART footprints. The relative differences of the footprints calculated for the highest and lowest 10 % of eBC concentrations and ω o showed that the highest eBC concentrations and the lowest ω o were associated with air masses coming from South America, Australia, and Africa. Vertical simulations that take BC particle removal processes into account show that there would be essentially no BC particles arriving at Dome C from north of latitude 10 • S at altitudes < 1600 m. The main biomass-burning regions Africa, Australia, and Brazil are more to the south, and their smoke plumes have been observed at higher altitudes than that, so they can get transported to Antarctica.
The seasonal cycle of BC emissions from wildfires and agricultural burning and other fires in South America, Africa, and Australia were calculated from data downloaded from the GFED. The maximum total emissions were in August-September, but the peak of monthly average eBC concentrations is observed in November, 2-3 months later, not only at Dome C, but also at the South Pole and the coastal stations. This is considerably longer than the 30 d presented in an earlier study. If this peak eBC concentration is really due to the peak emissions from the above-mentioned fires in the surrounding continents it means that the aerosol from these fires remains in air for several months and gets mixed essentially over the entire Southern Hemisphere. This does not mean that pollution episodes from the surrounding continents could not be transported faster than in 2 months. Earlier studies have clearly shown that biomass-burning smoke can get transported from South America or Africa to Antarctica in 1-2 weeks. The explanation for the considerably longer transport time to Dome C, 2 months, is that it was obtained from the correlation of seasonal cycles of BC emissions and observed concentrations, so it is more like a statistical estimate of seasonal transport times. That suggests that on average the BC emitted in the biomass burning in the surrounding continents circulates the Southern Hemisphere and migrates slowly towards the Antarctic upper plateau.
The seasonal contributions of air masses from different regions were calculated as the fraction of time the FLEXPART trajectories spent over Antarctica, the surrounding oceans and continents. The BC emissions were multiplied by the fraction of time the trajectories spent over each of the three major source areas. In October-November, when the eBC concentrations at Dome C are the largest the 2-month-shifted residence-time-weighted BC emissions from South America are approximately an order of magnitude larger than from Africa and Oceania, suggesting that South American BC emissions are the largest contributor to eBC at Dome C.
The seasonal cycles of scattering coefficient at Dome C and at the SPO were compared with the seasonal cycles of secondary and primary marine aerosol emissions. The seasonal cycles of SSA emissions were simulated with the FLEX-SSA offline tool. The seasonal cycles of biogenic secondary aerosols were estimated from monthly average phytoplankton biomass concentrations obtained from the CALIOP satellite sensor data. The correlation coefficients between scattering coefficients measured at Dome C and phytoplankton biomass concentrations and a biogenic secondary aerosol emission proxy are high. It may take a month for the biogenic aerosol to be formed and get transported from sea level to Dome C. The scattering coefficients measured at the SPO correlated much better with the SSA emission fluxes in the Southern Ocean than the scattering coefficients measured at Dome C. The difference between the scattering coefficients at these sites may either be due to geographical locations, different size ranges measured by the instruments or both. Data availability. The Dome C aerosol physical measurement data are openly available in the Zenodo data repository (https://doi.org/10.5281/zenodo.6354922, Virkkula, 2022). The SPO data are available at the EBAS database (http://ebas.nilu.no, last access: 10 April 2022, Sheridan, 2022). The Dome C aerosol chemical composition data will be available upon request by writing to Rita Traversi, Silvia Becagli or Mirko Severi.
Author contributions. AV contributed to designing the measurements of aerosol physical properties at Dome C, processing and analyzing data, and writing most of the manuscript with contributions from all the coauthors. HG did all the FLEXPART modeling and processed GFED data and wrote the associated text. JB wrote the codes for data processing. TP reviewed and commented on the manuscript. MB, CL, and AL were responsible for the installation and long-term operation of the aerosol instruments at Dome C. SB, RT, and MS were responsible for the sampling and analyses of the PM 1 and PM 10 filter samples. VV contributed to designing and administering the measurements at Dome C. PS and EA were respon-sible for the SPO aerosol measurements. All the authors reviewed and commented on the manuscript.