It's $0, 012, 012 120 201, 012 120 201 120 201 012 201 012 120, ...$
To easily construct the Thue-Morse sequence, you can do the following:
Start with $0$;
Replace $(0\to01)$ and $(1\to10)$.
So you get: $$0, 0 1, 01 10, 0110 1001, 01101001 10010110,...$$
To easily construct the ternary version of Thue-Morse sequence, you can do the following:
Start with $0$;
Replace $(0\to012), (1\to120), (2\to201).$
So you get: $$0, 0 1 2, 012 120 201, 012120201 120201012 201012120$$
To easily construct the $4$-ary version of Thue-Morse sequence, you can do the following:
Start with $0$;
Replace $(0\to 0123), (1\to 1230), (2\to2301), (3\to3012).$
So you get: $$0, 0 1 2 3, 0123 1230 2301 3012, 0123123023013012 1230230130120123 2301301201231230 3012012312302301$$
To easily construct the $n$-ary version of Thue-Morse sequence, you can do the following:
Start with $0$;
Replace
$$(0\to012\cdots[n-2][n-1]n),
(1\to123\cdots[n-1]n0),$$
$$(2\to234\cdots n01),$$
$$\dots$$
$$(n\to n01\cdots[n-3][n-2][n-1]).$$
There's probably no "closed form" for any $s>1$. Curiously, though,
the sum can be evaluated for integers $s \leq 0$, in the following sense.
The function $\zeta_{TM}$ extends to an analytic function on
${\bf C} \backslash \{ 1 \}$, with a simple pole at $s=1$
of residue $1/2$, and taking rational values at integers $s \leq 0$,
starting
$\zeta_{TM}(0) = -1/4$,
$\zeta_{TM}(-1) = -1/24$,
$\zeta_{TM}(-2) = 0$,
$\zeta_{TM}(-3) = +1/240$,
and in general
$\zeta_{TM}(s) = \zeta(s) / 2$ for integers $s \leq 0$
(so in particular $\zeta_{TM}$ inherits the "trivial zeros" of $\zeta$ at
$s = -2, -4, -6, \ldots$).
It is more convenient to work with the Dirichlet series whose
$(n+1)^{-s}$ coefficient is $1 - 2 t_n = (-1)^{t_n}$, because
the generating function for $(-1)^{t_n}$, call it
$$
T(z) = \sum_{n=0}^\infty (-1)^{t_n} z^n,
$$
factors as an infinite product:
$$
T(z) = (1-z) (1-z^2) (1-z^4) (1-z^8) \cdots
= \prod_{m=0}^\infty \bigl( 1 - z^{2^m} \bigr).
$$
So define
$$
Z_{TM}(s) = \zeta(s) - 2 \zeta_{TM}(s)
= \sum_{n=0}^\infty \frac{(-1)^{t_n}}{(n+1)^s}.
$$
The usual Mellin-transform trick gives an integral formula:
$$
\Gamma(s) Z_{TM}(s)
= \sum_{n=0}^\infty (-1)^{t_n} \! \int_0^\infty x^{s-1} e^{-(n+1)x} \, dx
= \int_0^\infty x^{s-1} e^{-x} T(e^{-x}) \, dx.
$$
This gives an analytic continuation of $\Gamma(s) Z_{TM}(s)$ to the entire
complex plane, because $T(e^{-x})$ decays faster than any power of $x$
as $x \to 0$: each factor $1 - e^{-2^m x}$ of the infinite product
is $O_m(x)$ and in $(0,1)$. Since $\Gamma(s)$ has no zeros,
but does have simple poles at $s = 0, -1, -2, -3, \ldots$,
it follows that $Z_{TM}$ is an entire function with simple zeros
at the same $s$, and no other real zeros (the integral for
$\Gamma(s) Z_{TM}(s)$ is plainly positive for all real $s$).
Since $Z_{TM} = \zeta - 2 \zeta_{TM}$, we conclude that
$\zeta_{TM}(s) = \frac12 \zeta(s)$ at those $s$, as claimed.
[The integral formula can also be used to compute
$Z_{TM}(s)$, and thus also $\zeta_{TM}(s)$, to high precision;
for example, using gp's "intnum" function we find
$Z_{TM}(2) = 0.6931534522\ldots$ (this is not $\log 2$,
though it's quite close -- the difference is $\lt 10^{-5}$),
so $\zeta_{TM}(2) = (\zeta(2) - Z_{TM}(2)) / 2 = 0.4758903073\ldots$.]
Best Answer
The answer is of course yes, see the MaothOverflow post here.