. Harima-nada sea suffers from oriogotrophication/nori damages due to diatom bloom • 4 type phytoplankton were modeled in the benthic-pelagic ecosystem model • The model reproduced the lifecycle and blooms of harmful diatoms in Harima-nada • The physical-chemical-biological causes of blooms were revealed by the model • This phytoplankton modelling approach enables to apply harmful bloom problems Harmful diatom (Coscinodiscus waillesi and Eucampia zodiacus) blooms frequently occur in the western coastal area of Japan from late autumn to early spring, causing discoloration of edible seaweed nori (Pyropia yezoensis). One of the causes of the discoloration is thought to be nutrient depletion (oligotrophication) due to nutrient uptake by diatoms. The goal of our study was to elucidate the mechanisms underlying phytoplankton blooms and the effects of such blooms on comprehensive nutrient cycles in the benthic-pelagic ecosystem. We developed a new ecosystem model, EMAGIN-L.P., simultaneously representing phytoplankton life cycle and the benthic-pelagic ecosystem. The life cycles of four types of phytoplankton—dinoflagellates (DFL), small diatoms (SDA), Coscinodiscus wailesii (CDW), and Eucampia zodiacus (EZ)—were numerically modeled and incorporated into an existing ecosystem model, ECOHYM (Sohma et al., 2008a), representing the food web among the benthic and pelagic systems and carbon-nitrogen-phosphorus-oxygen-coupled cycles driven by biogeochemical and physical processes. The model was applied to the Harima-Nada Sea where the discoloration of nori frequently occurs. The model reproduced the differences in timing and scale of the four phytoplankton blooms, and each bloom could be explained by a combination of temperature suitability, light availability, nutrient limitation, and vertical mixing within the water column. Furthermore, through a sensitivity analysis, competition for nutrient uptake among different types of phytoplankton was observed in the model. That is, from June to October, when DFL and SDA blooms occur, the SDA nutrient uptake rate is high and SDA suppresses the DFL bloom, while from November to March, when CDW and EZ blooms occur, the nutrient uptake rate of EZ is high and the CDW bloom is slightly suppressed. The developed model representing the phytoplankton bloom mechanisms in the whole benthic-pelagic coupled ecosystem will be a useful tool for revealing the mechanisms underlying the nutrient cycle and the contribution of harmful diatom blooms to the discoloration of nori. The modeling approach presented in this paper (i.e., categorizing different types of major phytoplankton and incorporating their life cycle into the entire benthic-pelagic coupling ecosystem) can be an effective method to understand, predict, and manage the problems caused by harmful diatom blooms across different damage types, harmful diatom species, and ocean areas. [ABSTRACT FROM AUTHOR]