New Insights into Health Risk Assessments for Inhalational Exposure to Metal(loid)s: The Application of Aqueous Chemistry Modelling in Understanding Bioaccessibility from Airborne Particulate Matter

Aqueous modelling of chemical speciation in simulated lung fluid (SLF) enables a better understanding of the underlying chemical factors that influence metal(loid) inhalation bioaccessibility from airborne particulate matter. Such an approach can be used to supplement experimental techniques that are integral to the health risk assessment of metal(loid) exposure by inhalational routes. In this paper, we modelled the aqueous chemistry of airborne particulate-bound metal(loid)s (As, Cu, Mn, Pb and Zn) in a SLF based on Gamble’s solution (neutral pH). The modelling was performed using two software packages (Geochemist’s Workbench 14 and OLI Studio 9.5) and a total of five thermochemical databases (GWB Thermo, MINTEQ, PHREEQC, WATEQ4F and the default database for OLI Studio). Modelled results were compared with experimentally determined bioaccessibilities for the NIST 2710a standard reference material (SRM) and with literature-reported bioaccessibilities for NIST 1648a and BCR 038 SRMs. Whilst the models correctly describe the observed increase in bioaccessibility for more dilute solid/liquid extraction ratios, the performance of the models against the fractional bias of the mean (FBmean) and the normalised mean square error (NMSE) statistical metrics was generally outside the acceptance criteria. Findings from an analysis of the main aqueous chemical species predicted to be present in SLF indicate that carbonate and chloride complexes of Cu, Mn, Pb and Zn predominate, whilst free cations (for Cu, Mn and Zn) and hydroxides (for Cu) also play a role in solubilisation. Arsenic is not predicted to form significant complexes with the SLF components and is present in solution mainly as the HAsO42− ion and its conjugate acid, H2AsO4−. For modelled runs where glycine and citrate were present, significant increases in the bioavailability of Cu and Zn were predicted as a result of complexation with these ligands. An additional finding from our experimental bioaccessibility results for NIST 2710a was that the inclusion of the lung fluid surfactant dipalmitoylphosphatidylcholine (DPPC) in the SLF did not significantly affect the bioaccessibility. Our study provides useful insights into the likely aqueous- and solid-phase speciation of metal(loid)s in SLF and highlights that future developments in this area should consider the role of mineralogy and surface interactions.


Introduction
Inhalation of particulate matter (PM), especially the finer fractions, can contribute to adverse health outcomes, including increased incidence of allergies, asthma, bronchitis, emphysema and cancers [1]. Moreover, it is known that the metal and metalloid components of the particulates are also associated with these health effects [1,2]. Upon inhalation, fine particulate fractions, i.e., those with an aerodynamic diameter less than 4 µm, are most A harsher extraction environment for PTEs is provided by artificial lung fluids that simulate phagocytosis, i.e., the engulfment of particles by alveolar macrophages. The macrophage lysosomal fluids contain proteases and bacterial enzymes and are at a much lower pH than the extracellular lung fluids, typically pH 4.5 [14]. These conditions are simulated in artificial lysosomal fluids (ALFs) by the inclusion of glycine, citrate, pyruvate, tartrate, and lactate [7,14,15].
There are contrasting literature opinions as to whether the more conservative ALF should be used to assess bioaccessibility [14,15] or whether uncertainty about the proportion of time that a particle is resident in lysosomal fluid means that SLF is a more realistic extraction fluid [4]. The choice is important because bioaccessibility is generally much higher in ALF compared to SLF and, in some studies, reported to be as high as 100% for Pb, Cu and As [15]. For example, bioaccessibility of Pb in three standards, BCR 723, NIST 2710a and NIST 1648a, ranged from 51 to 76% in ALF compared to 7.1 to 11.9% in modified Gamble's solution [14]. Nevertheless, the difference was less pronounced for Mn, Cu and Zn, while for Cd and Co, the extractabilities were similar for the two fluids and, in some cases, higher for modified Gamble's solution.
There has been much literature debate about the appropriate S/L (g/mL) to use in simulated lung fluid extractions [7], and several authors have calculated the likely range of these values [8,14]. Some studies have opted to carry out extractions over a wide range of values, for example Caboche et al. [8], who used a range of 1/30 to 1/50,000 and Pelfrêne et al., who used a range of 1/1000 to 1/10,000 [14]. This approach is a sensible one because the ratio will be dependent on a wide range of factors, including concentration of airborne particulates and deposition location. However, the chosen ratio has been shown to influence bioaccessibility, for example by Caboche et al., who found significantly lower proportions of bioaccessible Pb and Zn at an S/L ratio of 1/30 compared to 1/500 [8].
The lower bioaccessibilities at more concentrated extraction conditions have been ascribed to the saturation of species in solution and a restricted exposure of particulate surface to the extractant medium [8].
Another significant factor affecting bioaccessibility is the possibility of occlusion of the particulate surface because of the formation of surface complexes, which may limit the extent of dissolution of particle components. Wragg and Klinck [21] and Kastury et al. [20] have both demonstrated that coatings (2-5 µm thick) of Pb phosphates form as a result of Pb complexing with phosphate from the neutral SLFs (Pb-Si-P-Al-O in the case of Wragg and Klinck [21] in Gamble's solution and O-Pb-Ca-P-Si-Al-Fe in the case of Kastury et al. in Hatch's solution [20]). The same coatings are not observed with the ALF, where such complexation is thought unlikely to occur at the lower pH [20].
A more thorough understanding of the aqueous chemistry of PTEs in SLF may help in the prediction of bioaccessibility. Therefore, in this paper, we compare the performance of two aqueous chemistry modelling packages and a total of five thermochemical databases in predicting the equilibrium concentrations of PTE components of the particles when subjected to SLFs, i.e., to compare modelled against observed bioaccessibilities. Data for observed bioaccessibilities are provided by our own experimental determination for NIST 2710a and literature bioaccessibilities values for NIST 1648a and BCR 038 [8]. Moreover, we analyse the aqueous-and solid-phase chemical species that are predicted to be present, and the insights that this provides for PTE bioaccessibility. There are literature examples of aqueous chemistry modelling being applied to a range of related environmental scenarios, such as in the prediction of speciation and solubility of metals and metalloids during leaching from smelter fly ash [19,[22][23][24] and mine tailings [16,17], and in assessing the impact of acid rain on weathering and Al 3+ release in Swedish conifer forests [25]. Specific applications involving bioaccessibility modelling have included the simulation of speciation and solubility of PTEs from mine tailings in simulated lung and gut fluids [16] and the dissolution of uranium from particulates in simulated lung fluids [3]. Most of these studies have used the PHREEQC geochemical modeler [26], with either the MINTEQ [3,17,19,22,23,27] or WATEQ4F [18] thermochemical databases [27,28]. Geochemist's Workbench (GWB) has also been used, with its built-in Thermo thermochemical database [16]. The MINTEQ (3383 aqueous species, 709 minerals, 12 gases [27]), PHREEQC (187 aqueous species, 57 minerals, 8 gases [26]) and WATEQ4F (361 aqueous species, 310 minerals, 8 gases [28]) thermochemical databases are all available through the GWB 14 interface, along with the built-in Thermo database (647 aqueous species, 624 minerals, 10 gases [29]), and are all used in this paper. In addition, we compare predicted bioaccessibilities using another geochemical modelling package, OLI Studio, and its built-in thermochemical framework (80 inorganic elements of the periodic table, and their associated aqueous species, as well as over 8000 organic species) [30]. As far as we are aware, this package has not been applied to the prediction of bioaccessibility, though it has been used in environmental applications [31][32][33].

Experimental Determination of Inhalation Bioaccessibility of As, Cu, Mn, Pb and Zn from NIST 2710a
The observed bioaccessibility of PTEs from NIST 2710a was determined for triplicate samples using the Gamble's solution formulation of Caboche et al.: NaCl, 6400 mg L −1 ; CaCl 2 ·2H 2 O, 255 mg L −1 ; Na 2 HPO 4 , 150 mg L −1 ; NaHCO 3 , 2700 mg L −1 ; Na 3 citrate·2H 2 O, 160 mg L −1 ; NH 4 Cl, 118 mg L −1 ; Glycine, 190 mg L −1 [8]. However, a modification was that the bioaccessibility determinations were made in the presence and absence of DPPC (0.02%) in order to test the effect of lung fluid surfactants on bioaccessibility; these are thought to aid wettability of particles and prevent aggregation [8]. Appropriate sample weights of powdered standard reference material (SRM) to achieve the desired range of S/L (g mL −1 ) ratios (1/30, 1/500, 1/2000 and 1/4000) were mixed with a fixed volume of 30 mL of SLF. The extraction solutions were then placed on an end-over-end shaker at 37 • C, and 40 rotations per minute, for 24 h. The extract was centrifuged at 14,600 rpm for 20 min. The supernatant was then removed, and 100 µL of concentrated nitric acid added and the resultant solution refrigerated at 4 • C until analysed. For analysis, 9 mL of the acidified extract was placed in a tube and 1 mL of internal standard (Sc) added.
The total PTE concentration in NIST 2710a was determined using microwave digestion according to the method of Okorie et al. [34]. A soil sample (0.5 g) was accurately weighed into a PTFE (Teflon) vessel pre-cleaned with 10% HNO 3 solution. Then, 13 mL of aqua regia was added to the sample in a fume hood. The vessels were sealed and placed in the microwave carousel for microwave digestion. The microwave-assisted digestion was achieved using a closed vessel microwave-accelerated reaction system (MARS 5, CEM Corporation). This instrument was set at a power of 440 Watts for 40 min, with a cooling time of 20 min. The maximum pressure and temperature allowed for the XP-1500 vessel ranged from 100 to 600 psi and 240 to 300 • C, respectively. After cooling, the digest was filtered through the filter paper into a 50 mL volumetric flask and made up to the mark with deionised water. Finally, this solution was transferred into a 50 mL screw-cap Sarstedt tube ready for ICP-AES analysis. For analysis of the samples, aliquots of the NIST 2710a extracts, blanks and calibration solutions were prepared by adding a Sc solution (1 ppm) as an internal standard and diluted with 1% HNO 3 solution to ensure that the solutions are of similar matrix with the standard calibration solutions. Approximately 0.1-5 mL of the digested solution was diluted in a 10 mL tube with 1 mL of internal standard added. The volume was made up to 10 mL with 1% HNO 3 solution before the ICP-AES analysis.
The ICP-AES instrument (PerkinElmer Optima-8000, Beaconsfield, UK) was calibrated with a range of concentrations-0.01, 0.1, 1, 2, 5, 10, 20 mg L −1 -and a calibration blank prepared from a multi-element calibration solution. A linear calibration curve was achieved with good correlation coefficient (R 2 > 0.9999) for each of the analytes. Prior to analysis of the sample, a suitable wavelength for each element was selected to avoid chemical interferences. All elemental concentrations were reported as mg kg −1 , dry weight. The ICP-AES operating conditions for the plasma were as follows: RF power, 1300 watts; RF frequency, 40 MHz; nebuliser flow, 0.55 L min −1 ; auxiliary flow, 0.2 L min −1 ; plasma flow, 15 L min −1 ; equilibrium time, 15 s; plasma aerosol type, wet; nebuliser start-up conditions, gradual; plasma view, axial. The peristaltic pump conditions were as follows: sample flow rate, 1.5 L min −1 ; sample flush time, 30 s; wash rate, 2.00 mL min −1 ; wash time, 60 s; read delay, 2 s; replicates, 3; background correction, 2 points (for all elements).

Aqueous Chemistry Simulation
Aqueous chemistry simulations were carried out using two software packages: GWB 14 and OLI Studio 9.5. For both packages, and for all thermochemical databases used, a basic aqueous chemistry simulation was used, whereby it is assumed that all species are initially in the aqueous state, and the calculation is then made of the likely solid states (precipitates) that form, and thus the equilibrium aqueous species that remain. Whilst the models can also simulate surface adsorption and redox reactions, neither of these options were used.

Aqueous Chemistry Simulation Using GWB 14
For the GWB 14 aqueous chemistry modelling, spreadsheets of concentrations for the 'basis' species (the initial aqueous species present prior to any calculation of complexes and precipitates [35]) were prepared for each of the thermochemical datasets used: Thermo, PHREEQC, WATEQ4F and MINTEQ. Some of these thermochemical datasets had a reduced number of available chemical species; for example, WATEQ4F and PHREEQC do not include Cr or Hg species and, PHREEQC does not include As Cr, Hg or Ni species. Furthermore, only MINTEQ contains the glycine and citrate species that are used in Gamble's solution, and none contain the surfactant DPPC. The calculation of the concentration of the initial aqueous species in 1 L of water was made on the basis of (a) the components of the Gamble's extraction solution; (b) the concentrations given in the certificate of analysis for the SRMs (NIST 2710a was used for our extractability work; while NIST 1648a and BCR 038 were used to compare with the extractability data of Caboche et al. [8]); and (c) the S/L ratio being modelled (1/30, 1/500, 1/2000 and 1/4000 to compare with the experimental extraction results obtained in the present paper, with the addition of 1/20,000 and 1/50,000 to compare with the literature extractability data [8] In addition, separate modelling runs were carried out using MINTEQ in the absence and presence of citrate and glycine. The specified physical conditions were pH 7.4 and 37 • C, corresponding to the experimental extraction conditions used to replicate the situation in the human lung.
Once specified, the equilibrium concentrations of the aqueous species were calculated using the GWB React program. Prior to running React, three of the aqueous species needed to be substituted with mineral species in order to constrain the model [35]: quartz replacing H 4 SiO 4 , hematite replacing Fe 3+ and gibbsite replacing Al 3+ . The charge balance species was set to Cl − , and precipitation reactions were enabled. For the main set of modelling runs, the final pH of the solution was calculated by the React program, though a more limited set of simulations was also carried out in which the pH was set at 6.0, 6.8, 7.1 and 7.4. This was performed to allow the investigation of the effect of pH on bioaccessibility at the physiological pHs for particulate-lung fluid interactions, given the possibility of local changes in pH as a result of particle deposition.

Aqueous Chemistry Simulation Using OLI Studio 9.5
For modelling using OLI, a similar procedure was carried out to that specified for GWB 14 in terms of defining the initial aqueous model, though because of the difference in inputting chemical substances in OLI 9.5 compared to GWB 14, hydroxides of the metals were used as inputs. All of the species listed in Section 2.2.1 were included in the model, with the addition of Ni(OH) 2 and Cd(OH) 2 and Sb(OH) 5 , which were not available in the GWB thermochemical databases, but were in OLI Studio. To allow direct comparison to the GWB results, the pH was fixed at the value calculated from the GWB runs for the specified S/L ratio. To achieve the desired pH, HCl was specified as the titrant. The default aqueous chemistry thermochemical database for OLI 9.5 was used, and was run with and without citrate and glycine, to allow comparison with the MINTEQ database used in GWB 14.

Calculation of Modelled Bioccessibilty
Outputs from both GWB 14 and OLI Studio 9.5 give aqueous and whole-system elemental totals, allowing a straightforward calculation of the % of the aqueous species present in solution compared to the theoretical maximum. For the extractability work carried out in the present paper, comparisons were made with As, Cu, Mn, Pb, and Zn. For comparisons with the work of Caboche et al. [8], only Zn and Pb were used, as these were the only two PTEs for which extractability data were determined over a range of S/L ratios. In addition to the elemental totals in the solid and liquid phases, GWB 14 and OLI 9.5 generate a full aqueous-and solid-phase speciation of the elements.

Evaluation of Model Performance
To obtain a more objective consideration of the performance of the models, two statistical tests were used: fractional bias of the mean (FB mean ), Equation (1), which gives an indication of the extent of systematic bias, and normalised mean square error (NMSE), Equation (2), which quantifies the random bias. Both of these tests are used to evaluate environmental models in the literature [36,37]. In Equations (1) and (2), C obs and C pred are the observed and predicted concentrations, respectively. For both FB mean and NMSE, a perfect model would give values of zero. However, in evaluative terms, FB mean (absolute) values of <0.3 (i.e., less than 30% deviation from the actual mean) and NMSE values <1.5 are generally considered acceptable [37].  Table 2 details the measured bioaccessibility of As, Cu, Mn, Pb and Zn in the SLF, both in the presence and absence of DPPC surfactant. In both cases, the lowest bioaccessibilities are associated with the highest S/L ratios, i.e., the most concentrated solutions, as has previously been observed [8,14,17]; this is likely to be due to species becoming saturated at these higher concentrations and also the possibility of a reduced surface area for the exposed particulate [8]. The difference in % bioaccessibility at an S/L ratio of 1/30 compared to the more dilute solutions is more significant for our study compared to other studies, for example that by Caboche et al. [8], though this may simply reflect the different SRMs used. Where % bioaccessibility data are present for NIST 2710a from other studies [14], the S/L ratios used are not equivalent, and so a direct comparison cannot be made. Nevertheless, Pelfrêne et al. did observe a significant S/L ratio-dependent extraction profile for many metals between the S/L ratios 1/1000 and 1/10,000, most notably for Pb (2.2-11.8%), Zn (13.2-25.9%) and Mn (31.8-40.9%), though Cu ranged in the opposite direction (49.4-36.9%) [14].
We found that the general order of extractability for NIST 2710a (from highest to lowest) is Cu > Mn~As > Zn > Pb. Table 3 shows that our results are in very good agreement with previous inhalation bioaccessibility studies on this SRM [14,15], particularly for the study of Pelfrêne et al., where two variants of Gamble's solution were used, with 'Modified Gambles' referring to the addition of DPPC [14]. The study of Meza-Figueroa et al. [15] reported lower extractabilities for all elements compared to both our study and that of Pelfrêne et al. [14]. This might be explained by the different SLF compositions used, particularly the absence of glycine in the Meza-Figueroa et al. study [15]. It has previously been reported that glycine can form complexes with several different metals ions, thus aiding extractability [15].   [14], who compared SLF extractabilities in Gambles solution, with and without DPPC present. For the five PTEs investigated in our study, Pelfrêne et al. [14] previously found no statistical difference in % bioaccessibility for Pb or Mn in any of the three SRMs. Small, but statistically significant differences were observed with Cu (lower when DPPC present) and Zn (higher when DPPC present) for NIST 2710a. Nevertheless, Li et al. [38] found that the Pb bioaccessibility range for PM 2.5 sampled in Nanjing reduced from 23-43% to 5.6-18% if DPCC was omitted from Gamble's solution. The addition of DPCC is thought to aid wettability of the particles and prevent aggregation [5,8,39].

Comparison of Modelled and Experimentally Determined Bioaccessibilities
For our modelling predictions using the five thermochemical databases, we have made comparisons with both our own experimentally determined bioaccessibilities for NIST 2710a (Montana soil), Figure 1, and with bioaccessibilities reported in the literature for NIST 1648a (urban atmospheric PM), and BCR 038 (fly ash), Figure 2. For two of the databases, MINTEQ and OLI, the inclusion of glycine and citrate as species allowed a comparison to be made of bioaccessibility in the presence and absence of these metal ion chelators, that are integral components of the Gamble's extraction solution [5,8]. In addition to the extractability profiles presented in Figures 1 and 2, we have summarised in Table 4 the most significant aqueous species that are predicted to be present for NIST 2710a at an S/L ratio of 1/500 g mL −1 for each of the thermochemical databases. The top five most abundant species (if present) are presented in order, though the first two or three species will generally account for~95% or more of the aqueous phase for each of the species. A detailed analysis of the concentrations of the species has not been presented, as this will vary according to pH, the S/L ratio and elemental concentration for the SRMs.
Nevertheless, the list of species gives a useful comparison of the components of the SLF that solubilise the PTEs under investigation. Table 4 also shows the likely precipitating species.  Figure 3. Influence of pH on observed against modelled PTE% bioaccessibilities for NIST 2710a. Modelling was carried out only for MINTEQ (+gly. and citr.) and OLI (+gly. and citr.), at the pHs specified in the legend.

Conclusions
In this paper, we have evaluated whether aqueous chemical modelling can allow a better understanding of the underlying chemical factors influencing PTE inhalation bioaccessibility from airborne particulate matter. Such an approach can be used to supplement experimental techniques that are integral to the health risk assessment of PTE exposure by inhalational routes. This is important because experimental bioaccessibilities for As, Cu, Pb, Mn and Zn in neutral SLF vary considerably, and we need to understand the reasons for this variability. We have shown that aqueous chemical modelling may provide useful insights into the components of SLF that aid the solubilisation of the PTEs.
One of the most important insights arising from the work is that, when modelled with five different aqueous chemical thermochemical databases, there was a consensus amongst the models that at the lowest S/L ratios (greatest dilutions), As, Cu, Pb, Mn and Zn could be up to 100% bioaccessible. Mineralogy, surface occlusion, surface complexation and adsorption are the likely explanations for the variability in observed bioaccessibility, and for the discrepancy between modelled and experimental results. The models used do allow for surface adsorption and specific minerals to be included; however, in most cases, the mineralogy of modelled particles will be unknown and so it will not be possible to include these aspects in the model. Nevertheless, a conservative approach to risk assessment might be to assume 100% bioaccessibility on the basis of the uncertainty of how the S/L ratio and residence time will vary during the lifetime of the particle in the respiratory system, and the difficulty in determining a detailed mineralogy for the PTEs of concern.
From the evaluation of the modelled results using the statistical performance measures, FBmean and NMSE, none of the models stood out in terms of meeting the performance criteria across the range of PTEs, though some did meet both of the performance criteria for individual PTEs. These were the GWB Thermo database for Cu and WATQ4F for Zn for our NIST 2710a experimental results, and Zn modelled using MINTEQ for the Caboche et al. [8] comparison. OLI also appears to give a reasonable performance for Pb Figure 3. Influence of pH on observed against modelled PTE% bioaccessibilities for NIST 2710a. Modelling was carried out only for MINTEQ (+gly. and citr.) and OLI (+gly. and citr.), at the pHs specified in the legend.

Conclusions
In this paper, we have evaluated whether aqueous chemical modelling can allow a better understanding of the underlying chemical factors influencing PTE inhalation bioaccessibility from airborne particulate matter. Such an approach can be used to supplement experimental techniques that are integral to the health risk assessment of PTE exposure by inhalational routes. This is important because experimental bioaccessibilities for As, Cu, Pb, Mn and Zn in neutral SLF vary considerably, and we need to understand the reasons for this variability. We have shown that aqueous chemical modelling may provide useful insights into the components of SLF that aid the solubilisation of the PTEs.
One of the most important insights arising from the work is that, when modelled with five different aqueous chemical thermochemical databases, there was a consensus amongst the models that at the lowest S/L ratios (greatest dilutions), As, Cu, Pb, Mn and Zn could be up to 100% bioaccessible. Mineralogy, surface occlusion, surface complexation and adsorption are the likely explanations for the variability in observed bioaccessibility, and for the discrepancy between modelled and experimental results. The models used do allow for surface adsorption and specific minerals to be included; however, in most cases, the mineralogy of modelled particles will be unknown and so it will not be possible to include these aspects in the model. Nevertheless, a conservative approach to risk assessment might be to assume 100% bioaccessibility on the basis of the uncertainty of how the S/L ratio and residence time will vary during the lifetime of the particle in the respiratory system, and the difficulty in determining a detailed mineralogy for the PTEs of concern.
From the evaluation of the modelled results using the statistical performance measures, FB mean and NMSE, none of the models stood out in terms of meeting the performance criteria across the range of PTEs, though some did meet both of the performance criteria for individual PTEs. These were the GWB Thermo database for Cu and WATQ4F for Zn for our NIST 2710a experimental results, and Zn modelled using MINTEQ for the Caboche et al. [8] comparison. OLI also appears to give a reasonable performance for Pb and Zn. Modelled results are for the thermochemical databases described in the legend. No experimental bioaccessibilities were determined for S/L ratios of 1/20,000 or 1/50,000, though modelling was carried out at these S/L ratios so as to be consistent with those used by Caboche et al. [8]. Note that for As, all the curves are overlaid at 100%.

Conclusions
In this paper, we have evaluated whether aqueous chemical modelling can allow a better understanding of the underlying chemical factors influencing PTE inhalation bio- Figure 3. Influence of pH on observed against modelled PTE% bioaccessibilities for NIST 2710a. Modelling was carried out only for MINTEQ (+gly. and citr.) and OLI (+gly. and citr.), at the pHs specified in the legend.

Conclusions
In this paper, we have evaluated whether aqueous chemical modelling can allow a better understanding of the underlying chemical factors influencing PTE inhalation bioac-

Comparison of Modelled and Predicted Bioaccessibilities for NIST 2710a
From Figure 1, we can see that the bioaccessibilities predicted by the different thermochemical databases vary considerably and have a mixed performance with respect to correlation with experimentally determined values. Nevertheless, the experimentally determined trend of increasing bioaccessibility with lower S/L ratios is replicated by the modelling, though different models seem to work better for specific PTEs. Thus, OLI Studio seems to describe Pb bioaccessibility well, whereas WATQ4F would appear to work better for Zn, and PHREEQC for Mn. Of the literature explanations for the relationship between the S/L ratio and bioaccessibility, i.e., saturation of chemical species at the highest S/L ratios and reduced surface area for extraction [8], the modelling can only account for the former, though this seems to be a dominant process.
FB mean and NMSE were applied to predicted vs. actual bioaccessibility data for the individual metals. The results of the analysis are shown in Table 5. The values obtained for both FB mean and NMSE show that the models generally perform outside the acceptability criteria, though there are exceptions, such as the performance of the GWB Thermo database for Cu and WATQ4F for Zn, where both FB mean and NMSE have values within the performance criteria. OLI also appears to give a reasonable performance for Pb and Zn (NMSE conforms to the performance criteria, both with and without glycine/citrate, though the latter performs better). Perhaps the most important observation from the results of modelling studies shown in Figure 1 is that many of the models predict very high bioaccessibilities at the most dilute S/L ratios for all PTEs, with several predicting high bioaccessibilities at medium S/L ratios. Thus, for As, all models predict 100% bioaccessibility across all S/L ratios. For Cu and Zn all models predict 100% bioaccessibility at the lowest S/L ratio. For Mn, all models except WATEQ4F; and for Pb, all models except WATEQ4F and MINTEQ predict high bioaccessibilities at low S/L ratios.
From Table 4, we can see that for modelling runs conducted in the absence of glycine and citrate, the main solubilising SLF components are carbonate and chloride across all species apart from As. Free ions are also predicted to play a role in the solubility of Cu, Mn and Zn, and hydroxide complexes are, additionally, predicted to aid the solubility of Cu. In simulated leaching experiments of fly ash components into water, Ettler et al. found that chloro complexes tend to dominate the solubility of all metals studied (Pb, Cu, Cd and Zn), with free ions also significant for Cu and Zn [19,22]. Additional contributions to aqueous solubility were provided by a number of minor soluble sulphate and carbonate complexes, though the pHs were lower than in our study [19,22]. From the species data in Table 4, the metals will precipitate out of solution as carbonates, and phosphates. Arsenic solubility appears unaffected by the SLF components, but are instead dependent on the presence of the HAsO 4 2− ion and its conjugate acid, H 2 AsO 4 − (the pKa of the latter, arsenic acid, is 6.97 [40], so a similar pH to the SLF).
The major difference between the modelling and experimental extraction conditions is the presence of the metal chelating species glycine and citrate in the SLF: all experimental extractions were carried out in the presence of these species, whereas for the modelling, only the MINTEQ and OLI databases allowed the effect of these species to be modelled. We can see from Figure 1, for example, that the presence of glycine/citrate shows a significant increase in the bioavailability of Cu. Both MINTEQ and OLI have a large, predicted increase in bioaccessibility at high S/L ratios for Cu and, to a smaller extent, Zn. From Table 4, it is clear that for both OLI and MINTEQ, at an S/L ratio of 1/500 (pH~7.0), the diglycinate complex of Cu is the most significant aqueous Cu species; and for Zn, the Zn-citrate complexes make a significant contribution to the aqueous stability (though the Zn 2+ ion and carbonate complexes are also significant components). Manganese is also predicted to be solubilised by citrate, though the effect is marginal at an S/L ratio of 1/500. The results of the statistical evaluation in Table 5 show a better performance for the prediction of copper concentrations when citrate and glycine are included in OLI and MINTEQ, though this improvement is not evident for any of the other PTEs studied.
The discrepancy between the predicted and experimental bioaccessibilities for many of the PTEs is likely to be explained by the underlying mineralogy of the PTEs components in NIST 2710a, as well as the possibility of surface complexation and surface occlusion. Many authors have discussed the role of mineralogy as a determining factor in bioaccessibility and have identified the minerals that are likely to be associated with the highest bioaccessibilities, such as chlorides and sulphates [3,8,[16][17][18][19][20][21][22][23]41]. However, for most particulate samples, a detailed mineralogy will not be available and so cannot be reflected in the modelling, unlike for specific sources such as in the modelling work carried out by Thomas et al. [16] and Helser et al. [17]. Moreover, notwithstanding the mineralogy, we do know from experimental studies that much higher bioaccessibilities are possible when using the lower pH ALF [14,15], even up to 100% for Cu, Pb and As [15], and so it is likely that other factors are also involved. There is good evidence of surface occlusion being a factor in SLFs, particularly for Pb containing particles (galena particles from mine waste), where, over the course of 24 h, a 2-5 µm thick surface coating of lead phosphate develops, through intermediate Si-P-Al-O phases. A similar process has been observed by Kastury et al. in neutral Hatch's solution, where a 3-4 µm thick O-Pb-Ca-P-Si-Al-Fe forms on PM 10 particles taken from areas of mining/smelting [20]. Such surface occlusion will inhibit further interaction between the particles and the SLF. Table 4 shows that metal phosphates are likely to precipitate from solution, including for Pb (MINTEQ), supporting the plausibility that such reactions might also occur at the particle surface. In addition, there is the potential of surface adsorption of released PTE species, including specific surface complexation processes such as between iron minerals (oxides and oxohydroxysulphates) and AsO 4 3− and Pb 2+ [16].

Comparison of Modelled and Predicted
Bioaccessibilities for NIST 1648a and BCR 038 Figure 2 shows similar trends to Figure 1 in terms of the modelling results, though in this case, for comparison to the literature study of Caboche et al. [8], only data for Pb and Zn are available. The comparison of modelled concentrations for NIST 1648a and BCR 038 allows us to observe the effect of particulate PTE concentration on bioaccessibility because of the approximately ten-fold difference in concentrations of Pb and Zn between the two standards (with NIST 1648 the higher of the two and having similar concentrations to NIST 2710a, Figure 1). In addition, taking Figures 1 and 2 together, we can see the wide range of reported experimental bioaccessibilities, particularly for Zn.
As with the data presented in Figure 1, most of the models predict very high bioaccessibilities for Pb and Zn at the lowest S/L ratio, with the exception of WATEQ4F and MINTEQ (with and without glycine/citrate) for Pb and WATEQ4F for the NIST 1648a Zn result. The default Thermo database in GWB 14 predicts 100% extractability at almost all S/L ratios and for all PTEs; however, we see that for Zn in Figure 2, the curve for the Thermo database is displaced downwards, especially for BCR 038 (to 79%), where the concentration of Zn is approximately one-tenth of that for NIST 2710a and NIST 1648a. The detailed speciation output from the runs indicates that this is due to the formation of a zinc chromate precipitate. Another anomaly is the drop in Zn concentration observed at an S/L ratio of 1/500 for NIST 1648a using the PHREEQC and OLI (without glycine/citrate) models: this is due to the predicted formation of Smithsonite (ZnCO 3 ). Table 6 reports the results of the FB mean and NMSE analysis for the Caboche et al. comparison alone, and for a combined analysis of the Pb and Zn results for BCR 038, NIST 1648a and NIST 2710a. As with the results for NIST 2710a alone ( Figure 1 and Table 5), most of the parameters exceed the accepted performance criteria, except for Zn modelled using MINTEQ (FB mean and NMSE both withing the performance criteria). Table 6. Performance metrics for bioaccessibility modelling of NIST 1648a and BCR 038 [8]. Values in bold indicate 'acceptable' performance [37]. The reasons for the discrepancies between predicted and experimentally determined bioaccessibilities in SLF will be the same as discussed for our work on NIST 2710a, though in the case of the comparisons with the results of Caboche et al. [8], we can see that such discrepancies extend to S/L ratios of 1/20,000 and 1/50,000 (our own work has a lowest S/L ratio of 1/4000).

Effect of pH
As the deposition of particulate matter in the lungs is likely to alter the local pH of the lung fluid, it was important to carry out a series of runs in which the effect of pH was investigated for the range of values observed in our modelling studies. Figure 3 shows bioaccessibility profiles for PTEs in NIST 2710a over the pH range of 6-7.4, modelled using MINTEQ (with glycine and citrate) and OLI. It is evident that whilst there is relatively little effect of pH on the extractabilities of As (both oxidation states were modelled, though only As(V) is shown), Cu and Zn, there is a more significant effect for Mn and Pb. For Pb, modelled using OLI, bioaccessibilities greater than 60% are predicted at S/L ratios of 1/50,000; however, lower pHs result in reduced bioaccessibilities (excluding S/L ratios of 1/30 and 1/500). For Pb modelled using MINTEQ, there is a less clear relationship, with the profiles appearing more similar than for OLI. From the speciation output from the OLI runs, the principle difference accounting for the increased aqueous concentration of Pb is the increased abundance of soluble Pb(CO 3 ) 2 2− at higher pHs. For Mn, modelling results for both OLI and MINTEQ show a higher bioaccessibility at lower pHs. The output from the OLI runs indicates that at the lower pHs, there is an increased abundance of the three principal aqueous species, Mn-citrate − , Mn 2+ and MnCl + , most probably due to pH-dependent competitive interactions with other PTEs and matrix components. We know from the literature that pH has the potential to play a major role in bioaccessibility, with experimental results showing high bioaccessibilities for Pb, As and Zn when using ALF (pH 4.5) [14,15]. In the case of ALF, this may be due to the increased formation of soluble metal chloride complexes at the lower pH, a phenomenon that also seems to partially account for the changes in Mn solubility.  Figure 3. Influence of pH on observed against modelled PTE% bioaccessibilities for NIST 2710a. Modelling was carried out only for MINTEQ (+gly. and citr.) and OLI (+gly. and citr.), at the pHs specified in the legend.

Conclusions
In this paper, we have evaluated whether aqueous chemical modelling can allow better understanding of the underlying chemical factors influencing PTE inhalation b accessibility from airborne particulate matter. Such an approach can be used to supp ment experimental techniques that are integral to the health risk assessment of PTE exp sure by inhalational routes. This is important because experimental bioaccessibilities As, Cu, Pb, Mn and Zn in neutral SLF vary considerably, and we need to understand t reasons for this variability. We have shown that aqueous chemical modelling may provi useful insights into the components of SLF that aid the solubilisation of the PTEs.
One of the most important insights arising from the work is that, when modell with five different aqueous chemical thermochemical databases, there was a consens amongst the models that at the lowest S/L ratios (greatest dilutions), As, Cu, Pb, Mn a Zn could be up to 100% bioaccessible. Mineralogy, surface occlusion, surface comple tion and adsorption are the likely explanations for the variability in observed bioacces bility, and for the discrepancy between modelled and experimental results. The mod used do allow for surface adsorption and specific minerals to be included; however, most cases, the mineralogy of modelled particles will be unknown and so it will not possible to include these aspects in the model. Nevertheless, a conservative approach risk assessment might be to assume 100% bioaccessibility on the basis of the uncertain of how the S/L ratio and residence time will vary during the lifetime of the particle in t respiratory system, and the difficulty in determining a detailed mineralogy for the PT of concern.
From the evaluation of the modelled results using the statistical performan measures, FBmean and NMSE, none of the models stood out in terms of meeting the perf mance criteria across the range of PTEs, though some did meet both of the performan Figure 3. Influence of pH on observed against modelled PTE% bioaccessibilities for NIST 2710a. Modelling was carried out only for MINTEQ (+gly. and citr.) and OLI (+gly. and citr.), at the pHs specified in the legend.

Conclusions
In this paper, we have evaluated whether aqueous chemical modelling can allow a better understanding of the underlying chemical factors influencing PTE inhalation bioaccessibility from airborne particulate matter. Such an approach can be used to supplement experimental techniques that are integral to the health risk assessment of PTE exposure by inhalational routes. This is important because experimental bioaccessibilities for As, Cu, Pb, Mn and Zn in neutral SLF vary considerably, and we need to understand the reasons for this variability. We have shown that aqueous chemical modelling may provide useful insights into the components of SLF that aid the solubilisation of the PTEs.
One of the most important insights arising from the work is that, when modelled with five different aqueous chemical thermochemical databases, there was a consensus amongst the models that at the lowest S/L ratios (greatest dilutions), As, Cu, Pb, Mn and Zn could be up to 100% bioaccessible. Mineralogy, surface occlusion, surface complexation and adsorption are the likely explanations for the variability in observed bioaccessibility, and for the discrepancy between modelled and experimental results. The models used do allow for surface adsorption and specific minerals to be included; however, in most cases, the mineralogy of modelled particles will be unknown and so it will not be possible to include these aspects in the model. Nevertheless, a conservative approach to risk assessment might be to assume 100% bioaccessibility on the basis of the uncertainty of how the S/L ratio and residence time will vary during the lifetime of the particle in the respiratory system, and the difficulty in determining a detailed mineralogy for the PTEs of concern.
From the evaluation of the modelled results using the statistical performance measures, FB mean and NMSE, none of the models stood out in terms of meeting the performance criteria across the range of PTEs, though some did meet both of the performance criteria for individual PTEs. These were the GWB Thermo database for Cu and WATQ4F for Zn for our NIST 2710a experimental results, and Zn modelled using MINTEQ for the Caboche et al. [8] comparison. OLI also appears to give a reasonable performance for Pb and Zn. Influence of pH on observed against modelled PTE% bioaccessibilities for NIST 2710a. Modelling was carried out only for MINTEQ (+gly. and citr.) and OLI (+gly. and citr.), at the pHs specified in the legend.

Conclusions
In this paper, we have evaluated whether aqueous chemical modelling can allow a better understanding of the underlying chemical factors influencing PTE inhalation bioaccessibility from airborne particulate matter. Such an approach can be used to supplement experimental techniques that are integral to the health risk assessment of PTE exposure by inhalational routes. This is important because experimental bioaccessibilities for As, Cu, Pb, Mn and Zn in neutral SLF vary considerably, and we need to understand the reasons for this variability. We have shown that aqueous chemical modelling may provide useful insights into the components of SLF that aid the solubilisation of the PTEs.
One of the most important insights arising from the work is that, when modelled with five different aqueous chemical thermochemical databases, there was a consensus amongst the models that at the lowest S/L ratios (greatest dilutions), As, Cu, Pb, Mn and Zn could be up to 100% bioaccessible. Mineralogy, surface occlusion, surface complexation and adsorption are the likely explanations for the variability in observed bioaccessibility, and for the discrepancy between modelled and experimental results. The models used do allow for surface adsorption and specific minerals to be included; however, in most cases, the mineralogy of modelled particles will be unknown and so it will not be possible to include these aspects in the model. Nevertheless, a conservative approach to risk assessment might be to assume 100% bioaccessibility on the basis of the uncertainty of how the S/L ratio and residence time will vary during the lifetime of the particle in the respiratory system, and the difficulty in determining a detailed mineralogy for the PTEs of concern.
From the evaluation of the modelled results using the statistical performance measures, FB mean and NMSE, none of the models stood out in terms of meeting the performance criteria across the range of PTEs, though some did meet both of the performance criteria for individual PTEs. These were the GWB Thermo database for Cu and WATQ4F for Zn for our NIST 2710a experimental results, and Zn modelled using MINTEQ for the Caboche et al. [8] comparison. OLI also appears to give a reasonable performance for Pb and Zn.
The models did predict the experimentally observed increase in bioaccessibility for more dilute S/L ratios, and the importance of glycine and citrate in solubilising metals (OLI and MINTEQ). The modelling confirms that pH is an important determinant of bioaccessibility, particularly for Pb and Mn over a narrow range of possible pH values that might occur as a result of particulate matter interaction with lung fluid, especially at the highest S/L ratios. Our experimental bioaccessibility determinations using NIST 2710a showed that DPPC was not a significant factor affecting bioaccessibility.
Overall, this study provides insights into the likely aqueous-and solid-phase PTE speciation that occurs in a basic model that excludes the influence of mineralogy and surface interactions. We propose to extend the work to include aspects of mineralogy and surface adsorption. As far as we are aware, this is the first comparative study to have been carried out on such a wide range of aqueous chemistry thermochemical models.