The three-dimensional eco-hydrodynamic model was applied to estimate the autochthonous COD caused by production of phytoplankton in Jinhae Bay. A residual current was simulated, using a hydrodynamic model, to have a sightly complicated pattern in the inner part of the bay, ranging from 0.001 to 5 cm/s. In the outer part of the bay, the simulated current flowed out to the south sea with a southward flow at a maximum of 25 cm/s. The results of the ecological model simulation of COD levels showed high concentrations, exceeding 4 mg/L, in the inner bay of Masan, an area of wastewater discharge, and lower levels, approaching less than 1 mg/L, closer to the outer part of the bay. The simulation results of Autochthonous COD by two methods using ecological modeling, showed high ratio over 70% of total COD. Therefore, it is more important to consider nutrients than organic matters in the region for control COD standard.

The three-dimensional eco-hydrodynamic model was applied to estimate the physical process in terms of COD (chemical oxygen demand) and net supply(or decomposition) rate of COD in Kamak Bay to find proper management plan for oxygen demanding organic matters. The estimation results of the physical process in terms of COD showed that transportation of COD is dominant in surface level while accumulation of COD is dominant in bottom level. In the case of surface level, the net supply rate of COD was 0～0.50 mg/m2/day. The net decomposition rate of COD was 0～0.04 mg/m2/day in middle level(3～6m) and 0.05～0.15 mg/m2/day in bottom level(6m～bottom). These results indicates that the biological decomposition and physical accumulation of COD are occurred predominantly at the northern part of bottom level. Therefore, it is important to consider both allochthonous and autochthonous oxygen demanding organic matters in the region.

A three-dimensional ecological model (EMT-3D) was applied to Nonylphenol in Tokyo Bay. EMT-3D was calibrated with data obtained in the study area. The simulated results of dissolved Nonylphenol were in good agreement with the observed values, with a correlation coefficient(R) of 0.7707 and a coefficient of determination (R2) of 0.5940. The results of sensitivity analysis showed that biodegradation rate and bioconcentration factor are most important factors for dissolved Nonylphenol and Nonylphenol in phytoplankton, respectively. In the case of Nonylphenol in particulate organic carbon, biodegradation rate and partition coefficient were important factors. Therefore, the parameters must be carefully considered in the modeling. The mass balance results showed that standing stocks of Nonylphenol in water, in particulate organic carbon and in phytoplankton are 8.60×105 g, 2.19×102 g and 3.78×100 g, respectively. With respect to the flux of dissolved Nonylphenol, biodegradation in the water column, effluent to the open sea and partition to particulate organic carbon were 6.02×103 g/day, 6.02×102 g/day and 1.02×101 g/ day, respectively.