A zero dimensional box model (PNCMjzb) with six state variables (ammonium, nitrate, dissolved organic nitrogen, phytoplankton, zooplankton and detritus) was developed to study nitrogen cycling in the Jiaozhou Bay pelagic ecosystem. The dominant processes within these compartments are considered with nitrogen as flow currency. Phytoplankton and zooplankton are treated as separate state variables, assuming that the species composition was dominated by two or three species the dynamic constants of which are similar and that they represent the entire plankton community. The microbial loop has not been integrated explicitly in the model. The turnover of bacteria is included implicitly in processes such as detritus decomposition, DON remineralization, pelagic nitrification and denitrification. The model is driven by two forcing variables, viz. water temperature and light intensity. Historical data from the1980s and 1990s were compiled and used for model calibration. In this paper (part I), the consideration of every main compartment in the model is interpreted in detail. And the applied equations and parameters are presented. The main results from the simulations together with discussion about phytoplankton dynamics and primary production in Jiaozhou Bay are presented in the next paper (part II). A zero dimensional box model (PNCMjzb) with six state variables (ammonium, nitrate, dissolved organic nitrogen, phytoplankton, zooplankton and detritus) was developed to study nitrogen cycling in the Jiaozhou Bay pelagic ecosystem. The dominant processes within these compartments are considered with nitrogen as flow currency. Phytoplankton and zooplankton are treated as separate state variables, assuming that the species composition was dominated by two or three species the dynamic constants of which are similar and that they represent the entire plankton community. The microbial loop has not been integrated explicitly in the model. The turnover of bacteria is included implicitly in processes such as detritus decomposition, DON remineralization, pelagic nitrification and denitrification. The model is driven by two forcing variables, viz. water temperature and light intensity. Historical data from the1980s and 1990s were compiled and used for model calibration. In this paper (part I), the consideration of every main compartment in the model is interpreted in detail. And the applied equations and parameters are presented. The main results from the simulations together with discussion about phytoplankton dynamics and primary production in Jiaozhou Bay are presented in the next paper (part II ).