We update our approximate parametrizations of the three-loop splitting functions for the evolution of unpolarized parton densities in perturbative QCD. The new information taken into account is given by the additional Mellin moments recently calculated by Retey and Vermaseren. The inclusion of these constraints reduces the uncertainties of our approximations considerably and extends their region of applicability by about one order of magnitude to lower momentum fractions x.PACS: 12.38. Bx, 13.60.Hb In order to achieve a high accuracy of the predictions of perturbative QCD for hard processes, the calculations need to transcend the standard next-to-leading order (NLO) approximation. For processes with initial-state hadrons, the next-to-next-to-leading order (NNLO) expressions include the three-loop splitting functions. The computation of these functions is under way [1], but will not be completed in the near future [2].Partial results have already been obtained [3][4][5][6][7][8][9], however, most notably the five lowest even-integer moments for the flavour non-singlet combination entering electromagnetic deep-inelastic scattering (DIS) [3], and four moments for the singlet splitting functions [4]. In refs. [10,11] we have demonstrated that this information -due to the smoothening effect of the ubiquitous convolution with the initial parton densities and the small size of the corrections -is fully sufficient for momentum fractions x > ∼ 0.1 and leaves only small uncertainties down to x ≃ 10 −3 at scales above about 10 GeV 2 . We have provided approximate parametrizations for the three-loop splitting functions, including quantitative estimates of their residual uncertainties. These results have already been applied [12] to structure functions in DIS and Drell-Yan cross sections at hadron colliders, for which the subprocess cross sections have been computed up to NNLO [13,14].Very recently the fixed-moment calculations of refs. [3,4] have been extended, using improved computing resources, up to the twelfth moment [2]. For the first time also (odd) moments of the three-loop valence splitting functions have been obtained there. These results provide a severe check of our approximation procedure, as the latter led to rather tight predictions for the (tenth and) twelfth moments. This test is passed by the results of refs. [10,11]. In this letter we update these parametrizations by including the moments of ref.[2] in the derivation. As a result the residual uncertainties are greatly reduced, and the region of safe applicability is extended by about one order of magnitude in x, an improvement most relevant for applications to structure functions at HERA [12].Our notations for the parton densities and splitting functions are as follows: the nonsinglet combinations of quark and antiquark densities, q i andq i , are given byN f stands for the number of effectively massless flavours. The corresponding splitting functions are denoted by P ± NS and P V NS ≡ P − NS + P S NS . The latter function, P S NS , occurs for t...