Description: Methane concentrations in the water column and emissions to the atmosphere were determined for three tropical coastal lagoons surrounded by mangrove forests on the Yucatán Peninsula, Mexico. Surface water dissolved methane was sampled at different seasons over a period of 2 years in areas representing a wide range of salinities and anthropogenic impacts. The highest surface water methane concentrations (up to 8378 nM) were measured in a polluted canal associated with Terminos Lagoon. In Chelem Lagoon, methane concentrations were typically lower, except in the polluted harbor area (1796 nM). In the relatively pristine Celestún Lagoon, surface water methane concentrations ranged from 41 to 2551 nM. Methane concentrations were negatively correlated with salinity in Celestún, while in Chelem and Terminos high methane concentrations were associated with areas of known pollution inputs, irrespective of salinity. The diffusive methane flux from surface lagoon water to the atmosphere ranged from 0.0023 to 15 mmol CH4 m-2 d-1. Flux chamber measurements revealed that direct methane release as ebullition was up to 3 orders of magnitude greater than measured diffusive flux. Coastal mangrove lagoons may therefore be an important natural source of methane to the atmosphere despite their relatively high salinity. Pollution inputs are likely to substantially enhance this flux. Additional statistically rigorous data collected globally are needed to better consider methane fluxes from mangrove-surrounded coastal areas in response to sea level changes and anthropogenic pollution in order to refine projections of future atmospheric methane budgets.

Description: The benthic environment is a crucial component of marine systems in the provision of ecosystem services, sustaining biodiversity and in climate regulation, and therefore important to human society. With the contemporary increase in computational power, model resolution and technological improvements in quality and quantity of benthic data, it is necessary to ensure that benthic systems are appropriately represented in coupled benthic-pelagic biogeochemical and ecological modelling studies. In this paper we focus on five topical challenges related to various aspects of modelling benthic environments: organic matter reactivity, dynamics of benthic-pelagic boundary layer, microphytobenthos, biological transport and small-scale heterogeneity, and impacts of episodic events. We discuss current gaps in their understanding and indicate plausible ways ahead. Further, we propose a three-pronged approach for the advancement of benthic and benthic-pelagic modelling, essential for improved understanding, management and prediction of the marine environment. This includes: (A) development of a traceable and hierarchical framework for benthic-pelagic models, which will facilitate integration among models, reduce risk of bias, and clarify model limitations; (B) extended cross-disciplinary approach to promote effective collaboration between modelling and empirical scientists of various backgrounds and better involvement of stakeholders and end-users; (C) a common vocabulary for terminology used in benthic modelling, to promote model development and integration, and also to enhance mutual understanding.

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: 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.

Description: Takahe seep, located on the Opouawe Bank, Hikurangi Margin, is characterized by a well-defined subsurface seismic chimney structure ca. 80,500 m2 in area. Sub-seafloor geophysical data based on acoustic anomaly layers indicated the presence of gas hydrate and free gas layers within the chimney structure. Reaction-transport modeling was applied to porewater data from 11 gravity cores to constrain methane turnover rates and benthic methane fluxes in the upper 10 m. Model results show that methane dynamics were highly variable due to transport and dissolution of ascending gas. The dissolution of gas (up to 3761 mmol m−2 yr−1) dwarfed the rate of methanogenesis within the simulated sediment column (2.6 mmol m−2 yr−1). Dissolved methane is mainly consumed by anaerobic oxidation of methane (AOM) at the base of the sulfate reduction zone and trapped by methane hydrate formation below it, with maximum rates in the central part of the chimney (946 and 2420 mmol m−2 yr−1, respectively). A seep-wide methane budget was constrained by combining the biogeochemical model results with geophysical data and led to estimates of AOM rates, gas hydrate formation and benthic dissolved methane fluxes of 3.68 × 104 mol yr−1, 73.85 × 104 mol yr−1and 1.19 × 104 mol yr−1, respectively. A much larger flux of methane probably escapes in gaseous form through focused bubble vents. The approach of linking geochemical model results with spatial geophysical data put forward here can be applied elsewhere to improve benthic methane turnover rates from limited single spot measurements to larger spatial scales.

Description: The sediment-water interface is an important site for material exchange in marine systems and harbor unique microbial habitats. The flux of nutrients, metals, and greenhouse gases at this interface may be severely dampened by the activity of microorganisms and abiotic redox processes, leading to the “benthic filter” concept. In this study, we investigate the spatial variability, mechanisms and quantitative importance of a microbially-dominated benthic filter for dissolved sulfide in the Eastern Gotland Basin (Baltic Sea) that is located along a dynamic redox gradient between 65 and 173 m water depth. In August-September 2013, high resolution (0.25 mm minimum) vertical microprofiles of redox-sensitive species were measured in surface sediments with solid-state gold-amalgam voltammetric microelectrodes. The highest sulfide consumption (2.73–3.38 mmol m−2 day−1) occurred within the top 5 mm in sediments beneath a pelagic hypoxic transition zone (HTZ, 80–120 m water depth) covered by conspicuous white bacterial mats of genus Beggiatoa. A distinct voltammetric signal for polysulfides, a transient sulfur oxidation intermediate, was consistently observed within the mats. In sediments under anoxic waters (〉140 m depth), signals for Fe(II) and aqueous FeS appeared below a subsurface maximum in dissolved sulfide, indicating a Fe(II) flux originating from older sediments presumably deposited during the freshwater Ancylus Lake that preceded the modern Baltic Sea. Our results point to a dynamic benthic sulfur cycling in Gotland Basin where benthic sulfide accumulation is moderated by microbial sulfide oxidation at the sediment surface and FeS precipitation in deeper sediment layers. Upscaling our fluxes to the Baltic Proper; we find that up to 70% of the sulfide flux (2281 kton yr−1) toward the sediment-seawater interface in the entire basin can be consumed at the microbial mats under the HTZ (80–120 m water depth) while only about 30% the sulfide flux effuses to the bottom waters (〉120 m depth). This newly described benthic filter for the Gotland Basin must play a major role in limiting the accumulation of sulfide in and around the deep basins of the Baltic Sea.

Description: 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 (〈200 m) and 37 Gmol yr−1 by slope sediments (200–2000 m). The predicted deep-sea flux (〉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.

Description: Redox-sensitive mobilization of nutrients from sediments strongly affects the eutrophic state of the central Baltic Sea; a region associated with the spread of hypoxia and almost permanently anoxic and sulfidic conditions in the deeper basins. Ventilation of these basins depends on renewal by inflow of water enriched in oxygen (O2) from the North Sea, occurring roughly once per decade. Benthic fluxes and water column distributions of dissolved inorganic nitrogen species, phosphate (PO43-), dissolved inorganic carbon (DIC), sulfide (HS-), and total oxygen uptake (TOU) were measured along a depth gradient in the Eastern Gotland Basin (EGB). Campaigns were conducted during euxinic conditions of the deep basin in Aug./Sept. 2013 and after two inflow events in July/Aug. 2015 and March 2016 when O2 concentrations in deep waters reached 60 μM. The intrusion of O2-rich North Sea water into the EGB led to an approximate 33 and 10% reduction of the seabed PO43- and ammonium (NH4+) release from deep basin sediments. Post-inflow, the deep basin sediment was rapidly colonized by HS- oxidizing bacteria tentatively assigned to the family Beggiatoaceae, and HS- release was completely suppressed. The presence of a hypoxic transition zone (HTZ) between 80 and 120 m water depth was confirmed not only for euxinic deep-water conditions during 2013 but also for post-inflow conditions. Because deep-water renewal did not ventilate the HTZ, where PO43- and NH4+ fluxes were highest, high seabed nutrient release there was relatively unchanged. Extrapolation of the in situ nutrient fluxes indicated that, overall, the reduction in PO43- and NH4+ release in response to deep-water renewal can be considered as minor, reducing the internal nutrient load by 2 and 12% only, respectively. Infrequent inflow events thus have a limited capacity to sustainably reduce internal nutrient loading in the EGB and mitigate eutrophication.