Experimental and diagnostic protocol for the physical component of the CMIP6 Ocean Model Intercomparison Project (OMIP)

The Ocean Model Intercomparison Project (OMIP) aims to provide a framework for evaluating, understanding, and improving the ocean and sea-ice components of global climate and earth system models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses these aims in two complementary manners: (A) by providing an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing, (B) by providing a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) offering details for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows that of the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II have become the standard method to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP (Scenario MIP), as well as the ocean-sea ice OMIP simulations. The bulk of this paper offers scientific rationale for saving these diagnostics.

OMIP contribution to CMIP6: experimental and diagnostic protocol for the physical component of the Ocean Model Intercomparison Project

The Ocean Model Intercomparison Project (OMIP) is an endorsed project in the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses CMIP6 science questions, investigating the origins and consequences of systematic model biases. It does so by providing a framework for evaluating (including assessment of systematic biases), understanding, and improving ocean, sea-ice, tracer, and biogeochemical components of climate and earth system models contributing to CMIP6. Among the WCRP Grand Challenges in climate science (GCs), OMIP primarily contributes to the regional sea level change and near-term (climate/decadal) prediction GCs. OMIP provides (a) an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing; and (b) a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) detailing methods for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II (Interannual Forcing) have become the standard methods to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP, HighResMIP (High Resolution MIP), as well as the ocean/sea-ice OMIP simulations.

In-situ based reanalysis of the global ocean temperature and salinity with ISAS: variability of the heat content and steric height

The In Situ Analysis System (ISAS) was developed to produce gridded fields of temperature and salinity that preserve as much as possible the time and space sampling capabilities of the Argo network of profiling floats. Since the first global re-analysis performed in 2009, the system has evolved and a careful delayed mode processing of the 2002-2012 dataset has been carried out using version 6 of ISAS and updating the statistics to produce the ISAS13 analysis. This last version is now implemented as the operational analysis tool at the Coriolis data centre. The robustness of the results with respect to the system evolution is explored through global quantities of climatological interest: the Ocean Heat Content and the Steric Height. Estimates of errors consistent with the methodology are computed. This study shows that building reliable statistics on the fields is fundamental to improve the monthly estimates and to determine the absolute error bars. The new mean fields and variances deduced from the ISAS13 re-analysis and dataset show significant changes relative to the previous ISAS estimates, in particular in the southern ocean, justifying the iterative procedure. During the decade covered by Argo, the intermediate waters appear warmer and saltier in the North Atlantic and fresher in the Southern Ocean than in WOA05 long term mean. At inter-annual scale, the impact of ENSO on the Ocean Heat Content and Steric Height is observed during the 2006-2007 and 2009-2010 events captured by the network.

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.

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.

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.

Inertial dissipation method applied to derive turbulent fluxed over the ocean during the Surface of the Ocean, Fluxes and Interactions with the Atmosphere Atlantic Stratocumulus Transition Experiment (SOFIA/ASTEX) and Structure des Echanges Mer-Atmosphere, Proprietes des Heterogeneites Oceaniques: Recherche Experimentale (SEMAPHORE) experiments with low to moderate wind speed

The transfer coefficients for momentum and heat have been determined for 10 m neutral wind speeds (U-10n) between 0 and 12 m/s using data from the Surface of the Ocean, Fluxes and Interactions with the Atmosphere (SOFIA) and Structure des Echanges Mer-Atmosphere, Proprietes des Heterogeneites Oceaniques: Recherche Experimentale (SEMAPHORE) experiments. The inertial dissipation method was applied to wind and pseudo virtual temperature spectra from a sonic anemometer, mounted on a platform (ship) which was moving through the turbulence held. Under unstable conditions the assumptions concerning the turbulent kinetic energy (TKE) budget appeared incorrect. Using a bulk estimate for the stability parameter, Z/L (where Z is the height and L is the Obukhov length), this resulted in anomalously low drag coefficients compared to neutral conditions. Determining Z/L iteratively, a low rate of convergence was achieved. It was concluded that the divergence of the turbulent transport of TKE was not negligible under unstable conditions. By minimizing the dependence of the calculated neutral drag coefficient on stability, this term was estimated at about -0.65Z/L. The resulting turbulent fluxes were then in close agreement with other studies at moderate wind speed. The drag and exchange coefficients for low wind speeds were found to be C-en x 10(3) = 2.79U(10n)(-1) + 0.66 (U-10n < 5.2 m/s), C-en x 10(3) = C-hn x 10(3) = 1.2 (U-10n greater than or equal to 5.2 m/s), and C-dn x 10(3) = 11.7U(10n)(-2) + 0.668 (U-10n < 5.5 m/s), which imply a rapid increase of the coefficient values as the wind decreased within the smooth flow regime. The frozen turbulence hypothesis and the assumptions of isotropy and an inertial subrange were found to remain valid at these low wind speeds for these shipboard measurements. Incorporation of a free convection parameterization had little effect.

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

The aim of this study is to clarify the role of the Southern Ocean storms on interior mixing and meridional overturning circulation. A periodic and idealized numerical model has been designed to represent the key physical processes of a zonal portion of the Southern Ocean located between 70 and 40° S. It incorporates physical ingredients deemed essential for Southern Ocean functioning: rough topography, seasonally varying air–sea fluxes, and high-latitude storms with analytical form. The forcing strategy ensures that the time mean wind stress is the same between the different simulations, so the effect of the storms on the mean wind stress and resulting impacts on the Southern Ocean dynamics are not considered in this study. Level and distribution of mixing attributable to high-frequency winds are quantified and compared to those generated by eddy–topography interactions and dissipation of the balanced flow. Results suggest that (1) the synoptic atmospheric variability alone can generate the levels of mid-depth dissipation frequently observed in the Southern Ocean (10−10–10−9 W kg−1) and (2) the storms strengthen the overturning, primarily through enhanced mixing in the upper 300 m, whereas deeper mixing has a minor effect. The sensitivity of the results to horizontal resolution (20, 5, 2 and 1 km), vertical resolution and numerical choices is evaluated. Challenging issues concerning how numerical models are able to represent interior mixing forced by high-frequency winds are exposed and discussed, particularly in the context of the overturning circulation. Overall, submesoscale-permitting ocean modeling exhibits important delicacies owing to a lack of convergence of key components of its energetics even when reaching Δx =  1 km.