IOSVRM method for earth gravity field estimation using GOCE satel-lite gravimetric observations and comparison with Tikhonov and TSVD methods

نوع مقاله : مقاله پژوهشی‌

نویسندگان

1 Assistant Professor, Department of Geomatics, Thechnical School of Marand, University of Tabrız, Tabrız, Iran

2 Assistant Professor, Department of Civil Engineering, Engineering Faculty, Behbahan Khatam Alanbia University of Technology, Behbahan, Iran

چکیده

The ill-posed problems could be seen anywhere in our daily lives. An ill-posed problem is a problem that there are no uniqueness solutions (there is no solution or two or more solutions for the problem) or the solutions are unstable; i.e., an arbitrarily small error in the observation may lead to extremely large errors in the solutions. The main difficulty in solving ill-posed problems is instability of their solutions with respect to small variations of input data. A regularized estimation of an ill-posed problem is always biased; thus, it's worth obtaining the solution from different methods for reliable evaluation of the uncertainty in our
estimation. The regularization methods, such as Tikhonov's method, are used to obtain stable solutions for solving ill-posed problems. In Tikhonov's regularization method, a scalar
quantity is used as the stabilization parameter to solve ill-posed problems; whereas, In the Optimally Scaled Vector Regularization Method (OSVRM), a vector is used as the
stabilization parameter. In this paper, a comparison has been made between the results of Tikhonov, TSVD, and OSVRM methods in terms of accuracy of the results for estimation of the earth gravity field from GOCE satellite data. The RMSE of the results of Tikhonov, EGM96 model, and the OSVRM method – that use a vector instead of scalar as
regularization parameter – in the order of 10-5, 10-9, and 10-13, respectively. It is seen that the results obtained from the OSVRM method are much better compared to the Tikhonov method and EGM96 model for solving linear ill-posed problems. On the other hand, a
significant improvement has been achieved in the stability and accuracy of numerical results for linear problems solution.
 
 

کلیدواژه‌ها


عنوان مقاله [English]

IOSVRM method for earth gravity field estimation using GOCE satel-lite gravimetric observations and comparison with Tikhonov and TSVD methods

نویسندگان [English]

  • Abolfazl Shahamat 1
  • Lotfollah Emadali 2
1 Assistant Professor, Department of Geomatics, Thechnical School of Marand, University of Tabrız, Tabrız, Iran
2 Assistant Professor, Department of Civil Engineering, Engineering Faculty, Behbahan Khatam Alanbia University of Technology, Behbahan, Iran
چکیده [English]

The ill-posed problems could be seen anywhere in our daily lives. An ill-posed problem is a problem that there are no uniqueness solutions (there is no solution or two or more solutions for the problem) or the solutions are unstable; i.e., an arbitrarily small error in the observation may lead to extremely large errors in the solutions. The main difficulty in solving ill-posed problems is instability of their solutions with respect to small variations of input data. A regularized estimation of an ill-posed problem is always biased; thus, it's worth obtaining the solution from different methods for reliable evaluation of the uncertainty in our
estimation. The regularization methods, such as Tikhonov's method, are used to obtain stable solutions for solving ill-posed problems. In Tikhonov's regularization method, a scalar
quantity is used as the stabilization parameter to solve ill-posed problems; whereas, In the Optimally Scaled Vector Regularization Method (OSVRM), a vector is used as the
stabilization parameter. In this paper, a comparison has been made between the results of Tikhonov, TSVD, and OSVRM methods in terms of accuracy of the results for estimation of the earth gravity field from GOCE satellite data. The RMSE of the results of Tikhonov, EGM96 model, and the OSVRM method – that use a vector instead of scalar as
regularization parameter – in the order of 10-5, 10-9, and 10-13, respectively. It is seen that the results obtained from the OSVRM method are much better compared to the Tikhonov method and EGM96 model for solving linear ill-posed problems. On the other hand, a
significant improvement has been achieved in the stability and accuracy of numerical results for linear problems solution.
 
 

کلیدواژه‌ها [English]

  • Tikhonov stabilization method
  • regularization parameter
  • singular value decomposition
  • OSVRM method
  • GOCE satellite
  1. Ardalan, A. A., Safari, A. and Allahtavakoli, Y., 2008, On the optimum method for estimation of regularization parameter of downward continuation in the problem of geoid computation without Stokes formula,  Journal Of The Earth And Space Physics, Vol. 34 (3), pages 57-78

    Chein-Shan Liu (2012); “Optimally scaled vector regularization method to solve ill-posed linear problems”; Applied Mathematics and Computation, Volume 218, Issue 21, Pages 10602-10616

    Hadamard, J., 1923, Lectures on Cauchy’s Problem in Linear Partial Differential Equations. Yale University Press.

    Hansen, P. C., 1992, Numerical tools for analysis and solution of Fredholm integral equations of first kind, Inverse Problems 8:849-872., Printed in the UK.

    Per Christian Hansen And Dianne Prost O’leary, 1993, The Use Of The L-Curve In The Regularization Of Discrete Ill-Posed Problems; SIAM j. Sci. Comput. Vol. 14, No. 6, pp. 1487-1503, November 1993

    1. Kusche, J., Klees, R. (2002): “Regularization of gravity field estimation from satellite gravity gradients”; Journal of Geodesy, volume76, pages359–368

    Liu, C.-S., Hong, H. K., Atluri, S. N., 2010, Novel algorithms based on the conjugate gradient method for inverting ill-conditioned matrices and a new regularization method to solve ill-posed linear systems, CMES: Comput. Model. Eng. Sci. 60, 279–308.

    Liu, C.-S., 2011, A revision of relaxed steepest descent method from the dynamics on an invariant manifold, CMES: Comput. Model. Eng. Sci. 80, 57–86

    M.S. Petrovskaya · A.N. Vershkov (2006): “Non-singular expressions for the gravity gradients in the local north-oriented and orbital reference frames”; Journal of Geodesy, volume 80, pages117–127.