Knowns and unknowns of plastic waste flows in the Netherlands.

Plastic entering the environment is a growing threat for ecosystems. We estimate the annual mass of known Dutch plastic waste generated and littered and where it ends up. We use two methods: (1) a material flow analysis of plastic waste separately collected from 13 economic sectors (including households, industry and imports) and estimate the amount sent to processing plants or exported and (2) a mismanagement model from observations of litter (on Dutch beaches and riverbanks) plus estimates of inadequately managed exported plastic scraps entering the environment abroad. In 2017 (the most recent complete data set available), an estimate of 1990 (±111) kilotonnes [kt] of plastic waste was separately collected. The top three plastic waste generating sectors (74% of the total) were households, clothing and textiles, and importation. Our mismanagement model estimates that 4.3-21.2 kt enters the environment annually; almost all of which occurs in foreign countries after inadequate management of imported Dutch waste. We highlight unknowns, including the source and/or destination of imported (623 kt) and exported (514 kt) plastics, plastics in non-household mixed waste streams and the plastic fraction of some separately collected waste, for example, e-waste. Our results stress the need for improved monitoring and reporting of plastic waste. Beyond the Netherlands, our recommendations could also help other high-income countries' decision-makers reach their circular economy goals.

Improved Model-Data Agreement With Strongly Eddying Ocean Simulations in the Middle-Late Eocene.

Model simulations of past climates are increasingly found to compare well with proxy data at a global scale, but regional discrepancies remain. A persistent issue in modeling past greenhouse climates has been the temperature difference between equatorial and (sub-)polar regions, which is typically much larger in simulations than proxy data suggest. Particularly in the Eocene, multiple temperature proxies suggest extreme warmth in the southwest Pacific Ocean, where model simulations consistently suggest temperate conditions. Here, we present new global ocean model simulations at 0.1° horizontal resolution for the middle-late Eocene. The eddies in the high-resolution model affect poleward heat transport and local time-mean flow in critical regions compared to the noneddying flow in the standard low-resolution simulations. As a result, the high-resolution simulations produce higher surface temperatures near Antarctica and lower surface temperatures near the equator compared to the low-resolution simulations, leading to better correspondence with proxy reconstructions. Crucially, the high-resolution simulations are also much more consistent with biogeographic patterns in endemic-Antarctic and low-latitude-derived plankton, and thus resolve the long-standing discrepancy of warm subpolar ocean temperatures and isolating polar gyre circulation. The results imply that strongly eddying model simulations are required to reconcile discrepancies between regional proxy data and models, and demonstrate the importance of accurate regional paleobathymetry for proxy-model comparisons.

Influence of Particle Size and Fragmentation on Large-Scale Microplastic Transport in the Mediterranean Sea.

Microplastic particles move three-dimensionally through the ocean, but modeling studies often do not consider size-dependent vertical transport processes. In addition, microplastic fragmentation in ocean environments remains poorly understood, despite fragments making up the majority of microplastic pollution in terms of the number of particles and despite its potential role in mass removal. Here, we first investigate the role of particle size and density on the large-scale transport of microplastics in the Mediterranean Sea and next analyze how fragmentation may affect transport and mass loss of plastics. For progressively smaller particle sizes, microplastics are shown to be less likely to be beached and more likely to reach open water. Smaller particles also generally get mixed deeper, resulting in lower near-surface concentrations of small particles despite their higher total abundance. Microplastic fragmentation is shown to be dominated by beach-based fragmentation, with ocean-based fragmentation processes likely having negligible influence. However, fragmentation remains a slow process acting on decadal time scales and as such likely does not have a major influence on the large-scale distribution of microplastics and mass loss over periods less than 3 years.

Investigating microscale patchiness of motile microbes under turbulence in a simulated convective mixed layer.

Microbes play a primary role in aquatic ecosystems and biogeochemical cycles. Spatial patchiness is a critical factor underlying these activities, influencing biological productivity, nutrient cycling and dynamics across trophic levels. Incorporating spatial dynamics into microbial models is a long-standing challenge, particularly where small-scale turbulence is involved. Here, we combine a fully 3D direct numerical simulation of convective mixed layer turbulence, with an individual-based microbial model to test the key hypothesis that the coupling of gyrotactic motility and turbulence drives intense microscale patchiness. The fluid model simulates turbulent convection caused by heat loss through the fluid surface, for example during the night, during autumnal or winter cooling or during a cold-air outbreak. We find that under such conditions, turbulence-driven patchiness is depth-structured and requires high motility: Near the fluid surface, intense convective turbulence overpowers motility, homogenising motile and non-motile microbes approximately equally. At greater depth, in conditions analogous to a thermocline, highly motile microbes can be over twice as patch-concentrated as non-motile microbes, and can substantially amplify their swimming velocity by efficiently exploiting fast-moving packets of fluid. Our results substantiate the predictions of earlier studies, and demonstrate that turbulence-driven patchiness is not a ubiquitous consequence of motility but rather a delicate balance of motility and turbulent intensity.

Limited Lateral Transport Bias During Export of Sea Surface Temperature Proxy Carriers in the Mediterranean Sea.

Some lipid-biomarker-based sea surface temperature (SST) proxies applied in the modern Mediterranean Sea exhibit large offsets from expected values, generating uncertainties in climate reconstructions. Lateral transport of proxy carriers along ocean currents prior to burial can contribute to this offset between reconstructed and expected SSTs. We perform virtual particle tracking experiments to simulate transport prior to and during sinking and derive a quantitative estimate of transport bias for alkenones and glycerol dibiphytanyl glycerol tetraethers (GDGTs), which form the basis of the UK' 37 and TEX86 paleothermometers, respectively. We use a simple 30-day surface advection scenario and sinking speeds appropriate for the export of various proxy carriers (6, 12, 25, 50, 100, 250, 500, and 1000 md-1). For the assessed scenarios, lateral transport bias is generally small (always <0.85°C) within the Mediterranean Sea and does not substantially contribute to uncertainties in UK' 37- or TEX86-based SSTs.

Identifying Marine Sources of Beached Plastics Through a Bayesian Framework: Application to Southwest Netherlands.

Beaches are thought to be a large reservoir for marine plastics. To protect vulnerable beaches, it is advantageous to have information on the sources of this plastic. Here, we develop a universally applicable Bayesian framework to map sources of plastic arriving on a specific beach. In this framework, we combine Lagrangian backtracking simulations of drifting particles with estimates of plastic input from coastlines, rivers and fisheries. The advantage over traditional Lagrangian simulations is that the Bayesian framework can consider information on known sources, and thus facilitates spatiotemporal source attribution for plastic arriving at the specified beach. We show that the main sources for our target beach in southwest Netherlands are the east coast of the UK, the Dutch coast, the English Channel (fisheries) and the Thames, Seine, Rhine and Trieux (rivers). We also show that floating time is a major uncertainty in source attribution using backtracking.

Simulating Lagrangian Subgrid-Scale Dispersion on Neutral Surfaces in the Ocean.

To capture the effects of mesoscale turbulent eddies, coarse-resolution Eulerian ocean models resort to tracer diffusion parameterizations. Likewise, the effect of eddy dispersion needs to be parameterized when computing Lagrangian pathways using coarse flow fields. Dispersion in Lagrangian simulations is traditionally parameterized by random walks, equivalent to diffusion in Eulerian models. Beyond random walks, there is a hierarchy of stochastic parameterizations, where stochastic perturbations are added to Lagrangian particle velocities, accelerations, or hyper-accelerations. These parameterizations are referred to as the first, second and third order "Markov models" (Markov-N), respectively. Most previous studies investigate these parameterizations in two-dimensional setups, often restricted to the ocean surface. On the other hand, the few studies that investigated Lagrangian dispersion parameterizations in three dimensions, where dispersion is largely restricted to neutrally buoyant surfaces, have focused only on random walk (Markov-0) dispersion. Here, we present a three-dimensional isoneutral formulation of the Markov-1 model. We also implement an anisotropic, shear-dependent formulation of random walk dispersion, originally formulated as a Eulerian diffusion parameterization. Random walk dispersion and Markov-1 are compared using an idealized setup as well as more realistic coarse and coarsened (50 km) ocean model output. While random walk dispersion and Markov-1 produce similar particle distributions over time when using our ocean model output, Markov-1 yields Lagrangian trajectories that better resemble trajectories from eddy-resolving simulations. Markov-1 also yields a smaller spurious dianeutral flux.

Global Modeled Sinking Characteristics of Biofouled Microplastic.

Microplastic debris ending up at the sea surface has become a known major environmental issue. However, how microplastic particles move and when they sink in the ocean remains largely unknown. Here, we model microplastic subject to biofouling (algal growth on a substrate) to estimate sinking timescales and the time to reach the depth where particles stop sinking. We combine NEMO-MEDUSA 2.0 output, that represents hydrodynamic and biological properties of seawater, with a particle-tracking framework. Different sizes and densities of particles (for different types of plastic) are simulated, showing that the global distribution of sinking timescales is largely size-dependent as opposed to density-dependent. The smallest particles we simulate (0.1 µm) start sinking almost immediately around the globe and their trajectories take the longest time to reach their first sinking depth (relative to larger particles). In oligotrophic subtropical gyres with low algal concentrations, particles between 1 and 0.01 mm do not sink within the simulation time of 90 days. This suggests that in addition to the comparatively well-known physical processes, biological processes might also contribute to the accumulation of floating plastic (of 1-0.01 mm) in subtropical gyres. Particles of 1 µm in the gyres start sinking largely due to vertical advection, whereas in the equatorial Pacific they are more dependent on biofouling. The qualitative impacts of seasonality on sinking timescales are small, however, localized sooner sinking due to spring algal blooms is seen. This study maps processes that affect the sinking of virtual microplastic globally, which could ultimately impact the ocean plastic budget.

Resolution dependency of sinking Lagrangian particles in ocean general circulation models.

Any type of non-buoyant material in the ocean is transported horizontally by currents during its sinking journey. This lateral transport can be far from negligible for small sinking velocities. To estimate its magnitude and direction, the material is often modelled as a set of Lagrangian particles advected by current velocities that are obtained from Ocean General Circulation Models (OGCMs). State-of-the-art OGCMs are strongly eddying, similar to the real ocean, providing results with a spatial resolution on the order of 10 km on a daily frequency. While the importance of eddies in OGCMs is well-appreciated in the physical oceanographic community, other marine research communities may not. Further, many long term climate modelling simulations (e.g. in paleoclimate) rely on lower spatial resolution models that do not capture mesoscale features. To demonstrate how much the absence of mesoscale features in low-resolution models influences the Lagrangian particle transport, we simulate the transport of sinking Lagrangian particles using low- and high-resolution global OGCMs, and assess the lateral transport differences resulting from the difference in spatial and temporal model resolution. We find major differences between the transport in the non-eddying OGCM and in the eddying OGCM. Addition of stochastic noise to the particle trajectories in the non-eddying OGCM parameterises the effect of eddies well in some cases (e.g. in the North Pacific gyre). The effect of a coarser temporal resolution (once every 5 days versus monthly) is smaller compared to a coarser spatial resolution (0.1° versus 1° horizontally). We recommend to use sinking Lagrangian particles, representing e.g. marine snow, microplankton or sinking plastic, only with velocity fields from eddying Eulerian OGCMs, requiring high-resolution models in e.g. paleoceanographic studies. To increase the accessibility of our particle trace simulations, we launch planktondrift.science.uu.nl, an online tool to reconstruct the surface origin of sedimentary particles in a specific location.

Closing the Mediterranean Marine Floating Plastic Mass Budget: Inverse Modeling of Sources and Sinks.

Estimates of plastic inputs into the ocean are orders of magnitude larger than what is found in the surface waters. This can be due to discrepancies in the sources of plastic released into the ocean but can also be explained by the fact that it is not well-known what the most dominant sinks of marine plastics are and on what time scales these operate. To get a better understanding on possible sources and sinks, an inverse modeling methodology is presented here for a Lagrangian ocean model, estimating floating plastic quantities in the Mediterranean Sea. Field measurements of plastic concentrations in the Mediterranean are used to inform parametrizations defining various sources of marine plastics and removal of plastic particles because of beaching and sinking. The parameters of the model are found using inverse modeling, by comparison of model results and measurements of floating plastic concentrations. Time scales for the sinks are found, and likely sources of plastics can be ranked in importance. A new mass balance is made for floating plastics in the Mediterranean: for 2015, there is an estimated input of 2100-3400 tonnes, and of plastics released since 2006, about 170-420 tonnes remain afloat in the surface waters, 49-63% ended up on coastlines, and 37-51% have sunk down.

Influence of Barotropic Tidal Currents on Transport and Accumulation of Floating Microplastics in the Global Open Ocean.

Floating plastic debris is an increasing source of pollution in the world's oceans. Numerical simulations using models of ocean currents give insight into the transport and distribution of microplastics in the oceans, but most simulations do not account for the oscillating flow caused by global barotropic tides. Here, we investigate the influence of barotropic tidal currents on the transport and accumulation of floating microplastics, by numerically simulating the advection of virtual plastic particles released all over the world's oceans and tracking these for 13 years. We use geostrophic and surface Ekman currents from GlobCurrent and the currents caused by the four main tidal constituents (M 2 , S 2 , K 1 , and O 1 ) from the FES model. We analyze the differences between the simulations with and without the barotropic tidal currents included, focusing on the open ocean. In each of the simulations, we see that microplastic accumulates in regions in the subtropical gyres, which is in agreement with observations. The formation and location of these accumulation regions remain unaffected by the barotropic tidal currents. However, there are a number of coastal regions where we see differences when the barotropic tidal currents are included. Due to uncertainties of the model in coastal regions, further investigation is required in order to draw conclusions in these areas. Our results suggest that, in the global open ocean, barotropic tidal currents have little impact on the transport and accumulation of floating microplastic and can thus be neglected in simulations aimed at studying microplastic transport in the open ocean.

Phytoplankton thermal responses adapt in the absence of hard thermodynamic constraints.

To better predict how populations and communities respond to climatic temperature variation, it is necessary to understand how the shape of the response of fitness-related rates to temperature evolves (the thermal performance curve). Currently, there is disagreement about the extent to which the evolution of thermal performance curves is constrained. One school of thought has argued for the prevalence of thermodynamic constraints through enzyme kinetics, whereas another argues that adaptation can-at least partly-overcome such constraints. To shed further light on this debate, we perform a phylogenetic meta-analysis of the thermal performance curves of growth rate of phytoplankton-a globally important functional group-controlling for environmental effects (habitat type and thermal regime). We find that thermodynamic constraints have a minor influence on the shape of the curve. In particular, we detect a very weak increase of maximum performance with the temperature at which the curve peaks, suggesting a weak "hotter-is-better" constraint. Also, instead of a constant thermal sensitivity of growth across species, as might be expected from strong constraints, we find that all aspects of the thermal performance curve evolve along the phylogeny. Our results suggest that phytoplankton thermal performance curves adapt to thermal environments largely in the absence of hard thermodynamic constraints.

Early Last Interglacial ocean warming drove substantial ice mass loss from Antarctica.

The future response of the Antarctic ice sheet to rising temperatures remains highly uncertain. A useful period for assessing the sensitivity of Antarctica to warming is the Last Interglacial (LIG) (129 to 116 ky), which experienced warmer polar temperatures and higher global mean sea level (GMSL) (+6 to 9 m) relative to present day. LIG sea level cannot be fully explained by Greenland Ice Sheet melt (∼2 m), ocean thermal expansion, and melting mountain glaciers (∼1 m), suggesting substantial Antarctic mass loss was initiated by warming of Southern Ocean waters, resulting from a weakening Atlantic meridional overturning circulation in response to North Atlantic surface freshening. Here, we report a blue-ice record of ice sheet and environmental change from the Weddell Sea Embayment at the periphery of the marine-based West Antarctic Ice Sheet (WAIS), which is underlain by major methane hydrate reserves. Constrained by a widespread volcanic horizon and supported by ancient microbial DNA analyses, we provide evidence for substantial mass loss across the Weddell Sea Embayment during the LIG, most likely driven by ocean warming and associated with destabilization of subglacial hydrates. Ice sheet modeling supports this interpretation and suggests that millennial-scale warming of the Southern Ocean could have triggered a multimeter rise in global sea levels. Our data indicate that Antarctica is highly vulnerable to projected increases in ocean temperatures and may drive ice-climate feedbacks that further amplify warming.

Micro- and Nanoplastics in Alpine Snow: A New Method for Chemical Identification and (Semi)Quantification in the Nanogram Range.

We present a new method for chemical characterization of micro- and nanoplastics based on thermal desorption-proton transfer reaction-mass spectrometry. The detection limit for polystyrene (PS) obtained is <1 ng of the compound present in a sample, which results in 100 times better sensitivity than those of previously reported by other methods. This allows us to use small volumes of samples (1 mL) and to carry out experiments without a preconcentration step. Unique features in the high-resolution mass spectrum of different plastic polymers make this approach suitable for fingerprinting, even when the samples contain mixtures of other organic compounds. Accordingly, we got a positive fingerprint of PS when just 10 ng of the polymer was present within the dissolved organic matter of snow. Multiple types of microplastics (polyethylene terephthalate (PET), polyvinyl chloride, and polypropylene carbonate), were identified in a snowpit from the Austrian Alps; however, only PET was detected in the nanometer range for both snowpit and surface snow samples. This is in accordance with other publications showing that the dominant form of airborne microplastics is PET fibers. The presence of nanoplastics in high-altitude snow indicates airborne transport of plastic pollution with environmental and health consequences yet to be understood.

Laboratory Measurements of the Wave-Induced Motion of Plastic Particles: Influence of Wave Period, Plastic Size and Plastic Density.

The transport of plastic particles from inland sources to the oceans garbage patches occurs trough coastal regions where the transport processes depend highly on wave-induced motions. In this study, experimental measurements of the plastic particles wave-induced Lagrangian drift in intermediate water depth are presented investigating the influence of the wave conditions, particle size and density on the motion of relatively large plastic particles. A large influence of the particle density is observed causing particles to float or sink for relative densities lower and larger than water respectively. The measured net drift of the floating particles correlates well with theoretical solutions for particle Stokes drift, where the net drift is proportional to the square of the wave steepness. Floating particles remain at the free water surface because of buoyancy and no evidence of any other influence of particle inertia on the net drift is observed. Nonfloating particles move close to the bed with lower velocity magnitudes than the floating particles' motion at the free surface. The drift of nonfloating particles reduces with decreasing wave number, and therefore wave steepness.

Circadian clock helps cyanobacteria manage energy in coastal and high latitude ocean.

The circadian clock coordinates cellular functions over the diel cycle in many organisms. The molecular mechanisms of the cyanobacterial clock are well characterized, but its ecological role remains a mystery. We present an agent-based model of Synechococcus (harboring a self-sustained, bona fide circadian clock) that explicitly represents genes (e.g., kaiABC), transcripts, proteins, and metabolites. The model is calibrated to data from laboratory experiments with wild type and no-clock mutant strains, and it successfully reproduces the main observed patterns of glycogen metabolism. Comparison of wild type and no-clock mutant strains suggests a main benefit of the clock is due to energy management. For example, it inhibits glycogen synthesis early in the day when it is not needed and energy is better used for making the photosynthesis apparatus. To explore the ecological role of the clock, we integrate the model into a dynamic, three-dimensional global circulation model that includes light variability due to seasonal and diel incident radiation and vertical extinction. Model output is compared with field data, including in situ gene transcript levels. We simulate cyanobaceria with and without a circadian clock, which allows us to quantify the fitness benefit of the clock. Interestingly, the benefit is weakest in the low latitude open ocean, where Prochlorococcus (lacking a self-sustained clock) dominates. However, our attempt to experimentally validate this testable prediction failed. Our study provides insights into the role of the clock and an example for how models can be used to integrate across multiple levels of biological organization.

Environmental versus operational drivers of drifting FAD beaching in the Western and Central Pacific Ocean.

In an effort to increase purse seine fishing efficiency for tropical tunas, over 30,000 drifting Fish Aggregating Devices (dFADs) are deployed every year by fishers in the Western and Central Pacific Ocean (WCPO). The use of dFADs also impacts ecosystems, in particular through marine pollution and dFAD beaching. This paper presents the first estimate of dFAD beaching events in the WCPO (>1300 in 2016-2017) and their distribution. Lagrangian simulations of virtual dFADs, released subject to contrasting deployment distributions, help us determine the relative importance of operational versus environmental drivers of dFADs drifting to beaching areas. The highest levels of beaching, occurring on Papua New Guinea and Solomon Islands, are likely a result of the prevailing westward oceanic circulation and subsequent local processes driving dFADs towards land. Similarly, high beaching rates in Tuvalu appear to be due to the general circulation of the WCPO. In contrast, beaching in Kiribati Gilbert Islands appear to be more strongly related to dFAD deployment strategy. These findings indicate that reducing beaching events via changes in deployment locations may be difficult. As such, management approaches combining dFAD deployment limits, the use of biodegradable dFADs, recoveries at-sea close to sensitive areas and/or beached dFAD removal should be considered.

The Role of Ekman Currents, Geostrophy, and Stokes Drift in the Accumulation of Floating Microplastic.

Floating microplastic in the oceans is known to accumulate in the subtropical ocean gyres, but unclear is still what causes that accumulation. We investigate the role of various physical processes, such as surface Ekman and geostrophic currents, surface Stokes drift, and mesoscale eddy activity, on the global surface distribution of floating microplastic with Lagrangian particle tracking using GlobCurrent and WaveWatch III reanalysis products. Globally, the locations of microplastic accumulation (accumulation zones) are largely determined by the Ekman currents. Simulations of the North Pacific and North Atlantic show that the locations of the modeled accumulation zones using GlobCurrent Total (Ekman+Geostrophic) currents generally agree with observed microplastic distributions in the North Pacific and with the zonal distribution in the North Atlantic. Geostrophic currents and Stokes drift do not contribute to large-scale microplastic accumulation in the subtropics, but Stokes drift leads to increased microplastic transport to Arctic regions. Since the WaveWatch III Stokes drift and GlobCurrent Ekman current data sets are not independent, combining Stokes drift with the other current components leads to an overestimation of Stokes drift effects and there is therefore a need for independent measurements of the different ocean circulation components. We investigate whether windage would be appropriate as a proxy for Stokes drift but find discrepancies in the modeled direction and magnitude. In the North Pacific, we find that microplastic tends to accumulate in regions of relatively low eddy kinetic energy, indicating low mesoscale eddy activity, but we do not see similar trends in the North Atlantic.

Risk assessment of plastic pollution on marine diversity in the Mediterranean Sea.

Plastic marine pollution is an increasing threat to global marine diversity. Quantifying this threat is particularly difficult and complex, especially when evaluating multiple species with different ecological requirements. Here, we examine the semi-enclosed basin of the Mediterranean Sea where the inputs of plastic pollution and its impact on marine diversity are still widely unknown. Eighty-four species from six taxonomic classes were evaluated to assess the risk of ingesting plastic marine debris, integrating inter-specific factors such as plastic exposure rates and life history traits (e.g., motility, habitat, and body size). Species were modelled within a spatial context to identify and estimate their exposure to plastic ingestion across the Mediterranean Sea using literature data, species distribution maps and plastic dispersion models. Our approach identified hotspots for the risk of plastic ingestion across multiple taxa in the Mediterranean Sea, highlighting that coastal species are at higher risk of ingesting plastic in the marine environment than open-sea species. The plastic exposure analysis indicated that species with larger home ranges were more at risk of exposure with increased distances while local species were more likely to be exposed to plastic closer to the centre of their home range location. The approach used in this study can be applied to support management and mitigation efforts throughout the Mediterranean Sea and in other geographic regions to minimize the impact of plastic pollution on marine diversity.

Plastics in sea surface waters around the Antarctic Peninsula.

Although marine plastic pollution has been the focus of several studies, there are still many gaps in our understanding of the concentrations, characteristics and impacts of plastics in the oceans. This study aimed to quantify and characterize plastic debris in oceanic surface waters of the Antarctic Peninsula. Sampling was done through surface trawls, and mean debris concentration was estimated at 1,794 items.km-2 with an average weight of 27.8 g.km-2. No statistical difference was found between the amount of mesoplastics (46%) and microplastics (54%). We found hard and flexible fragments, spheres and lines, in nine colors, composed mostly of polyurethane, polyamide, and polyethylene. An oceanographic dispersal model showed that, for at least seven years, sampled plastics likely did not originate from latitudes lower than 58°S. Analysis of epiplastic community diversity revealed bacteria, microalgae, and invertebrate groups adhered to debris. Paint fragments were present at all sampling stations and were approximately 30 times more abundant than plastics. Although paint particles were not included in plastic concentration estimates, we highlight that they could have similar impacts as marine plastics. We call for urgent action to avoid and mitigate plastic and paint fragment inputs to the Southern Ocean.