Nonlinear least squares estimations have been widely applied in positioning. However, nonlinear least squares estimations are generally biased. As the Gauss-Newton method has been widely applied to obtain a nonlinear least squares solution, we propose an iterative procedure for obtaining unbiased estimations with this method. The characteristics of the linearization error are discussed and a systematic error source of the linearization error needs to be removed to guarantee the unbiasedness. Both the geometrical condition and the statistical condition for unbiased nonlinear least squares estimations are revealed. It is shown that for long-distance observations of high precision, or for a positioning configuration with the lowest Geometric Dilution Of Precision (GDOP), the nonlinear least squares estimations tend to be unbiased; but for short-distance cases, the bias in the nonlinear least squares solution should be estimated to obtain unbiased values by removing the bias from the nonlinear least squares solution. The proposed results are verified by the Monte Carlo method and this shows that the bias in nonlinear least squares solution of short-distance distances cannot be ignored.