Suppose that $\mathbf{b}:M\to S^2$ is given and let $B\subset M\times S^2\times S^2\times S^2$ consist of the oriented orthonormal frames $(x,e_1,e_2,e_3)$ that satisfy $e_3 = \mathbf{b}(x)$. Then the projection $x:B\to M = \mathbb{R}^3$ makes $B$ into a circle bundle over $M$, and we have the usual first structure equations
$$
\mathrm{d}x = e_i\,\omega_i = e_1\,\omega_1 + e_2\,\omega_2 + e_3\,\omega_3
\qquad\text{and}\qquad
\mathrm{d}e_i = e_j\,\omega_{ji}\,,
$$
where $\omega_{ij}=-\omega_{ji}$ and the $\omega_i$ satisfy the second structure equations
$$
\mathrm{d}\omega_i = -\omega_{ij}\wedge\omega_j
\qquad\text{and}\qquad
\mathrm{d}\omega_{ij} = - \omega_{ik}\wedge\omega_{kj}\,.
$$
By construction, the $1$-forms $\omega_1,\omega_2,\omega_3,\omega_{12}$ are a basis for the $1$-forms on $B$, while we must have functions $a_{ir}$ (where $i=1,2$ and $r=1,2,3$) such that
$$
\omega_{i3} = a_{ir}\,\omega_r\,
$$
Now, the choice of a frame field $(\mathbf{t},\mathbf{n},\mathbf{b}):M\to S^2\times S^2\times S^2$ with the above properties, i.e.,
$$
\mathrm{d}\mathbf{t}\equiv 0
\quad \text{modulo}\quad \mathbf{n},\omega_2,\omega_3\,,
$$
is really a choice of a section of $B$ that falls into the locus where
$$
\omega_{31} = -\omega_{13} = -a_{1r}\,\omega_r \equiv -a_{11}\,\omega_1 \equiv 0
\quad\text{modulo}\quad \omega_2,\omega_3\,,
$$
i.e, in the locus where $a_{11}=0$. Thus, we need to know how $a_{11}$ varies on the fibers of $B\to M$ in order to understand its zero locus.
Thus, whether there exists a foliation of the desired kind depends at least on a certain inequality in the first derivatives of the $2$-plane field being satisfied.
Added Remark: A 'geometric' interpretation of this obstruction can be made as follows: On $B$ (which is determined by the unit vector field $\mathbf{b}:M\to S^2$) quadratic form $$ Q = a_{11}\,{\omega_1}^2 + (a_{12}{+}a_{21})\,\omega_1{\circ}\omega_2 + a_{22}\,{\omega_2}^2 $$ turns out to be the pullback under the submersion $x:B\to M$ of a quadratic form $Q_{\mathbf{b}}$ on the $2$-plane field $D = \mathbf{b}^\perp\subset TM$, and what one finds is that, if $\mathbf{t}:M\to S^2$ is actually a section of $D$ such that its integral curves have their principal normals also lying in $D$, then $\mathbf{t}$ must be a null vector for the quadratic form $Q_{\mathbf{b}}$. Thus, if $Q_{\mathbf{b}}$ is definite (positive or negative), then there is no solution, while if $Q_{\mathbf{b}}$ is of split type, there are two distinct null directions at each point, and these define the two possible line fields with the desired property of their principal normals.
The interesting case is what happens when $Q_{\mathbf{b}}$ vanishes identically, which, of course, can happen, but is a very special case (in particular, the plane field $D$ has to be Frobenius, but even this is not sufficient). Then there is much more that can happen with the foliations.