First, an introduction.
Consider the surface of the Earth (whose mass is $M$): the geometry of the spacetime nearby can be described to a first approximation by the Schwarzschild Metric.
If we wish to compute the 4-acceleration of an observer which is stationary with respect to the $(r, \theta, \varphi)$ coordinates, we must do so with the relativistic expression $$ a^\mu = u^\nu \nabla_\nu u^\mu, $$ where $u^\mu \propto (1, 0, 0, 0)$ is the four-velocity of the observer. Doing the computation, we get that the only nonzero component of the acceleration is
$$ a^r = \frac{GM}{r^2}, $$ where $r$ is the (fixed) value of the radial coordinate, so the magnitude of the acceleration is $$ \left|a\right| = \frac{1}{\sqrt{1 - 2GM / r}} \frac{GM}{r^2}. $$
This may be used to explain why we perceive a force downwards: we are in a frame which is accelerating upwards, as per the equivalence principle.
In the computation the relevant Christoffel symbol is $$ \Gamma^r_{tt} = \left( 1 - \frac{2GM}{r}\right) \frac{GM}{r^2}, $$ which describes how the perturbation to the metric component $g_{tt} = - (1 + 2 \Phi)$ varies, that is, how much faster time passes at different elevations.
Now, to my question: is there an intuitive way to explain how this result for the acceleration comes about?
Mathematically it is clear to me how the fact that time passes at different rates at different elevations implies that stationary observers with respect to the coordinates are accelerating, but I would have a hard time explaining it without introducing the whole differential-geometric framework.
Ideally, I'm looking for a graphical explanation which could be understood without a lot of math, but a different approach to the computation could also be useful.
Edit: Such an argument might be only possible in the weak field limit: it is definitely fair for an answer to only apply in the limit $r \gg GM $.