Euler's Choice:
When Euler crafted the zeta function, he knew that $\zeta(1)$ diverged, so he made $\zeta(1)$ undefined.
When he crafted the zeta generating function using the Bernoulli numbers, he allowed that function to return $-\frac{1}{2}$ for $\zeta(0)$. I believe he could foresee no consequences in doing this.
Riemann's Oversight:
When Riemann modified the zeta function to handle the complex numbers, he defined $\zeta(1)$ to return ComplexInfinity. This would allow the reciprocal to return a $0$ to match the $0$ returned by $\sum_{n=1}^{\infty}\frac{\mu(n)}{n^s}$ when $s=1$. His oversight was to not fix $\zeta(0)$ to also return Complex Infinity for the same reason as he did for $\zeta(1)$.
My Testing:
I modified Mathematica to tweak $\zeta(0)$:
Unprotect[Zeta];
Zeta[0] := ComplexInfinity;
Protect[Zeta];
Replace the middle line with:
Zeta[0]=.;
to restore Zeta back to the original.
So far, my testing has not broken anything. However, we all know that zeta is found in many functions, so we should do exhaustive testing.
A Conjecture:
I think that the Critical Line being on Re $\frac{1}{2}$ is a consequence of a prime being a factor of exactly $\frac{1}{2}$ of the square-free numbers.
The Question:
What do you think the consequences would be if the definition of $\zeta(0)$ is changed?