A passive dosing method to determine fugacity capacities and partitioning properties of leaves

The capacity of leaves to take up chemicals from the atmosphere and water in ﬂ uences how contaminants are transferred into food webs and soil. We provide a proof of concept of a passive dosing method to measure leaf/polydimethylsiloxane partition ratios ( K leaf/PDMS ) for intact leaves, using polychlorinated biphenyls (PCBs) as model chemicals. Rhododendron leaves held in contact with PCB-loaded PDMS reached between 76 and 99% of equilibrium within 4 days for PCBs 3, 4, 28, 52, 101, 118, 138 and 180. Equilibrium K leaf/PDMS extrapolated from the uptake kinetics measured over 4 days ranged from 0.075 (PCB 180) to 0.371 (PCB 3). The K leaf/PDMS data can readily be converted to fugacity capacities of leaves ( Z leaf ) and subsequently leaf/water or leaf/air partition ratios ( K leaf/water and K leaf/air ) using partitioning data from the literature. Results of our measurements are within the variability observed for plant/air partition ratios ( K plant/air ) found in the literature. Log K leaf/air from this study ranged from 5.00 (PCB 3) to 8.30 (PCB 180) compared to log K plant/air of 3.31 (PCB 3) to 8.88 (PCB 180) found in the literature. The method we describe could provide data to characterize the variability in sorptive capacities of leaves that would improve descriptions of uptake of chemicals by leaves in multimedia fate models.


Introduction
Leaves are a sink for semi-volatile organic chemicals in the atmosphere and a vector for transfer of air pollutants to soil with litter 1,2 and into the terrestrial food chain. 3 Studies comparing forested areas with nearby clearings have found that deposition of some organic pollutants to forested soils was signicantly higher because of the shedding of leaves and waxes loaded with chemicals. 2,4 Enhanced deposition of chemicals under forest canopies is known as the forest lter effect and is thought to be most pronounced for chemicals with octanol/air partition ratios (K OA ) between 10 7 and 10 11 and air/water partition ratios (K AW ) greater than 10 À6 . 5 Multimedia models indicate that the forest lter effect can have a substantial impact on the fate of chemicals in the environment, in particular in boreal regions characterized by cool temperatures and dense forest coverage. 6,7 Understanding the forest lter effect requires knowledge of the fugacity capacities of leaves (Z leaf ), which are important parameters in multimedia fate and transport models such as CoZMo-POP 8 and BETR. 9 By denition, leaf/air and leaf/water partition ratios (K leaf/air and K leaf/water ) are the ratios Z leaf /Z air and Z leaf /Z water respectively. Available measurements indicate that plant/air partition ratios (K plant/air ) can differ by up to 3 orders of magnitude between plant species. 2,5,10,11 An interesting observation is that the partition ratios are plant species-dependent in the case of whole leaves, but one regression curve could be tted to a wide range of literature data for partitioning to isolated cuticles, which are the waxy outermost part of the leaves. 12 There is a need for consistent data on how species diversity inuences the partitioning of semi-volatile pollutants to leaves, and how partitioning to whole leaves relates to partitioning to isolated cuticles.
In this study we illustrate the use of sheets of polydimethylsiloxane (PDMS) as passive dosing devices for determining leaf/PDMS partition ratios (K leaf/PDMS ) of whole leaves for polychlorinated biphenyls (PCBs). PDMS is a convenient reference phase to calculate Z leaf and subsequently K leaf/air and K leaf/water . 13 We further calculate cuticle/water partition ratios from our measurements that are consistent with measurements derived from isolated cuticles reported in the literature.

Material and methods
All native and labeled standards used in this study were purchased from Larodan (Solna, Sweden). Vendor specications and purities for the solvents can be found in Table s1. † PDMS sheets of the type SSP-M823 with a thickness of 610 mm were purchased from Shielding Solutions Ltd (Essex, U.K.), the European distributor for Specialty Silicone Products, Inc.

Experimental setup
PDMS disks 18 mm in diameter were punched from the sheets, soaked in 200 mL of methanol for at least a week and then soaked overnight in acetone. The disks were loaded with PCBs in a ask by adding 100 mL of methanol spiked with 1 mg of each of the PCB congeners 3, 4, 28, 52, 101, 118, 138 and 180 dissolved in isooctane. In total 8 mL of isooctane was spiked to the loading solution. The PDMS disks and the loading solution were le to equilibrate in the dark for a week, then 50 to 100 mL aliquots of water from a MilliQ purication unit (Merck, Darmstadt, Germany) were successively added over the course of several hours up to a nal volume of 400 mL to force the chemicals into the PDMS as described by Birch et al. 14 The PDMS disks were then transferred to MilliQ water and le for one hour to remove traces of methanol, aer which they were dried overnight in a fume hood to remove traces of isooctane. No measurable amounts of isooctane were found in the PDMS disks aer this drying step in a trial study.
Rhododendron leaves (R. ponticum L.) were picked from a bush on the Frescati campus of Stockholm University. Rhododendron leaves were chosen because of their year-round availability and their rigid structure. Care was taken to select leaves that had similar exposure to sunlight. All leaves were wiped clean using wet paper tissues and blotted dry. In total, 15 leaves were collected. Disks 18 mm in diameter were punched from the leaves and shuffled randomly to avoid a possible bias of the leaf source. A subset of 17 disks was dried in the oven at 60 C for 3 days to determine the dry weight (DW) of the leaves. The volume of the leaf disks was measured by volume displacement of water using a separate batch of 54 disks.
Our passive dosing apparatus is similar to that of Kim et al., 12 which is a modied version of the setup used by Mayer et al. 15 and Trapp et al. 16 to measure the kinetics of chemical transport through different matrices. Leaf disks were sandwiched between two loaded PDMS disks (Fig. 1) and the assembly was pressed together with glass plates held by metal clamps. Blanks consisting of cleaned PDMS disks in contact with leaves were kept in a closed jar to avoid cross-contamination from the loaded PDMS. Leaves were sampled in triplicates before the leaves and PDMS were brought into contact (at 0 hours) and aer 2, 4, 8, 24, 48 and 96 hours contact time. The temperature and relative humidity in the fume hood were not monitored during this experiment but were later measured to be 20 AE 0.5 C and 32.1 AE 4.8%.

Measurement of total solvent-extractable organic matter
The extractable organic matter of a batch of Rhododendron leaves was measured using the modied Jensen extraction method (mod. 2). 17 The extracts were then le in the fume hood until they reached a constant weight (within 24-48 h).

Extraction of PCBs from leaves
Aer terminating exposure, the leaf disks were immediately transferred to pre-weighed microcentrifuge tubes containing 1 mL of acetonitrile, spiked with 100 ng of 13 C labeled PCBs 28, 101, 138 and 180, and a 5 mm diameter stainless steel ball was added. The disks were then homogenized using a Mini G 1600 tissue homogenizer (SPEX, Stanmore, U.K.) and extracted using QuEChERS (Quick-Easy-Cheap-Effective-Rugged-Safe, Phenomenex, California, U.S.A.). A detailed description of the method can be found in Text s1. † The nal extracts in acetonitrile were solvent-exchanged to isooctane, reduced to 0.5 mL using a gentle stream of nitrogen, and spiked with 50 ng of PCB 53 before being analyzed. The concentrations in the leaves are given on a dry volume basis (ng m À3 ), with the volume of a batch of dry leaves measured by the displacement of water.

PDMS extraction
The top and bottom PDMS disks of each sandwich were pooled, transferred to a test tube containing 4 mL of acetone and 100 ng of 13 C labeled PCBs and le overnight. The extract was then transferred to a new vial and solvent-exchanged to isooctane. The nal extract was reduced to 0.5 mL using a gentle stream of nitrogen, spiked with 50 ng of PCB 53 and then diluted by a factor 10 before analysis.

Instrumental analysis
The samples were analyzed on a gas chromatograph (Trace 1310, Thermo Scientic, U.S.) coupled to a single quadrupole mass spectrometer (ISQ LT, Thermo Scientic, U.S.). The ions used to quantify each of the analytes can be found in Table s2 † and an overview of the GC/MS program in Text s2. † The analyte

View Article Online
concentrations were quantied using the labeled internal standards in all samples except the leaf samples collected between 0 and 48 hours. In these samples we used PCB 53 as an internal standard and corrected the concentrations with the average recoveries of the labeled standards determined in the 96 hour samples.
Quality assurance/quality control All the blanks in this study (n ¼ 8 for the leaves to 13 for the PDMS) underwent the same extraction procedure as the actual samples with at least one blank processed at every sampling time point. Method quantication limits (MQLs) were calculated using the average concentration of the analyte in the blanks + 9 times the standard deviation (Table s3 †). When the concentration in the blank was too low to be measured, the lowest concentration in the calibration curve (4.84 AE 2.01 pg mL À1 , stddev.) was used. In those cases where analytes in the blanks were only detected in one of the replicates, we calculated the MQL conservatively as ten times the concentration found in the blanks. As the concentration in the blanks did not show any correlation with time, all blanks were combined to derive MQLs.

Results
Unless specied otherwise, all uncertainties provided in this section are standard deviations and all concentration ratios (C leaf /C PDMS ) are given on a dry volume basis.
Quality assurance/quality control The MQLs ranged from 27 to 197 ng g À1 DW for the leaves and 10 to 886 ng g À1 for the PDMS (Table s3 †). All the samples in this study had concentrations > MQL. The recoveries of the labeled PCBs from the QuEChERSbased leaf extraction ranged from 69 to 80% (Table s4 †) which is 10% lower than the recoveries reported for a similar extraction method for PCBs from catsh. 18 The recoveries of the labeled internal standards from the PDMS ranged from 105% to 120% for the extractions from blank PDMS and from 135% to 155% for the extractions from loaded PDMS. One possible explanation for the high recoveries from loaded PDMS is the dilution by a factor 10 of the loaded PDMS samples before analysis. Peak areas of the volumetric standard PCB 53 were roughly 20 times lower in diluted samples than in the standards and the undiluted samples. Therefore the apparently high recoveries may reect an artefact associated with extrapolation from our calibration series which is based on a series of replicates with a specic concentration of labeled internal standards and PCB 53. The recoveries of labeled PCBs from PDMS do not affect the concentration measurements used to determine uptake kinetics in the passive dosing experiments since all analytes were quantied using labeled internal standards.

Water loss from leaves
The leaves gradually dried out during the experiment (Fig. s1 †). While fresh leaves weighed on average 65 mg per disk, their average weight was reduced to 27 mg (a loss of 58%) at the end of the experiment, which is comparable to the DW (29 mg).

Leaf/PDMS partitioning
The concentrations of individual PCBs in the loaded PDMS disks varied by a factor of 2 to 4 between the lowest and highest measured values. No correlation was found between exposure time and concentration of the analytes (Pearson, P > 0.7).
Due to the variability in the PCB concentrations in the PDMS disks, the concentrations in leaves were normalized by calculating the concentration ratio between each individual leaf and the pooled PDMS disks it was exposed to. This calculation was done on a volume/volume basis.
A one-phase association curve tted the data well for all congeners, with R 2 values ranging from 0.86 to 0.91 (Fig. 2). This association curve was used to estimate the leaf/PDMS concentration ratio at equilibrium (K leaf/PDMS ). According to the curve ts, the chemicals had reached between 76% (AE17, PCB 4) and 99% (AE4, PCB 3) of equilibrium aer 4 days. The estimated time to reach 95% of equilibrium ranged between 3 and 6 days for all congeners except PCB 4 (8 days). Upon reaching equilibrium, between 2 and 6% of the analytes in each PDMS/leaf/PDMS system would have been transferred to the leaves (Table s5 †).
Z leaf (Table 1) was estimated by multiplying K leaf/PDMS with Z PDMS obtained from measurements of the PDMS/water partition ratios (K PDMS/water ) multiplied by Henry's law constants (H) from the literature measured at 25 C. 13,19,20 K leaf/water was then calculated by dividing Z leaf by Z water (K PDMS/water Â H) and K leaf/air by dividing Z leaf with Z air (1/RT ¼ 0.00040 mol m À3 Pa À1 , with R being the gas constant and T the temperature in Kelvin).

PDMS loading
The nal ratio of loading solution to water was two times higher in our experiments than in Birch et al. 14 Some of the variability in the PCB concentrations in the loaded PDMS may have been caused by using isooctane as solvent for the standards, as upon the addition of water to the loading solution some of the isooctane was forced into the PDMS. Furthermore, the addition of water caused the PDMS to oat on the loading solution, thereby exposing some disks differently than others. We therefore advise the use of more polar solvents such as methanol or acetone for the stock solutions of the chemicals and adding lower amounts of water, as suggested by Booij et al. 21 Leaf/PDMS partitioning The time required to approach equilibrium between the PDMS and the leaves was similar for the mono-to hepta-chlorinated congeners and corresponded well with those from a comparable setup for the uptake of PAHs in isolated cuticles. 12 That study showed that up to 4 days were needed to reach equilibrium between the cuticle and the PDMS for all but the highest molecular weight PAH congeners investigated in their study. Slower equilibration kinetics for the high molecular weight

View Article Online
PAHs were attributed by the authors to higher molar volumes of those compounds. It has recently been shown that there can be considerable differences in the mass transfer kinetics in these types of studies, possibly caused by a difference in pressure exerted on the system by the use of certain clamps or magnets to keep the setup together. 22 While unlikely to affect the leaf/PDMS partition ratios, it is possible that the mass transfer kinetics into the leaves were affected by the pressure on the so leaf tissue which could alter the structure of the leaves. No damage was observed on the leaves in this study in which we used glass sheets to distribute the pressure from the clamps more evenly across the leaves. While the variability between replicate measurements in this study is much larger than that observed in Mayer et al. 15 and Trapp et al., 16 it is comparable with that in the study of Kim et al. for isolated cuticles. 12 A quantitative comparison with this study is not possible however as the variability can only be estimated from Fig. 2 in their study. 12 The observation that the uptake conforms to rst-order kinetics and had similar rate constants for the different PCB congeners (Table 1) could be an indication that equilibrium with the whole leaf was not reached with this experimental setup. Plant uptake models based on a two-compartment approach have shown that there is a compartment in some leaves that can respond rapidly to changes in atmospheric concentrations and a second compartment which reacts more slowly. 23,24 Our rst-order kinetic uptake model assumes implicitly that there is only one compartment, and it ts the data well. We see no evidence of two-compartment uptake in our kinetic curves and all the ts have high R 2 values. Based on the literature data cited above, 23,24 it is possible that we measure kinetics of uptake into the fast responding compartment only.
The studies by Wild et al. 25 and subsequently Li and Chen 26 and Li and colleagues 27 have shown that phenanthrene can be transferred across the cuticle and into the epidermis within 24-48 hours which is within the timeframe of our experiment. Fig. 2 The leaf/PDMS concentration ratios (on a volume/volume basis) of the analytes over time. Extrapolation to equilibrium (K leaf/PDMS ) was done using a one-phase association curve in Graphpad Prism of the form Y ¼ Y max (1 À e (Àkt) ). The broken blue lines show the fitted curves while the dotted black lines show Y max ¼ K leaf/PDMS (AEstd error, shaded area). No difference was found between the measured k values (P > 0.15). Values for k and K leaf/PDMS are reported in Table 1. Table 1 Rate constants (k in h À1 ), log K leaf/PDMS and values for Z leaf , log K leaf/water and log K leaf/air for dry leaves. Values of k and Z leaf are given with standard deviations which show the uncertainty of the extrapolation using the nonlinear regression. For log K leaf/PDMS , log K leaf/water and log K leaf/air , the uncertainty was propagated from the input data and is given as a range of AE1 standard deviation of the measured value (in brackets). A more detailed description of the error propagation can be found in Text s3

View Article Online
In a recent study in our lab, 22 the use of the method presented in Mayer et al. 15 and Trapp et al. 16 to measure the mass transfer kinetics of PAHs and PCBs through leaves gave an indication that this process is very slow with mass transfer coefficients of 1.65 Â 10 À6 to 2.17 Â 10 À6 m h À1 for the transfer of uorene and phenanthrene respectively through leaves of a Hydrangea species while mass transfer through Rhododendron leaves was too slow to result in quantiable amounts of analytes penetrating through the leaf and accumulating in the acceptor PDMS on the other side within 48 hours. In combination with this study, it seems that while hydrophobic organic chemicals are quickly taken up by the leaf, their transfer through the entire leaf is slow. K leaf/water was converted to cuticle/water partition ratios (K cuticle/water ) by dividing by the volume fraction of the leaves that is cuticle (4.7%), assuming a cuticle thickness of 4.15 mm (average from 3 Rhododendron species). 28 In this case it meant applying a correction factor of 21.3 (1/0.047). Values for K cuticle/water were within a factor 2.5 of K cuticle/water estimated from regressions found in the literature for a variety of plant species (Fig. 3). 12,29 Fitting a regression curve through the entire set of values for log K cuticle/water (n ¼ 75) against the log K ow of the respective chemicals provides a good t (R 2 ¼ 0.95) and gives the impression that it is possible to estimate K cuticle/water from K ow measurements with the regression: Aside from K cuticle/water , there is a wide range of literature available on plant/air partitioning from both eld studies and models. Differences between plant species and the applied methods result in a range of 3 orders of magnitude of measured K plant/air (Fig. 4). Converting our K leaf/PDMS to K plant/air using partition ratios from the literature resulted in data points that lie within the range of existing data for clover, plantain, hawk's beard and yarrow measured using a fugacity meter, 10 and that differed by roughly an order of magnitude from data derived from deposition measurements, 2,5 or empirical regressions. 14 Recently, Vorkamp et al. 30 demonstrated that equilibrium sampling of PCBs in indoor air using silicone-coated baking paper is feasible. The K leaf/PDMS data from the present study can provide conversion factors for equilibrium sampling concentrations from silicone-based air sampling to equilibrium partitioning concentrations in leaves. More specically, the K leaf/PDMS data can be used (i) to estimate PCB concentrations in leaves based on equilibrium sampling measurements in air, (ii) for consistency checking between equilibrium sampling in air and leaf monitoring data or (iii) for checking the equilibrium status of PCBs in leaves based on parallel equilibrium sampling in air and measurements of concentrations in leaves. The fugacity capacity of leaves can also be used directly in studies of the uptake of pollutants in plants at contaminated sites 31 and in studies that use leaves as passive samplers for indoor air 32 or as a screening tool to assess spatial variability of semi-volatile chemicals in air. 33 To give lipid-normalized concentrations, the concentrations in the leaves were normalized to the leaves' total solventextractable organic matter (determined using the Jensen extraction) 17 as an indicator of lipid content. There was poor agreement between these so-derived lipid/PDMS partition ratios (K lipid/PDMS ) with measurements of K lipid/PDMS for olive oil and a wide range of animal lipids from the literature (Fig. s2 †). [34][35][36][37] The lower K lipid/PDMS values for the Rhododendron leaves are consistent with a certain fraction of the extractable organic matter not being available for partitioning. For instance, the crystalline nature of some cuticular waxes could make them less available for partitioning.
The solvent-extractable organic matter of leaves normally does not contain the depolymerizable lipids such as cutin (and for some species also cutan), which have been identied as potentially responsible for a major fraction of the sorptive capacities of leaves for HOCs. 38 If the depolymerizable lipids reached equilibrium in our experiments and were not extracted with other lipids, then the K lipid/PDMS values that we measure for the leaves would be highly overestimated. However, our measured K lipid/PDMS are below those found in the literature. A possible explanation is that the depolymerizable lipids are part of the slowly responding compartment observed in other studies, and did not reach equilibrium with the PDMS in our experiments. Another possible explanation for the poor agreement could be that the partitioning properties of the extractable organic matter of foliage and the other lipids studied in the literature are fundamentally different from each other.
Our experimental setup provides a straightforward way of measuring the fugacity capacity of leaves, which can then readily be applied to calculate K leaf/water , K leaf/air and other partition ratios of interest. While no observation of a twocompartment system was made, it is possible that our setup only reached equilibrium with the cuticle and the inner compartments of the leaf that are easy to access. The impact of the loss of water from the leaves on the fugacity capacity of these leaves for HOCs is most likely very limited. However, future studies should test if these results can be extrapolated to leaves with a higher water content than the Rhododendron leaves, which would require simple modications of the experimental setup such as keeping moist tissues underneath the setup, or keeping the setup in a closed environment with a beaker of water.