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: 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: 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: Geochemical data (CH4, SO42−, I−, Cl−, particulate organic carbon (POC), δ13C-CH4, and δ13C-CO2) are presented from the upper 30 m of marine sediment on a tectonic submarine accretionary wedge offshore southwest Taiwan. The sampling stations covered three ridges (Tai-Nan, Yung-An, and Good Weather), each characterized by bottom simulating reflectors, acoustic turbidity, and different types of faulting and anticlines. Sulfate and iodide concentrations varied little from seawater-like values in the upper 1–3 m of sediment at all stations; a feature that is consistent with irrigation of seawater by gas bubbles rising through the soft surface sediments. Below this depth, sulfate was rapidly consumed within 5–10 m by anaerobic oxidation of methane (AOM) at the sulfate-methane transition. Carbon isotopic data imply a mainly biogenic methane source. A numerical transport-reaction model was used to identify the supply pathways of methane and estimate depth-integrated turnover rates at the three ridges. Methane gas ascending from deep layers, facilitated by thrusts and faults, was by far the dominant term in the methane budget at all sites. Differences in the proximity of the sampling sites to the faults and anticlines mainly accounted for the variability in gas fluxes and depth-integrated AOM rates. By comparison, methane produced in situ by POC degradation within the modeled sediment column was unimportant. This study demonstrates that the geochemical trends in the continental margins offshore SW Taiwan are closely related to the different geological settings.

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: We present sedimentary geochemical data and in situ benthic flux measurements of dissolved inorganic nitrogen (DIN: NO3−, NO2−, NH4+) and oxygen (O2) from 7 sites with variable sand content along 18°N offshore Mauritania (NW Africa). Bottom water O2 concentrations at the shallowest station were hypoxic (42 μM) and increased to 125 μM at the deepest site (1113 m). Total oxygen uptake rates were highest on the shelf (−10.3 mmol O2 m−2 d−1) and decreased quasi-exponentially with water depth to −3.2 mmol O2 m−2 d−1. Average denitrification rates estimated from a flux balance decreased with water depth from 2.2 to 0.2 mmol N m−2 d−1. Overall, the sediments acted as net sink for DIN. Observed increases in δ15NNO3 and δ18ONO3 in the benthic chamber deployed on the shelf, characterized by muddy sand, were used to calculate apparent benthic nitrate fractionation factors of 8.0‰ (15εapp) and 14.1‰ (18εapp). Measurements of δ15NNO2 further demonstrated that the sediments acted as a source of 15N depleted NO2−. These observations were analyzed using an isotope box model that considered denitrification and nitrification of NH4+ and NO2−. The principal findings were that (i) net benthic 14N/15N fractionation (εDEN) was 12.9 ± 1.7‰, (ii) inverse fractionation during nitrite oxidation leads to an efflux of isotopically light NO2− (−22 ± 1.9‰), and (iii) direct coupling between nitrification and denitrification in the sediment is negligible. Previously reported εDEN for fine-grained sediments are much lower (4–8‰). We speculate that high benthic nitrate fractionation is driven by a combination of enhanced porewater–seawater exchange in permeable sediments and the hypoxic, high productivity environment. Although not without uncertainties, the results presented could have important implications for understanding the current state of the marine N cycle.

Description: A simple earth system model is developed to simulate global carbon and phosphorus cycling over the late Quaternary. It is focused on the geological cycling of C and P via continental weathering, volcanic and metamorphic degassing, hydrothermal processes and burial at the seabed. A simple ocean model is embedded in this geological model where the global ocean is represented by surface water, thermocline and deep water boxes. Concentrations of dissolved phosphorus, dissolved inorganic carbon, and total alkalinity are calculated for each box. The partial pressure of CO2 in the atmosphere (pCO(2A)) is determined by exchange processes with the surface ocean and the continents. It serves as key prognostic model variable and is assumed to govern surface temperatures and global sea-level. The model is formulated as autonomous system, in which the governing equations have no explicit time-dependence. For certain parameter values, the model does not converge towards a steady-state but develops stable self-sustained oscillations. These free oscillations feature pCO(2A) minima and maxima consistent with the ice-core record when vertical mixing in the ocean is allowed to vary in response to pCO(2A)-controlled temperature change. A stable 100-kyr cycle with a rapid transition from glacial to interglacial conditions is obtained when additional non-linear equations are applied to calculate deep ocean mixing, iron fertilization and the depth of organic matter degradation as function of pCO(2A)-controlled surface temperature. The delta C-13 value of carbon in the ocean/atmosphere system calculated in these model runs is consistent with the benthic delta C-13 record. However, the simulated C-13 depletion in the glacial ocean is not driven by the decline in terrestrial carbon stocks but by sea-level change controlling the rates of organic carbon burial and weathering at continental margins. The pCO(2A)-and delta C-13 oscillations develop without any form of external Milankovitch forcing. They are induced and maintained by sea-level change generating persistent imbalances in the marine carbon and phosphorus budgets. Stable oscillations are also obtained when sea-level change is allowed to lag temperature with a realistic time scale for ice-sheet adjustment

Description: An empirical function is derived for predicting the rate-depth profile of particulate organic carbon (POC) degradation in surface marine sediments including the bioturbated layer. The rate takes the form of a power law analogous to the Middelburg function. The functional parameters were optimized by simulating measured benthic O2 and NO3− fluxes at 185 stations worldwide using a diagenetic model. The novelty of this work rests with the finding that the vertically-resolved POC degradation rate in the bioturbated zone can be determined using a simple function where the POC rain rate is the governing variable. Although imperfect, the model is able to fit 71 % of paired O2 and NO3− fluxes to within 50% of measured values. It further provides realistic geochemical concentration-depth profiles, NO3− penetration depths and apparent first-order POC mineralization rate constants. The model performs less well on the continental shelf due to the high heterogeneity there. When applied to globally resolved maps of rain rate, the model predicts a global denitrification rate of 182 ± 88 Tg yr−1 of N and a POC burial rate of 107 ± 52 Tg yr−1 of C with a mean carbon burial efficiency of 6.1%. These results are in very good agreement with published values. Our proposed function is conceptually simple, requires less parameterization than multi-G type models and is suitable for non-steady state applications. It provides a basis for more accurately simulating benthic nutrient fluxes and carbonate dissolution rates in Earth system models.