An investigation into the use of < 38 µm fraction as a proxy for < 10 µm road dust particles.

It is well documented that a large portion of urban particulate matters is derived from road dust. Isolating particles of RD which are small enough to be inhaled, however, is a difficult process. In this study, it is shown for the first time that the < 38 µm fraction of road dust particles can be used as a proxy for road dust particles < 10 µm in bioaccessibility studies. This study probed similarities between the < 10 and < 38µm fractions of urban road dust to show that the larger of the two can be used for analysis for which larger sample masses are required, as is the case with in vitro analysis. Road dust, initially segregated to size < 38 µm using sieves, was again size segregated to < 10 µm using water deposition. Both the original < 38 µm and the separated < 10 µm fractions were then subject to single particle analysis by SEM-EDX and bulk analysis by ICP-OES for its elemental composition. Dissolution tests in artificial lysosomal fluid, representative of lung fluid, were carried out on both samples to determine % bioaccessibility of selected potentially harmful elements and thus probe similarities/differences in in vitro behaviour between the two fractions. The separation technique achieved 94.3% of particles < 10 µm in terms of number of particles (the original sample contained 90.4% as determined by SEM-EDX). Acid-soluble metal concentration results indicated differences between the samples. However, when manipulated to negate the input of Si, SEM-EDX data showed general similarities in metal concentrations. Dissolution testing results indicated similar behaviour between the two samples in a simulated biological fluid.

Introduction

Urban air quality is of significant importance to the majority of us living or working in our cities worldwide. Short-term events such as the 1952 London smog episode have the ability to dominate headlines because of the 4000 deaths due to acute exposure to air pollution, 12,000 when integrated over the whole of the winter (Whittaker et al. 2004). Chronic exposure to air pollution, to which we are all subject, provides equally shocking statistics. The European Environment Agency estimates that poor air quality is responsible for 467,000 premature deaths per year amongst Europeans (EEA 2016). The World Health Organization estimates that air pollution is responsible for 8% of lung cancer deaths, 5% of cardiopulmonary deaths and 3% of respiratory infection deaths (WHO 2009).

Such deleterious effects on the population are often attributed to particulate matter (PM). Emission sources of PM are numerous and commonly include industrial and commercial output, vehicular emissions and wear, abrasion of road surfaces, etc. (Wang et al. 2016). Another large source of PM is believed to be resuspended road dust (RD), which can account for as much as 74% of total suspended particles by mass (Hien et al. 1999; Harrison et al. 1997) and has also been observed to be the largest and second largest contributor to PM10 and PM2.5, respectively (Landis et al. 2017). Emphasis on RD as possibly the most significant contributor to PM will continue to grow because of difficulty regulating this substance and its many contributing materials, unlike vehicle exhaust emissions (Padoan et al. 2017).

Because of the range and nature of sources, PM should be considered unique for any given location. For that reason, PM is a topic which has been studied extensively, although recent consensus amongst the scientific community indicates that the micro-chemical structure of PM is responsible for the degree of toxicity (Rohr et al. 2011).

Due to the nature of PM, and particularly respirable PM10, it can be difficult to analyse micro-chemical structure, without specialised techniques such as micro-Raman spectroscopy (Worobiec et al. 2010) or surface-enhanced Raman spectroscopy (SERS) (Tian et al. 2014). There are further difficulties in obtaining a suitable sample mass with a view to carry out robust bulk in vitro studies, since common methods for modelling either the respiratory or the gastrointestinal route typically use upwards of 0.5 g per sample, including replicates (Colombo et al. 2008; Wragg et al. 2011). Collection on filters using high-volume particulate matter (PM10) samplers typically yields samples of less than 0.1 g (Shahsavani et al. 2012; Behrooz et al. 2017).

The aim of this study is to see if < 38 µm fraction of RD is analogous to the < 10 µm fraction, and thus assess the possibility of using the larger, more abundant size fraction of RD as a proxy for metals in airborne respirable particles originating from RD. This will be carried out using metal concentrations as a marker. A water deposition technique can be used to separate the RD fraction of < 10 µm from the original < 38 µm sample. The < 38 µm sample will also be soaked in deionised water to simulate the conditions to which the < 10 µm fraction was exposed, to allow for both samples to lose some water-soluble components in the process. Both fractions the < 10 µm and the < 38 µm will then be acid digested to compare their elemental profiles and exposed to artificial lung fluid to compare mobility of metals within the body. RD used in the study was collected in 2014 from Oxford Road, Manchester. The collection technique and handling of RD have previously been described by Brown et al. (2015).

Methods

<10 µm particle separation

Size separation was carried out using a sedimentation-in-water process based on the principles of Stokes’ law, which gives the velocity of a sphere in a viscous fluid under the force of gravity. Equation 1 is Stokes’ law, where V is velocity of a particle; g is acceleration due to gravity; a is the diameter of a particle (cm); d1 is the density of a particle; d2 is the density of water; µ is the viscosity of water (Haynes 2018). Mean density of particles was measured to be 3.2 g/cm3, as determined by a water displacement technique similar to Blake and Hartage (1986).

Stokes’ law

The separation technique is presented in Fig. 1. The process uses a sealed 100 cm3 measuring cylinder. The distance between the 100 cm3 mark on the cylinder and the 80 cm3 mark was measured to be 3.1 cm. The time taken for mean density (d1) particles of 10 µm to travel that distance is 259 s, as determined by Eq. 1.

Fig. 1

Separation schematic (https://chemix.org/)

A 2 g sample of RD was accurately weighed into the measuring cylinder and then made up to the 100 cm3 mark with deionised water (18 MΩ). The cylinder and contents were thoroughly agitated by hand and then placed on a bench, and this is shown as step a. in Fig. 1. A stopwatch was then used to time 259 s, the amount of time taken for particles larger than 10 µm to travel 3.1 cm, shown as step b. in Fig. 1. A 20 cm3 pipette was then inserted into the cylinder to siphon off the top 20 cm3 of water. This process removes suspended particles < 10 µm from the cylinder, shown as step c. in Fig. 1. The water in the pipette, containing the < 10 µm particles, is then aspirated into a clean centrifuge vial, shown as step d. in Fig. 1. The centrifuge vial and contents were then centrifuged at 2000 rpm for 15 min, and the supernatant was then discarded, shown as step e. in Fig. 1. The particles remaining in the centrifuge tube were transferred to a clean weighing boat and left to dry at room temperature in a clean brown paper bag. This process is repeated 20 times for each sample with another 20 cm3 of deionised water added to the measuring cylinder before step a. to replace the 20 cm3 removed at step c. on the previous iteration. Twenty extractions were carried out as preliminary testing indicated this was the amount needed for the water in the pipette to appear to be visibly clear of suspended particles. The method described here is based on one described by Boisa et al. (2014).

Particle size analysis

A Zeiss Supra 40VP field emission computer-controlled scanning electron microscope with energy-dispersive X-ray microanalysis (CC-SEM–EDX) was used to determine particle size distribution (PSD) of both samples once mounted by dispersion onto silver foil. A backscattered electron detector was used at 1000 × magnification to acquire an image of the stub surface; the heavy elemental background from the silver foil enables a high contrast black and white image to be obtained. The smallest measurable diameter was 0.5 µm as particles smaller than this are indistinguishable from imperfections in the silver foil surface. EDAX Genesis software is used to control the instrument, generating a series of images spiralling out from the centre of the stub, essentially mapping it. The program recognises particles based on their dark colour relative to the background, and greyscale sensitivity is programmed by the operator beforehand. The program measures the x and y ferets of a particle based on a number of pixels.

Particle size can be represented by average diameter, calculated from the two-dimensional projection of a given particle. Equation 2, as used by Potgieter-Vermaak et al. (2012), shows how average particle diameter is calculated for each particle where Dp is average particle diameter, Dmax is maximum feret diameter, Dmin is minimum feret diameter.

Average particle diameter

CC-SEM–EDX clustering

Clustering was carried out using the same instrumentation technique as described in “Particle size analysis” section. Each particle identified is exposed to an acquisition of 15 s at 25 kV giving metal concentrations as percentage of each particle based on signal peaks at discrete KeV values indicative of either Kα emission or Lα emission. Elemental data were then separated into clusters based on mole fractions per element, as described by Kandler et al. (2011). Clusters for these data are defined in Table 1, along with mole fraction criteria. AE refers to total mole fraction; in each case, the contributions of C, O and Ag have been subtracted.

Table 1

SPA clusters and criteria

Cluster	Criteria
Silicates	Si/AE > 0.2, Mg/Si < 1.33, Al/Si < 1.33, Fe/Si < 0.5, Ti/Si < 0.5, Na/Si < 0.7
Calcites	Ca/AE > 0.5, Mg/Ca < 0.33, Si/Ca < 0.5, S/Ca < 0.25, P/Ca < 0.15
Ti-rich	Ti/AE > 0.3, Mg/Ti < 1, Al/Ti < 1, Fe/Ti < 1, Si/Ti < 1, S/Ti < 1, Na/Ti < 1
Fe-rich	Fe/AE > 0.15, Si/Fe < 1.05, Ti/Fe < 1.3
Al-rich—not aluminosilicates	Al/AE > 0.15, Si/AE < 0.2, Ti/Al < 1.3
Trace containing	(Cu + Cr + Mn + Ni + Pb + Sn + V+Zn)/AE > 0.05

For the CC-SEM–EDX clustering, another aliquot of the < 38 µm fraction was used which was not exposed to water in the manner described in “Particle size analysis”. This was done in order to minimise any loss of water-soluble metals. Particles have been clustered in accordance with Table 1, with extra criteria added to group particles with a diameter of < 10 µm.

Acid-soluble metals determination

250 mg ± 10 mg of each < 10 µm separated sample was analysed for its acid-soluble metal content by dissolution in aqua regia (7.5 cm3 HCl, 2. cm3 HNO3); three extractions were carried out per sample. This process was facilitated with a 1600 W CEM Mars 5 Microwave equipped with Teflon digestion tubes. Analysis was carried out using a Varian Vista MPX ICP-OES with SeaSpray nebuliser. Initial system stability checks were carried out with a 5 mg cm−3 solution of Mn during the torch alignment whereby the upper values must be in excess of 300,000 counts per second, and sample analysis exclusively used the axial position. A four-point linear calibration was achieved for each element; a calibration coefficient of > 0.99 for all elements was determined. Method verification was carried out using a certified reference material (CRM), river clay sediment LGC6139. 500 mg ± 10 mg of CRM was weighed out in triplicate, digested using aqua regia and analysed using ICP-OES. Determined concentrations of available metals were within 4% of the certified mean concentrations. Five replicates were carried out per sample.

Metal dissolution determination

For the extraction method, 150 mg ± 10 mg of the < 10 µm particles was measured into a PET centrifuge tube along with 15 cm3 of artificial lysosomal fluid (ALF). ALF (Colombo et al. 2008) was selected as an appropriate representation of the lung environment cells encounter following contact with the macrophage immune response (Wiseman and Zereni 2014). The mixture was then placed in an incubator shaker for 24 h. This was deemed an appropriate length of time to expose the samples to ALF, as residence time for particles in the lung is still somewhat disputed by the scientific community. Despite this, the aim of this paper is to investigate the use of larger particles as a proxy for respirable particles so extraction time should not be important. The analysis method for the ALF-extracted samples was the same as for the acid-extracted samples. Al, As, Cr, Cu, Fe, Mn, Ni, Pb, Sn and Zn were quantified in each sample. These particular metals were selected to give an insight into the effective change for metals deemed to be of both crustal and anthropogenic origins (Alomary and Belhadj 2007; Gunawardana et al. 2012).

Results

Particle size was determined using CC-SEM–EDX. A total of 1500 particles were analysed for each of the two fractions. Figure 2 shows particle size distribution as determined by CC-SEM–EDX for the < 10 µm and < 38 µm fractions. The figure divides particles into bins based on particle size, starting with 0.5–1.5 µm and then 1.5–2.5 µm, etc., presented along the x-axis as average particle diameter (Dp). The primary y-axis shows the number of particles in each bin divided by the total number of particles (n/n). Presenting particle size distribution data in this manner is common for particulate matter (Chen et al. 1997; Lin et al. 2005; Potgieter-Vermaak et al. 2012; Yue et al. 2013) and RD (McKenzie et al. 2008); these studies tend to log the number of particles in each bin. The authors felt that logging the number of particles in each bin failed to give a clearer visual representation of the particle size distribution. The secondary y-axis of Fig. 2 shows a cumulative distribution of particles, by summing the number of particles in each preceding bin, as a percentage of the total number of particles (N(Dp)/n). The cumulative distribution is included to give an indication of the type of distribution observed. Both the < 38 µm and the < 10 µm fractions showed somewhat of a log-normal distribution; however, there are too many smaller particles to fit this distribution with statistical confidence. It is observed that the number of particles smaller than 1.5 µm dropped off in the < 38 µm fraction, but carried on increasing in the < 10 µm fraction. It is possible that this was caused by the separation technique failing to remove some particles between 1.5 and 10 µm from the < 38 µm, therefore under representing them in the < 10 µm fraction.

Fig. 2

Particle size distributions from each fraction, where n/n indicates the fraction of particles present in that size range and N(Dp)/n indicates the cumulative distribution

Geometric median and geometric standard deviation are calculated using equations given by Potgieter-Vermaak et al. (2012). The geometric median particle size in the < 10 µm fraction is 1.7 µm (geometric standard deviation 2.3 µm); 94.3% of particles were smaller than 10 µm. In the < 38 µm fraction, the geometric median particle size was 3.5 µm (geometric standard deviation 2.3 µm); 90.4% of particles were smaller than 10 µm. These results show that a good separation efficiency was achieved with this technique and that a significant portion of the particles larger than 10 µm have been removed. The results fail to report particles smaller than 0.5 µm as these particles could not be quantified by the analytical techniques used in this study. The distribution of particles observed here in the < 10 µm fraction was similar to that reported by McKenzie et al. (2008). In both cases, particles were measured down to a size of 0.5 µm, albeit using different instrumentations. Interestingly, results from this study showed a greater exponential skew towards the smallest particles. This could possibly be due to differences in sites or the vacuum collection method used by Mckenzie et al. (2008). The large number of small particles observed here in the < 38 µm fraction supports the concept that RD inherently contains a large number of small particles (Kong et al. 2012). This is a crucial observation with respect to the aim of this study; a significant portion of road dust particles, by number of particles, are capable of entering the thoracic cavity.

Table 2 shows cluster data for the original < 38 µm fraction, and for particles within this sample that have a diameter of < 10 µm, the portion within each cluster observed is presented as a percentage, in order to normalise the number of particles of a given observed size. A large number of particles were found to match the criteria for more than one cluster. These mixed clusters are also defined in Table 2.

Table 2

Abundance of clusters present in the 1500 particles of < 38 µm and < 10 µm fraction

Cluster	Clusters present in particles < 38 µm (%)	Clusters present in particles < 10 µm (%)
Silicates	57.7	56.3
Silicates and calcites mixed	0.3	0.4
Silicates and trace mixed	6.9	7.3
Calcites	0.4	0.4
Calcites and Al-rich mixed	4.2	4.3
Calcites, Al-rich and trace mixed	7.1	7.3
Ti-rich	0.2	0.1
Fe-rich	2.5	2.3
Fe-rich and trace mixed	0.3	0.3
Al-rich	1.6	1.3
Al-rich and trace mixed	8.2	7.4
Trace	7.2	7.3
Undefined	3.4	5.2

Table 2 shows markedly similar clustering data between the two size fractions. The largest absolute difference between the two fractions was the silicates cluster, where the < 38 µm fraction contained 1.4% more silicates than the < 10 µm fraction. Calcites, Ti-rich, Fe-rich, Al-rich and trace clusters were all within 0.3% of each other with respect to the different size fractions. There was a slightly more noticeable difference between the mixed clusters, where the differences between the size fractions ranged up to 0.8%. This could suggest that in both samples particles are commonly present as aggregated collections of particles, in agreement with previous research on the topic (Thorpe et al. 2007).

The undefined row of Table 2 indicates that particles fit none of the clustering criteria. In most circumstances, this classification was attributed to highly carbonaceous particles, containing trace amounts of various other elements. It is observed that there is a significantly larger portion of undefined particles in the < 10 µm size fraction, which could be the result of a larger fraction of soot/vehicle exhaust particles.

Bulk elemental concentration

In order to further verify if the < 38 µm can act as a predictor for bioaccessible behaviour of smaller particles, the bulk elemental concentration of the particles was investigated. Results from the acid digestion experiment indicated a significant difference in metal concentrations observed between the < 38 and < 10 µm fractions. Elemental concentrations were increased in the < 10 µm fraction for each element quantified. Table 3 represents the data from both fractions and includes the % increase from < 38 to < 10 µm.

Table 3

Acid-soluble metal concentrations

	Al	As	Cr	Cu	Fe	Mn	Ni	Pb	Sn	Zn
<10 µm (mg/kg)	212,700 ± 7.9%	49.7 ± 8.6%	674.8 ± 5.1%	7224 ± 0.9%	397,000 ± 1.1%	6180 ± 0.8%	258.1 ± 6.7%	3290 ± 0.8%	535.5 ± 1.4%	11,240 ± 0.9%
<38 µm (mg/kg)	140,100 ± 3.6%	30.0 ± 11%	636.8 ± 3%	4576 ± 1.6%	355,400 ± 1.5%	5788 ± 1.6%	209.7 ± 3.8%	2320 ± 0.4%	461.3 ± 1.6%	7564 ± 1.6%
Increase (%)	51.8	66.3	6.0	57.8	11.7	6.8	23.1	41.8	16.1	48.6

Elemental concentrations in the < 10 µm fraction showed the trend Fe > Al > Zn > Cu > Mn > Pb > Cr > Sn > Ni > As. The same pattern was observed in the < 38 µm fraction, with the exception of Mn, which was more abundant than copper. Albuquerque et al. (2017) concluded a similar pattern of metal concentration of PM10, Zn > Cu > Pb > Mn > Cr > Ni > As, the notable difference being the increased Pb concentration relative to Mn. Other studies investigating metal concentrations of PM10 also indicated a similar order, again with Pb being more abundant than Mn (Lim et al. 2010). It is worth noting that specific metal concentrations in terms of absolute concentration and relative to other metals are very dependent on site. It is also worth acknowledging that the literature quoting concentrations of metals within airborne PM10 will generally be reported in units of mass per volume (of air), as collection on filters is the most common method of quantification.

Such large differences in metal concentrations between the fractions may appear to conclusively disprove any possibility that the two samples may be comparable. However, it is documented that the largest constituent of RD is Si-containing particles, and therefore, it is expected that the concentration of Si increases in larger grain fractions (Brown et al. 2015). It is feasible that the difference between the two samples may be the significantly larger concentration of Si in the < 38 µm fraction. If that was the only reason, one would expect a fairly constant difference across the elements, which is not observed. On the other hand, as we report on acid-soluble elemental profiles and not complete digestion of the sample, we can also attribute some of these differences to selective or incomplete digestion of some elements, thereby explaining the variation in differences observed. The larger < 38 µm fraction will evidently show lower solubility of the silicate-containing particles than the < 10 µm fraction which has a larger surface area-to-volume ratio. It is also well known that the larger fraction of dust is more silicate-rich than the smaller fractions. Furthermore, the differential dissolution between the two fractions could also mean that some of the metals are bonded in a different matter, for example absorbed on to the surface of other particles, resulting in higher yields for the smaller fraction.

An alternative approach to assessing the difference between the aqua regia-soluble metals in the samples, we can refer back to the SEM–EDX work originally carried out as a means of ascertaining particle size distribution. Recalculating the percentage weight (%Wt) composition for each particle once C, O and Si have been removed, we can compare quantities of each remaining element present in each sample. A full suite of crustal and trace elements quantified by SEM–EDX were left in the %Wt equation to increase accuracy. Distribution of each element also quantified by ICP-OES was found to be skewed left (Al, As, Cr, Cu, Fe, Mn, Ni, Pb, Sn and Zn). Therefore, the Mann–Whitney U test was used to compare each of them between the < 10 and < 38 µm faction, the null hypothesis (H0) being that the samples are statistically not the same. With the exception of Fe, Zn and As, it was found that we can accept the alternative hypothesis for all other elements. In the case of As, the SEM–EDX data indicated too few As-containing particles in either fraction to perform a Mann–Whitney U test. With regard to Fe and Zn, the Mann–Whitney U test indicated that these elements were more enriched in the < 10 µm fraction.

ALF-soluble metals’ results

Since the elemental concentrations of the two fractions, as determined by SEM, showed significantly similar profiles, the last step was to test the bioaccessible profiles of the two fractions. Concentrations of ALF leachate were normalised to mg/kg, in the same way as the acid digested samples were. For simplicity, the data can then be represented as bioaccessibility (%) for each element using Eq. 3. Results for each metal are presented in Table 4.

Table 4

Metal leachate profiles from both fractions displayed as percentage bioaccessibility

	Al	As	Cr	Cu	Fe	Mn	Ni	Pb	Sn	Zn
< 10 µm	1.74 ± 0.08	10.79 ± 0.1	5.14 ± 0.08	8.39 ± 0.03	3.11 ± 0.08	7.61 ± 0.06	9.93 ± 0.07	7.56 ± 0.05	5.03 ± 0.1	10.08 ± 0.05
< 38 µm	1.56 ± 0.04	12.35 ± 0.12	3.33 ± 0.04	5.51 ± 0.03	2.24 ± 0.02	5.64 ± 0.02	8.98 ± 0.04	6.26 ± 0.03	4.57 ± 0.02	8.49 ± 0.02
% Underestimation	10.3	− 14.5	35.2	34.3	28.0	25.9	9.6	17.2	9.1	15.8

Bioaccessibility (%)

Standard deviations for each metal are given as calculated by propagation of errors (Miller and Miller 2005). The results do seem to indicate a strong similarity between the leaching behaviours of the two samples across the full range of metals. There is an average difference of just 1.35% between any given elements with respect to fraction size. The largest absolute difference is 2.88% for Cu; the smallest is 0.18% for Al. This appears to support the elemental composition results obtained using SEM–EDX, as well as the clustering data, which suggested similar associations of the elements. It is, therefore, not unreasonable to say that the results suggest that the chemical composition of the two samples is similar in their leaching behaviour in the ALF. Vastly different leaching behaviours between the two samples may suggest that the elements reported on here are differently associated with each other or other elements between the two samples. In general, it is observed that the < 38 µm fraction will underestimate the % bioaccessibility on average by 20%, with notable underestimation of 26–35% for Mn, Fe, Cu and Cr. In comparison with other studies, the bioaccessibility in ALF results achieved here for both size fractions is significantly lower than that collected by Wiseman and Zereni (2014), who achieved bioaccessibility around an order of magnitude higher across the suite of metals analysed, generally between 60 and 80% bioaccessible. It should be noted, however, that the particles analysed in this study were PM10, PM2.5 and PM1, which would be expected to leach metals more readily. Pelfrene et al. (2017) carried out ALF dissolution tests on a range of certified reference materials. Results indicated large variability in dissolution of metals between each of the reference materials; % bioaccessibility of Mn was found to be 5.5% (± 0.1), 44.3% (± 0.2) and 46.8% (± 2.6) amongst the three analysed reference materials; this would indicate that the chemical composition of a specific sample affects the ability of metals to leach into simulated biological fluids. The similarity in bioaccessibilities of the two samples analysed in this study indicates that the metals may be present in the same chemical composition and thus are representative of one another.

Conclusions

The aim of this study was to assess the possibility of using larger size fractions of urban particles to act as a proxy for respirable particles, in terms of metal concentrations and their in vitro profile. The data show that there are a large number of similarities between the leaching behaviours of the two fractions and that similarities also exist in the clustering profile. Despite this, the absolute metal concentrations between the fractions are quite dissimilar with regard to their aqua regia-soluble behaviour. One explanation could be that silica-rich particles are more concentrated in large particles and thus dilute metal concentrations therein. Another explanation could be that metals within the smaller particles are more efficiently solubilised. Similarities in leaching behaviour in artificial lysosomal fluid and elemental association as determined by clustering CC-SEM–EDX data, as presented in this study, suggest that the larger fraction is a precursor to the smaller. We would like to propose that the larger fraction has a significant presence of aggregates which ultimately disintegrate to form the smaller fraction, providing particles of similar chemical composition. This conclusion is touched upon in discussion “CC-SEM–EDX clustering” section.

Referring back to the aim, as outlined in introduction section of this paper, it does appear that the < 38 µm fraction could be used as at least a first approximation of metal enrichment in respirable particles for the purpose of in vitro tests. It therefore enables researchers to easily use a larger sample size and thus provide more representative results. Further work on the topic using intra-particle analytical techniques, to elucidate the chemical composition of metals, for example Raman spectroscopy, could broaden the conclusions drawn here.