In this paper, we obtain gradient continuity estimates for viscosity solutions of ΔNpu=f in terms of the scaling critical L(n,1)-norm of f, where ΔNp is the normalized p-Laplacian operator. Our main result corresponds to the borderline gradient continuity estimate in terms of the modified Riesz potential ̃𝕀fq. Moreover, for f∈Lm with m>n, we also obtain C1,α estimates. This improves one of the regularity results in [A. Attouchi, M. Parviainen and E. Ruosteenoja, C1,α regularity for the normalized p-Poisson problem, J. Math. Pures Appl. (9) 108(4) (2017) 553–591], where a C1,α estimate was established depending on the Lm-norm of f under the additional restriction that p>2 and m>max(2,n,p2). We also mention that differently from the approach in the above paper, which uses methods from divergence form theory and nonlinear potential theory, the method in this paper is more non-variational in nature, and it is based on separation of phases inspired by the ideas in [L. Wang, Compactness methods for certain degenerate elliptic equations, J. Differential Equations 107(2) (1994) 341–350]. Moreover, for f continuous, our approach also gives a somewhat different proof of the C1,α regularity result.