The aim of this research is to develop the parameterization of subgrid-scale oceanic processes using the results of a fully three-dimensional large eddy simulation. The advanced parameterization scheme will be incorporated into an ocean general circulation model or a coupled atmosphere-ocean general circulation model to project the interannual variability in ocean mixed layer processes and other processes. By resolving the uncertainties of conventional parameterization schemes, we can expect significant improvement in the ability to model future climate change.