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.

The three-dimensional eco-hydrodynamic model was applied to estimate the physical process in terms of nutrients and net uptake(or regeneration) rate of nutrients in Kamak Bay for scenario analysis to find proper management plan. The estimation results of the physical process in terms of nutrients showed that transportation of nutrients is dominant in surface level while accumulation of nutrients is dominant in bottom level. In the case of dissolved inorganic nitrogen, the results showed that the net uptake rate was 0～60 mg/m2/day in surface level(0～3m), and the net regeneration rate was 0.0～10.0 mg/m2/day in middle level(3～6m) and above 10 mg/m2/day in bottom level(6m～below). In the case of dissolved inorganic phosphorus, the net uptake rate was 0.0～3.0 mg/m2/day in surface level, and the net regeneration rate was 0.5～1.5 mg/m2/day in middle level and 1.0～3.0 mg/m2/day in bottom level. These results indicates that net uptake and transport of nutrients are occurred predominantly at the surface level and the net generation and accumulation are dominant at bottom level. Therefore, it is important to consider the re-supplement of nutrients due to regeneration of bottom water.

It is very important to interprete and simulate the variation of phytoplankton maximum region for the prediction and control of red tide. This study was composed of two parts, first, the hydrodynamic simulation such as residual current and salinity diffusion, and second, the ecological simulation such as phytoplankton distribution according to freshwater discharge and pollutant loads. Without the Nakdong river discharge, residual current was stagnated in inner side of this estuary, and surface distribution of salinity was over 25psu. On the contrary, with summer mean discharge, freshwater stretched very far outward and some waters flowed into Chinhae Bay through the Kadok channel, and low salinity extended over coastal sea and salinity front occurred. From the result of contributed physical process to phytoplankton biomass, the accumulation was occurred at the west part of this estuary and the Kadok channel with the Nakdong river discharge. When more increased input discharge, the accumulation band was transported to outer side of this estuary. The frequently outbreak of red tide in this area is caused by accumulation of physical processes. The phytoplankton maximum region located inner side of this estuary without the Nakdong river discharge and with mean discharge of winter, but it was moved to outer side when mean discharge of the Nakdong river was increased. The variation of input concentration from the land loads was not largely influenced on phytoplankton biomass and location of maximum region. When discharge was increased, phytoplankton maximum region was transferred to inner side of the Kadok channel. On the other hand, when discharge was decreased, phytoplankton maximum region was transferred to inner side of this estuary and chlorophyll a contents increased to over 20㎍/L. Therefore, if any other conditions are favorable for growth of phytoplankton, decrease of discharge causes to increase of possibility of red tide outbreak.

It is noted that the red tides and the oxygen-deficient water mass are extensively developed in Masan Bay during summer. The nutrients mass balance was calculated in Masan Bay, using the three-dimensional numerical hydrodynamic model and the material cycle model. The material cycle model was calibrated with the data obtained on the field of the study area in June 1993. The nutrients mass balance calculated by the combination of the residual currents and material cycle model results showed nutrients of surface and middle levels to be transported from the inner part to the outer part of Masan Bay, and nutrients of bottom level to be transported from outer part to inner part of Masan Bay. The uptake rate of DIN in the box A1(surface level of inner part) was found to be 337.5㎎/㎥ ·day, the largest value in all 9 boxes and that of DIP was found to be 18.6㎎/㎥·day in box A1, and the regeneration rate of DIN was found to be 78.2㎎/㎥· day in the box A3(bottom level of inner part), and that of DIP was found to be 18.6㎎/㎥· day in box A1. The regenerations of DIN and DIP in the water column of the entire Bay were found to be 7.66ton/day and 760㎏/day, respectively. And the releases of DIN and DIP from the sediments of the entire Bay were found to be 2.86ton/day and 634㎏/ day, respectively. The regeneration rate was 2.5 times as high as the release rate in DIN, and 1.2 times in DIP. The results of mass balance calculation showed not only the nutrients released from the sediments but the nutrients regenerated in water column to be important in the control and management of water quality in Masan Bay.