We build up a quantitative second-order Sobolev estimate of lnw for positive p-harmonic functions w in Riemannian manifolds under Ricci curvature bounded from below and also for positive weighted p-harmonic functions w in weighted manifolds under the Bakry-Émery curvature-dimension condition.