Non-conservative processes change the potential vorticity (PV) of the upper ocean, and later, through the subduction of surface waters into the interior, affect the general ocean circulation. Here we focus on how boundary layer turbulence, in the presence of submesoscale horizontal buoyancy gradients, generates a source of potential vorticity at the ocean surface through a balance known as the Turbulent ThermalWind. This source of PV injection at the submesoscale can be of similar magnitude to PV fluxes from the wind and surface buoyancy fluxes, and hence can lead to a net injection of PV onto outcropped isopycnals even during periods of surface buoyancy loss. The significance of these dynamics is illustrated using a high-resolution realistic model of the North Atlantic Subtropical Mode Water (18° water), where it is demonstrated that injection of PV at the submesoscale reduces the rate of mode-water PV removal by a factor of ~ 2, and shortens the annual period of mode water formation by ~3 weeks, relative to air-sea fluxes alone. Submesoscale processes thus provide a direct link between small-scale boundary layer turbulence and the gyre-scale circulation, through their effect on mode water formation, with implications for understanding the variability and biogeochemical properties of ocean mode waters globally.