A direct determination of the World Ocean barotropic circulation

The time-mean Argo float displacements and the World Ocean Atlas 2009 temperature–salinity climatology are used to obtain the total, top to bottom, mass transports. Outside of an equatorial band, the total transports are the sum of the vertical integrals of geostrophic- and wind-driven Ekman currents. However, these transports are generally divergent, and to obtain a mass conserving circulation, a Poisson equation is solved for the streamfunction with Dirichlet boundary conditions at solid boundaries. The value of the streamfunction on islands is also part of the unknowns. This study presents and discusses an energetic circulation in three basins: the North Atlantic, the North Pacific, and the Southern Ocean. This global method leads to new estimations of the time-mean western Eulerian boundary current transports maxima of 97 Sverdrups (Sv; 1 Sv ≡ 106 m3 s−1) at 60°W for the Gulf Stream, 84 Sv at 157°E for the Kuroshio, 80 Sv for the Agulhas Current between 32° and 36°S, and finally 175 Sv for the Antarctic Circumpolar Current at Drake Passage. Although the large-scale structure and boundary of the interior gyres is well predicted by the Sverdrup relation, the transports derived from the wind stress curl are lower than the observed transports in the interior by roughly a factor of 2, suggesting an important contribution of the bottom torques. With additional Argo displacement data, the errors caused by the presence of remaining transient terms at the 1000-db reference level will continue to decrease, allowing this method to produce increasingly accurate results in the future.

A kinetic energy budget and internal instabilities in the Fine Resolution Antarctic Model

An energy analysis of the Fine Resolution Antarctic Model (FRAM) reveals the instability processes in the model. The main source of time-mean kinetic energy is the wind stress and the main sink is transfer to mean potential energy. The wind forcing thus helps maintain the density structure. Transient motions result from internal instabilities of the Bow rather than seasonal variations of the forcing. Baroclinic instability is found to be an important mechanism in FRAM. The highest values of available potential energy are found in the western boundary regions as well as in the Antarctic Circumpolar Current (ACC) region. All subregions with predominantly zonal flow are found to be baroclinically unstable. The observed deficit of eddy kinetic energy in FRAM occurs as a result of the high lateral friction, which decreases the growth rates of the most unstable waves. This high friction is required for the numerical stability of the model and can only be made smaller by using a finer horizontal resolution. A grid spacing of at least 10-15 km would be required to resolve the most unstable waves in the southern part of the domain. Barotropic instability is also found to be important for the total domain balance. The inverse transfer (that is, transfer from eddy to mean kinetic energy) does not occur anywhere, except in very localized tight jets in the ACC. The open boundary condition at the northern edge of the model domain does not represent a significant source or sink of eddy variability. However, a large exchange between internal and external mode energies is found to occur. It is still unclear how these boundary conditions affect the dynamics of adjacent regions.

How essential are Argo observations to constrain a global ocean data assimilation system?

Observing system experiments (OSEs) are carried out over a 1-year period to quantify the impact of Argo observations on the Mercator Ocean 0.25° global ocean analysis and forecasting system. The reference simulation assimilates sea surface temperature (SST), SSALTO/DUACS (Segment Sol multi-missions dALTimetrie, d'orbitographie et de localisation précise/Data unification and Altimeter combination system) altimeter data and Argo and other in situ observations from the Coriolis data center. Two other simulations are carried out where all Argo and half of the Argo data are withheld. Assimilating Argo observations has a significant impact on analyzed and forecast temperature and salinity fields at different depths. Without Argo data assimilation, large errors occur in analyzed fields as estimated from the differences when compared with in situ observations. For example, in the 0–300 m layer RMS (root mean square) differences between analyzed fields and observations reach 0.25 psu and 1.25 °C in the western boundary currents and 0.1 psu and 0.75 °C in the open ocean. The impact of the Argo data in reducing observation–model forecast differences is also significant from the surface down to a depth of 2000 m. Differences between in situ observations and forecast fields are thus reduced by 20 % in the upper layers and by up to 40 % at a depth of 2000 m when Argo data are assimilated. At depth, the most impacted regions in the global ocean are the Mediterranean outflow, the Gulf Stream region and the Labrador Sea. A significant degradation can be observed when only half of the data are assimilated. Therefore, Argo observations matter to constrain the model solution, even for an eddy-permitting model configuration. The impact of the Argo floats' data assimilation on other model variables is briefly assessed: the improvement of the fit to Argo profiles do not lead globally to unphysical corrections on the sea surface temperature and sea surface height. The main conclusion is that the performance of the Mercator Ocean 0.25° global data assimilation system is heavily dependent on the availability of Argo data.

Oceanic control of multidecadal variability in an idealized coupled GCM

Idealized ocean models are known to develop intrinsic multidecadal oscillations of the meridional overturning circulation (MOC). Here we explore the role of ocean–atmosphere interactions on this low-frequency variability. We use a coupled ocean–atmosphere model set up in a flat-bottom aquaplanet geometry with two meridional boundaries. The model is run at three different horizontal resolutions (4°, 2° and 1°) in both the ocean and atmosphere. At all resolutions, the MOC exhibits spontaneous variability on multidecadal timescales in the range 30–40 years, associated with the propagation of large-scale baroclinic Rossby waves across the Atlantic-like basin. The unstable region of growth of these waves through the long wave limit of baroclinic instability shifts from the eastern boundary at coarse resolution to the western boundary at higher resolution. Increasing the horizontal resolution enhances both intrinsic atmospheric variability and ocean–atmosphere interactions. In particular, the simulated atmospheric annular mode becomes significantly correlated to the MOC variability at 1° resolution. An ocean-only simulation conducted for this specific case underscores the disruptive but not essential influence of air–sea interactions on the low-frequency variability. This study demonstrates that an atmospheric annular mode leading MOC changes by about 2 years (as found at 1° resolution) does not imply that the low-frequency variability originates from air–sea interactions.

Comparative phylogeography of the western Indian Ocean reef fauna

Assessing patterns of connectivity at the community and population levels is relevant to marine resource management and conservation. The present study reviews this issue with a focus on the western Indian Ocean (WIO) biogeographic province. This part of the Indian Ocean holds more species than expected from current models of global reef fish species richness. In this study, checklists of reef fish species were examined to determine levels of endemism in each of 10 biogeographic provinces of the Indian Ocean. Results showed that the number of endemic species was higher in the WIO than in any other region of the Indian Ocean. Endemic species from the WIO on the average had a larger body size than elsewhere in the tropical Indian Ocean. This suggests an effect of peripheral speciation, as previously documented in the Hawaiian reef fish fauna, relative to other sites in the tropical western Pacific. To explore evolutionary dynamics of species across biogeographic provinces and infer mechanisms of speciation, we present and compare the results of phylogeographic surveys based on compilations of published and unpublished mitochondrial DNA sequences for 19 Indo-Pacific reef-associated fishes (rainbow grouper Cephalopholis argus, scrawled butterflyfish Chaetodon meyeri, bluespot mullet Crenimugil sp. A, humbug damselfish Dascyllus abudafur/Dascyllus aruanus, areolate grouper Epinephelus areolatus, blacktip grouper Epinephelus fasciatus, honeycomb grouper Epinephelus merra, bluespotted cornetfish Fistularia commersonii, cleaner wrasse Labroides sp. 1, longface emperor Lethrinus sp. A, bluestripe snapper Lutjanus kasmira, unicornfishes Naso brevirosris, Naso unicornis and Naso vlamingii, blue-spotted maskray Neotrygon kuhlii, largescale mullet Planiliza macrolepis, common parrotfish Scarus psicattus, crescent grunter Terapon jarbua, whitetip reef shark Triaenodon obesus) and three coastal Indo-West Pacific invertebrates (blue seastar Linckia laevigata, spiny lobster Panulirus homarus, small giant clam Tridacna maxima). Heterogeneous and often unbalanced sampling design, paucity of data in a number of cases, and among-species discrepancy in phylogeographic structure precluded any generalization regarding phylogeographic patterns. Nevertheless, the WIO might have been a source of haplotypes in some cases and it also harboured an endemic clade in at least one case. The present survey also highlighted likely cryptic species. This may eventually affect the accuracy of the current checklists of species, which form the basis of some of the recent advances in Indo-West Pacific marine ecology and biogeography.

The pianosa contourite depositional system (northern Tyrrhenian sea): drift morphology and plio-quaternary stratigraphic evolution

The Pianosa Contourite Depositional System (CDS) is located in the Corsica Trough (Northern Tyrrhenian Sea), a confined basin dominated by mass transport and contour currents in the eastern flank and by turbidity currents in the western flank. The morphologic and stratigraphic characterisation of the Pianosa CDS is based on multibeam bathymetry, seismic reflection data (multi-channel high resolution mini GI gun, single-channel sparker and CHIRP), sediment cores and ADCP data. The Pianosa CDS is located at shallow to intermediate water depths (170 to 850 m water depth) and is formed under the influence of the Levantine Intermediate Water (LIW). It is 120 km long, has a maximum width of 10 km and is composed of different types of muddy sediment drifts: plastered drift, separated mounded drift, sigmoid drift and multicrested drift. The reduced tectonic activity in the Corsica Trough since the early Pliocene permits to recover a sedimentary record of the contourite depositional system that is only influenced by climate fluctuations. Contourites started to develop in the Middle-Late Pliocene, but their growth was enhanced since the Middle Pleistocene Transition (0.7–0.9 Ma). Although the general circulation of the LIW, flowing northwards in the Corsica Trough, remained active all along the history of the system, contourite drift formation changed, controlled by sediment influx and bottom current velocity. During periods of sea level fall, fast bottom currents often eroded the drift crest in the middle and upper slope. At that time the proximity of the coast to the shelf edge favoured the formation of bioclastic sand deposits winnowed by bottom currents. Higher sediment accumulation of mud in the drifts occurred during periods of fast bottom currents and high sediment availability (i.e. high activity of turbidity currents), coincident with periods of sea level low-stands. Condensed sections were formed during sea level high-stands, when bottom currents were more sluggish and the turbidite system was disconnected, resulting in a lower sediment influx.

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.