Mixed Layer formation and restratification in presence of mesoscale and submesoscale turbulence

Recent realistic high resolution modeling studies show a net increase of submesoscale activity in fall and winter when the mixed layer depth is at its maximum. This submesoscale activity increase is associated with a reduced deepening of the mixed layer. Both phenomena can be related to the development of mixed layer instabilities, which convert available potential energy into submesoscale eddy kinetic energy and contribute to a fast restratification by slumping the horizontal density gradient in the mixed layer. In the present work, the mixed layer formation and restratification was studied by uniformly cooling a fully turbulent zonal jet in a periodic channel at different resolutions, from eddy resolving (10 km) to submesoscale permitting (2 km). The effect of the submesoscale activity, highlighted by these different horizontal resolutions, was quantified in terms of mixed layer depth, restratification rate and buoyancy fluxes. Contrary to many idealized studies focusing on the restratification phase only, this study addresses a continuous event of mixed layer formation followed by its complete restratification. The robustness of the present results was established by ensemble simulations. The results show that, at higher resolution, when submesoscale starts to be resolved, the mixed layer formed during the surface cooling is significantly shallower and the total restratification almost three times faster. Such differences between coarse and fine resolution models are consistent with the submesoscale upward buoyancy flux, which balances the convection during the formation phase and accelerates the restratification once the surface cooling is stopped. This submesoscale buoyancy flux is active even below the mixed layer. Our simulations show that mesoscale dynamics also cause restratification, but on longer time scales. Finally, the spatial distribution of the mixed layer depth is highly heterogeneous in the presence of submesoscale activity, prompting the question of whether it is possible to parameterize submesoscale effects and their effects on the marine biology as a function of a spatially-averaged mixed layer depth.

Detection of mesoscale thermal fronts from 4km data using smoothing techniques: Gradient-based fronts classification and basin scale application

In order to optimize frontal detection in sea surface temperature fields at 4 km resolution, a combined statistical and expert-based approach is applied to test different spatial smoothing of the data prior to the detection process. Fronts are usually detected at 1 km resolution using the histogram-based, single image edge detection (SIED) algorithm developed by Cayula and Cornillon in 1992, with a standard preliminary smoothing using a median filter and a 3 × 3 pixel kernel. Here, detections are performed in three study regions (off Morocco, the Mozambique Channel, and north-western Australia) and across the Indian Ocean basin using the combination of multiple windows (CMW) method developed by Nieto, Demarcq and McClatchie in 2012 which improves on the original Cayula and Cornillon algorithm. Detections at 4 km and 1 km of resolution are compared. Fronts are divided in two intensity classes (“weak” and “strong”) according to their thermal gradient. A preliminary smoothing is applied prior to the detection using different convolutions: three type of filters (median, average and Gaussian) combined with four kernel sizes (3 × 3, 5 × 5, 7 × 7, and 9 × 9 pixels) and three detection window sizes (16 × 16, 24 × 24 and 32 × 32 pixels) to test the effect of these smoothing combinations on reducing the background noise of the data and therefore on improving the frontal detection. The performance of the combinations on 4 km data are evaluated using two criteria: detection efficiency and front length. We find that the optimal combination of preliminary smoothing parameters in enhancing detection efficiency and preserving front length includes a median filter, a 16 × 16 pixel window size, and a 5 × 5 pixel kernel for strong fronts and a 7 × 7 pixel kernel for weak fronts. Results show an improvement in detection performance (from largest to smallest window size) of 71% for strong fronts and 120% for weak fronts. Despite the small window used (16 × 16 pixels), the length of the fronts has been preserved relative to that found with 1 km data. This optimal preliminary smoothing and the CMW detection algorithm on 4 km sea surface temperature data are then used to describe the spatial distribution of the monthly frequencies of occurrence for both strong and weak fronts across the Indian Ocean basin. In general strong fronts are observed in coastal areas whereas weak fronts, with some seasonal exceptions, are mainly located in the open ocean. This study shows that adequate noise reduction done by a preliminary smoothing of the data considerably improves the frontal detection efficiency as well as the global quality of the results. Consequently, the use of 4 km data enables frontal detections similar to 1 km data (using a standard median 3 × 3 convolution) in terms of detectability, length and location. This method, using 4 km data is easily applicable to large regions or at the global scale with far less constraints of data manipulation and processing time relative to 1 km data.

Influence of the sea state on Mediterranean heavy precipitation: a case-study from HyMeX SOP1

Sea state can influence the turbulent air–sea exchanges, especially the momentum flux, by modifying the sea-surface roughness. The high-resolution non-hydrostatic convection-permitting model MESO-NH is used here to investigate the impact of a more realistic representation of the waves on heavy precipitation during the Intense Observation Period (IOP) 16a of the first HyMeX Special Observation Period (SOP1). Several quasi-stationary mesoscale convective systems developed over the western Mediterranean region, two of them over the sea, and resulted in heavy precipitation on the French and Italian coasts on 26 October 2012. Three different bulk parametrizations are tested in this study: a reference case (NOWAV) without any wave effect, a parametrization taking into account theoretical wave effects (WAV) and a last one with realistic wave characteristics from the MFWAM analyses (WAM). Using a realistic wave representation in WAM significantly increases the roughness length and the friction velocity with respect to NOWAV and WAV. The three MESO-NH sensitivity experiments of the IOP16a show that this surface-roughness increase in WAM generates higher momentum fluxes and directly impacts the low-level dynamics of the atmosphere, with a slowdown of the 10 m wind, when and where the wind speed exceeds 10 m s−1 and the sea state differs from the idealized one. The turbulent heat fluxes are not significantly influenced by the waves, these fluxes being controlled by the moisture content rather than by the wind speed in the simulations. Although the convective activity is globally well reproduced by all the simulations, the difference in the low-level dynamics of the atmosphere influences the localization of the simulated heavy precipitation. Objective evaluation of the daily rainfall amount and of the 10 m wind speed against the observations confirms the positive impact of the realistic wave representation on this simulation of heavy precipitation.

QuasiGeostrophic Diagnosis of Mixed-Layer Dynamics Embedded in a Mesoscale Turbulent Field

A quasigeostrophic model is developed to diagnose the three-dimensional circulation, including the vertical velocity, in the upper ocean from high-resolution observations of sea surface height and buoyancy. The formulation for the adiabatic component departs from the classical surface quasigeostrophic framework considered before since it takes into account the stratification within the surface mixed layer that is usually much weaker than that in the ocean interior. To achieve this, the model approximates the ocean with two constant stratification layers: a finite-thickness surface layer (or the mixed layer) and an infinitely deep interior layer. It is shown that the leading-order adiabatic circulation is entirely determined if both the surface streamfunction and buoyancy anomalies are considered. The surface layer further includes a diabatic dynamical contribution. Parameterization of diabatic vertical velocities is based on their restoring impacts of the thermal wind balance that is perturbed by turbulent vertical mixing of momentum and buoyancy. The model skill in reproducing the three-dimensional circulation in the upper ocean from surface data is checked against the output of a high-resolution primitive equation numerical simulation