Search

Carbon dioxide (CO

Highlights • We present a 5 myr record of biogeochemical cycling in a Cretaceous upwelling area. • A novel quantitative approach for the evaluation of Fe speciation proxies was applied. • Ferruginous proxy signature reflects intense chemical weathering rather than anoxia. • Water column redox conditions evolved from oxic to nitrogenous to euxinic before OAE2. • Smaller seawater nitrate inventory facilitated sedimentary H2S release and euxinia. Abstract Oceanic Anoxic Events (OAEs) in Earth's history are regarded as analogues for current and future ocean deoxygenation, potentially providing information on its pacing and internal dynamics. In order to predict the Earth system's response to changes in greenhouse gas concentrations and radiative forcing, a sound understanding of how biogeochemical cycling differs in modern and ancient marine environments is required. Here, we report proxy records for iron (Fe), sulfur and nitrogen cycling in the Tarfaya upwelling system in the Cretaceous Proto-North Atlantic before, during and after OAE2 (∼93 Ma). We apply a novel quantitative approach to sedimentary Fe speciation, which takes into account the influence of terrigenous weathering and sedimentation as well as authigenic Fe (non-terrigenous, precipitated onsite) rain rates on Fe-based paleo-redox proxies. Generally elevated ratios of reactive Fe (i.e., bound to oxide, carbonate and sulfide minerals) to total Fe (FeHR/FeT) throughout the 5 million year record are attributed to transport-limited chemical weathering under greenhouse climate conditions. Trace metal and nitrogen isotope systematics indicate a step-wise transition from oxic to nitrogenous to euxinic conditions over several million years prior to OAE2. Taking into consideration the low terrigenous sedimentation rates in the Tarfaya Basin, we demonstrate that highly elevated FeHR/FeT from the mid-Cenomanian through OAE2 were generated with a relatively small flux of additional authigenic Fe. Evaluation of mass accumulation rates of reactive Fe in conjunction with the extent of pyritization of reactive Fe reveals that authigenic Fe and sulfide precipitation rates in the Tarfaya Basin were similar to those in modern upwelling systems. Because of a smaller seawater nitrate inventory, however, chemolithoautotrophic sulfide oxidation with nitrate was less efficient in preventing hydrogen sulfide release into the water column. As terrigenous weathering and sediment flux determine how much authigenic Fe is required to generate an anoxic euxinic or ferruginous proxy signature, we emphasize that both have to be taken into account when interpreting Fe-based paleo-redox proxies.

Periodic changes in sediment composition are usually ascribed to insolation forcing controlled by Earth’s orbital parameters. During the Cretaceous Thermal Maximum at 97–91 Myr ago (Ma), a 37–50-kyr-long cycle that is generally believed to reflect obliquity forcing dominates the sediment record. Here, we use a numerical ocean model to show that a cycle of this length can be generated by marine biogeochemical processes without applying orbital forcing. According to our model, the restricted proto-North Atlantic and Tethys basins were poorly ventilated and oscillated between iron-rich and sulfidic (euxinic) states. The Panthalassa Basin was fertilized by dissolved iron originating from the proto-North Atlantic. Hence, it was less oxygenated while the proto-North Atlantic was in an iron-rich state and better oxygenated during euxinic periods in the proto-North Atlantic. This redox see-saw was strong enough to create significant changes in atmospheric $$p_{\mathrm {CO}_2}$$ . We conclude that most of the variability in the mid-Cretaceous ocean–atmosphere system can be ascribed to the internal redox see-saw and its response to external orbital forcing. An internal redox see-saw between the Panthalassa Basin and the proto-North Atlantic can explain cyclic changes in the sediment record throughout the Cretaceous Thermal Maximum 97 to 91 million years ago, according to simulations with a numerical ocean model.

Highlights • Amazon shelf sediments take up seawater potassium (K) due to reverse weathering. • Amazon shelf sediments release terrigenous phosphorus (P) during resuspension. • Updated estimates of sedimentary K uptake and P release on Amazon shelf are presented. • Sedimentary K uptake on Amazon shelf corresponds to 13% of global riverine K input. • Sedimentary P release is ~5 times higher than dissolved P discharge of Amazon River. Abstract In this study, we identify and quantify processes that lead to sedimentary potassium (K) sequestration and phosphorus (P) release on the Amazon shelf. To this end, seven short sediment cores were recovered from the Amazon shelf during R/V Meteor cruise M147. All of the sediment cores investigated in this study are characterized by elevated K to aluminum (Al) ratios compared to Amazon riverine suspended matter, which indicates that seawater K+ is incorporated into the solid phase on the entire Amazon shelf. Pore water silica (Si) profiles are characterized by irregularly increasing concentrations and plateaus, thus, deviating from the asymptotic shape that is typically found in continental margin sediments. At one site, a dissolved Si plateau coincides with a K+ minimum suggesting that these solutes are incorporated into authigenic minerals, a process referred to as reverse weathering. Previous flux estimates for elements that participate in reverse weathering on the Amazon shelf were derived from pore water diffusive fluxes, reaction rates estimated from sediment incubations and solid phase extractions. In this study, we took an alternative approach, which is based on the concentration difference between shelf sediments and river suspended particles. The resulting K flux due to reverse weathering of 1.7 ∙ 1011 mol yr−1 is in agreement with previous estimates and corresponds to 13% of the global riverine dissolved K+ input. Previous studies demonstrated that Amazon riverine particulate P is partly solubilized on the Amazon shelf. However, these results are exclusively based on sediment data close to the river mouth and no distinction between terrestrial and marine sediment components was made. Here, we quantify P release from Amazon shelf sediments by comparing terrestrial P concentrations in shelf sediments with P concentrations in river suspended particles. The resulting solubilized P flux of 2.2 ∙ 1010 mol yr−1 is about five to six times higher than previous estimates and about seven times the Amazon riverine dissolved P discharge. The magnitudes of the presented fluxes imply that the alteration of riverine shelf sediments significantly affects the mean concentrations of dissolved K+ and P in the global ocean.

Highlights • MeBo drilling in Danube fan down to 147 m recovered limnic to marine deposits. • Molecular and stable isotope characterization of light hydrocarbons, CO2, and H2O. • H and O isotopic compositions of pore water reflect paleoclimate variations. • Isotope relations prove microbial carbonate reduction as major methanogenic pathway. • Control of δ2H–CH4 by δ2H–H2O may lead to misinterpretation of methanogenic paths. Abstract We report on the geochemistry of light hydrocarbons and pore water in sediments down to 147 m below seafloor (mbsf), at two sites within the gas hydrate stability field of the Danube deep-sea fan, Black Sea. Sediments were drilled with MARUM-MeBo200 and comprise the transition from limnic to the recent marine stage. Stable C/N ratios (mean 5.1 and 5.6) and δ13C-Corg values (mean −25.8‰ V-PDB) suggest relatively uniform bulk organic matter compositions. In contrast, pore water δ2H and δ18O values varied considerably from approx. −120‰ to −30‰ V-SMOW and from −15‰ to −3‰ V-SMOW, respectively. These data pairs plot close to the ‘Global Meteoric Water Line’ and indicate paleo temperature variations. Depletions of pore water in 2H and 18O below 40 mbsf indicate low temperatures and likely reflect conditions during (the) last glacial period(s). Methane was much more abundant than the only other hydrocarbons found in notable concentrations, ethane and propane ((C1/(C2+C3) ≥20,000). Relatively constant δ13C–CH4 (~−70‰ V-PDB) and δ13C–C2H6 (~−52‰ V-PDB) values with depth indicate that methane and ethane are predominantly of microbial origin and that their formation was not limited by carbon availability. In contrast, δ2H–CH4 values varied in a large range (approx. −310 to −240‰ V-SMOW) with depth and positively correlated with trends observed for δ2H–H2O. Isotope separations (Δδ13C(CH4–CO2), Δδ2H(CH4–H2O)) substantiate that microbial carbonate reduction (CR) is the prevalent methanogenic pathway throughout the sediments irrespective of their geochemical history. Remarkably, in δ13C–CH4 – δ2H–CH4 diagrams widely used, samples characterized by δ2H–CH4 values more negative than approx. −250‰ plot out of the field assigned for pure CR. We conclude that assignments of microbial methanogenic pathways based on classical interpretations of δ13C–CH4 – δ2H–CH4 pairs can lead to misinterpretations, as severe 2H-depletions of methane formed through microbial CR can result from 2H-depletions of the pore water generated during low-temperature climatic periods.

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.

A diagenetic model is used to simulate the diagenesis and burial of particulate organic carbon (C org ) 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 C org with respect to organic P (C org :P org ) and total reactive P (C org :P reac ), burial efficiencies of C org and P org , and inorganic carbon-to-phosphorus regeneration ratios. If P org is preferentially mineralized relative to C org during aerobic respiration, as many previous studies suggest, then the simulated P org pool is found to be completely depleted. A modified model that incorporates the redox-dependent microbial synthesis of polyphosphates and P org (termed the microbial P pump) allows preferential mineralization of the bulk P org pool relative to C org 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 P org away from the aerobic sediment layers where mineralization rates are highest, thereby mitigating diffusive PO 4 3− 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 O 2 levels. The results also help to explain C org :P org ratios in the geological record and the persistence of P org in ancient marine sediments.

We have developed and employed an Earth system model to explore the forcings of atmospheric pCO2 change and the chemical and isotopic evolution of seawater over the last glacial cycle. Concentrations of dissolved phosphorus (DP), reactive nitrogen, molecular oxygen, dissolved inorganic carbon (DIC), total alkalinity (TA), 13C-DIC, and 14C-DIC were calculated for 24 ocean boxes. The bi-directional water fluxes between these model boxes were derived from a 3-D circulation field of the modern ocean (Opa 8.2, NEMO) and tuned such that tracer distributions calculated by the box model were consistent with observational data from the modern ocean. To model the last 130 kyr, we employed records of past changes in sea-level, ocean circulation, and dust deposition. According to the model, about half of the glacial pCO2 drawdown may be attributed to marine regressions. The glacial sea-level low-stands implied steepened ocean margins, a reduced burial of particulate organic carbon, phosphorus, and neritic carbonate at the margin seafloor, a decline in benthic denitrification, and enhanced weathering of emerged shelf sediments. In turn, low-stands led to a distinct rise in the standing stocks of DIC, TA, and nutrients in the global ocean, promoted the glacial sequestration of atmospheric CO2 in the ocean, and added 13C- and 14C-depleted DIC to the ocean as recorded in benthic foraminifera signals. The other half of the glacial drop in pCO2 was linked to inferred shoaling of Atlantic meridional overturning circulation and more efficient utilization of nutrients in the Southern Ocean. The diminished ventilation of deep water in the glacial Atlantic and Southern Ocean led to significant 14C depletions with respect to the atmosphere. According to our model, the deglacial rapid and stepwise rise in atmospheric pCO2 was induced by upwelling both in the Southern Ocean and subarctic North Pacific and promoted by a drop in nutrient utilization in the Southern Ocean. The deglacial sea-level rise led to a gradual decline in nutrient, DIC, and TA stocks, a slow change due to the large size and extended residence times of dissolved chemical species in the ocean. Thus, the rapid deglacial rise in pCO2 can be explained by fast changes in ocean dynamics and nutrient utilization whereas the gradual pCO2 rise over the Holocene may be linked to the slow drop in nutrient and TA stocks that continued to promote an ongoing CO2 transfer from the ocean into the atmosphere.

Previous studies have suggested that weathering and benthic phosphorus (P) fluxes, triggered by climate warming, can increase the oceanic P inventory on millennial time scales, promoting ocean productivity and deoxygenation. In this study, we assessed the major uncertainties in projected P inventories and their imprint on ocean deoxygenation using an Earth system model of intermediate complexity for a business-as-usual carbon dioxide (CO2) emission scenario until year 2300 and subsequent linear decline to zero emissions until year 3000. Model results suggest a large spread in the simulated oceanic P inventory due to uncertainties in (1) assumptions for weathering parameters, (2) the representation of bathymetry on slopes and shelves in the model bathymetry, (3) the parametrization of benthic P fluxes and (4) the representation of sediment P inventories. Our best estimate for changes in the global ocean P inventory by the year 5000 caused by global warming amounts to +30 % compared to pre-industrial levels. Weathering, benthic and anthropogenic fluxes of P contributed +25 %, +3 % and +2 % respectively. The total range of oceanic P inventory changes across all model simulations varied between +2 % and +60 %. Suboxic volumes were up to 5 times larger than in a model simulation with a constant oceanic P inventory. Considerably large amounts of the additional P left the ocean surface unused by phytoplankton via physical transport processes as preformed P. Nitrogen fixation was not able to adjust the oceanic nitrogen inventory to the increasing P levels or to compensate for the nitrogen loss due to increased denitrification. This is in contrast to palaeo reconstructions of large-scale deoxygenation events. We suggest that uncertainties in P weathering, nitrogen fixation and benthic P feedbacks need to be reduced to achieve more reliable projections of oceanic deoxygenation on millennial timescales.

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.

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.

An edited version of this paper was published by AGU. Copyright 2018 American Geophysical Union. Source at https://doi.org/10.1002/2017GC007288. During expedition MARIA S. MERIAN MSM57/2 to the Svalbard margin offshore Prins Karls Forland, the seafloor drill rig MARUM‐MeBo70 was used to assess the landward termination of the gas hydrate system in water depths between 340 and 446 m. The study region shows abundant seafloor gas vents, clustered at a water depth of ∼400 m. The sedimentary environment within the upper 100 m below seafloor (mbsf) is dominated by ice‐berg scours and glacial unconformities. Sediments cored included glacial diamictons and sheet‐sands interbedded with mud. Seismic data show a bottom simulating reflector terminating ∼30 km seaward in ∼760 m water depth before it reaches the theoretical limit of the gas hydrate stability zone (GHSZ) at the drilling transect. We present results of the first in situ temperature measurements conducted with MeBo70 down to 28 mbsf. The data yield temperature gradients between ∼38°C km−1 at the deepest site (446 m) and ∼41°C km−1 at a shallower drill site (390 m). These data constrain combined with in situ pore‐fluid data, sediment porosities, and thermal conductivities the dynamic evolution of the GHSZ during the past 70 years for which bottom water temperature records exist. Gas hydrate is not stable in the sediments at sites shallower than 390 m water depth at the time of acquisition (August 2016). Only at the drill site in 446 m water depth, favorable gas hydrate stability conditions are met (maximum vertical extent of ∼60 mbsf); however, coring did not encounter any gas hydrates.

Carbon cycling and fluid seepage in marine sediments over the late Quaternary were investigated at a now-extinct pockmark located in a mega-pockmark field in the SW Xisha Uplift (NW South China Sea). Measured particulate organic carbon (POC) content, and porewater sulfate (SO4 2-), dissolved inorganic carbon (DIC) concentrations and δ34S-SO4 2- distributions were used to constrain a non-steady-state reaction-transport model and quantify POC mineralization rates as well as estimate the time when fluid flow ceased at the investigated pockmark. An increase in POC content and δ34S-SO4 2- and a decrease in sulfate concentrations in the upper ca. 2m at the pockmark and a reference core implied an increase in the flux and reactivity of organic matter during the early Holocene around 10kyrB. P. caused by enhanced primary productivity during the strengthened southwestern summer monsoon. These features were simulated with the model assuming a Holocene increase in POC flux and reactivity. Subsequently, starting from an initial condition reminiscent of a modern active cold seep (Hydrate Ridge), hindcast simulations showed that fluid seepage at the pockmark ceased ca. 39kyr ago. This corresponds to a relative sea level high-stand, which is believed to be associated with gas hydrate stabilization. The non-steady-state model presented in this contribution can also be used to constrain the time when fluid seepage ceased at other presently extinct cold seeps when suitable sediment and porewater data are available.

Literature data on benthic dissolved iron (DFe) fluxes (µmol m−2 d−1), bottom water oxygen concentrations (O2BW, μM), and sedimentary carbon oxidation rates (COX, mmol m−2 d−1) from water depths ranging from 80 to 3700 m were assembled. The data were analyzed with a diagenetic iron model to derive an empirical function for predicting benthic DFe fluxes: inline image where γ (= 170 µmol m−2 d−1) is the maximum flux for sediments at steady state located away from river mouths. This simple function unifies previous observations that COX and O2BW are important controls on DFe fluxes. Upscaling predicts a global DFe flux from continental margin sediments of 109 ± 55 Gmol yr−1, of which 72 Gmol yr−1 is contributed by the shelf ( 2000 m) of 41 ± 21 Gmol yr−1 is unsupported by empirical data. Previous estimates of benthic DFe fluxes derived using global iron models are far lower (approximately 10–30 Gmol yr−1). This can be attributed to (i) inadequate treatment of the role of oxygen on benthic DFe fluxes and (ii) improper consideration of continental shelf processes due to coarse spatial resolution. Globally averaged DFe concentrations in surface waters simulated with the intermediate-complexity University of Victoria Earth System Climate Model were a factor of 2 higher with the new function. We conclude that (i) the DFe flux from marginal sediments has been underestimated in the marine iron cycle and (ii) iron scavenging in the water column is more intense than currently presumed.

Carbon cycling in Peruvian margin sediments (11° S and 12° S) was examined at 16 stations from 74 m on the inner shelf down to 1024 m water depth by means of in situ flux measurements, sedimentary geochemistry and modeling. Bottom water oxygen was below detection limit down to ca. 400 m and increased to 53 μM at the deepest station. Sediment accumulation rates and benthic dissolved inorganic carbon fluxes decreased rapidly with water depth. Particulate organic carbon (POC) content was lowest on the inner shelf and at the deep oxygenated stations (< 5%) and highest between 200 and 400 m in the oxygen minimum zone (OMZ, 15–20%). The organic carbon burial efficiency (CBE) was unexpectedly low on the inner shelf (< 20%) when compared to a global database, for reasons which may be linked to the frequent ventilation of the shelf by oceanographic anomalies. CBE at the deeper oxygenated sites was much higher than expected (max. 81%). Elsewhere, CBEs were mostly above the range expected for sediments underlying normal oxic bottom waters, with an average of 51 and 58% for the 11° S and 12° S transects, respectively. Organic carbon rain rates calculated from the benthic fluxes alluded to a very efficient mineralization of organic matter in the water column, with a Martin curve exponent typical of normal oxic waters (0.88 ± 0.09). Yet, mean POC burial rates were 2–5 times higher than the global average for continental margins. The observations at the Peruvian margin suggest that a lack of oxygen does not affect the degradation of organic matter in the water column but promotes the preservation of organic matter in marine sediments.

Shallow gas migration along hydrocarbon wells constitutes a potential methane emission pathway that currently is not recognized in any regulatory framework or greenhouse gas inventory. Recently, the first methane emission measurements at three abandoned offshore wells in the Central North Sea (CNS) were conducted showing that considerable amounts of biogenic methane originating from shallow gas accumulations in the overburden of deep reservoirs were released by the boreholes. Here, we identify numerous wells poking through shallow gas pockets in 3-D seismic data of the CNS indicating that about one-third of the wells may leak, potentially releasing a total of 3–17 kt of methane per year into the North Sea. This poses a significant contribution to the North Sea methane budget. A large fraction of this gas (∼42%) may reach the atmosphere via direct bubble transport (0–2 kt yr–1) and via diffusive exchange of methane dissolving in the surface mixed layer (1–5 kt yr–1), as indicated by numerical modeling. In the North Sea and in other hydrocarbon-prolific provinces of the world shallow gas pockets are frequently observed in the sedimentary overburden and aggregate leakages along the numerous wells drilled in those areas may be significant.

Methane seepage from the upper continental slopes of Western Svalbard has previously been attributed to gas hydrate dissociation induced by anthropogenic warming of ambient bottom waters. Here we show that sediment cores drilled off Prins Karls Foreland contain freshwater from dissociating hydrates. However, our modeling indicates that the observed pore water freshening began around 8 ka BP when the rate of isostatic uplift outpaced eustatic sea-level rise. The resultant local shallowing and lowering of hydrostatic pressure forced gas hydrate dissociation and dissolved chloride depletions consistent with our geochemical analysis. Hence, we propose that hydrate dissociation was triggered by postglacial isostatic rebound rather than anthropogenic warming. Furthermore, we show that methane fluxes from dissociating hydrates were considerably smaller than present methane seepage rates implying that gas hydrates were not a major source of methane to the oceans, but rather acted as a dynamic seal, regulating methane release from deep geological reservoirs., Methane seepage from continental slopes has been attributed to gas hydrate dissociation induced by anthropogenic bottom water warming. Here, the authors show that hydrates dissociated before the Anthropocene when the isostatic rebound induced by deglaciation of the Arctic ice sheet outpaced eustatic sea-level rise.

Monthly time-series data (1998–2009) of bottom water oxygen, nitrate and nitrite concentrations from the outer shelf (150 m water depth) in the oxygen minimum zone offshore Peru were coupled to a layered biogeochemical sediment model to investigate benthic-pelagic coupling over multi-annual time scales. The model includes the mineralization of four reactive pools of particulate organic carbon (POC) with lifetimes of 0.13, 1.3, 20, and 1700 year that were constrained using empirical data. Total POC rain rates to the seafloor were derived from satellite based estimates of primary production. Solute fluxes and concentrations in sediment porewater showed highly dynamic behavior over the course of a typical year. Conversion of fixed N to N2 by denitrification varied from 1.1 mmol m−2 d−1 of N in winter to 1.8 mmol m−2 d−1 of N in summer with a long term mean N loss for the shelf of 1.5 mmol m−2 d−1 of N. Fixed N loss across the whole time-series agreed very well with a previously-published vertically-integrated sediment model for coupling the benthic and pelagic N cycle in regional and global models. Dissimilatory nitrate reduction to ammonium (DNRA) emerges as a major process in the benthic N cycle, producing on average 1.9 mmol m−2 d−1 of ammonium: more than twice the rate of ammonification of organic nitrogen. The model predicts sulfide emissions from the sediment of up to 1 mmol m−2 d−1 when POC rain rate exceeds 20 mmol m−2 d−1, in agreement with past observations of benthic sulfide fluxes and sulfide plume distributions in the water column. This study demonstrates that sediments on the Peruvian shelf are not static repositories that are independent of changes taking place in the water column. Our results strongly suggest the shelf sediments must exert an important feedback on biogeochemical processes in the overlying waters, and should be considered in regional model studies.

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

The isotopic composition of Phanerozoic marine sediments provides important information about changes in seawater chemistry. In particular, the radiogenic strontium isotope (87Sr/86Sr) system is a powerful tool for constraining plate tectonic processes and their influence on atmospheric CO2 concentrations. However, the 87Sr/86Sr isotope ratio of seawater is not sensitive to temporal changes in the marine strontium (Sr) output flux, which is primarily controlled by the burial of calcium carbonate (CaCO3) at the ocean floor. The Sr budget of the Phanerozoic ocean, including the associated changes in the amount of CaCO3 burial, is therefore only poorly constrained. Here, we present the first stable isotope record of Sr for Phanerozoic skeletal carbonates, and by inference for Phanerozoic seawater (δ88/86Srsw), which we find to be sensitive to imbalances in the Sr input and output fluxes. This δ88/86Srsw record varies from ∼0.25‰ to ∼0.60‰ (vs. SRM987) with a mean of ∼0.37‰. The fractionation factor between modern seawater and skeletal calcite Δ88/86Srcc-sw, based on the analysis of 13 modern brachiopods (mean δ88/86Sr of 0.176±0.016‰, 2 standard deviations (s.d.)), is -0.21‰ and was found to be independent of species, water temperature, and habitat location. Overall, the Phanerozoic δ88/86Srsw record is positively correlated with the Ca isotope record (δ44/40Casw), but not with the radiogenic Sr isotope record ((87Sr/86Sr)sw). A new numerical modeling approach, which considers both δ88/86Srsw and (87Sr/86Sr)sw, yields improved estimates for Phanerozoic fluxes and concentrations for seawater Sr. The oceanic net carbonate flux of Sr (F(Sr)carb) varied between an output of -4.7x1010mol/Myr and an input of +2.3x1010mol/Myr with a mean of -1.6x1010mol/Myr. On time scales in excess of 100Myrs the F(Sr)carb is proposed to have been controlled by the relative importance of calcium carbonate precipitates during the “aragonite” and “calcite” sea episodes. On time scales less than 20Myrs the F(Sr)carb seems to be controlled by variable combinations of carbonate burial rate, shelf carbonate weathering and recrystallization, ocean acidification, and ocean anoxia. In particular, the Permian/Triassic transition is marked by a prominent positive δ88/86Srsw-peak that reflects a significantly enhanced burial flux of Sr and carbonate, likely driven by bacterial sulfate reduction (BSR) and the related alkalinity production in deeper anoxic waters. We also argue that the residence time of Sr in the Phanerozoic ocean ranged from ∼1Myrs to ∼20Myrs.

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.

Most of the marine biotic crises that occurred during the hot Mesozoic era have been linked to episodes of extreme warmth(1,2). Others, however, may have occurred during cooler intervals that interrupted Cretaceous greenhouse warmth(3-5). There are some indications of cooling in the late Aptian(6-8) (116-114 Myr ago), but it has not been definitively linked to biotic crisis. Here we assess the timing and magnitude of late Aptian cooling and its association with biotic crises using a suite of geochemical and micropalaeontological assessments from a marine sediment core from the North Atlantic Ocean as well as global biogeochemical modelling. Sea surface temperatures derived from the TEX86 proxy suggest that surface waters cooled by about 5 degrees C during the two million years, coincident with a positive delta C-13 excursion of approximately 2 parts per thousand in carbonates and organic carbon. Surface productivity was enhanced during this period, but the abundance of planktonic foraminifera and nannoconid phytoplankton declined. Our simulations with a biogeochemical model indicate that the delta C-13 excursion associated with the cooling could be explained by the burial of about 812,000 gigatons of carbon over 2.5 million years. About 50% of the this carbon burial occurred in the Atlantic, Southern and Tethys ocean basins. We conclude that global cooling during greenhouse conditions can cause perturbations to marine ecosystems and biogeochemical cycles at scales comparable to those associated with global warming

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.

Seafloormethane emission from the Quepos Slide on the submarine segment of the Costa Rica fore-arc margin was estimated by extrapolating flux measurements from individual seeps to the total area covered by bacterial mats. This approach is based on the combination of detailed mapping to determine the abundance of seeps and the application of a numerical model to estimate the amount of benthic methane fluxes. Model results suggest that the majority of the studied seeps transport rather limited amount of methane (on average: *177 lmol cm-2 a-1) into the water column due to moderate upward advection, allowing for intense anaerobic oxidation of methane (AOM; on average: 53 % of the methane flux is consumed). Depth-integrated AOM rates (56–1,538 lmol CH4 cm-2 a-1) are comparable with values reported from other active seep sites. The overall amount of dissolved methane released into the water column from the entire area covered by bacterial mats on the Quepos Slide is estimated to be about 0.28 9 106 mol a-1. This conservative estimate which relies on rather accurate determinations of seafloor methane fluxes emphasizes the potential importance of submarine slides as sites of natural methane seepage; however, at present the global extent of methane seepage from submarine slides is largely unknown.

Highlights • Very high rates of dissimilatory nitrate reduction to ammonium by Thioploca. • Non-steady state model predicts Thioploca survival on intracellular nitrate reservoir. • Ammonium release by Thioploca may be coupled to pelagic N loss by anammox. • Thioploca may contribute to anammox long after bottom water nitrate disappearance. • Model indicates that benthic foraminifera account for 90% of benthic N2 production. Abstract Benthic N cycling in the Peruvian oxygen minimum zone (OMZ) was investigated at ten stations along 12oS from the middle shelf (74 m) to the upper slope (1024 m) using in situ flux measurements, sediment biogeochemistry and modelling. Middle shelf sediments were covered by mats of the filamentous bacteria Thioploca spp. and contained a large ‘hidden’ pool of nitrate that was not detectable in the porewater. This was attributed to a biological nitrate reservoir stored within the bacteria to oxidize sulfide to sulfate during ‘dissimilatory nitrate reduction to ammonium’ (DNRA). The extremely high rates of DNRA on the shelf (15.6 mmol m-2 d-1 of N), determined using an empirical steady-state model, could easily supply all the ammonium requirements for anammox in the water column. The model further showed that denitrification by foraminifera may account for 90% of N2 production at the lower edge of the OMZ. At the time of sampling, dissolved oxygen was below detection limit down to 400 m and the water body overlying the shelf had stagnated, resulting in complete depletion of nitrate and nitrite. A decrease in the biological nitrate pool was observed on the shelf during fieldwork concomitant with a rise in porewater sulfide levels in surface sediments to 2 mM. Using a non-steady state model to simulate this natural anoxia experiment, these observations were shown to be consistent with Thioploca surviving on a dwindling intracellular nitrate reservoir to survive the stagnation period. The model shows that sediments hosting Thioploca are able to maintain high ammonium fluxes for many weeks following stagnation, potentially sustaining pelagic N loss by anammox. In contrast, sulfide emissions remain low, reducing the economic risk to the Peruvian fishery by toxic sulfide plume development.

Current de-oxygenation of the oceans is associated with severe habitat loss and distinct changes in the species composition of bentho-pelagic communities. We investigated the distributions of epibenthic megafauna across the Peruvian OMZ (11°S) at water depths ranging from ∼80 to 1000 m water depth using sea floor images. Likely controls of distributions were adressed by combining the abundances of major groups with geochemical parameters and sea-floor topography. In addition to bottom-water oxygen levels and organic-carbon availability, particular emphasis is laid on the effects of local hydrodynamics. Beside the occurrence of microbial mats at the shelf and upper slope, distinct zones of highly abundant megafauna, dominated by gastropods (900 ind. m−2), ophiuroids (140 ind. m−2), and pennatulaceans (20 ind. m−2), were observed at the lower boundary of the OMZ. Their distribution extended from 460 m water depth (O2 levels < 2 μM), where gastropods were abundant, to 680 m (O2 ∼6 μM) where epifaunal abundances declined sharply. Bottom water O2 represents a major factor that limits the ability of metazoans to invade deeply into the OMZ where they could have access to labile organic carbon. However, depending on their feeding mode, the distribution of organisms appeared to be related to local hydrodynamics caused by the energy dissipation of incipient internal M2 tides affecting the suspension, transport and deposition of food particles. This was particularly evident in certain sections of the investigated transect. At these potentially critical sites, energy dissipation of internal tides is associated with high bottom shear stress and high turbulences and coincides with elevated turbidity levels in the benthic boundary layer, increased Zr/Al-ratios, low sedimentation rates as well as a shift in the grain size towards coarser particles. In or near such areas, abundant suspension-feeding organisms, such as ophiuroids, pennatulaceans, and tunicates were present, whereas deposit-feeding gastropods were absent. The influence of local hydrodynamic conditions on the distribution of epibenthic organisms has been neglected in OMZ studies, although it has been considered in other settings.

The accumulation of methane hydrate in marine sediments is controlled by a number of physical and biogeochemical parameters including the thickness of the gas hydrate stability zone (GHSZ), the solubility of methane in pore fluids, the accumulation of particulate organic carbon at the seafloor, the kinetics of microbial organic matter degradation and methane generation in marine sediments, sediment compaction and the ascent of deep-seated pore fluids and methane gas into the GHSZ. Our present knowledge on these controlling factors is discussed and new estimates of global sediment and methane fluxes are provided applying a transport-reaction model at global scale. The modeling and the data evaluation yield improved and better constrained estimates of the global pore volume within the modern GHSZ ( ≥ 44 × 1015 m3), the Holocene POC accumulation rate at the seabed (~1.4 × 1014 g yr−1), the global rate of microbial methane production in the deep biosphere (4−25 × 1012 g C yr−1) and the inventory of methane hydrates in marine sediments ( ≥ 455 Gt of methane-bound carbon).

Pore water and solid phase data for redox-sensitive metals (Mn, Fe, V, Mo and U) were collected on a transect across the Peru upwelling area (11°S) at water depths between 78 and 2025 m and bottom water oxygen concentrations ranging from ~0 to 93 µM. By comparing authigenic mass accumulation rates and diffusive benthic fluxes, we evaluate the respective mechanisms of trace metal accumulation, retention and remobilization across the oxygen minimum zone (OMZ) and with respect to oxygen fluctuations in the water column related to the El Nino Southern Oscillation (ENSO). Sediments within the permanent OMZ are characterized by diffusive uptake and authigenic fixation of U, V and Mo as well as diffusive loss of Mn and Fe across the benthic boundary. Some of the dissolved Mn and Fe in the water column re-precipitate at the oxycline and shuttle particle-reactive trace metals to the sediment surface at the lower and upper boundary of the OMZ. At the lower boundary, pore waters are not sufficiently sulfidic as to enable an efficient authigenic V and Mo fixation. As a consequence, sediments below the OMZ are preferentially enriched in U which is delivered via both in situ pre-cipitation and lateral supply of U-rich phosphorites from further upslope. Trace metal cycling on the Peruvian shelf is strongly affected by ENSO-related oxygen fluctuations in bottom water. During periods of shelf oxygenation, surface sediments receive particulate V and Mo with metal (oxyhydr)oxides that derive from both terrigenous sources and precipitation at the retreating oxycline. After the recurrence of anoxic conditions, metal (oxyhydr)oxides are reductively dissolved and the hereby liberated V and Mo are authigenically removed. This alternation between supply of particle-reactive trace metals during oxic periods and fixation during anoxic periods leads to a preferential accumulation of V and Mo compared to U on the Peruvian shelf. The decoupling of V, Mo and U accumulation is further accentuated by the varying susceptibility to re-oxidation of the different authigenic metal phases. While authigenic U and V are readily re-oxidized and recycled during periods of shelf oxygenation, the sequestration of Mo by authigenic pyrite is favored by the transient occurrence of oxidizing conditions.Our findings reveal that redox-sensitive trace metals respond in specific manner to short-term oxygen fluctuations in the water column. The relative enrichment patterns identified might be useful for the reconstruction of past OMZ extension and large-scale redox oscillations in the geological record.

This study investigates the biogeochemical processes that control the benthic fluxes of dissolved nitrogen (N) species in Boknis Eck - a 28 m deep site in the Eckernförde Bay (southwestern Baltic Sea). Bottom water oxygen concentrations (O2-BW) fluctuate greatly over the year at Boknis Eck, being well-oxygenated in winter and experiencing severe bottom water hypoxia and even anoxia in late summer. The present communication addresses the winter situation (February 2010). Fluxes of ammonium (NH4+), nitrate (NO3-) and nitrite (NO2-) were simulated using a benthic model that accounted for transport andbiogeochemical reactions and constrained with ex situ flux measurements and sediment geochemical analysis. The sediments were a net sink for NO3- (-0.35 mmol m-2 d-1 of NO3-), of which 75% was ascribed to dissimilatory reduction of nitrate to ammonium (DNRA) by sulfide oxidizing bacteria, and 25% to NO3- reduction to NO2- by denitrifying microorganisms. NH4+ fluxes were high (1.74 mmol m-2d-1 of NH4+), mainly due to the degradation of organic nitrogen, and directed out of the sediment. NO2-fluxes were negligible. The sediments in Boknis Eck are, therefore, a net source of dissolved inorganic nitrogen(DIN = NO3- + NO2- + NH4+) during winter. This is in large part due to bioirrigation, which accounts for 76% of the benthic efflux of NH4+, thus reducing the capacity for nitrification of NH4+. The combined rate of fixed N loss by denitrification and anammox was estimated at 0.08 mmol m-2 d-1 of N2, which is at the lower end of previously reported values. A systematic sensitivity analysis revealed that denitrification and anammox respond strongly and positively to the concentration of NO3- in the bottomwater (NO3-BW).Higher O2-BW decreases DNRA and denitrification but stimulates both anammox and the contribution ofanammox to total N2 production (%Ramx). A complete mechanistic explanation of these findings is provided. Our analysis indicates that nitrification is the geochemical driving force behind the observed correlation between %Ramx and water depth in the seminal study of Dalsgaard et al. (2005). Despite remaining uncertainties, the results provide a general mechanistic framework for interpreting the existing knowledge of N-turnover processes and fluxes in continental margin sediments, as well as predicting the types of environment where these reactions are expected to occur prominently.

Benthic nitrogen (N) cycling was investigated at six stations along a transect traversing the Peruvian oxygen minimum zone (OMZ) at 11 °S. An extensive dataset including porewater concentration profiles and in situ benthic fluxes of nitrate (NO3–), nitrite (NO2–) and ammonium (NH4+) was used to constrain a 1–D reaction–transport model designed to simulate and interpret the measured data at each station. Simulated rates of nitrification, denitrification, anammox and dissimilatory nitrate reduction to ammonium (DNRA) by filamentous large sulfur bacteria (e.g. Beggiatoa and Thioploca) were highly variable throughout the OMZ yet clear trends were discernible. On the shelf and upper slope (80 – 260 m water depth) where extensive areas of bacterial mats were present, DNRA dominated total N turnover (less-than-or-equals, slant 2.9 mmol N m–2 d–1) and accounted for greater-or-equal, slanted 65 % of NO3– + NO2– uptake by the sediments from the bottom water. Nonetheless, these sediments did not represent a major sink for dissolved inorganic nitrogen (DIN = NO3– + NO2– + NH4+) since DNRA reduces NO3– and, potentially NO2–, to NH4+. Consequently, the shelf and upper slope sediments were recycling sites for DIN due to relatively low rates of denitrification and high rates of ammonium release from DNRA and ammonification of organic matter. This finding contrasts with the current opinion that sediments underlying OMZs are a strong sink for DIN. Only at greater water depths (300 – 1000 m) did the sediments become a net sink for DIN. Here, denitrification was the major process (less-than-or-equals, slant 2 mmol N m–2 d–1) and removed 55 – 73 % of NO3– and NO2– taken up by the sediments, with DNRA and anammox accounting for the remaining fraction. Anammox was of minor importance on the shelf and upper slope yet contributed up to 62 % to total N2 production at the 1000 m station. The results indicate that the partitioning of oxidized N (NO3–, NO2–) into DNRA or denitrification is a key factor determining the role of marine sediments as DIN sinks or recycling sites. Consequently, high measured benthic uptake rates of oxidized N within OMZs do not necessarily indicate a loss of fixed N from the marine environment.

This study provides new estimates for the global offshore methane hydrate inventory formed due to microbial CH4 production under Quaternary and Holocene boundary conditions. A multi-1D model for particular organic carbon (POC) degradation, gas hydrate formation and dissolution is presented. The novel reaction-transport model contains an open three-phase system of two solid compounds (organic carbon, gas hydrates), three dissolved species (methane, sulfates, inorganic carbon) and one gaseous phase (free methane). The model computes time-resolved concentration profiles for all compounds by accounting for chemical reactions as well as diffusive and advective transport processes. The reaction module builds upon a new kinetic model of POC degradation which considers a down-core decrease in reactivity of organic matter. Various chemical reactions such as organic carbon decay, anaerobic oxidation of methane, methanogenesis, and sulfate reduction are resolved using appropriate kinetic rate laws and constants. Gas hydrates and free gas form if the concentration of dissolved methane exceeds the pressure, temperature, and salinity-dependent solubility limits of hydrates and/or free gas, with a rate given by kinetic parameters. Global input grids have been compiled from a variety of oceanographic, geological and geophysical data sets including a new parameterization of sedimentation rates in terms of water depth. We find prominent gas hydrate provinces offshore Central America where sediments are rich in organic carbon and in the Arctic Ocean where low bottom water temperatures stabilize methane hydrates. The world’s total gas hydrate inventory is estimated at 0.82 x 10sup13 m3 - 2.10 x 10sup15 m3 CH4 (at STP conditions) or, equivalently, 4.18–995 Gt of methane carbon. The first value refers to present day conditions estimated using the relatively low Holocene sedimentation rates; the second value corresponds to a scenario of higher Quaternary sedimentation rates along continental margins. Our results clearly show that in-situ POC degradation is at present not an efficient hydrate forming process. Significant hydrate deposits in marine settings are more likely to have formed at times of higher sedimentation during the Quaternary or as a consequence of upward fluid transport at continental margins.

The question whether large scale glaciations on Antarctica were possible in a late Mesozoic greenhouse climate such as the Late Cretaceous is an intriguing one. The most recent years have provided an increasing number of studies investigating the growth and decay of paleo-continental ice sheets on Antarctica possibly large enough to affect sea level. Since the outcome of these studies doesn't provide a basis for a conclusive decision we have performed a number of model runs using an Atmospheric General Circulation Model (AGCM) to test whether large volumes of snow might have accumulated even under Late Cretaceous greenhouse conditions. By varying orbital parameters as well as topography, and atmospheric CO2 concentrations our models indicate the possibility of an Antarctic ice shield build-up large enough to drive sea level fluctuations on the order of tens of meters within ~ 20,000 years. This is supported under the assumption of pCO2 levels < 800 ppm, low insolation, and elevated topography. The growth of a major Antarctic ice sheet would be possible on reasonable time scales. To accumulate about half the present day snow/ice volume which is required to explain the documented shifts in oxygen isotopes our model results suggest a time span between 20,000 and 80,000 years for these ice volumes to accumulate.

International audience; Four seep sites located within an ∼20 km2 area offshore Georgia (Batumi seep area, Pechori Mound, Iberia Mound, and Colkheti Seep) show characteristic differences with respect to element concentrations, and oxygen, hydrogen, strontium, and chlorine isotope signatures in pore waters, as well as impregnation of sediments with petroleum and hydrocarbon potential. All seep sites have active gas seepage, near surface authigenic carbonates and gas hydrates. Cokheti Seep, Iberia Mound, and Pechori Mound are characterized by oil-stained sediments and gas seepage decoupled from deep fluid advection and bottom water intrusion induced by gas bubble release. Pechori Mound is further characterized by deep fluid advection of lower salinity pore fluids. The Pechori Mound pore fluids are altered by mineral/water reactions at elevated temperatures (between 60 and 110 °C) indicated by heavier oxygen and lighter chlorine isotope values, distinct Li and B enrichment, and K depletion. Strontium isotope ratios indicate that fluids originate from late Oligocene strata. This finding is supported by the occurrence of hydrocarbon impregnations within the sediments. Furthermore, light hydrocarbons and high molecular weight impregnates indicate a predominant thermogenic origin for the gas and oil at Pechori Mound, Iberia Mound, and Colkheti Seep. C15+ hydrocarbons at the oil seeps are allochtonous, whereas those at the Batumi seep area are autochthonous. The presence of oleanane, an angiosperm biomarker, suggests that the hydrocarbon source rocks belong to the Maikopian Formation. In summary, all investigated seep sites show a high hydrocarbon potential and hydrocarbons of Iberia Mound, Colkheti Seep, and Pechori Mound are predominantly of thermogenic origin. However, only at the latter seep site advection of deep pore fluids is indicated.

Cretaceous anoxic events may have been triggered by massive volcanic CO2 degassing as large igneous provinces (LIPs) were emplaced on the seafloor. Here, we present a comprehensive modeling study to decipher the marine biogeochemical consequences of enhanced volcanic CO2 emissions. A biogeochemical box model has been developed for transient model runs with time-dependent volcanic CO2 forcing. The box model considers continental weathering processes, marine export production, degradation processes in the water column, the rain of particles to the seafloor, benthic fluxes of dissolved species across the seabed, and burial of particulates in marine sediments. The ocean is represented by twenty-seven boxes. To estimate horizontal and vertical fluxes between boxes, a coupled ocean–atmosphere general circulation model (AOGCM) is run to derive the circulation patterns of the global ocean under Late Cretaceous boundary conditions. The AOGCM modeling predicts a strong thermohaline circulation and intense ventilation in the Late Cretaceous oceans under high pCO2 values. With an appropriate choice of parameter values such as the continental input of phosphorus, the model produces ocean anoxia at low to mid latitudes and changes in marine δ13C that are consistent with geological data such as the well established δ13C curve. The spread of anoxia is supported by an increase in riverine phosphorus fluxes under high pCO2 and a decrease in phosphorus burial efficiency in marine sediments under low oxygen conditions in ambient bottom waters. Here, we suggest that an additional mechanism might contribute to anoxia, an increase in the C:P ratio of marine plankton which is induced by high pCO2 values. According to our AOGCM model results, an intensively ventilated Cretaceous ocean turns anoxic only if the C:P ratio of marine organic particles exported into the deep ocean is allowed to increase under high pCO2 conditions. Being aware of the uncertainties such as diagenesis, this modeling study implies that potential changes in Redfield ratios might be a strong feedback mechanism to attain ocean anoxia via enhanced CO2 emissions. The formation of C-enriched marine organic matter may also explain the frequent occurrence of global anoxia during other geological periods characterized by high pCO2 values.

This study combines sediment geochemical analysis, in situ benthic lander deployments and numerical modeling to quantify the biogeochemical cycles of carbon and sulfur and the associated rates of Gibbs energy production at a novel methane seep. The benthic ecosystem is dominated by a dense population of tube-building ampharetid polychaetes and conspicuous microbial mats were unusually absent. A 1D numerical reaction-transport model, which allows for the explicit growth of sulfide and methane oxidizing microorganisms, was tuned to the geochemical data using a fluid advection velocity of 14 cm yr−1. The fluids provide a deep source of dissolved hydrogen sulfide and methane to the sediment with fluxes equal to 4.1 and 18.2 mmol m−2 d−1, respectively. Chemosynthetic biomass production in the subsurface sediment is estimated to be 2.8 mmol m−2 d−1 of C biomass. However, carbon and oxygen budgets indicate that chemosynthetic organisms living directly above or on the surface sediment have the potential to produce 12.3 mmol m−2 d−1 of C biomass. This autochthonous carbon source meets the ampharetid respiratory carbon demand of 23.2 mmol m−2 d−1 to within a factor of 2. By contrast, the contribution of photosynthetically-fixed carbon sources to ampharetid nutrition is minor (3.3 mmol m−2 d−1 of C). The data strongly suggest that mixing of labile autochthonous microbial detritus below the oxic layer sustains high measured rates of sulfate reduction in the uppermost 2 cm of the sulfidic sediment (100–200 nmol cm−3 d−1). Similar rates have been reported in the literature for other seeps, from which we conclude that autochthonous organic matter is an important substrate for sulfate reducing bacteria in these sediment layers. A system-scale energy budget based on the chemosynthetic reaction pathways reveals that up to 8.3 kJ m−2 d−1 or 96 mW m−2 of catabolic (Gibbs) energy is dissipated at the seep through oxidation reactions. The microorganisms mediating sulfide oxidation and anaerobic oxidation of methane (AOM) produce 95% and 2% of this energy flux, respectively. The low power output by AOM is due to strong bioenergetic constraints imposed on the reaction rate by the composition of the chemical environment. These constraints provide a high potential for dissolved methane efflux from the sediment (12.0 mmol m−2 d−1) and indicates a much lower efficiency of (dissolved) methane sequestration by AOM at seeps than considered previously. Nonetheless, AOM is able to consume a third of the ascending methane flux (5.9 mmol m−2 d−1 of CH4) with a high efficiency of energy expenditure (35 mmol CH4 kJ−1). It is further proposed that bioenergetic limitation of AOM provides an explanation for the non-zero sulfate concentrations below the AOM zone observed here and in other active and passive margin sediments.

A simple prognostic tool for gas hydrate (GH) quantification in marine sediments is presented based on a diagenetic transport-reaction model approach. One of the most crucial factors for the application of diagenetic models is the accurate formulation of microbial degradation rates of particulate organic carbon (POC) and the coupled formation of biogenic methane. Wallmann et al. (2006) suggested a kinetic formulation considering the ageing effects of POC and accumulation of reaction products (CH4, CO2) in the pore water. This model is applied to data sets of several ODP sites in order to test its general validity. Based on a thorough parameter analysis considering a wide range of environmental conditions, the POC accumulation rate (POCar in g/m2/yr) and the thickness of the gas hydrate stability zone (GHSZ in m) were identified as the most important and independent controls for biogenic GH formation. Hence, depth-integrated GH inventories in marine sediments (GHI in g of CH4 per cm2 seafloor area) can be estimated as: GHI = a · POCar · GHSZb · exp (– GHSZc/POCar/d) + e with a = 0.00214, b = 1.234, c = –3.339, d = 0.3148, e = –10.265. The transfer function gives a realistic first order approximation of the minimum GH inventory in low gas flux (LGF) systems. The overall advantage of the presented function is its simplicity compared to the application of complex numerical models, because only two easily accessible parameters need to be determined.

Two volcanic ash layer dominated marine gravity cores were recovered during research cruise RV Sonne 173/3 along the Pacific margin offshore Nicaragua. Thick ash layers dominate the recovered sediments. Pore water samples of these – one incoming plate and one slope site – sediment cores have been analysed for light rare earth elements, alkalinity, sulphate, phosphate, ammonia, calcium, and manganese contents. The data provide a systematic look at changes in REE during diagenesis proceeding from open ocean sediments to highly reducing near-shore sediments. The pore water patterns of all the dissolved species mentioned above suggest that volcanic glass alteration constitutes a major source for REE release in marine sediments when low rates of anaerobic degradation of organic matter, like typical for incoming plate/open ocean sediments, prevail. Under these conditions the release of REE results in fluxes of at least 0.04 nmol La cm− 2 yr− 1 and 0.09 nmol Ce cm− 2 yr− 1.

The sediment temperature distribution at mud volcanoes provides insights into their activity and into the occurrence of gas hydrates. If ambient pressure and temperature conditions are close to the limits of the gas hydrate stability field, the sediment temperature distribution not only limits the occurrence of gas hydrates, but is itself influenced by heat production and consumption related to the formation and dissociation of gas hydrates. Located in the Sorokin Trough in the northern Black Sea, the Dvurechenskii mud volcano (DMV) was in the focus of detailed investigations during the M72/2 and M73/3a cruises of the German R/V Meteor and the ROV Quest 4000 m in February and March 2007. A large number of in-situ sediment temperature measurements were conducted from the ROV and with a sensor-equipped gravity corer. Gas hydrates were sampled in pressurized cores using a dynamic autoclave piston corer (DAPC). The thermal structure of the DMV suggests a regime of fluid flow at rates decreasing from the summit towards the edges of the mud volcano, accompanied by intermittent mud expulsion at the summit. Modeled gas hydrate dissociation temperatures reveal that the gas hydrates at the DMV are very close to the stability limits. Changes in heat flow due to variable seepage rates probably do not result in changes in sediment temperature but are compensated by gas hydrate dissociation and formation. (C) 2009 Elsevier Ltd. All rights reserved.

Vodyanitskii mud volcano is located at a depth of about 2070 m in the Sorokin Trough, Black sea. It is a 500-m wide and 20-m high cone surrounded by a depression, which is typical of many mud volcanoes in the Black Sea. 75 kHz sidescan sonar show different generations of mud flows that include mud breccia, authigenic carbonates, and gas hydrates that were sampled by gravity coring. The fluids that flow through or erupt with the mud are enriched in chloride (up to similar to 650 mmol L-1 at similar to 150-cm sediment depth) suggesting a deep source, which is similar to the fluids of the close-by Dvurechenskii mud volcano. Direct observation with the remotely operated vehicle QUEST revealed gas bubbles emanating at two distinct sites at the crest of the mud volcano, which confirms earlier observations of bubble-induced hydroacoustic anomalies in echosounder records. The sediments at the main bubble emission site show a thermal anomaly with temperatures at similar to 60 cm sediment depth that were 0.9 degrees C warmer than the bottom water. Chemical and isotopic analyses of the emanated gas revealed that it consisted primarily of methane (99.8%) and was of microbial origin (delta D-CH4 = -170.8 parts per thousand (SMOW), delta C-13-CH4 = -61.0 parts per thousand (V-PDB), delta C-13-C2H6 = -44.0 parts per thousand (V-PDB)). The gas flux was estimated using the video observations of the ROV. Assuming that the flux is constant with time, about 0.9 +/- 0.5 x 106 mol of methane is released every year. This value is of the same order-of-magnitude as reported fluxes of dissolved methane released with pore water at other mud volcanoes. This suggests that bubble emanation is a significant pathway transporting methane from the sediments into the water column. (C) 2009 Elsevier Ltd. All rights reserved.

In 2000, the Carbon Dioxide in the Ocean working group of the North Pacific Marine Science Organization (PICES) performed an international experiment on the intercalibration of the measurements of the total alkalinity in seawater using certified reference materials (CRM). Taking part in this experiment, Russian specialists presented the method by Bruevich. The results of the intercalibration showed that the alkalinity values obtained by Bruevich’s method using modern burettes, an Na2CO3 reactant of high purity as a standard to ascertain the acid titre, and corrections for the acid density and for the weights of the acid and seawater samples in vacuum are in agreement with the standard within ±1 μ M/kg.

Porewater data from vent sites of the northeastern shelf off Sakhalin Island, Sea of Okhotsk, exhibit bottom-water concentrations down to a sediment depth of up to 300 cm. Below this depth, solute concentrations rapidly change due to methanogenesis and anaerobic methane oxidation (AMO). The profile shapes suggest an irrigation-like process that mixes on a meter scale. At these sites active gas emanation into the overlying water column and near-surface gas hydrates are commonly observed. We propose that methane gas bubbles rise through the soft surface sediments and cause mixing of the porewater. Mathematically, the bubble-induced irrigation can be described by eddy diffusion enhancing the diffusive transport of solutes by several orders of magnitude. A 3-D numerical transport-reaction model was developed to investigate the parameters defining the mixing process, such as bubble rise velocity, tube size, tube distribution in the sediment, and ebullition frequency. Model consistency with the field data requires eddy diffusivities ⩾1 × 105 cm2/a, tube densities of >4 tubes/m2 (equivalent to a tube spacing of

Future warming is predicted to shift the Earth system into a mode with progressive increase and vigour of extreme climate events possibly stimulating other mechanisms that invigorate global warming. This study provides new data and modelling investigating climatic consequences and biogeochemical feedbacks that happened in a warmer world not, vert, similar 112 Myr ago. Our study focuses on the Cretaceous Oceanic Anoxic Event (OAE) 1b and explores how the Earth system responded to a moderate not, vert, similar 25,000 yr lasting climate perturbation that is modelled to be less than 1 °C in global average temperature. Using a new chronological model for OAE 1b we present high-resolution elemental and bulk carbon isotope records from DSDP Site 545 from Mazagan Plateau off NW Africa and combine this information with a coupled atmosphere–land–ocean model. The simulations suggest that a perturbation at the onset of OAE 1b caused almost instantaneous warming of the atmosphere on the order of 0.3 °C followed by a longer (not, vert, similar 45,000 yr) period of not, vert, similar 0.8 °C cooling. The marine records from DSDP Site 545 support that these moderate swings in global climate had immediate consequences for African continental supply of mineral matter and nutrients (phosphorous), subsequent oxygen availability, and organic carbon burial in the eastern subtropical Atlantic, however, without turning the ocean anoxic. The match between modelling results and stratigraphic isotopic data support previous studies [summarized in Jenkyns, H.C., 2003. Evidence for rapid climate change in the Mesozoic–Palaeogene greenhouse world. The Royal Society, 361: 1885–1916.] in that methane emission from marine hydrates, albeit moderate in dimension, may have been the trigger for OAE 1b, though we can not finally rule out alternative mechanisms. Following the hydrate mechanism a total of 1.15 × 1018 g methane carbon (δ13C = − 60 ‰), equivalent to about 10% to the total modern gas hydrate inventory, generated the δ13Ccarb profile recorded in the section. Modelling suggests a combination of moderate-scale methane pulses supplemented by continuous methane emission at elevated levels over not, vert, similar 25,000 yr. The proposed mechanism, though difficult to finally confirm in the geological past, is arguably more likely to occur in a warmer world and apparently perturbs global climate and ocean chemistry almost instantaneously. This study shows that, once set-off, this mechanism can maintain Earth's climate in a perturbed mode over geological time leading to pronounced changes in regional climate.

Controversy over the oxygen isotope composition of seawater began in the 1950's, since which time there has been no agreement over whether the oxygen isotope composition of the oceans has changed over time. Resolving this uncertainty would allow the δ18O values of demonstrably well preserved marine authigenic precipitates to be used to reconstruct surface climate trends back to early Archean times and would help towards a more quantitative description of Earth's global water cycle on geological time scales. Isotopic studies of marine carbonate and silica reveal a trend of increasing δ18O values with decreasing age since the Archean. This trend has been interpreted by some to reflect a progressive increase in seawater δ18O through time; however, it is generally accepted on the basis of ophiolite studies and theoretical considerations that seawater δ18O cannot change significantly because of the buffering effects of ocean crust alteration at mid-ocean ridges. As a result many alternative interpretations have been proposed, including meteoric alteration; warmer paleoclimates; higher seawater pH; salinity stratification and biased sampling. Here we review these interpretations in the light of an updated compilation of marine carbonate δ18O from around the world, covering the Phanerozoic and Precambrian rock records. Recent models of the geological water cycle demonstrate how long-term trends in chemical weathering and hydrothermal circulation can indeed influence the O-isotope composition of the global ocean to the extent necessary to explain the carbonate δ18O trend, with residual variation attributed to climatic fluctuations and post-depositional alteration. We present the further development of an existing model of the geological water cycle. In the model, seawater δ18O increased from about − 13.3‰ to − 0.3‰ over a period of 3.4 Ga, with average surface temperatures fluctuating between 10 °C to 33 °C, which is consistent with known biological constraints. Similar temperature variations are also obtained, although with higher starting seawater δ18O composition, when more conservative approaches are used that take into account the systematic effects of diagenetic alteration on mean calcite δ18O values. In contrast to much published opinion, the average δ18O value of ocean crust in the model remained relatively unchanged throughout all model runs. Invariable ophiolite δ18O values can, therefore, not be used as a definitive argument against changing seawater δ18O. The most likely explanation for the long-term trend in seawater δ18O invokes two stepwise increases in the ratio of high- to low-temperature fluid/rock interactions. An initial increase may have occurred close to the Archean–Proterozoic boundary, but a possibly more significant increase took place near the Proterozoic–Phanerozoic boundary. Possible explanations for extremely low seawater δ18O during the Archean include higher continental weathering rates caused by a combination of higher atmospheric pCO2 (necessarily combined with high CO2 outgassing rates), a greater abundance of relatively easily weathered volcanic rocks in greenstone belts and partial emergence of spreading ridges. The second increase may have been caused by the suppression of low-temperature overprinting of ocean crust alteration by the formation of a thick sediment cover on ridge flanks due to the emergence of shelly plankton at the beginning of the Phanerozoic. Postulated increases in spreading ridge depths since the Archean would also have enhanced the efficiency of vertical heat flux and changed the depth at which hydrothermal fluids boil, both of which would favour high- over low-temperature interactions with time.

Pore fluids from the Green Canyon Block in the northern Gulf of Mexico show distinct differences with respect to element concentrations and oxygen, hydrogen, and strontium isotope signatures. The shallowest of the three investigated sites (GC185 or Bush Hill at 540 m water depth) is interpreted as a seafloor piercing mud mound and the two deeper areas (GC415 East and West at 950 and 1050 m water depth) as gas vent and oil seep sites. All three locations accommodate near-surface gas hydrates and the sediment surface is populated with chemosynthetic communities. They are characterized by a distinct increase in salinity with depth. However, the origin of this increasing salinity is different for the GC415 sites and Bush Hill and the depth source of the fluids is considerably different for all sites. The more saline fluids of the GC415 sites result from the dissolution of halite by formation water from two different sources. The fluids of GC415 East have most likely a deeper origin (early Cenozoic or even Mesozoic) and experienced elevated temperatures leading to mineral/water reactions including mineral transformations (e.g. smectite to illite transformation) and dissolution (e.g. feldspar dissolution). This process is expressed by the heavier oxygen isotope values and distinct Li, Sr, and Ca enrichments. The fluids of GC415 West have a shallower origin which is expressed by a smaller enrichment in Li, Sr, and Ca and lighter oxygen isotopes. The fluids from Bush Hill are less saline and its fluid signature indicates intensive water/mineral interaction. Oxygen and hydrogen isotope values as well as Na/Cl and Br/Cl molar ratios suggest that the salt enrichment was caused by phase separation under sub-critical conditions. A simple heat flow model simulation suggests that sub-critical phase separation may have occurred at a depth of ∼ 1650 m at ∼ 350 °C.

The Hg distribution and some mineralogical-geochemical features of bottom sediments up to a depth of 10 m in the Deryugin Basin showed that the high and anomalous Hg contents in the Holocene deposits are confined to a spreading riftogenic structure and separate fluid vents within it. The accumulations of Hg in the the sediments were caused by its fluxes from gas and low-temperature hydrothermal vents under favorable oceanological conditions in the Holocene. The two mainly responsible for the high and anomalous Hg contents are infiltration (fluxes of hydrothermal or gas fluids from the sedimentary cover) and plume (Hg precipitation from water plumes with certain hydrochemical conditions forming above endogenous sources). The infiltration anomalies of Hg were revealed in the following environments: (1) near gas vents on the northeastern Sakhalin slope, where high Hg contents are associated only with Se and were caused by the accumulation of gases ascending from beneath the gas hydrate layer; (2) in the area of inferred occasionally operating low-temperature hydrothermal seeps in the central part of the Deryugin Basin, in which massive barite chimneys, hydrothermal Fe-Mn crusts, and anomalous contents of Mn, Ba, Zn, and Ni in sediments develop.

We developed and employed an earth system model to explore the forcings of atmospheric pCO2 change and the chemical and isotopic evolution of seawater over the last glacial cycle. Concentrations of dissolved phosphorus, reactive nitrogen, molecular oxygen, dissolved inorganic carbon (DIC), total alkalinity (TA), 13C-DIC and 14C-DIC were calculated for 24 ocean boxes. The bi-directional water fluxes between these model boxes were derived from a 3-D circulation field of the modern ocean (Opa 8.2, NEMO) and tuned such that tracer distributions calculated by the box model were consistent with observational data from the modern ocean. To model the last 130 kyr, we employed records of past changes in sea-level, ocean circulation, and dust deposition. According to the model, about half of the glacial pCO2 drawdown may be attributed to marine regressions. The glacial sea-level low-stands implied steepened ocean margins, a reduced burial of particulate organic carbon, phosphorus, and neritic carbonate at the margin seafloor, a decline in benthic denitrification, and enhanced weathering of emerged shelf sediments. In turn, they led to a distinct rise in the standing stocks of DIC, TA, and nutrients in the global ocean, promoted the glacial sequestration of atmospheric CO2 in the ocean, and added 13C- and 14C-depleted DIC to the ocean as recorded in benthic foraminifera signals. The other half of the glacial drop in pCO2 was linked to reduced deep ocean dynamics, a shoaling of Atlantic meridional overturning circulation, and a rise in iron fertilization. The increased transit time of deep waters in the glacial ocean led to significant 14C depletions with respect to the atmosphere. The deglacial rapid and stepwise rise in atmospheric pCO2 was induced by upwelling both in the Southern Ocean and subarctic North Pacific and promoted by a drop in dust-borne iron discharge to the Southern Ocean. The deglacial sea-level rise led to a gradual decline in nutrient, DIC, and TA stocks, a slow change due to the large size and extended residence times of dissolved chemical species in the ocean. Thus, the rapid deglacial rise in pCO2 was dominated by fast changes in ocean dynamics and reduced dust deposition whereas the gradual pCO2 rise over the Holocene may be linked to the slow drop in nutrient and TA stocks that continued to promote an ongoing CO2 transfer from the ocean into the atmosphere.

As a result of extensive hydrocarbon exploration, the North Sea hosts several thousand abandoned wells; many believed to be leaking methane. However, how much of this greenhouse gas is emitted into the water column and ultimately reaches the atmosphere is not known. Here, we investigate three abandoned wells at 81–93 m water depth in the Norwegian sector of the North Sea, all of which show gas seepage into the bottom water. The isotopic signature of the emanating gas points towards a biogenic origin and hence to gas pockets in the sedimentary overburden above the gas reservoirs that the wells were drilled into. Video-analysis of the seeping gas bubbles and direct gas flow measurements resolved initial bubble sizes ranging between 3.2 and 7.4 mm in diameter with a total seabed gas flow between 1 and 19 tons of CH4 per year per well. Estimated total annual seabed emissions from all three wells of ∼24 tons are similar to the natural seepage rates at Tommeliten, suggesting that leaky abandoned wells represent a significant source of methane into North Sea bottom waters. However, the bubble-driven direct methane transport into the atmosphere was found to be negligible (