Transport and retention of micro-polystyrene in coarse riverbed sediments: effects of flow velocity, particle and sediment sizes

Riverbed sediments have recently been found to be an important reservoir for microplastics. But the hydrogeological factors that control the abundance of microplastics are complex and conceptual frameworks priorising the parameters affecting their transport and retention during deep riverbed filtration are still missing. In this study a series of saturated column experiments was conducted to investigate the vertical distribution patterns of secondary polystyrene fragments (100–2000 μm) in dependence on their particle size, grain size of the sediment, seepage velocity and duration of infiltration flow. The columns with a length of 50 cm were operated with flow velocities between 1.8 m d −1 and 27 m d −1 . Invasive samples obtained after the experiments were density separated and then depth profiles of microplastic concentrations were retrieved using fluorescence imaging analysis. Most polystyrene particles were retained in the upper 20 cm and 15 cm of the medium gravel and coarse sand sediments, respectively. Through the high particle retention riverbed sediments can act as a temporary sink or long term retention site for the transport of microplastic particles (MPPs) from streams to oceans. A small fraction of particles ranging from 100 to 500 μm in size was observed down to infiltration depths of 50 cm suggesting that MPPs at the pore scale have the potential to be advectively transferred via hyporheic exchange or induced bank filtration into coarse riverbed sediments and alluvial aquifers. MPP abundance over column depth follows an exponential relationship with a filter coefficient that was found to depend significantly on the flow rate, MPP and sediment grain size, as indicated by multiple linear regression (R 2 = 0.92). The experimentally derived empirical relation allows to estimate particle abundances of initially negatively buoyant MPP in riverbed sediments by surface water infiltration.

Rivers play an important role in the transport of plastic into lakes, seas, and oceans [10].Microplastics concentrations in rivers vary over a wide range, and may exceed several thousands of particles per cubic meter of river water in urbanized areas [17,38,39].The particle density, shape, and size as well as as physicochemical or biotic interactions determine whether MPP can be buoyant, neutral or sink in surface waters [19,37,46].For polystyrene particles (100-1000 μm), for instance, the modelled mean residence time in the epilimnion (2-18 m) of a stratified lake was estimated to be less than 1 day.In contrast, for the smallest microplastic particles (<= 30 μm) the mean residence times were orders of magnitudes higher [12].As a result, especially the smallest fraction of MPP tend to float in the water column and may never reach bed sediments.In stream ecosystems hyporheic exchange is considered as an important mechanism for the exchange of water and energy between water column and bed sediments [27] and for the delivery of microplastic particles to the streambed [11].
A certain number of MPP sink onto, infiltrate and accumulate in riverbed sediments [2,29].Recent studies examined the spatial and temporal distribution of microplastics in top riverine sediments worldwide, often with highest abundance of PE, PP, and PS [18].In sediment samples from the Thames River (Canada), Corcoran et al. [7] found the highest number of MPPs in samples with the finest grain size and with the highest portion of organic debris.Apart from that, MPPs can be remobilized into the surface water, transported and reemerged into the water column via hyporheic flow paths [11] or further relocated into the aquifer.In bank filtrate of an sandy to gravelly, alluvial aquifer (adjacent the Lower Havel River, North-East Germany) and in raw drinking waters exclusively from groundwater (Oldenburg, North-West Germany) only low concentrations of 4 MPP m −3 (100-700 μm) [25] and of 0.7 MPP m −3 (50-150 μm in size) [24] were detected, respectively.But studies analysing even smaller particles found much higher concentrations of 197 MPP m −3 (MPPs > 5 μm, South-West Germany) [30] and of 3.6 * 10 6 MPP m −3 (MPPs > 1 μm, Czech Republic) [31] in raw drinking water.Detected particles were mostly identified as PE, PA, PET, PVC, PS or epoxy resin.
Microplastic abundances in riverbed sediments are difficult to detect and so far only a few studies are available.Microplastic abundances can exhibit a clear decreasing trend with increasing depth.Fraser et al. [14] highlighted that the mean abundance of microplastic (50-5000 μm) in Tonglu and Fuyang riverbed sediments (China), that were collected with a gravity corer, significantly decreased from 290 MPP kg −1 (surface sediments), via 130 MPP kg −1 (5-10 cm) to 60 MPP kg −1 (10-15 cm).PE, PS, and PET fragments were the dominant polymer types present in sediment samples and accounted for 86-98% of the total abundance.And 63% of MPPs at 10-15 cm were larger than 100 μm.
Frei et al. [15] found highest MPP accumulation for particle sizes from 20 to 500 μm in the upper 10 cm (exceeding 50,000 MPP kg −1 ) of five freeze cores extracted from the Roter Main River sediments (Germany), consisting of high fractions of sand, medium to fine gravels and cobbles.MPP concentrations generally decreased with increasing depth (60 cm) except for the 10-20 cm segment where the lowest particle abundance (≈ 10,000 MPP kg −1 ) was detected.Large PS particles (500-5000 μm) were found down to a depth of 60 cm below the streambed surface.
However, also other types of microplastic distribution in riverbeds can occur, for example, microplastic concentrations were found to increase from shallow to deep layers in Qinhuai River sediments (China) mainly composed by silt [28].The mean MPP concentration across 15 sites gradually increased with increasing depth from 20 MPP kg −1 (0-10 cm) to 47 MPP kg −1 (41-50 cm) when analysing 10 cm thick segments.The average number of MPP < 2 mm in each layer across 15 sites increased gradually with depth, whereas the number of largest microplastics (> = 2 mm) decreased with the increasing depth [28].Thus, the riverbed could represent an important reservoir for microplastics and may act as a long-term pollutant sink, provided that the previously deposited microplastics are not fully remobilised or substantially carried into aquifers.
Soil column or lysimeter experiments can be very helpful in order to better understand distribution patterns of MPPs as observed in the environment and to obtain a better grasp on their transport processes in unsaturated and saturated sediments.Column experiments are commonly conducted with well-characterized granular media (e.g., glass beads, gravel or sand) to explore the physical and chemical influences on the transport and retention of microplastic particles in porous media via monitoring changes in the column effluent microplastic concentration as a function of time or detecting microplastic retention profiles.
For saturated conditions as in riverbed sediments, studies so far focused on transport of micron sized plastic PS spheres (< 10 μm), thus a primary MPP.For them, particle mobility decreased and deposition increased for increasing ionic strength, lower flow rates, and larger ratio of the colloid to the median grain diameter [3,4,9].Furthermore, the aggregation kinetics of MPP is an important factor in controlling their environmental fate and transport, whereby the aggregation of PS particles among themselves in freshwater environments (i.e.lakes and rivers) is unlikely [8,21].Dong et al. [9] found that 2 μm sized PS were not retained in sand columns operated with deionized water.The ionic strength and diameter of the collector and MPP were the dominant factors affecting the retention of plastic colloids [22], however, just for glass beads as porous medium.
Most of the current experimental works use small size spherical particles, while MPPs in the environment have various shapes and sizes that may affect their transport.Recently some studies also focused their observations on larger microplastics, such as fragments [43,45].Microplastic fragments could possibly even be detected within undisturbed sediment cores by tomographic methods [41,42].Concordantly, particles > 100 μm were found to be mobile in unsaturated laboratory experiments with glass beads; the infiltration depth increased with smaller particleto-grain size ratios [45].According to these authors, 100 μm sized MPPs can penetrate to an average depth of 0.13 m in coarse sand and fine gravel, whereas 1000 μm sized MPPs are only detectable to 0.01 m in coarse sand and 0.05 m in fine gravel.
Despite an increasing number of studies the source, transport and fate of MPPs in streambed sediments remain unclear [44].To gain a better understanding of the transport and retention of irregularly shaped microplastics in water saturated sediments we investigated and statistically analysed the abundance, distribution and retention of secondary polystyrene particles (100-2000 μm) in column experiments, with regard to dependence on flow rate, duration of infiltration event, particle size and sediment grain size.

Sediment and microplastic particles
As an analogue to coarse-grained river sediments, medium gravel (d 10 /d 50 /d 90 : 5.60/6.60/7.80mm) and coarse sand (d 10 /d 50 /d 90 : 1.13/1.51/1.90mm) was used in the column experiments as porous media (Fig. 1a).Before use, the rounded, pH neutral, organic free sediments (provided by supplier Steinlando GmbH & Co. KG according to DIN 53200) were rinsed with tap water in stainless steel sieves with a mesh size of 2 mm (medium gravel) or 1 mm (coarse sand) until the water was clear and then they were dried at 100 °C for 24 h.
Fragments of solid, fluorescent polystyrene (PS), that offer a simple and reliable detection with fluorescence microscopy [23], were produced by cryomilling (Retsch UZM 100) red fluorescent polystyrene (ρ ≈ 1.06 g cm −3 ) granulate (Magic Pyramid) and fractionated into the size classes 100-250 μm, 250-500 μm, 500-1000 μm and 1000-2000 μm using stainless steel sieves (for particle identification and size quantification see section 2.4).The resulting fragments were completely stained with the fluorescent dye.In order to remove any grinding dust adhering to the particles, the PS were rinsed with deionized water and 30% ethanol solution, dried at 60 °C and stored in the dark in amber glass bottles until use.The excitation and emission wavelength of these particles were at 550 ± 25 nm and 605 ± 70 nm, respectively.

Column preparation and operation
Stainless steel columns with a length of 50 cm and an inner diameter of 8 cm (volume ≈2.5 dm 3 ) were used for the laboratory experiments.The rinsed sediment was gradually filled into the column with a metal funnel while lightly tapping the column with a rubber mallet and saturating from below with tap water.Through the tapping and saturation, the sediment was compacted and air inclusions could escape.
A sealing ring (aramid fibers/NBR binder), a 10 μm stainless steel mesh and a 2 cm thick brass spacer ring that was filled with medium gravel were put on top of the packed column (Fig. 1c) and covered with the sealing cap with a 3/8″ two-way ball valve for PTFE hose connections.The lower sealing cap was connected with a PTFE hose with an inner diameter of 11 mm to the inlet vessel.The column experiments were operated using tap water (chemical parameters are listed in Table S2 in the supporting information).The total porosity of the sediment was calculated as the ratio of pore to column volume.The pore volume was estimated via the weight difference between the dry and the water saturated column divided by the density of water.
To mimic surface water infiltration into riverbed sediments the columns were operated in downward flow direction driven by a free gradient instead of the much more common used upward flow direction driven by a peristaltic pump [8,43].This ensures that the direction of gravity forces is in the direction of flow.Before starting the column experiments each column was flushed with at least three pore volumes of water (≈3000 ml in total) and the flow rate was adjusted by changing the pressure difference between in-and outlet (Δh, Fig. 1d).An indirect determination of the column flow rate (Q [ml min −1 ]) was performed via repeated measurements of time needed to collect 100 ml of effluent.
Once the flow rate had reached equilibrium, the flow was interrupted, the column was opened and drained by only a few millilitres, and the PS particles were added to the sediment surface.The Petri dishes used as well as tweezers and spatulas were rinsed with deionized water over the sediment surface to release any remaining MPP while at the same time refilling the column by the previously drained water volume.The column was closed again and the experiment was started with the afore established flow rate.Despite the risk of getting hydraulic artefacts in the column, this experimental approach was preferred over adding particles to the water reservoir and maintaining the predetermined flow rate, as this alternative would lead to experimental difficulties such as keeping the MPPs in dispersion, fixing a constant MPP concentration in the reservoir or avoiding possible MPP adhesion in the pipes.
A series of 11 column tests was carried out including -five different flow rates (Q) ranging between 4 ml min −1 and 40 ml min −1 in columns filled with medium gravel (mGr).These columns were named according to their flow rate and the used sediment (e.g.Q4_mGr or Q13_mGr), -two different Q of 2 ml min −1 and 4 ml min −1 in columns filled with coarse sand (cSa) (e.g.Q4_cSa), and -different durations.If not mentioned otherwise experiments were conducted for a duration until about eight pore volumes (PV) passed the columns.
If the duration was adapted to shorter times, the column names were extended by the duration of the experiment (e.g.Q10_mGr_4PV or Q10_mGr_2PV).
The flow rate was measured regularly approximately every 15 minutes.Despite all efforts and adjustments in the experimental setup the flow rates slightly varied over time for each experiment (Table 1).At the end of each experiment the columns were drained while maintaining the afore established outflow rate (though for the first columns Q38_mGr and Q20_mGr in a less stringent way than for the others).
The predefined exchange conditions agree well with water fluxes observed in surface water groundwater systems, ranging from tens of meters per day of hyporheic exchange flow in gravelly mountain streams (e.g.[48]), via a few meters per day in flood affected bank filtration systems in gravelly aquifers (e.g.[47]), to about 1 meter per day in bank filtration systems in sandy to gravelly aquifers during base flow conditions (e.g.[26]).

Sediment sampling and PS extraction
After the column transport experiments the sediment was sampled as 2.5 cm segments for the topmost 5 cm of the column and in intervals of 5 cm for the rest of the column from 5 to 50 cm.The sediment was carefully removed with a stainless-steel spoon down to 15 cm for all columns filled with medium gravel and down to 20 cm for all columns filled with coarse sand.To sample the sediment below that depth it was pushed out of the column using an individually manufactured sediment extraction pushing device (Fig. 1e).All sediment samples were transferred to 250 ml (2.5 cm intervals) or 500 ml (5 cm intervals) wide neck Erlenmeyer flasks, dried at 60 °C for at least 48 h and weighted (Table S3 in Supporting Information).
In order to separate microplastic particles from the denser matrix components we used a sodium chloride solution (NaCl) with a density of at least 1.2 g cm −3 .The Erlenmeyer flasks were approximately half-filled with the density solution, stirred for 2 min with a stainless-steel spatula (I) and treated for 2 min in an ultrasonic bath (II) to destroy soil agglomerates (steps I and II were repeated once more).Afterwards, the Erlenmeyer flasks were filled with the density solution and left for 24 h in the dark.Subsequently, all floating compounds were decanted or captured by direct vacuum filtration.The fluorescent particles were clearly visible under normal laboratory conditions.This also applies for the smallest fraction of 100-250 μm.Visual inspection of the remaining sediment after density separation revealed no visible evidence of fluorescent MPPs, thus confirming that the overwhelming majority of particles were successfully separated.
During the decantation, the Erlenmeyer flask was rotated slowly so that particles adhering to the inner wall of the flask were flushed out and the eluate was poured through a vacuum filtration unit on glass fibre (GF) filters (ROTILABO ® Typ CR263, Ø: 47 mm, mesh size: 0.7 μm) with a hand vacuum pump.The GF filters were regularly replaced in dependence of the captured number of PS and other particles to avoid particle overlays and touches while maintaining a manageable number of filters.In average 5.2 filters were used per sediment sample.In order to reduce the number of filters, this method used for Q38_mGr and Q20_mGr was replaced by vacuum filtration.
For direct vacuum filtration (PC 3001 VARIO select) the flotate in the Erlenmeyer flask was sucked off with a silicone hose up to the lower end of the flask neck and refilled with NaCl solution.This process was repeated three times.After the third filling, the separation solution was sucked off up to approx. 1 cm above the sediment.All floating compounds were captured through a filtration funnel on a 10 μm stainless steel mesh (Sartorius) and transferred to a GF filter (ROTILABO ® Typ CR263).In average 2.4 filters were used per sediment sample.All GF filters were transferred to glass petri dishes, covered with aluminium foil and placed in the drying oven at 60 °C.Average particle count for the filters of the first sediment layer (0-2.5 cm, highest MPP occurrence) was 286 ± 110.
To avoid sample contamination all materials were rinsed intensively with deionized water and ethanol solution and placed in the ultrasonic device.Two blank filters were processed within each column test (n = 20).For that 250 g of the sediment used for the column filling was transferred to a 250 ml Erlenmeyer flask and processed in the same way as the column sediment samples.

PS detection and characterisation
To detect the fluorescent PS particles all sample filters were scanned using the fluorescence imaging system ChemiDoc (Bio-Rad Laboratories, Inc).The MPPs on the filters were excited for 0.003 s using green coaxial illumination (green epi illumination) and one image per filter was recorded in the range of 602 ± 50 nm using an integrated CCD camera allowing low-noise recording with a pixel resolution of 32.7 μm.
For the characterization of particle dimensions, each recorded image was evaluated using ImageJ [33,40].Image analyses included Non-local Means Denoising [5], Thresholding, Watershed analyses and Particle analyses (Rasband, 1997(Rasband, -2018) ) in order to calculate the properties of each particle region, such as area, circumference, minimum Feret diameter and circularity.If the number of neighbouring pixels with a relevant fluorescence signal Table 1 Overview of experimental conditions.Q mean is the average flow rate, std. is the standard deviation of Q and v a is the average seepage velocity of each column ( v a = Qmean θ •A column ).The volume of the columns is ≈1000 ml * estimated value (average porosity of all gravel filled columns)

experiment porosity θ [−]
Q mean [ml min −1 ] std(Q) [ml min −1 ] v a [m d −1 ] std (v a ) [m d −1  was larger than two they were classified as one microplastic particle.Consequently, the limit of detection was roughly 100 μm.Then the minimum Feret diameter was evaluated, for which ImageJ selects the minimum distance between two parallel tangents of the particle positioned in any arbitrary angle around the x-or y-axis (minimum calliper around the particle).The minimum Feret diameter was chosen because it is expected to be the minimum size that could just allow MPPs to penetrate through the sediment pore throats.Furthermore, the reliability of the analytical MPP detection and characterisation was verified with a recovery test of 80 (20 per size class) randomly selected particles.The dimensions of those particles were additionally measured using fluorescence phase contrast microscope (KEYENCE Biozero BZ) to identify the linearity and offset between both methods (results are presented in section 3.1.).
Based on the minimum Feret diameter the relative MPP abundance (C L /C Tot ) per size class (L) and depth segment (z) was calculated as: where C L is the MPP concentration [MPP kg −1 ], C Tot is the input particle concentration in the column [MPP kg −1 ], n z, L is the number of particles of a certain segment and size class and m z is the dry mass of the sediment segment taken (Table S3 in Supporting Information).

Microplastic particle retention
The relative particle abundance versus depth after several PV were flushed through (retention profile) can be described by first-order reaction kinetics [16]: where λ L [cm −1 ] is the filter coefficient of a clean filter medium in respect to a certain particle size class (L).According to that, particle abundance in a porous medium decreases exponentially with depth (z).For the shallowest sampling interval (0-2.5 cm) z was assigned to be 0 because a large amount of MPPs were present on top of sediment whereas for all other sampling depth intervals z was ascribed to the midpoint of the segment.As particle retention is strongly dependent on particle size eq. 2 strictly applies only for a suspension uniform in size, i.e. monodisperse (particles of a certain size class).The filter coefficient was evaluated by fitting (1) eq. 2 to the measured retention profiles by minimizing the root mean square error (C L /C Tot > 0.023, below this threshold, the particle abundance does not necessarily decrease exponentially with increasing depth, but rather remains constant for low values with increasing depth, see section 3.2).
To relate the λ L to our set of experimental conditions (predictors, such as flow velocity, sediment grain and MPP size) we applied a generalized multiple linear regression model (GLM) (implemented in the R-package 'stats') [6,34]: where minFeret L is the average minimum Feret diameter of all particles within a size class, v a is the average seepage velocity, and d 50 is the mean porous media grain size.In order to obtain a dimensionless form of the predictors, they were divided by their respective mean values (mean).The coefficients β i [i = 0...3], representing the individual effect of each predictor variable on λ L , are adjusted to minimize the error term (ε) based on the sum of squared errors.For the GLM the response variable λ L is transformed with a link function (λ L * ) to accommodate for non-normal conditional distributions [13].An appropriate link function enhances the prerequisites (normal distribution of error terms, homoscedasticity, and absence of multicollinearity) of multiple linear regression models.The obtained regression coefficients allow the interpretation of the predictors' influence on the filter coefficient.

Microplastic particle detection and laboratory blanks
The recovery test of 80 randomly selected particles (20 per size class) proved the high accuracy of the analytical particle identification and quantification procedure.All fluorescent PS on the filters were correctly assigned as particles (recovery = 100%) and the false discovery rate was only 1.3% (precision = 98.7%), i.e. just on one filter a pixel group wrongly indicated the presence of a MPP 130 μm in size (constituting a single false positive).The false positive identification of a filter area seemed a result of fluorescence scattering from larger MPPs when not fully removed by the application of Non-Local Means Denoising and Thresholding.This implies that such false positive detects will not occur if there are no real MPPs on a filter.Furthermore, all identified particles were detected with an accurate size (no systematic overor underestimation).The relation between minimum (3) Feret Diameter derived with the automated particle identification and quantification procedure and individual size estimates with the fluorescence phase contrast microscope was found to be linear with a slope of 0.95 (R 2 = 0.96) (Fig. S4 in Supporting Information).Recovery rates were determined for five of the column experiments carried out (Tab.S1 in Supporting Information).Average MPP recovery rate was 84%.In total 12 fluorescent particles with an average size of 260 μm were found on 8 out of the 20 blank controls.A maximum of two particles were found in the size classes between 100 and 250 μm and 250-500 μm on any blank control.These particles could have adhered to any of the experimental columns (e.g.valves) or materials used for MPP extraction (e.g., pipes or sealings) and thus been transferred to the blanks from previously processed columns.This is more fluorescent particles on the blanks than desirable, but has to be seen in relation to the high number of detected particles on each filter.As this constitutes only 0.11% of the average number of MPPs on a filter, it follows that a contamination of the analysed samples by fluorescent PS is rather negligible.

Gravel packed columns
The relative particle abundances found in the first 2.5 cm strongly varied with the flow rate.Relative abundances of microplastic particles between 100 and 250 μm ranged between 0.47 and 0.86 for the flow rates of 38 ml min −1 and 4 ml min −1 , respectively.Generally, the abundance of larger particles in 2.5 cm and 5 cm is higher compared to smaller particles whereas this relation inverted with increasing depth (Fig. 2).All MPPs larger than 1000 μm were retained in the uppermost 5 cm (Table 2).Particles in between 500 and 1000 μm penetrated to an average depth of 20 cm whereas smaller particles were partly found over the entire length of the experimental columns (Table 2).A minor fraction of particles < 500 μm were relocated to depths below 20 cm of the gravel packed columns for flow rates equal and higher than 13 ml min −1 (Fig. 2a&b, Fig. 3c).For those conditions the relative MPP abundance from 20-50 cm remained quite constant with average values of 0.023 (73 MPP kg −1 ) and 0.012 (26 MPP kg −1 ) for size classes of 100-250 μm and 250-500 μm, respectively (Fig. 2& 3).
The maximum ratio of MPP size to sediment grain size required to allow for MPP infiltration below 20 cm in our experimental columns was about 0.08, which is somewhat smaller than the relation of 0.11 proposed by Waldschläger and Schüttrumpf [45] inferred from unsaturated experiments with spherical glass beads; i.e. the maximum MPP size that was found to infiltrate into sediments in our saturated experimental columns was just slightly lower (by a factor of 0.77) compared to Waldschläger and Schüttrumpf [45].Even if saturated and unsaturated c) Low Q Fig. 2 PS retention profiles for the gravel packed columns operated with flow rates (Q) of about 38 ml min −1 (a, Q38_mGr), 20 ml min −1 (b, Q20_ mGr) and 4 ml min −1 (c, Q4_mGr) flow conditions could not be exactly compared to each other such differences might result from deviations in water fluxes through the experimental columns (0.7 ml min −1 cm −2 in our study vs. 15.6 ml min −1 cm −2 in Waldschläger and Schüttrumpf [45]).

Sand packed columns
Particles above 500 μm were retained in the uppermost 2.5 cm of the sandy column fillings (Table 2, Fig. S5 in Supporting Information).Smaller particles infiltrated to an average depth of 15 cm.Average particle abundances of size classes of 100-250 μm and 250-500 μm at a depth of 15 cm were below 4.5 × 10 −3 (8 MPP kg −1 ) and 2.1 × 10 −3 (4 MPP kg −1 ), respectively.The maximum particle penetration depths in the sand columns were found to be higher than those of Tumwet et al. [43] who found maximum penetration depth of 2.5 cm for PVC of 125-200 μm, comparable sediment sizes (1.6-2.0 mm) and applied water fluxes (0.7 ml min −1 cm −2 , ultrapure water).Under these experimental conditions, spherical LDPE particles did not even migrate from its initial position.c) Many PV Fig. 3 PS retention profiles for the gravel packed columns operated with comparable flow rates of about 10 ml min −1 (a & b) and 13 ml min −1 (c), but with different duration of infiltration that is pore volumes (PV) passed through.Applied flow volumes were about 2 PV (a, Q10_mGr_2PV), 4 PV (b, Q10_mGr_4PV) and 8 PV (c, Q13_mGr)

Variation in duration of infiltration flow
We observed comparable spatial patterns of particle abundances for the gravel packed columns when operated with different pore volumes, that is different duration of the infiltration flow.Again, the particle abundance of all size classes strongly decreased with depth.The highest number of particles of all size classes were found in the first 2.5 cm of the sediment (Fig. 3).Slight differences were observed in the maximum particle penetration depths (Table 2), which can also be caused by not exactly identical exchange flows besides the variations in the applied flow volume.Thus, the retention of MPPs was found not to be dependent on extending the flow duration beyond 2 PV.Due to the slightly varying flow throughout the column experiments the effect of infiltration duration on the particle retention profiles could not be further quantified.However, for the studied range of variation in flow durations the experimental results so far indicate that extending the duration had a negligible effect on the retention profiles (Fig. 3), compared to the variation in column filling material and flow velocity.The reason seems to be that a steady particle retention profile quickly established within the very first exchanged pore volumes, and not much further transport of the microplastic particles occurred when already retained in the column, at least for the stable hydrochemical and hydrodynamic conditions of our experiment.

Influence of experimental parameters on the filter coefficient
The filter coefficient (λ L ) was optimised for each individual experiment in order to achieve the best fit between theoretical (eq.2) and measured particle retention profiles.The optimized filter coefficient, describing the local decrease in particle abundance with depth, ranged between 0.18 (Q38_mGr, 100-250 μm) and 1.0 (Q2_cSa, 250-500 μm) (Table 3, Fig. 4).Low values of λ L represent a low filtering effect of the sediment, i.e.MPPs are substantially relocated to greater depth, whereas high values of λ L correspond to a high filtering effect of the sediment, i.e.MPPs are retained already in the first centimetres of the sediment (Fig. 4).
The filter coefficient was found to -continuously increase with increasing MPP size, -continuously decrease with increasing flow rate, and -substantially decrease with increasing sediment particle size.
These general relations did apply with exception for the comparison between Q2_cSa and Q4_cSa.Here the filter coefficient for the flow rate of 4 ml min −1 is slightly larger than for 2 ml min −1 , not smaller.It may be Q4_cSa that deviates from the general relations of the filter coefficients due to slight unknown differences in the column sediment structure affecting the particle transport (experimental error).For Q4_cSa the particle abundance at 5-10 cm and 15-20 cm is slightly higher than the particle abundance in 2.5-5 cm depths (Fig. S5 in Supporting Information).This peculiarity also affects the value and the quality of the estimated filter coefficient for this column.The goodness of fit between theoretical and experimental retention profile for Q4_cSA, size class 100-250 μm (R 2 = 0.40) was substantially lower than for all our remaining experiments (average R 2 for all model fits = 0.86, Tab.S6 in Supporting Information).The good model fits furthermore verify that MPP abundance decreases exponentially with increasing depth, as is typical for a particle filtration.
The average difference in filter coefficient between the size class of 100-250 μm and 250-500 μm of all gravel filled and sand filled columns was 0.05 cm −1 and 0.24 cm −1 , respectively (Table 3).The filter coefficient of all gravel packed columns increased on average by 0.015 per decrease of 1 m d −1 (1.4 ml min −1 ) in flow velocity.These results extend studies previously limited to very small, primary PS spheres (< 10 μm) [3,4,9] to more realistic circumstances and show that particle transport of PS fragments strongly depend on the hydrodynamics, and the colloid and collector grain sizes.We also examined more systematically how the MPP retention (quantified by the filter coefficient) is related to dimensionless flow velocity, MPP and sediment grain sizes.GLM regression was performed using the log transformation as link function.Fig. 5 depicts the effect of the seepage velocity for different MPP size and sediment particle size on log(λ L ) [λ L in cm −1 ] .The applied GLM revealed significant effects of all predictors on the filter coefficient (p-value < 0.01, R 2 = 0.92): Because the estimated beta values hardly differ, it can be concluded that the influence of each predictor on the filter coefficient is comparable magnitude (d 50 has a slightly greater influence on the retention behavior of MPPs than v a ).The minFeret L is specific for a certain MPP size class (L) and therefore linearly influences the offset of the multiple GLM but with an overall significant effect on the filter  S7 & S8 in the Supporting Information).
Equation 4 allows to estimate retention of microplastic fragments in riverbed sediments and alluvial aquifers mainly consisting of medium sands and large gravels.In other words, by an approximate estimate of the filter coefficient (eq.4), the retention profiles of MPPs in rough porous media (eq.2) and the maximum depth of particle appearance in riverbeds can be calculated based on the prevailing environmental conditions (flow velocity, MPP and median sediment particle size).
In summary, while experimental sediment columns provide a mechanistic insight into microplastic transport, their findings should be interpreted with caution -Natural riverbeds often exhibit a wide range of sediment types, including various particle sizes and organic content that are not accounted for in the current experiments.-River flow is dynamic and influenced by various factors such as channel morphology, flow velocity, and turbulences.In laboratory settings, flow conditions are typically simplified, e.g.no lateral forces on top of the infiltration point were present that can cause horizontal particle movement or particle remobilization) -Natural riverbeds are teeming with diverse microbial and macroinvertebrate communities.Their activities, such as bioturbation and biofilm formation, can significantly alter microplastic transport patterns.-River systems experience seasonal fluctuations in flow rates, sedimentation patterns, and MPP loads.
These aspects can explain heterogeneous retention profiles of MPPs found in natural riverbed sediments [14] and the occurrence of large particles in relatively high streambed depth [15] that could not be explained by MPP transport with infiltrating surface water alone.So far eq. 4 is only valid for PS fragments in the applied range of seepage velocity, MPP size, and sediment particle size.In the experiments only organic free, homogeneous sediment materials with a narrow grain size distributions were considered.Extrapolations beyond that range be treated with caution.
Nevertheless, this study provides a mechanistic approach to estimate retention profiles and advective transport distance of MPPs in riverbed sediments caused by surface water infiltration.Our results show that river sediments represent both, a potential sink of MPP due to the high retention within the first cm of sediment but also a potential pathway of pore scale microplastics (100-500 μm) from surface waters to hyporheic zones and aquifers.The much larger proportion of MPPs are, at least temporarily, immobilized in the streambeds.Overall, this expands our ability to understand transport behaviour and distribution of MPP in river and lake bottom sediments.

Conclusions
Plastic pollution is a global and pervasive problem in freshwater ecosystems.The results of these column studies highlight that MPPs are generally retained in saturated sediments, although in dependence of flow velocity, sediment grain and MPP size.In the gravel packed columns all particles above 1000 μm were retained on the uppermost 5 cm.MPPs in between 500 and 1000 μm could infiltrate to an average depth of 20 cm whereas MPP below 500 μm in size were relocated to depth of 50 cm for high flow velocities.In the sand packed columns MPPs above 500 μm were retained in the uppermost 2.5 cm whereas MPPs below 500 μm infiltrated to an average depth of 15 cm.Nevertheless, some MPPs below 500 μm were transported down the full length of the column and thus could pass larger distances than that.This implies that MPP < 500 μm on the pore scale can indeed be advectively transferred into coarse grained sediments, can remain mobile and thus pose a potential threat to ground-and drinking water analogous to dissolved substances in induced bank filtration [26].
In general, it seems that already within 2 PV of water flow the major part of MPP retention has been occurring and afterwards no substantial relocation of MPP occurs, at least if conditions remain stable and no further MPP input happens.We observed an exponential decrease in MPP abundance with increasing depth, which significantly depends on flow sediment grain and MPP size and can be well described by first-order deep filtration theory.We therefore assume that MPP distributions in natural riverbed sediments, which deviate from this theoretical behavior, are influenced by a variety of dynamic processes such as variable sedimentation rates, remobilization, and sediment rearrangements.However, there is still an urgent need to further mechanistically unravel the role of environmental factors affecting the transport.

Fig. 1 a
Fig. 1 a Top view of columns (Ø = 8 cm) filled with medium gravel, d 50 = 6.60 mm (left) and coarse sand, d 50 = 1.51 mm (right); b Fragments of fluorescent polystyrene particles within size class 500-1000 μm (left) and 250-500 μm (right); c Top view of columns covered with brass spacer ring filled with gravel (left) as well as sealing ring and 10 μm stainless steel mesh (right); d Sketch of the column setup (not to scale; column Ø = 8 cm and L = 50 cm); and e Sketch of the self-made sediment extraction pushing device coefficient.Just considering a GLM with a single predictor revealed a significant effect of v a (R 2 = 0.64) and d 50 (R 2 = 0.47) on λ L * (details of the GLM's are shown in Table

Fig. 4
Fig.4 Comparison of PS retention profiles of the particle Size Class from 100 to 250 μm.Note y-axis is log scaled.Below a particle abundance of 0.023 (area highlighted in light grey), the observed change in particle abundance with increasing depth partly deviates from the exponential model assumption

Fig. 5
Fig.5 Relation between flow rate and filter coefficient (λ L ).The bars depict the standard deviation of the averaged flow velocities.Due to better readability, the generally quite small standard error of λ L was not displayed.These values are listed in TableS6in the supporting information

Table 2
Maximum penetration depth (depth in which particles are not present for the first time when screening downwards) for each experiment depending on particle size class

Table 3
Overview of optimized filter coefficient values of each column and particle size class, as far as applicable.The corresponding R 2 and standard error is listed in Tab.S6 in the Supporting Information may have been influenced by experimental errors (e.g.heterogeneous sediment layers) *