[step:Convert volume equality into equality of the averaged Jacobian quotient]Let $d\sigma_p$ denote the Riemannian measure on the unit sphere $S_pM\subset T_pM$ induced by the [inner product](/page/Inner%20Product) $g_p$; under any linear isometry $T_pM\cong \mathbb{R}^n$, its total mass is $\mathcal{H}^{n-1}(S^{n-1})$. Define the model area density
\begin{align*}
w_k:(0,\rho_k)&\to (0,\infty),\\
t&\mapsto \mathcal{H}^{n-1}(S^{n-1})\operatorname{sn}_k(t)^{n-1},
\end{align*}
and define the averaged quotient
\begin{align*}
A:(0,\infty)&\to [0,\infty),\\
t&\mapsto \frac{1}{\mathcal{H}^{n-1}(S^{n-1})}\int_{S_pM} q(t,\xi)\,d\sigma_p(\xi).
\end{align*}
The cut-time function is measurable and $J$ is smooth on the regular polar-coordinate domain, so the extended quotient $q$ is a nonnegative measurable function on $(0,\infty)\times S_pM$. Hence [Tonelli's Theorem](/page/Tonelli%27s%20Theorem) applies to all iterated integrals below. Since $q(\cdot,\xi)$ is nonincreasing for every $\xi$, integrating the inequality $q(t_1,\xi)\geq q(t_2,\xi)$ for $0<t_1<t_2<R$ over $(S_pM,d\sigma_p)$ gives that $A$ is nonincreasing. The [geodesic polar-coordinate volume formula](/page/Geodesic%20Polar%20Coordinates) gives, for every $\rho\in(0,\rho_k)$,
\begin{align*}
\operatorname{Vol}_g(B(p,\rho))
&=
\int_0^\rho\int_{S_pM} q(t,\xi)\operatorname{sn}_k(t)^{n-1}\,d\sigma_p(\xi)\,d\mathcal{L}^1(t)\\
&=
\int_0^\rho A(t)w_k(t)\,d\mathcal{L}^1(t).
\end{align*}
Also,
\begin{align*}
V_k(\rho)=\int_0^\rho w_k(t)\,d\mathcal{L}^1(t).
\end{align*}
Hence
\begin{align*}
\frac{\operatorname{Vol}_g(B(p,\rho))}{V_k(\rho)}
=
\frac{\int_0^\rho A(t)w_k(t)\,d\mathcal{L}^1(t)}
{\int_0^\rho w_k(t)\,d\mathcal{L}^1(t)}.
\end{align*}
Because $R<\rho_k$, the weight $w_k(t)$ is positive for every $t\in(0,R)$. We use the following monotone weighted-average fact: if $A$ is nonincreasing on $(0,R)$ and $w_k>0$, then
\begin{align*}
\rho\mapsto \frac{\int_0^\rho A(t)w_k(t)\,d\mathcal{L}^1(t)}{\int_0^\rho w_k(t)\,d\mathcal{L}^1(t)}
\end{align*}
is nonincreasing. Indeed, for $0<\rho_1<\rho_2<R$, every value of $A$ on $(0,\rho_1)$ is at least every value of $A$ on $(\rho_1,\rho_2)$ up to null sets, so adding the outer shell can only decrease the weighted average.
Let
\begin{align*}
C=\frac{\operatorname{Vol}_g(B(p,r))}{V_k(r)}=\frac{\operatorname{Vol}_g(B(p,R))}{V_k(R)}.
\end{align*}
The equality of the averages over $(0,r)$ and $(0,R)$ implies that the weighted average of $A$ over the shell $(r,R)$ is also $C$; this follows by subtracting the identity for the weighted integral over $(0,r)$ from the identity over $(0,R)$. Because $A$ is nonincreasing, for every $s\in(0,r)$ and $t\in(r,R)$ one has $A(s)\geq A(t)$. Suppose there were $s_0\in(0,r)$ with $A(s_0)>C$. Then $A(s)\geq A(s_0)>C$ for every $s\in(0,s_0)$, and since $w_k>0$ on $(0,R)$ the weighted average of $A$ over $(0,r)$ would be strictly larger than $C$. Suppose instead there were $t_0\in(r,R)$ with $A(t_0)<C$. Then $A(t)\leq A(t_0)<C$ for every $t\in(t_0,R)$, and since $w_k>0$ on $(0,R)$ the weighted average of $A$ over $(r,R)$ would be strictly smaller than $C$. Thus $A(s)\leq C\leq A(t)$ for all $s\in(0,r)$ and $t\in(r,R)$, while monotonicity gives the reverse inequality across the cut point. It follows that $A(t)=C$ for every $t\in(0,R)$.[/step]