Benthic silicon cycling in the Arctic Barents Sea: a reaction–transport model study

• Published in: Biogeosciences
• Main Authors: J. P. J. Ward, K. R. Hendry, S. Arndt, J. C. Faust, F. S. Freitas, S. F. Henley, J. W. Krause, C. März, A. C. Tessin, R. L. Airs
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2022
• Subjects: Ecology; QH540-549.5; Life; QH501-531; Geology; QE1-996.5; Arctic; Arctic Ocean; Barents Sea; Phytoplankton; Sea ice
• Online Access: https://doi.org/10.5194/bg-19-3445-2022; https://doaj.org/article/f40d67c44f3d461ea3b1a82a455fef10

Over recent decades the highest rates of water column warming and sea ice loss across the Arctic Ocean have been observed in the Barents Sea. These physical changes have resulted in rapid ecosystem adjustments, manifesting as a northward migration of temperate phytoplankton species at the expense of silica-based diatoms. These changes will potentially alter the composition of phytodetritus deposited at the seafloor, which acts as a biogeochemical reactor and is pivotal in the recycling of key nutrients, such as silicon (Si). To appreciate the sensitivity of the Barents Sea benthic system to the observed changes in surface primary production, there is a need to better understand this benthic–pelagic coupling. Stable Si isotopic compositions of sediment pore waters and the solid phase from three stations in the Barents Sea reveal a coupling of the iron (Fe) and Si cycles, the contemporaneous dissolution of lithogenic silicate minerals (LSi) alongside biogenic silica (BSi), and the potential for the reprecipitation of dissolved silicic acid (DSi) as authigenic clay minerals (AuSi). However, as reaction rates cannot be quantified from observational data alone, a mechanistic understanding of which factors control these processes is missing. Here, we employ reaction–transport modelling together with observational data to disentangle the reaction pathways controlling the cycling of Si within the seafloor. Processes such as the dissolution of BSi are active on multiple timescales, ranging from weeks to hundreds of years, which we are able to examine through steady state and transient model runs. Steady state simulations show that 60 % to 98 % of the sediment pore water DSi pool may be sourced from the dissolution of LSi, while the isotopic composition is also strongly influenced by the desorption of Si from metal oxides, most likely Fe (oxyhydr)oxides (FeSi), as they reductively dissolve. Further, our model simulations indicate that between 2.9 % and 37 % of the DSi released into sediment pore waters is subsequently ...