We prove that the Cheng–Yau gradient estimate on positive harmonic functions on
manifolds with nonnegative Ricci curvature is globally stable under certain
perturbations of the metric. In some cases, one only needs the condition
Ricci,
with
and
sufficiently small.