Empirical Lagrangian parametrization for wind-driven mixing of buoyant particles at the ocean surface

Turbulent mixing is a vital component of vertical particulate transport, but ocean global circulation models (OGCMs) generally have low-resolution representations of near-surface mixing. Furthermore, turbulence data are often not provided in OGCM model output. We present 1D parametrizations of wind-driven turbulent mixing in the ocean surface mixed layer that are designed to be easily included in 3D Lagrangian model experiments. Stochastic transport is computed by Markov-0 or Markov-1 models, and we discuss the advantages and disadvantages of two vertical profiles for the vertical diffusion coefficient Kz. All vertical diffusion profiles and stochastic transport models lead to stable concentration profiles for buoyant particles, which for particles with rise velocities of 0.03 and 0.003 m s−1 agree relatively well with concentration profiles from field measurements of microplastics when Langmuir-circulation-driven turbulence is accounted for. Markov-0 models provide good model performance for integration time steps of Δt≈30 s and can be readily applied when studying the behavior of buoyant particulates in the ocean. Markov-1 models do not consistently improve model performance relative to Markov-0 models and require an additional parameter that is poorly constrained

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.Peer ReviewedPostprint (published version

Modelling size distributions of marine plastics under the influence of continuous cascading fragmentation

Field studies have shown that plastic fragments make up the majority of plastic pollution in the oceans in terms of abundance. How quickly environmental plastics fragment is not well understood, however. Here, we study this fragmentation process by considering a model which captures continuous fragmentation of particles over time in a cascading fashion. With this cascading fragmentation model we simulate particle size distributions (PSDs), specifying the abundance or mass of particles for different size classes. The fragmentation model is coupled to an environmental box model, simulating the distributions of plastic particles in the ocean, coastal waters, and on the beach. We compare the modelled PSDs to available observations, and use the results to illustrate the effect of size-selective processes such as vertical mixing in the water column and resuspension of particles from the beach into coastal waters. The model quantifies the role of fragmentation on the marine plastic mass budget: while fragmentation is a major source of (secondary) plastic particles in terms of abundance, it seems to have a minor effect on the total mass of particles larger than 0.1 mm. Future comparison to observed PSD data allow us to understand size-selective plastic transport in the environment, and potentially inform us on plastic longevity

Modelling carbon export mediated by biofouled microplastics in the Mediterranean Sea

Marine microplastics can be colonized by biofouling microbial organisms, leading to a decrease in microplastics' buoyancy. The sinking of biofouled microplastics could therefore represent a novel carbon export pathway within the ocean carbon cycle. Here, we model how microplastics are biofouled by diatoms, their consequent vertical motion due to buoyancy changes, and the interactions between particle-attached diatoms and carbon pools within the water column. We initialize our Lagrangian framework with biogeochemical data from NEMO-MEDUSA-2.0 and estimate the amount of organic carbon exported below 100 m depth starting from different surface concentrations of 1-mm microplastics. We focus on the Mediterranean Sea that is characterized by some of the world's highest microplastics concentrations and is a hotspot for biogeochemical changes induced by rising atmospheric carbon dioxide levels. Our results show that the carbon export caused by sinking biofouled microplastics is proportional to the concentration of microplastics in the sea surface layer, at least at modeled concentrations. We estimate that, while current concentrations of microplastics can modify the natural biological carbon export by < 1%, future concentrations projected under business-as-usual pollution scenarios may lead to carbon exports up to 5% larger than the baseline (1998–2012) by 2050. Areas characterized by high primary productivity, that is, the Western and Central Mediterranean, are those where microplastics-mediated carbon export results to be the highest. While highlighting the potential and quantitatively limited occurrence of this phenomenon in the Mediterranean Sea, our results call for further investigation of a microplastics-related carbon export pathway in the global ocean

Empirical Lagrangian parametrization for wind-driven mixing of buoyant particles at the ocean surface

Turbulent mixing is a vital component of vertical particulate transport, but ocean global circulation models (OGCMs) generally have low-resolution representations of near-surface mixing. Furthermore, turbulence data are often not provided in OGCM model output. We present 1D parametrizations of wind-driven turbulent mixing in the ocean surface mixed layer that are designed to be easily included in 3D Lagrangian model experiments. Stochastic transport is computed by Markov-0 or Markov-1 models, and we discuss the advantages and disadvantages of two vertical profiles for the vertical diffusion coefficient Kz. All vertical diffusion profiles and stochastic transport models lead to stable concentration profiles for buoyant particles, which for particles with rise velocities of 0.03 and 0.003gmgs-1 agree relatively well with concentration profiles from field measurements of microplastics when Langmuir-circulation-driven turbulence is accounted for. Markov-0 models provide good model performance for integration time steps of "t≈30gs and can be readily applied when studying the behavior of buoyant particulates in the ocean. Markov-1 models do not consistently improve model performance relative to Markov-0 models and require an additional parameter that is poorly constrained

Perceptual multivariate visualisation of volumetric Lagrangian fluid-flow processes

Lagrangian flow data in oceanography are highly complex, encompassing not only the underpinning Eulerian, advective, vectorial flow fields and the three-dimensional position coordinates of tracer particles but also supplementary trajectory information such as interaction radii of particles, lifecycle source-to-sink information and biochemical process data. Visualising all those data cooperatively in its three-dimensional context is a prime challenge, as it demands to present all relevant information to enable a contextual analysis of the flow process while preventing the most commonly-occurring perceptual issues of clutter, colourisation conflicts, artefacts and the lack of spatial references in fluid-flow applications. In this article, we present visualisation design approaches for 4D spatio-temporal data in their context and introduce a novel colour-mapping approach for 3D velocity tensors. The employed visualisation approach is evaluated towards perceptual adequacy and efficacy with respect to algebraic visualisation design and on an oceanographic case study. The technical and perceptual elements have further implications and applications for still-picture and animated volumetric visualisation design in related applications of the natural sciences, such as geological flow mapping

Inferring source regions and supply mechanisms of iron in the Southern Ocean from satellite chlorophyll data

Highlights • Shelf sediment iron source concentrated around coastal margins. • No large iron flux from sediments on shallow submerged plateaus in the open ocean. • Horizontal advection of iron more important than upwelling of iron at ocean fronts. • Western boundary currents major supply mechanism of iron for Sub-Antarctic Zone. Abstract Primary productivity is limited by the availability of iron over large areas of the global ocean. Changes in the supply of iron to these regions could have major impacts on primary productivity and the carbon cycle. However, source regions and supply mechanisms of iron to the global oceans remain poorly constrained. Shelf sediments are considered one of the largest sources of dissolved iron to the global ocean, and a large shelf sediment iron flux is prescribed in many biogeochemical models over all areas of bathymetry shallower than 1000 m. Here, we infer the likely location of shelf sediment iron sources in the Southern Ocean, by identifying where satellite chlorophyll concentrations are enhanced over shallow bathymetry (2 mg m−3 are only found within 50 km of a continental or island coastline. These results suggest that sedimentary iron sources only exist on continental and island shelves. Large sedimentary iron fluxes do not seem present on seamounts and submerged plateaus. Large chlorophyll blooms develop where the western boundary currents detach from the continental shelves, and turn eastward into the Sub-Antarctic Zone. Chlorophyll concentrations are enhanced along contours of sea surface height extending off the continental shelves, as shown by the trajectories of virtual water parcels in satellite altimetry data. These analyses support the hypothesis that bioavailable iron from continental shelves is entrained into western boundary currents, and advected into the Sub-Antarctic Zone along the Dynamical Subtropical Front. Our results indicate that upwelling at fronts in the open ocean is unlikely to deliver iron to the ocean surface from deep sources. Finally, we hypothesise how a reduction in sea level may have altered the distribution of shelf sediment iron sources in the Southern Ocean and increased export production over the Sub-Antarctic Zone during glacial intervals

Atlantic multi-decadal oscillation covaries with Agulhas leakage

The interoceanic transfer of seawater between the Indian Ocean and the Atlantic, ‘Agulhas leakage’, forms a choke point for the overturning circulation in the global ocean. Here, by combining output from a series of high-resolution ocean and climate models with in situ and satellite observations, we construct a time series of Agulhas leakage for the period 1870–2014. The time series demonstrates the impact of Southern Hemisphere westerlies on decadal timescales. Agulhas leakage shows a correlation with the Atlantic Multi-decadal Oscillation on multi-decadal timescales; the former leading by 15 years. This is relevant for climate in the North Atlanti

Sea surface slope as a proxy for Agulhas Current strength

The linear relation between the strength of the Agulhas Current at nominal latitude 34°S and the gradient in sea level height anomaly across the current is investigated in a 1/10° resolution regional numerical ocean model. Our results show that the strength of the current can be estimated with reasonable accuracy using altimeter data, once it has been calibrated using in-situ transport measurements. Three years of transport measurements provide a calibration with worst-case correlation R = 0.78. In that case the errors in proxy transport have a standard deviation of 9.8 Sv, compared to a 20.2 Sv standard deviation of the transport time series itself. From these results we conclude that the design of the Agulhas Current Timeseries (ACT) experiment, a three-year deployment of moorings across the Agulhas Current and along a TOPEX/Jason altimeter ground track, will likely produce a good quality multi-decadal time series of Agulhas Current strength

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 mm and 10 μm do not sink within the 90‐day simulation time. 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 mm–10 μm) 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, localised 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