The hydrography and stratification of the subpolar North Atlantic is highly variable, with convection activating and deactivating across parts of the Labrador and Irminger seas. Likely consequential for the Atlantic Meridional Overturning Circulation (AMOC), this variability is examined in an eddy-rich ocean model hindcast spanning 1958–2021 and in 1950–2050 simulations with four climate models, spanning differences in ocean resolution (eddy-rich or eddy-permitting), code and implementation. Stratification of the Labrador and Irminger seas is quantified with the Potential Energy Anomaly (PEA) in the upper 1000 m of the water column. Monthly PEA anomalies are evaluated alongside corresponding anomalies of sea surface temperature, salinity, and density. For 30-year windows, moving correlations between PEA and surface properties are obtained over the 100-year simulations to characterize the evolving relationships. As climate change progresses, stratification in three of the four models is increasingly associated with variable surface salinity, in both regions. Lagrangian analyses of surface flow pathways in the decades preceding 1990 and 2040 are undertaken for one of the models in which surface salinity grows in influence. The subpolar presence of low-salinity Arctic waters and high-salinity subtropical Atlantic waters are found to increase and decrease respectively by 2040. Furthermore, in three of the four models, surface density compensation associated with correlation of surface temperature and salinity anomalies is progressively replaced by combined surface warming and freshening, lowering surface density, and strengthening stratification. The extent of these model-dependent changes and processes are of consequence for the projected fate of the AMOC by the mid twenty-first century.