Depth dependence of westward-propagating North Atlantic features diagnosed from altimetry and a numerical 1/6° model

International audienceA 1/6° numerical simulation is used to investigate the vertical structure of westward propagation between 1993 and 2000 in the North Atlantic ocean. The realism of the simulated westward propagating signals, interpreted principally as the signature of first-mode baroclinic Rossby waves (RW), is first assessed by comparing the simulated amplitude and zonal phase speeds of Sea Level Anomalies (SLA) against TOPEX/Poseidon-ERS satellite altimeter data. Then, the (unobserved) subsurface signature of RW phase speeds is investigated from model outputs by means of the Radon Transform which was specifically adapted to focus on first-mode baroclinic RW. The analysis is performed on observed and simulated SLA and along 9 simulated isopycnal displacements spanning the 0-3250 m depth range. Simulated RW phase speeds agree well with their observed counterparts at the surface, although with a slight slow bias. Below the surface, the simulated phase speeds exhibit a systematic deceleration with increasing depth, by a factor that appears to vary geographically. Thus, while the reduction factor is about 15-18% on average at 3250 m over the region considered, it appears to be much weaker (about 5-8%) in the eddy-active Azores Current, where westward propagating structures might be more coherent in the vertical. In the context of linear theories, these results question the often-made normal mode assumption of many WKB-based theories that the phase speed is independent of depth. Alternatively, these results could also suggest that the vertical structure of westward propagating signals may significantly depend on their degree of nonlinearity, with the degree of vertical coherence possibly increasing with the degree of nonlinearity

ASSESSING THE REALISM OF OCEAN SIMULATIONS AGAINST HYDROGRAPHY AND ALTIMETRY

Future climate prediction systems will include ocean models at eddy-admitting to eddy-resolving resolution, i.e. ° on the horizontal or finer. The development and calibration of such models requires the us

Can we map the interannual variability of the whole upper Southern Ocean with the current database of hydrographic observations?

International audienceWith the advent of Argo floats, it now seems feasible to study the interannual variations of upper ocean hydrographic properties of the historically undersampled Southern Ocean. To do so, scattered hydrographic profiles often first need to be mapped. To investigate biases and errors associated both with the limited space-time distribution of the profiles and with the mapping methods, we colocate the mixed-layer depth (MLD) output from a state-of-the-art 1/12° DRAKKAR simulation onto the latitude, longitude, and date of actual in situ profiles from 2005 to 2014. We compare the results obtained after remapping using a nearest neighbor (NN) interpolation and an objective analysis (OA) with different spatiotemporal grid resolutions and decorrelation scales. NN is improved with a coarser resolution. OA performs best with low decorrelation scales, avoiding too strong a smoothing, but returns values over larger areas with large decorrelation scales and low temporal resolution, as more points are available. For all resolutions OA represents better the annual extreme values than NN. Both methods underestimate the seasonal cycle in MLD. MLD biases are lower than 10 m on average but can exceed 250 m locally in winter. We argue that current Argo data should not be mapped to infer decadal trends in MLD, as all methods are unable to reproduce existing trends without creating unrealistic extra ones. We also show that regions of the subtropical Atlantic, Indian, and Pacific Oceans, and the whole ice-covered Southern Ocean, still cannot be mapped even by the best method because of the lack of observational data

Intrinsic variability of the Atlantic Meridional Overturning Circulation at interannual-to-multidecadal time scales

The low-frequency variability of the Atlantic meridional overturning circulation (AMOC) is investigated from 2, ¼°, and ° global ocean–sea ice simulations, with a specific focus on its internally generated (i.e., “intrinsic”) component. A 327-yr climatological ¼° simulation, driven by a repeated seasonal cycle (i.e., a forcing devoid of interannual time scales), is shown to spontaneously generate a significant fraction R of the interannual-to-decadal AMOC variance obtained in a 50-yr “fully forced” hindcast (with reanalyzed atmospheric forcing including interannual time scales). This intrinsic variance fraction R slightly depends on whether AMOCs are computed in geopotential or density coordinates, and on the period considered in the climatological simulation, but the following features are quite robust when mesoscale eddies are simulated (at both ¼° and ° resolutions); R barely exceeds 5%–10% in the subpolar gyre but reaches 30%–50% at 34°S, up to 20%–40% near 25°N, and 40%–60% near the Gulf Stream. About 25% of the meridional heat transport interannual variability is attributed to intrinsic processes at 34°S and near the Gulf Stream. Fourier and wavelet spectra, built from the 327-yr ¼° climatological simulation, further indicate that spectral peaks of intrinsic AMOC variability (i) are found at specific frequencies ranging from interannual to multidecadal, (ii) often extend over the whole meridional scale of gyres, (iii) stochastically change throughout these 327 yr, and (iv) sometimes match the spectral peaks found in the fully forced hindcast in the North Atlantic. Intrinsic AMOC variability is also detected at multidecadal time scales, with a marked meridional coherence between 35°S and 25°N (15–30 yr periods) and throughout the whole basin (50–90-yr periods)

Stochastic variability of oceanic flows above topography anomalies

International audienceWe describe a stochastic variability mechanism which is genuinely internal to the ocean, i.e. not due to fluctuations in atmospheric forcing. % The key ingredient is the existence of closed contours of bottom topography surrounded by a stirring region of enhanced eddy activity. This configuration leads to the formation of a robust but highly variable vortex above the topography anomaly. The vortex dynamics integrates the white noise forcing of oceanic eddies into a red noise signal for the large scale volume transport of the vortex. The strong interannual fluctuations of the transport of the Zapiola anticyclone ($\sim 100 \ Sv$) in the Argentine basin are argued to be partly due to such eddy-driven stochastic variability, on the basis of a $310$ years long simulation of a comprehensive global ocean model run driven by a repeated-year forcing

Oceanic hindcast simulations at high resolution suggest that the Atlantic MOC is bistable

All climate models predict a freshening of the North Atlantic at high latitude that may induce an abrupt change of the Atlantic Meridional Overturning Circulation (hereafter AMOC) if it resides in the bistable regime, where both a strong and a weak state coexist. The latter remains uncertain as there is no consensus among observations and ocean reanalyses, where the AMOC is bistable, versus most climate models that reproduce a mono-stable strong AMOC. A series of four hindcast simulations of the global ocean at 1/12° resolution, which is presently unique, are used to diagnose freshwater transport by the AMOC in the South Atlantic, an indicator of AMOC bistability. In all simulations, the AMOC resides in the bistable regime: it exports freshwater southward in the South Atlantic, implying a positive salt advection feedback that would act to amplify a decreasing trend in subarctic deep water formation as projected in climate scenarios

Variability and trends in Southern Ocean eddy activity in 1/12° ocean model simulations

The response of eddy kinetic energy (EKE) to the strengthening of Southern Hemisphere winds occurring since the 1950s is investigated with a global ocean model having a resolution of 1/12° in the Antarctic Circumpolar Current domain. The simulations expose regional differences in the relative importance of stochastic and wind-related contributions to inter-annual EKE changes. In the Pacific and Indian sectors the model captures the EKE variability observed since 1993 and confirms previous hypotheses of a lagged response to regional wind stress anomalies. Here, the multi-decadal trend in wind stress is reflected in an increase in EKE typically exceeding 5 cm2 sec-2 decade-1. In the western Atlantic EKE variability is mostly stochastic, is weakly correlated with wind fluctuations, and its multi-decadal trends are close to zero. The non-uniform distribution of wind-related changes in the eddy activity could affect the regional patterns of ocean circulation and biogeochemical responses to future climate change

The North Atlantic subpolar gyre in four high resolution models

The authors present the first quantitative comparison between new velocity datasets and high-resolution models in the North Atlantic subpolar gyre [1/10° Parallel Ocean Program model (POPNA10), Miami Isopycnic Coordinate Ocean Model (MICOM), ° Atlantic model (ATL6), and Family of Linked Atlantic Ocean Model Experiments (FLAME)]. At the surface, the model velocities agree generally well with World Ocean Circulation Experiment (WOCE) drifter data. Two noticeable exceptions are the weakness of the East Greenland coastal current in models and the presence in the surface layers of a strong southwestward East Reykjanes Ridge Current. At depths, the most prominent feature of the circulation is the boundary current following the continental slope. In this narrow flow, it is found that gridded float datasets cannot be used for a quantitative comparison with models. The models have very different patterns of deep convection, and it is suggested that this could be related to the differences in their barotropic transport at Cape Farewell. Models show a large drift in watermass properties with a salinization of the Labrador Sea Water. The authors believe that the main cause is related to horizontal transports of salt because models with different forcing and vertical mixing share the same salinization problem. A remarkable feature of the model solutions is the large westward transport over Reykjanes Ridge [10 Sv (Sv ≡ 106 m3 s−1) or more