A new global dataset of annually averaged ultrafine particle (UFP) concentrations at the Earth's surface for the years 2015–2017 has been developed through numerical simulations using the ECHAM/MESSy Atmospheric Chemistry model (EMAC). We present total and size-resolved concentrations along with their interannual variability. Size distributions of emitted particles from the contributing source sectors have been derived based on literature reports. The model results of UFP concentrations are evaluated using particle size distribution and particle number concentration measurements from available datasets and the literature. While we obtain reasonable agreement between the model results and observations (logarithmic-scale correlation of r=0.76 for non-remote, polluted regions), the highest values of observed, street-level UFP concentrations are systematically underestimated, whereas in rural environments close to urban areas the model generally overestimates observed UFP concentrations. As the relatively coarse global model does not resolve concentration gradients in urban centres and industrial UFP hotspots, high-resolution data of anthropogenic emissions are used to account for such differences in each model grid box, obtaining UFP concentrations with unprecedented horizontal resolution at the Earth's surface. This observation-guided downscaling further improves the agreement with observations, leading to an increase in the logarithmic-scale correlation between observed and simulated UFP concentrations to r=0.84 in polluted environments (and 0.95 in all regions), a decrease in the root mean squared logarithmic error (from 0.57 to 0.43), and removal of discrepancies associated with air quality and population density gradients within the model grid boxes. The model results are made publicly available for studies on public health and other impacts of atmospheric UFPs, as well as for intercomparison with other regional and global models and datasets.