In December 2019, the 13th revision of the International Geomagnetic Reference Field (IGRF) was released by the International Association of Geomagnetism and Aeronomy (IAGA) Division V Working Group V-MOD. This revision comprises two new spherical harmonic main eld models for epochs 2015.0 (DGRF-2015) and 2020.0 (IGRF-2020) and a model of the predicted secular variation for the interval 2020.0 to 2025.0 (SV-2020-2025). The models were produced from candidates submitted by fteen international teams. These teams were led by the British Geological Survey (UK), China Earthquake Administration (China), Universidad Complutense de Madrid (Spain), University of Colorado Boulder (USA), Technical University of Denmark (Denmark), GFZ German Research Centre for Geosciences (Germany), Institut de physique du globe de Paris (France), Institut des Sciences de la Terre (France), Pushkov Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation (Russia), Kyoto University (Japan), University of Leeds (UK), Max Planck Institute for Solar System Research (Germany), NASA Goddard Space Flight Center (USA), University of Potsdam (Germany), and Universite de Strasbourg (France). The candidate models were evaluated individually and compared to all other candidates as well to the mean, median and a robust Huber-weighted model of all candidates. These analyses were used to identify, for example, the variation between the Gauss coecients or the geographical regions where the candidate models strongly di ered. The majority of candidates were suciently close that the di erences can be explained primarily by individual modeling methodologies and data selection strategies. None of the candidates were so di erent to warrant their exclusion from the nal IGRF-13. The IAGA V-MOD task force thus voted for two approaches: the median of the Gauss coecients of the candidates for the DGRF-2015 and IGRF-2020 models and the robust Huber-weighted model for the predictive SV-2020-2025. In this paper, we document the evaluation of the candidate models and provide details of the approach used to derive the nal IGRF-13 products.