Articles | Volume 18, issue 4
Biogeosciences, 18, 1451–1461, 2021
Biogeosciences, 18, 1451–1461, 2021

Research article 25 Feb 2021

Research article | 25 Feb 2021

Bioturbation has a limited effect on phosphorus burial in salt marsh sediments

Bioturbation has a limited effect on phosphorus burial in salt marsh sediments
Sebastiaan J. van de Velde1,2, Rebecca K. James3, Ine Callebaut4, Silvia Hidalgo-Martinez5, and Filip J. R. Meysman5,6 Sebastiaan J. van de Velde et al.
  • 1Bgeosys, Geoscience, Environment & Society, Université Libre de Bruxelles, Brussels, Belgium
  • 2Operational Directorate Natural Environment, Royal Belgian Institute of Natural Sciences, Brussels, Belgium
  • 3Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, the Netherlands
  • 4Analytical, Environmental & Geo-Chemistry, Vrije Universiteit Brussel, Brussels, Belgium
  • 5Department of Biology, Universiteit Antwerpen, Wilrijk, Belgium
  • 6Department of Biotechnology, Technical University of Delft, Delft, the Netherlands

Correspondence: Sebastiaan J. van de Velde ( and Filip J. R. Meysman (


It has been hypothesized that the evolution of animals during the Ediacaran–Cambrian transition stimulated the burial of phosphorus in marine sediments. This assumption is centrally based on data compilations from marine sediments deposited under oxic and anoxic bottom waters. Since anoxia excludes the presence of infauna and sediment reworking, the observed differences in P burial are assumed to be driven by the presence of bioturbators. This reasoning however ignores the potentially confounding impact of bottom-water oxygenation on phosphorus burial. Here, our goal is to test the idea that bioturbation increases the burial of organic and inorganic phosphorus (Porg and Pinorg, respectively) while accounting for bottom-water oxygenation. We present solid-phase phosphorus speciation data from salt marsh ponds with and without bioturbation (Blakeney salt marsh, Norfolk, UK). In both cases, the pond sediments are exposed to oxygenated bottom waters, and so the only difference is the presence or absence of bioturbating macrofauna. Our data reveal that the rate of Porg and Pinorg burial are indistinguishable between bioturbated and non-bioturbated sediments. A large terrestrial fraction of organic matter and higher sedimentation velocity than generally found in marine sediments (0.3 ± 0.1 cm yr−1) may partially impact these results. However, the absence of a clear effect of bioturbation on total P burial puts into question the presumed importance of bioturbation for phosphorus burial.

1 Introduction

The evolution of animals near the Ediacaran–Cambrian transition ( 542 Myr ago) was a major evolutionary event in Earth’s history (Mangano and Buatois2017; Meysman et al.2006; Wood et al.2019). Early benthic animals developed the ability to burrow (the so-called “burrowing revolution”;  Meysman et al.2006), which profoundly changed the geochemical cycling and burial of elements in the seafloor (Aller1977; McIlroy and Logan1999; Meysman et al.2006). Benthic fauna affect the seafloor in two separate ways: by the reworking of solid-phase particles (bio-mixing) and by the flushing of burrows (bio-irrigation), lumped together under the term “bioturbation” (Kristensen et al.2012). Bio-mixing and bio-irrigation can have distinct effects on organic carbon mineralization and early diagenesis (Kostka et al.2002; van de Velde and Meysman2016). For instance, bio-irrigation can promote aerobic respiration by flushing oxygenated bottom waters into deeper anoxic horizons (Archer and Devol1992; van de Velde and Meysman2016), whereas bio-mixing transports fresh organic from the sediment–water interface into the anoxic zone, thus stimulating anaerobic mineralization pathways (Berner and Westrich1985; van de Velde and Meysman2016).

Bioturbation has been proposed to play a key role in the sedimentary phosphorus (P) cycle. Via bio-mixing, bioturbating organisms transport P that is adsorbed on iron oxides from the oxic zone at the sediment surface into deeper sedimentary layers (Slomp et al.1996). In the anoxic zone of the sediment, these iron oxides are reduced, and P is released further away from the sediment–water interface (SWI). Accordingly, P is retained longer in the sediment, which could stimulate the precipitation and eventual burial of inorganic P (Pinorg) minerals, such as apatite (Slomp et al.1996). Via bio-irrigation, benthic fauna can increase the availability of oxygen in the sediment (Volkenborn et al.2019), which could stimulate the production of microbial polyphosphate compounds within the sediment column (Dale et al.2016). Microbial polyphosphates are generated during the breakdown of organic matter under oxic conditions (Diaz et al.2008) and can be converted into more refractory organic P (Porg; as for example phosphate esters or phosphonates) or inorganic P minerals during diagenesis, which would then constitute a permanent Porg burial sink (Berner et al.1993; Van Cappellen and Ingall1994; Diaz et al.2008; Goldhammer et al.2010; Ingall and Jahnke1997). Overall, bioturbation could increase the burial of organic and inorganic P in the sediment. Accordingly, it has been hypothesized that the rise of bioturbation at the Ediacaran–Cambrian boundary increased the burial of P in marine sediments (Boyle et al.2014; Dale et al.2016; van de Velde et al.2018). Marine P burial plays a key role in the long-term functioning of the Earth system because P is considered the long-term limiting nutrient for marine primary productivity (Van Cappellen and Ingall1996). If more P becomes buried, photosynthesis would decrease, thereby limiting the overall burial of Corg in the seafloor, subsequently increasing atmospheric CO2 and lowering atmospheric O2 (Bergman et al.2004; Berner1982). Hence, the rise of bioturbation may have increased atmospheric CO2 concentrations and decreased atmospheric O2 concentrations, thus inducing warmer climatic conditions and more widespread ocean anoxia (Boyle et al.2014; van de Velde et al.2018).

However, other factors play an important role for the burial of P in marine sediments. High sedimentation velocities stimulate the preservation of organic matter (Canfield1994), causing less mineralization of organic matter and consequently more burial of Porg associated with organic matter. Additionally, bottom-water redox conditions are important for P burial; sediments underlying oxic bottom waters generally show higher burial rates of Porg than sediments underlying anoxic bottom waters (Ingall and Jahnke1997; Ruttenberg2014; Slomp and Van Cappellen2007), as sediments underlying oxic bottom waters likely receive an extra source of Porg in the form of microbial polyphosphates (Diaz et al.2008). It is difficult to disentangle the confounding effect or sediment accumulation, redox conditions and bioturbation because within the present-day seafloor, most sediments deposited under oxic conditions generally also experience bioturbation (Levin et al.1991). As a result, the proposed effect of bioturbation on P burial has not been extensively verified under field or laboratory conditions.

To be able to adequately tease apart the impact of bioturbation from other confounding factors, one requires sites that have oxygenated waters and a similar input of organic matter but no bioturbation. While these conditions are extremely rare in the modern seafloor, they are found in salt marsh ponds along the North Sea coast of Norfolk (UK), which contain sediments with overlying oxygenated water that are either bioturbated or non-bioturbated (Antler et al.2019; Hutchings et al.2019). These two different sediment types can be found in neighbouring ponds, less than 5 m apart, and no systematic difference in local sediment input, organic matter supply or other boundary conditions has been found between the two pond types (Antler et al.2019; Hutchings et al.2019; van de Velde et al.2020a). This remarkable biogeochemical dichotomy between the pond sediments has recently been attributed to alternative stable states, in which small initial differences between ponds are amplified through non-linear positive feedbacks in the sedimentary iron–sulfur cycle (van de Velde et al.2020a).

Whatever the cause of the biogeochemical dichotomy, the important aspect here is that by comparing the geochemistry of the two oxygenated pond types, we can single out the effect of burrowing fauna on sediment biogeochemistry. These ponds within the Norfolk salt marsh complex, hence, provide a unique environment to study the impact of bioturbation on the burial of organic and inorganic phosphorus, without the confounding effect of bottom-water oxygenation. To this end, we collected solid-phase phosphorus data during three separate visits and quantified the burial rates of P in the bioturbated and non-bioturbated ponds.

2 Materials and methods

2.1 Field site

The Blakeney salt marsh (Fig. 1) is part of a larger salt marsh complex along the North Sea coast of East Anglia (UK). The higher, vegetated marsh hosts several shallow, water-filled ponds with a surface area of  50–500 m2 and a water depth of 10–20 cm (Fig. 1;  van de Velde et al.2020a). These ponds show a conspicuous dichotomy in terms of their sediment geochemistry and belong to either one of two end-member types. Pond sediments are either heavily bioturbated, and the solid phase is rich in iron oxides, or sediments are non-bioturbated, and the pore water is rich in hydrogen sulfide (Antler et al.2019; Hutchings et al.2019; van de Velde et al.2020a). The bioturbated ponds are colonized by large macrofauna, mostly Nereis and Arenicola, at high densities ( 1000 organisms m−2; Antler et al.2019), while the non-bioturbated sediments do not show signs of macrofauna or burrows. Water column concentrations of oxygen, dissolved inorganic carbon (DIC) and nutrients are not statistically different between pond types, and both types of ponds have the same sedimentation flux (0.9 ± 0.1 kg m−2 yr−1), indicating they receive a similar input of detrital minerals. Furthermore, all ponds display similar inputs of organic matter, thus suggesting that the only important difference between the ponds is the presence of burrowing fauna (van de Velde et al.2020a). Depth profiles of 137Cs show a well-defined peak in the non-bioturbated ponds, suggesting they have been undisturbed by fauna for at least 60 years (van de Velde et al.2020a).

Figure 1Overview of the field site in the Blakeney salt marsh system, UK. Sampled ponds are outlined in white (non-bioturbated) and blue (bioturbated), with ponds where sediment cores taken denoted with the sampling year. Coordinates of the sampled ponds are given in Table A1. (b, c) Map data © Google Maps 2020.

2.2 Sediment sampling and analysis

Sediment cores were collected on three separate visits (October 2015, August 2016 and August 2018). A total of 12 ponds were examined in total (sampling sites are indicated in Fig. 1). During each sampling campaign, two replicate sediment cores were collected from each sampled pond. Core sectioning was done at 0.5 cm resolution from 0 to 3 cm depth, at 1 cm resolution between 3 and 8 cm depth, and in 2 cm slices from 8 to 22 cm depth. Sediment sections were collected in 50 mL centrifuge tubes (polypropylene; TPP Techno Plastic Products AG, Switzerland). In 2015 and 2016, sediment cores were processed under anaerobic conditions in a glove bag with N2, freeze-dried and stored in a sealed aluminium bag under N2 atmosphere for later solid-phase analysis. In 2018, cores were immediately sliced in open air in the field and subsequently freeze-dried and stored under room conditions in 50 mL centrifuge tubes. This difference in sampling procedure reflects the subsequent analyses. In 2015 and 2016, the solid phase was analysed for different phosphorus fractions following the SEDEX (sequential extraction method) phosphorus extraction (Ruttenberg1992; Slomp et al.1996). In 2018, due to the higher number of samples, we decided to use a simpler and faster extraction method that only differentiates between inorganic and organic phosphorus fractions (Bowman1989; Olsen and Sommers1982).

Sediment samples from all three campaigns were analysed for carbon and nitrogen content. To this end, freeze-dried solid-phase samples were ground to a fine powder and analysed by an Interscience Flash 2000 organic element analyser (precision < 5 %) for determination of particulate organic carbon (POC) and total nitrogen (TN). Before analysis, samples for POC were first acidified with 0.1M HCl to remove the inorganic carbon (Nieuwenhuize et al.1994). Concentrations of POC and TN are expressed as mass % of dry sediment. Organic matter Corg:Ntot was calculated as the molar ratio of POC over TN. The POC and TN results have been presented previously in van de Velde et al. (2020a).

The SEDEX procedure used in 2015 and 2016 separates total sedimentary P into five fractions: exchangeable P (Pexch), P associated with iron (PFe), authigenic P (Pauth), detrital P (Pdet) and P associated with organic matter (Porg). All extractions were performed on a subsample of 300 mg, under room temperature and under constantly agitated conditions (the extraction procedure is detailed in Table A2). Inorganic phosphorus (Pinorg) is calculated as

(1) P inorg = P exch + P Fe + P auth + P det

and expressed as µmol per gram dry sediment. The extraction procedure used in 2018 separates the sediment phosphorus content into two fractions: total phosphorus (Ptot) and inorganic phosphorus (Pinorg). Organic phosphorus is then calculated as the difference between Ptot and Pinorg (Bowman1989; Olsen and Sommers1982). Extractions were performed on a subsample of 1 g, under room temperature and under constantly agitated conditions (extraction procedure detailed in Table A2).

2.3 Burial fluxes and solid-phase inventories

Burial fluxes of solid-phase species were calculated based on the sedimentation flux (Jsed) and the concentration of the solid component at the bottom of the sediment column (Csolid).

(2) J burial = J sed C solid

The sedimentation flux was previously determined based on 210Pb and 137Cs dating and was statistically indistinguishable between pond types (0.9 ± 0.1 kg m−2 yr−1; van de Velde et al.2020a).

Solid-phase inventories were calculated by integrating measured concentrations over the first 20 cm of the sediment cores:

(3) INV = ρ solidphase x down x up ( 1 - ϕ x ) C solid d x ,

where ϕx is the porosity at depth x and ρsolidphase is the solid-phase density (previously determined to be 2.2 g cm−3; van de Velde et al.2020a). The porosity depth profile was determined from the water content and solid-phase density, considering the salt content of the pore water. The water content of the sediment was determined as the difference in sediment weight before and after freeze-drying. Porosity profiles have been presented previously in van de Velde et al. (2020a).

2.4 Statistics

The measured values of Corg:Ntot and Corg:Porg, calculated burial fluxes of Porg and Pinorg, and solid-phase inventories were averaged over the duplicate cores from each pond, and one-way ANOVA tests were used to test for significant differences between bioturbated and non-bioturbated ponds. Due to the different methods used for P extraction in 2015 and 2016 relative to 2018, the inclusion of the sampling year as a random effect was tested through model evaluation using AIC (Akaike information criterion). Sampling year was not significant for all dependent factors, and so it was excluded from subsequent analyses. Residuals were tested for normality and homoscedasticity, and all but Corg:Porg fulfilled these assumptions. The values of Corg:Porg were log-transformed before analysis.

3 Results

3.1 Bioturbation and organic phosphorus burial

Individual depth profiles of solid-phase variables (particulate organic carbon, POC; total nitrogen, TN; total phosphorus, TP; particulate organic phosphorus, POP; and particulate inorganic phosphorus, PIP) show variation within the sediments of a given pond type, indicating spatial heterogeneity (Fig. 2). Still, averaged depth profiles of the six non-bioturbated and the six bioturbated ponds reveal that non-bioturbated sediments have generally a stronger down-core gradient (Fig. 2a–j), which is expected as sediment bio-mixing by bioturbating fauna acts to erase solid-phase gradients (van de Velde and Meysman2016). Differences in POC, TN and POP inventories in the top 20 cm between bioturbated and non-bioturbated cores are insignificant (Fig. 2a, d, g, h; Table A4). Consistently, Corg:Ntot and Corg:Porg values are not significantly different between bioturbated and non-bioturbated pond types (p>0.1; Table A4). Values of Corg:Ntot range from 9 to 15, with an average value of 12 (Fig. 2k). These values are slightly higher than expected for sediments from fully marine settings (Corg:Ntot< 10; Burdige2006) but are consistent with sediments from temperate salt marshes, which generally have Corg:Ntot values of 10 or higher (Spivak et al.2018). Similarly, Corg:Porg values averaged around 500 in both the non-bioturbated and bioturbated sediment cores (Fig. 2g). These Corg:Porg values are higher than expected for marine sediments with similar sedimentation rates underlying oxygenated waters (Corg:Porg of ∼200; Slomp and Van Cappellen2007) and are more representative of low-oxygen and anoxic marine environments (Corg:Porg= 300–700; Slomp and Van Cappellen2007). However, since the pond waters were oxygenated at the time of sampling (van de Velde et al.2020a), the elevated Corg:Porg values most likely reflect the contribution of plant material from the surrounding marsh (Corg:Porg> 500; Table A3), which is substantially elevated above the Corg:Porg of marine plankton ( 106; Redfield1934).

Figure 2Vertical solid-phase profiles of (a, b) particulate organic carbon (POC), (c, d) total nitrogen (TN), (e, f) total phosphorus (TP), (g, h) particulate organic phosphorus (POP) and (i, j) particulate inorganic phosphorus (PIP). Individual cores are plotted as a light-grey lines, and the average vertical profile over all cores is plotted as black dots. Values are averaged inventories and errors are 1 standard deviation. (k) Boxplots of the Corg:Ntot and Corg:Porg of the particulate organic matter fraction.


Overall, our organic P data show two main findings: (i) bioturbated sediments have similar amounts of Porg (Pinv=0.067± 0.005 mmol cm−2; Fig. 3a) compared with non-bioturbated sediments (Pinv=0.072± 0.004 mmol cm−2; Figs. 23a), and (ii) differences in Porg burial in bioturbated sediments (17 ± 4 µmol m−2 d−1) are not statistically significant from non-bioturbated sediments (19 ± 2 µmol m−2 d−1) (Fig. 3b; p>0.1; Table A4).

Figure 3Boxplots of (a) inventories and (b) burial rates of inorganic and organic P fractions. The data of 2015 are not included in (a), since some depth samples were not analysed for P fractionation.


3.2 Bioturbation and inorganic phosphorus burial

Our solid-phase analyses show that bioturbated sediments (203 ± 70 µmol cm−2) contain significantly more particulate inorganic phosphorus than non-bioturbated sediments (113 ± 17 µmol cm−2) (p<0.1; Figs. 2i, j, 3a; Table A4).

SEDEX extractions, performed on sediment cores collected in 2015 and 2016, show that the large difference in Pinorg is caused by the much higher Pexch and PFe contents in the bioturbated sediments (Fig. 4; Pexch+ PFe is up to 90 mmol P cm−2 higher in the bioturbated sediments; Table A5). Phosphorus adsorbed onto iron oxide minerals is transported from the sediment–water interface to the deeper sedimentary layers by the downward mixing of benthic fauna (Slomp et al.1996). Indeed,  50% of the Fe minerals in the bioturbated sediments at the Blakeney salt marsh are in oxidized form, whereas < 10 % are in the non-bioturbated sediments (see van de Velde et al.2020a, for an extended discussion of the Fe–S cycle at the field site).

The increased inventory of Pexch and PFe at depth seemingly does not lead to more precipitation of Pauth (Table A5), since we do not find the detectable formation of authigenic apatite, the concentrations of which are negligible throughout all sediment cores in both ponds (<1µmol g−1). While Pinorg concentrations are high in the top layers, they strongly decrease with depth; thus Pinorg is not efficiently buried. Overall, we find that the burial rate of Pinorg is 38 ± 3 µmol m−2 d−1 in bioturbated sediments and 36 ± 3 µmol m−2 d−1 in the non-bioturbated sediment, and this difference is not significant (p>0.1; Table A4).

Figure 4Vertical solid-phase profiles of non-bioturbated (a) and bioturbated (f) particulate inorganic phosphorus (PIP). Values are depth-averaged concentrations. (b–e, g–j) Solid-phase phosphorus speciation (SEDEX) of cores collected in 2015 and 2016. Pexch: exchangeable phosphorus, PFe: iron-bound phosphorus, Pauth: authigenic phosphorus, Pdet: detrital phosphorus, Porg: organic phosphorus. Note that the Pauth fraction is too small to be visible in the figure.


4 Discussion

By comparing non-bioturbated and bioturbated sediment cores collected from an East Anglian salt marsh, this study provides a field verification for the hypothesis that bioturbation stimulates the burial of inorganic and organic P minerals. Our field data show that a bioturbated sediment contains more inorganic P, which exists mainly in the form of iron-associated P. Nevertheless, the accumulation of inorganic P occurs principally in the top layers of the sediment. At 20 cm, the Pinorg levels in bioturbated sediments decrease to similar levels as in non-bioturbated sediments (Fig. 4). We find no evidence for the formation of authigenic P (Fig. 4), which is at odds with previous hypotheses (Slomp et al.1996; Zhao et al.2020). Potentially, the high concentrations of Fe2+ in the pore water of the bioturbated sediments (up to 300 µM; van de Velde et al.2020a) promoted the formation of vivianite over apatite (Ruttenberg2014). Because vivianite is extracted in the PFe fraction (Nembrini et al.1983), we cannot separate between P associated with iron oxides and P in the form of vivianite. The precipitation of vivianite in the bioturbated sediments could potentially stimulate the burial of Pinorg, but this effect is not readily seen in our data. As a consequence, we find that the burial of Pinorg was not significantly different between bioturbated and non-bioturbated ponds, which does not support the idea that bioturbation stimulates the burial of inorganic P minerals. Our data are hence not congruent with previous modelling and field studies that have suggested that bioturbating fauna stimulates the burial of authigenic apatite (Slomp et al.1996; Zhao et al.2020). These studies argued that bioturbators mix iron oxides, on which phosphate is adsorbed, further away from the sediment–water interface (SWI), increasing the retention time of phosphate in the sediment and subsequently stimulating the precipitation and eventual burial of inorganic phosphorus minerals (Slomp et al.1996; Dale et al.2016). Bioturbators also flush their burrows however, which removes phosphate from the pore water (Dale et al.2016). Consequently, depending on the animal community (irrigation versus mixing), the pore water can become undersaturated with respect to apatite, leading to lower precipitation rates and burial of apatite. The difference between oxic bioturbated and oxic non-bioturbated sediments inferred from the model study of Dale et al. (2016) is only a few µmol m−2 d−1, which is likely undetectable in a field study like ours (Fig. 3). Consequently, our results suggest that the effect of bioturbation on Pinorg minerals is small, consistent with diagenetic modelling results (Dale et al.2016).

We find that the burial of Porg is indistinguishable between bioturbated and non-bioturbated sediments (Fig. 2g). Accordingly, our field data do not support the idea that bioturbation increases the burial of organic phosphorus. It should be noted however that the C:N:P ratios of the organic matter in the salt marsh pond sediments suggest a large contribution of plant material from the surrounding marsh (Fig. 2), which has elevated C:N:P ratios compared to marine organic matter (Table A3). Because terrestrial material is less easily degraded, the addition of plant material to the sediment organic matter pool dilutes the signal of the marine organic matter (Ruttenberg2014), which could explain why we see no effect of bioturbation. Another possible reason is that the sedimentation velocity in Norfolk salt marsh ponds is relatively high (0.3 cm yr−1, compared to a mean value of 0.1 cm yr−1 for shelf sediments; Burwicz et al.2011). High sedimentation velocities stimulate the preservation of organic matter (Canfield1994) and would potentially further dilute the effect of bioturbation on Porg burial. Nevertheless, the POC profiles suggest an appreciable amount of organic matter is being degraded in the salt marsh pond sediments (Fig. 2). Indeed, organic matter mineralization rates based on POC and nutrient profiles range from 6 to 38 mmol C m−2 d−1, which is comparable to rates observed in shallow marine sediments (Burdige2007) and thus indicates that the Blakeney salt marsh site shows a comparable metabolic activity to marine sediments. The observation then still remains; bioturbation seems to have no significant impact on Porg burial.

Bioturbation is believed to impact Porg burial by stimulating the formation of polyphosphates, which can act as an intermediate for the formation of apatite or more refractory organic P compounds (Berner et al.1993; Goldhammer et al.2010). Polyphosphates are formed during aerobic respiration of organic matter (Diaz et al.2008), and because bio-irrigation stimulates aerobic respiration in the sediments (Archer and Devol1992; van de Velde and Meysman2016), it has been hypothesized that more polyphosphates could be formed in bioturbated sediments (Dale et al.2016). However, other studies have shown that polyphosphates and the subsequent conversion to apatite also occur under anoxic conditions (Goldhammer et al.2010), so it is questionable whether periodic influshing of oxygenated water would have a large impact on polyphosphate formation. Unfortunately, our analysis method does not allow us to differentiate between bulk organic P and microbial polyphosphates to investigate whether polyphosphate formation plays a role in both bioturbated and non-bioturbated sediments. This would be an interesting avenue for future research.

A number of caveats need to be taken into account when interpreting our results, as observations made in the surface sediments of the Norfolk salt marsh ponds are not necessarily representative of burial signals made several metres deep in coastal and shelf sediments. As noted above, the Corg:Porg ratio of buried organic matter is high, as it is influenced by the deposition of the salt marsh vegetation with high C:P ratios, and therefore, caution is required when extrapolating our results to other coastal and shelf sediments. A difference in timescale is another caveat to take into account. Based on an average sedimentation velocity of 0.3 cm yr−1 measured via radionuclide dating (van de Velde et al.2020a), our samples represent only 60 years of sediment accumulation, and they hence indicate that, after ∼60 years of early diagenesis, there is no discernible difference in burial of P minerals. Diagenesis later in the burial history could still affect long-term P burial. Data reported in the literature are collected from geological formations (see e.g. Ingall et al.1993) or deep drill cores (see e.g. Slomp et al.2004). These samples have undergone several hundreds to thousands of years of diagenesis before being sampled and analysed. However, we cannot see a mechanism via which bioturbation could have an imprint on the sediment that eventually leads to differential diagenesis of P minerals later in their burial history. If bioturbation does not already create a difference in P burial in the early phase of diagenesis, it is highly unlikely that will happen afterwards. So the difference in timescale cannot easily explain why P burial rates are similar in bioturbated and non-bioturbated sediments.

Overall, more field studies and laboratory experiments are required to verify our results. Nevertheless, the absence of a clear impact of bioturbation on P burial questions whether the modelled effects of the evolution of bioturbation in the early Cambrian, if they are solely based on inferred effects of bioturbation on P burial, do not overestimate the true impact of the burrowing revolution. Our study highlights that assumptions made in Earth system models should be firmly grounded in field observations, before accurate inferences about large-scale questions can be made, such as the impact of the burrowing revolution on the composition of the atmosphere and climate.

Appendix A

Table A1Coordinates and type of the ponds sampled in the 2015, 2016 and 2018 field campaigns in the Blakeney salt marsh. See Fig. 1c in the main text for relative geographical location of the ponds.

Download Print Version | Download XLSX

(Ruttenberg1992; Slomp et al.1996)(Bowman1989; Olsen and Sommers1982)

Table A2Specifics for the sequential extractions of phosphorus. After the extraction, the sample was centrifuged (2500 g for 10 min), and the supernatant was filtered (0.45 µm cellulose acetate). ICP-OES: inductively coupled plasma optical emission spectrometry.

Download Print Version | Download XLSX

Table A3Summary of salt marsh plant properties. Data were first published in van de Velde et al. (2020a).

Download Print Version | Download XLSX

Table A4Results from one-way ANOVA analysis and calculated 95 % confidence intervals. Residuals were tested for normality and heterogeneity; all passed these assumptions except C:P, which was log-transformed. df: degrees of freedom.

Download Print Version | Download XLSX

Table A5Inventories of the individual phosphorus fractions. Porg and Pinorg are reported as mean ± 1 SD; SEDEX fractions are given as a range (due to the low number of replicates) from the 2016 cores (the 2015 cores were omitted because some depth samples were not analysed for P fractionation).

Download Print Version | Download XLSX

Data availability

All data presented in this paper are available from the VLIZ data repository (, van de Velde et al.2020b).

Author contributions

SJV conceived the hypothesis. SJvdV and FJRM organized the field sampling. All authors contributed to the field sampling. SJvdV and IC performed the SEDEX extractions. SJvdV analysed the data. RKJ did the statistical analyses. SJvdV and FJRM wrote the paper with input from all co-authors.

Competing interests

The authors declare that they have no conflict of interest.


The authors would like to thank Tom Van der Spriet from the University of Antwerp for the phosphorus extraction of the 2018 sediment samples.

Financial support

This research has been supported by the Fonds Wetenschappelijk Onderzoek (grant no. G038819N), the Universiteit Antwerpen (grant no. TOPBOF), the Netherlands Organisation for Scientific Research (grant no. 016.VICI.170.072) and the Belgian Federal Science Policy Office (grant no. FED-tWIN2019-prf-008).

Review statement

This paper was edited by Marilaure Grégoire and reviewed by Peter Kraal and one anonymous referee.


Aller, R. C.: The influence of macrobenthos on chemical diagenesis of marine sediments, PhD thesis, Yale University, New Haven, Connecticut, 600 pp., 1977. a

Antler, G., Mills, J. V., Huthings, A., Redeker, K., and Turchyn, A. V.: The sedimentary carbon-sulfur-iron interplay – a lesson from East Anglian salt marsh sediments, Front. Earth Sci., 7, 140,, 2019. a, b, c, d

Archer, D. and Devol, A. H.: Benthic oxygen fluxes on the Washington shelf and slope: A comparison of in situ microelectrode and chamber flux measurements, Limnol. Oceanogr., 37, 614–629,, 1992. a, b

Bergman, N. M., Lenton, T. M., and Watson, A. J.: COPSE: A new model of biogeochemical cycling over phanerozoic time, Am. J. Sci., 304, 397–437,, 2004. a

Berner, R. A.: Burial of Organic Carbon and Pyrite sulfur in the modern ocean: Its geochemical and environmental sinificance, Am. J. Sci., 282, 451–473, 1982. a

Berner, R. A. and Westrich, J. T.: Bioturbation and the early diagenesis of carbon and sulfur, Am. J. Sci., 285, 193–206,, 1985 a

Berner, R. A., Ruttenberg, K. C., Ingall, E. D., and Rao, J.-L.: The Nature of Phosphorus Burial in Modern Marine Sediments, in: Interactions of C, N, P and S Biogeochemical Cycles and Global Change, edited by: Wollast, R., Mackenzie, F. T., and Chou, L., Springer, Berlin and Heidelberg, Germany, 365–378, 1993. a, b

Bowman, R. A.: A sequential extraction procedure with concentrated sulfuric acid and dilute base for soil organic phosphorus, Soil Sci. Soc. Am. J., 53, 362–366,, 1989. a, b, c

Boyle, R. A., Dahl, T. W., Dale, A. W., Zhu, M., Brasier, M. D., Canfield, D. E., and Lenton, T. M.: Stabilization of the coupled oxygen and phosphorus cycles by the evolution of bioturbation, Nat. Geosci., 7, 671–676,, 2014. a, b

Burdige, D. J.: Geochemistry of Marine Sediments, Princeton University Press, Princeton , USA, 2006. a

Burdige, D. J.: Preservation of Organic Matter in Marine Sediments: Controls, Mechanisms, and an Imbalance in Sediment Organic Carbon Budgets?, Chem. Rev., 107, 467–485,, 2007. a

Burwicz, E. B., Rüpke, L. H., and Wallmann, K.: Estimation of the global amount of submarine gas hydrates formed via microbial methane formation based on numerical reaction-transport modeling and a novel parameterization of Holocene sedimentation, Geochim. Cosmochim. Ac., 75, 4562–4576,, 2011. a

Canfield, D. E.: Factors influencing organic carbon preservation in marine sediments, Chem. Geol., 114, 315–329, 1994. a, b

Dale, A. W., Boyle, R. A., Lenton, T. M., Ingall, E. D., and Wallmann, K.: A model for microbial phosphorus cycling in bioturbated marine sediments: Significance for phosphorus burial in the early Paleozoic, Geochim. Cosmochim. Ac., 189, 251–268,, 2016. a, b, c, d, e, f, g

Diaz, J., Ingall, E., Benitez-Nelson, C., Paterson, D., De Jonge, M. D., Mcnulty, I., and Brandes, J. A.: Marine Polyphosphate: A Key Player in Geologic Phosphorus Sequestration, Science, 320, 652–656,, 2008. a, b, c, d

Goldhammer, T., Brüchert, V., Ferdelman, T. G., and Zabel, M.: Microbial sequestration of phosphorus in anoxic upwelling sediments, Nat. Geosci., 3, 557–561,, 2010. a, b, c

Hutchings, A. M., Antler, G., Wilkening, J., Basu, A., Bradbury, H. J., Clegg, J. A., Gorka, M., Lin, C. Y., Mills, J. V., Pellerin, A., Redeker, K., Sun, X., and Turchyn, A. V.: Creek dynamics determine pond subsurface geochemical heterogeneity in East Anglian (UK) salt marshes, Front. Earth Sci., 7, 41,, 2019. a, b, c

Ingall, E. and Jahnke, R.: Influence of water-column anoxia on the elemental fractionation of carbon and phosphorus during sediment diagenesis, Mar. Geol., 139, 219–229,, 1997. a, b

Ingall, E. D., Bustin, R. M., and Van Cappellen, P.: Influence of water column anoxia on the burial and preservation of carbon and phosphorus in marine shales, Geochim. Cosmochim. Ac., 57, 303–316, 1993. a

Kostka, J. E., Gribsholt, B., Petrie, E., Dalton, D., Skelton, H., and Kristensen, E.: The rates and pathways of carbon oxidation in bioturbated saltmarsh sediments, Limnol. Oceanogr., 47, 230–240,, 2002. a

Kristensen, E., Penha-Lopes, G., Delefosse, M., Valdemarsen, T., Quintana, C. O., and Banta, G. T.: What is bioturbation? the need for a precise definition for fauna in aquatic sciences, Mar. Ecol. Prog. Ser., 446, 285–302,, 2012. a

Levin, L. A., Huggett, C. L., and Wishner, K. F.: Control of deep-sea benthic community structure by oxygen and organic-matter gradients in the eastern Pacific Ocean, J. Mar. Res., 49, 763–800,, 1991. a

Mángano, G. M. and Buatois, L. A.: The Cambrian revolutions: Trace-fossil record, timing, links and geobiological impact, Earth-Science Rev., 173, 96–108,, 2017. a

McIlroy, D. and Logan, G. A.: The impact of bioturbation on infaunal ecology and evolution during the Proterozoic-Cambrian transition, Palaios, 14, 58–72,, 1999. a

Meysman, F. J. R., Middelburg, J. J., and Heip, C. H. R.: Bioturbation: a fresh look at Darwin’s last idea., Trends Ecol. Evol., 21, 688–95,, 2006. a, b, c

Nembrini, G. P., Capobianco, J. A., Viel, M., and Williams, A. F.: A Mössbauer and chemical study of the formation of vivianite in sediments of Lago Maggiore (Italy), Geochim. Cosmochim. Ac., 47, 1459–1464,, 1983. a

Nieuwenhuize, J., Maas, Y. E. M., and Middelburg, J. J.: Rapid analysis of organic carbon and nitrogen in particulate materials, Mar. Chem., 45, 217–224, 1994. a

Olsen, S. R. and Sommers, L. E.: Phosphorus, in: Methods of Soil Analysis, Part 2, edited by: Page, A. L., Miller, R. H., and Keeney, D. R., American Society of Agronomy, Inc. Madison., American Society of Agronomy, Soil Science Society of America, Madison, USA, 403–430, 1982. a, b, c

Redfield, A. C.: On the proportions of organic derivations in sea water and their relation to the composition of plankton, in: James Johnstone Memorial Volume, edited by: Daniel, R. J., University Press of Liverpool, Liverpool, UK, 177–192, 1934. a

Ruttenberg, K. C.: Development of a sequential extraction method for different forms of phosphorus in marine sediments, Limnol. Oceanogr., 37, 1460–1482, 1992. a, b

Ruttenberg, K. C.: The Global Phosphorus Cycle, 2nd Edn., Elsevier, 499–558,, 2014. a, b, c

Slomp, C. P. and Van Cappellen, P.: The global marine phosphorus cycle: sensitivity to oceanic circulation, Biogeosciences, 4, 155–171,, 2007. a, b, c

Slomp, C. P., Epping, E. H. G., Helder, W., and Van Raaphorst, W.: A key role for iron-bound phosphorus in authigenic apatite formation in North Atlantic continental platform sediments, J. Mar. Res., 54, 1179–1205, 1996a.  a, b, c, d, e, f, g, h

Slomp, C. P., Thomson, J., and De Lange, G. J.: Controls on phosphorus regeneration and burial during formation of eastern Mediterranean sapropels, Mar. Geol., 203, 141–159,, 2004. a

Spivak, A. C., Gosselin, K. M., and Sylva, S. P.: Shallow ponds are biogeochemically distinct habitats in salt marsh ecosystems, Limnol. Oceanogr., 63, 1622–1642,, 2018. a

Van Cappellen, P. and Ingall, E. D.: Benthic phosphorus regeneration, net primary production, and ocean anoxia: A model of the coupled marine biogeochemical cycles of carbon and phosphorus, Paleoceanography, 9, 677–692, 1994. a

Van Cappellen, P. and Ingall, E. D.: Redox stabilization of the Atmosphere and Oceans by Phosphorus-Limited Marine Productivity, Science, 271, 493–496,, 1996. a

van de Velde, S. and Meysman, F. J. R.: The influence of bioturbation on iron and sulphur cycling in marine sediments: a model analysis, Aquat. Geochem., 22, 469–504,, 2016. a, b, c, d, e

van de Velde, S., Mills, B., Meysman, F. J., Lenton, T. M., and Poulton, S. W.: Early Palaeozoic ocean anoxia and global warming driven by the evolution of shallow burrowing, Nat. Commun., 9, 2554,, 2018. a, b

van de Velde, S., Hidalgo-Martinez, S., Callebaut, I., Antler, G., James, R., Leermakers, M., and Meysman, F.: Burrowing fauna mediate alternative stable states in the redox cycling of salt marsh sediments, Geochim. Cosmochim. Ac., 276, 31–49,, 2020a. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o

van de Velde, S., James, R. K., Callebaut, I., Hidalgo-Martinez, S., and Meysman, F. J. R.: Sedimentary and plant carbon, nitrogen and phosphorus measurements of Blakeney salt marsh ponds, Marine Data Archive,, 2020b. a

Volkenborn, N., Woodin, S. A., Wethey, D. S., and Polerecky, L.: Bioirrigation, in: Encyclopedia of Ocean Sciences, edited by: Cochran, J. K., Bokuniewicz, H. J., and Yager, P. L., Academic Press, 663–670,, 2019. a

Wood, R., Liu, A. G., Bowyer, F., Wilby, P. R., Dunn, F. S., Kenchington, C. G., Cuthill, J. F. H., Mitchell, E. G., and Penny, A.: Integrated records of environmental change and evolution challenge the Cambrian Explosion, Nat. Ecol. Evol., 3, 528–538,, 2019. a

Zhao, M., Zhang, S., Tarhan, L. G., Reinhard, C. T., and Planavsky, N.: The role of calcium in regulating marine phosphorus burial and atmospheric oxygenation, Nat. Commun., 11, 2232,, 2020. a, b

Short summary
Some 540 Myr ago, animal life evolved in the ocean. Previous research suggested that when these early animals started inhabiting the seafloor, they retained phosphorus in the seafloor, thereby limiting photosynthesis in the ocean. We studied salt marsh sediments with and without animals and found that their impact on phosphorus retention is limited, which implies that their impact on the global environment might have been less drastic than previously assumed.
Final-revised paper