[step:Define the Busemann functions and protected levels]
Let $\mathcal{R}_o$ denote the set of all unit-speed geodesic rays $\rho: [0,\infty) \to M$ with $\rho(0)=o$. For each $\rho \in \mathcal{R}_o$, define the Busemann function
\begin{align*}
b_\rho: M &\to \mathbb{R}, \\
x &\mapsto \lim_{t\to\infty}\bigl(t-d_g(x,\rho(t))\bigr),
\end{align*}
where $d_g: M \times M \to [0,\infty)$ is the Riemannian distance induced by $g$. For fixed $x\in M$, define $f_x:[0,\infty)\to\mathbb{R}$ by $f_x(t)=t-d_g(x,\rho(t))$. If $0\le s\le t$, then the triangle inequality gives
\begin{align*}
d_g(x,\rho(t))\le d_g(x,\rho(s))+d_g(\rho(s),\rho(t))=d_g(x,\rho(s))+t-s,
\end{align*}
so $f_x(s)\le f_x(t)$. Thus $f_x$ is nondecreasing. The triangle inequality also gives
\begin{align*}
t-d_g(o,x)\le d_g(x,\rho(t))\le t+d_g(o,x),
\end{align*}
because $d_g(o,\rho(t))=t$, and hence
\begin{align*}
-d_g(o,x)\le f_x(t)\le d_g(o,x).
\end{align*}
Therefore $f_x$ is monotone and bounded, so the defining limit for $b_\rho(x)$ exists. Define, for each $a\ge 0$,
\begin{align*}
C_a := \bigcap_{\rho\in \mathcal{R}_o}\{x\in M : b_\rho(x)\ge -a\}.
\end{align*}
[/step]