Global observations of tropospheric nitrogen dioxide (NO 2 ) columns have been shown to be feasible from space, but consistent multi-sensor records do not yet exist, nor are they covered by planned activities on the international level. Harmonised, multi-decadal records of NO 2 columns and their associated uncertainties can provide crucial information how the emissions and concentrations of nitrogen oxides evolve over time. Here we describe the development of a new, community best practice NO 2 retrieval algorithm based on a synthesis of existing approaches. Detailed comparisons of these approaches led us to implement an enhanced spectral fitting method for NO 2 , a 1° × 1° TM5-MP data assimilation scheme to estimate the stratospheric background, and improve air mass factor calculations. Guided by the needs expressed by data users, producers, and WMO GCOS guidelines, we incorporated detailed per-pixel uncertainty information in the data product, along with easily traceable information on the relevant quality aspects of the retrieval. We applied the improved QA4ECV NO 2 algorithm on the most actual level-1 data sets to produce a complete 22-year data record that includes GOME (1995-2003), SCIAMACHY (2002–2012), GOME-2(A) (2007 onwards) and OMI (2004 onwards). The QA4ECV NO 2 spectral fitting recommendations and TM5-MP stratospheric column and air mass factor approach are currently also applied to S5P-TROPOMI. The uncertainties in the QA4ECV tropospheric NO 2 columns amount to typically 40 % over polluted scenes. First validation results of the QA4ECV OMI NO 2 columns and their uncertainties over Tai’an, China in June 2006 suggests little bias (−27thinsp;%) and better precision than suggested by uncertainty propagation. We conclude that our improved QA4ECV NO 2 long-term data record is providing valuable information to quantitatively constrain emissions, deposition, and trends in nitrogen oxides on a global scale.