Regarding $T^{\mu\nu};_{\mu}=0$ in general relativity In Physics we say we have a conserved entity $\vec{P}$ (to use a common example) if we can write:
$$\frac{\partial}{\partial x^{\mu}}P^{\mu}=P^{\mu},_{\mu}=0\tag{1}$$
Where "," denotes the partial derivative. Such a formulation allows us to apply Gauss's law and say that the change of our entity within some volume is equal to the rate of flux of that entity (momentum here) through a bounding surface. This can be simply written:
$$\vec{\nabla}\cdot\vec{P}=0\tag{2}$$
The situation is changed in General Relativity, as we take covariant derivatives such that the above expression becomes:
$$P^{\mu};_{\mu}=0\tag{3}$$
Many claim that this essentially rules out a general integral conservation law in General Relativity (Moller, theory of relativity pg 337; Wheeler, gravitation and inertia pg 44; Fock space time and gravitation 163). However; a vector can be written in terms of arbitrary basis $\hat{e}$ as:
$$\vec{P}=P^{\mu}\hat{e}_{\mu}\tag{4}$$
And the covariant derivative can be simply written as:
$$\vec{P};_{\alpha}=\frac{\partial P^{\mu}}{\partial x^{\alpha}}\hat{e}_{\mu}+P^{\mu}\frac{\partial\hat{e}_{\mu}}{\partial x^{\alpha}}\tag{5}$$
Where the latter term (involving derivatives of the basis) is most commonly re-written in terms of Christoffel symbols. This expression can clearly be written in the form:
$$\left(\vec{P}\right),_{\alpha}=\left(P^{\mu}\hat{e}_{\mu}\right),_{\alpha}\tag{6}$$
The same is true of of the standard covariant expression for the stress energy tensor in General Relativity:
$$T^{\mu\nu};_{\mu}=0\tag{7}$$
Which can be written as:
$$\left(T^{\mu\nu}\hat{e}_{\mu}\hat{e}_{\nu}\right),_{\mu}=0\tag{8}$$
I don't see why this isn't a conservation law, as it only deals with the partial derivative, could someone please point me in the right direction? much appreciated
 A: Switching to invariant language solves nothing. What you write as $(\vec{P})_{,\alpha}$ is still a covariant derivative.
On a general manifold, we have no canonical way of differentiating, we need a connection. The metric provides one and we use that.
In flat space, we have a natural way of differentiating, but if we use a basis that itself changes throughout space, the component of the derivative is not the derivative of the components, because we need to factor in the rate of change of the basis fields too, hence why we talk about "covariant derivative" in flat-space context. But that is a coordinate-based phenomenon, in the invariant sense, you still use the same derivative you'd use otherwise. In particular, if in curves spacetime, you'd calculate something like $\vec{P}_{[,\mu,\nu]}$, it would give nonzero result (provided you do it properly and contract the derivative-index with a basis vector one you do the first derivative). (Mathematicians would call this $\nabla^2_{X,Y}Z-\nabla^2_{Y,X}Z=R(X,Y)Z$.)
This doesn't change the fact that the main problem is that a vector $v_x\in T_xM$ at $x$ and a vector $v_y\in T_yM$ live in different vector spaces, operations between them don't make mathematical sense.
So if you'd calculate the integral of a vectorial (or tensorial) object, the operation $$ \int_M\vec{P}\ d\mu $$ would correspond to $$ \lim_{\Delta x^\mu\rightarrow 0}\sum_{i}\vec{P}(x_i)\sqrt{-g}(x_i)\Delta x^0\cdot...\cdot\Delta x^3, $$ where the index $i$ labels a partition of the manifold $M$ into small coordinate-rectangles, and it is impossible to add up $\vec{P}(x_i)$ and $\vec{P}(x_j)$ for $i\neq j$, since they are located in different spaces. No vector-valued integrals in differential geometry, sorry.
A: *

*The problem seems to be OP's eq. (5), which incorrectly uses partial instead of covariant derivatives on the right-hand side. 

*If we write a $(r,s)$ tensor 
$$T^{\mu_1\ldots\mu_r}{}_{\nu_1\ldots\nu_s}$$ in index-free notation
$$T~=~ \frac{\partial}{\partial x^{\mu_1}}\otimes \cdots \otimes \frac{\partial}{\partial x^{\mu_r}} ~T^{\mu_1\ldots\mu_r}{}_{\nu_1\ldots\nu_s} ~ dx^{\nu_1}\otimes \cdots \otimes dx^{\nu_s},$$
it doesn't mean that covariant derivatives may be replaced with partial derivatives. 

*Concerning how covariant derivatives act on basis elements, the reader may find my Phys.SE answer here interesting.
A: This isn't the root reason why you don't, in general get integral forms of conservation equations in general relativity.  The root reason why is that integral conservation laws set a time derivative of some quantity equal to an integral of a field.  
But, in a fully covariant spacetime, you can mess this up by choosing a crazy time coordinate. 
So, if you're going to use the Gauss's law trick, and set the time derivative of some charge equal to the integral over some timelike or null stack of two-surfaces that bound a stack of 3-spaces, the timelike or null derivative of the charge will have to be an isometry of the 3-space, so that you don't get any geometrically induced weirdness.  
At asymptotic infinity, any timelike vector will do, because you're asymptoting to Minkowski space.  But this isn't the only case where this is true.  The event horizon of an isolated black hole also has these properties, with the null tangent providing the role of "time."
