Bottom RedOx Model (BROM, v.1.0): a coupled benthic-pelagic model for simulation of seasonal anoxia and its impact
Abstract. Interaction between seawater and benthic sediments plays an important role in global biogeochemical cycling. Benthic fluxes of chemical elements (C, N, P, O, Si, Fe, Mn, S) directly affect redox state and acidification (i.e. pH and carbonate saturation), which in turn determine the functioning of the benthic and pelagic ecosystems. The redox state of the near bottom layer can change and oscillate in many regions responding to the supply of organic matter, physical regime and coastal discharge. The goal of this work was to develop a model that captures key biogeochemical processes occurring at the bottom boundary layer and sediment–water interface and analyze the changes that result from seasonal variability in redox conditions in the water column. We used a modular approach allowing the model to be coupled to existing hydrophysical models in 1-D, 2-D or 3-D. The model is capable to simulate seasonality in production and respiration of organic matter as well as in mixing, that leads to variation of redox conditions in the bottom boundary layer. Production and reduction of organic matter and varying redox conditions in the bottom boundary layer affect the carbonate system and lead to changes in pH and alkalinity. Bacteria play a significant role in the fate of organic matter due to chemosynthesis (autotrophs) and consumption of organic matter (heterotrophs). Changes in the bottom boundary layer redox conditions modify the distribution of nutrients (N and P) and redox metals (Mn and Fe). The model can be used for analyzing and interpreting data on sediment-water exchange, and estimating the consequences of forcing such as climate change, external nutrient loading, ocean acidification, carbon storage leakages, and point-source metal pollution.