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 interrelationships 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 to the strength of the subpolar AMOC and gyre circulation, and it is also linked to the subtropical AMOC, although the strength of this relationship is model-dependent and affected by the inclusion of the Ekman component. 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 to the subtropical AMOC might be more realistic.