A Gravitational Portrait of Earth from 500 Kilometers Up ✎

Original: "High-dimensional experiments for the downward continuation using the LRFMP algorithm"
arXiv:2308.04167 · 2023-08-08 · arXiv non-exclusive · independent analysis · 2 min
A new algorithm, LRFMP, processes half a million satellite measurements for the first time to recover Earth's gravitational field with a resolution of 20 kilometers.

Setting up an orchestra without a conductor: the algorithm itself decides which instruments (basis functions) and at what volume (weights) should play to reproduce Earth's gravitational symphony.

Two GRACE-FO satellites fly one behind the other at an altitude of 500 kilometers, continuously measuring how the distance between them changes. When a mountain range or an underground water reservoir passes beneath them, the gravitational field shifts slightly, and the satellites either draw closer or drift apart. In December 2022, they transmitted 529,832 such points, roughly one for every 20 kilometers along the orbit. That's like measuring gravity at intervals the length of Manhattan.

The problem is that measurements taken at altitude must be used to recover the field at the surface. This is called downward continuation, and it is exponentially ill-posed: the slightest noise in the data above turns into huge errors below, and fine details of the field are simply lost.

Previously, the LRFMP algorithm (Learning Regularized Functional Matching Pursuit) had only been tested on artificial grids of 12,500 points, and it was unclear whether it could handle real-world data volumes. The authors derived a closed-form formula for the inner products of basis functions, replacing slow numerical series. It's as if the orchestra, instead of tuning its instruments each time, received ready-made sheet music that can be played immediately. After 1,600 iterations, the relative data error was about 5%, close to the noise level, while the relative root-mean-square error at the surface was about 17.5%. Errors concentrate in the Andes, the Himalayas, and the Pacific Ring of Fire, where high-frequency details of the field are not visible from such an altitude.

Comparison with previous work shows that this is the first time the method has been scaled to a realistic amount of data. It can now compete with neural network approaches, yet it remains transparent: each selected dictionary element has a physical meaning.

This is the first test of LRFMP on real satellite geometry with half a million points, and it opens the door to practical use in geodesy. The method can be applied to monitor mass redistribution associated with melting glaciers and changes in groundwater storage. The next step is to move from already processed coefficients to raw ranging measurements, which are what the satellites directly record. Then Earth's gravitational symphony will sound even clearer.

Each GRACE-FO satellite completes approximately 496 orbits around Earth in December 2022, crossing the equator twice per orbit. That gives 1984 points at the equator, corresponding to a resolution of about 20 km: as if you were measuring gravity at intervals the length of Manhattan.

In Arthur C. Clarke's novel 'The Fountains of Paradise', engineers build a space elevator that requires precise knowledge of Earth's gravitational field. The method described here for recovering the potential from satellite data could be part of such technology.

\langle K(x,\cdot),K(\tilde{x},\cdot)\rangle_{H^2(\Omega)} = \frac{1}{64\pi}\sum_{k=0}^5 \binom{5}{k}2^k (q\frac{d}{dq})^k \phi(q)
The inner product is expressed as a finite sum of derivatives of a simple function, which speeds up computations
Scientists
Albert EinsteinRobert H. DickeFritz ZwickyJohannes KeplerSubrahmanyan Chandrasekhar
Concepts
inverse problem GPS GPCP laser ranging gravitational superposition experiments satellite imagery gps free navigation distance gravitational wave cosmology cosmic distance ladder gravimeter teleparallel gravity linear inverse model geopotential height multi-messenger astronomy cosmological simulation time reversal symmetry sky survey tianqin lsst