Regularization parameter estimation for underdetermined problems by the χ 2 principle with application to 2D focusing gravity inversion