Sinking timescales of virtual particles subject to biofouling are more dependent on the initial size of particles than initial densityLow algal presence in subtropical gyres and minimal biofouling can contribute to particles between 1 and 0.01 mm remaining at the surfaceModeled biofouled particles of 0.1 microns start sinking almost immediately and show a global median sinking timescale of one day
Estimates suggest that one hundred times more plastic enters the ocean annually than is found at the sea surface (Van Sebille et al., 2015). It is still unknown where the plastic goes and how much resides in the other ocean reservoirs, including the sea floor (Courtene‐Jones et al., 2020), the water column (Choy et al., 2019; Koelmans et al., 2017; Ye & Andrady, 1991), the shorelines (Lebreton et al., 2019) and marine biota (F. Galgani et al., 2015; Lusher et al., 2017). By improving our understanding of the ocean plastic budget, stakeholders including policy‐makers, conservationists, and scientists can focus their attention on the reservoir with the highest concentrations. The impact of plastic pollution on aquatic species as well as human health (e.g., through consumption of seafood (Smith et al., 2018) and sea salt (Kim et al., 2018)) can then also be further investigated.
Many recent studies on the transport of plastic in the ocean have focused on what occurs at the surface (Van Sebille et al., 2020). The three‐dimensional pathway of plastic, however, is still rather poorly understood (Mountford & Morales Maqueda, 2019). Around half of the mass of plastic produced worldwide is positively buoyant and therefore should float upon entry in the ocean (Andrady, 2011; Geyer et al., 2017). Some studies have shown, however, that plastic types with initial densities lower than seawater have been found in the water column and on the ocean floor (Courtene‐Jones et al., 2020; Hidalgo‐Ruz et al., 2012; Schwarz et al., 2019). When the density of a particle exceeds the density of its surrounding seawater, it begins to sink. A few processes could induce such sinking of microplastic by increasing its density (Van Sebille et al., 2020), including entrainment in marine snow (Michels et al., 2018; Porter et al., 2018; Zhao et al., 2018), fecal pellet sinking (Cole et al., 2016; Kvale et al., 2020), migration after ingestion by marine life (Botterell et al., 2019; Lusher et al., 2017), suspension in inorganic particles or metals (e.g., ambient toxins; Richard et al. (2019)) and algal growth or biofouling (Amaral‐Zettler et al., 2020; Kooi et al., 2017; Rummel et al., 2017). Here, we focus on the latter process since biofouling is recognized as one of the main potential sinks of surface microplastics and requires improved modeled estimates (Egger, Sulu‐Gambari, et al., 2020; Nguyen et al., 2020; Van Melkebeke et al., 2020).
Biofouling is the attachment and accumulation of biological organisms on the surface of a submerged object (Long et al., 2015; Ye & Andrady, 1991). The presence of plastic provides a new substrate for microbial communities to thrive (also known as the “plastisphere”; Wright et al. (2020); Zettler et al. (2013)). This can have a negative impact on other species in the euphotic zone via (1) reduced light penetration, (2) impaired oxygen and carbon dioxide production and (3) a thickened barrier for gaseous exchanges at the air‐sea interface. Upon sinking, the plastic induces a downward flux of organic carbon and hence can also affect the ocean’s biological and carbon pump (L. Galgani & Loiselle, 2020).
In this study, we estimate the sinking characteristics of particles on a global scale and therefore focus on a numerical modeling approach, since microplastic biofouling observations are still sparse and regional, as well as logistically hard to execute. The following are examples of locations where observational studies have been conducted; Biscayne Bay (Ye & Andrady, 1991), the Bay of Bengal (Artham et al., 2009), Cape Town (Fazey & Ryan, 2016), Germany (Kaiser et al., 2017) and Monterey Bay (Choy et al., 2019). Fazey and Ryan (2016) show that 17 days are required for 50% of low‐density polyethylene sheets of 5 × 5 × 0.1 mm3 to sink due to biofouling, and their smaller samples lose buoyancy much faster than their larger ones (49 days for 50 × 50 × 0.1 mm3).
Kooi et al. (2017) propose a model (hereafter referred to as the Kooi model) to estimate the sinking timescale of biofouled microplastic particles depending on their size, density, biofilm growth, and fixed depth profiles for water density, light, salinity, temperature, and viscosity. Their results suggest that under idealized hydrodynamic conditions, due to a trade‐off between the collision frequency with algae and surface‐to‐volume ratio, the smallest particles (0.1 μm) start sinking immediately and the larger particles (0.1–10 mm) take around 30 days to start sinking. Their study explores the vertical transport of microplastic at a fixed latitude and longitude, hence the motivation behind our study to investigate the three‐dimensional transport while using more realistic hydrodynamic and biological conditions that can vary temporally and spatially.
The aim of our work is to generate a global map of the sinking characteristics of initially buoyant virtual particles subject to biofouling. To improve understanding of the fate of floating microplastic in the ocean, the sinking characteristics we explore are (1) the residence time at the surface before sinking (referred to as the sinking timescale) and (2) the time it takes to sink to a depth when its vertical velocity is no longer oriented downwards (referred to as its first sinking depth). All processes that could affect transport below the surface after the virtual particle reaches its first sinking depth are therefore beyond the scope of this current work. Furthermore, three‐dimensional advection (horizontal and vertical) is included in this study, however, any effects of turbulence on the sinking characteristics are not included.
We explore the sinking response of the virtual particles to the following parameters: particle size, particle density, three‐dimensional advection, global seawater properties (temperature, salinity and density) and biological properties (algal concentrations and growth), the mixed layer depth (MLD) and seasonality. The initial particle size and density are defined at the start of each run and the other parameters vary dynamically in time and space, affecting the particles’ size and density due to the growth of the biofilm. Regarding the particles' size, “microplastic” is now commonly defined among the scientific community as <5 mm, and here we define it as spherical particles that have a radius between 1 mm and 0.1 μm, which almost covers the same size range used in Kooi et al. (2017). We also use three initially buoyant densities (representing commonly used polymer types). We have therefore further developed the one‐dimensional (vertical) Kooi model, and present a Lagrangian analysis of the global three‐dimensional sinking characteristics of biofouled microplastic. Implications of this work include bringing us closer to understanding the potential sinking fate of “missing” floating plastic, due to size‐selective sinking of particles as a result of biofouling and large‐scale advection, that can vary in different locations of the global ocean.
The three‐dimensional transport of microplastic is simulated using the Lagrangian particle‐tracking framework, OceanParcels Version 2.1.6 (Delandmeter & van Sebille, 2019). We use the NEMO‐MEDUSA‐2.0 ORCA0083‐N06 output (Yool et al., 2013), hereafter MEDUSA, for the hydrodynamic data and biogeochemical data, which has a 1/12° global horizontal resolution and 75 depth levels. The data include the three‐dimensional ocean velocity fields, temperature, salinity, phytoplankton concentrations and primary productivity (the latter two are converted to algal concentrations and biofilm growth, explained in Section 2.2). We use five‐day averages from 2000 to 2009 (available from http://opendap4gws.jasmin.ac.uk/thredds/nemo/root/catalog.html). It is therefore important to note that we are investigating the effects of large‐scale advective features (such as upwelling, downwelling, convergence, and divergence). Previous work has shown that Lagrangian experiments using temporal means of up to nine days do not exhibit any significant degradation in these large‐scale flow characteristics (Qin et al., 2014).
Particles are released on a global 2° × 2° grid (9,620 particles) from 70°S to 80°N. The full latitudinal range is not used due to interpolation challenges when using the tripolar MEDUSA grid. The particles are released at 0.6 m (the surface depth in MEDUSA) and we simulate the particle trajectories for 90 days forward in time. The Lagrangian framework provides spatial and temporal interpolation of the fields following the C‐grid interpolation scheme (Delandmeter & van Sebille, 2019). We use the three‐dimensional fourth‐order Runge‐Kutta method with an integration time step of 30 s and the three‐dimensional position of each particle is stored every 12 h (180 time steps per simulation). Five sizes are chosen for the microplastic with a radius from 1000 to 0.1 μm (decreasing by one order of magnitude). Three different densities representing different initially buoyant plastic types are used; expanded polystyrene (EPS = 30 kgm−3), polypropylene (PP = 840 kgm−3) and low‐density polyethylene (LDPE = 920 kgm−3 ). We assume that all virtual microplastic particles are spherical and initially pristine, as in the Kooi model. It should be noted that our results are subject to the parametrizations and design of the model and that for example, the assumption of microplastic being present at the sea surface as clean and spherical must be taken into consideration when interpreting the results. To optimize efficiency (and reduce intensive computing time), we choose to run simulations for the “most typical year” from the 10 years mentioned above; 2004 (see Section S1 for a description on how we determine that year).
Two regional simulations are additionally generated in order to evaluate the effects of removing the biofouling process and removing advection on the sinking of the particles, where 25 particles are released on a 10° x 10° grid (with a 2° x 2° horizontal resolution). One region is in the algal‐rich zone of the equatorial Pacific (EqPac), from 4°S‐4°N and 140°W–148°W. The other region is in the algal‐scarce North Pacific subtropical gyre (NPSG), from 28°N‐36°N and 135°W–143°W. This latter region has been chosen since it is shown to have the highest concentrations of marine debris in the NPSG from Manta trawl samples (Egger, Nijhof, et al., 2020). We also choose to run this second simulation set‐up in the NPSG for 1000 days, from January 1, 2004, in order to test whether a longer simulation time allows for virtual microplastic particles to leave the accumulation zone. Lastly, these 1000 days simulations are repeated with horizontal advection only and vertical advection only, to isolate the effects of the different components of the three‐dimensional velocity fields.
The equations governing the biofouling and sinking of microplastic in this study are based on the Kooi model (see Kooi et al. (2017) for the detailed method). Our code to generate the results can be found here: https://doi.org/10.5281/zenodo.4543145. The sinking velocity, V s (ms−1), is dependent on the density difference between the combined plastic plus biofilm and its surrounding seawater, as well as the buoyancy of the virtual microplastic (i.e., the particle’s size and density):
The attached algal growth, dA/dt (no. m−2s−1), is dynamically estimated as:
It should be mentioned at this stage that the Kooi equations produce unrealistic results for the sinking velocities of 10 mm size particles (which reach a V s value of 20 ms−1). Hence, we only focus on particles of a 1 mm particle radius (that have a maximum V s of 0.2 ms−1 ) and smaller. Furthermore, within the context of the Kooi model, any effects of turbulence are not included in this study. This has been decided due to the fact that a wind‐driven mixing term proposed by Kukulka et al. (2012) used in previous studies (e.g., Wichmann et al. (2019)) does not hold for negatively buoyant particles, which is the case for biofouled particles that sink in our study. Parametrizations of mixing (such as wind‐driven mixing or sub‐grid processes) must therefore be further developed in future work and have been left out of the current study. Furthermore, the advantage of this set‐up for our process study is to separate the effects of advection from wind‐driven mixing. We would expect that if one were to include vertical mixing to our study, particles with a smaller sinking velocity than that of the mixing term would exhibit shorter sinking timescales (for example, the smallest particles) and could sink deeper.
Our results are represented using two quantities. The first is the sinking timescale, T s, defined as the number of days until the vertical velocity (V s from Equation 1) plus the vertical advection from MEDUSA is downward, within the simulation time of 90 days. Since the particles are all initially positively buoyant, we force any rising particles to remain at 0.6 m (the ‘surface depth’ in MEDUSA), with a vertical velocity of 0 m s−1 . This is to avoid algal attachment due to oscillations between 0 and 0.6 m (see Equation S17) that are meaningless since we treat 0.6 m as our surface. Simulations for the four seasons are first run separately, between December 2003 and November 2004 inclusive (released on December 1st, 2003 for DJF, December‐January‐February; and so on for MAM, March‐April‐May, JJA, June‐July‐August and SON, September‐October‐November). The four T s maps produced have been averaged at each initial release location (for only those simulations where particles do sink within the 90 days) in order to produce one global map of T s estimations. These T s maps are therefore plotted where the particles are first released (as opposed to the location where the particles first sink, since it would be harder to interpret particles that overlap spatially). The first sinking depth, Z s, for each individual season is defined as the depth each particle reaches until its vertical velocity is no longer downwards. We do this to highlight the initial effects of biofouling on the sinking of particles that are originally found at the surface, rather than the final position of the particles. If a submerged particle remains at a certain depth, or ascends, that part of the trajectory is not included in this study. Furthermore, if a particle continuously sinks throughout the 90 days and its vertical velocity remains downwards, its depth at 90 days is recorded as Z s. We therefore capture each particle’s trajectory as a function of depth and time until Z s, its first sinking depth.
The globally mapped sinking timescale (T s ) as a result of biofouling and three‐dimensional advection is estimated for different virtual microplastic sizes and densities (Figure 1). A particle sinks as soon as its density exceeds that of its surrounding seawater density and vertical advection is downwards. As explained in Kooi et al. (2017), the sinking timescale is a trade‐off between the particle’s size and the surface‐to‐volume ratio. The larger particles are more likely to encounter ambient algal cells, hence the collision rate is higher and the density can increase sooner. On the other hand, the smaller the object, the greater the relative surface area, meaning that the attachment of a very small number of ambient algae can cause it to sink. This concept can be visualized by the results for the smallest particles simulated here (0.1 μm in Figures 1m, 1n, and 1o) where particles start sinking almost immediately globally (median global T s = 1 day for all three densities; 30, 840 and 920 kgm‐3 ). The second‐smallest particles (1 μm in Figures 1j, 1k, and 1l) take slightly longer in some regions than the smallest particles, and then the three largest sizes (1, 0.1 and 0.01 mm in Figures 1a–1i) produce almost identical global T s maps. For these three largest sized particles, the median global T s is 40–43 days for the density particle type, LDPE (Figures 1c, 1f and 1i), 35–38 days for PP (Figures 1b, 1e and 1h) and 33–36.5 days for the least dense, EPS (Figures 1a, 1d and 1g; although these differences are not statistically significant, with a standard deviation of ±18 days for all nine global maps). Our results mirror the one dimensional Kooi et al. (2017) results in that an asymptotically shaped relationship between the log of increasing particle sizes versus the sinking timescale is reached in almost all locations.
Another result to highlight is what occurs in the five subtropical gyres. For the larger sizes in this study (1 mm–10 μm in Figures 1a–1i) the subtropical gyres are distinguished from the other regions by the fact that the particles do not sink within the simulation time of 90 days (black patches). This behavior can be attributed to the subtropical gyres being largely oligotrophic zones (Morel et al., 2010), implying that there is low biological activity (lighter green patches in Figure 3, left column). This leads to low algal collision frequencies such that the density of the particles does not exceed that of its surrounding seawater. On the other hand, since less algal attachment is required to increase the density of smaller particles (1 μm in Figure 1j–1l and 0.1 μm in Figure 1m–1o), the small particles in the subtropical gyres do sink within the 90 days (albeit later than in the other oceanic regions for 1 μm particles). Even when extending this to a regional simulation of 1000 days, only 2 of the 25 particles sink for LDPE particles of 1 mm (Figure S3). This suggests that as a result of the biofouling process, larger microplastic particles can be more present at the surface than smaller particles in subtropical gyres. This follows the size‐selective removal theory proposed by Cózar et al. (2014) and further supported by Egger, Nijhof, et al. (2020). Furthermore, for all particle sizes, the equatorial regions show the shortest sinking timescales of <10 days, which coincides with regions of higher algal concentrations (darker green patches in Figure 3, left column; regardless of the season). These subtropical and equatorial results prompt us to run regional simulations to better understand T s in regions of extreme high and low algal concentrations and with different large‐scale physical processes, such as upwelling and downwelling (see Section 3.4).
Lastly, the qualitative differences in T s as a function of size are more distinguishable than the differences in T s for different densities of microplastic. The most sensitive particle size to density differences is 1 μm; in the subtropical gyres the 1 μm particles of 30 and 840 kgm−3 take longer to sink (up to twice as long, if they sink at all within 90 days), relative to 920 kgm−3 particles (Figures 1j and 1k relative to 1l). The general mapped patterns of T s for the three densities are otherwise very similar globally for the other size classes and hence our results demonstrate that the global modeled distribution of sinking timescales of microplastic is very similar for different types of initially buoyant plastic (including EPS, PP, LDPE). This is somewhat expected due to the fact that the sizes can range over many orders of magnitude (five in this case) and the densities do not even reach two orders of magnitude (due to the definition of buoyant plastic types being below the average density of seawater: 1025 kgm−3).
For the Z s results, we use one of the seasons as an example; boreal spring (March‐April‐May; MAM). The deeper Z s (m) and sooner Z s timescale (days) for the particles of all sizes released in the northern hemisphere (red trajectories in Figure 2) can be explained by the deeper mixed layer depth and larger algal blooms during the spring in the northern Atlantic. To clarify, the Z s timescale includes both T s (days) as well as the time it takes from leaving the surface to reaching Z s. The effect of seasonality on both T s and Z s is explored further in Section 3.3.
Although as stated in the previous section, T s is almost identical globally for the three largest sizes, the time it takes to reach the first sinking depth, Z s , is not identical for those three sizes (1–0.01 mm in Figures 2a–2i). Rather, it gradually takes longer to reach Z s with a decrease in particle size. For example, the largest‐sized virtual LDPE particles can sink from the surface to around 100 m within a few days (Figure 2c), whereas the smallest‐sized virtual particles can take 40 days to reach a similar depth (Figure 2o), and the latter start sinking almost immediately around the globe. Once again, as with the T s results, the difference in Z s for the different initial densities is very small, therefore all further simulations have used the density of LDPE; 920 kgm−3.
Sinking characteristics of particles can be influenced by biological and physical features, such as ambient algal concentrations and the mixed layer depth (Figure 3 and Figure 4). Since we release the virtual particles at 0.6 m, we show surface algal concentrations to illustrate the surrounding algae available for initial attachment and growth that influences when particles start sinking (i.e., T s). The MLD is shown since waters maintain generally homogeneous properties from the surface to the base of the mixed layer. With a deeper mixed layer, for example, a deeper constant seawater density profile from the sea surface could result in sinking particles having a deeper Z s since sudden changes in surrounding density are not encountered (Figure 4).
In this study, we distinguish the subtropical gyres as the regions that maintain the lowest surface algal concentrations across all seasons (which are consistent with previous work defining the convergence zones; e.g., Figure 4d in Maximenko et al. (2012)). Using one initial particle size (1 μm) to isolate effects of seasonality on the sinking timescale (Figures 3b, 3d, 3f, and 3h), the regions with higher algal concentrations generally correspond to shorter T s (equatorial regions, Southern Ocean and northern Atlantic and Pacific), while lower concentrations correspond to a longer T s (subtropical gyres) across the seasons (this is also seen for other sizes, e.g. 0.1 mm in Figure S1). Focusing on one season, the spring is known to influence algal concentrations due to algal blooms in MAM in the northern hemisphere and SON in the southern hemisphere (Figures 3c and 3g, respectively). During boreal spring, the North Atlantic and North Pacific show the lowest T s across all seasons, with some particles sinking within 10 days in the subtropical gyres, while the southern hemisphere subtropical gyres show T s values of up to 90 days or more (Figure 3d). Similarly, during austral spring (SON) the southern hemisphere subtropical gyres show comparatively lower T s values relative to the other seasons, reducing the patches with T s > 90 days (Figure 3h).
Regarding the Z s results as a function of seasonality (Figure 4), the spring months show consistency with T s results. During MAM, 1 μm particles released in the northern hemisphere reach their Z s sooner (Figure 4d), and during SON, particles released in the southern hemisphere reach their Z s sooner (Figure 4h). This is also seen for all other sizes (for example, 0.1 mm particles in Figure S2). The maximum Z s reached for particles released in the northern and southern hemisphere is also deeper for their respective spring months relative to the winter or summer months. This once again suggests that during spring, particles sink deeper with higher surface algae present (Figure 3c and 3g, respectively) and with deeper MLDs (Figure 4c and 4g, respectively).
Although T s differences have been mentioned for different seasons, it is important to highlight that the results across the seasons are generally fairly similar. We can therefore further justify averaging the four seasons for each of the sizes in Figure 1. Furthermore, other physical and biological parameters from MEDUSA have also been analyzed, such as the euphotic layer depth and the total primary productivity (which is converted to algal growth for our model), however, they produce very minor global distribution differences across the four seasons (not shown), hence the focus on MLD and algal concentrations.
We now isolate the effects of advection and biofouling on the sinking characteristics by comparing regional results with and without these processes. Since the subtropical gyres are a recurring region of interest in this study, we choose to compare T s in the oligotrophic North Pacific subtropical gyre (NPSG) to T s in the algal‐rich equatorial Pacific (EqPac). Furthermore, the particle size that deserves a more in‐depth analysis is the 1 μm, since particles sink in the NPSG within 90 days (Figure 5a), contrary to particles larger than 1 μm (Figure 1).
In the NPSG simulation without biofouling while keeping advection (Figure 5b), T s slightly increases relative to simulations with both biofouling and advection (Figure 5a). The simulation lacking advection while keeping biofouling (Figure 5c) results in almost all of the 25 particles staying at the surface throughout the 90 days. This suggests that three‐dimensional advection is the process with the largest impact on T s for 1 μm in the NPSG. When repeating these simulations for 1000 days, and running one simulation with horizontal advection only and the other with vertical advection only, we see that vertical advection specifically plays a key role in transporting 1 μm particles below the surface in this downwelling zone (Figure S5d vs. Figure S4d). Since algal concentrations are low throughout all seasons (Figure 3) in the NPSG, the fact that removing biofouling has a lower impact on T s is somewhat expected. The Z s results further confirm this; those 1 μm particles that do sink within the 90 days of the MAM simulation, have a very shallow Z s (<2 m in Figure 6a), and removing biofouling has a very minor impact on this first sinking depth (Figure 6b). If particles of all sizes are fixed to the northern edge of the NPSG (without horizontal advection), the particles can be exposed to more algae than those that are advected toward the center of the gyre, and hence can sink after > 700 days (Figure S5).
In the EqPac, although the presence of algae is consistently high throughout the year (Figure 3), T s is almost not affected by the removal of biofouling or advection for 1 μm particles (T s < 10 days in Figures 5d–5f). When analyzing the impacts of these processes on Z s (for MAM), however, if local advection in this upwelling region is isolated (and biofouling is removed), very shallow first sinking depths are produced (Z s < 2 m; Figure 6e) whereas Z s can otherwise reach 20‐50 m when advection (mainly upwelling here) is removed (Figure 6f).
In this study, we simulate the effects of biofouling on the global sinking timescales of initially floating virtual microplastic (with a radius of 1 mm to 0.1 μm). Sinking of buoyant particles in the marine environment depends on the density difference between the particle and the surrounding seawater, as well as downward vertical advection for smaller particles. Therefore, particle properties combined with the global temporal and spatial variability of physical and biological properties can result in different sinking timescales in different regions of the ocean. The factors that have been analyzed here include the initial radius and density of microplastic, the seasonality of algal concentrations and mixed layer depth and the presence of three‐dimensional advection. We have therefore expanded on the Kooi et al. (2017) study that models the theoretical effects of biofouled particles in the vertical dimension and now map three‐dimensional, global sinking characteristics.
Our results suggest that the sinking of virtual particles subject to biofouling is largely dependent on the ambient algal concentrations as well as the initial size of the particle (rather than the initial density). The 1 μm particles seem to be the most sensitive to different particles’ initial densities, where the sinking timescale in subtropical gyres can half for the densest particles (920 kgm−3) relative to the least dense (30 kgm−3). The smallest virtual particle size analyzed here (0.1 μm) shows an almost immediate sinking timescale in the global ocean (with a median of less than a day). The smaller the particles, however, the longer they take to sink to their first sinking depth. Spring blooms (in both hemispheres) can lead to localized regions with higher surface algal concentrations and in turn shorter particle sinking timescales. Despite this, the overall qualitative distribution of T s is quite similar across seasons. We also hypothesize that the combination of high algal concentrations and deeper mixed layers can reduce the time particles take to reach Z s.
Although the median global T s matches the one dimensional results from Kooi et al. (2017) (1 day for 0.1 μm to 30–40 days for 10 μm to 1 mm), the novel aspect of our study is the global map we provide of T s. It allows us to explore different regions with different biological and physical features. The most prominent result in the global T s maps of particles between 10 μm and 1 mm is the lack of sinking after 90 days in low‐productive (oligotrophic) regions where particles have a low encounter rate with algae. Even for a simulation of 1000 days in the oligotrophic North Pacific subtropical gyre, only 2 out of 25 particles of a 1 mm size manage to sink. A more detailed analysis into the smaller particles suggests that T s and Z s of 1 μm particles are more dependent on advection than biofouling in subtropical gyres. In other regions with higher algal concentrations (e.g., the equatorial Pacific) T s is shorter than in the NPSG. Upwelling can keep particles from having a deep Z s in the EqPac, and therefore a deeper Z s occurs when biofouling is dominant. Our results hint to the fact that particles between 10 μm and 1 mm might be accumulating in the five subtropical gyres due to both physical reasons as well as biological properties (i.e., the lack of algae to make them sink) as also implied by Egger, Nijhof, et al. (2020). This supports the size‐selective removal theory of surface particles as proposed by Cózar et al. (2014) and seen in observations in Fazey and Ryan (2016), where smaller particles are found at a lower concentration at the sea surface than expected since they sink sooner.
In this study, we focus on implementing the one dimensional Kooi model in a more realistic biogeochemical and hydrodynamic setting by using temporally and spatially varying output data from MEDUSA. Both the choice of particle transport equations (e.g., drag parameterizations depending on shape, size, surface roughness, etc.) and environmental parameters (e.g., mixing strength) put constraints on the particle size range it is valid for (e.g., Monroy et al. (2017); De la Fuente et al. (2020)). Further work is required to determine the effects of these parameters as well as replacing the “pristine, spherical” particles with different shapes or fragmentation properties (e.g., DiBenedetto et al. (2019)) once they are further parametrized in the literature. Although our aim is to design a process study, the sinking timescales we present are subject to the particle release strategy we have chosen. In the real ocean, clean microplastic would rarely be found in the open ocean, and it could already be weathered, fragmented and biofouled before reaching the subtropical gyres, for example. Including turbulence in the upper layers of the water column would also have an impact on when the particles start sinking and how deep they sink (i.e., probably even sooner and deeper for smaller particles). It is also important to mention that a key result in the Kooi et al. (2017) study is the potential oscillatory behavior of microplastic as a result of biofouling and theoretical defouling below the euphotic layer depth, and this process has yet to be explored further.
To conclude, we have theoretically characterized the effects of biofouling on changing the density of initially buoyant microplastics by estimating the time it takes for the particles to sink, if at all, on a global scale and for a range of plastic sizes and densities. Improving parametrizations of marine microplastic biofouling is crucial to being one step closer to modeling the global ocean plastic budget by understanding biofouling’s role on the fate of the missing 99% of floating marine plastic debris.
DL, CK and EvS are part of the “Tracking Of Plastic In Our Seas” (TOPIOS) project, supported through funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement no. 715386). Simulations were carried out on the Dutch National e‐Infrastructure with the support of SURF Cooperative (project no. 16371 and 2019.034). CEJ is supported by a Skempton Scholarship. CL acknowledges support from the Swiss National Science Foundation under grant 174124. We would like to thank Andrew Yool for initial clarifications of MEDUSA output and Matthias Egger for fruitful discussions about the size‐selective removal theory. There is no conflict of interest among the co‐authors.
The NEMO‐MEDUSA‐2.0 ORCA0083‐N06 output for the hydrodynamic data and biogeochemical data used in this study is available from: http://opendap4gws.jasmin.ac.uk/threads/nemo/root/catalog.html. The code to run the simulations is available on Zenodo at: https://doi.org/10.5281/zenodo.4543145.