Marine biogeochemical processes can strongly interact with processes occurring in adjacent ice and sediments. This is especially likely in areas with shallow water and frequent ice cover, both of which are common in the Arctic. Modeling tools are therefore required to simulate coupled biogeochemical systems in ice, water, and sediment domains. We developed a 1D sympagic-pelagic-benthic transport model (SPBM) which uses input from physical model simulations to describe hydrodynamics and ice growth and modules from the Framework for Aquatic Biogeochemical Models (FABM) to construct a user-defined biogeochemical model. SPBM coupled with a biogeochemical model simulates the processes of vertical diffusion, sinking/burial, and biogeochemical transformations within and between the three domains. The potential utility of SPBM is demonstrated herein with two test runs using modules from the European regional seas ecosystem model (ERSEM) and the bottom-redox model biogeochemistry (BROM-biogeochemistry). The first run simulates multiple phytoplankton functional groups inhabiting the ice and water domains, while the second simulates detailed redox biogeochemistry in the ice, water, and sediments. SPBM is a flexible tool for integrated simulation of ice, water, and sediment biogeochemistry, and as such may help in producing well-parameterized biogeochemical models for regions with strong sympagic-pelagic-benthic interactions.A biogeochemical model suitable for the Arctic should take into account the specific conditions of this region, such as the seasonal to permanent ice cover and the presence of shelf areas. Thus, the model should preferably combine processes occurring in three domains: ice, water column, and sediments. Each of these domains has some specific features and modeling challenges:Ice. The Arctic ice-algal primary production is a significant part of the total primary production of the Arctic region [7]. Photosynthetic microorganisms extend the production season, provide a winter and early spring food source, and contribute to organic carbon export to depth [8]. A modeling study [9] estimated an average Arctic ice-algal primary production of 21.7 Tg C year −1 , which equates to roughly 5% of total pelagic primary production [10] for this area. Other authors [7] estimated sea ice-algal production accounting for 5-10% of total Arctic and Southern Ocean primary productivity. Another modeling study [11] suggested that under a mild climate change scenario the sea ice community around Greenland may become generally more productive while pelagic phytoplankton productivity may decrease. It is therefore desirable to include the ice domain in biogeochemical modeling studies of the Arctic region. There are three main approaches to implement ice algae behavior according to the place where algae live in the ice column [12,13]: in the bottom layer of an ice column with fixed thickness, in the bottom layer of an ice column with variable thickness, or in any layer of an ice column. Recent research suggests that ice-algal models ...