Intraseasonal variability of mixed layer depth in the tropical Indian Ocean

In this paper, we use an observational dataset built from Argo in situ profiles to describe the main large-scale patterns of intraseasonal mixed layer depth (MLD) variations in the Indian Ocean. An eddy permitting (0.25A degrees) regional ocean model that generally agrees well with those observed estimates is then used to investigate the mechanisms that drive MLD intraseasonal variations and to assess their potential impact on the related SST response. During summer, intraseasonal MLD variations in the Bay of Bengal and eastern equatorial Indian Ocean primarily respond to active/break convective phases of the summer monsoon. In the southern Arabian Sea, summer MLD variations are largely driven by seemingly-independent intraseasonal fluctuations of the Findlater jet intensity. During winter, the Madden-Julian Oscillation drives most of the intraseasonal MLD variability in the eastern equatorial Indian Ocean. Large winter MLD signals in northern Arabian Sea can, on the other hand, be related to advection of continental temperature anomalies from the northern end of the basin. In all the aforementioned regions, peak-to-peak MLD variations usually reach 10 m, but can exceed 20 m for the largest events. Buoyancy flux and wind stirring contribute to intraseasonal MLD fluctuations in roughly equal proportions, except for the Northern Arabian Sea in winter, where buoyancy fluxes dominate. A simple slab ocean analysis finally suggests that the impact of these MLD fluctuations on intraseasonal sea surface temperature variability is probably rather weak, because of the compensating effects of thermal capacity and sunlight penetration: a thin mixed-layer is more efficiently warmed at the surface by heat fluxes but loses more solar flux through its lower base.

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.

Mixed-layer salinity budget in the tropical Indian Ocean: seasonal cycle based only on observations

The mixed-layer salinity (MLS) budget in the tropical Indian Ocean is estimated from a combination of satellite products and in situ observations over the 2004-2012 period, to investigate the mechanisms controlling the seasonal MLS variability. In contrast with previous studies in the tropical Indian Ocean, our results reveal that the coverage, resolution, and quality of available observations are now sufficient to approach a closed monthly climatology seasonal salt budget. In the South-central Arabian Sea and South-western Tropical Indian Ocean (SCAS and STIO, respectively), where seasonal variability of the MLS is pronounced, the monthly MLS tendency terms are well captured by the diagnostic. In the SCAS region, in agreement with previous results, the seasonal cycle of the MLS is mainly due to meridional advection driven by the monsoon winds. In the STIO, contrasting previous results indicating the control of the meridional advection over the seasonal MLS budget, our results reveal the leading role of the freshwater flux due to 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

Optimizing PAT data transmission: assessing the accuracy of temperature summary data to estimate environmental conditions

Pop-up archival tags (PAT) provide summary and high-resolution time series data at predefined temporal intervals. The limited battery capabilities of PATs often restrict the transmission success and thus temporal coverage of both data products. While summary data are usually less affected by this problem, as a result of its lower size, it might be less informative. We here investigate the accuracy and feasibility of using temperature at depth summary data provided by PATs to describe encountered oceanographic conditions. Interpolated temperature at depth summary data was found to provide accurate estimates of three major thermal water column structure indicators: thermocline depth, stratification and ocean heat content. Such indicators are useful for the interpretation of the tagged animal's horizontal and vertical behaviour. The accuracy of these indicators was found to be particularly sensitive to the number of data points available in the first 100 m, which in turn depends on the vertical behaviour of the tagged animal. Based on our results, we recommend the use of temperature at depth summary data as opposed to temperature time series data for PAT studies; doing so during the tag programming will help to maximize the amount of transmitted time series data for other key data types such as light levels and depth.

Characterizing, modelling and understanding the climate variability of the deep water formation in the North-Western Mediterranean Sea

Observing, modelling and understanding the climate-scale variability of the deep water formation (DWF) in the North-Western Mediterranean Sea remains today very challenging. In this study, we first characterize the interannual variability of this phenomenon by a thorough reanalysis of observations in order to establish reference time series. These quantitative indicators include 31 observed years for the yearly maximum mixed layer depth over the period 1980–2013 and a detailed multi-indicator description of the period 2007–2013. Then a 1980–2013 hindcast simulation is performed with a fully-coupled regional climate system model including the high-resolution representation of the regional atmosphere, ocean, land-surface and rivers. The simulation reproduces quantitatively well the mean behaviour and the large interannual variability of the DWF phenomenon. The model shows convection deeper than 1000 m in 2/3 of the modelled winters, a mean DWF rate equal to 0.35 Sv with maximum values of 1.7 (resp. 1.6) Sv in 2013 (resp. 2005). Using the model results, the winter-integrated buoyancy loss over the Gulf of Lions is identified as the primary driving factor of the DWF interannual variability and explains, alone, around 50 % of its variance. It is itself explained by the occurrence of few stormy days during winter. At daily scale, the Atlantic ridge weather regime is identified as favourable to strong buoyancy losses and therefore DWF, whereas the positive phase of the North Atlantic oscillation is unfavourable. The driving role of the vertical stratification in autumn, a measure of the water column inhibition to mixing, has also been analyzed. Combining both driving factors allows to explain more than 70 % of the interannual variance of the phenomenon and in particular the occurrence of the five strongest convective years of the model (1981, 1999, 2005, 2009, 2013). The model simulates qualitatively well the trends in the deep waters (warming, saltening, increase in the dense water volume, increase in the bottom water density) despite an underestimation of the salinity and density trends. These deep trends come from a heat and salt accumulation during the 1980s and the 1990s in the surface and intermediate layers of the Gulf of Lions before being transferred stepwise towards the deep layers when very convective years occur in 1999 and later. The salinity increase in the near Atlantic Ocean surface layers seems to be the external forcing that finally leads to these deep trends. In the future, our results may allow to better understand the behaviour of the DWF phenomenon in Mediterranean Sea simulations in hindcast, forecast, reanalysis or future climate change scenario modes. The robustness of the obtained results must be however confirmed in multi-model studies.

Mean structure of the North Atlantic subtropical permanent pycnocline from in-situ observations

In the north Atlantic subtropical gyre, the oceanic vertical structure of density is characterized by a region of rapid increase with depth. This layer is called the permanent pycnocline. The permanent pycnocline is found below a surface mode water ,which are ventilated every winter when penetrated locally by the mixed layer. Assessing the structure and variability of the permanent pycnocline is of a major interest in the understanding of the climate system because the pycnocline layer delimits important heat and anthropogenic reservoir. Moreover, the heat content structure translate into changes in the large scale stratification feature, such as the permanent pycnocline. We developed a new objective algorithm for the characterization of the large scale structure of the permanent pycnocline (OAC-P). Argo data have been used with OAC-P to provide a detailed description of the mean structure of the North-Atlantic subtropical pycnocline (e.g.: depth, thickness, temperature, salinity, density, potential vorticity). Results reveal a surprisingly complex structure with inhomogeneous properties. While the classical bowl shape of the pycnocline depth is captured, much more complex pycnocline structure emerges at the regional scale. In the southern recirculation gyre of the Gulf Stream Extension, the pycnocline is deep, thick, the maximum of stratification is found in the middle on the layer and follow an isopycnal surface. But local processes influence and modify this textbook description and the pycnocline is characterized by a vertically asymmetric structure and gradients in thermohaline properties. T/S distribution along the permanent pycnocline depth is complex and reveals a diversity of water masses resulting from mixing of different source waters. We will present the observed mean structure of the North-Atlantic subtropical permanent pycnocline and relate it to physical processes that constraint it.

Dating of sediments in the Biscay bay: Implication for pollution chronology

Important historical informations on the temporal changes of anthropogenic pollution in marine environment can be assessed using sediment analysis. Dating is a crucial prerequisite to reconstruct pollution events, to calculate fluxes, and thus to allow comparison between different sites. This work presents estimates of accumulation rates of sediments in the Bay of Biscay. Fives cores were collected during RIKEAU 2002 cruise on board o/v Thalia in order to study temporal changes in PAH and organohalogens compounds content of sediment. We compare chronostratigraphic estimates on cores derived from the natural radionuclide 210Pb in excess with estimates from the known times of introduction of the artificial radionuclide 137Cs to the environment. 210Pb, 226Ra and 137Cs were measured directly by non-destructive gamma spectrometry using a well type γ-detector. Total 210Pb and 226Ra activities vary from 30 to 150 mBq g-1, and 20 to 36 mBq g-1 respectively; 137Cs presents lower levels (< 5 mBq g-1). Profiles of 210Pb in three cores present a well mixed layer, from 2-3 to 10 cm, in the uppermost sediments, followed by an exponential decrease of activities, suitable for the determination of sedimentation rates. Under constant flux and sedimentation rate assumptions, vertical accretion rates derived from 210Pb present a large range from nearly 0.1 cm yr-1 up to almost 0.3 cm yr-1. Differences are mainly due to relative position of studied cores regarding the muddy patch. Although the moderate level of 137Cs limits the accuracy of this dating method, profiles of 137Cs with depth strengthen mean rates derived from 210Pb data. The implication of this dating on pollutant inputs in sediments of the Bay of Biscay is briefly discussed.