Well, two and a half years later, I stumbled across an answer. The following is a paraphrasing of Petersen's "Riemannian Geometry," page 40.
We work on a Riemannian manifold $(M,\langle \cdot, \cdot \rangle)$ with Levi-Civita connection $\nabla$. Let $\text{Ric}$ denote the symmetric $(1,1)$-tensor $\text{Ric}(v) = \sum R(v,e_i)e_i$, where $\{e_i\}$ is an orthonormal basis. Let $S$ denote the scalar curvature, i.e. $S = \text{tr}(\text{Ric}) = \sum \langle \text{Ric}(e_i), e_i \rangle$. We note that $\nabla S = dS$.
Prop: $dS = 2\,\text{div}(\text{Ric})$.
Proof: Let $\{E_i\}$ be an orthonormal frame at $p \in M$ with $\nabla E_i|_p = 0$. Let $X$ be a vector field with $\nabla X |_p = 0$. From the second Bianchi identity:
$$\begin{align*}
dS(X)(p) = XS(p) & = X \sum \langle \text{Ric}(E_i), E_i \rangle \\
& = X \sum \langle R(E_i, E_j)E_j, E_i\rangle \\
& = \sum \langle \nabla_X [R(E_i, E_j)E_j], E_i \rangle \\
& = \sum \langle (\nabla_X R)(E_i, E_j)E_j, E_i \rangle \\
& = -\sum \langle (\nabla_{E_j}R)(X, E_i)E_j, E_i \rangle - \sum \langle (\nabla_{E_i}R)(E_j,X)E_j, E_i \rangle \\
& = -\sum (\nabla_{E_j}R)(X, E_i, E_j, E_i) - \sum (\nabla_{E_i}R)(E_j, X, E_j, E_i) \\
& = \sum (\nabla_{E_j}R)(E_j, E_i, E_i, X) + \sum (\nabla_{E_j} R)(E_i, E_j, E_j, X) \\
& = 2 \sum (\nabla_{E_j}R)(E_j, E_i, E_i, X) \\
& = 2 \sum \nabla_{E_j}(R(E_j, E_i, E_i, X) \\
& = 2 \sum \nabla_{E_j} \langle \text{Ric}(E_j), X \rangle \\
& = 2 \sum \nabla_{E_j} \langle \text{Ric}(X), E_j \rangle \\
& = 2 \sum \langle \nabla_{E_j}(\text{Ric}(X)), E_j) \rangle \\
& = 2 \sum \langle (\nabla_{E_j}\text{Ric})(X), E_j \rangle \\
& = 2\,\text{div}(\text{Ric})(X)_p.
\end{align*}$$
There seems to be some disagreement in the literature about whether we assume $\lambda$ is constant from the outset when we write $\text{Ric}_{ij} = \lambda g_{ij}$. Here is a hint to a proof that $\lambda$ is constant using the classical second Bianchi identity
$$R_{ijk\ell;m}+R_{ij\ell m;k}+R_{ijmk;\ell}=0\,.$$
Working with $g_{ij}(p) = \delta_{ij}$, consider $\sum\limits_{i,k,\ell} \delta_{k\ell} R_{iki\ell}$.
Best Answer
Here's the idea: define a covariant $5$-tensor $T$ by $$T_{ijklm} = R_{ijkl;m} + R_{ijlm;k} + R_{ijmk;l}.$$ The differential Bianchi identity says that $T$ is the zero tensor. First contract $T$ on its first and fourth indices, yielding a $3$-tensor $U$ which is also zero: $$0 = U_{jkm}= T_{ijk}{}^i{}_m = R_{jk;m} - R_{jm;k} + R_{ijmk;}{}^{i}.$$ Then contract $U$ on its first two indices: $$0 = U_{j}{}^j{}_m = S_{;m} - R_{jm;}{}^{j} - R_{im;}{}^{i}.$$