The following argument results in a conclusion that I find strange, and makes me suspect there is something wrong with the reasoning.
First, consider a timelike geodesic $\gamma$ with normalized tangent vector $k^a$. Then consider a set $x^a_{(i)}$ of $(d-1)$ orthonormal spacelike vectors at some point $p$ on $\gamma$ which are also orthogonal to $k^a$ there:
$g_{ab} x^a_{(i)} x^b_{(j)} = \delta_{ij}$ and $g_{ab} x^a_{(i)} k^b = 0$.
The $d$ vectors $(k^a, x^a_{(i)})$ form an orthonormal frame at $p$, and by parallel transporting the $x^a_{(i)}$ along $\gamma$, they yield an orthonormal frame at every point on $\gamma$.
Now, introduce a Fermi coordinate system in a neighborhood of $\gamma$. That is, identify points on $\gamma$ by their proper time $\tau$, and for each $\tau$, introduce the proper distances $y^i$ along spatial geodesics fired off of $\gamma$ in the directions $x^a_{(i)}$. The Fermi coordinates are $(\tau, y^i)$.
$k^a$ and $x^a_{(i)}$ can be thought of as coordinate basis vector fields in a neighborhood of $\gamma$: $k^a = (\partial/\partial \tau)^a$ and $x^a_{(i)} = (\partial/\partial y^i)^a$. They therefore all commute:
$[k,x_{(i)}]^a = 0 = [x_{(i)},x_{(j)}]^a$.
Since the $x^a_{(i)}$ are parallel transported along $\gamma$, we have $k^a \nabla_a x^b_{(i)} = 0$, and therefore the vanishing commutator implies $x^a_{(i)} \nabla_a k^b = 0$. But we also had $k^a \nabla_a k^b = 0$ from the geodesic equation. Since the $(k^a, x^a_{(i)})$ form a basis, these conditions together imply that $\nabla_a k^b = 0$.
I find the result $\nabla_a k^b = 0$ strange and possibly incorrect, but I can't find an error in the reasoning above. Is this argument actually correct, or is there a flaw somewhere?