With the rapid advancement of satellite remote sensing technology, many scientists and organizations, including NASA, ESA, NAOC, and Roscosmos, observe and study significant changes in the geomagnetic field, which has greatly promoted research on the geomagnetic field and made it an important research direction in Earth system science. In traditional geomagnetic field research, tesseroid cells face degradation issues in high-latitude regions and accuracy limitations. To overcome these limitations, this paper introduces the Discrete Global Grid System (DGGS) to construct a geophysical model, achieving seamless global coverage through multi-level grid subdivision, significantly enhancing the processing capability of multi-source and multi-temporal spatial data. Addressing the challenges of the lack of analytical solutions and clear integration limits for DGGS cells, a method for constructing shape functions of arbitrary isoparametric elements is proposed based on the principle of isoparametric transformation, and the shape functions of isoparametric DGGS cells are successfully derived. In magnetic vector forwarding, considering the potential error amplification caused by Poisson’s formula, the DGGS grid is divided into six regular triangular sub-units. The triangular superconvergent point technique is adopted, and the positions of integration points and their weight coefficients are accurately determined according to symmetry rules, thereby significantly improving the calculation accuracy without increasing the computational complexity. Finally, through the forward modeling algorithm based on tiny tesseroid cells, this study comprehensively compares and analyzes the computational accuracy of the DGGS-based magnetic vector forwarding algorithm, verifying the effectiveness and superiority of the proposed method and providing new theoretical support and technical means for geophysical research.
Read full abstract