In (mathematical) physics in order to compute path integrals one often makes an infinite dimensional change of variables and uses infinite Jacobian as a purely formal expression. This step is done in formal analogy with the finite dimensional case (e.g. in the Faddeev-Popov construction).
Next one often needs to compute this Jacobian. In some situations it is a "determinant" of a self-adjoint elliptic differential operator on a compact manifold. A conventional procedure to compute it is to use the $\zeta$-function regularization. This procedure seems to me to be quite arbitrary. While I am not an expert, I have never seen in literature any motivation of it except that it gives the right answer for finite dimensional transformations. To compare with, in classical analysis there are many different methods to sum divergent series, but they may lead to different answers.
QUESTIONS. 1) Are there other conventional methods which are used in explicit computations of infinite Jacobians?
2) Are there any nice pleausible general properties of the $\zeta$-function regularized determinants which distingush them from other methods?
ADDED: This method leads to the identity mentioned in the answer below by Zurab Silagadze $$1+2+3+4+\dots=-\frac{1}{12}$$ which seems to me completely counter-intuitive.