“…The BFM was recently developed with a carbonate system (Cossarini et al, 2015) and revised phytoplankton nutrient uptake processes (Lazzari et al, 2016). Moreover, the BFM has been implemented in the Mediterranean Sea in a wide range of applications, including 1D and 3D configurations aimed at studying the interaction between optics and biogeochemistry (Terzić et al, 2021(Terzić et al, , 2019, conducting climate scenario simulations (Lazzari et al, 2014), estimating carbon sequestration services (Melaku Canu et al, 2015), simulating high-resolution coastal dynamics in marginal seas (Cossarini et al, 2017), analysing temporal scales of multi-decadal variability (Di Biagio et al, 2019), and studying CDOM (chromophoric-dissolved organic matter) spatiotemporal variability (Lazzari et al, 2021). 3DVar-Bio is the data assimilation scheme for the correction of phytoplankton biomass and nutrient concentrations (i.e.…”