During the Last Glacial Maximum, the northwestern North Atlantic constituted a major conduit for Labrador and Greenland ice sheet meltwaters. Vertical density gradients in its upper water masses have been reconstructed by combining information from transfer functions based on dinocysts and from oxygen isotope measurements (δ18O) in planktonic foraminifera. Transfer functions yield temperature and salinity and thus potential density (σθ) for the warmest (August) and coldest (February) months in the photic zone. The δ18O values in different size fractions of epipelagic (Globigerina bulloides) and mesopelagic (Neogloboquadrina pachyderma left-coiled (Npl)) foraminifera allow us to assess σθ gradients through the pycnocline between surface and intermediate waters, based on the calibration of a σθ versus δ18O relationship from transfer function reconstructions. The size and density of Npl shells provide further constraints on these σθ gradients. The results show the development of a very strong pycnocline during the LGM with a difference of about 3 (summer) to 1.5 (winter) σθ units between surface and underlying waters. They indicate conditions unfavorable for vertical convection and support the hypothesis of the spreading of a shallow, low-salinity buoyant layer over the northern North Atlantic. This layer depicted a strong E-W gradient, with maximum seasonal contrast and minimum absolute σθ values westward.