Multiscale Brazil nut effects in bioturbated sediment

Scientific Reports volume 12, Article number: 11450 (2022) Cite this article

916 Accesses

2 Citations

16 Altmetric

Metrics details

Size segregation in granular materials is a universal phenomenon popularly known as the Brazil nut effect (BNE), from the tendency of larger nuts to end on the top of a shaken container. In nature, fast granular flows bear many similarities with well-studied mixing processes. Instead, much slower phenomena, such as the accumulation of ferromanganese nodules (FN) on the seafloor, have been attributed to the BNE but remain essentially unexplained. Here we document, for the first time, the BNE on sub-millimetre particles in pelagic sediment and propose a size segregation model for the surface mixed layer of bioturbated sediments. Our model explains the size distribution of FN seeds, pointing to a uniform segregation mechanism over sizes ranging from < 1 mm to > 1 cm, which does not depend on selective ingestion by feeding organisms. In addition to explaining FN nucleation, our model has important implications for microfossil dating and the mechanism underlying sedimentary records of the Earth's magnetic field.

If a container filled with mixed nuts is shaken, size segregation occurs, with the larger Brazil nuts ending on the top1,2. This counterintuitive phenomenon is known as the Brazil nut effect (BNE). The BNE appears in processes involving granular mixing or flow3. In simple terms, it is caused by the ability of small particles to infiltrate voids that develop preferentially beneath large particles when the interlocking structure of granular materials is disrupted during shaking or shearing1,4,5,6,7. Despite the simplicity of this principle, the BNE depends in a surprisingly complex manner on how the local structure of granular materials gets disrupted, the material cohesion, and the relative density of the constituents8.

The BNE occurs also in geological transport processes: for instance, riverbeds get stabilized by the surface accumulation or large pebbles during bedload transport9. Fast dynamics in geophysical mass flows bears many similarities with well-studied industrial granular mixing processes10,11. Much slower phenomena, such as the natural lifting of buried archaeological artifacts12, the migration of coarse debris to the sediment surface13, and the accumulation of ferromanganese nodules (FN) on the seafloor14, remain essentially unexplained, despite having been ascribed to a form of bioturbation-driven BNE, in which burrowing organisms push aside particles that are too large to be ingested13. The extreme slowness of this ‘biological pumping’ machine prevents a direct observation of the BNE, so that its existence is usually inferred by exclusion of alternate explanations, as it is the case for the relative scarcity of buried FN15,16,17. While FN nuclei18 are large enough to remain on the sediment surface and grow, we are left with the question whether objects of the size of microfossils would also be affected by the BNE. This is an important issue for radiocarbon dating, since an upward offset would make them older than the surrounding sediment, as indeed sometimes observed19,20, contrary to the negative age offset caused by the preferential dissolution of weaker shells, known as the Barker effect21.

Here we document, for the first time, the BNE occurrence on sub-millimetre microtektite fragments (Fig. 1) that have been deposited ~ 788 ka ago22 in a pelagic sediment in the Indian Ocean. Because this was an instantaneous event on the geological timescale, distinct depth distributions for different microtektite size classes represent the impulse response generated by the combined action of sediment mixing and size segregation in the surface mixed layer (SML). The observed impulse responses have been modelled with a size segregation mechanism based on the shear-induced BNE. Our model predicts the correct minimum size of FM nuclei and the microfossil age offsets required to reconcile observed discrepancies.

Examples of microtektite fragments found in the Indian Ocean core MD90-0961.

Sediment mixing in the SML is caused by bioturbation (Fig. 2a). Mixing is represented mathematically by a stochastic process in which individual sediment particles perform a biased random walk23, until they reach the bottom of the SML, where the randomizing action of bioturbation ceases (Fig. 2b). In case of random walks with finite waiting time and jump length distributions23, the concentration C of a conservative tracer (in our case, microtektites) inside a SML where each volume element undergoes sufficient (i.e., ~ 25) bioturbation events24 before being definitively buried, is governed by a simplified version of the diffusion–advection equation obtained by neglecting the porosity gradient25,26:

where t and z are the time and the depth below the sediment–water interface, respectively, \({D}_{\mathrm{t}}\) is the diffusion coefficient of the tracer particles, \({v}_{\mathrm{b}}\) the bulk burial velocity, and \({v}_{\mathrm{t}}\) an additional upward velocity of the tracer particles, due for instance to bio-advection or size segregation (Fig. 2c). The upper boundary condition \(\left({v}_{\mathrm{b}}-{v}_{\mathrm{t}}\right)C-{D}_{\mathrm{t}}{\partial }_{z}C={F}_{\mathrm{t}}/{\varphi }_{\mathrm{s}}{\rho }_{\mathrm{s}}\) at z = 0, where \({\varphi }_{\mathrm{s}}\) is the volume fraction of solids and \({\rho }_{\mathrm{s}}\) their density, is controlled by the incoming tracer flux \({F}_{\mathrm{t}}\). A microtektite input event is then described by \({F}_{\mathrm{t}}={\Phi }_{\mathrm{t}}\delta (t)\) where \(\delta (t)\) is the Dirac impulse, and \({\Phi }_{\mathrm{t}}\) the microtektite fluence.

Particle transport within the SML. (a) Burrowing and reworking by benthic organisms. (b) Individual particles (e.g., a microtektite fragment) perform a biased random walk starting from the sediment surface, until disturbances cease below the SML. (c) Under certain conditions, the random walk in (b) is governed by a diffusion–advection equation, where the advective velocity is the sum of burial, bioadvective, and size segregation velocities. (d) Depth-dependent concentration C of tracer particles that are small enough to get buried (e.g., microtektite fragments, solid line) and of ferromanganese nodules large enough to stay indefinitely in the top part or the SML (dashed line). The curve below the SML indicates the microtektite distribution inside the historical layer resulting from an instantaneous deposition process.

Bioturbation intensity declines with depth, so that \({D}_{\mathrm{t}}\) and \({v}_{\mathrm{t}}\) are unknown functions of z. In practice, different depth-dependent diffusion models yield similar fits to experimental data25,27, which means that the SML can be represented by an equivalent homogeneous layer with thickness L and constant \({D}_{\mathrm{t}}\), \({v}_{\mathrm{b}}\) and \({v}_{\mathrm{t}}\). Solution of Eq. (1) with \(C\left(0,z\right)=\delta (z)\) yields the impulse response \(\mathcal{I}\left(t\right)=C(t,L)\) of the system, which can be converted to a depth-dependent concentration profile over \(z>L\) using the age model of the sediment (Fig. 2d). The microscopic equivalent to the impulse response is a Wiener process with constant drift, starting at \(\left(t,z\right)=(\mathrm{0,0})\) and ending at \(\left(t,z\right)=({t}_{L},L)\), where \({t}_{L}\) is the transit (or escape) time28 with probability density function \(\mathcal{I}\left(t\right)\). The age T of particles found at depth \(z>L\) is a stochastic variable related to \({t}_{L}\) by \(T={t}_{L}+{t}_{\mathrm{b}}\), where \({t}_{\mathrm{b}}\) is the burial time from the bottom of the SML derived from the age model. The stochasticity of T is an important factor affecting single specimen dating29.

Solution examples (Fig. 3a) show how \({v}_{\mathrm{t}}\) increases the time needed to cross the SML, due to the reduced or inverted tracer velocity \({v}_{\mathrm{b}}-{v}_{\mathrm{t}}\). The mean transit time \(\langle {t}_{L}\rangle\), defined as the expectation of \(\mathcal{I}\left(t\right)\), diverges above a critical \({v}_{\mathrm{t}}/{v}_{\mathrm{b}}\) threshold (Fig. 3b). This threshold is close to 1 when advection is the dominant transport mechanism of the tracer particles across the SML. Diffusion ensures a non-negligible probability to escape the SML even if \({v}_{\mathrm{t}}>{v}_{\mathrm{b}}\), yielding a higher \({v}_{\mathrm{t}}/{v}_{\mathrm{b}}\) threshold that depends on the inverse Péclet number \(G={D}_{\text{s}}/L{v}_{\mathrm{b}}\) of the bulk sediment, where \({D}_{\mathrm{s}}\) is the bulk diffusion coefficient. In all cases, \(\mathcal{I}\left(t\right)\) becomes dramatically skewed as the threshold is approached, converging to a uniform distribution over t > 0. This means that size segregation tends to redistribute large particles above the stratigraphic depth corresponding to their deposition age. The grain size dependence of \(\mathcal{I}\left(t\right)\) has obvious consequences for dating. While G affects the skewness of \(\mathcal{I}\left(t\right)\), and thus the stratigraphic age of individual particles, but not the mean age26—since \(\langle {t}_{L}\rangle =L/{v}_{\mathrm{b}}\) for \({v}_{\mathrm{t}}=0\)—size segregation increases the apparent age of larger particles with respect to the bulk, individually and on average, until a meaningful stratigraphic relation is lost.

Effect of advection on distribution and age of sediment particles. (a) Modelled impulse response of tracer particles, with \({t}_{0}=L/{v}_{\mathrm{b}}\) being the mean transit time of regular sediment particles through the SML, for selected ratios between the tracer segregation velocity \({v}_{\mathrm{t}}\) and the burial velocity \({v}_{\mathrm{b}}\). (b) Transit time enhancement as a function of \({v}_{\mathrm{t}}/{v}_{\mathrm{b}}\) for selected values of the inverse Péclet number \({G}_{\mathrm{s}}\) of the bulk sediment.

The expected concentration of microtektites belonging to a given size interval \([{s}_{1},{s}_{2}]\) is governed by the model function

where \({g}_{\mathrm{t}}\) is the empirical grain size distribution determined from microtektite counts over all depths (Supplementary Figure S1), \({z}_{0}\) the depth in sediment corresponding to the time of the deposition event, and \(\mathcal{I}\) the impulse response obtained from the solution of Eq. (1). While \({v}_{\mathrm{b}}\) is derived from the age model of the sediment core, \({\Phi }_{\mathrm{t}}\), \({D}_{\mathrm{t}}\), \(L\), and \({v}_{\mathrm{t}}\) must be determined by fitting microtektite profiles for different size classes. Because conventional diffusion–advection models based on \({v}_{\mathrm{t}}=0\) provide good fits of stable tracer concentration profiles30, parameter estimates obtained with Eq. (2) cannot be entirely significant, so that additional constraints need to be applied to the size dependencies of \({D}_{\mathrm{t}}\) and \({v}_{\mathrm{t}}\).

Segregation of large particles in the SML can be driven either directly by the BNE, or indirectly by biogenic graded bedding resulting from the selective transport of finer particles through ingestion31,32,33, burrow lining34, infills35, and resuspension36. While graded bedding has been observed in sediments dominated by single benthic organisms, it is not a typical feature of regularly deposited sediment37,38. Furthermore, size segregation resulting from graded bedding is not expected to depend on particle size above the maximum dimension of ingestible particles, while, as shown later, this dependence is required to explain the formation of FN. Plastic deformation of sediment around burrowing benthic organisms39,40 is a possible BNE driving mechanism, because the deformation field includes a vertical gradient of horizontal displacement around burrow tips, which is analogous to the horizontal shearing used in many BNE experiments6,41. In this case, both the diffusivity and the segregation velocity are proportional to the shear rate6,42. Size segregation might also be driven by microbial-induced bubble formation in organic-rich sediments43.

Experiments with sorted glass beads44, which share with microtektites the lack of preferential ingestion by feeding organism31, indicate that the size dependence of tracer diffusivity is governed by a power law of the form \({D}_{\mathrm{t}}\propto {s}^{-q}\) with q ≈ 0.52. Percolation of smaller particles through random media also displays a power-law dependence on particle size45. Accordingly, we assume that the diffusion of large tracer particles with size \(s\) in a sediment with mean grain size \({s}_{0}\) is given by \({D}_{\mathrm{t}}={{D}_{\mathrm{s}}(s/{s}_{0})}^{-q}\). Granular mixing experiments show that the advective velocity of large grains is proportional to \(s/{s}_{0}-{\psi }_{\mathrm{c}}\), where \({\psi }_{\mathrm{c}}\approx\) 2.8 is a critical size ratio threshold in binary mixtures4,8. Therefore, we model the segregation velocity as \({v}_{\mathrm{t}}=(s/{s}_{0}-{\psi }_{\mathrm{c}}){\beta }_{0}{D}_{\mathrm{s}}\), with \({\beta }_{0}\) being an unknown coefficient that expresses the segregation efficiency of bioturbation. Most granular mixing studies have been performed with cohesionless particles, which are a poor analogue to fine-grained sediment. Experiments with wetted grains show that cohesiveness tends to reduce size segregation46, as long a clumping is prevented, although this effect is much less pronounced in case of non-spherical particles47. Cohesive forces tend to suppress the ability of small particles to infiltrate voids, as they cannot freely fall. However, the biasing effect of gravity, which is the primary cause of size segregation, does not cease. A similar effect has been observed for the torque experienced by magnetic particles in a cohesive, bioturbated sediment, in presence of a weak magnetic field48,49. In this case, the resulting magnetic alignment was found to be proportional to the ratio between the magnetic torque and the torques that resist particle rotation. With these considerations in mind, effects of sediment cohesiveness are entirely accounted by \({\beta }_{0}\).

A Poisson regression model has been used to fit microtektite counts (Fig. 4) for three size classes, using the above models for \({D}_{\mathrm{t}}(s)\) and \({v}_{\mathrm{t}}(s)\). Model residuals are generally compatible with counting uncertainties estimated with bootstrapping, up to few exceptions that might be explained by sediment heterogeneity. The size segregation parameters \({\beta }_{0}\) and q are both significantly different from zero at a > 99.4% confidence level (Table 1). Estimates of \({D}_{\mathrm{s}}\) and L are comprised within the typical ranges obtained from radioactive tracers for similar sediments50. The power-law exponent q ≈ 0.25 for the size dependence of \({D}_{\mathrm{s}}\) (Table 1) is smaller than the value obtained by Wheatcroft44 for 10–300 µm glass beads, possibly because most microtektites are too large to be ingested.

Microtektite distribution fits. (a) Concentration of three microtektite size classes in core MD90-0961 (dots), and Poisson fit with Eq. (2) (solid lines). The dashed lines represent the impulse response (rescaled to the same maximum value of the fits) for sediment particles of average size, as predicted by the same model. Notice the offset of large microtektites with respect to the impulse response. (b) Model residuals corresponding to the fits in (a) (dots), and 90% confidence intervals obtained from a bootstrapping simulation of microtektite counts (error bars).

Tektite profiles illustrate how size segregation offsets the age distribution of buried objects (solid lines in Fig. 4a), relative to that of regular sediment particles (dashed lines in Fig. 4a). The predicted mean age offsets for the three size classes of Fig. 4 (~ 0.16 kyr for 0.05–0.2 mm, ~ 0.87 kyr for 0.2–0.45 mm, and ~ 2.9 kyr for 0.45–0.9 mm) are comparable with the maximum positive offsets reported for radiocarbon ages20,21. The role of size segregation in the generation of positive age offsets increases with the Péclet number, producing a ‘runaway’ effect when the net burial velocity \({v}_{\mathrm{b}}-{v}_{\mathrm{t}}\) in the SML vanishes (Fig. 3). For this reason, sediments with low deposition rates are expected to be particularly prone to age offsets caused by size segregation. For instance, decreasing \({v}_{\mathrm{b}}\) to 1.5 cm/kyr for a sediment with same properties as MD90-0961 would increase the age offset of a 0.5 mm object from ~ 0.73 to ~ 10 kyr. Large positive offsets caused by the BNE can explain foraminifera age differences that cannot arise from selective dissolution alone21,51,52. On the other hand, size segregation is less sensitive to changes of the diffusion coefficient: for instance, doubling \({D}_{\mathrm{s}}\) increases the age offset of the 0.5 mm object of the above example to ~ 1.7 kyr, because the resulting increase of \({v}_{\mathrm{t}}\) is partially compensated by a decrease of Pe.

In case of a stationary flux of large tracer particles, the BNE produces a concentration gradient within the SML with a similar dependence on \({v}_{\mathrm{t}}/{v}_{\mathrm{b}}\) as \(\langle {t}_{L}\rangle\) (Supplementary Figure S2): this is because conservation of the vertical tracer flux requires a decrease of the net burial velocity to be compensated by a higher concentration. Hence, the interpretation of foraminiferal concentration variations within the SML might be biased by the BNE. For instance, the concentration of G. bulloides tests in sediments of the Oman margin53,54, which have been used to reconstruct the Indian summer monsoon during the last ~ 2000 years, might increase by 9–40% in the uppermost ~ 6 cm, if the same segregation parameters of Table 1 are assumed along with bioturbation data representative of the oxygen minimum zone in the northwest Arabian Sea (i.e., \({D}_{\mathrm{s}}\approx\) 150 cm2/kyr, L ≈ 6 cm, and \({v}_{\mathrm{b}}=\) 3–20 cm/kyr54,55).

Size segregation at sub-millimetre scales has important paleomagnetic implications because it requires a reorganization of the sediment microstructure assimilable with a truly diffusive process, which causes a reorientation of magnetic carriers. While non-local transport in the SML is hardly distinguishable from true diffusion23, the two processes affect sedimentary records of the Earth magnetic field in a drastically different manner. Upward conveyor belt transport removes material at depth and redeposits it on the sediment surface, where a so-called detrital remanent magnetization (DRM) is acquired by partial alignment of suspended particles in the magnetic field. Buried sediment, which is not affected by bioturbation in this model, would carry an intact DRM coeval with deposition age56. Local disruption of the sediment structure, on the other hand, erases the existing DRM and replaces it with a post-depositional magnetization (PDRM) younger than the age of deposition. Conventional PDRM models assume that this magnetization is acquired below the SML during early diagenesis57; however, laboratory experiments have shown that PDRM acquisition can be driven by bioturbation through the rotational component of diffusion48,49. Diffusive sediment mixing introduces a delay of the order of \(L/{v}_{\text{b}}\) in magnetostratigraphic records not affected by diagenesis. This delay is compatible with observed offsets between magnetic mineral and 10Be records of the Matuyama-Brunhes field reversal58, if appropriated estimates of L for marine sediments50 are used.

If FN nucleation and microtektite segregation have the same origin, segregation parameters estimated from microtektite profiles can be used to predict the minimum size of FN nuclei, which is 1–5 mm18. A simple growth model assumes that the size \(s\left(T\right)={s}_{\mathrm{n}}+\gamma T\) of a FN of age \(T\) increases linearly in time from the initial nucleus size \({s}_{\mathrm{n}}\) at a constant growth rate \(\gamma \approx\) 1–5 mm/Myr59. In case of a stationary flux \({F}_{0}\) of seeds with initial grain size distribution \({n}_{0}(s)\), the size distribution \(n(s)\) of growing seeds at the sediment–water interface is given by

where \(r(s)\) is the ratio between the tracer concentrations at z = L and z = 0, respectively, obtained from the steady-state solution of Eq. (1). Solutions of Eq. (3) with the size segregation parameters of Table 1 and SML properties typical of FN fields at a water depth of ~ 4000 m50,60 predict minimum seed sizes of ~ 2–3 mm (Fig. 5), which are comparable with observed ones18. Under these conditions, an increase of \({v}_{\mathrm{b}}\) from 0.5 to 0.8 cm/kyr, a decrease of \({D}_{\mathrm{s}}\) from 22 to 8 cm2/kyr, or a decrease of \({\beta }_{0}\) from 0.075 to 0.064 m−1 are sufficient to suppress the growth of ~ 2.3 mm seeds, confirming that a minimum degree of bioturbation and a sufficiently small sedimentation rate are necessary physical conditions for FN growth17. Finally, a linear dependence of \({v}_{\mathrm{t}}\) on size up to at least ~ 4 mm, well above the limit of selective transport by benthic organisms, is required to reduce the sinking probability of growing FN nuclei so that they can grow by several cm (dashed line in Fig. 5). This demonstrates that FN are kept at the sediment surface by the BNE, rather than other mechanisms like the selective removal of fine-grained sediment by bottom currents17.

Steady-state distribution of growing ferromanganese nodule seeds predicted by Eq. (3) with size segregation parameters from microtektite fits of core MD90-0961 (Table 1), a nodule grow rate of 5 mm/Myr, and L = 6 cm, \({D}_{\mathrm{s}}\) = 22 cm2/kyr and \({v}_{\mathrm{b}}\) = 0.5 cm/kyr as representative SML properties at a water depth of ~ 4000 m. Results are shown for lognormal size distributions of depositing seeds (solid lines) with σ = 0.1 and three values of µ (numbers in mm). The dashed line represents the case of µ = 2.3 mm when the linear increase of \({v}_{\mathrm{t}}(s)\) stops at s = 4 mm.

We report, for the first time, the size segregation of 0.05–0.9 mm microtektite fragments in a pelagic sediment from the Indian Ocean. The depth distribution of these fragments can be explained by a bioturbation-driven BNE in the SML. As a result, large particles experience an upwardly directed segregation velocity \({v}_{\mathrm{t}}\) relative to the bulk sediment, which increases linearly with particle size. Above a sediment-dependent size threshold (e.g., ~ 1 mm in core MD90-0961), \({v}_{\mathrm{t}}\) exceeds the burial velocity, and the probability of burial below the SML becomes small. This has two important consequences: (1) buried objects that are much larger than the mean grain size of sediment, such as microfossils, are significantly older than their stratigraphic age and tend to lose any relation with stratigraphy for sizes above the \({v}_{\mathrm{t}}={v}_{\mathrm{b}}\) threshold, and (2) > 1 mm particles tend to remain on the sediment surface for long times, serving as seeds for the growth of FN under favourable conditions. In the latter case, continuous growth further decreases the burial probability, explaining the scarcity of buried nodules. A single empirical model for the size dependence of segregation velocity and diffusivity, derived from experiments on granular mixing, explains our microtektite counting results and correctly predicts the minimum size of FN seeds, despite the > 10 orders of magnitude difference between bioturbation and laboratory time scales. Plastic deformation of sediment associated with burrowing is the most likely BNE driving mechanism.

The BNE has important implications for the fundamental understanding of phenomena that depends on sediment micromechanics, such as FN growth and paleomagnetic records, and for the interpretation of foraminifera ages and concentration variations. Large positive age offsets caused by the BNE can explain foraminifera age differences that cannot arise from selective dissolution alone. Furthermore, the BNE produces a concentration gradient within the SML, which might affect the interpretation of recent climatic variations. The effect of physical sediment properties like cohesiveness on size segregation needs to be investigated to assess the role played by the BNE in the redistribution of large particles in sediment, beyond the single example presented here.

Core MD90-0961 (5°03.71′ N, 73°52.57′ E) was collected during the SEYMAMA research cruise of the R/V Marion Dufresne in 1990. The 45-m long core was retrieved on the eastern margin of the Chagos-Maldive-Laccadive Ridge at a water depth of 2450 m and is composed of calcareous nannofossil ooze with abundant foraminifera. Typical microtektite concentrations amount to few counts per sample (~ 3 g); therefore, counts from three sampling campaigns (Supplementary Tables S1, S2, S3) have been gathered into three size classes with a total of 137, 287, and 49 counts, respectively (Table 2). Sediment preparation details are given in the Supplementary Information.

The diffusion–advection equation does not have a simple analytical solution in case of size segregation. A series expansion of of \(\mathcal{I}\left(t\right)\) was obtained with the solution approach of Guinasso and Schink26, taking the modified boundary conditions into account (Supplementary Methods).

The significance of the size segregation parameters \({\beta }_{0}\) and \(q\) have been tested using the ratio \(\Lambda\) between the likelihoods of the null hypothesis \({H}_{0}\) that segregation does not occur (\({\beta }_{0}\) = 0 or \({\beta }_{0}\) = q = 0), and of the full model \({H}_{1}\), assuming that microtektite counts are governed by Poisson statistics (Supplementary Methods).

All data supporting our findings are given in the Supplementary Tables.

Metzger, M. J., Remy, B. & Glasser, B. J. All the Brazil nuts are not on top: Vibration induced granular size segregation of binary, ternary and multi-sized mixtures. Powder Technol. 205, 42–51 (2011).

Article CAS Google Scholar

Gajjar, P. et al. Size segregation of irregular granular materials captured by time-resolved 3D imaging. Sci. Rep. 11, 8351 (2021).

Article CAS Google Scholar

Rosato, A. D., Blackmore, D. L., Zhang, N. & Lan, Y. A perspective on vibration-induced size segregation of granular materials. Chem. Eng. Sci. 57, 265–275 (2002).

Article CAS Google Scholar

Jullien, R., Meakin, P. & Pavlovitch, A. Three-dimensional model for particle-size segregation by shaking. Phys. Rev. Lett. 69, 640 (1992).

Article ADS CAS PubMed Google Scholar

Duran, J., Rajchenbach, J. & Clément, E. Arching effect model for particle size segregation. Phys. Rev. Lett. 70, 2431 (1993).

Article ADS CAS PubMed Google Scholar

Van der Vaart, K. et al. Underling asymmetry within particle size segregation. Phys. Rev. Lett. 114, 238001 (2015).

Article ADS PubMed CAS Google Scholar

Li, Z. et al. Microscopic structure and dynamics study of granular segregation mechanism by cyclic shear. Sci. Adv. 7, eabe8737 (2021).

Article ADS PubMed PubMed Central Google Scholar

Kudrolli, A. Size separation in vibrated granular matter. Rep. Prog. Phys. 67, 209–247 (2004).

Article ADS Google Scholar

Ferdowsi, B., Ortiz, C. P., Houssais, M. & Jerolmack, D. J. River-bed armouring as a granular segregation phenomenon. Nat. Commun. 8, 1363 (2017).

Article ADS PubMed PubMed Central CAS Google Scholar

Johnson, C. G. et al. Grain-size segregation and levee formation in geophysical mass flows. J. Geophys. Res. 117, F01032 (2012).

ADS Google Scholar

Gray, J. M. N. T., Gajjar, P. & Kokelaar, P. Particle-size segregation in dense granular avalances. C. R. Phys. 16, 73–85 (2015).

Article ADS CAS Google Scholar

Luria, D., Fantalkin, A., Zilberman, E. & Ben-Dor, E. Identifying the Brazil nut effect in archaeological site formation processes. Med. Geosci. Rev. 2, 267–281 (2020).

Article Google Scholar

McCave, I. Biological pumping upwards of the coarse fraction of deep-sea sediments. J. Sediment. Res. 58, 148–158 (1988).

Article Google Scholar

Piper, D. Z. & Fowler, B. New constraint on the maintenance of Mn nodules at the sediment surface. Nature 286, 880–883 (1980).

Article ADS CAS Google Scholar

Finney, B., Heath, G. R. & Lyle, M. Growth rates of manganese-rich nodules at MANOP Site H (Eastern North Pacific). Geochim. Cosmochim. Acta 48, 911–919 (1984).

Article ADS CAS Google Scholar

Pattan, J. N. & Parthiban, G. Do manganese nodules grow or dissolve after burial? Results from the Central Indian Ocean Basin. J. Asian Earth Sci. 30, 696–705 (2007).

Article ADS Google Scholar

Dutkiewicz, A., Judge, A. & Müller, R. D. Environmental predictors of deep-sea polymetallic nodule occurrence in the global ocean. Geology 48, 293–297 (2020).

Article ADS CAS Google Scholar

Von Stackelberg, U. Growth history of manganese nodules and crusts of the Peru Basin. In Manganese Mineralization: Geochemistry and Mineralogy of Terrestrial and Marine Deposits (eds. Nicholson, K., Hein, J. R., Bfihn, B. & Dasgupta S.) (Geological Society Special Publication No. 119, 1997).

Thomson, J. et al. Radiocarbon age offsets in different-sized carbonate components of deep-sea sediments. Radiocarbon 37, 91–101 (1995).

Article CAS Google Scholar

Brown, L., Cook, G. T., MacKenzie, A. B. & Thomson, J. Radiocarbon age profiles and size dependency of mixing in Northeast Atlantic sediments. Radiocarbon 43, 929–937 (2001).

Article CAS Google Scholar

Broecker, W. & Clark, E. Radiocarbon-age differences among coexisting planktic foraminifera shells: The Barker Effect. Paleoceanography 26, 2222 (2011).

ADS Google Scholar

Jourdan, F., Nomade, S., Wingate, M. T., Eroglu, E. & Deino, A. Ultraprecise age and formation temperature of the Australasian tektites constrained by 40Ar/39Ar analyses. Meteorit. Planet. Sci. 54, 2573–2591 (2019).

Article ADS CAS Google Scholar

Meysman, F. J. R., Malyuga, V. S., Boudreau, B. P. & Middelburg, J. J. A generalized stochastic approach to particle dispersal in soils and sediments. Geochim. Cosmochim. Acta 72, 3460–3478 (2008).

Article ADS CAS Google Scholar

Reed, D. C., Huang, K., Boudreau, B. P. & Meysman, F. J. R. Steady-state tracer dynamics in a lattice-automaton model of bioturbation. Geochim. Cosmochim. Acta 70, 5855–5867 (2007).

Article ADS CAS Google Scholar

Meysman, F. J. R., Malyuga, V. S., Boudreau, B. P. & Middelburg, J. J. The influence of porosity gradients on mixing coefficients in sediments. Geochim. Cosmochim. Acta 71, 961–973 (2007).

Article ADS CAS Google Scholar

Guinasso, N. & Schink, D. Quantitative estimates of biological mixing rates in abyssal sediments. J. Geophys. Res. 80, 3032–3043 (1975).

Article ADS Google Scholar

Deleersnijder, E., Beckers, J.-M. & Delhez, E. The residence time of settling particles in the surface mixed layer. Environ. Fluid Mech. 6, 25–42 (2006).

Article Google Scholar

Dominé, M. First passage time distribution of a Wiener process with drift concerning two elastic barriers. J. Appl. Probab. 33, 164–175 (1996).

Article MathSciNet MATH Google Scholar

Lougheed, B. C., Ascough, P., Dolman, A. M., Löwemark, L. & Metcalfe, B. Re-evaluating 14C dating accuracy in deep-sea sediment archives. Geochronology 2, 17–31 (2020).

Article ADS Google Scholar

Glass, B. P. & Köberl, C. Australasian microtektites and associated impact ejecta in the South China Sea and the Middle Pleistocene supereruption of Toba. Meteorit. Planet. Sci. 41, 305–326 (2006).

Article ADS CAS Google Scholar

Fornes, W. L., DeMaster, D. J. & Smith, C. R. A particle introduction experiment in Santa Catalina Basin sediments: Testing the age-dependent mixing hypothesis. J. Mar. Res. 59, 97–112 (2001).

Article Google Scholar

Shull, D. H. & Yasuda, M. Size-selective downward particle transport by cirratulid polychaetes. J. Mar. Res. 59, 453–473 (2001).

Article Google Scholar

Guieb, R. A., Jumars, P. A. & Self, R. F. L. Adhesive-based selection by a tentacle-feeding polychaete for particle size, shape and bacteria coating in silt and sand. J. Mar. Res. 62, 260–281 (2004).

Article Google Scholar

Schieber, J. & Wilson, R. D. Burrows without a trace—How meioturbation affects rock fabrics and leaves a record of meiobenthos activity in shales and mudstones. PalZ 95, 767–791 (2021).

Article Google Scholar

Herringshaw, L. G., Sherwood, O. A. & McIlroy, D. Ecosystem engineering by bioturbating polychaetes in event bed microcosms. Palaios 25, 46–58 (2010).

Article ADS Google Scholar

Zhang, L., Shang, J., He, W., You, B. & Fan, C. The role of tubificid worms (Limnodrilus hoffmeisteri) in sediment resuspension: a microcosm study. Ann. Limnol. Int. J. Limnol. 50, 253–260 (2014).

Article Google Scholar

Carlin, J. et al. Variability in shelf sedimentation in response to fluvial sediment supply and coastal erosion over the past 1,000 years in Monterey Bay, CA, United States. Front. Earth Sci. 7, 113 (2019).

Article ADS Google Scholar

Martín, J., Puig, P., Masqué, P., Palanques, A. & Sánchez-Gómez, A. Impact of bottom trawling on deep-sea sediment properties along the flanks of a submarine canyon. PLoS ONE 9, e104536 (2014).

Article ADS PubMed PubMed Central CAS Google Scholar

Dorgan, K. M. The biomechanics of burrowing and boring. J. Exp. Biol. 218, 176–183 (2015).

Article PubMed Google Scholar

Du Clos, K. T. Visualizing subsurface burrowing by the polychaete Alitta virens with particle imaging velocimetry. Limnol. Oceanogr. Methods 12, 703–712 (2014).

Article Google Scholar

Guillard, F., Forterre, Y. & Pouliquen, O. Scaling laws for segregation forces in dense sheared granular flows. J. Fluid Mech. 807, R1 (2016).

Article ADS MathSciNet MATH CAS Google Scholar

Rognon, P. & Macaulay, M. Shear-induced diffusion in dense granular fluids. Soft Matter 17, 5271 (2021).

Article ADS CAS PubMed Google Scholar

Srivastava, A., Kikuchi, K. & Ishikawa, T. Microbial Brazil nut effect. Soft Matter 17, 10428 (2021).

Article ADS CAS PubMed Google Scholar

Wheatcroft, R. A. Experimental tests for particle size-dependent bioturbation in the deep ocean. Limnol. Oceanogr. 37, 90–104 (1992).

Article ADS Google Scholar

Anderson, A. N., Crawford, J. W. & McBratney, A. B. On diffusion in fractal soil structures. Soil Sci. Soc. Am. J. 64, 19–24 (2000).

Article ADS CAS Google Scholar

Jarray, A., Shi, H., Scheper, B. J., Habibi, M. & Luding, S. Cohesion-driven mixing and segregation of dry granular media. Sci. Rep. 9, 13480 (2019).

Article ADS PubMed PubMed Central CAS Google Scholar

Zhao, L.-L., Li, Y.-W., Yang, X.-D., Jiao, Y. & DEM Hou, Q.-F. study of size segregation of wet particles under vertical vibration. Adv. Powder Technol. 30, 1386–1399 (2019).

Article Google Scholar

Egli, R. & Zhao, X. Natural remanent magnetization acquisition in bioturbated sediment: General theory and implications for relative paleointensity reconstructions. Geochem. Geophys. Geosyst. 16, 995–1016 (2015).

Article ADS Google Scholar

Zhao, X., Egli, R., Gilder, S. A. & Müller, S. Microbially assisted recording of the Earth's magnetic field in sediment. Nat. Commun. 7, 10673 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Solan, M. et al. Worldwide measurements of bioturbation intensity, ventilation rate, and the mixing depth of marine sediments. Sci. Data 6, 1–6 (2019).

Article Google Scholar

Paull, C. K., Hills, S. J., Thierstein, H. R., Bonani, G. & Wölfli, W. 14C offsets and apparently non-synchronous δ18O stratigraphies between nannofossil and foraminiferal pelagic carbonates. Quat. Res. 35, 274–290 (1991).

Article Google Scholar

Sepulcre, S., Durand, N. & Bard, E. Large 14C age offsets between the fine fraction and coexisting planktonic foraminifera in shallow Caribbean sediments. Quat. Geochronol. 38, 61–74 (2017).

Article Google Scholar

Anderson, D. M., Overpeck, J. T. & Gupta, A. K. Increase in the Asian Southwest Monsoon during the past four centuries. Science 297, 596–599 (2002).

Article ADS CAS PubMed Google Scholar

Anderson, D. M., Baulcomb, C. K., Duvivier, A. K. & Gupta, A. K. Indian summer monsoon during the last two millennia. J. Quat. Sci. 25, 911–917 (2010).

Article Google Scholar

Smith, C. R., Levin, L. A., Hoover, D. J., McMurtry, G. & Gage, J. D. Variations in bioturbation across the oxygen minimum zone in the northwest Arabian Sea. Deep-Sea Res. II(47), 227–257 (2000).

ADS Google Scholar

Katari, K., Tauxe, L. & King, J. A reassessment of post-depositional remanent magnetism: Preliminary experiments with natural sediments. Earth Planet. Sci. Lett. 183, 147–160 (2000).

Article ADS CAS Google Scholar

Roberts, A. P., Tauxe, L. & Heslop, D. Magnetic paleointensity stratigraphy and high-resolution Quaternary geochronology: Successes and future challenges. Quat. Sci. Rev. 61, 1–16 (2013).

Article ADS Google Scholar

Suganuma, Y. et al. 10Be evidence for delayed acquisition of remanent magnetization in marine sediments: Implications for a new age for the Matuyama-Brunhes boundary. Earth Planet. Sci. Lett. 296, 443–450 (2010).

Article ADS CAS Google Scholar

Conrad, T., Heinb, J. R., Paytan, A. & Clague, D. A. Formation of Fe-Mn crusts within a continental margin environment. Ore Geol. Rev. 87, 25–40 (2016).

Article Google Scholar

Soetaert, K. et al. Modeling 210Pb-derived mixing activity in ocean margin sediments: Diffusive versus nonlocal mixing. J. Mar. Res. 54, 1207–1227 (1996).

Article Google Scholar

Download references

This study was supported by the ERC advanced grant 339899 "EDIFICE". We thank S. van de Velde for his constructive review.

Tatiana Savranskaia

Present address: Helmholtz-Zentrum Potsdam Deutsches GeoForschungsZentrum GFZ, Telegrafenberg, 14473, Potsdam, Germany

Institut de Physique Du Globe de Paris, CNRS, Université de Paris, 75005, Paris, France

Tatiana Savranskaia, Ramon Egli & Jean-Pierre Valet

Zentralanstalt für Meteorologie und Geodynamik (ZAMG), 1190, Wien, Austria

Ramon Egli

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

T.S., R. E. and J.-P. V. conceived and planned the experiments. T.S. and R. E. wrote the original draft. T.S. carried out sample preparation and experiments. R.E. designed the mathematical model and performed numerical calculations. J.-P. V. supervised the project. All authors critically revised the manuscript, agree to be fully accountable for ensuring the integrity and accuracy of the work, and read and approved the final manuscript.

Correspondence to Tatiana Savranskaia or Ramon Egli.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Savranskaia, T., Egli, R. & Valet, JP. Multiscale Brazil nut effects in bioturbated sediment. Sci Rep 12, 11450 (2022). https://doi.org/10.1038/s41598-022-14276-w

Download citation

Received: 03 April 2022

Accepted: 03 June 2022

Published: 06 July 2022

DOI: https://doi.org/10.1038/s41598-022-14276-w

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.