Simultaneous Measurements of PM2.5- and PM10-bound Benzo(a)pyrene in a Coastal Urban Atmosphere in Poland: Seasonality of Dry Deposition Fluxes and Influence of Atmospheric Transport

Parallel 24-h measurements of PM2.5and PM10-bound benzo(a)pyrene (BaP) were conducted in a coastal urban region in northern Poland in 2019. The aim was to examine distribution profiles, dry deposition fluxes and perform a trajectory analysis to investigate potential source regions of pollutants. The mean ratio of PM2.5/PM10 was 0.70 ± 0.23, suggesting significant heterogeneity in PM composition during the study period, both in fine and coarse fraction. The measurements revealed a strong seasonal variation in the concentration and dry deposition fluxes of BaP. The peak BaP concentrations of 6.14 ng m–3 and 6.65 ng m–3, respectively in fine and coarse PM fractions were observed in December. They were directly attributed to local emission from major industrial (i.e., coal-fired power plants, refinery, solid waste burning, metal recycling plants) and residential sources (domestic heating units), and to the meteorological situation that led to the reduced dissipation of air pollutants and their transport outside the study domain. During the period from October to November, high Pb/Zn ratios in PM2.5 and PM10 were also observed, suggesting significant contribution of industrial emission (coal combustion) and local nonexhaust traffic emissions, including tyre abrasion in the case of Zn and the resuspension of road dust rich in Pb. The highest daily BaP deposition fluxes of 2652 ng m–2 (PM2.5) and 5561 ng m–2 (PM10) were reported in winter, while the lowest values of daily BaP deposition (2 ng m–2) were found during spring and summer measurements. Additionally, the three-dimensional FLEXTRA transport model was used to identify the impact of different sources on air masses over the study domain. FLEXTRA results showed that the air masses transported from distant S to W regions greatly affected the chemical composition of PM2.5 and PM10 during autumn and winter.


INTRODUCTION
The role of atmosphere in the distribution and transformation of polycyclic aromatic hydrocarbons (PAHs) is crucial. Polycyclic aromatic hydrocarbons are emitted mainly from anthropogenic and industrial sources, including energy production, incomplete combustion of fossil fuels, coal and oil, iron/steel manufacturing, cement production, non-ferrous metal smelting, municipal waste incineration, vehicular emission (tyre/brake wear), petroleum refining, coke production and lowtemperature processes such as wood burning (Slezakova et al., 2007;Martiellini et al., 2012;Pereira et al., 2017;Ma et al., 2018;Nguyen et al., 2018;Yu et al., 2018;Hsu et al., 2019). The parent PAHs, such as naphthalene (Nap), acenaphthene (Ace), acenaphthylene (Acy), phenanthrene (Phe), fluorine (Flu), anthracene (Ant), benz(a)anthracene (BaA), chrysene (Chry), pyrene (Pyr), current environmental models and interdisciplinary field experiments. Therefore, taking into consideration the above facts, in this study, we performed daily simultaneous measurements of PM2.5 and PM10, at the coastal research platform, to examine the chemistry of a target and carcinogenic PAH -the BaP, with a special emphasis on its transformation and transport in the urban coastal atmosphere during four seasons. This first and in-depth case study of the seasonal cycle of BaP in fine and coarse particulate matter over the southern Baltic Sea has been undertaken in response to current and emerging scientific problems such as identification of major PAH sources associated with shipping activities and typical urban facilities, establish toxicity of PM2.5bound BaP and to determine conditions that influence BaP distribution and deposition within a study domain. Moreover, in this study we compare our results with those from other worldwide locations to consolidate recent knowledge on the atmospheric BaP chemistry and to better understand the relationships between toxic compounds, air quality and coastal environment.

Details of the Experiment and Aerosol Collection
The 24-hour fine (PM2.5) and coarse (PM10) particulate matter samples for the comprehensive chemical analysis were collected from January/March to December 2019, at the coastal site in Gdynia (54°52'55.1''N,18°56'13.5''E), Poland. Gdynia is a city in Pomorskie Province. Its area is 135.14 km 2 and it has up to 247 478 inhabitants. Gdynia can be characterised by high industrial, urban and traffic emissions. There is one of the biggest international harbours with shipbuilding industry (3,678 of incoming ships in 2016, BIP Gdynia). This city is a well-situated receptor area, suitable to investigate local and regional urban plumes, aerosol properties and long-range transport of pollutants originating from the surrounding countries, e.g., Germany, Belarus, Ukraine, Lithuania, Russia, Czech Republic, Slovakia and from Scandinavian countries. This region provides an excellent opportunity to examine the impact of shipping emission on aerosol chemistry, air quality and atmospheric processes. There are several dominant inland anthropogenic sources and shipping lanes within the study domain. The Google map of the sampling site is shown in Fig. 1, with major local emission sources such as coal-fired power plants, high traffic roads, domestic and commercial sectors, municipal solid waste recycling units, petrochemical plants and the refinery. In the vicinity of the sampling site there is also a large port area for various types of ships, i.e., small container ships, passenger vessels, cargo ships, etc.
In this study, the sampling instruments were mounted on the roof of Gdynia Aquarium. This building is located less than 100 metres from the shoreline. The sampling site was comprehensively equipped with a temperature-controlled set of automatic (PNS-18T Comde Derenda GmbH, Germany) and manual (LVS 3.1 Comde Derenda GmbH, Germany) samplers with a flow rate of 2.30 m 3 h -1 . The samplers contained meteorological sensors to register parameters such as air temperature, relative humidity, pressure, wind speed and direction. All the scientific instruments were installed on a special platform, about 4 m above flat-roof buildings, within the area clear of any obstacles, buildings and trees, in accordance with the European requirements for aerosol monitoring stations. After each daily measurement, environmental samples were immediately transported to the laboratory for subsequent analytical procedure. In this study, a total of 566 PM10 and PM2.5 samples were obtained and examined for BaP distribution and seasonal variability. In addition, an impact of local/regional contaminant sources located within the study domain was assessed as well as biological risk related to atmospheric BaP over Gdynia. This experimental BaP-focused project was consistent with all criteria and requirements for basic research. The backward trajectory calculations were based on FLEXTRA air mass transport model. It was used to identify the impact of different sources on air masses over the study domain ( Fig. 1), which we further discuss in the following sections.
During the entire sampling period, fine particulate matter concentrations ranged from 1.72 µg m -3 (July 29 th ) to 167 µg m -3 (February 17 th ), with the mean ± standard deviation of 24.1 ± 22.5 µg m -3 . Daily PM10 mass concentrations ranged from 3.45 µg m -3 to 65.5 µg m -3 (Table 1). However, it is important to highlight that PM10 was sampled from March to December. In both cases, higher values of mass PM concentrations were observed for the cold study period as compared to the warm period. This finding was mainly due to significant contribution of local  high-temperature processes in residential (i.e., coal combustion, domestic heating units) and commercial sectors (coal-fired power plants) during the cold study period. The observations also confirmed that sampling days with high PM mass concentrations coincided with the specific meteorological situation (i.e., lower height of the mixing layer, lower temperature, higher relative humidity) and local chemistry-climate processes (NOx, SO2, CO). Meteorological conditions over the study domain exhibited a large variation during the whole study period (Table 1). Specifically, the monthly mean air temperature ranged between -2.0°C (January) and 20.5°C (July). Regarding the ambient relative humidity (Rh), its mean value was 65%, the monthly peak was observed in January (mean of 75%) whereas the minimum in June (mean of 60%). The southerly winds blowing from industrial areas were much frequent during winter months, while for spring and summer observations strong westerly and southwesterly winds were predominant. An important meteorological feature of the study domain is that coastal winds from the north-to-east sectors are registered during spring and summer months. Moreover, the haze-fog episodes are typical of spring months.

Quantitative Analysis of Benzo(a)pyrene in PM2.5 and PM10 Samples
The fine (PM2.5) and coarse (PM10) fractions of atmospheric particulate matter were collected on quartz microfiber filters (QMA, 47 mm in diameter, Whatman). In this study, the 3-step pretreatment procedure was adopted to each filter before sampling, i.e., it was baked in a muffle furnace at 550°C for 8h, conditioned in a borosilicate glass desiccator under constant temperature and humidity (24°C and 40%) for 24 h, then weighted at the microbalance with 10 µg precision. Immediately after field sampling, the filters were packed in aluminum foil, conditioned for 24 h in a glass desiccator, weighted and then folded again in aluminum foil and kept in a zipped bag at -20°C until the main analysis. All handling works with blanks, procedural and environmental filters were performed in a clean room equipped with HEPA units.
To determine the concentration of BaP in air samples, the quartz microfiber filters were extracted using Accelerated Solvent Extractor (ASE350 model, Dionex, UK). The accelerated extraction was applied under controlled pressure, temperature and time, which is recommended by U.S. EPA and described in method TO-13A. Such approach has been previously applied in many similar environmental studies . Here, dichloromethane and n-hexane in the ratio of 1:1 was used as a basic solvent. Briefly, 15 mL of DCM:n-hexane solvent (Merck) was added to the tightly-capped ASE cell containing one quarter of a filter. The extraction program was as follows: ASE350 system was heated to 100°C for 5 min, the pressure was maintained at 1500 psi and remained static for 15 min, then the system was cooled down to the lab temperature. Each filter was extracted 3 times in the same manner. After extraction, the raw aliquots were concentrated to about 1 mL in a rotatory evaporator at the bath temperature of 35°C. The sample extracts were then transferred to a pre-baked glass column filled with glass wood, 1 cm of anhydrous Na2SO4, 2 g of 2% activated silica gel and 1 cm of anhydrous Na2SO4 on top, and eluted with 20 mL of DCM and hexane mixture (1:1, v/v). Finally, the extracts were concentered via rotary evaporator in a vacuum to 1 mL and then gently evaporated to 0.5 mL under nitrogen stream in a 35°C water bath. After condensation, the extracts were transferred into glass chromatographic ampules and stored in a refrigerator until HPLC-FLD analysis.
Concentrations of benzo(a)pyrene in PM2.5 and PM10 were quantitatively determined by Shimadzu modular Prominence HPLC system (Japan) coupled with fluorescence detection (excitation, λ = 290 nm, emission, λ = 430 nm). Data were proceeded through the use of Shimadzu Image v.5 software. The 5-ring PAH congener was quantified based on a several-point standard calibration curve (0.05; 0.02; 0.03; 0.05; 0.10; 0.15 µg mL -1 ). The calibration curve was linear with a correlation coefficient (R 2 ) of 0.9999 for BaP. The PAHs standard solution in acetonitrile (100 µg BaP mL -1 , Chiron) was used to prepare sub-standards and determine the analysed congener in environmental air samples. The chromatographic separation of PAHs was obtained through the use of a Kinetex LC column (150 mm × 4.6 mm i.d., particle size of 3.5 µm, Phenomenex) and 50%:50% acetonitrile -deionized water as a mobile phase. The injection volume was 25 µL. The mobile phase flow was 0.5 mL min -1 . The method detection limit (MDL) for the measured BaP, expressed as three times the laboratory blanks, was on average 0.01 ng m -3 . In addition, the precision and accuracy of the analytical procedure were controlled using certified reference materials, i.e., ERM-CZ120 Fine dust. In general, 10 pre-cleaned QMA filters were loaded with 0.5 mg of CRM for PAHs quantification of recovery. The average recovery for the BaP analysis ranged between 88% and 105%. Field and procedural blanks were measured following the same procedure as in the case of environmental filters. The results from duplicate blanks represented less than 5% of the concentrations quantified in samples.

Calculation of BaP Dry Deposition Fluxes
The deposition fluxes of BaP (Fdry) for each sampling day (ng m -2 day -1 ) were calculated using the following equation: where CBaP is analyte concentration (ng m -3 ) in the particulate phase of PM2.5 or PM10, and Vd is deposition velocity of particles during a sampling day (m s -1 ). Here, it was assumed that Vd for fine particles was 0.5 cm s -1 (Rohde et al., 1980;Poor et al., 2004), while for coarse particles Vd was established at 1.0 cm s -1 . The monthly deposition fluxes of BaP were obtained by summing up the daily values, separately for PM2.5 and PM10.

Air Mass Trajectories
In order to investigate the sources of PAH congeners in PM2.5 and PM10 over Gdynia, the air mass trajectory analysis was performed. The air mass transport models have been widely used in many aerosol-oriented studies (Forster et al., 2001;Chang et al., 2011;Halse et al., 2011;Aalto et al., 2015;Moroni et al., 2015;ChooChuay et al., 2020). In the present study, 24-h backward trajectories (BTs) were computed using a three-dimensional kinematic FLEXTRA air mass trajectory model (NILU and Institute of Meteorology and Geophysics, Vienna) with numerical meteorological data provided by European Centre for Medium Range Weather Forecast (ECMWF). This Lagrangian particle dispersion model has been described by Stohl et al. (1995). Specifically, the meteorological data used for each FLEXTRA trajectory simulation was analysed with a spatial resolution of 1.25 degree and temporal resolution of 6 hours. Three-dimensional trajectories were then calculated using the vertical wind provided by ECMWF. Stohl and Seibert (1998) have pointed out that accuracy of the trajectories is typically of the order of 20% of the travel distance, however uncertainties can also be much larger for individual cases. In the present study, the BTs were calculated for the period of 12 th January-31 th December, 2019 (6:00 UTC). In total, 340 BTs were calculated for the study domain, at three different arrival heights: 500, 1000 and 1500 m, to elucidate conditions below and over the boundary layer. In this study, nine separate air mass trajectory sectors were designated over the study domain. The sectors were marked as: Nnorth, NE -north-east, E -east, SE -south-east, S -south, SW -south-west, W -west, NWnorth-west and L -local. Furthermore, the BTs and their contribution to the total number of air mass trajectories in each season were assigned to monthly deposition fluxes. The monthly mean deposition fluxes of PM2.5-and PM10-bound BaP were calculated for each season (spring, summer, autumn, winter) and discussed in Section 3.6.

Monthly Distribution of the PM2.5 to PM10 Ratio
As shown in Fig. 2, the mean mass ratios of PM2.5∕PM10 in Gdynia were significantly lower in June 2019 (0.35) than in the other months of the study period. It was found that the range of the PM2.5/PM10 ratio was quite similar in August and September. The mean values of PM2.5/PM10 observed in December and November were 0.78 and 0.92, respectively, and they were almost 1.4-1.6 times higher than the values reported in March.
It has been previously reported that lower PM2.5 to PM10 ratio may indicate significant contribution of local sources of coarse airborne particles originating from road and construction dust or sea salt aerosol and soil dust emissions (Slezakova et al., 2007;Upadhyay et al., 2011;Fang et al., 2017;Siudek and Frankowski, 2018). On the other hand, the opposite trend may be mostly attributed to larger input of the primary PM2.5 from vehicle emissions with a significant contribution from local coal burning, including power generation plants or residential heating (Xu et al., 2017). For example, Permadi et al. (2018) have found that the observed values of the PM2.5/PM10 ratio ranged between 0.74 and 0.83 at the mixed urban sites in southeast Asia, suggesting predominant influence of road traffic, residential cooking and other urban activities. Ma and Jia (2016) have reported the annual mean PM2.5/PM10 ratios for three megacities in China, i.e., Beijing, Shanghai and Guangzhou, to be 0.63 ± 0.17, 0.69 ± 0.13 and 0.67 ± 0.06, respectively, indicating that the populations of fine and coarse particles were more heterogeneous, and pointing to the large variability in contribution of primary and secondary PM sources. Hsu et al. (2019) have found the mean PM2.5/PM10 ratio of 0.54 ± 0.35 at a single site in Taipei city, indicating that the role of fine particle fraction was predominant in shaping the urban air quality of that region. In the present study, the mean ratio of PM2.5/PM10 was 0.70 ± 0.23, suggesting significant heterogeneity in the PM composition both for fine and coarse fractions. The analysis of the PM2.5/PM10 ratio in Gdynia also revealed a seasonal trend of particulate matter and differences in emission sources during the study period. Similar findings have been demonstrated by Xu et al. (2017). Table 2 shows temporal variations of BaP concentrations measured in PM2.5 (from January to December) and PM10 (from March to December) in Gdynia. Seasonally, the BaP concentrations in PM2.5 and PM10 were much higher during the cold study period (mean values of 1.59 ng m -3 and 1.857 ng m -3 , respectively) as compared to those registered in the warm study period (0.10 ng m -3 and 0.06 ng m -3 , respectively). The peak concentrations of BaP in fine (6.14 ng m -3 ) and coarse (6.65 ng m -3 ) PM fractions were observed in December 2019 (Table 2), i.e., in the period of heavy PAHs emission from local anthropogenic sources, i.e., coal-fired power plants, domestic heating units etc., and the specific meteorological situation which limited dissipation and transport of air pollutants outside the study domain. In this study, a simple statistical screening of possible relationships between selected variables was conducted. Correlations between BaP and meteorological parameters were different (Spearman coefficient, p < 0.05). In particular, BaP showed a higher negative correlation (-0.690) with temperature. This suggests that the ambient air temperature has significant impact on PM2.5-bound BaP processes. A more detailed analysis for BaP/T correlations in different seasons revealed the following trend: autumn (-0.706) > winter (-0.564) > spring (-0.432) > summer (0.164). Previous studies in urban environments have mentioned that low temperature in a cold study period (autumn-winter) together with greater emission of particulate-phase PAHs from the residential/commercial heating sector and stagnant conditions can significantly enhance high concentration of 5-ring congeners in PM2.5 and PM10, while the increase in temperature during the summer period favours vaporization of PAH from particle to gas phase and/or its rapid photochemical degradation Wang et al., 2017). In relation to air pressure, PM2.5-bound BaP showed relatively low correlations (-0.116) during the whole sampling period. In contrast, the correlation for BaP versus relative humidity was weaker and positive (0.336). This result suggests that relative humidity can be identified as an important factor for BaP processes in the atmosphere over the Baltic coastal region. The similar dependencies have been reported by Wang et al. (2016).

Seasonal Variation of the BaP Concentration in Gdynia
Our results are in general agreement with some previous studies carried out in Poland (Staniszewska et al., 2013;Wiśniewska et al., 2019), western Europe (Amodio et al., 2009;Pietrogrande et al., 2011;Martellini et al., 2012;Polachova et al., 2020) and Asia Hsu et al., 2019;Wang et al., 2016;Wang et al., 2017;Urbančok et al., 2017;Gong et al., 2018;Nguyen et al., 2018). Additionally, relatively high BaP concentrations in both PM phases during winter study period can be explained by frequent occurrence of a thermal inversion layer (Hsu et al., 2019) as well as low and/or below zero air temperature that directly favours rapid multiphase chemical reactions of organic compounds in the water phase of atmospheric particles, i.e., condensation and adsorption (Hsu et al., 2019). Other factors that can explain elevated BaP concentrations in ambient air during winter are: low concentrations of ozone, reduced photochemical activity (Pietrogranade et al., 2011) and high loadings of atmospheric particulate matter .
As can be seen in Table 2, relatively low and less variable concentrations of BaP were registered between May and September for PM2.5, and between March and September in the case of PM10. This can be directly linked to significant changes in the local emission of PAHs and increasing role of atmospheric decomposition of 5-ring congeners (i.e., photo-, thermal-and heterogeneous chemical reactions with reactive gases). Several studies have highlighted the significant influence of high ambient temperature and strong oxidation capacity coupled with low relative humidity on the reduction processes of atmospheric BaP in summer (Hsu et al., 2019). Furthermore, based on gas/particulate partitioning of PAHs in four seasons, Nguyen et al. (2018) have found that particulate 5-ring PAHs exhibited the lowest contribution in summer, whereas higher concentrations of particulate PAHs were observed in winter. Similar atmospheric transformations of PAHs could explain the observed BaP profiles in this study during the warm season.

Seasonal Variation of the BaP Deposition Fluxes in Gdynia
As shown in Fig. 3, high deposition fluxes of BaP in PM2.5 and PM10 were predominant in winter. Specifically, the monthly mean PM2.5-bound BaP deposition was as follows: 614 ng m -2 (Nov) > 459 ng m -2 (Dec) > 244 ng m -2 (Oct), suggesting the increase of BaP importance due to intensive PAHs emission from local combustion sources and the observed severe pollution events in December 2019. In comparison, during the cold study period, the median BaP deposition in PM10 ranged between 352 ng m -2 in October and 1100 ng m -2 in November (Fig. 3). The highest BaP deposition flux of 2652 ng m -2 day -1 in PM2.5 was found on 16 th December, while the lowest values (1 ng m -2 day -1 ) were observed in April and May.
As for PM10, the peak BaP deposition flux of 5561 ng m -2 day -1 was reported on 18 th December, while the minimum BaP deposition (2 ng m -2 day -1 ) was found in April and August. In Gdynia, the monthly mean BaP deposition fluxes in PM2.5 substantially decreased from 365 ng m -2 in January to 85 ng m -2 in April, indicating a decrease in contribution of particulate PAHs during the period from spring to autumn (Fig. 3). In addition, the seasonal variability of BaP deposition fluxes at the coastal site in Gdynia was likely affected by meteorological conditions during the measurements (i.e., lower marine/planetary boundary layer height in winter vs. summer, stagnant conditions with low wind speed and high humidity that enhanced accumulation of air pollutants, seasonal differences in formation/decomposition processes in the coastal atmosphere, air temperature and relative humidity).

Comparison of the BaP Concentrations and Deposition Fluxes in Gdynia and Other Sites
The mean PM10-and PM2.5-bound benzo(a)pyrene concentrations measured in Gdynia in 2019 were 0.50 and 0.64 ng m -3 , respectively. These BaP levels were slightly lower than those observed by Staniszewska et al. (2013) and Lewandowska et al. (2018), mostly due to lower frequency of severe pollution events during the heating season in 2019 as compared to previous PAHsoriented measurements in 2008 and 2012.
A more detailed comparison of regional variations in BaP concentrations and deposition fluxes is presented in Table 3. A review study by Ma et al. (2018) have shown that atmospheric PAH concentrations across the Asian domain are the highest of all types of the considered measurement sites and that carcinogenic benzo(a)pyrene caused serious environmental problems. For example, the urban sites in Asia showed very variable concentrations of BaP in the ambient air, from < 1 ng m -3 measured in Hiroshima, Japan (Kakimoto et al., 2002) to 16.6 ng m -3 in Kathmandu, Nepal , with a mean value of 4.39 ng m -3 in Chengdu, China (Shi et al., 2015), which was much higher than those obtained in suburban sites such as New Brunswick (Lee et al., 2006). Moreover, Niu et al. (2017) showed extremely high mean values of BaP in PM2.5 within Beijing-Tianjin-Hebei region (6.87-18.30 ng m -3 , Table 3). The mean BaP levels in other large urban/industrial regions of Asia, such as in Seoul (2.55 ng m -3 , Lee et al., 2006), were much higher compared to Longyearbyen in Norway (0.89 ng m -3 , Drotikova et al., 2020). It was also found that BaP concentrations in TPM at urban traffic/residential sites in Europe (Pietrogrande et al., 2011;Martellini et al., 2012;Siudek, 2018;Šišovic et al., 2012), North America  and Australia (Mishra et al., 2016) were much lower than in similar urban areas of India (Ray et al., 2017) and China (Shi et al., 2015;Ma et al., 2018). Such large spatial differences in the TPM-bound BaP concentrations are a combined effect, which can be linked to different sampling periods, impact of major local/regional PAH emission sources (increasing or decreasing trends, or absence of any trends), patterns in precipitation height, number of episodes with extremely high deposition fluxes of atmospheric PAH congeners.
As presented in Table 3, mean BaP concentrations in some rural locations of Asia, i.e., Ulsan, South Korea (0.16 ng m -3 , Nguyen et al., 2018), Changhua, Taiwan (0.187 ng m -3 , Chen et al., 2016), were lower than those at the coastal site in Gdynia. Season-depended features can significantly affect the annual and multi-year variability in BaP levels in a study region (Šišovic et al., 2012;Ma et al., 2018;Siudek, 2018). For example, Liu et al. (2017) have shown some significant temporal trends in atmospheric PAH concentrations and deposition values for a group of monitoring PAHs in the US. They have found that in the 1991-2005 period, the mean BaP concentration level in the urban ambient air was 0.35 ng m -3 . In addition, the areas in the U.K. included in a national program (1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005) showed annual mean BaP levels of 0.25 ng m -3 (Meijer et al., 2008), while almost 10 times higher annual PM-bound BaP concentrations were identified in the mainland of China, including large cities such as Beijing, Chengdu, Dalian, Gangzhou, Nanchang, Shihezi,  Kumning, Lanzhou, Lhasa, Harbin, Xi'an (Ma et al., 2018) and Korea (Park et al., 2002). Recent long-term measurements at a European background site in Košetnice (Czech Republic) showed a large range for BaP concentrations in total particulate matter (< LOQ-9.1 ng m -3 , Shahpoury et al., 2015), with a mean level of 0.4 ng m -3 , which was almost 4.5 times lower than in the city of Poznań, central Poland (Siudek, 2018).
The mean annual PM2.5-BaP deposition flux in Gdynia amounted to 294 ng m -2 , which was much lower than that observed in Izmit Bay (1500 ng m -2 month -1 , Pekey et al., 2007), Cardiff and Manchester (6600 and 9000 ng m -2 month -1 , respectively, Hallsal et al., 1997). Arellano et al. (2018) have shown that BaP deposition fluxes at the European high-altitude stations such as Lochnagar (UK), Pyrenees (Spain) and Tyrolean Alps (Austria) are relatively low and vary between 20.0 and 61.0 ng m -2 month -1 , while in the Tatra Mountains the BaP deposition flux does not exceed 280 ng m -2 month -1 . In contrast, Gigliotti et al. (2005) have reported that BaP deposition values at some coastal sites in the North America tended to be one order of magnitude lower than the values from the Western Mediterranean Sea (Lipiatou et al., 1997). The monthly BaP fluxes measured in the present study were relatively low compared to previous measurements by Lewandowska et al. (2018), largely due to some significant decrease in wintertime emission of PAHs from local sources in the sampling year of 2019.
According to the recent national inventory report (IIR, 2020), PAH emissions in Poland have shown a significant decreasing trend between 2010 (309.19 t) and 2018 (231.14 t). The largest national PAH source is fuel combustion in non-industrial combustion plants (mostly residential/commercial heating sector), contributing 81-82% to the total difference in the period of 2005-2018(IIR, 2020. It should be noted that a large decline of 24% in the emission from residential sector have been reported for the period of 1990-2018, mainly due to replacing low-efficiency coal-fired boilers for more efficient devices and the implementation of thermomodernisation strategy (IIR, 2020). The second important anthropogenic origin of atmospheric PAHs is associated with industrial processes related to metal industry and other solvent and product usage as the dominant subsectors, followed by waste incineration and open burning of waste (IIR, 2020). The fugitive emission from fuels and field burning of agricultural residues are minor contributors to the national PAH emission.
In relation to the study region, it has been established that the BaP emission from residential/ commercial sector (domestic heating, coal-fired power plants) amounted to 7 018.1 kg y -1 , which was 98.3% of the total anthropogenic emission at the provincial-scale in 2019 (Voivodeship raport, 2020). It is noteworthy that the BaP emission from industrial and transport sectors was 111.4 kg and 11.2 kg, respectively. The percentage of other sources was very low (< 0.01%, 0.27 kg y -1 ).

Comparison of the Pb/Zn Ratios in PM2.5 and PM10
The geochemical data from observations in Poland related to the variability in metalliferous particles in inhalable fraction are still limited and more results are needed to consider the multiyear trends in atmospheric trace metal contents and their bioavailability (Siudek, 2021). In Gdynia, the monthly-averaged Pb/Zn ratios both in PM2.5 and PM10 displayed a quite similar trend, indicating the consistency of their possible local/regional sources with industrial hotspots (Fig. 4).
For example, the highest monthly median Pb/Zn ratios of 0.37 and 0.35 for PM2.5 and PM10, respectively, were observed in December, while the lowest mean Pb/Zn ratios in the fine fraction were found in February (Fig. 4). It has to be noted that the Pb/Zn ratios in PM2.5 revealed larger variability of concentration ranges as compared to PM10, in particular for the period from September to December. This finding may suggest contribution from local non-exhaust traffic emissions, including tyre abrasion of Zn and resuspension of road dust rich in Pb. However, other anthropogenic sources (i.e., solid waste burning, metal recycling plants, industrial processes) are also possible. In Gdynia, the PM2.5 and PM10-bound Zn had different seasonal profile, suggesting the impact of local/regional sources related to natural soil resuspension and important role of industrial emission (coal-combustion). As presented in Fig. 1, the local Zn sources are associated with petrochemical refinery and plants, vehicular exhaust emission, and also with soil and road dust resuspension within the study domain and the surroundings. Amato et al. (2014) have found that urban roads are characterised by large content of crustal species such as Al, Ca, K, Ti, Fe and Mn with significant amount of Pb. They have also reported that the road dust emission reveals clear seasonality with maximum values in summer and minimum in winter (Amato et al., 2014). Additionally, Garg et al. (2000) have reported that brake wear dust of motor vehicles contains large amount of Pb. In Gdynia, the local vehicle-related sources seemed to be more evident with winds (< 2 m s -1 ) coming from the west and southwest.
Previous studies have highlighted that during winter measurements, local atmosphere can be more polluted and, consequently, the concentration of metallic elements such as Zn and Pb in PM samples can be significantly higher (Alastuey et al., 2016), reflecting greater influence of local anthropogenic sources (i.e., coal combustion in power plants, shipping emission, waste incineration) and industrial activities (i.e., metallurgical processes). In Gdynia, we observed that this effect, similarly as in other polluted regions, was enhanced by thermal inversion layer and thick atmospheric mixing layer present in the cold study period, and this mitigated the transport of polluted particles up to higher altitudes leading to the significant impact of urban plumes on the study domain.

Seasonal Variation of the BaP Deposition in Relation to Air Mass Transport Sectors
It is well known that both local and regional emission of air pollutants combined with their atmospheric transport play a crucial role in depositional processes . In this study, one of the priority issues was to examine the differences in dry deposition fluxes of BaP in relation to air mass trajectories. Table 4 shows that PM2.5 and PM10-bound BaP deposition fluxes over Gdynia vary seasonally depending on the air mass trajectory clusters, with maximum values reported in winter and minimum values in summer.
Specifically, in winter, the mean PM2.5-bound BaP deposition fluxes associated with southern trajectories (SE-S-SW) were much higher than those attributed to local and northern trajectories (NW-N-NE). The daily mean deposition fluxes of PM2.5-bound BaP from southern sectors during  (Table 4). As reported by Lewandowska et al. (2018), when winter air masses originated from distant southern areas of Poland and Europe, they could transport higher loads of PAH-rich particles and therefore the deposition fluxes of air pollutants in Gdynia could be higher. Furthermore, most recent PAHsoriented studies by Siudek (2018) have demonstrated a significant contribution of particulatephase pollutants originating from industrial/urban and non-local sources in central Poland during autumn and winter study periods. However, it should be highlighted that the domestic heating from local residential sector had a predominant role in the total budget of PM. Previous studies have reported that the emissions from coal and wood burning in domestic systems of urban areas are much higher in winter than in a warm study period due to higher supply for energy. Here, daily mean PM2.5-bound deposition fluxes in spring were highest for the SW sector (343.28 ng m -2 ) followed by N (272.76 ng m -2 ), S (226.95 ng m -2 ) and local clusters (223.15 ng m -2 ), while the lowest deposition fluxes were observed for the NE cluster (76.60 ng m -2 , Table 4). In contrast, daily mean PM10-bound BaP deposition fluxes of 590.66 ng m -2 were found for the local cluster. This was at least 2 times higher as compared to the other clusters (range: 76.80-255.07 ng m -2 ), indicating that there were significant differences in the deposition of BaP during spring measurements. The results of PM2.5-bound and PM10-bound BaP showed that the summer deposition fluxes were below 60 ng m -2 day -1 for all clusters, and the highest values were associated with the local cluster (local traffic emission, shipping emission). In comparison with summer measurements, the values of BaP deposition fluxes in autumn months were much higher for both size-mode particles. Specifically, the highest mean daily PM10-bound BaP deposition flux of 1229.96 ng m -2 was measured for the S cluster, while the highest deposition of BaP in PM2.5 (528.35 ng m -2 ) was found for the SE cluster. As mentioned in some previous studies, due to higher air temperatures in summer, the condensation of most volatile PAHs can be significantly reduced and, as a consequence, the deposition fluxes are typically lower compared to other seasons (Arellano et al., 2018). Some similar trend was observed in this study.

CONCLUSIONS
Simultaneous measurements of PM2.5 and PM10 samples in a coastal urban region of northern Poland were performed for the first time in 2019 to determine the seasonality of PM2.5 and PM10bound benzo(a)pyrene, the target congener of five-ring PAHs. The concentrations of BaP were quantitatively determined by HPLC-FLD. The seasonal variability of BaP concentrations and deposition fluxes was evident. The highest values were observed in winter, being the result of large emission from local residential sector (domestic heating units) and major industrial sources (i.e., coal-fired power plants, refinery, solid waste burning, metal recycling plants) as well as meteorological situation leading to the reduced dissipation of air pollutants and their limited transport outside the study domain. In contrast, summer measurements exhibited considerably lower BaP concentrations and deposition fluxes in particulate matter over this coastal urban region, which was likely related to a decrease in local coal combustion (commercial and residential), vaporization of PAHs from particle to gas phase and/or its rapid photochemical degradation. It was also observed that the local traffic emission played a significant role. The large-scale atmospheric transport of air pollutants from distant industrially-impacted regions may influence the coastal aerosol population in the southern Baltic Sea, especially in winter.
The limitation of this study is that we have not addressed the influence of precipitation amount on the deposition fluxes of BaP during the sampling period, as considered by Staniszewska et al. (2013). Due to lack of routine PM10 measurements in January and February 2019, the obtained values of the dry deposition flux in coarse fraction for winter 2019, in all selected sectors, can be underestimated.