Modeling submerged biofouled microplastics and their vertical trajectories

Abstract. The fate of (micro)plastic particles in the open ocean is controlled by physical and biological processes. Here, we model the effects of biofouling on the subsurface vertical distribution of spherical, virtual plastic particles with radii of 0.01–1 mm. For the physics, four vertical velocity terms are included: advection, wind-driven mixing, tidally induced mixing, and the sinking velocity of the biofouled particle. For the biology, we simulate the attachment, growth and loss of algae on particles. We track 10,000 particles for one year in three different regions with distinct biological and physical properties: the low productivity region of the North Pacific Subtropical Gyre, the high productivity region of the Equatorial Pacific and the high mixing region of the Southern Ocean. The growth of biofilm mass in the euphotic zone and loss of mass below the euphotic zone result in the oscillatory behaviour of particles, where the larger (0.1–1.0 mm) particles have much shorter average oscillation lengths (  5000 m). Our results show that the vertical movement of particles is mainly affected by physical (wind-induced mixing) processes within the mixed layer and biological (biofilm) dynamics below the mixed layer. Furthermore, positively buoyant particles with radii of 0.01–1.0 mm can sink far below the euphotic zone and mixed layer in regions with high near-surface mixing or high biological activity. This work can easily be coupled to other models to simulate open-ocean biofouling dynamics, in order to reach a better understanding of where ocean (micro)plastic ends up.