Measurement of ocean temperatures using instruments carried by Antarctic fur seals

The study aimed to test the utility of instruments deployed on marine mammals for measuring physical oceanographic variation and, using this method, to examine temperature variation in the coastal waters around South Georgia. There was a significant correlation between temperature measurements made using a towed undulating oceanographic recorder (UOR) and concurrent measurements from time-depth recorders (TDRs) fitted to lactating Antarctic fur seals foraging from the coast of South Georgia. Congruence was found at horizontal spatial scales from 0.01°×0.01° to 0.5°×0.5° (degrees of latitude and longitude), and at a vertical scale of 10 m. However, there was no significant correlation between temperature measured by TDRs in the top 5 m and sea surface temperature (SST) measured by satellite remote sensing. TDR data provided information about temperature variation vertically through the water column, and through time. The UOR data were used to recalibrate the TDR data in order to correct for the slow response time of the TDR thermistor relative to the speed of seal movements through the water column. Seasonal temperature variation was apparent, and temperatures also varied between regions, and with bathymetry. These results were consistent with the current interpretation of the coastal oceanography around South Georgia. In particular, the relationship between on- and off-shelf waters showed larger amounts of warmer surface water in a region in which more run-off was to be expected. The study also showed that Antarctic fur seals concentrate their activity in regions of colder, and presumably oceanic, water. Such instrumented animals could provide near real time data for assimilation into ocean models.     

Analysis of the contribution of ice algae to the ice-covered ecosystem in Lake Saroma by means of a coupled ice–ocean ecosystem model

A new coupled ice–ocean ecosystem model that links the pelagic and ice ecosystems was used to clarify the role of ice algae in ice-covered ocean ecosystems. The model was applied to Lake Saroma (Hokkaido, Japan) in 1992. Comparison of the model's results with observational data confirmed that the model reproduced the behavior of the ecosystem with acceptable accuracy during the period from winter to spring. The primary production of the ice algae is effectively transported into the pelagic system by means of physical releasing effects: brine convection, ice melting and freezing, and diffusion generated at the bottom of the ice. Ice algae released from the ice are rapidly exported because of their high sinking speed and the shallow depth of Lake Saroma. For this reason, the zooplankton in Lake Saroma cannot graze these released algae. However, zooplankton actively graze the ice algae living along the bottom of the ice. These results show that, before their release, ice algae play an important role as a food source for overwintering zooplankton. A sensitivity analysis revealed a positive correlation between the sinking speed of the released ice algae and the magnitude of the spring bloom by pelagic phytoplankton, and that the time when secondary production becomes active is an important factor in the linkage between these two algal populations.     

Mesoscale subduction at the Almeria–Oran front: Part 1: Ageostrophic flow

This paper presents a detailed diagnostic analysis of hydrographic and current meter data from three, rapidly repeated, fine-scale surveys of the Almeria–Oran front. Instability of the frontal boundary, between surface waters of Atlantic and Mediterranean origin, is shown to provide a mechanism for significant heat transfer from the surface layers to the deep ocean in winter. The data were collected during the second observational phase of the EU funded OMEGA project on RRS Discovery cruise 224 during December 1996. High resolution hydrographic measurements using the towed undulating CTD vehicle, SeaSoar, traced the subduction of Mediterranean Surface Water across the Almeria–Oran front. This subduction is shown to result from a significant baroclinic component to the instability of the frontal jet. The Q-vector formulation of the omega equation is combined with a scale analysis to quantitatively diagnose vertical transport resulting from mesoscale ageostrophic circulation. The analyses are presented and discussed in the presence of satellite and airborne remotely sensed data; which provide the basis for a thorough and novel approach to the determination of observational error.     

海洋船舶的轨迹估计算法

准确估计和预测海洋导航中机动船舶的轨迹是改善海上安全和保安的重要工具。因此,许多常规海洋导航系统和船舶交通管理和报告服务都为此目的配备了雷达设施。然而,船舶操纵轨迹的预测的准确性主要取决于船舶位置,速度和加速度的估计的良好性。因此,本研究提出了一个基于曲线运动模型的机动海洋船只模型,其基于用于相同目的的线性位置模型的测量。此外,还假定与白高斯噪声相关联的系统状态和测量模型。扩展卡尔曼滤波器被提出作为用于估计位置,速度和加速度的自适应滤波器算法,用于预测机动的海洋船舶轨迹。最后,提出的模型被实施,并获得成功的计算结果,在本研究中的海洋导航中的船舶操纵轨迹的预测。     

基于小波的SPIHT雷达图像压缩算法

雷达在现代航海导航中占据着很重要的地位,在综合船桥系统和船舶交通管理系统里,雷达图像需要通过网络传送给其他导航系统如ECDIS或远程图像监视系统。由于雷达图像数据量非常大,为了保证图像的实时性,通过网络传输雷达图像时,先要对雷达图像进行压缩。雷达天线的扫描速度大约为20 RPM,所以一幅雷达图像的数据压缩、数据传输和数据重建必须在3 s内完成,这对雷达图像的压缩和重建提出了很高的要求。论文采用基于小波变换的SPIHT算法对雷达图像进行压缩.在追求高压缩比率情况下,保证了较低的均方误差和较高的信噪比,最后用C 实现了算法并给出了实验对比数据。     

Quality control of ocean temperature and salinity profiles — Historical and real-time data

A new automated quality control system for oceanic temperature and temperature–salinity profiles is presented. Substantial development was needed for some of the quality control algorithms although the checks were based on documented procedures used elsewhere where appropriate. A new automated ship track check was developed: the results of an undetected position error can be very damaging to ocean analyses. Also important is a check against a gridded background, this can be a climatology but near the surface it is advantageous to use an estimate that is evolving over time. Bayesian probability theory is used in the background check and the associated check against nearby observations (buddy check). The system was used to process archive data for 1956–2004. As a by-product monthly model-free objective analyses for this period were produced. Versions of the system are used for near-real time ocean analysis and for initialising both short-range ocean forecasts and seasonal atmosphere–ocean forecasts. The main features of the oceanic observing systems are presented along with quality control statistics, examples of errors that can occur and some additional problematic cases.     

High-resolution synthetic monitoring by a 4-dimensional variational data assimilation system in the northwestern North Pacific

A data assimilation system is applied to the integrated monitoring of oceanic states in the northwestern North Pacific by combining a high resolution ocean general circulation model with an adjoint method. A comparison of assimilation results with observations shows that the system is better able to represent synoptic features of ocean circulation than do models or data alone. Furthermore, meso-scale features associated with frontal structures and eddies, which are often seen in the Kuroshio and Oyashio extension regions and the Sea of Japan, are better defined in the assimilation results. These features suggest that our 4D-VAR high-resolution data assimilation system is capable of providing time series data which satisfy the model physics and fit the observations, and hence the ocean state derived from our system has greater information content than that obtained from earlier methods.     

Preoperational ocean forecasting in the southeastern Mediterranean Sea: Implementation and evaluation of the models and selection of the atmospheric forcing

Within the framework of several local and international programs, a quasi-operational ocean-forecasting system for the Southeastern Mediterranean Sea has been established and evaluated through a series of preoperational tests. The Princeton Ocean Model (POM) is used for simulating and predicting the hydrodynamics while the Wave Model (WAM) is used for predicting surface waves. Both models were set up to allow varying resolution and multiple nesting. In addition, POM was set up to be easily relocatable to allow rapid deployment of the model for any region of interest within the Mediterranean Sea. A common requirement for both models is the need for atmospheric forcing. Both models require time varying wind or wind stress. In addition, the hydrodynamic model requires initial conditions as well as time dependent surface heat fluxes, fresh water flux, and lateral boundary conditions at the open boundaries. Several sources of atmospheric forcing have been assessed based on their availability and their impact on the quality of the ocean models' forecasts. The various sources include operational forecast centers, other research centers, as well as running an in-house regional atmospheric model. For surface waves, higher spatial and temporal resolution of the winds plays a central role in improving the forecasts in terms of significant wave height and the timing of various high wave events. For the hydrodynamics, using the predicted wind stress and heat fluxes directly from an atmospheric model can potentially produce short range ocean forecasts that are nearly as good as hindcasts forced with gridded atmospheric analyses. Finally, a high-resolution, nested version of the model has shown to be stable under a variety of forcing conditions and time scales, thus indicating the robustness of the selected nesting strategy. For the southeastern corner of the Mediterranean, at forecast lead times of up to 4 days the high-resolution model shows improved skill over the coarser resolution driving model when compared to satellite derived sea surface temperatures. Most of the error appears to be due to the analysis error inherent in the initial conditions.     

VHF海上安全信息自动播发系统在VTS应用现状与工作建议

海上安全信息是船舶航行及作业安全的重要保障。VTS管辖水域通常船舶交通流密集、生产繁忙,为更好地满足船舶信息需求、维护辖区水域安全,VTS中心开始在传统人工服务手段的基础上探索实现海上安全信息的自动播发。然而,目前VTS人员对于VHF自动播发语音信息的性质存在认知差异,自动播发系统的使用存在较大的随意性。该文对VHF自动播发语音信息的性质、VHF海上安全信息自动播发系统应用需求及应用现状进行了分析,在此基础上提出了相关工作建议。     

Variability of the western Galician upwelling system (NW Spain) during an intensively sampled annual cycle. An EOF analysis approach

The key features of the western Galician shelf hydrography and dynamics are analyzed on a solid statistical and experimental basis. The results allowed us to gather together information dispersed in previous oceanographic works of the region. Empirical orthogonal functions analysis and a canonical correlation analysis were applied to a high-resolution dataset collected from 47 surveys done on a weekly frequency from May 2001 to May 2002. The main results of these analyses are summarized bellow. Salinity, temperature and the meridional component of the residual current are correlated with the relevant local forcings (the meridional coastal wind component and the continental run-off) and with a remote forcing (the meridional temperature gradient at latitude 37°N). About 80% of the salinity and temperature total variability over the shelf, and 37% of the residual meridional current total variability are explained by two EOFs for each variable. Up to 22% of the temperature total variability and 14% of the residual meridional current total variability is devoted to the set up of cross-shore gradients of the thermohaline properties caused by the wind-induced Ekman transport. Up to 11% and 10%, respectively, is related to the variability of the meridional temperature gradient at the Western Iberian Winter Front. About 30% of the temperature total variability can be explained by the development and erosion of the seasonal thermocline and by the seasonal variability of the thermohaline properties of the central waters. This thermocline presented unexpected low salinity values due to the trapping during spring and summer of the high continental inputs from the River Miño recorded in 2001. The low salinity plumes can be traced on the Galician shelf during almost all the annual cycle; they tend to be extended throughout the entire water column under downwelling conditions and concentrate in the surface layer when upwelling favourable winds blow. Our evidences point to the meridional temperature gradient acting as an important controlling factor of the central waters thermohaline properties and in the development and decay of the Iberian Poleward Current.     

Parameterisation of clastic sediments including benthic structures

The sediment transport processes in the south-western Baltic Sea are predicted by means of a numerical model in the project DYNAS. There are two sediment parameters that influence the results of modelling remarkably: critical shear stress velocity and bottom roughness. This paper presents the way how to parameterise these factors and extrapolate them into the investigation area. The critical shear stress velocity is parameterised basing on grain size data, combining approximations after Hjulström [Hjulström, F., 1935: Studies in the morphological activity of rivers as illustrated by the river Fyris. Geological Institution of University of Uppsala: Bulletin (25): 221–528.], Shields [Shields, A., 1936: Anwendung der Ähnlichkeits-Mechanik und der Turbulenzforschung auf die Geschiebebewegung. Mitteilungen der Preussischen Versuchsanstalt für Wasserbau und Schiffahrt (26): 26 pp.] and Bohling [Bohling, B., 2003: Untersuchungen zur Mobilität natürlicher und anthropogener Sedimente in der Mecklenburger Bucht. unpublished doctoral thesis, Mathematisch-Naturwissenschaftliche Fakultät, Ernst-Moritz-Arndt-Universität Greifswald/Germany, 156 pp.]. The roughness length, in the case of absence of macro zoo-benthos and their structures, is parameterised basing on grain size too employing Soulsby [Soulsby, R.L., 1997: Dynamics of Marine Sands: a Manual for Practical Applications. London, Thomas Telford Publications. 249 pp.], Nielsen [Nielsen, P., 1983: Analytical determination of nearshore wave height variation due to refraction shoaling and friction. Coastal Engineering 7, 233–251.] and Yalin [Yalin, M.S., 1977: Mechanics of Sediment Transport. Pergamon Press, New York. 298 pp.]. No equivalent simple parameterisations for biologically caused bed roughness exist. Here, findings of Friedrichs [Friedrichs, M., 2004: Flow-induced effects of macro zoo-benthic structures on the near-bed sediment transport. Dissertation, Universität Rostock, 80 S.] and estimations by the DYNAS biologists group were combined in order to derive roughness lengths from abundance measurements of four previously selected key species which represent the originators of the dominating benthic structures at the sea floor in the south-western Baltic Sea. Critical shear stress velocity and bed roughness are known at few sample sites only. They were extrapolated into the larger investigation area using a proxy-target concept. The mean near bottom milieu (bathymetry, median grain size, salinity, oxygen) which was derived using results from numerical modelling serves as the proxy. Since the milieu parameters are measured at the sampling sites for which the target parameters have been determined, a combined hierarchical and supervised classification was employed to transfer the local knowledge into the unknown investigation area.     

杭州湾跨海大桥VTS软科学研究

船舶交通管理系统(VTS)是一项牵涉面非常广、极复杂的系统工程，开展软科学研究，是合理规划和科学管理的前提和依据。杭州湾跨海大桥桥区水域建立船舶交通管理系统，条件、情况比较特殊复杂，应该更加重视和加强实现系统功能所必须依赖的软科学研究。     

Validation of a hybrid optimal interpolation and Kalman filter scheme for sea surface temperature assimilation

A hybrid data assimilation scheme designed for operational assimilation of satellite sea surface temperatures (SST) into an ocean model has been developed and validated against in-situ observations. The scheme consists of an optimal interpolation (OI) part and a greatly simplified Kalman filter (KF) part.The OI is performed only in the longitudinal and latitudinal directions. A climatological field is used as a background field for the interpolation. It is constructed by fitting daily averages of satellite SST to the annual mean, annual, and semiannual harmonics in a 20 km by 20 km grid. The background error covariance is approximated by a spatially varying two-dimensional exponential covariance model. The parameters of the covariance model are fitted to the deviations of the satellite data from the background field using data from a full year.The simplified KF uses ocean model forecasts as a background field. It is based on the assumption that it is possible to neglect horizontal SST covariances in the filter and that the typical time scale for vertical mixing in the mixed layer is much shorter than the average time between observations. We therefore assume that the error variance in a column of water is evenly spread out throughout the mixed layer. The result of these simplifications is a computationally very efficient KF.A one year validation of the scheme is performed for year 2001 using an operational eddy resolving ocean model covering the North Sea and the Baltic Sea. It is found that assimilation of sea surface temperature data reduces the model root mean square error from 1.13 °C to 0.70 °C. The hybrid scheme is found to reduce the root mean square error slightly more than the simplified KF without OI to 0.66 °C. The inclusion of spatially varying satellite error variances does not improve the performance of the scheme significantly.     

Wave data assimilation with the Kalman filter

A new data assimilation method for ocean waves is presented, based on an efficient low-rank approximation to the Kalman filter. Both the extended Kalman filter and a truncated second-order filter are implemented. In order to explicitly estimate past wind corrections based on current wave measurements, the filter is extended to a fixed-lag Kalman smoother for the wind fields. The filter is tested in a number of synthetic experiments with simple geometries. Propagation experiments with errors in the boundary condition showed that the KF was able to accurately propagate forecast errors, resulting in spatially varying error correlations, which would be impossible to model with time-independent assimilation methods like OI. An explicit comparison with an OI assimilation scheme showed that the KF also is superior in estimating the sea state at some distance from the observations. In experiments with errors in the driving wind, the modeled error estimates were also in agreement with the actual forecast errors. The bias in the state estimate, which is introduced through the nonlinear dependence of the waves on the driving wind field, was largely removed by the second-order filter, even without actually assimilating data. Assimilation of wave observations resulted in an improved wave analysis and in correction of past wind fields. The accuracy of this wind correction depends strongly on the actual place and time of wave generation, which is correctly modeled by the error estimate supplied by the Kalman filter. In summary, the KF approach is shown to be a reliable assimilation scheme in these simple experiments, and has the advantage over other assimilation methods that it supplies explicit dynamical error estimates.     

Surface chlorophyll in the Black Sea over 1978–1986 derived from satellite and in situ data

Absolute values of chlorophyll a concentration and its spatial and seasonal variations in the Black Sea were assessed by using satellite CZCS and in situ data. Since the satellite CZCS had operated for the 1978–1986 period, CZCS data was used for assessing the past state of the Black Sea just before the onset of drastic changes observed in late 1980s. The approach used for the calculation of the absolute values of chlorophyll a concentration from CZCS data was based on the direct comparison of in situ chlorophyll a data and those of CZCS and by applying the algorithm developed for the transformation of CZCS data into chlorophyll a values. CZCS Level 2 data related with pigment concentration having a spatial resolution of 1 km at nadir were used. The daily Level 3 files were derived by binning Level 2 values into 4-km grid cells and the monthly and seasonal Level 3 files were created by averaging the daily Level 3 files over the corresponding period. In situ chlorophyll a data were obtained by spectrophotometric and fluorometric methods in 15 scientific cruises over the 1978–1986 period. Total number of ship-measured data used for the comparison with those CZCS values was 590.Chlorophyll a concentration (Chl) was derived from CZCS values (C) with regression equations Chl=kC; the coefficient of transformation k was calculated from six different data sets by taking into account distinctions between subregions and seasons. The reasons for difference in the k values have been analyzed.Statistical comparison of the chlorophyll a values measured in situ and those derived from CZCS data was based on log-transformed data and gave the following results: regression SLOPE=0.842, regression INTERCEPT=−0.081, coefficient of determination (R²)=0.806, root–mean–square ERROR=0.195. The mean monthly chlorophyll a distributions derived from CZCS data over 1978–1986 have been constructed and the mean seasonal chlorophyll a values in different regions have been calculated and analyzed. The significant difference in chlorophyll concentration between the western shelf regions and the open part of the Black Sea has been demonstrated, especially in warm season. At almost all seasons, the highest chlorophyll concentration is observed in the western interior shelf region which is under strong influence of Danube. The summer mean chlorophyll concentration in this region is 18 times higher than that in the open parts and about nine times higher than in the eastern shelf region. The greatest seasonal variations are observed in the open part of the Black Sea: chlorophyll concentration in cold season is four to six times higher than in summer and three to five times higher than in April and October. To the contrary, in the western interior shelf regions, the concentration is higher in May–October (about twice than that in November–March). Seasonal variations in the western outer shelf regions are smoothed out as compared with both the western interior shelf and the open regions.     

1/32° real-time global ocean prediction and value-added over 1/16° resolution

A 1/32° global ocean nowcast/forecast system has been developed by the Naval Research Laboratory at the Stennis Space Center. It started running at the Naval Oceanographic Office in near real-time on 1 Nov. 2003 and has been running daily in real-time since 1 Mar. 2005. It became an operational system on 6 March 2006, replacing the existing 1/16° system which ceased operation on 12 March 2006. Both systems use the NRL Layered Ocean Model (NLOM) with assimilation of sea surface height from satellite altimeters and sea surface temperature from multi-channel satellite infrared radiometers. Real-time and archived results are available online at http://www.ocean.nrlssc.navy.mil/global_nlom. The 1/32° system has improvements over the earlier system that can be grouped into two categories: (1) better resolution and representation of dynamical processes and (2) design modifications. The design modifications are the result of accrued knowledge since the development of the earlier 1/16° system. The improved horizontal resolution of the 1/32° system has significant dynamical benefits which increase the ability of the model to accurately nowcast and skillfully forecast. At the finer resolution, current pathways and their transports become more accurate, the sea surface height (SSH) variability increases and becomes more realistic and even the global ocean circulation experiences some changes (including inter-basin exchange). These improvements make the 1/32° system a better dynamical interpolator of assimilated satellite altimeter track data, using a one-day model forecast as the first guess. The result is quantitatively more accurate nowcasts, as is illustrated by several model-data comparisons. Based on comparisons with ocean color imagery in the northwestern Arabian Sea and the Gulf of Oman, the 1/32° system has even demonstrated the ability to map small eddies, 25–75 km in diameter, with 70% reliability and a median eddy center location error of 22.5 km, a surprising and unanticipated result from assimilation of altimeter track data. For all of the eddies (50% small eddies), the reliability was 80% and the median eddy center location error was 29 km. The 1/32° system also exhibits improved forecast skill in relation to the 1/16° system. This is due to (a) a more accurate initial condition for the forecast and (b) better resolution and representation of critical dynamical processes (such as upper ocean – topographic coupling via mesoscale flow instabilities) which allow the model to more accurately evolve these features in time while running in forecast mode (forecast atmospheric forcing for the first 5 days, then gradually reverting toward climatology for the remainder of the 30-day forecast period). At 1/32° resolution, forecast SSH generally compares better with unassimilated observations and the anomaly correlation of the forecast SSH exceeds that from persistence by a larger amount than found in the 1/16° system.     

Operational surface drift prediction using linear and non-linear hyper-ensemble statistics on atmospheric and ocean models

Multimodel super-ensemble forecasts, which exploit the power of an optimal local combination of individual models usually show superior forecasting skills when compared to individual models because they allow for local correction and/or bias removal. Deterministic approaches to the problem of surface drift are often limited by strong assumptions on the underlying physics. A new approach based on linear and non-linear optimization is proposed, using hyper-ensemble deduced statistics to forecast at short time scale Lagrangian drifts from combined atmospheric and ocean operational models and local observations that were made available during the MREA04 field experiment along the West coast of Portugal. Optimization methods are based on a training/forecast cycle. The performance and the limitations of the hyper-ensembles and the individual models are discussed. Results suggest that our statistical methods reduce the position errors significantly for 12 to 48 h forecasts and hence compete with pure deterministic approaches.     

The Mercator global ocean operational analysis system: Assessment and validation of an 11-year reanalysis

This paper presents Prototype Système 2 Global (PSY2G), the first Mercator global Ocean General Circulation Model (OGCM) to assimilate along-track sea level anomaly (SLA) satellite data. Based on a coarse resolution ocean model, this system was developed mainly for climatic purposes and will provide, for example, initial oceanic states for coupled ocean-atmosphere seasonal predictions. It has been operational since 3 September 2003 and produces an analysis and a two-week forecast for the global ocean every week. The PSY2G system uses an incremental assimilation scheme based on the Cooper and Haines [Cooper, M., Haines, K., 1996. Data assimilation with water property conservation. J. Geophys. Res., 101, 1059-1077.] lifting–lowering of isopycnals. The SLA increment is obtained using an optimal interpolation method then the correction is partitioned into baroclinic and barotropic contributions. The baroclinic ocean state correction consists of temperature, salinity and geostrophic velocity increments and the barotropic correction is a barotropic velocity increment. A reanalysis (1993–2003) was carried out that enabled the PSY2G system to perform its first operational cycle. All available SLA data sets (TOPEX/Poséïdon, ERS2, Geosat-Follow-On, Jason1 and Envisat) were assimilated for the 1993–2003 period. The major objective of this study is to assess the reanalysis from both an assimilation and a thermodynamic point of view in order to evaluate its realism, especially in the tropical band which is a key region for climatic studies. Although the system is also able to deliver forecasts, we have mainly focused on analysis. These results are useful because they give an a priori estimation of the qualities and capabilities of the operational ocean analysis system that has been implemented. In particular, the reanalysis identifies some regional biases in sea level variability such as near the Antarctic Circumpolar Current, in the eastern Equatorial Pacific and in the Norwegian Sea (generally less than 1 cm) with a small seasonal cycle. This is attributed to changes in mean circulation and vertical stratification caused by the assimilation methodology. But the model's low resolution, inaccurate physical parameterisations (especially for ocean–ice interactions) and surface atmospheric forcing also contribute to the occurrence of the SLA biases. A detailed analysis of the thermohaline structure of the ocean reveals that the isopycnal lifting–lowering tends to diffuse vertically the main thermocline. The impact on temperature is that the surface layer (0–200 m) becomes cooler whereas in deeper waters (from 500 to 1500 m), the ocean becomes slightly warmer. This is particularly true in the tropics, between 30°N and 30°S. However it can be demonstrated that the assimilation improves the variability in both surface currents and sub-surface temperature in the Equatorial Pacific Ocean.     

Inter-comparing five forecast operational systems in the North Atlantic and Mediterranean basins: The MERSEA-strand1 methodology

The methodology to achieve a real time inter-comparison of five state-of-the-art operational forecast systems for the North Atlantic and Mediterranean basins is presented. All systems provide analysis and near real-time prediction of the three-dimensional ocean through Opendap servers. A standard set of diagnostics called metrics, is described. Definition and examples of metrics are given. An inter-comparison of the five systems is conducted over a 1 year period using those metrics. It is shown that the methodology developed allows a successful inter-comparison. It has been adopted by the GODAE community. It is also shown that the systems are consistent with the current knowledge of the ocean circulation and climatologies. Systems are deficient in the representation of specific water masses characteristics as Mode waters. Data assimilation of vertical profiles of temperature and salinity solve such deficiencies. Metrics also allow a monitoring of the system's North Atlantic overturning stream function and will allow detecting any changes in the coming year system's thermohaline circulation.