Modelling turbulent vertical mixing sensitivity using a 1-D version of NEMO

Abstract

Abstract. Through two numerical experiments, a 1-D vertical model called NEMO1D was used to investigate physical and numerical turbulent-mixing behaviour. The results show that all the turbulent closures tested (k + l from Blanke and Delecluse, 1993, and two equation models: generic length scale closures from Umlauf and Burchard, 2003) are able to correctly reproduce the classical test of Kato and Phillips (1969) under favourable numerical conditions while some solutions may diverge depending on the degradation of the spatial and time discretization. The performances of turbulence models were then compared with data measured over a 1-year period (mid-2010 to mid-2011) at the PAPA station, located in the North Pacific Ocean. The modelled temperature and salinity were in good agreement with the observations, with a maximum temperature error between −2 and 2 • C during the stratified period (June to October). However, the results also depend on the numerical conditions. The vertical RMSE varied, for different turbulent closures, from 0.1 to 0.3 • C during the stratified period and from 0.03 to 0.15 • C during the homogeneous period. This 1-D configuration at the PAPA station (called PAPA1D) is now available in NEMO as a reference configuration including the input files and atmospheric forcing set described in this paper. Thus, all the results described can be recovered by downloading and launching PAPA1D. The configuration is described on the NEMO site (http://www.nemo-ocean.eu/Using-NEMO/ Configurations/C1D_PAPA). This package is a good starting point for further investigation of vertical processes.

Validation of an ocean shelf model for the prediction of mixed-layer properties in the Mediterranean Sea west of Sardinia

“…This is also in agreement with the findings of Reffray et al (2015). Hence, it was conjectured that the parameterisation of the vertical transport of heat and/or momentum was not adequate.…”

“…An indicator of why the k-ω parameterisation performed better than the other closure schemes is possibly found in the publication of Reffray et al (2015). Here, a one-dimensional model implemented in a three-dimensional circulation model was used to investigate physical and numerical turbulentmixing behaviour.…”

Do the Amazon and Orinoco freshwater plumes really matter for hurricane‐induced ocean surface cooling?

“…The vertical diffusion coefficients are given by a Generic Length Scale (GLS) scheme with a k-e turbulent closure [Umlauf and Burchard, 2003], complemented with the type ''A'' full equilibrium form of Canuto et al [2001] stability functions. A full description of the scheme is given in Reffray et al [2015] and the scheme has already proven its efficiency in regional configuration [e.g., Maraldi et al, 2013].…”

Dissipation of the energy imparted by mid-latitude storms in the Southern Ocean

“…A thirdorder upstream biased scheme (UP3) is used for both tracer and momentum advection, with no explicit diffusion. The vertical diffusion coefficients are given by a generic length scale (GLS) scheme with a k −ε turbulent closure (Reffray et al, 2015). Bottom friction is linear with a bottom drag coefficient of 1.5 × 10 −3 m s −1 .…”