Comparison of Atmospheric Travel Distances of Several Pahs Calculated by Two Fate and Transport Models (the Tool and Elpos) with Experimental Values Derived from a Peat Bog Transect

Multimedia fate and transport models are used to evaluate the long range transport potential (LRTP) of organic pollutants, often by calculating their characteristic travel distance (CTD). We calculated the CTD of several polycyclic aromatic hydrocarbons (PAHs) and metals using two models: the OECD P OV & LRTP Screening Tool (The Tool), and ELPOS. The absolute CTDs of PAHs estimated with the two models agree reasonably well for predominantly particle-bound congeners, while discrepancies are observed for more volatile congeners. We test the performance of the models by comparing the relative ranking of CTDs with the one of experimentally determined travel distances (ETDs). ETDs were estimated from historical deposition rates of pollutants to peat bogs in Eastern Canada. CTDs and ETDs of PAHs indicate a low LRTP. To eliminate the high influence on specific model assumptions and to emphasize the difference between the travel distances of single PAHs, ETDs and CTDs were analyzed relative to the travel distances of particle-bound compounds. The ETDs determined for PAHs, Cu, and Zn ranged from 173 to 321 km with relative uncertainties between 26% and 46%. The ETDs of two metals were shorter than 325 those of the PAHs. For particle-bound PAHs the relative ETDs and CTDs were similar, while they differed for Chrysene.


Introduction
In the atmosphere, some semi-volatile organic chemicals (SVOCs) are transported long distances before being deposited [1].To evaluate whether the environmental release of a man-made chemical may need to be curtailed through international agreements like the Stockholm Convention on Persistent Organic Pollutants, a number of compound-specific hazard indicators is often determined and compared with threshold criteria [2].One of these indicators is a contaminant's spatial range or long-range transport potential (LRTP) [3], which can be characterized by the characteristic travel distance (CTD) [4].The CTD is defined as the distance from a source where the concentration of a chemical is reduced by 63%.A fraction of the chemical can travel a multiple of the CTD.
Several models such as the "Environmental Long-range Transport and Persistence of Organic Substances Model" (ELPOS) [5] and the P OV & LRTP Screening Tool ("The Tool") [6] have been developed and applied to estimate the CTD of various chemicals.Both are steady-state fugacity models that estimate the CTD based on a compound's equilibrium partitioning properties, transformation rates in environmental media, and environmental properties, such as wind speed.The Tool has been designed to be as basic as possible while including all essential information [7], whereas ELPOS is more complex and includes more compartments and additional options.Both models can be used to assess a chemical's LRTP by comparing and ranking its CTD against that of identified persistent organic pollutants [8].By comparing the results of various models the influence of different model assumptions and design features can become apparent [9].Both models are not intended to simulate real world conditions.However, to increase confidence in the results and to identify possible model failure, their estimates should qualitatively agree with empirically derived data.To date, experimental data that facilitate the evaluation of these models with respect to the estimated CTDs are sparse [10,11].
Polycyclic aromatic hydrocarbons (PAHs) are emitted by natural combustion processes, such as forest fires, but anthropogenic emissions dominate, specifically through the incomplete combustion of carbon-based fuels.They undergo atmospheric long-range transport [12].Their saturation vapor pressures and, therefore, their gas/particle partitioning characteristics span a wide range [13].Based on their gas/particle partitioning behavior in the rural atmosphere of Eastern Canada, PAHs can be divided in three groups [14]: g-PAHs that are almost exclusively present in the gas phase (e.g., phenanthrene (Phen) and fluoranthene (Flt)), c-PAHs (c for changing partitioning behavior) that are present in the gas phase in summer and particle-bound in winter (e.g., benzo[a]anthracene (B[a]A) and chrysene (Chry), and p-PAHs that are almost exclusively particle-bound all year long (e.g., benzo[b+k]fluoranthene (B[b+k]F), benzo[a]pyrene (B[a]P), and indeno [1,2,3-cd]pyrene(Ind)).To empirically determine LRTP, measurements at variable distances from a distinct contaminant source generating measurable concentrations have to be done [5].The rates of decline of log-transformed concentrations over large latitudinal gradients in North America have been used to determine empirical travel distance (ETDs) of pesticide chemicals [10,11].A general problem with such an approach is that additional sources along a long transect may affect the results.Moreover, substantial changes in temperature and precipitation along latitudinal gradients might complicate data interpretation due to the temperature dependence of gas/particle partitioning and atmospheric reaction rates and the dependence of deposition rates on the magnitude and temporal variability of precipitation.This study aimed at using a shorter, regional transect downwind of a dominant point source to determine ETDs of a range of PAHs, which were then used to evaluate the CTDs derived with The Tool and ELPOS.
CTDs not only scale linearly with wind speed, but also rely on several simplifying and generalizing assumptions and are based on a simplistic representation of the environment.Therefore, absolute values given by the models are not directly comparable to ETDs.To facilitate a comparison, both ETDs and CTDs were normalized to those of inert aerosol particles and Ind.Therefore, we compare relative rather than absolute values.Because most metals, including zinc (Zn) and copper (Cu), undergo atmospheric transport exclusively bound to particles [15] and are not degraded in the atmosphere, they can be used to estimate the atmospheric transport distance of particles.As Ind is also mainly bound to particles, we used it as a second proxy for particle-bound substances when normalizing travel distances.
We selected four ombrotrophic bogs along a ~400 km transect downwind of a large mining and smelting operation in Eastern Canada (Greater Sudbury area, Ontario; Figure 1) to derive the temporal trend of atmospheric deposition of PAHs, Zn and Cu.The smelters in the Greater Sudbury area have been a known source of air pollution, especially for copper and sulphur-based acid precipitation, to the surrounding area for the last century [16].Smelters emit PAHs and metals to the atmosphere as a consequence of wood or coal firing and roasting processes [15,17,18].Ombrotrophic peat bogs have previously been used as archives of the atmospheric deposition of several classes of organic and inorganic contaminants [19][20][21][22].The combination of continuous growth and limited bioturbation with atmospheric wet and dry deposition as the only source of contaminants makes them well suited to derive the historical atmospheric deposition trend of non-polar contaminants.Historical depositions were used here as they represent a long-term average.More importantly, the Greater Sudbury area was the dominant regional source in the past, whereas today there are many potential sources of PAHs.As deposition rates of pollutants correlate to their air concentrations [23], the decline of deposition rates over distance can be used to derive the ETD.

Sampling and Analysis
Peat cores were sampled at four bogs located downwind of the Greater Sudbury area (Canada; Figure 1) in summer 2008.Coordinates of the sampling locations are given in Table S1 in the Supplementary Material.At each site, three peat cores were sampled in hollows with a box corer and cut into 5 cm segments in the field.Samples were transported to the laboratory, frozen, freeze dried, and milled to a fine powder.For determination of PAHs, recovery standards (mixture of 7 deuterium-substituted PAHs) were added to 5 g (dry weight) of sample before extraction by accelerated solvent extraction with hexane.Extracts were purified by column chromatography (5 g silica, 3 g alox 15% deactivated) followed by size exclusion chromatography.The concentrations of 12 PAHs (10 EPA-PAHs: phenanthrene, fluoranthene, pyrene (Pyr), benzo[a]anthracene, chrysene, benzo[b]fluoranthene, benzo[k]fluoranthene, benzo[a]pyrene, indeno[1,2,3-cd]pyrene, and benzo[ghi]perylene (B[ghi]P), and additionally benzo[e]pyrene (B[e]P), and benzo[j]fluoranthene (B[j]F)) were determined as described in [22].Due to the incomplete chromatographic separation of B[b]F and B[k]F, the sum of the compounds is reported as B[b+k]F.Because these isomers have similar physico-chemical properties (less than 10% difference) [24], the ETD of B[b+k]F was compared to the CTD of B[b]F.A mass of 200 mg of the milled peat samples was digested by a microwave assisted acid dissolution technique [25] and Cu and Zn were analyzed by ICP-MS.The peat cores were dated using 210 Pb (t 1/2 = 22.26 years) by applying the constant rate of supply model [26].

Calculation of ETDs
The temporal resolution that can be achieved with peat cores is relatively low [22].Individual segments of 5 cm span periods of 20 to 30 years.In contrast to previous work [20] we did not calculate 10-year-averaged deposition rates, as this procedure might reduce the absolute magnitude of reconstructed deposition rates.Instead, we based our analyses on the segments showing the maximum PAH deposition rate.Such maxima have previously been suggested as chrono-stratigraphic markers in peat profiles by Turetsky et al. [27].In all peat cores PAH deposition rates peaked around 1880-1930, coinciding with the onset of intensive smelter operations in Sudbury.At that time, the smelters were one of the main sources of PAHs in the study region as the population was small and other sources, such as traffic, were not that relevant yet.
For each compound, the maximum deposition rates from each of the three replicate cores per bog were ln-transformed and linearly regressed against the distance of the sampling site from the Inco Superstack in Sudbury.In cases where the decline of the deposition rates was significant (p < 0.1) the negative reciprocal of the slope m (km −1 ) of the regression line corresponds to the ETD for the respective compound, and the ETD consequently corresponds to the distance where the air concentration of a compound declined by 63% (1/e).The error of the slope of the regression was used to estimate the uncertainty of the ETDs.The software package R [28] was used for these calculations.
Wind speed is an important parameter constraining atmospheric transport of contaminants, but reliable data on average wind speed are difficult to obtain for large areas and longer periods.
To facilitate the comparison of measured and calculated transport distances, we eliminated this source of uncertainty by normalizing the PAH transport distances to those of aerosol particles.As a proxy for the transport distance of particles, we used the average transport distance of the metals Cu and Zn or the PAH Ind. Subsequently, we refer to these normalized distance as r metal -ETD and r Ind -ETD.

Model Description
The Tool and ELPOS have both been designed to distinguish between chemicals with low and high LRTP [29].They are steady state fugacity models (level III, no equilibrium, inter-compartment transport [30]) that assess the environmental hazard of organic chemicals using the metrics of CTD and overall persistence (P OV ) under different emission scenarios [7].With only three generic compartments (air, soil and water), The Tool requires five input parameters for each compound: the air-water (K AW ) and octanol-water (K OW ) partition coefficients as well as degradation half-lives (t 1/2 ) in air, water, and soil.A Monte Carlo (MC) analysis routine estimates the 95% confidence interval of P OV and CTD.In this study we used dispersion factors of 5 and 10 for partition coefficients and half-lives, respectively, comparable to the factors of 3 and 10 used previously [8].ELPOS (Version 2.2) has eight compartments: air, water, sediment, four soil compartments (natural, agricultural, industrial, urban), and plants [5], necessitating seven chemical input parameters: K AW, K OW and the organic carbon partition coefficient (K OC ); as well as t 1/2 in air, water, soil, and sediment.Instead of estimating the K OC from K OW , we used the option to supply it directly as well as the K OA -approach to estimate the particle-bound fraction in air.
In addition to the CTD, both models also provide the fraction of a compound bound to particles in air as output information.As degradation rates for pollutants bound to particles are generally slow compared to the gas phase [31] and reliable estimates are difficult to obtain, both models assume that compounds associated with particles do not degrade [29].
Some of the default environmental input parameters of both models were adjusted to reflect the characteristics of Southern Ontario (fraction of model region covered by water = 10%, average water depth = 10 m, yearly average precipitation = 1008 mm, average wind speed = 13.1 km•h −1 taken from Canadian Climate Normals for North Bay-A, WMO ID: 71731, period 1971-2000 [32]).More detailed information on the model input parameters is provided as Supporting Material (Tables S2 and S3).Three temperature scenarios reflect the annual average (4 °C), summer (13 °C), and winter (limited to 0 °C as partition behavior to frozen compartments might differ and is not taken into account by the model calculations) conditions in Southern Ontario [32].
Partition coefficients (K OW and K AW ) at 25 °C for the PAHs were taken from Ma et al. [24].The K OC values at 25 °C were estimated using a poly-parameter linear free energy relationship (pp-LFER) [33] and solute descriptors from Absolv online [34].While ELPOS has a built-in routine to adjust these coefficients to a user-supplied temperature, The Tool requires the user to supply temperature-adjusted input information.The K AW values were adjusted using the temperature dependent pp-LFERs from Goss et al. [35]; K OA values from Ma et al. [24] were adjusted using enthalpies of vaporization, derived as described in Goss and Schwarzenbach [36].Finally, temperature-adjusted K OW were obtained as the product of K AW and K OA .We preferred predicted values since consistency between the properties of different PAHs is more important than accuracy of single properties.Degradation half-lives in air reflecting the reaction of PAHs in the gas phase with hydroxyl radicals were estimated with the Atmospheric Oxidation Program (AOPWIN v1.92) using default settings for a 12-h day as done by others [10,11].However, AOPWIN does not estimate activation energies, which are necessary to calculate degradation rate constants at different temperatures [29].Temperature dependent rate constants for the reaction with hydroxyl radicals are only available for Phen and Flt [37].The temperature dependence is so small that in case of The Tool degradation of all PAHs in air can be assumed to not depend on temperature [29].However, the variability in daytime length and OH radical concentration with season was taken into account.Therefore we used the averaged OH concentrations calculated for the atmospheric boundary layer (1000 and 900 hPa) by Spivakovsky et al. [38].Average daylight hours were taken from [39].Degradation half-lives in water, soil and sediment were estimated with EPISuite v4.0 [40] and retrieved through the Level III Fugacity Model implemented in this software.In ELPOS degradation rates are temperature adjusted by an implemented routine.By setting the temperature to 25 °C and using the manually adjusted values calculated for the Tool the implemented temperature adjustment was bypassed.
Whereas the two models were not designed to estimate CTDs of metals, it is possible to use them to calculate the CTD of exclusively particle-bound substances.As the metals for which we were able to determine ETDs (Cu, Zn) do not partition into the gas phase, they can be considered surrogates for particle-bound organic contaminants.The CTDs for such particle-bound contaminants was estimated by assigning partition coefficients that ensured almost exclusive partitioning to particles (log K OW = 10, log K AW = −9, yielding log K OA of 19).In addition, we assigned a very long environmental half-life (t 1/2 = 10 10 h).Additional information on this choice of parameters is available as Supporting Material (Figure S1).All model input parameters are listed in Tables S4 and S5.Table S6 summarizes the temperature-adjusted chemical properties.As for the ETDs, relative CTDs were calculated by normalizing the CTDs of the PAHs to the CTD of Ind (r Ind -CTD) and of exclusively particle-bound compounds (r particle -CTD).
The sensitivity of the CTD to changes in input values was tested by a sensitivity analysis as described in Wania and Dugani [41].A sensitivity of 1 indicates proportional changes of input values and CTDs.A sensitivity <1 indicates that the CTD does not respond linearly to changes of the input value.We analyzed the sensitivity of temperature (ELPOS only), wind speed, height of the air compartment, aerosol deposition velocity, precipitation rate, rain scavenging ratio, half-life times in all compartments, and partition coefficients for the 4 °C scenario by varying the values by 10%.

Empirically Determined Travel Distances
The maximum deposition of PAHs and the two metals was measured in peat segments dated to the period 1880-1930 (Figure S2 in Supporting Material, data given in Table S7).The increase in deposition around 1880 coincides with the onset of intense smelting operations in Sudbury.In most peat profiles, we observed declining deposition rates from 1920 onwards, when blast furnaces and roast yards were replaced by multi-hearth roasters and reverberatory furnaces for smelting [42].More detailed data and interpretation is available in Thuens et al. [22].The maximum deposition rates of Chry, B[b+k]F, B[a]P and Ind, as well as of Cu and Zn declined exponentially (p < 0.1) with distance from Sudbury (Figure 2).Together with the prevailing westerly winds, this finding is a powerful indication that the dominant source of atmospheric contamination to the studied bogs during the time of peak deposition was located in the Greater Sudbury area.For Phen, Flt, Pyr, and B[a]A, this decline was not significant (see Figure S3 in Supporting Material); for B[ghi]P, B[e]P, and B[j]P no ETD could be calculated, as concentrations were below the level of quantification in several samples.Therefore, these compounds are not considered in the following discussion.
The ETDs of PAHs, Cu, and Zn are shown in Table 1.They range from 173 to 321 km with relative uncertainties between 26% and 46%.The ETDs of the two metals were shorter than those of the PAHs.The discrepancy between the two metal ETDs reflects the overall uncertainty of our approach.Their average ETD (198 ± 45 km) was used to calculate the r metal -ETD of PAHs.Using data from several studies reporting metal deposition rates at variable distances from smelters [43][44][45], we calculated ETDs for Cu, Zn, Pb and Ni between 21 km and 67 km.Average uncertainty was 17% (Table S8 in Supporting Material).These ETDs are considerably lower than those we have obtained here.In these previous studies, samples were taken within a distance of up to 125 km from the smelters.It is known that metals such as Ni, Zn, and Cu are associated with both smaller and larger aerosol particles in the atmosphere [46].As larger particles are deposited more rapidly than smaller ones [47], large particles will dominate the deposition close to the emission source, while fine particles will be deposited at longer distances.This might cause divergent ETDs depending on the spatial scale and resolution of a particular study: the apparent ETD may be short in the immediate vicinity of the smelter, where deposition of coarse particles dominates, and long at larger distances, where mainly fine particles are deposited.This interpretation would explain the longer metal ETDs observed in our study, as the distances of the sites we analyzed are outside the range influenced by larger particles.As coarse particles are not subject to atmospheric long-range transport, we conclude that the metal ETDs determined in this study are more representative of the travel distances of the particles relevant for long-range transport of organic chemicals.

The Tool
With the default settings of The Tool, CTDs up to 2860 km were calculated for particles at 4 °C.After adjusting the settings and input parameters to the specific conditions in the study region the estimated CTDs were lower by a factor of 7.5.This reduction is mainly caused by changes in the partition coefficients and changes in size of the individual compartments, while wind speed contributed to a minor extent.The 95% confidence intervals of the absolute CTDs were larger for g-PAHs than p-PAHs (e.g., factor of 46 for B[a]A, 40 for Phen, 9 for Chry and 2 to 3 for the p-PAHs) indicating a higher uncertainty of the CTD for B[a]A and the g-PAHs.Details are summarized in Table S9 (Supporting Material).The confidence intervals cover ranges from almost 10 to 5500 km, indicating the high uncertainty of the model results.
Using The Tool, the CTD of particles was longer than the CTDs of the PAHs except for the CTD of Phen and Flt in winter.Consequently, most r particle -CTDs are <1 (Figure 3A).The CTD for Phen in winter was greater than all other CTDs by a factor of 4. The p-PAHs had almost the same CTD as the particles, while those of g-and c-PAHs were substantially shorter.B[a]A had the shortest CTD.As the CTD of Ind was almost the same as the one of particles, r Ind -CTDs were equal to the r particle -CTDs.
The CTD of particles was identical for the three temperature scenarios.The CTDs of all PAHs decreased with increasing temperature, which can be attributed to a higher proportion of these compounds being in the gas phase at higher temperatures and the gas phase degradation half-lives in air being shorter because of longer days and higher OH concentrations in summer, when temperatures are higher.Conversely, the high CTD of Phen and Flt for the 0 °C scenario is caused by the long half-lives in air, which are a result of short days and low OH concentrations in winter.The Tool predicted the c-PAHs to be mainly particle-bound in winter and equally distributed between gas and particle phase in summer (13 °C), whereas the p-PAHs were almost exclusively particle-bound at all temperatures.This is in good agreement with measured gas-particle distributions [14].
The CTDs of all studied compounds were directly proportional to wind speed (sensitivity of 1, Figure 3B).For the g-PAHs, half-life in air was the only other parameter with high sensitivity.The CTDs of the less volatile PAHs were sensitive to the particle deposition velocity, the precipitation rate, and the height of the air compartment while the sensitivity to the half-life in air was lower (Figure 3B).This finding reflects the importance of the removal of c-and p-PAHs from the atmosphere by dry or wet deposition.

ELPOS
With default settings, CTDs from 295 km (Flt) to 700 km were calculated with ELPOS (particles, 4 °C).With region-specific settings, CTDs between 142 km (B[a]A) and 517 km (particles) were estimated for the annual average temperature of 4 °C (Table S10 in Supporting Material).This reduction was caused by a lower height of the air compartment and a higher precipitation in the study area.The r particle -CTDs for all PAHs covered the range from 0.10 to 4.6 under the three temperature scenarios (Figure 3C).B[a]A had the shortest r particle -CTD due to its high degradation rate in air in summer, Phen had the longest in winter.
The temperature dependence of the chemical properties was substantial, e.g., K AW were up to 90% lower and half-lives 200% longer in the winter scenario compared to summer (Table S6).This effect had a high impact on the calculated CTDs.With rising temperatures, the CTDs of all PAHs diminished due to increased reactive loss from air.The relation between temperature and CTD was almost linear, with coefficients of determination for the regression of CTD against temperature >0.98 (Figure 3).The effect of temperature was largest for B[a]A and B[a]P and the CTD in winter double that in summer.
The particle-bound fraction of B[a]P and Ind estimated by ELPOS was less than 100% (given in Table S10 in Supporting Material) for all temperature scenarios and decreased with increasing temperature.This contrasted with measurements in Southern Ontario [14] where these PAH could not be detected in the gas phase.The overestimation of the gaseous fraction by ELPOS might cause an overestimation of B[a]P and Ind removal by degradation.
The estimated CTD of all compounds was directly proportional to wind speed (Figure 3D).The half-life of the volatile compounds in air was another highly sensitive parameter, whereas the particle-bound substances were more sensitive to precipitation and the height of the air compartment.The dry particle deposition velocity only substantially influenced the p-PAHs and particles.The sensitivity of the CTD against changes in temperature was ≤0.21 and similar regarding all compounds other than B[b]F and particles.

Model Comparison
When the same degradation half-lives and temperature-adjusted partition coefficients were used in both models, i.e., the internal temperature correction of ELPOS was not used, the CTDs calculated by both models agreed reasonably well (Figure 4A) and were significantly correlated (R 2 = 0.92, p < 0.00015).S11 and S12; the lines in sub-figures B and D are shown for visual support in identifying the individual parameters).The assignment of PAHs to groups according to their atmospheric partition behavior is indicated in brackets.depovelo = dry particle deposition velocity.
However, most ELPOS-CTDs were above the 1:1-line, i.e., ELPOS systematically estimated longer CTDs than The Tool.CTDs normalized to that of the particles agreed very well between both models, with all compounds falling close to the 1:1-line (Figure 4B).Clearly, normalization eliminates systematic difference between the models.The agreement for p-and c-PAHs was even better after normalizing the CTDs to that of Ind (Figure 4C).
When the in-built temperature-adjustment of ELPOS was used, the agreement between the models diminished.On the one hand, a substantial increase of CTDs calculated using ELPOS was observed for Phen and Flt (Figure 4A) with the temperature adjustment.This finding was not surprising as the temperature-adjustment of the atmospheric half-lives in ELPOS substantially increased the half-lives with decreasing temperatures (for 0 °C the increase was 200%, see Table S6), while the half-lives were not adjusted in The Tool.Moreover, the sensitivity of the two gas-phase PAHs Phen and Flt against changes in half-lives was high in both models (Figure 3B,D).For this reason, these two PAHs should indeed be most affected by the discrepancy of the degradation rates in the gas phase.ELPOS applies a default activation energy of 30 kJ•mol −1 for the temperature adjustment of the atmospheric degradation rate for all compounds.Using 10 kJ•mol −1 instead, the increase of t 1/2 of Phen in air at lower temperatures was reduced by 46%.The experimentally determined activation energy of Phen and Flt is less than 10 kJ•mol −1 [37].Therefore, the temperature adjustment in ELPOS might underestimate atmospheric degradation of PAHs at low temperature.On the other hand, the CTDs of the less volatile PAHs generally decreased when the implemented temperature adjustment was used.This is because the default temperature dependence of the partitioning coefficients implemented in ELPOS was larger than the dependence obtained manually, i.e., the PAHs partitioned more into the particulate phase in ELPOS than this was the case for the adjusted input values we used in The Tool.When comparing the r particle -CTDs, the pattern was basically the same (Figure 4B).With the built-in adjustment, however, the r particle -CTDs of the less volatile PAHs calculated by ELPOS were now consistently smaller than those with the manual temperature adjustment.Overall, the major discrepancies between the two models can be seen for the more volatile compounds, while the r particle -CTDs for the less volatile PAHs were comparatively consistent.This finding is important for the comparison with measured ETDs, as these are only available for the less volatile compounds and thus for compounds with a lower overall deviation between the two models.
The sensitivity of both models against changes in the individual input parameters was similar for the g-PAHs, and only wind speed and t 1/2 in air had a strong influence.However, while this half-life still strongly influenced the CTD of the c-PAHs in ELPOS, i.e., sensitivity was close to 1, the influence was substantially lower in The Tool (Figure 3B,D).This result can be attributed to a substantially lower fraction of B[a]A and Chry in the atmospheric gas phase in The Tool (20%) than in ELPOS (84%) (see Supporting Material for details).As Su et al. [14] found them to be 80% in the gas phase in summer and 40% in winter, The Tool obviously under-predicts the gaseous fraction.
The sensitivity of these compounds against K AW in The Tool was higher than that of the other PAHs and than found in ELPOS.This observation is most likely due to its sensitivity to K OA , which was governing the gas/particle partitioning.Sensitivity to K OA , which was calculated from the K AW , thus indicates which compounds primarily partition from gas to particle phase and are crucially affected by the partitioning coefficient.In the Tool the c-PAHs, whereas in ELPOS the p-PAHs were affected.For the predominantly particle-bound PAHs, processes determining the fate of atmospheric particles, such as precipitation and dry particle deposition velocity, became more influential in both models, although the dry deposition velocity was more influential in The Tool.

Comparison of CTDs and ETDs
The absolute ETDs of PAHs derived from the peat transect were between 281 and 321 km, whereas the absolute CTDs with regional settings were between 44 and 480 km (average temperature scenario, 4 °C).The latter are comparable to the CTD calculated for aldrin (~200 km (The Tool) and 95 km (ELPOS), using the input values implemented in both models), a pesticide banned by the Stockholm Convention on Persistent Organic Pollutants [48], and which has-compared to most banned substances-a relatively low LRTP [49].
The direct comparison between measured ETDs and calculated CTDs was impeded by the fact that the Greater Sudbury area is a point source on the regional scale.Even without any degradation and deposition, the PAH concentrations would thus decline with distance from the source because of radial dilution.Such a dilution is not implemented into either model, and instead a line source and wind that is blowing in a constant direction are assumed.However, by normalizing the measured and calculated transport distances to a common reference, i.e., non-reactive particles, the comparison between measurements and model becomes feasible as the effects of radial dilution, wind speed, and a changing wind direction are eliminated.Moreover, when using a particle-bound substance as a reference, differences in parameters that affect the deposition of atmospheric particles between the real environment and the model settings are eliminated as well, and the relative comparison of a given set of compounds in terms of atmospheric travel distance is possible.
As both models are intended to calculate the transport potential and to compare compounds relative to each other, rather than to reproduce real transport distances in the environment, it is more meaningful to compare relative CTDs with relative ETDs instead of absolute values.This comparison is illustrated in Figure 5.When normalized to particles, the r particle -CTDs of both models agreed very well, and the r particle -CTDs of the three p-PAHs were close to 1, whereas the r particle -CTD of Chry was substantially lower.However, the measured r metal -ETD for all four compounds was much higher and around 1.5 for all compounds.Thus, the agreement between measured and modeled relative transport distances was not very good.This outcome might be attributed to a different atmospheric behavior of the particles that metals and PAHs are associated with, or to some inaccuracy of the estimated particle-CTDs.Although it is very likely that PAHs and metals in historic samples shared the same dominant source in the Greater Sudbury area, they might be associated with different aerosol size classes.The particles contained differ in their atmospheric transport behavior and metals associated with somewhat larger particles than PAHs may consequently be removed from the atmosphere faster.Alternatively, PAHs might not be as reactive as the models assume.The importance of these factors for the discrepancy between r metal -ETDs and r particle -CTDs of PAHs cannot be decided based on the currently available knowledge.To further reduce uncertainty, the CTD of PAHs can be normalized to the travel distance of Ind (or any other of the p-PAHs) instead of metals (Figures 4C and 5B), thereby referencing the transport distance to a compound with identical emission behavior and that is association with particle of similar size.This procedure yielded a very good agreement between measurement and model for B[b]F and B[a]P, while the estimated r Ind -CTDs for Chrysene (ELPOS: 0.38; The Tool: 0.42) were much lower than the actually measured r Ind -ETD (0.94).
The finding that the measured r-ETD of Chry was almost identical to the r-ETDs of the other three PAHs was unexpected, as Chry is substantially more volatile (see Supporting Material, Table S4 for partition coefficients).Both Chry and Ind have similar atmospheric half-lives but Chry is the compound with the higher volatility and it should thus be removed faster if the assumption is correct that particle-bound PAHs are degraded much slower than gas phase PAHs [29].The discrepancy between r-ETD and r-CTD can have two reasons: (i) The volatility of Chry might be overestimated in the models which assume an ideal partitioning behavior of PAHs.For 0 °C ELPOS assumed 24% of Chry to be bound to particles, while Su et al. [14] found about 60% of Chry bound to particles in winter.Harner and Bidleman [50] supposed that PAHs might be bound within particles during their formation in combustion processes and thus might be non-exchangeable with the atmosphere.Moreover, PAHs initially present on the surface of particles may become less accessible for reactions, partitioning and volatilization due to a coating of the particles by secondary aerosol components [17,31]; (ii) The actual atmospheric half-life of Chry might be much longer than that used for estimating its CTD.However, in order to reach a r Ind -CTD equal to the r Ind -ETD, in ELPOS a nearly 10-fold increase in t 1/2 would be necessary, while it is even impossible to reach such a CTD in The Tool by increasing the atmospheric half-life.
Two previous studies aimed at comparing measured and calculated transport distances of organic pollutants.Shen et al. [10] found a good agreement between ETDs of pesticides with high LRTP determined from latitudinal gradients of air concentrations and CTDs calculated by ELPOS and two other models.In contrast, ETDs for pesticides derived from water concentrations in North American lakes were tenfold higher than the CTDs calculated by four models.This difference could be reduced by lowering the hydroxyl radical concentration and introducing dry weather conditions [11].Unfortunately, both studies did not measure metals or other predominantly particle-bound reference substances to allow for the calculation of relative travel distances.

Conclusions
In this study the comparison between modeled and measured transport distances was limited to only a few compounds, as the peat profiles did not reveal statistically significant ETDs for more PAH congeners.Keeping this in mind, several general conclusions can be drawn: -The absolute CTDs of PAHs estimated with the two models agreed reasonably well for predominantly particle-bound congeners, while discrepancies were observed for the more volatile congeners (Figure 4).This could be attributed to the temperature-adjustment of input parameters (in particular, K AW and half-life in air) in ELPOS, which relies on generic assumptions, most importantly with respect to the temperature dependence of the atmospheric degradation half-life.Since the CTDs of rather volatile compounds in both models are highly sensitive to the half-life in air (Figure 3), the quality of the predictions would improve with better data on the temperature dependence of the atmospheric degradation half-lives of PAHs and the implementation of this dependence in both models.-Peat archives sampled from bogs at variable distance from a single dominant point source can be used to experimentally determine travel distances of pollutants.-Particles deposit onto peat and, therefore, peat archives can provide information about the relative transport distances of particle-bound compounds.Most other passive samplers collect only compounds in the gas phase.-Data from peat archives are suitable to experimentally evaluate predicted CTDs of PAHs if data are normalized to a reference compound with similar behavior, even if the source characteristics (point vs. linear source) differ between model and reality.However, in the present study, the results from the peat archive limited the number of compounds available for such an evaluation.This conclusion might also apply to other non-polar atmospheric contaminants measured in peat archives.-The underestimation of the r particle -CTD of Chry might be due to a "non-ideal" partitioning behavior, i.e., Chry associated to aerosol particles might not freely partition to the gas phase.More information on such processes is needed to verify this hypothesis experimentally.Such a process could also apply to other contaminants that are already particle-bound during their emission.Compounds predominantly emitted as gaseous compounds would be less affected.

Figure 1 .
Figure 1.Locations of sampled peat bogs in Ontario, Canada.

Figure 2 .
Figure 2. Average maximum deposition rates of (A) PAHs and (B) metals as a function of distance from the Greater Sudbury source area.The error bars represent the standard deviation of the three replicate peat cores.

Figure 3 .
Figure 3. Characteristic Travel Distance (CTD) of selected PAHs normalized to the CTD of particles estimated with (A) The Tool; and (C) ELPOS for the three temperature scenarios, and sensitivity of the CTD at 4 °C to the model input parameters estimated with (B) The Tool; and (D) ELPOS (note: data not presented here available in Supporting Material, TablesS11 and S12; the lines in sub-figures B and D are shown for visual support in identifying the individual parameters).The assignment of PAHs to groups according to their atmospheric partition behavior is indicated in brackets.depovelo = dry particle deposition velocity.

Figure 4 .
Figure 4. Comparison of characteristic travel distances predicted by The Tool and ELPOS for the average temperature scenario (4 °C).Results from ELPOS were obtained both with the built-in temperature correction of partition coefficients and reaction rates (ELPOS auto T-adj) and with the parameters that were temperature-adjusted manually as in The Tool (ELPOS man.T-adj).(A) Absolute CTDs; (B) CTDs relative to the particle-CTD; (C) CTDs relative to Ind.The 1:1-line is shown as dotted line.

Table 1 .
Empirical travel distances (ETDs) with standard error (SE), significance level (p), coefficient of determination (R 2 ), and relative experimentally determined travel distances (r metal -ETD and r Ind -ETD) with SE.