[guided]The only geometric input needed is the behaviour of Ricci curvature under multiplying the metric by a positive constant. We prove it through the Levi-Civita connection rather than through indexed coordinate components.
Fix $t_0 \in I$, and define a new metric $h$ on $M$ by
\begin{align*}
h=\lambda g(t_0).
\end{align*}
Let $\mathfrak{X}(M)=\Gamma(TM)$ be the space of smooth vector fields on $M$, and choose arbitrary vector fields $X,Y,Z \in \mathfrak{X}(M)$. Write $\nabla^g$ for the Levi-Civita connection of $g(t_0)$ and $\nabla^h$ for the Levi-Civita connection of $h$. The Koszul formula characterizes the Levi-Civita connection from the metric and Lie brackets of vector fields. Applied to $h$, it says
\begin{align*}
2h(\nabla^h_XY,Z)=X(h(Y,Z))+Y(h(Z,X))-Z(h(X,Y))-h(X,[Y,Z])+h(Y,[Z,X])+h(Z,[X,Y]).
\end{align*}
Now use the fact that $h=\lambda g(t_0)$ with $\lambda$ constant in space. Every term on the right-hand side acquires exactly one factor of $\lambda$; for instance,
\begin{align*}
X(h(Y,Z))=X(\lambda g(t_0)(Y,Z))=\lambda X(g(t_0)(Y,Z)).
\end{align*}
The same constant factor appears in the remaining five terms. Therefore the right-hand side of the Koszul formula for $h$ is
\begin{align*}
\lambda\bigl(2g(t_0)(\nabla^g_XY,Z)\bigr),
\end{align*}
where the expression in parentheses is the Koszul formula for the metric $g(t_0)$. The left-hand side also contains the metric $h$, hence
\begin{align*}
2h(\nabla^h_XY,Z)=2\lambda g(t_0)(\nabla^h_XY,Z).
\end{align*}
Combining the two identities and cancelling $2\lambda$, which is allowed because $\lambda>0$, gives
\begin{align*}
g(t_0)(\nabla^h_XY,Z)=g(t_0)(\nabla^g_XY,Z).
\end{align*}
Since $Z$ was arbitrary and $g(t_0)$ is nondegenerate, the two vector fields paired against $Z$ must be equal:
\begin{align*}
\nabla^h_XY=\nabla^g_XY.
\end{align*}
This proves that multiplying a metric by a positive constant does not change its Levi-Civita connection.
Next define $R^h$ and $R^g$ to be the curvature operators of $\nabla^h$ and $\nabla^g$. For a connection $\nabla$, its curvature operator is the map determined by
\begin{align*}
R(X,Y)Z=\nabla_X\nabla_YZ-\nabla_Y\nabla_XZ-\nabla_{[X,Y]}Z.
\end{align*}
Because this expression uses only the connection and the Lie bracket, the equality $\nabla^h=\nabla^g$ implies
\begin{align*}
R^h(X,Y)Z=R^g(X,Y)Z
\end{align*}
for all smooth vector fields $X,Y,Z$. The Ricci tensor is obtained by taking the trace of the linear endomorphism $V \mapsto R(V,X)Y$. Since the underlying curvature operators are equal, their traces are equal, and hence
\begin{align*}
\operatorname{Ric}(h)(X,Y)=\operatorname{Ric}(g(t_0))(X,Y).
\end{align*}
The vector fields $X$ and $Y$ were arbitrary, so this proves the global tensor equality
\begin{align*}
\operatorname{Ric}(h)=\operatorname{Ric}(g(t_0)).
\end{align*}[/guided]