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.

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.

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.

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.

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.

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.

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.

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.

Drift in ocean currents impacts intergenerational microbial exposure to temperature.

Microbes are the foundation of marine ecosystems [Falkowski PG, Fenchel T, Delong EF (2008) Science 320(5879):1034-1039]. Until now, the analytical framework for understanding the implications of ocean warming on microbes has not considered thermal exposure during transport in dynamic seascapes, implying that our current view of change for these critical organisms may be inaccurate. Here we show that upper-ocean microbes experience along-trajectory temperature variability up to 10 °C greater than seasonal fluctuations estimated in a static frame, and that this variability depends strongly on location. These findings demonstrate that drift in ocean currents can increase the thermal exposure of microbes and suggests that microbial populations with broad thermal tolerance will survive transport to distant regions of the ocean and invade new habitats. Our findings also suggest that advection has the capacity to influence microbial community assemblies, such that regions with strong currents and large thermal fluctuations select for communities with greatest plasticity and evolvability, and communities with narrow thermal performance are found where ocean currents are weak or along-trajectory temperature variation is low. Given that fluctuating environments select for individual plasticity in microbial lineages, and that physiological plasticity of ancestors can predict the magnitude of evolutionary responses of subsequent generations to environmental change [Schaum CE, Collins S (2014) Proc Biol Soc 281(1793):20141486], our findings suggest that microbial populations in the sub-Antarctic (∼40°S), North Pacific, and North Atlantic will have the most capacity to adapt to contemporary ocean warming.

Isolation by environment in the highly mobile olive ridley turtle (Lepidochelys olivacea) in the eastern Pacific.

Spatial and temporal scales at which processes modulate genetic diversity over the landscape are usually overlooked, impacting the design of conservation management practices for widely distributed species. We examine processes shaping population divergence in highly mobile species by re-assessing the case of panmixia in the iconic olive ridley turtle from the eastern Pacific. We implemented a biophysical model of connectivity and a seascape genetic analysis based on nuclear DNA variation of 634 samples collected from 27 nesting areas. Two genetically distinct populations largely isolated during reproductive migrations and mating were detected, each composed of multiple nesting sites linked by high connectivity. This pattern was strongly associated with a steep environmental gradient and also influenced by ocean currents. These findings relate to meso-scale features of a dynamic oceanographic interface in the eastern tropical Pacific (ETP) region, a scenario that possibly provides different cost-benefit solutions and selective pressures for sea turtles during both the mating and migration periods. We reject panmixia and propose a new paradigm for olive ridley turtles where reproductive isolation due to assortative mating is linked to its environment. Our study demonstrates the relevance of integrative approaches for assessing the role of environmental gradients and oceanographic currents as drivers of genetic differentiation in widely distributed marine species. This is relevant for the conservation management of species of highly mobile behaviour, and assists the planning and development of large-scale conservation strategies for the threatened olive ridley turtles in the ETP.

Threat of plastic pollution to seabirds is global, pervasive, and increasing.

Plastic pollution in the ocean is a global concern; concentrations reach 580,000 pieces per km(2) and production is increasing exponentially. Although a large number of empirical studies provide emerging evidence of impacts to wildlife, there has been little systematic assessment of risk. We performed a spatial risk analysis using predicted debris distributions and ranges for 186 seabird species to model debris exposure. We adjusted the model using published data on plastic ingestion by seabirds. Eighty of 135 (59%) species with studies reported in the literature between 1962 and 2012 had ingested plastic, and, within those studies, on average 29% of individuals had plastic in their gut. Standardizing the data for time and species, we estimate the ingestion rate would reach 90% of individuals if these studies were conducted today. Using these results from the literature, we tuned our risk model and were able to capture 71% of the variation in plastic ingestion based on a model including exposure, time, study method, and body size. We used this tuned model to predict risk across seabird species at the global scale. The highest area of expected impact occurs at the Southern Ocean boundary in the Tasman Sea between Australia and New Zealand, which contrasts with previous work identifying this area as having low anthropogenic pressures and concentrations of marine debris. We predict that plastics ingestion is increasing in seabirds, that it will reach 99% of all species by 2050, and that effective waste management can reduce this threat.

Anticipating changes to future connectivity within a network of marine protected areas.

Continental boundary currents are projected to be altered under future scenarios of climate change. As these currents often influence dispersal and connectivity among populations of many marine organisms, changes to boundary currents may have dramatic implications for population persistence. Networks of marine protected areas (MPAs) often aim to maintain connectivity, but anticipation of the scale and extent of climatic impacts on connectivity are required to achieve this critical conservation goal in a future of climate change. For two key marine species (kelp and sea urchins), we use oceanographic modelling to predict how continental boundary currents are likely to change connectivity among a network of MPAs spanning over 1000 km of coastline off the coast of eastern Australia. Overall change in predicted connectivity among pairs of MPAs within the network did not change significantly over and above temporal variation within climatic scenarios, highlighting the need for future studies to incorporate temporal variation in dispersal to robustly anticipate likely change. However, the intricacies of connectivity between different pairs of MPAs were noteworthy. For kelp, poleward connectivity among pairs of MPAs tended to increase in the future, whereas equatorward connectivity tended to decrease. In contrast, for sea urchins, connectivity among pairs of MPAs generally decreased in both directions. Self-seeding within higher-latitude MPAs tended to increase, and the role of low-latitude MPAs as a sink for urchins changed significantly in contrasting ways. These projected changes have the potential to alter important genetic parameters with implications for adaptation and ecosystem vulnerability to climate change. Considering such changes, in the context of managing and designing MPA networks, may ensure that conservation goals are achieved into the future.

Risk analysis reveals global hotspots for marine debris ingestion by sea turtles.

Plastic marine debris pollution is rapidly becoming one of the critical environmental concerns facing wildlife in the 21st century. Here we present a risk analysis for plastic ingestion by sea turtles on a global scale. We combined global marine plastic distributions based on ocean drifter data with sea turtle habitat maps to predict exposure levels to plastic pollution. Empirical data from necropsies of deceased animals were then utilised to assess the consequence of exposure to plastics. We modelled the risk (probability of debris ingestion) by incorporating exposure to debris and consequence of exposure, and included life history stage, species of sea turtle and date of stranding observation as possible additional explanatory factors. Life history stage is the best predictor of debris ingestion, but the best-fit model also incorporates encounter rates within a limited distance from stranding location, marine debris predictions specific to the date of the stranding study and turtle species. There is no difference in ingestion rates between stranded turtles vs. those caught as bycatch from fishing activity, suggesting that stranded animals are not a biased representation of debris ingestion rates in the background population. Oceanic life-stage sea turtles are at the highest risk of debris ingestion, and olive ridley turtles are the most at-risk species. The regions of highest risk to global sea turtle populations are off of the east coasts of the USA, Australia and South Africa; the east Indian Ocean, and Southeast Asia. Model results can be used to predict the number of sea turtles globally at risk of debris ingestion. Based on currently available data, initial calculations indicate that up to 52% of sea turtles may have ingested debris.

Strengthened currents override the effect of warming on lobster larval dispersal and survival.

Human-induced climate change is projected to increase ocean temperature and modify circulation patterns, with potential widespread implications for the transport and survival of planktonic larvae of marine organisms. Circulation affects the dispersal of larvae, whereas temperature impacts larval development and survival. However, the combined effect of changes in circulation and temperature on larval dispersal and survival has rarely been studied in a future climate scenario. Such understanding is crucial to predict future species distributions, anticipate ecosystem shifts and design effective management strategies. We simulate contemporary (1990s) and future (2060s) dispersal of lobster larvae using an eddy-resolving ocean model in south-eastern Australia, a region of rapid ocean warming. Here we show that the effects of changes in circulation and temperature can counter each other: ocean warming favours the survival of lobster larvae, whereas a strengthened western boundary current diminishes the supply of larvae to the coast by restricting cross-current larval dispersal. Furthermore, we find that changes in circulation have a stronger effect on connectivity patterns of lobster larvae along south-eastern Australia than ocean warming in the future climate so that the supply of larvae to the coast reduces by ~4% and the settlement peak shifts poleward by ~270 km in the model simulation. Thus, ocean circulation may be one of the dominant factors contributing to climate-induced changes of species ranges.

Studying an Agulhas ring's long-term pathway and decay with finite-time coherent sets.

Coherent sets in dynamical systems are regions in phase space that optimally "carry mass" with them under the system's evolution, so that these regions experience minimal leakage. The dominant tool for determining coherent sets is the transfer operator, which provides a complete description of Lagrangian mass transport. In this work, we combine existing transfer operator methods with a windowing scheme to study the spatial and temporal evolution of a so-called Agulhas ring: a large anticyclonic mesoscale eddy playing a key role in inter-ocean exchange of climate-relevant properties. Our focus is on ring decay over time and the windowing scheme enables us to study how the most coherent region (our estimate of the ring) varies in position and size over a period of more than two years. We compare the eddy-like structure and its spatio-temporal changes as revealed by our method and by a classical Eulerian approach.

The tropicalization of temperate marine ecosystems: climate-mediated changes in herbivory and community phase shifts.

Climate-driven changes in biotic interactions can profoundly alter ecological communities, particularly when they impact foundation species. In marine systems, changes in herbivory and the consequent loss of dominant habitat forming species can result in dramatic community phase shifts, such as from coral to macroalgal dominance when tropical fish herbivory decreases, and from algal forests to 'barrens' when temperate urchin grazing increases. Here, we propose a novel phase-shift away from macroalgal dominance caused by tropical herbivores extending their range into temperate regions. We argue that this phase shift is facilitated by poleward-flowing boundary currents that are creating ocean warming hotspots around the globe, enabling the range expansion of tropical species and increasing their grazing rates in temperate areas. Overgrazing of temperate macroalgae by tropical herbivorous fishes has already occurred in Japan and the Mediterranean. Emerging evidence suggests similar phenomena are occurring in other temperate regions, with increasing occurrence of tropical fishes on temperate reefs.

How well-connected is the surface of the global ocean?

The Ekman dynamics of the ocean surface circulation is known to contain attracting regions such as the great oceanic gyres and the associated garbage patches. Less well-known are the extents of the basins of attractions of these regions and how strongly attracting they are. Understanding the shape and extent of the basins of attraction sheds light on the question of the strength of connectivity of different regions of the ocean, which helps in understanding the flow of buoyant material like plastic litter. Using short flow time trajectory data from a global ocean model, we create a Markov chain model of the surface ocean dynamics. The surface ocean is not a conservative dynamical system as water in the ocean follows three-dimensional pathways, with upwelling and downwelling in certain regions. Using our Markov chain model, we easily compute net surface upwelling and downwelling, and verify that it matches observed patterns of upwelling and downwelling in the real ocean. We analyze the Markov chain to determine multiple attracting regions. Finally, using an eigenvector approach, we (i) identify the five major ocean garbage patches, (ii) partition the ocean into basins of attraction for each of the garbage patches, and (iii) partition the ocean into regions that demonstrate transient dynamics modulo the attracting garbage patches.

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.

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.

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.