As part of the French National Programme for Coastal Oceanography, this paper focuses on improvement of biogeochemical modelling in the: Bay of Seine (Eastern Channel), by introducing phosphorus in nutrient cycles. The Bay of Seine receives the Seine river, which exhibits very high nutrient concentrations, and this coastal zone constitutes a typical case of eutrophication in a river plume area. In terms of analyses, sequential extraction and analysis of sedimentary phosphorus were used in order, to measure various forms of particulate phosphorus in suspended matter and sediment (calcium-bound phosphate, Fe/Al-bound phosphate, exchangeable phosphate and organic phosphorus). In the modelling approach, the Bay is divided into 42 boxes and a two-layer, vertical thermohaline model is linked with the horizontal circulation scheme to take vertical stratification into account. The previous biological sub-model, with only nitrogen and silicon as nutrients, was improved by taking account of (a) various forms of bioavailable phosphorus (dissolved phosphate, exchangeable phosphate and organic phosphorus) and (b) experimental parameters, which govern the adsorption/desorption of phosphate at the solid-water interface. The ecological model accurately reproduces dissolved nutrient behaviour in the river plume, and concentrations of diatoms, flagellates and zooplankton are consistent with observed values. Simulated exchangeable phosphate and organic phosphorus in suspended matter are in the range of measured concentrations, whereas these particulate variables are simulated with less accuracy in the sediment. After evaluating the respective roles of nitrogen, phosphorus and silicon as limiting factors for phytoplanktonic growth, the model was used to estimate consequences of reduced phosphorus or nitrogen input for phytoplanktonic production in the Bay of Seine.