Abstract. Numerical models are a useful tool for studying marine ecosystems and associated biogeochemical processes in icecovered regions where observations are scarce. To this end, CSIB v1 (Canadian Sea-ice Biogeochemistry version 1), a new seaice biogeochemical model has been developed and embedded into the Nucleus for European Modelling of the Ocean (NEMO) modelling system. This model consists of a three-compartment (ice algae, nitrate, and ammonium) sea-ice ecosystem and a twocompartment (dimethylsulfoniopropionate and dimethylsulfide) sea-ice sulfur cycle which are coupled to pelagic ecosystem 5 and sulfur-cycle models at the sea ice-ocean interface. In addition to biological and chemical sources and sinks, the model simulates the horizontal transport of biogeochemical state variables within sea ice through a one-way coupling to a dynamicthermodynamic sea-ice model (LIM2). This paper describes technical aspects of implementing sea-ice biogeochemistry into NEMO and provides discussion on the results of several model experiments. Results of the reference simulation were evaluated by comparing the model outputs to observations and previous modelling studies. Additional simulations were conducted to 10 assess the model sensitivity to 1) the temporal resolution of the snowfall forcing data, 2) the representation of light penetration through snow, 3) advective and eddy-diffusive horizontal transport of sea-ice biogeochemical state variables, and 4) light attenuation by ice algae. The sea-ice biogeochemical model has been developed within the generic framework of NEMO to facilitate its use within different configurations and domains, and can be adapted for use with other NEMO-based submodels such as LIM3 and PISCES.