Use of High Performance Computing in Gravity Field Research

In the light of the three geoscientific satellite missions CHAMP, GRACE and GOCE the overall scientific aim is to achieve an automatism for the recovery of the Earth’s gravity field respectively the physical shape of the Earth, namely the geoid. Furthermore, an improved understanding of the spatial and temporal variations of the geoid is of great benefit for the study of the dynamics of the Earth’s lithosphere and upper mantle, global sea level variations, ocean circulation and ocean mass and heat transport, ice mass balance, the global water cycle and the interaction of these phenomena. This involves the determination of up to a hundred thousand unknown coefficients of the corresponding series expansion model from data sets which amount to several millions of observations provided by the satellites. The resulting system of equations which has to be solved for such an analysis cannot be evaluated without simplistic assumptions or within a satisfying time frame on personal computers due to hardware limitations. Consequently this challenging problem has to be tackled by means of high performance computing strategies. Only adoption of parallel programming standards such as MPI or OpenMP in conjunction with highly efficient numerical libraries allows for successfully accomplishing the demands of gravity field analysis. Indeed, the huge amount of data provided by satellite sensors, together with a high-resolution gravity field modeling, requires the determination of several ten thousands of unknown parameters and leads to the assignment that this problem is a true “challenge of calculus”.