Relationship between Indoor High Frequency Size Distribution of Ultrafine Particles and Their Metrics in a University Site

Exposure to ultrafine particles (UFPs size < 100 nm) in life and work environments can contribute to adverse health effects also in terms of health burden of related diseases over time. The choice of parameters which better characterize UFPs is challenging, due to their physical-chemical properties and their variable size. It is also strictly related to the availability of different instrumental techniques. In the present study we focus on real time high frequency (1 Hz) UFPs particle size distribution (PSD) and their relationship with total particle number concentration (TPNC) and mean particle diameter (Davg) as a contribution characterizing by size the human exposure to UFPs in an indoor site of the University of Rome “Sapienza” (Italy). Further considerations about UFPs contribution to nucleation mode (NM) and accumulation mode (AM) have been highlighted, also in order to investigate the contribution of polycyclic aromatic hydrocarbons (PAHs) surface-adsorbed on indoor air particles (pPAHs). High indoor TPNC values were registered during the rush hours (early morning and mid/late afternoon) according to the outdoor influences originated from anthropogenic activities. AM mainly contribute to the indoor TPNC during working days showing high correlation with pPAHs. These findings may provide useful indications in terms of occupational exposure to UFPs since there are many evidences that indoor exposures to such pollutants may be associated with adverse health effects also in working environments.


Introduction
Ambient particulate matter (PM) is considered one of the leading environmental health risk factors, associated with several health impacts worldwide [1]. Among the important and variable parameters of PM that may be linked with health consequences of exposures is the particle size [2]. Ambient aerosols generally include particles with aerodynamic diameter between 0.01 µm and 100 µm. Particles between 2.5 µm and 10 µm are defined as coarse particles, and particles between 2.5 µm and 0.1 µm, as fine particles, whereas particles with diameter less than 0.1 µm are defined as ultrafine particles (UFPs) or nanoparticles [3,4].
Many studies have shown that UFPs can contribute to adverse health effects in life and work environments [5][6][7][8][9][10][11], measurable also over time in terms of health burden of related diseases [12]. UFPs are very small compared with the cellular structures and so this may be important in the consequences they may cause to the lung [13]. A particular concern is their ability to reach the most distal lung regions (alveoli) and circumvent primary airway defenses. When inhaled, UFPs can pass through the respiratory tract with high efficiency down to the alveoli due to their small size. A fraction of UFPs penetrates the alveolar-capillary barrier and can thus be distributed throughout the body via the circulatory system [14]. Because of this property of UFPs, extra-pulmonary diseases related

UFPs Outdoor Sources
UFPs in indoor environments are influenced by both indoor and outdoor sources [17]. In this context, an overview of main outdoor source contributions is needed. With respect to their formation mechanism, UFPs are usually considered as belonging to two types. The smallest is the nucleation mode (NM, particle size < 30 nm), which includes UFPs newly formed through source condensation processes or atmospheric chemical reactions: usually formed by volatile precursors as exhaust dilutes and cools, they have short atmospheric life (minutes). Otherwise, particles with aerodynamic diameter greater than 30 nm tend to have longer atmospheric lifetimes. The higher number and larger surface area for a given mass concentration favor the occurrence of condensation and coagulation processes. In this way these particles grow into the accumulation mode (AM, typically in the size range 30 nm-1 µm, but may be larger under high relative humidity conditions) which contains primary and secondary particles with longer atmospheric life (days) [15,[18][19][20][21][22].
UFPs particles may be also classified by their source contribution in natural (e.g., marine aerosol, volcanic eruptions) [23,24] and anthropogenic, typically generated through combustion of biomass (e.g., tobacco smoking, wood burning, incense burning) or fossil fuels (e.g., coal, natural gas, diesel) [25]. In the last years other sources of anthropogenic UFPs-that may present also correlations with indoor environments-were introduced, referred as nanomaterials either deliberately produced (i.e., engineered nanomaterials) or not voluntarily derived by production processes (defined as incidental) [26][27][28].
Motor vehicles, especially those driven by diesel engines, have been indicated as a major source of ambient UFPs emission [15,29]. UFPs can be formed from the condensation in the gas phase of semi-volatile organic compounds (SVOCs) that bypass the filtration technology [30]. On-road submicron particle size distribution typically contains an AM with estimated number mean diameter of 100 nm, consisting of carbonaceous agglomerates, e.g., elemental carbon and organic carbon that includes polycyclic aromatic hydrocarbons (PAHs). Vehicular traffic is an important source of particle-bound PAHs (pPAHs) [31]. Diesel and spark ignition vehicles also contribute to the occurrence of NM and his formation is extremely dependent on the dilution conditions of the engine exhaust [32].
UFPs generated from combustion of gasoline and diesel differ in chemical composition from those generated by solid fuel combustion. Coal combustions can be considered as one of the main contribution to the anthropogenic budget of NM, in both rural and urban areas [33,34]. The exact size depends on the distance from the source and the available time for coagulation and growth, although they are well below the range in which direct effects of solar radiation through scattering and absorption are important (>300 nm) [35][36][37].
Other source of UFPs is represented by biomass burning, natural or accidental (e.g., wildfire burning) and anthropogenic (e.g., biomass fuels such as peat, trees, leaves, and grass). The smoldering phase of biomass burning is associated with long-term output of UFPs and particles in AM, showing a bimodal distribution, with the smallest mode at~10 nm and the dominant mode of the particle-size distribution was in the range of 29-52 nm [38,39].

UFPs Indoor Sources
If abundant literature is available concerning outdoor sources, indoor sources characterization is far less defined. Buildings provide moderate, yet incomplete protection against UFPs of outdoor origin. With regard to infiltration of ambient particles from outside, it is quantified by the infiltration factor, which, in turn, is a function of several parameters, such as building construction characteristics (e.g., cracks, internal sources), air exchange mechanisms, outdoor meteorological conditions and indoor air circulation [40,41].
Indoor UFPs derives also from the emissions of indoor sources [42][43][44] which include cooking, domestic heating, smoking [45], use of electric appliances and chemical products, abrasion of textile materials and presence of individuals as well as UFPs emitted from building and furnishings materials [41,46].
In working environments aforementioned sources are missing and those ascribable to the presence of people become of particular importance, such as particle resuspension, clothes abrasion and bioaerosol release (i.e., skin fragments) [47][48][49]. Consumer products such as cleaning agents, air fresheners and personal care products contain terpene species, such as limonene, α-pinene or α-terpinolene [50]. Emissions of such species and subsequent chemical reactions indoors can allow the formation of numerous multifunctional and sometimes harmful secondary pollutants [51] such as VOCs which have low enough vapour pressures such that they partition to the particle phase to form secondary organic aerosols (SOA) [52,53].
Laser printers are another of the major contributors of UFPs in office environments [54]. UFPs formation may happen through nucleation and condensation of VOCs or SVOCs released by the fuser, which comprises siloxanes and fluorinated compounds, or emitted from the chassis, which holds flame retardants, lubricants, and plasticizers, in a typical bimodal PSD, with a smaller mode <10 nm and a broader mode extending from~40 nm up to 100 nm [55].

UFPs Metrics
In this framework research questions are still open about UFPs metrics and indicators of interest that might be used to characterize potential emissions, air quality, exposure and health effects [15,20]. Exposure to UFPs is a big challenge, due to their physical-chemical properties and their variable size [4,56], which may influence the choice and availability of instrumental techniques. Major measurement parameters include: total particle number concentration (TPNC), particle size distribution (PSD, i.e., particle number concentration per each size bin), mass concentration and lung deposited surface area (LDSA). Usually, only some of these parameters are measured, and in addition not simultaneously [57]. A common measure of UFPs prevalence is the TPNC: UFPs dominate TPNC and have negligible particle mass concentrations. Further findings support the concept that LDSA is the dose measure that predicts pulmonary response, rather than mass [5,13].
In the present paper we focus on UFPs PSD dynamics and their relationship with TPNC and mean particle diameter (D avg ) as a contribution to characterize by size the human exposure to UFPs in a University classroom. In this view, we analyzed real time high frequency (1 Hz) data from a measurements campaign conducted in the winter period in a computer room (CR) at the University site of Rome "Sapienza" (Italy), to describe airborne UFPs by working and non-working days and by daily time slots. Different contributions of NM and AM are studied in detail to investigate the formation mechanisms of indoor UFPs. Further considerations about pPAHs concentration and PSD were proposed, in order to highlight the ultrafine fraction of such pollutants in indoor air.
The present study is a part of the INAIL project "VIEPI" (Integrated Evaluation of Indoor Particulate Exposure, activity no. 3 "Indoor and outdoor particle number concentration") [41] in which extensive experimental activities for chemical and biological characterization of PM, as well as investigation of infiltration/exfiltration mechanisms, were conducted with the final aim to increase knowledge on the relationships between meteorological parameters, indoor airflow, and PM concentration.

Measurement Strategy
High frequency (1 Hz) real time measurements were conducted in the CR of the Physics Department of Rome University "Sapienza" ("Enrico Fermi" building) from Saturday 23 November 2017 to Thursday 27 November 2017. The CR is located at the second floor, it has a volume of 450 m 3 and it includes about 50 workstations (each equipped with personal computer case, screen, keyboard and mouse) for a maximum capacity of 80 persons. Natural ventilation of the room is provided by the windows placed on the East side, while hot-water radiators are operating during the sampling period as heating appliances integrated by air conditioning devices, manually managed by the occupants (Figure 1). Laser printers are located in the internal hallway placed on the West side of the CR. Outside the CR faces East on the courtyard that, in south direction, is 50 m far away from a road characterized by high traffic volumes and vehicular emission. During the working days a succession of teaching activities, exercises and exams took place in the CR without interruption from 8 a.m. to 6 p.m. We estimated that the CR was occupied by about 42 students on average, with a turnover of students between different classes. Windows were generally opened for a few minutes at the end of each lesson [41].
The present study is a part of the INAIL project "VIEPI" (Integrated Evaluation of Indoor Particulate Exposure, activity no. 3 "Indoor and outdoor particle number concentration") [41] in which extensive experimental activities for chemical and biological characterization of PM, as well as investigation of infiltration/exfiltration mechanisms, were conducted with the final aim to increase knowledge on the relationships between meteorological parameters, indoor airflow, and PM concentration.

Measurement Strategy
High frequency (1 Hz) real time measurements were conducted in the CR of the Physics Department of Rome University "Sapienza" ("Enrico Fermi" building) from Saturday 23 November 2017 to Thursday 27 November 2017. The CR is located at the second floor, it has a volume of 450 m 3 and it includes about 50 workstations (each equipped with personal computer case, screen, keyboard and mouse) for a maximum capacity of 80 persons. Natural ventilation of the room is provided by the windows placed on the East side, while hot-water radiators are operating during the sampling period as heating appliances integrated by air conditioning devices, manually managed by the occupants (Figure 1). Laser printers are located in the internal hallway placed on the West side of the CR. Outside the CR faces East on the courtyard that, in south direction, is 50 m far away from a road characterized by high traffic volumes and vehicular emission. During the working days a succession of teaching activities, exercises and exams took place in the CR without interruption from 8 a.m. to 6 p.m. We estimated that the CR was occupied by about 42 students on average, with a turnover of students between different classes. Windows were generally opened for a few minutes at the end of each lesson [41]. In the present study data relating to non-working days (Saturday and Sunday) were grouped and indicated with P1, those relating to Monday are indicated with P2 and Tuesday, Wednesday and Thursday with P3. We reasonably assume that a transition of sources in indoor environments may happen on Monday, therefore the related data were processed separately.  In the present study data relating to non-working days (Saturday and Sunday) were grouped and indicated with P1, those relating to Monday are indicated with P2 and Tuesday, Wednesday and Thursday with P3. We reasonably assume that a transition of sources in indoor environments may happen on Monday, therefore the related data were processed separately.

Equipment
Two instruments are used in parallel in this study:

•
Fast mobility particle sizer (FMPS mod. 3091, TSI Inc., Shoreview, MN, USA)) that is based on the ionic mobility principle and it is used to obtain the normalized PSD dN/dlogDp (part/cm 3 ) in the size range 5.6-560 nm divided into 32 different channels and the TPNC, with 1 s time resolution (total flow 10 L/min). • PAS 2000 (EcoChemAnalytics, League City, TX, USA) to measure pPAHs surfaceadsorbed on carbon aerosol with aerodynamic diameter less than 1 µm, with a response time of 10 s in a measuring interval from 0 to 1000 ng/m 3 and a lower detection limit of 3 ng/m 3 (total flow 1 L/min).
Further specifications about such instruments and their affordability are reported elsewhere [59][60][61]. Both PAS 2000 and FMPS measurements occurred simultaneously at the same location inside the CR. The instrument probes were placed at about 1.5 m from the CR ground, far from the windows/doors, the heating sources, the workstations and blackboard (Figure 1).

Statistical Modelling Approach
TPNC (part/cm 3 ) is referred in the following, as the sum of particle number concentration (PNC) of all particles having size between 5.6 nm and 560 nm, while normalized PSD dN/dlogDp is the particle concentration of each of the 32 channels in the instrumental range 5.6-560 nm, generally used to compare data from instruments with different resolutions [62].
Based on PNC, the averaged diameter (D avg ) was calculated at each time (t) according to: where N j is the PNC related to the j th diameter d j , with j ranging from 5.6 nm to 560 nm (into the 32 different channels of FMPS). Principal component analysis (PCA) [63] was performed as qualitative tool on the original dataset composed of 288 rows (averaged data every half hour for six days) and 4 columns (TPNC, NM, AM and pPAHs), in order to explore independent factors (principal components) that will illustrate the variation of dependent variables such as pPAHs and TPNC during the six time slots (T1-T6) in the non-working and working days.
The relationship between high frequency distribution and metrics has been studied according to the NM and AM mode. Data analysis of the observed distributions based on conventional statistical approach and on the evaluation of high frequency distribution were performed using OriginPro ® 2020 (OriginLab Corp., Northampton, MA, USA). TPNC time series show everyday a similar trend that is particularly evident from Monday to Thursday: TPNC reach the highest values about at midnight, then they decrease to a relative minimum between 6 a.m. and 8 a.m., going up to a relative maximum between 8 a.m. and 12 p.m. and reaching another minimum value between 12 p.m. and 4 p.m. Then, the TPNC trend is increasing until it reaches the daily maximum (>20,000 part/cm 3 ) at about midnight. Thursday differs slightly from this trend showing the daily maximum value between 8 a.m. and 12 p.m.

TPNC and D avg Time Series
A similar trend is also appreciable during Saturday and Sunday even if the TPNC values are lower than the working days. The observed distribution during the whole campaign is shown in the contour plot of particle diameter and normalized PSD versus time (Figure 2b). In the red area it shows that, in the working days at about midnight and in the early morning hours, the PNC highest values range from 2000 part/cm 3 to 2600 part/cm 3 for the particle sizes in the interval 60-115 nm.
High TPNC values in CR are according to the outdoor during the rush hours (early morning and mid/late afternoon) [41,64,65]. In the winter period, outdoor TPNC peaks at rush hours (in the morning and in the evening) can be considered as the result of the motor vehicle emissions combined with the mixing layer height and low ambient temperature, which favors nucleation mechanisms, at least for particles with dimensions up to some tenths of a nanometer [66]. Furthermore the TPNC outdoor reflects the daily profiles of population mobility indicators already proposed by Gariazzo and Pelliccioni [67] for the city of Rome in their study on urban population mobility using mobile phone traffic data. A similar trend is also appreciable during Saturday and Sunday even if the TPNC values are lower than the working days. The observed distribution during the whole campaign is shown in the contour plot of particle diameter and normalized PSD versus time (Figure 2b). In the red area it shows that, in the working days at about midnight and in the early morning hours, the PNC highest values range from 2000 part/cm 3 to 2600 part/cm 3 for the particle sizes in the interval 60-115 nm.
High TPNC values in CR are according to the outdoor during the rush hours (early morning and mid/late afternoon) [41,64,65]. In the winter period, outdoor TPNC peaks at rush hours (in the morning and in the evening) can be considered as the result of the motor vehicle emissions combined with the mixing layer height and low ambient temperature, which favors nucleation mechanisms, at least for particles with dimensions up to some tenths of a nanometer [66]. Furthermore the TPNC outdoor reflects the daily profiles of population mobility indicators already proposed by Gariazzo and Pelliccioni [67] for the city of Rome in their study on urban population mobility using mobile phone traffic data.

Relationship between NM and AM and Their Contributions to UFPs
In Figure 3 PSD during working (P2, P3) and non-working days (P1) for the six 4-h time slots (T1-T6) are shown. The observed distributions are essentially composed by three modes. Since the peak of lower Aitken mode [19,21] at about 20 nm is negligible compared to the first one (centered at about 10 nm), two mono-modal distributions were used: the first distribution concerns particle diameters less than 30 nm and the second one refers to particle diameters greater than 30 nm, respectively referred to NM and AM as classified in literature [15,18,19].

Relationship between NM and AM and Their Contributions to UFPs
In Figure 3 PSD during working (P2, P3) and non-working days (P1) for the six 4-h time slots (T1-T6) are shown. The observed distributions are essentially composed by three modes. Since the peak of lower Aitken mode [19,21] at about 20 nm is negligible compared to the first one (centered at about 10 nm), two mono-modal distributions were used: the first distribution concerns particle diameters less than 30 nm and the second one refers to particle diameters greater than 30 nm, respectively referred to NM and AM as classified in literature [15,18,19]. In Table 1 are reported the PSD values for P1, P2 and P3 during the different time slots (T1-T6), for NM and AM distributions.   In Table 1 are reported the PSD values for P1, P2 and P3 during the different time slots  (T1-T6), for NM and AM distributions. The PSD mean values for NM distribution are 3700, 2855 and 3991 part/cm 3 for P1, P2 and P3 respectively. The standard deviation lowest value occurs for P1 (502 part/cm 3 ) indicating that P2 and P3 are more widespread than P1. It is worth noting that the coefficients of variation (CV) are respectively 0.14, 0.35 and 0.47 for P1, P2 and P3. These findings highlight that during the weekend the NM distributions are about the same and not dependent by the time slot, while during the working days higher variability occurred. For AM distributions, the PSD values are always greater than the NM mode. On Monday (P2) some specific features have to be mentioned, related to NM and AM distributions by time slot, highlighting a possible transition of sources in the indoor environment between non-working days (P1) and the other working days (P3).
To take into account the specific contribution (in percentage) of each particle size, PSD values were calculated (PSD % ) respect to the daily total concentration [68][69][70].
Starting from Figure 4, PSD % values for NM and AM distributions were reported in Table 2, where the sum of AM and NM values respectively for P1, P2 and P3, is equal to 100% (corresponding to the bi-modal distribution subtended area). Starting from Figure 4, PSD% values for NM and AM distributions were reported in Table 2, where the sum of AM and NM values respectively for P1, P2 and P3, is equal to 100% (corresponding to the bi-modal distribution subtended area).     The observed ratio between AM and NM (AM/NM) is 2.5 in the weekend (P1) and 6.2 in the working days (P2 and P3). Furthermore, high variability of PSD values are observed during the 24 h: AM/NM ratio in the working days nighttime (T1, T2 and T6) is 9.4 and in the working days daytime (T3, T4 and T5) is 2.8; similarly, in the weekend AM/NM ratio is 3.2 and 1.7 respectively for nighttime and daytime.
All the ratios indicate a strong differentiation between working and non-working days. The contribution of the AM is always higher than the NM, regardless of the time slots and the day of the week. In particular, the ratio AM/NM in the working days is about 3.0 and 1.6 times greater than the weekend, in the nighttime and in daytime respectively. Figure 5 shows these results for each time slots (T1-T6) and in P1, P2, P3. The high concentrations reported on working days are probably due to anthropogenic activities and they might favor the coagulation and aggregation phenomena of airborne particles in the AM that reaches its maximum value during the night.  NM and AM distributions are characterized by different and well distinguished average diameters (Table 3). In fact, the average peak diameter for NM in P1, P2 and P3 is 10.5 ± 0.1 nm and for AM is 81.4 ± 8.3 nm. The corresponding CV is 0.008 for NM and 0.102 for AM, thus indicating the greater variability of AM peak position than the NM. These data are consistent with the literature as NM can be mainly associated to fresh/primary particles (i.e., diesel exhaust particles from outdoor, heating or natural sources) while AM includes not only primary particles but also those that have aggregated or coming from different particulate sources (either indoor or outdoor) [32]. To sum up these findings confirm and support the previous TPNC and D avg time series results (Figure 2), showing that in the working days the great amount of particles in the nighttime refers to AM.
NM and AM distributions are characterized by different and well distinguished average diameters (Table 3). In fact, the average peak diameter for NM in P1, P2 and P3 is 10.5 ± 0.1 nm and for AM is 81.4 ± 8.3 nm. The corresponding CV is 0.008 for NM and 0.102 for AM, thus indicating the greater variability of AM peak position than the NM. These data are consistent with the literature as NM can be mainly associated to fresh/primary particles (i.e., diesel exhaust particles from outdoor, heating or natural sources) while AM includes not only primary particles but also those that have aggregated or coming from different particulate sources (either indoor or outdoor) [32].

pPAHs in the NM and AM
PCA results are shown in Figure 6. The first and the second principal components explain 71.53% and 26.31% of the total variance. PCA was carried out on data that were firstly unit variance scaled and mean centered. In the loading plot (Figure 6a) TPNC, pPAHs and AM are grouped together contributing to the same information (positively correlated) and their direction is coincident with the first principal component whereas NM is independent by other variables (direction orthogonal with respect to PNC, pPAHs and AM) and its direction is coincident with the second principal component. In the score plot ( Figure 6b) the data points referred to the weekend (empty circles), for every time slot, are grouped together in the opposite side of the first component (i.e., pPAHs concentration) whereas data points referred to the working days (full circles) mainly correlate with pPAHs concentration during the nighttime (blue tones) and with pPAHs and NM during the day time (red tones). Starting from the qualitative considerations allowed by the PCA, we proceed with the estimation of parameters that describe the correlations between the variables. Figure  7 shows TPNC and pPAHs time series, day by day; the correlation between pPAHs, TPNC, NM and AM distributions in the weekend and in the working days are shown respectively in (Figure 8a,b).   Starting from the qualitative considerations allowed by the PCA, we proceed with the estimation of parameters that describe the correlations between the variables. Figure 7 shows TPNC and pPAHs time series, day by day; the correlation between pPAHs, TPNC, NM and AM distributions in the weekend and in the working days are shown respectively in (Figure 8a,b).
Starting from the qualitative considerations allowed by the PCA, we proceed with the estimation of parameters that describe the correlations between the variables. Figure  7 shows TPNC and pPAHs time series, day by day; the correlation between pPAHs, TPNC, NM and AM distributions in the weekend and in the working days are shown respectively in (Figure 8a,b).

Figure 7.
Indoor time series of total particle number concentration (TPNC, black curve) and particle-bound polycyclic aromatic hydrocarbons (pPAHs, red curve). Correlation parameters are reported in Table 4 Table 4. Correlation parameters (slope, intercept and R 2 ) between particle-bound polycyclic aromatic hydrocarbons (pPAHs) and nucleation mode (NM), accumulation mode (AM), total particle number concentration (TPNC) respectively, in non-working and working days.

Non-Working Days
Working Days Correlation parameters are reported in Table 4. The linear correlations between pPAHs, TPNC and AM in the working days (pPAHs: mean = 79.98 ng/m 3 ; Std. Dev. = 26.96 ng/m 3 ; CV = 0.34) and in the weekend (pPAHs: mean = 30.77 ng/m 3 ; Std. Dev. = 21.67 ng/m 3 ; CV = 0.70) show about the same slope and pPAHs concentrations are in accordance with those reported by Srogi [71]. Table 4. Correlation parameters (slope, intercept and R 2 ) between particle-bound polycyclic aromatic hydrocarbons (pPAHs) and nucleation mode (NM), accumulation mode (AM), total particle number concentration (TPNC) respectively, in non-working and working days. These figures show that inside the CR particles bounded with PAHs are about 0.6% of TPNC. A non-linear correlation occurs between PAHs and NM. It is worth noting that PAHs can be likely bounded on particles having size diameter greater than 30 nm (AM).

Conclusions
In the present study we analyzed ultrafine PSD and their relationship with TPNC and D avg , as a contribution to characterize by size the indoor exposure to UFPs in a computer classroom at the site of University of Rome "Sapienza" (Italy). High frequency (1 Hz) data were collected in a weekly campaign including working and non-working days, for six different daily time slots, covering both night and daily time periods. UFPs contribution to NM and AM have been highlighted, including statistical analysis of pPAHs and TPNC correlations with aforementioned two size modes.
High TPNC values registered inside the CR during the rush hours (early morning and mid/late afternoon) are according to the outdoor influences originated from anthropogenic activities. This is confirmed by similar daily trends of TPNC observed in working and non-working days, even if the values in Saturday and Sunday were lower than in the working days.
It is noting that PSD analysis reveals a site-specific condition in which the AM contributions are always greater than the NM ones inside the CR, not already reported for other different work environments [69,72].
NM can be mainly associated to fresh/primary traffic exhaust particles, heating or natural sources from outdoor while AM includes both primary and aged particles coming from different indoor or outdoor sources [32,70].
During the weekend the NM distribution is about the same and not dependent by the time slot, while higher variability occurred during the working days in which the great amount of particles refers to AM, in particular during the nighttime.
AM mainly contribute to the TPNC inside the CR during working days showing high correlation with indoor pPAHs, probably associated to the accumulation of outdoor contributions seeing as inside the CR primary pPAHs sources are not usually present. These findings may have an impact in terms of occupational health since there are many evidences that indoor exposures to PAHs may be associated with adverse health effects also in working environments [73].
The analysis conducted on the high-frequency data of working/non-working days highlights the complexity of indoor dynamics, not only related to the presence/absence of internal work-related sources, but also influenced by the background characteristic of indoor environments and further investigation are needed in this view.
Finally, recent studies highlighted the important role of indoor air quality in SARS-CoV-2 transmission, highlighting that virus particles may have similar dynamics than other suspended UFPs (e.g., arising from road traffic, heating or other sources) [74][75][76]. In this context, although there is a need for more specific research exploring possible interactions between air pollutants in the ambient air and SARS-CoV-2 impact on human health, the proposed methodological approach may provide useful indications, with specific reference