The sensitivity of numerical simulation to vertical mixing parameterization schemes: a case study for the Yellow Sea Cold Water Mass

The vertical mixing parameterization scheme,by providing the eff ects of some explicitly missed physical processes and more importantly closing the energy budgets,is a critical model component and therefore imposes signifi cant impacts on model performance.The Yellow Sea Cold Water Mass(YSCWM),as the most striking and unique phenomenon in the Yellow Sea during summer,is dramatically aff ected by vertical mixing process during its each stage and therefore seriously sensitive to the proper choice of parameterization scheme.In this paper,a hindcast of YSCWM in winter of 2006 was implemented by using the Regional Ocean Modeling System(ROMS).Three popular parameterization schemes,including the level 2.5 Mellor-Yamada closure(M-Y 2.5),Generic Length Scale closure(GLS)and K-Profi le Parameterization(KPP),were tested and compared with each other by conducting a series of sensitivity model experiments.The infl uence of diff erent parameterization schemes on modeling the YSCWM was then carefully examined and assessed based on these model experiments.Although reasonable thermal structure and its seasonal variation were well reproduced by all schemes,considerable diff erences could still be found among all experiments.A warmer and spatially smaller simulation of YSCWM,with very strong thermocline,appeared in M-Y 2.5 experiment,while a spatially larger YSCWM with shallow mixed layer was found in GLS and KPP schemes.Among all the experiments,the discrepancy,indicated by core temperature,appeared since spring,and grew gradually by the end of November.Additional experiments also confi rmed that the increase of background diff usivity could eff ectively weaken the YSCWM,in either strength or coverage.Surface wave,another contributor in upper layer,was found responsible for the shrinkage of YSCWM coverage.The treatment of wave eff ect as an additional turbulence production term in prognostic equation was shown to be more superior to the strategy of directly increasing diff usivity for a coastal region.

The construction of high precision geostrophic currents based on new gravity models of GOCE and satellite altimetry data

The new gravity field models of gravity field and steady-state ocean circulation explorer(GOCE),TIM_R6 and DIR_R6,were released by the European Space Agency(ESA)in June 2019.The sixth generation of gravity models have the highest possible signal and lowest error levels compared with other GOCE-only gravity models,and the accuracy is significantly improved.This is an opportunity to build high precision geostrophic currents.The mean dynamic topography and geostrophic currents have been calculated by the 5th(TIM_R5 and DIR_R5),6th(TIM_R6 and DIR_R6)release of GOCE gravity field models and ITSG-Grace2018 of GRACE gravity field model in this study.By comparison with the drifter results,the optimal filtering lengths of them have been obtained(for DIR_R5,DIR_R6,TIM_R5 and TIM_R6 models are 1°and for ITSG-Grace2018 model is 1.1°).The filtered results show that the geostrophic currents obtained by the GOCE gravity field models can better reflect detailed characteristics of ocean currents.The total geostrophic speed based on the TIM_R6 model is similar to the result of the DIR_R6 model with standard deviation(STD)of 0.320 m/s and 0.321 m/s,respectively.The STD of the total velocities are 0.333 m/s and 0.325 m/s for DIR_R5 and TIM_R5.When compared with ITSG-Grace2018 results,the STD(0.344 m/s)of total geostrophic speeds is larger than GOCE results,and the accuracy of geostrophic currents obtained by ITSG-Grace2018 is lower.And the absolute errors are mainly distributed in the areas with faster speeds,such as the Antarctic circumpolar circulation,equatorial region,Kuroshio and Gulf Stream areas.After the remove-restore technique was applied to TIM_R6 MDT,the STD of total geostrophic speeds dropped to 0.162 m/s.