Abstract. The Subpolar North Atlantic (SPNA) is a region with prominent decadal variability that has experienced remarkable warming and cooling trends in the last few decades. These observed trends have been preceded by slow-paced increases and decreases in the Labrador Sea density (LSD), which are thought to be a precursor of large scale ocean circulation changes. This article analyses the inter-relationships between the LSD and the wider North Atlantic across an ensemble of coupled climate model simulations. In particular, it analyses the link between subsurface density and the deep boundary density, the Atlantic Meridional Overturning Circulation (AMOC), the Subpolar Gyre (SPG) circulation, and the upper ocean temperature in the eastern SPNA. All simulations exhibit considerable multidecadal variability in the LSD and the ocean circulation indices, which are found to be interrelated. LSD is strongly linked with the strength of subpolar AMOC and gyre circulation, and is also linked with the subtropical AMOC, although the strength of this relationship is model dependent. The connectivity of LSD with the subtropics is found to be sensitive to different model features, including: the mean density stratification in the Labrador Sea; the strength and depth of the AMOC; and the depth at which the LSD propagates southward along the western boundary. Several of these quantities can also be computed from observations, and comparison with these observation-based quantities suggests that models representing a weaker link with the subtropical AMOC may be more realistic. This would imply that RAPID AMOC measurements might not be adequate to represent decadal to multidecadal changes in the subpolar overturning circulation.