Ce travail de thèse a pour objectif de caractériser et quantifier les différences entre indices climatiques océaniques grande échelle simulés par une hiérarchie de neuf modèles numériques océaniques réalistes (cinq simulations libres et quatre réanalyses). Cette étude se focalise sur la région Atlantique Nord. Les indices climatiques suivants : transport méridien de volume (circulation méridienne d'overturning) et transport méridien de chaleur sont diagnostiqués dans les différents modèles. La période commune à huit des simulations s'étend sur neuf ans, l'étude de la variabilité de ces indices climatiques se focalise sur l'echelle interannuelle. L'influence de différents paramètres numériques : configuration et la résolution spatiales des modèles numérique libres, et l'influence de différents schémas d'assimilation de données séquentielle des réanalyses, sont évaluées au regard des observations et estimations disponibles de ces indices climatiques grande échelle. Il ressort de ce travail une importante diversité des solutions des modèles quant aux valeurs moyennes et aux variations de ces transports méridiens de volume et de chaleur. Les réanalyses océaniques étudiées ici, qui ont pourtant été contraintes vers un océan observé, peinent à converger vers des indices climatiques grande échelle cohérents, tant en moyenne qu'en variabilité interannuelle. Le manque de robustesse de ces réanalyses s'exprime par leur faible cohérence entre elles mais aussi au regard des estimations observationnelles disponibles, et illustre ainsi la difficulté de simuler des indices climatiques cohérents à grande échelle en assimilant avec des contraintes locales.