We developed a novel method to combine ocean data from observations and models, and applied it to produce gridded estimates of nutrients, oxygen, dissolved inorganic carbon and total alkalinity concentrations at latitudes >40° N and years 1980–2020. The new estimates showed improved accuracy and coverage relative to previous estimates, but highlighted remaining uncertainty in some poorly sampled regions. The work was largely motivated by a need for accurate input data for regional ocean models.