[guided]The only point needing care is the sentence “constant scaling leaves Ricci unchanged.” We prove it directly. Fix $t \in I$ and write $c:=a(t)>0$, so the metric at time $t$ is $\widetilde g:=c g_0$. Let $\nabla$ be the Levi-Civita connection of $g_0$, and let $\widetilde\nabla$ be the Levi-Civita connection of $\widetilde g$. Let $\mathfrak{X}(M)$ denote the vector space of smooth vector fields on $M$.
We compare the two connections using the Koszul formula for the Levi-Civita connection. This formula applies to the Levi-Civita connection of a Riemannian metric, and both $g_0$ and $\widetilde g$ are Riemannian metrics. For smooth vector fields $X,Y,Z \in \mathfrak{X}(M)$, the formula gives
\begin{align*}
2\widetilde g(\widetilde\nabla_XY,Z)=X\widetilde g(Y,Z)+Y\widetilde g(X,Z)-Z\widetilde g(X,Y)+\widetilde g([X,Y],Z)-\widetilde g([X,Z],Y)-\widetilde g([Y,Z],X).
\end{align*}
Since $\widetilde g=cg_0$ and $c$ is a constant on $M$, every term on the right has a common factor $c$. Therefore
\begin{align*}
2\widetilde g(\widetilde\nabla_XY,Z)=c\bigl(Xg_0(Y,Z)+Yg_0(X,Z)-Zg_0(X,Y)+g_0([X,Y],Z)-g_0([X,Z],Y)-g_0([Y,Z],X)\bigr).
\end{align*}
The expression in parentheses is $2g_0(\nabla_XY,Z)$ by the Koszul formula for $g_0$. Hence
\begin{align*}
2\widetilde g(\widetilde\nabla_XY,Z)=2c\,g_0(\nabla_XY,Z)=2\widetilde g(\nabla_XY,Z).
\end{align*}
Since $\widetilde g$ is nondegenerate, equality of inner products against every $Z$ implies
\begin{align*}
\widetilde\nabla_XY=\nabla_XY.
\end{align*}
The curvature endomorphism is built only from the connection and Lie brackets:
\begin{align*}
R(X,Y)Z=\nabla_X\nabla_YZ-\nabla_Y\nabla_XZ-\nabla_{[X,Y]}Z.
\end{align*}
Because the two connections agree, their curvature endomorphisms agree:
\begin{align*}
\widetilde R(X,Y)Z=R(X,Y)Z.
\end{align*}
It remains to check that the Ricci tensor as a $(0,2)$-tensor is unchanged, not merely rescaled. Let $p \in M$, and let $n:=\dim M$. Let $(e_1,\dots,e_n)$ be a $g_0$-orthonormal basis of $T_pM$. Since $\widetilde g= c g_0$, the basis $(c^{-1/2}e_1,\dots,c^{-1/2}e_n)$ is $\widetilde g$-orthonormal. We now use the standard orthonormal-frame contraction formula in the definition of Ricci curvature. Therefore, for $X,Y \in T_pM$,
\begin{align*}
\operatorname{Ric}(\widetilde g)(X,Y)=\sum_{i=1}^n \widetilde g(\widetilde R(c^{-1/2}e_i,X)Y,c^{-1/2}e_i)=\sum_{i=1}^n c\, g_0(c^{-1/2}R(e_i,X)Y,c^{-1/2}e_i)=\sum_{i=1}^n g_0(R(e_i,X)Y,e_i)=\operatorname{Ric}(g_0)(X,Y).
\end{align*}
Thus constant positive scaling leaves the Ricci tensor unchanged as a $(0,2)$-tensor. Applying this with $\widetilde g=g(t)$ gives
\begin{align*}
\operatorname{Ric}(g(t))=\operatorname{Ric}(g_0)=\kappa g_0.
\end{align*}[/guided]