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

. It has been hypothesized that the evolution of animals during the Ediacaran–Cambrian transition stimulated the burial of phosphorus in marine sediments. This assump-tion 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 reason-ing 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 (P org and P inorg , re-spectively) 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 P org and P inorg 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.

Abstract. 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 (P org and P inorg , 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 P org and P inorg 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.

Introduction
The evolution of animals near the Ediacaran-Cambrian transition (∼ 542 Myr ago) was a major evolutionary event in Earth's history (Mangano and Buatois, 2017;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 (Aller, 1977;McIlroy and Logan, 1999;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). Biomixing and bio-irrigation can have distinct effects on organic carbon mineralization and early diagenesis (Kostka et al., 2002;van de Velde and Meysman, 2016). For instance, bio-irrigation can promote aerobic respiration by flushing oxygenated bottom waters into deeper anoxic horizons (Archer and Devol, 1992;van de Velde and Meysman, 2016), whereas bio-mixing transports fresh organic from the sediment-water interface into the anoxic zone, thus stimulating anaerobic mineralization pathways (Berner and Westrich, 1985;van de Velde and Meysman, 2016).
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 (P inorg ) 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 (P org ; as for example phosphate esters or phosphonates) or inorganic P minerals during diagenesis, which would then constitute a permanent P org burial sink (Berner et al., 1993;Van Cappellen and Ingall, 1994;Diaz et al., 2008;Goldhammer et al., 2010;Ingall and Jahnke, 1997). 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 Ingall, 1996). If more P becomes buried, photosynthesis would decrease, thereby limiting the overall burial of C org in the seafloor, subsequently increasing atmospheric CO 2 and lowering atmospheric O 2 (Bergman et al., 2004;Berner, 1982). Hence, the rise of bioturbation may have increased atmospheric CO 2 concentrations and decreased atmospheric O 2 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 (Canfield, 1994), causing less mineralization of organic matter and consequently more burial of P org 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 P org than sediments underlying anoxic bottom waters (Ingall and Jahnke, 1997;Ruttenberg, 2014;Slomp and Van Cappellen, 2007), as sediments underlying oxic bottom waters likely receive an extra source of P org 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 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 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.

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 m 2 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 nonbioturbated, and the pore water is rich in hydrogen sulfide 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 nonbioturbated 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 ), indicat- Figure 1. Overview 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.
ing 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 137 Cs 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).

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 N 2 , freeze-dried and stored in a sealed aluminium bag under N 2 atmosphere for later solidphase 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 (Ruttenberg, 1992;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 (Bowman, 1989;Olsen and Sommers, 1982).
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 C org : N tot 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 (P exch ), P associated with iron (P Fe ), authigenic P (P auth ), detrital P (P det ) and P associated with organic matter (P org ). 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 (P inorg ) is calculated as P inorg = P exch + P Fe + P auth + P det (1) and expressed as µmol per gram dry sediment. The extraction procedure used in 2018 separates the sediment phosphorus content into two fractions: total phosphorus (P tot ) and inorganic phosphorus (P inorg ). Organic phosphorus is then calculated as the difference between P tot and P inorg (Bowman, 1989;Olsen and Sommers, 1982). Extractions were performed on a subsample of 1 g, under room temperature and under constantly agitated conditions (extraction procedure detailed in Table A2).

Burial fluxes and solid-phase inventories
Burial fluxes of solid-phase species were calculated based on the sedimentation flux (J sed ) and the concentration of the solid component at the bottom of the sediment column (C solid ).
The sedimentation flux was previously determined based on 210 Pb and 137 Cs 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: where φ x is the porosity at depth x and ρ solidphase is the solidphase 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).

Statistics
The measured values of C org : N tot and C org : P org , calculated burial fluxes of P org and P inorg , 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 C org : P org fulfilled these assumptions. The values of C org : P org were log-transformed before analysis.

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 downcore gradient (Fig. 2a-j), which is expected as sediment bio-mixing by bioturbating fauna acts to erase solid-phase gradients (van de Velde and Meysman, 2016). 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, C org : N tot and C org : P org values are not significantly different between bioturbated and non-bioturbated pond types (p > 0.1; Table A4). Values of C org : N tot 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 (C org : N tot < 10; Burdige, 2006) but are consistent with sediments from temperate salt marshes, which generally have C org : N tot values of 10 or higher (Spivak et al., 2018). Similarly, C org : P org values averaged around 500 in both the non-bioturbated and bioturbated sediment cores (Fig. 2g). These C org : P org values are higher than expected for marine sediments with similar sedimentation rates underlying oxygenated waters (C org : P org of ∼ 200; Slomp and Van Cappellen, 2007) and are more representative of low-oxygen and anoxic marine environments (C org : P org = 300-700; Slomp and Van Cappellen, 2007). However, since the pond waters were oxygenated at the time of sampling (van de Velde et al., 2020a), the elevated C org : P org values most likely reflect the contribution of plant material from the surrounding marsh (C org : P org > 500; Table A3), which is substantially elevated above the C org : P org of marine plankton (∼ 106; Redfield, 1934). Overall, our organic P data show two main findings: (i) bioturbated sediments have similar amounts of P org (P inv = 0.067 ± 0.005 mmol cm −2 ; Fig. 3a) compared with non-bioturbated sediments (P inv = 0.072 ± 0.004 mmol cm −2 ; Figs. 2, 3a), and (ii) differences in P org 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).
SEDEX extractions, performed on sediment cores collected in 2015 and 2016, show that the large difference in P inorg is caused by the much higher P exch and P Fe contents in the bioturbated sediments ( Fig. 4; P exch + P Fe 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 P exch and P Fe at depth seemingly does not lead to more precipitation of P auth (Table A5), since we do not find the detectable formation of authigenic apatite, the concentrations of which are neg-  ligible throughout all sediment cores in both ponds (< 1 µmol g −1 ). While P inorg concentrations are high in the top layers, they strongly decrease with depth; thus P inorg is not efficiently buried. Overall, we find that the burial rate of P inorg 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).

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 ironassociated P. Nevertheless, the accumulation of inorganic P occurs principally in the top layers of the sediment. At 20 cm, the P inorg 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 Fe 2+ 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 (Ruttenberg, 2014). Because vivianite is extracted in the P Fe 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 P inorg , but this effect is not readily seen in our data. As a consequence, we find that the burial of P inorg 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 P inorg minerals is small, consistent with diagenetic modelling results (Dale et al., 2016). We find that the burial of P org 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 (Ruttenberg, 2014), 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 (Canfield, 1994) and would potentially further dilute the effect of bioturbation on P org 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 (Burdige, 2007) 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 P org burial.
Bioturbation is believed to impact P org 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 Devol, 1992;van de Velde and Meysman, 2016), 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 C org : P org 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 longterm 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. Table A2. Specifics 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.

Fraction Extraction solution
Atmosphere Time Analysis SEDEX phosphorus extraction (Ruttenberg, 1992;Slomp et al., 1996) Table A4. Results 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.
One-way ANOVA P inorg burial rate P org burial rate C org : N tot C org : P org  Table A5. Inventories of the individual phosphorus fractions. P org and P inorg 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).