There is no contradiction in the three answers you are citing, just loose language. Only the third answer deals with your spare globally symmetric Lagrangian. The other two deal with QCD, its gauged (local) extension, involving gluons, which most texts write down: it involves an extra fermion bilinear/ gluon linear term, and a kinetic term for the gluons. It too is also globally, beyond being locally invariant, but the first two answers fuss the analog Noether's theorem for gauge theories which the preamble to your question excludes.
Noether's (first) theorem deals with the eight Lie-algebraic generators $F^\alpha, ~~~~\alpha=1,2,...,8$, which you may think of as eight linearly independent 3×3 matrices acting on the three a labels of the states you wrote down. (I assume you are not interested in the specific basis thereof, etc.),
$$
\psi_a\to \Bigl (\exp (i\theta^\alpha F_\alpha)\Bigr )_{ab} \psi_b .
$$
Summation over repeated indices is implied.
There is one "color rotation angle" parameter $\theta^\alpha$ for each of the eight generators written. For infinitesimally small such, the generic increment of the 3-spinors your wrote down is
$$
\delta \psi_a = i\theta^\alpha (F_\alpha)_{ab} \psi_b.
$$
These yield eight "color charges" through Noether's (first) theorem,
$$
Q^\alpha =i\int\!\! dx~~ \bar \psi_a (F^\alpha)_{ab} \gamma_0\psi_b,
$$
which you evidently checked are on-shell (~through the application of the equations of motion) time-conserved.
Commuting these eight charges with the 3-spinors yields the above $\delta \psi$ increments. Indeed, the two charges corresponding to the two Cartan generators, together with the Baryon charge, $\propto i\int\!\! dx~~ \bar \psi_a \gamma_0\psi_a$, also commuting with those, suffice to uniquely specify the independent components of each triplet $\psi_a$, some people call r,g,b for convenience. So: three colors, eight independent rotations thereof (conserved charges).
If/when you add gluon gauge fields, which transform among themselves through 8×8 (not 3×3) matrices, the expanded Lagrangian (QCD) is also globally SU(3) invariant, but through the magic of the gauge construction, it is, in addition, locally (gauge) invariant an issue outside the ambit of your question, if not your confusion.
The way the 8 Lie algebra generators and corresponding color charges compose among themselves (commutation) is best described through two roots, corresponding to two Cartan generators; you may think about them as commutator motions among these 8 charges.
So, in a way, your three questions harmonize, after all: it's just that they describe three different aspects of the elephant to the three blind men touching three different parts of her...
Clarification for comment question: Yes it does, only up to language. Algebraically, flavor SU(3) and color SU(3) are identical, so one uses labels and matrices of the two interchangeably, trusting the reader cannot be confused.
Since the same generators are involved, one automatically maps the three diagonal matrices' $\lambda_8; \lambda_3; {\mathbb I}$ diagonals, suitably normalized, to the state vectors $r=(1,0,0); g=(0,1,0); b=(0,0,1)$, virtually a change of basis. The idea is people learn their SU(3) from flavor (the eightfold way), and then instantly translate to color. As generators, they are conserved charges (the baryon number is the identity in that space of triplets), but not all conserved color charges.