 Solution

A steady state ${\mathbf{c}}_{\ast }$ satisfies $0=A{\mathbf{c}}_{\ast }-\mathbf{b}$, henceThe sign is positive on the right after rearranging the supplied equation. A solution exists exactly when $\mathbf{b}$ lies in the column space of $A$, which above was found to be the zero-sum vector subspace. Equivalently, multiply by the left null space vector $1^T$:This condition is necessary and sufficient for an algebraic steady state, since both the column space and that zero-sum vector subspace have dimension $3$. It expresses global conservation of mass: in the closed box model, total tracer production must balance total tracer removal. With the paper's convention $-\mathbf{b}$, a positive $b_i$ is a sink and a negative $b_i$ is a source. If the sum is nonzero, $dM/dt=-{\sum }_i{b}_i$ and the inventory changes linearly in time, preventing a steady state.

When the condition holds, all algebraic steady states form the affine solution space of a linear equation ${\mathbf{c}}_p+\alpha 1$, for any particular solution $A{\mathbf{c}}_p=\mathbf{b}$. The initial inventory selects one member through $1^T{\mathbf{c}}_{\ast }=1^T\mathbf{c}(0)$. For a physically admissible solution, its four concentrations must additionally be nonnegative; the compatibility condition alone imposes no such inequalities.