Study of Temporal Variations of Equivalent Black Carbon in a Coastal City in Northwest Spain Using an Atmospheric Aerosol Data Management Software

: Light-absorbing carbonaceous aerosols (including black carbon (BC)) pose serious health is-sues and play significant roles in atmospheric radiative properties. Two-year measurements (2015–2016) of aerosol light absorption, combined with measurements of sub-micrometric particles, were continu-ously conducted in A Coruña (northwest (NW) Spain) to determine their light absorption properties: absorption coefficients ( σ abs) and the absorption Ångström exponent (AAE). The mean and standard deviation of equivalent black carbon (eBC) during the period of study were 0.85 ± 0.83 µ g m − 3 , which are lower than other values measured in urban areas of Spain and Europe. High eBC concentrations found in winter are associated with an increase in emissions from anthropogenic sources in combination with lower mixing layer heights and frequent stagnant conditions. The pronounced diurnal variability suggests a strong influence from local sources. AAE had an average value of 1.26 ± 0.22 which implies that both fossil fuel combustion and biomass burning influenced optical aerosol properties. This also highlights biomass combustion in suburban areas, where the use of wood for domestic heating is encouraged, as an important source of eBC. All data treatment was gathered using SCALA © as atmospheric aerosol data management support software program.


Introduction
Atmospheric pollution supposes a risk factor for cancer and cardiovascular and respiratory diseases. However, atmospheric aerosol particles are poorly understood components of the atmosphere. They have a large diversity of sources, formation and transformation processes, which imply a large number of different aerosol species and, therefore, different properties. Hence, and taking into account their large variability, it is necessary to characterize them to know how they behave in the atmosphere in order to assess, for example, their effects on health and climate [1], and even in the degradation of building materials [2].
Black carbon (BC) is the fraction of the carbonaceous aerosol in the atmosphere that is characterized by its strong absorption of visible light [3] and also by its resistance to chemical transformation. BC is a good indicator of primary emissions and it is often used as an indicator for the efficiency of abatement initiatives [4,5]. The importance of BC determination also relies on the fact that BC and organic carbon are the constituents of particulate matter that are most likely to cause adverse health effects [6]. In this manner, the identification Taking all the previous information into account, the aim of this study was (i) to evaluate eBC levels and their temporal variability in a residential area in the NW of Spain, (ii) to assess the previous information with data related to Nt, and (iii) to identify eBC primary local and regional sources of emissions from fossil fuel combustion and biomass burning. For this purpose, (i) we carried out a sampling campaign from January 2015 to December 2016, and (ii) we developed new software features/functionalities for the SCALA © software application, which was the atmospheric aerosol data management software used to carry out this work.

Measurement Site
The study here presented was carried out at the University Institute of Research in Environmental Studies of the University of A Coruña (43 • 20 11" north (N,) 08 • 21 7" west (W); Galicia, Spain). Specifically, the Institute is located in Oleiros (Figure 1), in the northwest part of the autonomous region of Galicia (NW of Spain). It is a commuter town with a low population density characterized by (i) medium traffic density, especially in summer due to tourism caused by its beach area [25], and (ii) by single family homes with heating systems based on wood combustion. The nearest populations are two small villages at 2 and 4 km distance. The closest large population is the city of A Coruña (around 250,000 inhabitants) at a distance of 10 km. There are industries close to the sampling site that can influence air quality: three power plants (10 km W, 25 km southwest (SW) and 60 km northeast (NE)), one factory for the production and transformation of metals (9 km W), an industrial area (petrochemical refinery, primary aluminum manufacturing and graphite electrode manufacturing) (10 km W), one solid waste incinerator (25 km SW), several port areas (2-11 km NW) and one airport (4 km SW). There are also other significant activities in a radius of about 15 km: paper manufacturing and processing, food and beverage plants, ports, hospitals, funerary homes, printers and laundries.
The west of Spain lacks studies on eBC, on eBC in combination with Nt, and on source apportionment of eBC. However, there is a traditional use of wood in this region as an energy source for domestic heating in single family homes outside urban areas. In fact, the Galician Autonomous Government (in the northwest (NW) of Spain) has encouraged the use of biomass for this end.
Taking all the previous information into account, the aim of this study was (i) to evaluate eBC levels and their temporal variability in a residential area in the NW of Spain, (ii) to assess the previous information with data related to Nt, and (iii) to identify eBC primary local and regional sources of emissions from fossil fuel combustion and biomass burning. For this purpose, (i) we carried out a sampling campaign from January 2015 to December 2016, and (ii) we developed new software features/functionalities for the SCALA © software application, which was the atmospheric aerosol data management software used to carry out this work.

Measurement Site
The study here presented was carried out at the University Institute of Research in Environmental Studies of the University of A Coruña (43°20′11" north (N,) 08°21′7" west (W); Galicia, Spain). Specifically, the Institute is located in Oleiros (Figure 1), in the northwest part of the autonomous region of Galicia (NW of Spain). It is a commuter town with a low population density characterized by (i) medium traffic density, especially in summer due to tourism caused by its beach area [25], and (ii) by single family homes with heating systems based on wood combustion. The nearest populations are two small villages at 2 and 4 km distance. The closest large population is the city of A Coruña (around 250,000 inhabitants) at a distance of 10 km. There are industries close to the sampling site that can influence air quality: three power plants (10 km W, 25 km southwest (SW) and 60 km northeast (NE)), one factory for the production and transformation of metals (9 km W), an industrial area (petrochemical refinery, primary aluminum manufacturing and graphite electrode manufacturing) (10 km W), one solid waste incinerator (25 km SW), several port areas (2-11 km NW) and one airport (4 km SW). There are also other significant activities in a radius of about 15 km: paper manufacturing and processing, food and beverage plants, ports, hospitals, funerary homes, printers and laundries.  The prevailing wind direction in Galicia is northeast (NE) in summer and southwest (SW) in winter. The local wind pattern in the sampling site is mainly conditioned by the land-sea breeze because of its proximity to the sea (1 km). Northwest (NW) winds are dominant, although other wind directions are also recorded, and generally carry in relatively clean air from the sea. However, the main industrial sources (including A Coruña city) are located in this direction.

Instrumentation
eBC mass concentrations were measured using one multi-wavelength aethalometer (Magee Scientific aethalometer model AE33, Aerosol d.o.o, Slovenia), with a BGI, MiniPM ® Inlet at a flow rate of 5 L min −1 . Data were recorded with a 1 min time resolution. Filterbased light attenuation by the deposited aerosol particles was measured at 7 wavelengths (370, 470, 520, 590, 660, 880 and 950 nm). The eBC mass concentration was calculated using the measurement at the 950 nm wavelength with a MAC of 7.19 m 2 g −1 [9].
The system used to carry out the measurements of UFP was one SMPS, operated in the scanning mode. The SMPS consists of an Electrostatic Classifier (Model 3080, TSI United Kingdom) and a DMA (Model 3081, TSI United Kingdom) connected to a Water Condensation Particle Counter (Model 3785, TSI United Kingdom). Only the data of Nt were used in this study. The sampled aerosol was dried using a Nafion ® dryer. Additionally, a pre-impactor with 0.0514 cm nozzle was used. The SMPS was set with a sheath and a polydisperse aerosol flow rate of 10 and 11pm, respectively, to scan the size range between 10 and 289 nm (although this range is not properly UFP (<100 nm); it is called this way in this paper since the range is closer to the ultrafine range than to the fine range (<1 µm)). Every 5 min, the system sampled two 120 s scans per sample. Data were compensated for losses by diffusion and multiple charges with the Aerosol Instrument Manager ® software (version 9.0.0., TSI, Inc., St Paul, MN, USA).
Quality control protocols were considered for both instruments, including a weekly check for verifying a correct operation and data acquisition. Sampling flows were checked on a monthly basis for SMPS, and on an annual basis for aethalometer. Moreover, several tests were performed annually for this last instrument (e.g., leakage tests). Finally, SMPS national intercomparison exercises were performed in the framework of the Spanish Network of Environmental DMAs (REDMAAS) [26] in order to evaluate the accuracy, since there are no standards to verify the correct measurement of Nt.
The sampling campaign considered in this paper was performed during the years 2015 and 2016, with data coverage close to 75% for the aethalometer and close to 55% for the SMPS.
The source apportionment of eBC was performed using the two-component model described by Sandradewi et al. [8], using light absorption measurements at 470 nm and 950 nm [27,28], since eBC from fossil fuel has a weak dependence on wavelength (i.e., AAE~1), whereas eBC from biomass burning features a stronger absorption spectral dependence and shows enhanced absorption at a shorter wavelength (i.e., AAE > 1) [29,30]. The absorption at 470 nm was used, instead of the UV channel at 370 nm, to minimize the interferences introduced by types of organic compounds, based on the sensitivity of the aethalometer model due to different wavelength combinations carried out by Zotter et al. [31]. Even though 880 nm is considered the standard channel for eBC measurement by aethalometers, the 950 nm wavelength was used in this research according to the results obtained in the sensitivity of the aethalometer model using different pairs of wavelengths carried out by Zotter et al. [31]. Therefore, the analysis here conducted was based on the wavelength range 470-950 nm. For the purpose of this research, the software features of the first complete functional version of the web based software application called WEP-PROACLIM, which is an atmospheric aerosol data management program presented by Andrade-Garda et al. [32] and later known as SCALA © (Sampling Campaigns for Aerosols in the Low Atmosphere) (https://proaclim.udc.es) [33], had to be extended. That is to say, new functionalities (i.e., a new software increment following the software development process presented by Andrade-Garda et al. in [33,34]) were developed. Next, we will present the first version of SCALA © , then the new increment will be addressed. SCALA © is a web-based software application that was incrementally developed in a multidisciplinary way to integrally support the documentation and the management and analysis of atmospheric aerosol data from sampling campaigns. Therefore, this software application allows the comprehensive management of the sampling campaigns' life cycle (i.e., management of the profiles and processes involved in the start-up, beginning, development and ending of a campaign) and provides support for both intra-and inter-campaign data analysis. Thus, for example, the campaigns involving the different groups belonging to the Spanish REDMAAS are currently managed through SCALA © .
The need for SCALA © arises from the fact that such campaign management was traditionally approached in an eminently manual way. Thus, for example, date(s) and place(s) for a campaign were agreed upon after a huge exchange of mails and, or, telephone calls (or using Doodle in the best case), and data processing and analysis were performed by integrating and managing (through Dropbox or Google Drive in the best case) a handmade spreadsheet created by the environment technicians. Obviously, a great effort of format standardization was previously required for this integration. After that, to proceed with the data analysis, it was also necessary to have a skilled (error-prone) handling of, for example, a spreadsheet like Excel and/or a programming language like R [35]. In brief, there was no software application supporting, in a holistic way, the work related to all the sampling campaign activities. Accordingly, SCALA © was developed as a holistic solution, avoiding arbitrary errors in the manual handling of data (thus achieving efficacy) and optimizing the effort (thus achieving efficiency). As an example of the latter, it should be pointed out that the effort required to prepare the data files for uploading to the ACTRIS (Aerosols, Clouds and Trace gases Research InfraStructure network; the reference network in Europe [36]) website could be up to a month of work the first time; due to its complex format and, mainly, to the ad hoc work needed depending on each measuring instrument considered. The standardization of the data files for SCALA © takes less than thirty minutes, and much less even if, for example, macros are used. Moreover, the call for campaigns, the work during the campaigns and the data analysis have been fully integrated into a unique software application. Thus, once the data has been uploaded, the data analysis of both intra-and inter-campaigns is available immediately in SCALA © . SCALA © integrates (i) a public section (in Spanish), to promote and inform about related research activities or projects in the Spanish R&D Plan; and (ii) a private functional section (in English or Spanish, depending on the registered user's preferences), for managing and analyzing data and documentation from campaigns. There are five user profiles: administrator (technical administrator of the software application), group manager, campaign manager, technical user and external/guest/anonymous user (i.e., a non-registered user, who can request access). Next, the three main profiles are briefly described along with the most relevant features of SCALA © . More details are described by Andrade-Garda et al. [33].

•
Group manager: Responsible for a group of environmental technicians. She/he manages the members and the measuring instruments of her/his group, votes in surveys about the campaigns the group has been invited to, and sets the configuration for the group equipment in campaigns.

•
Campaign manager: The group manager that registers a campaign in SCALA © and, therefore, is responsible for managing the campaign life cycle (i.e., invitations to groups for participation in the campaign, campaign surveys for date and place agreement, campaign opening and campaign closing). • Technical user: A member belonging to a group of environmental technicians. Its main features are the following: Data files, documentation and incident uploads, and intraand inter-campaign analysis (data files selection, charts selection, and exportation options selection).
With regards to the above-mentioned new software increment for SCALA © , required to meet the new management needs arising from the new research presented in this paper, it was developed by a subset of the original multidisciplinary team following the same software development model used for the first version. This new increment mainly considered the possibility of adding aethalometers into the measurement instruments set managed by a group manager, along with the corresponding new functionalities regarding data file uploads, graphical representations and statistical analysis. Basically, and for what concerns here, this new software increment focused on the technical integration of SCALA © with R in order to (i) embed new plots (percentile rose plots, conditional probability function plots, and bivariate polar plots) into SCALA © and (ii) represent eBC concentrations from the aethalometers jointly with data already in SCALA © (i.e., particle number size distributions from SMPSs). The statistical elaborations here presented were embedded into SCALA © thanks to the mentioned technical integration with R; specifically, with R's openair [37] package. Briefly, SCALA © server actually invokes (i.e., "exec" in Computing) R features (in the SCALA © server), and the results coming from R are embedded into the web pages served to the user (web browser). In this manner, the R programming and testing activities are completely hidden from the user (environmental technician). Moreover, the user does not need to have R installed in her/his computer (computationally speaking, the client). All this (i) avoids computing and testing activities by non-technical users, (ii) avoids handling errors, (iii) automates and integrates the analysis into a single environment and (iv) reduces the required client's computational power.

eBC Mass Concentrations
The mean eBC value ± standard deviation obtained at this site during the period of study was 0.85 ± 0.83 µg m −3 . The measured values were lower than those measured in other sites around Europe (Table 1). In Figure 2a,b a combination of several variations is shown: hourly, weekly and monthly plots of eBC concentrations divided by season. With regard to the hourly variations, two peaks were detected during the day: one in the early hours and another in the evening, with differences in the mean concentrations (values were higher in the evening peak during winter than in the other seasons). These differences are likely to be related to the combination of anthropogenic emissions and boundary layer dynamics. An increase in eBC concentrations during the early hours can be explained by the increase of anthropogenic emissions and due to the formation of a nocturnal boundary layer (NBL), which favors the accumulation of pollutants. A decrease in eBC concentrations during the middle hours of the day is likely to be associated with a reduction of these emissions (from domestic heating and traffic) and due to a gradual increase in boundary layer height and wind speed. This fact enhances the vertical and horizontal diffusion of pollutants. A similar diurnal variation in eBC concentrations were observed in Ostrava [40] and urban sites of Barcelona and Granada [14]. A marked difference between workdays and weekends was also observed, pointing to traffic as the most important emission source of eBC in this area. domestic heating and traffic) and due to a gradual increase in boundary layer height and wind speed. This fact enhances the vertical and horizontal diffusion of pollutants. A similar diurnal variation in eBC concentrations were observed in Ostrava [40] and urban sites of Barcelona and Granada [14]. A marked difference between workdays and weekends was also observed, pointing to traffic as the most important emission source of eBC in this area.  In order to study the possible influence of near emission sources, data were represented as a percentile rose plot of eBC concentrations (Figure 3). The percentile intervals are shaded and shown by wind direction. The plot shows the variation by season and whether it is night-time or daylight hours (note that the open-air package estimates sunrise and sunset to give either daylight or night-time. This estimation, valid for dates between 1901 and 2099, is made using the location, date, time and astronomical algorithms to estimate the relative positions of the Sun and the location on the Earth's surface. For the sampling site here considered, for example, daylight hours in January range from 8 a.m. to 5 p.m. UTC). eBC concentrations are higher on average in the autumn and winter and when the wind comes from the south. The reason is most likely due to the presence of single family homes with heating systems based on wood combustion in this wind direction. eBC concentrations are somewhat higher at night-time in spring compared to daylight hours and higher during daylight hours in summer when the wind comes from the south-east. In this case, the influence of road traffic could be the main source of this pollutant in this wind direction. to 5 p.m. UTC). eBC concentrations are higher on average in the autumn and winter and when the wind comes from the south. The reason is most likely due to the presence of single family homes with heating systems based on wood combustion in this wind direction. eBC concentrations are somewhat higher at night-time in spring compared to daylight hours and higher during daylight hours in summer when the wind comes from the south-east. In this case, the influence of road traffic could be the main source of this pollutant in this wind direction. The Conditional Probability Function (CPF) is a commonly used approach when identifying atmospheric pollutant sources, since it provides directional information concerning major sources (i.e., it shows which wind direction and wind speed intervals are dominated by high concentrations and give the probability of this occurring) [43,44]. Figure 4 shows a CPF plot for eBC over the entire sampling period. It is clear that the highest concentrations (95th percentile eBC concentration: 2.5 µg m −3 ) are dominated by south and easterly wind directions. In this case, these directions correspond respectively to (i) a group of single family homes and (ii) one of the most transited roads near to the sampling site. There are very low conditional probabilities of these concentrations being experienced for other wind directions. That is to say, apportions coming from the closest large population city (A Coruña) or coming from the industrial area above-commented are not relevant, as well as the ones coming from the sea (NW winds). Very similar results are The Conditional Probability Function (CPF) is a commonly used approach when identifying atmospheric pollutant sources, since it provides directional information concerning major sources (i.e., it shows which wind direction and wind speed intervals are dominated by high concentrations and give the probability of this occurring) [43,44]. Figure 4 shows a CPF plot for eBC over the entire sampling period. It is clear that the highest concentrations (95th percentile eBC concentration: 2.5 µg m −3 ) are dominated by south and easterly wind directions. In this case, these directions correspond respectively to (i) a group of single family homes and (ii) one of the most transited roads near to the sampling site. There are very low conditional probabilities of these concentrations being experienced for other wind directions. That is to say, apportions coming from the closest large population city (A Coruña) or coming from the industrial area above-commented are not relevant, as well as the ones coming from the sea (NW winds). Very similar results are observed when representing the CPF plot over the sampling data segmented by the seasons of the year.  With the AE33 aethalometer, ultra-violet light absorption can be assessed. Thus, the independent effect of background wood smoke particles can be identified where wood burning is a major source of personal exposure to carbonaceous particles [45]. Wang et al. [46] commented that organic components of wood smoke particles absorb light at 370 nm more effectively than 880 nm. As commented above, this enhanced absorption (ΔC = BC With the AE33 aethalometer, ultra-violet light absorption can be assessed. Thus, the independent effect of background wood smoke particles can be identified where wood burning is a major source of personal exposure to carbonaceous particles [45]. Wang et al. [46] commented that organic components of wood smoke particles absorb light at 370 nm more effectively than 880 nm. As commented above, this enhanced absorption (∆C = BC (370 nm)-BC (880 nm)) can serve as an indicator of residential wood smoke particles. In this paper, the enhanced absorption was calculated for absorption at 470 and 950 nm, respectively. As can be seen in Figure 5a,b, the ∆C median values were substantially higher during cold seasons compared to warmer seasons, as was also reported by Wang et al. [46] (these authors found a difference between winter and summer of ∆C = 0.18 µg m −3 (0.24-0.06 µg m −3 )). During summer, a high dispersion of concentrations was found, probably due to episodes of biomass burning. The diurnal pattern showed an evening peak probably due to the use of biomass for domestic heating in the single family homes near to the sampling site.

UFP
To compare the levels of UFP at this site with those in other cities, the Nt concentrations were computed obtaining the mean ± standard deviation of 3825 ± 4023 particles cm −3 for the entire sampling period. These concentrations were lower than those reported in the studies summarized in Table 2, probably due to (i) the different diameter range measured and, (ii) the different characteristics of the sampling site. It is important to note that in most cases, the sampling period was shorter than the sampling period in this study.

UFP
To compare the levels of UFP at this site with those in other cities, the Nt concentrations were computed obtaining the mean ± standard deviation of 3825 ± 4023 particles cm −3 for the entire sampling period. These concentrations were lower than those reported in the studies summarized in Table 2, probably due to (i) the different diameter range measured and, (ii) the different characteristics of the sampling site. It is important to note that in most cases, the sampling period was shorter than the sampling period in this study.  [15] Nt is about, in real number concentration, 4500-5000 particles cm −3 in autumn and winter, and decreases to 2800-4000 particles cm −3 in spring and summer. Nt shows clear diurnal cycles, similar to eBC concentrations (Figure 6a,b) in autumn and winter, which points to traffic as the major source of this pollutant. This variation was also observed in two campaigns across Spain during 2012 and 2016 in the REDMAAS framework, one of the sampling sites being the University Institute of Research in Environmental Studies [1]; and also in the urban background sites in Leicester [15]. It is important to highlight the peak at 12 in spring and summer, as also has been reported in other papers [48], when the high insolation (i.e., high solar radiation) favors new particle formation. The high peak of Nt between the hours from 0 to 6 a.m. in winter (Figure 5b) could be explained as follows: the nocturnal reduction of the mixing layer height, along with the increase in the concentration of some gases and the decrease in temperature, can produce this increase in Nt. and also in the urban background sites in Leicester [15]. It is important to highlight the peak at 12 in spring and summer, as also has been reported in other papers [48], when the high insolation (i.e., high solar radiation) favors new particle formation. The high peak of Nt between the hours from 0 to 6 a.m. in winter (Figure 5b) could be explained as follows: the nocturnal reduction of the mixing layer height, along with the increase in the concentration of some gases and the decrease in temperature, can produce this increase in Nt. With regard to Nt concentrations, they are higher on average in the autumn and winter and when the wind comes from the south and east (Figure 7). But in summer during daylight hours Nt increases with different wind directions (east, northeast and southwest). Additionally, but in less proportion, Nt increases in spring, both during daylight hours and night hours, when wind comes from the northwest (from A Coruña city and its industrial area). With regard to Nt concentrations, they are higher on average in the autumn and winter and when the wind comes from the south and east (Figure 7). But in summer during daylight hours Nt increases with different wind directions (east, northeast and southwest). Additionally, but in less proportion, Nt increases in spring, both during daylight hours and night hours, when wind comes from the northwest (from A Coruña city and its industrial area).
With regard to Nt concentrations, they are higher on average in the autumn and winter and when the wind comes from the south and east (Figure 7). But in summer during daylight hours Nt increases with different wind directions (east, northeast and southwest). Additionally, but in less proportion, Nt increases in spring, both during daylight hours and night hours, when wind comes from the northwest (from A Coruña city and its industrial area).

Source Apportionment of eBC
The measured diurnal and hourly average cycles of the AAE during the period of study are shown in Figure 8a,b. AAE presents a strong diurnal variation during winter and autumn, with values in the evening hours until dawn consistent with an increase of biomass burning activity due to the use of the domestic wood stoves.
AAE shows a clear seasonal cycle, with the largest values during winter and the lowest during summer. Large AAE values during winter are related to an increase in the contribution of UV-absorbing particles, likely BrC from biomass burning emissions [39]. Furthermore, the large variability observed in AAE values indicate changes in the contribution of different absorbing aerosol types. This variability could be due to the impact of different sources of absorbing aerosols (i.e., fossil fuel (ff) and biomass burning (bb)) with a distinct absorption spectral dependence. The AAE average value was 1.26 ± 0.22 which is similar to the one obtained in a suburban area of Helsinki (1.25 ± 0.11) [39] and was slightly lower than the AAE average values found in Rome (1.41 ± 0.21) [19]. Thus, this supports the hypothesis that the highest eBC levels observed during cold seasons are strongly affected by wood combustion sources. The lowest AAE values obtained in summer suggests a low contribution from non-BC absorbing components to the absorption process, as Nicolás et al. commented [49] for a remote site in Spain.
Applying the model used in Sandradewi et al. [8], it was found that fossil fuel combustion was responsible for 69% of the eBC in 2015 and 64% in 2016. These values are slightly higher than the ones obtained in winter in three sites of Granada (52-68%) [50]. When applying this model, it is necessary to take into account, and should thus be considered carefully, that eBCff comprises traffic emissions as well as carbonaceous aerosols originating from fuel oil and natural gas combustion but excludes coal-burning organic aerosols. Another limitation might be the presence of mineral dust particles (notably containing iron oxides) [51,52], which also absorb light at near-UV wavelengths [53].

Source Apportionment of eBC
The measured diurnal and hourly average cycles of the AAE during the period of study are shown in Figure 8a,b. AAE presents a strong diurnal variation during winter and autumn, with values in the evening hours until dawn consistent with an increase of biomass burning activity due to the use of the domestic wood stoves.  The average diurnal cycles of eBCff and eBCbb are summarized in Figure 9a,b. The eBCff is practically the same as that of the total measured eBC and higher than eBCbb concentrations. As expected, higher contributions to black carbon from biomass burning were observed in colder periods, which can be explained by the higher intensity of domestic heating. The decrease in eBCbb concentrations during the warmer seasons is likely explained by the decrease in local and/or regional emissions from domestic heating by wood stoves. This evidently affects the eBCbb concentration levels especially at this site, wherein the wood burning emissions are dominated by local appliances during winter. However, a great variation in eBCbb concentrations was found in summer probably due to episodes of biomass burning fires. The eBCff did not show a distinct seasonal dependency. Higher eBCff on weekdays during autumn and winter could be due to trips made by car to second homes when there is good weather (i.e., more trips appear to go to schools, to work, shopping, etc.).
During weekdays, the previous figure shows that the eBCff concentration was characterized by two peaks, one in the morning and another in the late afternoon/early evening, coinciding with peak traffic times. On the weekend, these peaks are less pronounced. A relatively weak morning peak and a strong evening peak characterize the average diurnal cycle of eBCbb. In winter, the contribution of eBCbb is higher during night-time, which can be explained by the typical operating time pattern of domestic heating appliances in the coldest months. Moreover, Figure 9b shows wider shades during the early morning hours and the late evening of summer, probably due to (i) episodes of thermal inversion, which retain pollutants during these hours, and (ii) traditional summer barbecues around the beach and in the single family homes starting at 5 a.m. (UTC). This increase in some summertime sources of biomass burning was previously observed in London by Fuller et al. [54]. In the following Figures 10 and 11, the data are split by day and night-time and th season of the year, and here the trends can be observed. In the case of eBCff, south-ea and north-west wind directions contribute most to mean concentrations in autumn an spring, both during the day and at night. This could be explained by the prevailing win direction in autumn and winter from SW, where there is a road and also an industrial are In the following Figures 10 and 11, the data are split by day and night-time and the season of the year, and here the trends can be observed. In the case of eBCff, south-east and north-west wind directions contribute most to mean concentrations in autumn and spring, both during the day and at night. This could be explained by the prevailing wind direction in autumn and winter from SW, where there is a road and also an industrial area (see Figure 1). In the same seasons, NW is the dominant wind direction, which could bring pollutants from the city of A Coruña and its industrial area. No land-sea breeze effect was observed.   Once the previous analysis had been carried out with the entire sampling data, the daily values of both eBC and Nt were represented (not shown here) and the days in which the highest values of the two pollutants occurred were chosen. Once this selection of days was made, they were compared with the "classification of days with Saharan dust intrusion and with intrusion of particulate matter from biomass combustion provided by the Ministry for Ecological Transition" (https://www.miteco.gob.es/), to see if they coincided with any of those long-distance transport events. Two events were detected, and they are commented next through Figures 12 and 13. The data related to both events were not excluded in the overall analysis since it was confirmed that they do not affect the analysis.
As shown in Figure 12, we found a strong episode with high levels of the two pollutants during 4 February 2016, which coincides with an intense fresh biomass burning episode. Once the previous analysis had been carried out with the entire sampling data, the daily values of both eBC and Nt were represented (not shown here) and the days in which the highest values of the two pollutants occurred were chosen. Once this selection of days was made, they were compared with the "classification of days with Saharan dust intrusion and with intrusion of particulate matter from biomass combustion provided by the Ministry for Ecological Transition" (https://www.miteco.gob.es/), to see if they coincided with any of those long-distance transport events. Two events were detected, and they are commented next through Figures 12 and 13. The data related to both events were not excluded in the overall analysis since it was confirmed that they do not affect the analysis.
As shown in Figure 12, we found a strong episode with high levels of the two pollutants during 4 February 2016, which coincides with an intense fresh biomass burning episode. Another episode with high concentrations of eBC and Nt was detected from 26 October to 4 November 2016. In this case, dust and extremely fresh biomass burning episodes during those days were identified. This was confirmed by the high hourly levels found for eBCff and eBCbb (marked in yellow in Figure 13) during this period. With regard to Nt, and as happened with the above-mentioned episode, almost every hourly value are below 20,000 particles cm −3 during this episode. Another episode with high concentrations of eBC and Nt was detected from 26 October to 4 November 2016. In this case, dust and extremely fresh biomass burning episodes during those days were identified. This was confirmed by the high hourly levels found for eBCff and eBCbb (marked in yellow in Figure 13) during this period. With regard to Nt, and as happened with the above-mentioned episode, almost every hourly value are below 20,000 particles cm −3 during this episode.

Conclusions
A multi-wavelength aethalometer was used to study mass equivalent black carbon (eBC) characteristics at a suburban site in A Coruna during a two-year period. The eBC mass concentrations presented concentrations typical of other suburban sites, and lower than other urban sites in Europe. They show seasonal variations, significantly higher in winter and autumn than in spring and summer. This behavior is the result of meteorolog-

Conclusions
A multi-wavelength aethalometer was used to study mass equivalent black carbon (eBC) characteristics at a suburban site in A Coruna during a two-year period. The eBC mass concentrations presented concentrations typical of other suburban sites, and lower than other urban sites in Europe. They show seasonal variations, significantly higher in winter and autumn than in spring and summer. This behavior is the result of meteorological conditions during the studied period but can be also linked to the seasonal variations in emissions, especially higher in winter due to biomass burning. AAE shows a clear seasonal cycle, with the largest values during winter and the lowest during summer. The use of the aethalometer in areas where biomass burning could be an important source of eBC is very appealing, as these instruments are easy to operate and are often already installed in order to measure black carbon concentrations. Nt shows clear diurnal cycles, similar to eBC concentrations. The interpretation of measurements regarding eBC and Nt in this paper was made through the SCALA © software, an atmospheric aerosol data management software that was successfully extended to carry out this kind of study. The biomass burning source was identified (%eBCbb = 34, two-year data average), but it is not clear whether this source corresponds to open-air biomass burning of agricultural waste or biomass burning from domestic heating or both. Both practices are common in this area. Further studies are needed in order to assess the effects of these pollutants. Additionally, future pollution management strategies should take into account that biomass burning can be an important aerosol source during winter and spring seasons, since its human health effects may be different from those caused by diesel particles.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.