
==== Front
Environ Sci Pollut Res Int
Environ Sci Pollut Res Int
Environmental Science and Pollution Research International
0944-1344
1614-7499
Springer Berlin Heidelberg Berlin/Heidelberg

39129044
34460
10.1007/s11356-024-34460-x
Research Article
Modelling passive sampling of hydrophilic compounds under time-variable aqueous concentrations
http://orcid.org/0000-0002-9383-819X
Becker Benjamin benjamin.becker@bafg.de

1
Kochleus Christian 1
Spira Denise 1
Bachtin Julia 1
König Fabian 2
Meinecke Stefan 2
Möhlenkamp Christel 1
Booij Kees 3
1 https://ror.org/03kdvpr29 grid.425106.4 0000 0001 2294 3155 Federal Institute of Hydrology (BfG), Am Mainzer Tor 1, 56068 Koblenz, Germany
2 grid.425100.2 0000 0004 0554 9748 German Environment Agency (UBA), Schichauweg 58, 12307 Berlin, Germany
3 PaSOC, Greate Pierwei 25, 8821 LV Kimswerd, The Netherlands
Responsible Editor: Marcus Schulz

12 8 2024
12 8 2024
2024
31 39 5184451857
19 12 2023
19 7 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Passive sampling is a crucial method for evaluating concentrations of hydrophilic organic compounds in the aquatic environment, but it is insufficiently understood to what extent passive samplers capture the intermittent emissions that frequently occur for this group of compounds. In the present study, silicone sheets and styrene-divinyl benzene-reversed phase sulfonated extraction disks with and without a polyethersulfone membrane were exposed under semi-field conditions in a 31 m3 flume at three different flow velocities. Natural processes and spiking/dilution measures caused aqueous concentrations to vary strongly with time. The data were analyzed using two analytical models that account for these time-variable concentrations: a sampling rate model and a diffusion model. The diffusion model generally gave a better fit of the data than the sampling rate model, but the difference in residual errors was quite small (median errors of 19 vs. 25% for silicone and 22 vs. 25% for SDB-RPS samplers). The sampling rate model was therefore adequate enough to evaluate the time-integrative capabilities of the samplers. Sampler performance was best for SDB-RPS samplers with a polyethersulfone membrane, despite the occurrence of lag times for some compounds (0.1 to 0.4 days). Sampling rates for this design also spanned a narrower range (80 to 110 mL/day) than SDB-RPS samplers without a membrane (100 to 660 mL/day). The effect of biofouling was similar for all compounds and was consistent with a biofouling layer thickness of 150 µm.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11356-024-34460-x.

Keywords

Passive sampling
Diffusion model
Sampling rate model
Biofouling
Flow effects
Peak events
SDB-RPS
Silicone
German Environment AgencyFKZ: 3717 67 418 0 issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

Environmental monitoring of hydrophilic organic compounds in water can be done using grab sampling, continuous water sampling, and passive sampling. Emissions of hydrophilic compounds (for example polar pesticides and pharmaceuticals) often occur as peak concentration events, following rainfall and intermittent emissions of prescription drugs and illicit drugs (Gallé et al. 2020; Harman et al. 2011; Morosini et al. 2017; Ort et al. 2014). Passive sampling can be considered when the cost of high-frequency grab sampling and continuous water sampling is prohibitive. A current weakness of passive sampling of hydrophilic compounds is that the mechanisms of sampler-water exchange kinetics are not well understood, and that reported sampling rates (Rs) and sorption coefficients (Ksw) of these compounds show a large scatter, without having a clear relationship with flow velocity and temperature (Charriau et al. 2016; Harman et al. 2012).

Passive sampler calibration is commonly done by exposing samplers to constant aqueous concentrations (Cw) for time periods of typically 2 weeks, followed by data analysis with the sampling rate model1 Ns=mKswCw1-exp-RstmKsw

where Ns is the accumulated amount, m is the sampler mass, Ksw is the sorbent-water sorption coefficient in volume per mass units, and t is time. Rs is obtained from the initial slope of Ns vs. time, and Ksw from the plateau at long time scales (Huckins et al. 2006). The time-integrative window (TIW) is a useful parameter for summarizing the capability of passive samplers to yield time-integrated concentrations.2 TIW=mKswRs=mKswAko

where ko is the overall mass transfer coefficient, and A is the water-sampler exchange surface area. Sampling is essentially time-integrative when the exposure time is smaller than the TIW (Booij et al. 2007). The TIW is often multiplied by -ln(0.5) ≈ 0.7 to obtain the half-life time for equilibrium attainment, but sampler-water equilibrium is only defined when Cw is constant, and the factor 0.7 is of minor importance. Sampler designs can be manipulated to yield larger TIWs by choosing sorbents with higher Ksw, adding diffusion-limiting membranes that result in a smaller ko, and using a higher ratio of sorbent mass to surface area (Eq. 2). The use of diffusion-limiting membranes often results in lag times that originate from sorption of the target compounds to the membrane (Endo and Matsuura 2018; Vermeirssen et al. 2012).

The Rs model has also been used to model passive sampler response to peak concentration events, by evaluating if the Rs that is determined under constant Cw scenarios predicts the time-average Cw from repeated grab sampling (Bernard et al. 2018; Mazzella et al. 2008; Noro et al. 2019; Schreiner et al. 2020; Shaw and Mueller 2009). The overall conclusion from these studies is that the Cw from passive sampling and from repeated grab sampling agree within a factor of approximately 1.5, and that the presence of a polyethersulfone (PES) membrane improves the time-integrative capability of the sampler (Shaw and Mueller 2009).

Vermeirssen et al. (2012) used an empirical rate constant model for hydrophilic compound accumulation by POCIS and Chemcatchers. These authors then converted the rate constants to an equivalent Rs over an arbitrary time period of 30 d, using numerical methods. Schreiner et al. (2020) used a similar approach for Chemcatchers without a PES membrane. Analytical models for sampler-water exchange kinetics with time-variable Cw conditions have been missing so far, with the exception of the relatively simple case of Cw that varies linearly with time (Booij et al. 2003).

Tcaciuc et al. (2015) proposed to abandon Rs models because of their weak mechanistic basis, and to use diffusion models instead. These models are solutions of Fick’s second law of diffusion in the sorbent and allow modelling of partial rate control by sorbent, membrane (if present), and water boundary layer (WBL) (Booij 2021; Tcaciuc et al. 2015; Thompson et al. 2015). Analytical solutions of Fick’s second law have been applied in passive sampling research for constant Cw scenarios in infinite water volumes, and constant amount scenarios in finite water volumes.

Two approaches are available for applying diffusion models for scenarios with time-variable Cw and mixed rate control by sorbent, WBL, and membrane. Endo et al. (2019) used numerical integration of Fick’s second law for samplers with a PES membrane and a sorption phase of styrene–divinylbenzene reversed phase sulfonated polymers (SDB-RPS). Mikhailov and Özişik (1984, chapter 7.2) presented a series solution model that can be used for diffusive transport in a passive sampler sorbent, under conditions of time-variable Cw and partial rate control by the sorbent, but this model has not been applied so far in passive sampling research.

A present research challenge for passive sampling of hydrophilic compounds is to evaluate passive samplers for their capability of yielding time-integrated Cw estimates for realistic Cw scenarios, using sound mechanistic models. The aim of the present study was to assess passive sampler response to peak concentration events under semi-field conditions, to understand the effects of flow, biofouling, and presence/absence of a PES membrane on the exchange kinetics, and to evaluate the merits of a diffusion model and an Rs model for sampler-water exchange kinetics under time-variable Cw. The present study was carried out in three experimental phases (Table 1). The aim of Phase 1 was to evaluate uptake and release kinetics of silicone and SDB-RPS samplers, and to study the effect of flow velocity on the exchange kinetics of SDB-RPS samplers. The purpose of Phase 2 was to evaluate sampler response to three peak concentration events, including an evaluation of flow effects. Phase 3 aimed to study differences in exchange kinetics between SDB-RPS samplers with and without a PES membrane, and to evaluate the effect of biofouling on the exchange kinetics. Table 1 Summary of sampler deployment and retrieval. Spiking of Channel A was done at t = 1 day for Phase 1 and 3, and at t = 1, 3, and 6 days for Phase 2. U is the water flow velocity

Scenario	Treatment	
Phase 1	Spike Channel A at t = 1 day	
ChA	Deploy silicone and SDB-RPS samplers in Channel A at U = 9 cm/s at t = 0. Retrieve at regular time intervals up to 13 days	
ChA- > ChB	Deploy silicone and SDB-RPS samplers in Channel A at U = 9 cm/s at t = 0	
Transfer a subset of samplers from Channel A to Channel B at t = 2 days. Retrieve samplers up to 13 days	
Velocity effect I	U = 1 cm/s: deploy SDB-RPS samplers in Channel A at U = 1 cm/s. Transfer a subset to Channel B at t = 2 days. Retrieve at t = 8 days	
U = 18 cm/s: same procedure as for U = 1 cm/s	
U = 9 cm/s: use data from scenario ChA	
3 x [ChA <—> ChB]	Deploy silicone and SDB-RPS samplers in Channel A at t = 0, U = 9 cm/s	
Transfer samplers to Channel B at t = 2 days, and back to Channel A at t = 3 days	
Transfer samplers to Channel B at t = 4 days, and back to Channel A at t = 5 days	
Transfer samplers to Channel B at t = 6 days. Retrieve at t = 7 days	
Phase 2	Spike Channel A at t = 1, 3, and 6 days	
ChA	Deploy silicone and SDB-RPS samplers in Channel A at t = 0, U = 9 cm/s	
Decrease concentrations by partially draining/refilling Channel A after spiking	
Retrieve up to 17 days	
Velocity effect II	U = 1 cm/s: Deploy silicone and SDB-RPS samplers in Channel A, U = 1 cm/s. Retrieve samplers at t = 8 days,	
U = 18 cm/s: same procedure as for U = 1 cm/s	
U = 9 cm/s: use data from scenario ChA	
Phase 3	Spike Channel A at t = 1 day	
ChA	Deploy SDB-RPS samplers with and without PES membrane in Channel A at t = 0, U = 9 cm/s. Retrieve up to 18 days	
ChA- > ChB	Deploy SDB-RPS samplers with and without PES membrane in Channel A at t = 0, U = 9 cm/s	
Transfer samplers to Channel B at t = 2 days. Retrieve up to 18 days	
Fouling effect	Allow biofilm to develop on SDB-RPS samplers in Channel C for 2 weeks, and 4 weeks	
Deploy in Channel A at t = 0, U = 9 cm/s. Retrieve at t = 2 and 4 days	

Materials and methods

Materials

Translucent silicone (Specialty Silicone Products, SSP M823, thickness 250 µm) was obtained from Shielding Solutions, UK. Sheets were cut to a size of 55 × 90 mm (mass 1.35 g, exchange surface area 0.99 dm2) and were pre-extracted with ethyl acetate by soaking in solvent for 1 week, followed by 12 times 90 min microwave extraction with solvent replacement (Becker et al. 2020). The sheets were then transferred to a flask with methanol for 2 days, with solvent replacement after 1 day, followed by rinsing with bi-distilled water, drying with lint-free tissue, and storage at − 20 °C.

SDB-RPS extraction disks (thickness 660 µm, disk mass 314 mg, diameter 47 mm) were obtained from Affinisep, France. Ratios of porosity (φ) and squared tortuosity (θ) were measured as φ/θ 2 = 0.34, using the alabaster dissolution rate method (Booij et al. 2017). Disks were pre-extracted three times, by swirling gently with 20 mL methanol, followed by a 30 min rest period. Methanol was poured off and the process was repeated twice with distilled water. Disks were clamped between two stainless steel plates (bottom plate 100 × 70 mm, thickness 2 mm, upper plate 70 × 70 mm, thickness 2 mm, with a 40 mm diameter opening in the middle) (Figure S1-1). The effective sorbent mass was taken to be 314 mg × (40 mm/47 mm)2 = 227 mg. The exposed surface area was 0.126 dm2.

PES microporous membranes (Pall, Supor-200) had a thickness of 150 µm, nominal pore size of 0.2 µm, and diameter of 47 mm, and φ/θ 2 = 0.78 (Booij et al. 2017).

Exposure channels

Experiments were carried out in three channels from the artificial stream and pond system of the German Federal Environment Agency in Berlin-Marienfelde. Each channel consisted of two rows of 16 segments of 3 m length that were connected at the ends by U-turns, creating a closed loop with length of 104 m, width of 0.98 m, and water depth of 0.3 m (Figure S1-2). The channels were filled with water from a storage pond fed with manganese- and iron-free groundwater. The water volume in each channel was 31.35 m3. Screw pumps generated a flow of 9 cm/s, as measured with a Schiltknecht MiniWater 20, equipped with a MC20 MiniController and a MiniWater20 Micro sensor. Flow velocities of 18 and 1 cm/s were obtained by locally reducing channel width to 45 cm, and by inserting a secondary channel with an entrance grid, respectively (Figure S1-3). The channels contained a sediment layer of 15 cm thickness. Sediments in channel A and C were aged in the channel system over several years and can be regarded as largely natural in terms of carbon content and biology. Four segments contained macrophytes. Vegetation that had spread in other segments was removed before start of the experiments. Channel B was filled with fresh sediment and contained no plants.

Channel A was used for exposing passive samplers to elevated compound concentrations. Channel B was used for exposures in unspiked water. Channel C was used as a control channel. Peak concentration events were simulated by either transferring samplers back and forth between Channels A and B, or by partially draining Channel A followed by refilling with unspiked water.

Spiking of Channel A was done by injecting a spike solution (2 × 1 L) in methanol/water (10/90 v:v) at two opposite locations in the channel, using a valveless rotary piston pump (Ismatec, Switzerland) over a time period of 8 min. This resulted in a nearly homogeneous initial distribution of compounds over the length of the channel because spiked water travelled a distance of 43 m downstream of each spiking location during 8 min, as compared with 2 × 52 m for a complete water circulation. Spiked compounds were bentazon, carbamazepine, diclofenac, flufenacet, imidacloprid, metazachlor, nicosulfuron, pendimethalin, propiconazole, terbuthylazine, thiacloprid, and UV 326 (bumetrizole). Nominal concentrations (injected amount divided by channel volume) ranged between 50 and 1000 ng/L.

Sampler exposures

Exposure experiments were done during three time periods in 2019: 25 June–8 July (Phase 1), 27 August–13 September (Phase 2), and 14 October–1 November (Phase 3). Sampler types were silicone and SDB-RPS extraction disks with and without a PES membrane. Sampler types and exposure scenarios are summarized below and in Table 1.

Phase 1 was designed to allow a comparison of SDB-RPS and silicone sampler response to constant Cw and to peak concentration events. Exposure to constant Cw was only successful for nicosulfuron and terbuthylazine. The Cw of other compounds decreased with time, due to a variety of loss processes (e.g., sorption to channel components, biodegradation, photodegradation). Intentional peak concentration events were simulated by transferring samplers between the spiked Channel A to the unspiked Channel B. Samplers were mounted on racks that could be quickly transferred between the channels whenever needed (Figure S1-1).

Peak concentration events during Phase 2 were generated by spiking Channel A at t = 1, 3, and 6 days. Immediately after each spiking event, 15 m3 water was drained during 35 min, followed by refilling the channel at 31 L/min for 8 h. This draining/refilling operation was repeated two times more, with a refilling rate of 25 L/min for 10 h. Total time between start of the spiking and end of the draining/refilling procedure was 1.2 days. SDB-RPS and silicone samplers were deployed during Phase 2.

Exposure scenarios during Phase 3 included sampler exposure in Channel A only and transfer to Channel B after initial exposure in Channel A. SDB-RPS samplers with and without a PES membrane were deployed. The effect of fouling was evaluated by exposing SDB-RPS samplers that had developed a biofilm during 2 and 4 weeks in Channel C.

Daily averaged temperatures decreased by approximately 13 °C during each Phase. Average temperatures were 21 °C (Phase 1, 29 to 16 °C), 19 °C (Phase 2, 25 to 14 °C), and 11 °C (Phase 3, 16 to 1 °C).

Chemical analysis

Silicone sheets were extracted in a Soxhlet extractor (30 h, 30 mL methanol). SDB-RPS disks were shaken 3 times with 6 mL solvent (acetone, methanol, acetone) on a shaking Table (60 min, 120 RPM). Extracts were evaporated to 4 mL with a nitrogen stream. A 1 mL subsample was filtered through a PTFE filter (0.45 µm pore size), mixed with deuterated internal standards, and diluted with ultrapure water at a ratio of 1:10, 1:20, or 1:100 (v/v). Water samples were filtered and directly injected. For each batch of samples, an additional procedural blank was extracted to verify that the silicone sheets and SDB-RPS discs were not contaminated.

LC–MS/MS analysis was performed on an Agilent 1260 LC system coupled with a Sciex 4500 QTrap MS/MS System. Compounds were separated on a Phenomenex 100 × 2.1 mm Kinetex C18 column with 2.6 µm particle size and upstream security guard cartridge at 40 °C, using a gradient of 0.4 mM ammonium acetate in water and methanol. The tandem mass spectrometer was run in multiple reaction mode (MRM) with positive and negative electrospray ionization. Details of instrumental analysis, internal standards, detection limits, and quantitation limits are given in section S2.

Modelling of aqueous concentrations

Concentrations in water were modelled as an exponential decrease from the value immediately after spiking (Cw0) towards a constant level (Cw∞) as t→∞3 Cw=Cw0-Cw∞exp-αt+Cw∞

where α is a first-order decay constant and t is the time after spiking.

Accuracy of the Cw data was evaluated from the mass balance (MB) at the time of spiking.4 MB=Cw0Cw,pre-spike+Cw0,spike×100%

where (Cw,pre-spike) is the concentration immediately before spiking and Cw,0,spike is the increase in Cw, as calculated from spiked amount and channel volume.

Modelled Cw values (Eq. 3) were used for evaluating sampler-water exchange kinetics for all post-spike periods. Measured aqueous concentrations > LOQ were adopted for Channel B and for the pre-spike period in Channel A, but some inevitable uncertainty is associated with Cw values < LOQ, as Cw can have any value between 0 and LOQ in these cases. We chose for a conservative approach for dealing with this uncertainty. Concentrations below the LOQ were set to zero because the exposure water was aged, and isolated from direct and indirect sources. In some exceptional cases, there was evidence that Cw was > 0, as measurable amounts were detected in the samplers before the first spiking event. This was observed in Phase 2 for metazachlor, and in Phase 3 for metazachlor, nicosulfuron, and thiacloprid. After initial modelling of sampler-water exchange kinetics the pre-spike Cw of these compounds was adjusted to match the amounts in the samplers before spiking. For Channel B, metazachlor concentrations of 0.005 ng/mL were adopted because amounts in the samplers reached a plateau value in this channel during Phases 1 and 3. We verified that all hypothesized Cw values were indeed < LOQ.

Sampling rate model with time-variable Cw

The Rs model can be generalized to account for time-dependent Cw and nonzero initial concentrations (supplementary information, S3)5 Cst=exp-ketCs0+keKsw∫t0expket′Cwt′dt′

where Cs0 is the initial concentration in the sampler, and ke = Rs/(mKsw) is used to simplify notation. Evaluation of the integral is straightforward for the adopted Cw model (Eq. 3). The time domain was divided in separate stages (pre-spike and post-spike periods), each starting with its own t = 0. The final Cs from the previous stage was used as the Cs0 of the next stage. The most complex scenario was repeated sampler switching between Channels A and B during Phase 1, which comprised 7 stages (pre-spike exposure + three periods in Channel A + three periods in Channel B).

Combining the Cw model (Eq. 3) and the generalized Rs model (Eq. 5) yields the Rs model for exposures in Channel A.6 Cst=KswCw∞+Cs0-KswCw∞exp-ket+keKswCw0-Cw∞exp-αt-exp-ketke-α

Exposures at constant Cw were modeled by setting Cw0 = Cw∞. For nonlinear least squares analysis ke was replaced by Rs/(mKsw), because optimizing Ksw and ke would yield erroneous estimates of the standard errors of Rs.

Diffusion model with time-variable Cw

The diffusion model from Mikhailov and Özişik (1984) is more compactly written in terms of dimensionless time (τ)7 τ=DstL2

where L is the half-thickness for samplers that are exposed on both sides (silicone in the present study) or the thickness for samplers that are exposed on one-side, with an impermeable boundary at the other side (SDB-RPS in the present study). Model derivation for plane sheets is summarized in section S4. The time evolution of space-averaged concentrations is given by8 Csτ=∑n=1∞sin(μn)Fnμne-μn2τC~s,n0+μnsin(μn)Ksw′∫τ0Cw(τ′)eμn2τ′dτ′

where µn are solutions of9 cotμn-μnBi=0

The Biot number (Bi) is defined as10 Bi=kwLK′swDs

The term Fn in the denominator of Eq. 8 is a normalization factor.11 Fn=∫10cos2μnXdX=12+cosμnsin(μn)2μn

where X = x/L. The term C~s,n(0) between the square brackets in Eq. 8 is the integral transform of the initial concentration in the sorbent12 C~s,n0=∫10cos(μnX)Cs(X,0)dX

Substituting the Cw model (Eq. 3) into Eq. 8 yields13 Csτ=∑n=1∞sinμnFnμnC~s,n0e-μn2τ+μnsinμnKsw′Cw0-Cw∞-β+μn2e-βτ-e-μn2τ+Cw∞μn21-e-μn2τ

where β = αL2/Ds. Constant-Cw scenarios are covered by setting Cw0 = Cw∞.

The diffusion model has K’sw, Ds, and kw as adjustable parameters. For diffusion in porous media (SDB-RPS) this can be reduced to two parameters because Ds is inversely proportional to K’sw by (Booij 2021; Endo et al. 2019)14 Ds=φsDwθs2Ksw′

Dw at 25 °C was calculated from the McGowan molar volume (Schwarzenbach et al. 2016), and correction for exposure temperature was done using the relationship with dynamic viscosity (Dw ~ η−1.14) as suggested by Hayduk and Laudie (1974).

The time domain was divided into multiple stages that started at the time of sudden change in Cw (a new spike event or a transfer between Channels A and B), as outlined for the Rs model above. The initial concentration distribution for each next stage was set equal to the final concentration distribution for the previous stage.

Series resistance model

The series resistance model is useful for understanding the rate limiting steps and the effects of temperature, flow, and biofouling on the exchange kinetics. For single-phase polymers without biofouling layer this model takes the form (Booij 2021)15 1ko=1kw+δs(t)Ksw′Ds

where ko is the overall mass transfer coefficient (i.e., the proportionality constant between the flux and the effective concentration difference between sampler and water), kw is the mass transfer coefficient of the WBL, and Ds is the diffusion coefficient in the sorbent. The K’sw is the sorption coefficient in volume per volume units (concentration in bulk sorbent divided by concentration in water). The δs(t) is the effective thickness of the concentration boundary layer in the sorbent. The concentration gradient in the sorbent is infinite at t = 0, and hence δs(0) = 0. At increasing time, δs(t) gradually increases, and reaches a steady state value of approximately 0.3 to 0.4 times the half-thickness (L) when Dst/L2 > 0.1.

For porous sorbents, porous membranes, and biofouling, present knowledge suggests that the diffusion coefficient is inversely proportional to the sorption coefficient, which causes the product K’swDs to be independent of K’sw (Booij 2021; Booij et al. 2007)16 1ko=1kw+θm2dmφmDw+θs2δs(t)φsDw+θb2dbφbDw

where subscripts m and b refer to membrane and biofilm, ϕ is porosity, θ is tortuosity, and d is the thickness of membrane or biofilm. Here it is assumed that transport in membrane and biofilm is through the pore space only. The sorption coefficient K’sw still occurs implicitly in δs(t), because the concentration gradient in the sorbent is steeper for longer times when K’sw is larger.

Statistical methods

Best-fit parameters and their standard errors were obtained using weighted nonlinear least squares analysis in Microsoft Excel, using the Solver add-in (Billo 2001). A clear minimum in the residual sum of squares (SSQ) was not always observed. For the diffusion model this can occur when the kinetics is essentially controlled by the sorbent. A plot of SSQ vs. kw shows a plateau at large kw in that case, and only rises for kw values that are too small. The nature of the SSQ minimum (local minimum, global minimum, or plateau) was evaluated by running the model at fixed kw values between 0.1 and 1000 µm/s, optimizing the other parameters, and plotting SSQ vs. kw. When a plateau was found in the SSQ vs. kw plot, it was concluded that the simpler sorbent control model was applicable. When a minimum was found in the SSQ vs. kw plot, a partial F-test was used to determine if the complex model (adjustable kw) yielded a significantly better fit than the reduced model (fixed kw = 1000 µm/s).17 F=SSQreduced-SSQcomplexdfreduced-dfcomplexSSQcomplexdfcomplex

Values of kw were labelled as not significant when the F-value was smaller than the critical F-value (95% level of confidence, right-tailed F-test).

Similarly, for WBL-controlled kinetics, any Ds value that is larger than a certain threshold gives an equally good fit. Finally, when Cs is much smaller than the equilibrium concentration, then any high value of Ksw or K’sw will do. The significance of Ds and Ksw or K’sw was evaluated similarly to the significance of kw. The approach is further illustrated in section S5.

The partial F-test was also used to determine if statistically significant lag times occurred in the Rs model for SDB-RPS samplers that were covered by a PES membrane.

Results and discussion

Aqueous concentration modelling and compound selection

Modelling of Cw in Channel A was considered to be acceptable when a sufficient number of Cw values was above the LOQ, and when the mass balance at the time of spiking was between 80 and 120%. Bentazon, carbamazepine, imidacloprid, metazachlor, nicosulfuron (Phase 1 and 3 only), terbuthylazine (Phase 1 and 2 only), and thiacloprid, passed these criteria (Table 2). Mass balances for nicosulfuron during Phase 2 (120 to 160%) and for terbuthylazine in Phase 3 (57%) were outside the 80 to 120% window, but were still retained to allow comparison of sampler kinetics among the respective experimental Phases. Cw data for diclofenac, flufenacet, pendimethalin, propiconazole, and UV 326 did not meet the criteria for mass balance and detection frequency in any of the experimental Phases, and were excluded from further consideration (Tables S6-1, S6-2, and S6-3). The rapid disappearance of pendimethalin and UV 326 may be related to sorption to channel walls and sediment sections because of their relatively high logKow values of 5.20 and 5.55. The rapid disappearance of diclofenac may be related to degradation. The reason for low mass balances of flufenacet and propiconazole is unclear. Table 2 Parameters for aqueous concentration modelling of included compounds. (Model results for excluded compounds are given in Tables S6-1, S6-2, and S6-3)

	Phase 1	Phase 2	Phase 3	
Compound	Mass balance	α	n	Mass balance peak 1	Mass balance peak 2	Mass balance peak 3	α	n	Mass balance	α	n	
	(%)	(day−1)		(%)	(%)	(%)	(day−1)		(%)	(day−1)		
Bentazon	93	0.096	65	121	100	96	0.98	7	85	0.027	16	
Carbamazepine	82	0.017	65	96	100	98	2.3	32	97	0.56	16	
Imidacloprid	82	0.64	56	94	132	97	2.1	18	92	0.14	16	
Metazachlor	91	0.035	65	105	100	101	1.9	32	102	0.51	16	
Nicosulfuron	103	0	65	160	148	121	1.7	9	107	0	16	
Terbuthylazine	82	0	65	104	87	90	2.0	32	57	0	16	
Thiacloprid	80	0.028	65	94	95	92	1.9	29	87	0.043	16	

Silicone samplers

No accumulation in silicone was observed for bentazon, imidacloprid, and nicosulfuron. Accumulation of carbamazepine, metazachlor, terbuthylazine, and thiacloprid is characterized by low sorption coefficients (logKsw between 2.0 and 3.4) and large differences in Rs/A between compounds (0.1 to 33 L/[dm2 day])(Table 3). Amounts that are accumulated during a peak concentration event dissipate appreciable within a few days (Fig. 1, and Figure S7-1). TIW values are between 0.5 and 6 days, which makes silicone to be of little practical use for time-integrative sampling of these compounds, although this sampler is still interesting for research purposes. Table 3 Model results for silicone samplers. SE = standard error of the parameter in the column to the left. No accumulation was observed for bentazon, imidacloprid, and nicosulfuron. TIW is the time-integrative window (mKsw/Rs)

		Diffusion model		Rs model	
		logK'sw (mL/mL)	SE	Ds (µm2/s)	SE	kw (µm/s)	SE	s (%)	n	logKsw (mL/g)	SE	Rs/A (L/[dm2 day])	SE	s (%)	TIW (days)	
Carbamazepine	Phase 1	2.26	0.13	0.0056	0.0038	ns)1	-	16	24	2.04	0.04	0.32	0.05	40	3.7	
	Phase 2	2.28	0.01	0.0088	0.0010	0.4	0.2	13	25	2.13	0.03	0.36	0.04	19	4.0	
Metazachlor	Phase 1	2.88	0.02	0.15	0.04	7.2	3.1	21	24	2.87	0.02	15	1	24	0.5	
	Phase 2	3.10	0.06	0.012	0.004	ns)1	-	28	25	2.94	0.03	8.6	2.0	33	1.1	
Terbuthylazine	Phase 1	3.34	0.05	0.080	0.051	ns)1	-	18	24	3.30	0.02	33	5	20	0.6	
	Phase 2	3.48	0.02	0.032	0.007	ns)1	-	19	25	3.44	0.02	24	4	19	1.2	
Thiacloprid	Phase 1	3.10	0.18	0.000029	0.000020	ns)1	-	15	20	1.74	0.05	0.12	0.02	30	4.7	
	Phase 2	2.97	0.21	0.000044	0.000038	ns)1	-	20	25	1.76	0.06	0.10	0.01	26	6.0	
)1 Including kw as an adjustable parameter did not yield a significantly better fit

Fig. 1 Model fit and residual errors of thiacloprid sampling by silicone samplers in Phase 1 for the Rs model (left) and the diffusion model (right). Cw exposure scenarios were continuous exposure in Channel A (blue), exposure in Channel A for 2 days, followed by transfer to Channel B (amber), and three times switch between Channels A and B (green). Solid lines represent model fits

The diffusion model gave a better fit of the data than the Rs model (residual errors 13–28% vs. 19–40%, Table 3). No difference in residual errors was found for terbuthylazine (19%). Better performance of the diffusion model is supported by the more randomly distributed residual errors for this model (Fig. 1, Figure S7-1). No significant effect of flow (Phase 2) was observed for any of the compounds between 1 and 18 cm/s (Fig. 2, Figure S8-1).Fig. 2 Accumulation of thiacloprid by silicone samplers in Phase 2 for three flow velocities. Exposure in Channel A, sampling at t = 8 days

Including kw as an adjustable parameter in the diffusion model did not yield a better fit, except for carbamazepine in Phase 2 and metazachlor in Phase 1 (Table 3). This indicates insignificant rate control by the WBL. The kw estimate for carbamazepine in Phase 2 (0.4 ± 0.2 µm/s) is statistically significant (p = 0.03), but suspiciously small and likely erroneous. Significant flow effects would be expected for this compound if the kw estimate were realistic, but no flow effects were observed (Figure S8-1). Model results for fully sorbent-controlled kinetics yielded slightly higher residual errors (14% vs. 13%), and smaller Ds (0.0057 µm2/s) for carbamazepine. The kw estimate for metazachlor in Phase 1 (7.2 ± 3.1 µm/s) is significant (p = 0.02) and realistic, but no flow effects were observed for this compound for the Phase 2 data. We have no solid explanation for the statistical significance of kw for these two cases, but we consider that the relatively large standard errors of approximately 50% may be taken as a warning that the kw estimates are questionable after all. This further illustrates that conclusions on expected flow effects that are based on curve fitting need experimental verification by exposing samplers to different flow velocities.

SDB-RPS samplers without a membrane

Sorption coefficients for SDB-RPS samplers are substantially higher than for silicone (logKsw between 3.5 and 4.7; median values for Phases 1, 2, and 3, Table 4). Rs/A values range between 0.8 and 4.0 L/[dm2 day] (Rs between 0.1 and 0.5 L/day). The TIWs between 7 and 34 days indicate fair time-integrative capability for exposures of two weeks, except for bentazon (TIW = 10 days) and nicosulfuron (TIW = 7 days). Table 4 Model results for SDB-RPS samplers. SE = standard error of the parameter in the column to the left. The diffusion coefficient in the sorbent (Ds) follows from Eq. 14 and the optimized value of K’sw. TIW is the time-integrative window (mKsw/Rs)

Compound	Phase	Diffusion model		Rs model	
		logK'sw (mL/mL)	SE	kw (µm/s)	SE	s (%)	N	logKsw (mL/g)	SE	Rs/A (L/[dm2 day])	SE	s (%)	TIW (days)	
Bentazon	1	3.07	0.03	6.8	1.0	17	38	3.42	0.03	1.3	0.1	35	3.5	
	2	3.53	0.26	0.7	0.1	25	31	ns)1		0.38	0.02	27	 > 16	
	3	3.30	0.06	2.6	0.4	24	46	3.63	0.06	0.80	0.06	33	9.5	
Carbamazepine	1	5.01	0.07	8.9	0.6	18	38	4.37	0.05	5.2	0.3	23	8.0	
	2	5.08	0.08	6.6	0.6	13	36	4.43	0.05	4.0	0.2	15	12	
	3	5.48	0.16	3.4	0.2	20	53	4.63	0.09	2.6	0.09	20	30	
Imidacloprid	1	4.42	0.05	9.4	0.9	18	38	4.13	0.04	3.9	0.3	29	6.2	
	2	4.87	0.14	3.4	0.4	15	30	4.47	0.11	2.1	0.1	17	25	
	3	4.51	0.12	4.1	0.5	31	49	4.32	0.11	2.0	0.2	36	19	
Metazachlor	1	5.39	0.09	7.3	0.4	16	38	4.51	0.05	5.1	0.2	17	11	
	2	ns)1		3.9	0.2	14	36	ns)1		3.2	0.1	14	 > 16	
	3	5.23	0.15	3.8	0.3	23	53	4.45	0.08	2.7	0.1	23	19	
Nicosulfuron	1	3.01	0.05	1.9	0.2	21	38	3.30	0.04	0.75	0.06	32	4.7	
	2	3.17	0.13	1.9	0.2	21	38	ns)1	ns)1	0.45	0.04	32	 > 16	
	3	3.54	0.12	2.4	0.5	46	52	3.63	0.04	1.1	0.1	44	7.0	
Terbuthylazine	1	4.99	0.09	8.7	0.8	17	38	4.30	0.04	5.2	0.2	17	6.8	
	2	5.59	0.14	5.6	0.4	13	36	4.67	0.07	4.0	0.2	13	21	
	3	4.94	0.11	5.6	0.4	13	36	4.29	0.05	2.9	0.1	13	12	
Thiacloprid	1	5.20	0.07	7.9	0.5	16	38	4.45	0.05	5.1	0.2	19	9.9	
	2	5.83	0.23	4.5	0.3	13	36	4.84	0.15	3.4	0.1	14	37	
	3	5.50	0.18	3.5	0.2	22	53	4.69	0.12	2.6	0.1	22	34	
)1 Including K’sw or Ksw as an adjustable parameter did not yield a significantly better fit

The diffusion model yielded a better fit of the data for SDB-RPS samplers than the Rs model, but the difference in residual errors was typically only a few percent (Table 4). The largest difference was observed for bentazon (17 to 25% for the diffusion model versus 27 to 35% for the Rs model). The relatively large residual errors (44 to 46%) for nicosulfuron in Phase 3 result from appreciable scatter for t < 4 days (Figure S9-1). Residual errors were fairly randomly distributed, both for uptake in Channel A and for dissipation in Channel B (Fig. 3, Figure S9-1). This indicates that uptake and release kinetics are well-described by the same model parameter values, and that anisotropic exchange was not a major issue for the present data. Accumulated amounts with SDB-RPS are much higher than for silicone, due to the higher sorption coefficients of SDB-RPS. A similarity for these samplers is that accumulated compounds in SDB-RPS also dissipate on a time scale of a few days (Fig. 3, Figure S9-1).Fig. 3 Model plot (left) and residual errors (right) for thiacloprid sampling by SDB-RPS samplers in Phase 1, analyzed with the diffusion model. Results for the Rs model are nearly identical (Figure S7-1)

Amounts in SDB-RPS samplers increased with flow velocity for all exposure scenarios that were evaluated (Fig. 4, Figure S10-1). For fully WBL-controlled kinetics during the kinetic sampling stage at constant Cw it is expected that amounts follow a power law relationship with flow velocity: Ns ~ Un, with n = 1/2, because kw increases with the square root of U (Glanzmann et al. 2022; Stephens et al. 2005). Experimental values of n (obtained from the slopes of logNs versus logU) were 0.30 ± 0.15 (range 0 to 0.53, Table S10-1). The observation that slopes are generally smaller than 0.5 likely originates from partial equilibrium attainment and from the fact that the kinetics is only partially WBL-controlled.Fig. 4 Effect of flow on the accumulated amounts of thiacloprid in SDB-RPS samplers. Error bars represent standard deviations

The applicability of the square root relationship between kw and U was further evaluated as follows. First, the amounts at a flow velocity of 9 cm/s were calculated from the model parameters in Table 4. The amounts at U = 1 and 18 cm/s were then calculated by using a 3-times lower and √2 times higher kw, respectively. The ratio of calculated and measured amounts was 0.96 ± 0.18 (Fig. 5). The model slightly underestimated the amounts at 18 cm/s (amount ratio 0.87 ± 0.13), but overall gives a fair prediction of flow effects for the present velocity range.Fig. 5 Modelling of flow effects on the accumulation in SDB-RPS samplers. Vertical axis: ratio of modelled and experimental amounts at flow velocities of 1, 9, and 18 cm/s, adopting the relationship kw ~ U1/2

SDB-RPS samplers with a PES membrane

In SDB-RPS samplers with a PES membrane appreciable sorption to PES occurred for thiacloprid and terbuthylazine (Fig. 6, Figure S11-1). Sorption of bentazon, carbamazepine, imidacloprid, and metazachlor to PES was less than 20% at the end of the exposure. Nicosulfuron showed erratic responses (residual errors of 65%, Figure S11-1), and is excluded from further discussion of membrane effects. Amounts in samplers with a membrane were 2 to 3 times smaller than the amounts in samplers without a membrane, in line with observations by Shaw et al. (2009). An exception was bentazon, for which the amounts in both sampler types were approximately the same, which indicates again that the exchange kinetics for this compound is controlled by the sorbent.Fig. 6 Accumulation of thiacloprid in SDB-RPS (filled circles) and PES membrane (open circles), for uptake in Channel A (blue) and dissipation in Channel B (amber)

The presence of a PES membrane appreciably improved the time-integrative capabilities of SDB-RPS samplers. TIWs were between 17 and 140 days (Table 5). Rs/A values spanned a narrower range (0.6 to 0.9 L/[dm2 day]) than for samplers without a membrane (0.8 to 4.0 L/[dm2 day]). Minor lag times were observed for carbamazepine (0.1 days), terbuthylazine (0.4 days), thiacloprid (0.2 days). For the present set of compounds, the disadvantages of analyte sorption to PES are outweighed by the advantages of higher TIWs and more consistent Rs values. Table 5 Rs model results including lag times for SDB-RPS samplers with a PES membrane. SE = standard error of the parameter in the column to the left. TIW is the time-integrative window (mKsw/Rs)

	logKsw (mL/g)	SE	Lag time (days)	SE	Rs/A (L/[dm2 day])	SE	n	s (%)	TIW (days)	
Bentazon	3.78	0.05	ns)1		0.63	0.03	18	14	17	
Carbamazepine	ns)1		0.06	0.03	0.87	0.03	18	11	55)2	
Imidacloprid	ns)1		ns)1		0.71	0.03	16	15	53)2	
Metazachlor	ns)1		ns)1		0.77	0.02	20	10	71)2	
Terbuthylazine	4.05	0.08	0.38	0.06	0.91	0.06	15	8	22	
Thiacloprid	ns)1	ns	0.21	0.02	0.62	0.03	17	16	140)2	
)1 Including Ksw or lag time as an adjustable parameter did not yield a significantly better fit

)2 Based on the median logKsw from data for SDB-RPS disks without a membrane (Table 4)

The diffusion model yielded the same residual errors as the Rs model for all compounds except bentazon (Table S11-1). Residual errors were approximately 13% for carbamazepine, imidacloprid, and metazachlor, and were much larger for terbuthylazine (26%) and thiacloprid (38%). For the latter two compounds modeled amounts at short time scales were 80% higher than measured, and up to 40% lower at long time scales (Figure S11-1). This error structure is in line with the occurrence of lag times that originate from sorption to PES. This has been observed for sampling by Chemcatchers with SDB-RPS sorbent and PES membranes, for compounds with logKPES-w larger than 2 (Vermeirssen et al. 2012) or 3.5 (Estoppey et al. 2019). Data from Vermeirssen et al. (2012) indicate lag times for terbuthylazine of 0.5 to 4 days for Chemcatchers and 0.2 to 2 days for POCIS with HLB sorbent.

Lag times (tlag) were evaluated using the Rs model, by replacing t with t-tlag in Eq. 6, and ignoring amounts for t < tlag. Including lag times in the model yields a better fit only for carbamazepine, terbuthylazine, and thiacloprid (Table 5, Figure S11-2). Lag times < 0.4 days are minor when compared with typical exposure times of 2 weeks. Residual errors are slightly smaller for samplers with a PES membrane than for bare SDB-RPS samplers (12 vs 20%). Sampling rates are similar for all compounds (Rs = 95 ± 15 mL/day; mean ± standard deviation). This value is fairly well in line with a series resistance model for mixed WBL/membrane controlled kinetics (Eq. 16). Adopting φ/θ 2 = 0.78 for the PES membrane, Dw = 350 µm2/s at 11 °C during Phase 3, a membrane thickness of 150 µm, and assuming transport through the pore space only, results in a membrane conductivity of 0.78 × 350/150 = 1.8 µm/s (Eq. 16). Adopting further a WBL conductivity (kw) for Phase 3 of 3.6 µm/s (Table 4), yields a combined conductivity of (1/1.8 + 1/3.6)−1 = 1.2 µm/s = 1.0 L/(dm2 day). For a surface area of 0.126 dm2 the expected Rs would be 130 mL/day, which is 1.4 times higher than the experimental Rs. Further understanding of the experimental Rs can be gained when experimental values of porosity and tortuosity for the present membrane type are available.

kw estimates

Experimental kw values for SDB-RPS samplers were similar to the values that were calculated from water flow velocities and water properties, using the model from Glanzmann et al. (2022). (Fig. 7). The range of observed/predicted kw ratios was 0.4 to 1.7, excluding the likely erroneous value for bentazon in Phase 2. Some uncertainty in the observed kw values is associated with the 13 °C temperature decrease during all Phases, but kw is expected to decrease by only a factor of 1.3 over this temperature range (Section S5 from Booij and Chen 2018). Although the average observed kw agrees well with the Glanzmann et al. prediction, present data give no reason to prefer their recommended proportionality constant of 0.52 over the theoretical value of 0.664. Experimentally determined kw values from alabaster dissolution or dissipation of performance reference compounds from silicone samplers would be a valuable addition for future experiments (Booij et al. 2017; Glanzmann et al. 2022).Fig. 7 Estimates of kw from the diffusion model versus predicted values from Glanzmann et al. (2022). The solid line represents the 1:1 relationship. Error bars represent standard errors

Biofouling effects

Allowing biofouling to develop on SDB-RPS samplers in Channel C, prior to exposure in Channel A, caused an appreciable reduction of accumulated amounts (Fig. 8, Figure S12-1). This reduction was 13 ± 6% per week of pre-fouling (relative to the amount for samplers without pre-fouling, Fig. 9). The biofouling effect for terbuthylazine is less pronounced than for other compounds, because of the non-zero background concentrations in Channel C for this compound.Fig. 8 Effect of pre-fouling on accumulation of thiacloprid by SDB-RPS samplers in Channel A at t = 2 and 4 days, after pre-fouling for 0, 2, 4 weeks

Fig. 9 Reduction of accumulation by SDB-RPS samplers in Channel A at t = 2 and 4 days, after pre-fouling in Channel C

The relatively compound-independent effect of biofouling is in line with an earlier suggestion that the fouling layer can be modelled as an immobilized water layer (Eq. 16). A crude estimate of the effective thickness of this layer can be obtained as follows. The overall mass transfer coefficient for WBL and sorbent is ko = Rs/A, which is for Phase 3 equal to 2.1 ± 0.8 L/(dm2 day) = 2.4 ± 1.0 µm/s (based on Rs/A values for all compounds, Table 4). Four weeks of biofouling caused a reduction of approximately 50%, indicating that the biofouling resistance (θ 2db/[φ b Dw]) is equal to the combined resistances of WBL and sorbent. Adopting Dw = 350 µm2/s for 11 °C during Phase 3 yields an effective thickness of the fouling layer of θ b2db/φ b = 350/2.4 = 150 µm, which looks like a reasonable estimate in view of the low degree of visible fouling (Figure S12-2).

Implications

The analytical sampling rate model was adequate for evaluating passive sampler responses to time-variable Cw of the studied compounds. The best performance was observed for SDB-RPS disks with a PES membrane, despite the occurrence of analyte sorption to PES. Sampling rates spanned a relatively narrow range for all analytes, lag times were minor, and TIWs were highest. This conclusion is provisional because this sampler type was only evaluated in Phase 3 of the present study. Additional study with sampler exposures to multiple peak concentration events is needed for this sampler design. The time-integrative capability of SDB-RPS samplers without a PES membrane is fair, but Rs is more strongly compound dependent and harder to understand. Silicone samplers are of little use for the time-integrative sampling of the present compounds because of their short TIWs.

The diffusion model for exposures to time-variable Cw may be useful when a more fundamental understanding of sampler-water exchange mechanisms is needed. Knowledge of these mechanisms allows to assess the effect of flow and temperature on the exchange kinetics, using models for the temperature and flow dependency of kw, and the temperature dependency of K’sw and Ds. In the present study, the model parameters (kw and K’sw) were obtained by curve fitting of uptake and dissipation data. It would be interesting to test if calibration cost can be reduced by obtaining K’sw from batch sorption experiments and calculating kw from existing models.

The diffusion model is of limited use for practical application in routine monitoring. The better accuracy of this model only takes effect when exposure concentrations are either constant, or when time and duration of peak concentration events are known. This is usually not the case, but the diffusion model may be used to quantify the uncertainties that are associated with intermittent ambient concentrations, for example by comparing scenarios where peak concentrations occur at the beginning or at the end of the sampler exposure. The diffusion model may similarly be used to calculate the distribution of sampling rates that may be encountered, given the uncertainties in the occurrence of peak events and uncertainties in flow velocity and temperature at the exposure sites.

More importantly, the diffusion model may be used to identify sampler designs for which the transport resistance of the sorbent is much smaller than the combined resistance of membrane and WBL. In that case the numerical accuracy of the Rs model is the same as for the diffusion model, which greatly reduces the computational burden.

Supplementary Information

Below is the link to the electronic supplementary material.Supplementary file1 (DOCX 5437 KB)

Acknowledgements

We gratefully acknowledge Ronny Schmiediche and Stefan Loth for their experimental support at the stream channels. We thank Anja Duffek, Bilgin Karaoglan, Daniel Laska, and Etienne Vermeirssen for the input on the experimental design. We thank Etienne Vermeirssen for scientific support in the method development for the Chemcatcher. K. Booij has received financial support from the Federal Institute of Hydrology (BfG).

Author contribution

Study conceptualization and methodology were done by BB, DS, CK, SM, and FK. Experiments were performed by BB, DS, CK, SM, FK, and CM. Sample preparations and analysis were conducted by JB, CM, DS, BB, and CK. Stream channels were provided by the German Environment Agency (SM and FK). The manuscript was written by KB, BB, and CK. Critical revision of the manuscript was done by all authors. Data processing was performed by KB, CK, and BB. Data modelling was done by KB. Supervision and administration were conducted by SM and BB.

All authors read and approved the final manuscript.

Funding

This study was financed by the German Environment Agency (FKZ: 3717 67 418 0).

Data availability

The datasets used and/or analyzed during the current study are available from the corresponding author on request.

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

Becker B, Claus E, Möhlenkamp C, Bachtin J, Paschke A, Schäfer S (2020) Anwendung von Passivsammlern in Überwachungsprogrammen gemäß WRRL und MSRL - Identifizierung von Kontaminationsschwerpunkten, Referenzstandorten und neuen Schadstoffen (AnPassa). (No. UBA Texte 110/2020). Umweltbundesamt, Dessau-Roßlau. https://www.umweltbundesamt.de/publikationen/anwendung-von-passivsammlern-in
Bernard M Boutry S Tapie N Budzinski H Mazzella N Lab-scale investigation of the ability of polar organic chemical integrative sampler to catch short pesticide contamination peaks Environ Sci Pollut Res 2018 10.1007/s11356-018-3391-2
Bernard M, Boutry S, Tapie N, Budzinski H, Mazzella N (2018) Lab-scale investigation of the ability of polar organic chemical integrative sampler to catch short pesticide contamination peaks. Environ Sci Pollut Res. 10.1007/s11356-018-3391-210.1007/s11356-018-3391-2
Billo EJ Billo EJ Non-linear regression using the solver Excel for chemists: a comprehensive guide 2001 New York John Wiley & Sons Inc 223 238
Billo EJ (2001) Non-linear regression using the solver. In: Billo EJ (ed) Excel for chemists: a comprehensive guide. John Wiley & Sons Inc, New York, pp 223–238. 10.1002/0471220582.ch12
Booij K Passive sampler exchange kinetics in large and small water volumes under mixed rate control by sorbent and water boundary layer Environ Toxicol Chem 2021 40 1241 1254 10.1002/etc.4989 33492717
Booij K (2021) Passive sampler exchange kinetics in large and small water volumes under mixed rate control by sorbent and water boundary layer. Environ Toxicol Chem 40:1241–1254. 10.1002/etc.498933492717 10.1002/etc.4989
Booij K Chen S Review of atrazine sampling by polar organic chemical integrative samplers and Chemcatcher Environ Toxicol Chem 2018 37 1786 1798 10.1002/etc.4160 29687480
Booij K, Chen S (2018) Review of atrazine sampling by polar organic chemical integrative samplers and Chemcatcher. Environ Toxicol Chem 37:1786–1798. 10.1002/etc.416029687480 10.1002/etc.4160
Booij K Hofmans HE Fischer CV van Weerlee EM Temperature-dependent uptake rates of non-polar organic compounds by semipermeable membrane devices and low-density polyethylene membranes Environ Sci Technol 2003 37 361 366 10.1021/es025739i 12564909
Booij K, Hofmans HE, Fischer CV, van Weerlee EM (2003) Temperature-dependent uptake rates of non-polar organic compounds by semipermeable membrane devices and low-density polyethylene membranes. Environ Sci Technol 37:361–366. 10.1021/es025739i12564909 10.1021/es025739i
Booij K Vrana B Huckins JN Greenwood R Mills GA Vrana B Chapter 7, Theory, modelling and calibration of passive samplers used in water monitoring Passive Sampling Techniques in Environmental Monitoring 2007 Amsterdam Elsevier 141 169
Booij K, Vrana B, Huckins JN (2007) Chapter 7, Theory, modelling and calibration of passive samplers used in water monitoring. In: Greenwood R, Mills GA, Vrana B (eds) Passive Sampling Techniques in Environmental Monitoring. Elsevier, Amsterdam, pp 141–169
Booij K Maarsen NL Theeuwen M van Bommel R A method to account for the effect of hydrodynamics on polar organic compound uptake by passive samplers Environ Toxicol Chem 2017 36 1517 1524 10.1002/etc.3700 27925279
Booij K, Maarsen NL, Theeuwen M, van Bommel R (2017) A method to account for the effect of hydrodynamics on polar organic compound uptake by passive samplers. Environ Toxicol Chem 36:1517–1524. 10.1002/etc.370027925279 10.1002/etc.3700
Charriau A Lissalde S Poulier G Mazzella N Buzier R Guibaud G Overview of the Chemcatcher® for the passive sampling of various pollutants in aquatic environments Part A: principles, calibration, preparation and analysis of the sampler Talanta 2016 148 556 571 10.1016/j.talanta.2015.06.064 26653485
Charriau A, Lissalde S, Poulier G, Mazzella N, Buzier R, Guibaud G (2016) Overview of the Chemcatcher® for the passive sampling of various pollutants in aquatic environments Part A: principles, calibration, preparation and analysis of the sampler. Talanta 148:556–571. 10.1016/j.talanta.2015.06.06426653485 10.1016/j.talanta.2015.06.064
Endo S Matsuura Y Characterizing sorption and permeation properties of membrane filters used for aquatic integrative passive samplers Environ Sci Technol 2018 52 2118 2125 10.1021/acs.est.7b05144 29366322
Endo S, Matsuura Y (2018) Characterizing sorption and permeation properties of membrane filters used for aquatic integrative passive samplers. Environ Sci Technol 52:2118–2125. 10.1021/acs.est.7b0514429366322 10.1021/acs.est.7b05144
Endo S Matsuura Y Vermeirssen ELM Mechanistic model describing the uptake of chemicals by aquatic integrative samplers: comparison to data and implications for improved sampler configurations Environ Sci Technol 2019 53 1482 1489 10.1021/acs.est.8b06225 30608657
Endo S, Matsuura Y, Vermeirssen ELM (2019) Mechanistic model describing the uptake of chemicals by aquatic integrative samplers: comparison to data and implications for improved sampler configurations. Environ Sci Technol 53:1482–1489. 10.1021/acs.est.8b0622530608657 10.1021/acs.est.8b06225
Estoppey N Mathieu J Gascon Diez E Sapin E Delémont O Esseiva P de Alencastro LF Coudret S Folly P Monitoring of explosive residues in lake-bottom water using Polar Organic Chemical Integrative Sampler (POCIS) and chemcatcher: determination of transfer kinetics through Polyethersulfone (PES) membrane is crucial Environ Pollut 2019 252 767 776 10.1016/j.envpol.2019.04.087 31200202
Estoppey N, Mathieu J, Gascon Diez E, Sapin E, Delémont O, Esseiva P, de Alencastro LF, Coudret S, Folly P (2019) Monitoring of explosive residues in lake-bottom water using Polar Organic Chemical Integrative Sampler (POCIS) and chemcatcher: determination of transfer kinetics through Polyethersulfone (PES) membrane is crucial. Environ Pollut 252:767–776. 10.1016/j.envpol.2019.04.08731200202 10.1016/j.envpol.2019.04.087
Gallé T Frelat M Huck V Bayerle M Pittois D Braun C Quantitative use of passive sampling data to derive a complete seasonal sequence of flood event loads: a case study for maize herbicides in Luxembourg Environ Sci Process Impacts 2020 22 294 304 10.1039/C9EM00487D 31939971
Gallé T, Frelat M, Huck V, Bayerle M, Pittois D, Braun C (2020) Quantitative use of passive sampling data to derive a complete seasonal sequence of flood event loads: a case study for maize herbicides in Luxembourg. Environ Sci Process Impacts 22:294–304. 10.1039/C9EM00487D31939971 10.1039/C9EM00487D
Glanzmann V Booij K Reymond N Weyermann C Estoppey N Determining the mass transfer coefficient of the water boundary layer at the surface of aquatic integrative passive samplers Environ Sci Technol 2022 10.1021/acs.est.1c08088
Glanzmann V, Booij K, Reymond N, Weyermann C, Estoppey N (2022) Determining the mass transfer coefficient of the water boundary layer at the surface of aquatic integrative passive samplers. Environ Sci Technol. 10.1021/acs.est.1c0808810.1021/acs.est.1c08088
Harman C Reid M Thomas KV In situ calibration of a passive sampling device for selected illicit drugs and their metabolites in wastewater, and subsequent year-long assessment of community drug usage Environ Sci Technol 2011 45 5676 5682 10.1021/es201124j 21648435
Harman C, Reid M, Thomas KV (2011) In situ calibration of a passive sampling device for selected illicit drugs and their metabolites in wastewater, and subsequent year-long assessment of community drug usage. Environ Sci Technol 45:5676–5682. 10.1021/es201124j21648435 10.1021/es201124j
Harman C Allan IJ Vermeirssen ELM Calibration and use of the polar organic chemical integrative sampler - a critical review Environ Toxicol Chem 2012 31 2724 2738 10.1002/etc.2011 23012256
Harman C, Allan IJ, Vermeirssen ELM (2012) Calibration and use of the polar organic chemical integrative sampler - a critical review. Environ Toxicol Chem 31:2724–2738. 10.1002/etc.201123012256 10.1002/etc.2011
Hayduk W Laudie H Prediction of diffusion coefficients for nonelectrolytes in dilute aqueous solutions AIChE J 1974 20 611 615 10.1002/aic.690200329
Hayduk W, Laudie H (1974) Prediction of diffusion coefficients for nonelectrolytes in dilute aqueous solutions. AIChE J 20:611–615. 10.1002/aic.69020032910.1002/aic.690200329
Huckins JN Petty JD Booij K Monitors of organic chemicals in the environment: semipermeable membrane devices 2006 New York Springer
Huckins JN, Petty JD, Booij K (2006) Monitors of organic chemicals in the environment: semipermeable membrane devices. Springer, New York
Mazzella N Debenest T Delmas F Comparison between the polar organic chemical integrative sampler and the solid-phase extraction for estimating herbicide time-weighted average concentrations during a microcosm experiment Chemosphere 2008 73 545 550 10.1016/j.chemosphere.2008.06.009 18649919
Mazzella N, Debenest T, Delmas F (2008) Comparison between the polar organic chemical integrative sampler and the solid-phase extraction for estimating herbicide time-weighted average concentrations during a microcosm experiment. Chemosphere 73:545–550. 10.1016/j.chemosphere.2008.06.00918649919 10.1016/j.chemosphere.2008.06.009
Mikhailov MD Özişik MN Unified analysis and solutions of heat and mass diffusion 1984 New York John Wiley & Sons
Mikhailov MD, Özişik MN (1984) Unified analysis and solutions of heat and mass diffusion. John Wiley & Sons, New York
Morosini C, Marsoni M, Torretta V, Conti F, Ragazzi M, Rada EC, Cioca G (2017) Factors affecting spatial and temporal concentration variability of pharmaceuticals: comparison between two WWTPs. Sustainability 9. 10.3390/su9081466
Noro K Yabuki Y Banno A Tawa Y Nakamura S Validation of the application of a Polar Organic Chemical Integrative Sampler (POCIS) in non-steady-state conditions in aquatic environments J Water Environ Technol 2019 17 432 447 10.2965/jwet.19-057
Noro K, Yabuki Y, Banno A, Tawa Y, Nakamura S (2019) Validation of the application of a Polar Organic Chemical Integrative Sampler (POCIS) in non-steady-state conditions in aquatic environments. J Water Environ Technol 17:432–447. 10.2965/jwet.19-05710.2965/jwet.19-057
Ort C Nuijs ALN Berset J-D Bijlsma L Castiglioni S Covaci A Voogt P Emke E Fatta-Kassinos D Griffiths P Hernández F González-Mariño I Grabic R Kasprzyk-Hordern B Mastroianni N Meierjohann A Nefau T Östman M Pico Y Racamonde I Reid M Slobodnik J Terzic S Thomaidis N Thomas KV Spatial differences and temporal changes in illicit drug use in Europe quantified by wastewater analysis Addiction 2014 109 1338 1352 10.1111/add.12570 24861844
Ort C, Nuijs ALN, Berset J-D, Bijlsma L, Castiglioni S, Covaci A, Voogt P, Emke E, Fatta-Kassinos D, Griffiths P, Hernández F, González-Mariño I, Grabic R, Kasprzyk-Hordern B, Mastroianni N, Meierjohann A, Nefau T, Östman M, Pico Y, Racamonde I, Reid M, Slobodnik J, Terzic S, Thomaidis N, Thomas KV (2014) Spatial differences and temporal changes in illicit drug use in Europe quantified by wastewater analysis. Addiction 109:1338–1352. 10.1111/add.1257024861844 10.1111/add.12570
Schreiner VC Bakanov N Kattwinkel M Könemann S Kunz S Vermeirssen ELM Schäfer RB Sampling rates for passive samplers exposed to a field-relevant peak of 42 organic pesticides Sci Total Environ 2020 740 140376 10.1016/j.scitotenv.2020.140376 32927560
Schreiner VC, Bakanov N, Kattwinkel M, Könemann S, Kunz S, Vermeirssen ELM, Schäfer RB (2020) Sampling rates for passive samplers exposed to a field-relevant peak of 42 organic pesticides. Sci Total Environ 740:140376. 10.1016/j.scitotenv.2020.14037632927560 10.1016/j.scitotenv.2020.140376
Schwarzenbach RP Gschwend PM Imboden DM Environmental organic chemistry 2016 3 Hoboken John Wiley & Sons
Schwarzenbach RP, Gschwend PM, Imboden DM (2016) Environmental organic chemistry, 3rd edn. John Wiley & Sons, Hoboken
Shaw M Mueller JF Time integrative passive sampling: how well do Chemcatchers integrate fluctuating pollutant concentrations? Environ Sci Technol 2009 43 1443 1448 10.1021/es8021446 19350917
Shaw M, Mueller JF (2009) Time integrative passive sampling: how well do Chemcatchers integrate fluctuating pollutant concentrations? Environ Sci Technol 43:1443–1448. 10.1021/es802144619350917 10.1021/es8021446
Shaw M Eaglesham G Mueller JF Uptake and release of polar compounds in SDB-RPS Empore Disks; implications for their use as passive samplers Chemosphere 2009 75 1 7 10.1016/j.chemosphere.2008.11.072 19201009
Shaw M, Eaglesham G, Mueller JF (2009) Uptake and release of polar compounds in SDB-RPS Empore Disks; implications for their use as passive samplers. Chemosphere 75:1–7. 10.1016/j.chemosphere.2008.11.07219201009 10.1016/j.chemosphere.2008.11.072
Stephens BS Kapernick A Eaglesham G Mueller J Aquatic passive sampling of herbicides on naked particle loaded membranes: accelerated measurement and empirical estimation of kinetic parameters Environ Sci Technol 2005 39 8891 8897 10.1021/es050463a 16323791
Stephens BS, Kapernick A, Eaglesham G, Mueller J (2005) Aquatic passive sampling of herbicides on naked particle loaded membranes: accelerated measurement and empirical estimation of kinetic parameters. Environ Sci Technol 39:8891–8897. 10.1021/es050463a16323791 10.1021/es050463a
Tcaciuc AP Apell JN Gschwend PM Modeling the transport of organic chemicals between polyethylene passive samplers and water in finite and infinite bath conditions Environ Toxicol Chem 2015 34 2739 2749 10.1002/etc.3128 26109238
Tcaciuc AP, Apell JN, Gschwend PM (2015) Modeling the transport of organic chemicals between polyethylene passive samplers and water in finite and infinite bath conditions. Environ Toxicol Chem 34:2739–2749. 10.1002/etc.312826109238 10.1002/etc.3128
Thompson JM Hsieh CH Luthy RG Modeling uptake of hydrophobic organic contaminants into polyethylene passive samplers Environ Sci Technol 2015 49 2270 2277 10.1021/es504442s 25607420
Thompson JM, Hsieh CH, Luthy RG (2015) Modeling uptake of hydrophobic organic contaminants into polyethylene passive samplers. Environ Sci Technol 49:2270–2277. 10.1021/es504442s25607420 10.1021/es504442s
Vermeirssen ELM Dietschweiler C Escher BI van der Voet J Hollender J Transfer kinetics of polar organic compounds over polyethersulfone membranes in the passive samplers POCIS and Chemcatcher Environ Sci Technol 2012 46 6759 6766 10.1021/es3007854 22594693
Vermeirssen ELM, Dietschweiler C, Escher BI, van der Voet J, Hollender J (2012) Transfer kinetics of polar organic compounds over polyethersulfone membranes in the passive samplers POCIS and Chemcatcher. Environ Sci Technol 46:6759–6766. 10.1021/es300785422594693 10.1021/es3007854
