Imaging Earth's crustal magnetic field with satellite data: a regularized spherical triangle tessellation approach

We present a method for imaging the global crustal magnetic field at Earth's surface using a local basis representation and a minimum norm model regularization approach. The local basis consists of a spherical triangle tessellation (STT) parametrization of the radial component of the crustal field at Earth's reference spherical surface. The Green's function for Laplace's equation in spherical geometry with Neumann boundary conditions provides the necessary forward modelling scheme. We solve the inverse problem of estimating the crustal field from satellite magnetic observations by minimizing an objective function comprising a mean absolute deviation (L1-norm) measure of misfit plus a norm measuring model complexity. Both quadratic and entropy measures of field complexity are investigated. We report results from synthetic tests performed on a geophysically motivated scenario; these include a successful benchmark of the method and a comparison between quadratic and entropy regularization strategies. Applying our technique to real observations collected by the CHAMP, Ørsted and SAC-C satellites, we obtain stable images of the crustal magnetic field at Earth's surface that include sharp features with high amplitudes. We present details of two prototype crustal field models STT-CRUST-Q and STT-CRUST-E regularized using quadratic and entropy norms respectively; these provide a perspective complementary to that given by conventional spherical harmonic crustal field model

Imaging Earth's crustal magnetic field with satellite data: A regularized spherical triangle tessellation approach

ISSN:0956-540XISSN:1365-246

POGO satellite orbit corrections: an opportunity to improve the quality of the geomagnetic field measurements?

We present an attempt to improve the quality of the geomagnetic field measurements from the Polar Orbiting Geophysical Observatory (POGO) satellite missions in the late 1960s. Inaccurate satellite positions are believed to be a major source of errors for using the magnetic observations for field modelling. To improve the data, we use an iterative approach consisting of two main parts: one is a main field modelling process to obtain the radial field gradient to perturb the orbits and the other is the state-of-the-art GPS orbit modelling software BERNESE to calculate new physical orbits. We report results based on a single-day approach showing a clear increase of the data quality. That single-day approach leads, however, to undesirable orbital jumps at midnight. Furthermore, we report results obtained for a much larger data set comprising almost all of the data from the three missions. With this approach, we eliminate the orbit discontinuities at midnight but only tiny quality improvements could be achieved for geomagnetically quiet data. We believe that improvements to the data are probably still possible, but it would require the original tracking observations to be found.ISSN:1343-8832ISSN:1880-598