Description: Mud Volcanism and fluid seepage are widespread phenomena in the Gulf of Cadiz (SW Iberian Margin). In this seismically active region located at the boundary between the African and Eurasian plates, fluid flow is typically focused on deeply rooted active strike-slip faults. The geochemical signature of emanating fluids from various mud volcanoes (MVs) has been interpreted as being largely affected by clay mineral dehydration and recrystallization of Upper Jurassic carbonates. Here we present the results of a novel, fully-coupled 1D basin-scale reactive-transport model capable of simulating major fluid forming processes and related geochemical signatures by considering the growth of the sediment column over time, compaction of sediments, diffusion and advection of fluids, as well as convective and conductive heat flow. The outcome of the model is a realistic approximation to the development of the sediment pore water system over geological time scales in the Gulf of Cadiz. Combined with a geochemical reaction transport model for clay mineral dehydration and calcium carbonate recrystallization, we were able to reproduce measured concentrations of Cl, strontium and 87Sr/86Sr of emanating mud volcano fluids. These results support previously made qualitative interpretations and add further constraints on fluid forming processes, reaction rates and source depths. The geochemical signature at Porto MV posed a specific problem, because of insufficient constraints on non-radiogenic 87Sr/86Sr sources at this location. We favour a scenario of basement-derived fluid injection into basal Upper Jurassic carbonate deposits (Hensen et al., 2015). Although the mechanism behind such basement-derived flow, e.g. along permeable faults, remains speculative at this stage, it provides an additional source of low 87Sr/86Sr fluids and offers an idea on how formation water from the deepest sedimentary strata above the basement can be mobilized and eventually initiate the advection of fluids feeding MVs at the seafloor. The dynamic reactive-transport model presented in this study provides a new tool addressing the combined simulation of complex physical-geochemical processes in sedimentary systems. The model can easily be extended and applied to similar geological settings, and thus help us to provide a fundamental understanding of fluid dynamics and element recycling in sedimentary basins.

Description: The kinetics of particulate organic carbon (POC) mineralization in marine surface sediments is not well constrained. This creates considerable uncertainties when benthic processes are considered in global biogeochemical or Earth system circulation models to simulate climate-ocean interactions and biogeochemical tracers in the ocean. In an attempt to improve our understanding of the rate and depth distribution of organic carbon mineralization in bioturbated (0–10 cm) sediments, we parameterized a 1-D diagenetic model that simulates the reactivity of three discrete POC pools at global scale (a "multi-G" model). The rate constants of the three reactive classes (highly reactive, reactive, refractory) are fixed and determined to be 70 yr−1, 0.5 yr−1, and ~0.001 yr−1, respectively, based on the Martin curve model for pelagic POC degradation. In contrast to previous approaches, the reactivity of the organic material degraded in the seafloor is continuous with, and set by, the apparent reactivity of material sinking through the water column. The model is able to simulate a global database (185 stations) of benthic oxygen and nitrate fluxes across the sediment-water interface in addition to porewater oxygen and nitrate distributions and organic carbon burial efficiencies. It is further consistent with degradation experiments of fresh phytoplankton. We propose that an important yet mostly overlooked consideration in previous upscaling approaches is the proportion of the relative reactive POC classes reaching the seafloor in addition to their reactivity. The approach presented is applicable to both steady-state and non-steady state scenarios, and links POC degradation kinetics in sedimentary environments to water depth and the POC rain rate to the seafloor.

Description: Iron speciation was determined in hemiplegic sediments from a high productivity area to investigate systematically the early diagenetic reactivity of Fe. A combination of various leaching agents (1 M HCI, dithionite buffered in citrate/acetic acid, HF/H2SO4, acetic Cr(II)) was applied to sediment and extracted more than 80% of total Fe. Subsequent Fe species determination defined specific mineral fractions that are available for Fe reduction and fractions formed as products of Fe diagenesis. To determine the Fe speciation of (sheet) silicates we explored an extraction procedure (HF/H2SO4) and verified the procedure by application to standard rocks. Variations of Fe speciation of (sheet) silicates reflect the possible formation of Fe-bearing silicates in near surface sediments. The same fraction indicates a change in the primary input at greater depth, which is supported by other parameters. The Fe(II)/ Fe(III) -ratio of total sediment determined by extractions was compared with Mössbauer-spectroscopy ] at room temperature and showed agreement within 10%. M6ssbauer-spectroscopy indicates the occurrence of siderite in the presence of free sulfide and pyrite, supporting the importance of microenvironments during mineral formation. The occurrence of other Fe(II) bearing minerals such as ankerite (Ca-, Fe-, Mg-carbonate) can be presumed but remains speculative.

Description: Oxygen minimum zones (OMZs) that impinge on continental margins favor the release of phosphorus (P) from the sediments to the water column, enhancing primary productivity and the maintenance or expansion of low-oxygen waters. A comprehensive field program in the Peruvian OMZ was undertaken to identify the sources of benthic P at six stations, including the analysis of particles from the water column, surface sediments, and pore fluids, as well as in situ benthic flux measurements. A major fraction of solid-phase P was bound as particulate inorganic P (PIP) both in the water column and in sediments. Sedimentary PIP increased with depth in the sediment at the expense of particulate organic P (POP). The ratio of particulate organic carbon (POC) to POP exceeded the Redfield ratio both in the water column (202 ± 29) and in surface sediments (303 ± 77). However, the POC to total particulate P (TPP = POP + PIP) ratio was close to Redfield in the water column (103 ± 9) and in sediment samples (102 ± 15). This suggests that the relative burial efficiencies of POC and TPP are similar under low-oxygen conditions and that the sediments underlying the anoxic waters on the Peru margin are not depleted in P compared to Redfield. Benthic fluxes of dissolved P were extremely high (up to 1.04 ± 0.31 mmol m−2 d−1), however, showing that a lack of oxygen promotes the intensified release of dissolved P from sediments, whilst preserving the POC / TPP burial ratio. Benthic dissolved P fluxes were always higher than the TPP rain rate to the seabed, which is proposed to be caused by transient P release by bacterial mats that had stored P during previous periods when bottom waters were less reducing. At one station located at the lower rim of the OMZ, dissolved P was taken up by the sediments, indicating ongoing phosphorite formation. This is further supported by decreasing porewater phosphate concentrations with sediment depth, whereas solid-phase P concentrations were comparatively high.

Description: Highlights • An empirical diagenetic model includes redox-dependent P storage by microorganisms. • Sediment mixing and burrowing by animals strongly enhances P burial. • Supporting evidence for a decrease in oceanic P inventory in the early Paleozoic. Abstract A diagenetic model is used to simulate the diagenesis and burial of particulate organic carbon (Corg) and phosphorus (P) in marine sediments underlying anoxic versus oxic bottom waters. The latter are physically mixed by animals moving through the surface sediment (bioturbation) and ventilated by burrowing, tube-dwelling organisms (bioirrigation). The model is constrained using an empirical database including burial ratios of Corg with respect to organic P (Corg:Porg) and total reactive P (Corg:Preac), burial efficiencies of Corg and Porg, and inorganic carbon-to-phosphorus regeneration ratios. If Porg is preferentially mineralized relative to Corg during aerobic respiration, as many previous studies suggest, then the simulated Porg pool is found to be completely depleted. A modified model that incorporates the redox-dependent microbial synthesis of polyphosphates and Porg (termed the microbial P pump) allows preferential mineralization of the bulk Porg pool relative to Corg during both aerobic and anaerobic respiration and is consistent with the database. Results with this model show that P burial is strongly enhanced in sediments hosting fauna. Animals mix highly labile Porg away from the aerobic sediment layers where mineralization rates are highest, thereby mitigating diffusive PO43− fluxes to the bottom water. They also expand the redox niche where microbial P uptake occurs. The model was applied to a hypothetical shelf setting in the early Paleozoic; a time of the first radiation of benthic fauna. Results show that even shallow bioturbation at that time may have had a significant impact on P burial. Our model provides support for a recent study that proposed that faunal radiation in ocean sediments led to enhanced P burial and, possibly, a stabilization of atmospheric O2 levels. The results also help to explain Corg:Porg ratios in the geological record and the persistence of Porg in ancient marine sediments.

Description: Dissolved silicon isotope compositions have been analysed for the first time in pore waters (δ30SiPW) of three short sediment cores from the Peruvian margin upwelling region with distinctly different biogenic opal content in order to investigate silicon isotope fractionation behaviour during early diagenetic turnover of biogenic opal in marine sediments. The δ30SiPW varies between +1.1‰ and +1.9‰ with the highest values occurring in the uppermost part close to the sediment–water interface. These values are of the same order or higher than the δ30Si of the biogenic opal extracted from the same sediments (+0.3‰ to +1.2‰) and of the overlying bottom waters (+1.1‰ to +1.5‰). Together with dissolved silicic acid concentrations well below biogenic opal saturation, our collective observations are consistent with the formation of authigenic alumino-silicates from the dissolving biogenic opal. Using a numerical transport-reaction model we find that approximately 24% of the dissolving biogenic opal is re-precipitated in the sediments in the form of these authigenic phases at a relatively low precipitation rate of 56 μmol Si cm−2 yr−1. The fractionation factor between the precipitates and the pore waters is estimated at −2.0‰. Dissolved and solid cation concentrations further indicate that off Peru, where biogenic opal concentrations in the sediments are high, the availability of reactive terrigenous material is the limiting factor for the formation of authigenic alumino-silicate phases.

Description: During the last glaciation the Marmara Sea was isolated from the Mediterranean Sea because global sea level was below the depth of the Dardanelles sill. Prior to the postglacial reconnection to the Mediterranean Sea (~14.7 cal kyr BP), the surface waters of the Marmara Sea were brackish (Marmara Lake). Freshening of a previously saline Marmara Sea could have happened via spill-out of brackish to fresh water from the surface water of the Black Sea through the Bosphorus Strait. This hypothesis has not been tested against alternative possibilities (salt flushing by river run-off and precipitation). Here we use the dissolved Cl- and stable isotope composition (δO18 and δD) of Marmara Sea sediment pore water to estimate the salinity and stable isotope composition of Marmara Lake bottom water and to evaluate possible freshening scenarios. We use a transport model to simulate pore water Cl-, δO18 and δD in Marmara Sea sediments in the past 130 kyr, which includes the last interglacial (130-75 cal kyr BP), the last glacial (75-14 cal kyr BP) and the current postglacial period. Our results show that the bottom waters of the Marmara Lake were brackish (~4‰ salinity) and isotopically depleted (δO18~-10.2‰ and δD~-70‰, respectively) compared to modern seawater. Their salinity and stable isotope ratios show that they are a mixture of Mediterranean waters and Danube-like waters implying that the freshening took place via spill-out of freshwater through the Bosphorus. Our modelling approach indicates that the transit of fresh water from glacial Eurasia to the Mediterranean via the Marmara Sea started at least 50 cal kyr BP, was continuous throughout most of the last glaciation and persisted up to the post glacial reconnection to the Mediterranean through the Dardanelles sill (14.7 cal kyr BP). These results are consistent with previously published micropaleontological and geochemical investigations of sediment cores that indicate lacustrine conditions in the Marmara Sea from about 75 to 14.7 cal kyr BP.

Description: The accumulation of gas hydrates in marine sediments is essentially controlled by the accumulation of particulate organic carbon (POC) which is microbially converted into methane, the thickness of the gas hydrate stability zone (GHSZ) where methane can be trapped, the sedimentation rate (SR) that controls the time that POC and the generated methane stays within the GHSZ, and the delivery of methane from deep-seated sediments by ascending pore fluids and gas into the GHSZ. Recently, Wallmann et al. (2012) presented transfer functions to predict the gas hydrate inventory in diffusion-controlled geological systems based on SR, POC and GHSZ thickness for two different scenarios: normal and full compacting sediments. We apply these functions to global data sets of bathymetry, heat flow, seafloor temperature, POC input and SR, estimating a global mass of carbon stored in marine methane hydrates from 3 to 455 Gt of carbon (GtC) depending on the sedimentation and compaction conditions. The global sediment volume of the GHSZ in continental margins is estimated to be 60–67 × 1015 m3, with a total of 7 × 1015 m3 of pore volume (available for GH accumulation). However, seepage of methane-rich fluids is known to have a pronounced effect on gas hydrate accumulation. Therefore, we carried out a set of systematic model runs with the transport-reaction code in order to derive an extended transfer function explicitly considering upward fluid advection. Using averaged fluid velocities for active margins, which were derived from mass balance considerations, this extended transfer function predicts the enhanced gas hydrate accumulation along the continental margins worldwide. Different scenarios were investigated resulting in a global mass of sub-seafloor gas hydrates of ~ 550 GtC. Overall, our systematic approach allows to clearly and quantitatively distinguish between the effect of biogenic methane generation from POC and fluid advection on the accumulation of gas hydrate, and hence, provides a simple prognostic tool for the estimation of large-scale and global gas hydrate inventories in marine sediments.

Description: This study investigates the controls on organic carbon and molybdenum (Mo) accumulation in sediments deposited within the Western Interior Seaway across the Cenomanian–Turonian boundary interval (94.34–93.04 Ma) including Oceanic Anoxic Event 2 (OAE2). Carbon fluxes to the sediment–water interface (reflecting changes in primary productivity) and bottom-water oxygen concentrations (reflecting preservation effects) are reconstructed from field data and used to constrain a benthic model that simulates the geochemistry of unconsolidated sediments as they were deposited. The results show that increased availability of reactive iron prevents Mo sequestration as thiomolybdate (MoS42 −) during OAE2 (O2 ~ 105 μM) by (i) inhibiting sulfate reduction, and (ii) buffering any free sulfide that becomes available. In the post-OAE2 period (O2 ~ 50 μM), Mo accumulation is favored by a large reduction in iron flux. Importantly, this occurs in parallel with oxygenated bottom waters and high rates of aerobic carbon degradation in the surface sediments, implying that elevated Mo burial fluxes in ancient marine facies do not necessarily reflect euxinic or even anoxic conditions within the water column. Our findings suggest that both an increase in production and preservation lead to enrichment in organic carbon in the Western Interior Seaway. More generally, the results demonstrate that a careful consideration of the coupling between iron, carbon and oxygen cycles during the early stages of diagenesis is critical for interpreting geochemical proxies in modern and ancient settings.

Description: This study presents benthic data from 12 samplings from February to December 2010 in a 28 m deep channel in the southwest Baltic Sea. In winter, the distribution of solutes in the porewater was strongly modulated by bioirrigation which efficiently flushed the upper 10 cm of sediment, leading to concentrations which varied little from bottom water values. Solute pumping by bioirrigation fell sharply in the summer as the bottom waters became severely hypoxic (〈 2 μM O2). At this point the giant sulfide-oxidizing bacteria Beggiatoa was visible on surface sediments. Despite an increase in O2 following mixing of the water column in November, macrofauna remained absent until the end of the sampling. Contrary to expectations, metabolites such as dissolved inorganic carbon, ammonium and hydrogen sulfide did not accumulate in the upper 10 cm during the hypoxic period when bioirrigation was absent, but instead tended toward bottom water values. This was taken as evidence for episodic bubbling of methane gas out of the sediment acting as an abiogenic irrigation process. Porewater–seawater mixing by escaping bubbles provides a pathway for enhanced nutrient release to the bottom water and may exacerbate the feedback with hypoxia. Subsurface dissolved phosphate (TPO4) peaks in excess of 400 μM developed in autumn, resulting in a very large diffusive TPO4 flux to the water column of 0.7 ± 0.2 mmol m−2 d−1. The model was not able to simulate this TPO4 source as release of iron-bound P (Fe–P) or organic P. As an alternative hypothesis, the TPO4 peak was reproduced using new kinetic expressions that allow Beggiatoa to take up porewater TPO4 and accumulate an intracellular P pool during periods with oxic bottom waters. TPO4 is then released during hypoxia, as previous published results with sulfide-oxidizing bacteria indicate. The TPO4 added to the porewater over the year by organic P and Fe–P is recycled through Beggiatoa, meaning that no additional source of TPO4 is needed to explain the TPO4 peak. Further experimental studies are needed to strengthen this conclusion and rule out Fe–P and organic P as candidate sources of ephemeral TPO4 release. A measured C/P ratio of 〈 20 for the diffusive flux to the water column during hypoxia directly demonstrates preferential release of P relative to C under oxygen-deficient bottom waters. This coincides with a strong decrease in dissolved inorganic N/P ratios in the water column to ~ 1. Our results suggest that sulfide oxidizing bacteria could act as phosphorus capacitors in systems with oscillating redox conditions, releasing massive amounts of TPO4 in a short space of time and dramatically increasing the internal loading of TPO4 to the overlying water.