Here is the second of several animated models of bioturbation. It simulates the effect of mixing of a 5 cm thick layer by benthic organisms on a single 1 cm thick layer of yellow particles, followed by layers with 100% red particles.
The function [S,ABU] = BIOTURBATION_2(L,C,M) is a model for benthic mixing, where L is the number of sediment layers added to an initial pile of 10 layers, C is the number of particles per layer and M is the depth of homogeneous mixing. The function displays the mixed sediment column as yellow, red and blue pixels. It also returns the mixed sediment column as matrix and a vector containing the number of particles per layer after mixing.
Here is a plot of the sediment abundance after benthic mixing:
Berger, W.H., Heath, R.G. (1968) Vertical mixing in pelagic sediments. Journal of Marine Research 26, 134–143.
Trauth, M.H., Sarnthein, M., Arnold, M. (1997) Bioturbational mixing depth and carbon flux at the seafloor. Paleoceanography, 12, 517-526.
Trauth, M.H. (1998) TURBO: a dynamic-probabilistic simulation to study the effects of bioturbation on paleoceanographic time series. Computers and Geosciences, 24(5), 433-441.
Trauth, M.H. (2013) TURBO2: A MATLAB simulation to study the effects of bioturbation on paleoceanographic time series. Computers and Geosciences, 61, 1-10.