 Frontier letter
 Open Access
 Published:
Inflight scalar calibration and characterisation of the Swarm magnetometry package
Earth, Planets and Space volume 68, Article number: 129 (2016)
Abstract
We present the inflight scalar calibration and characterisation of the Swarm magnetometry package consisting of the absolute scalar magnetometer, the vector magnetometer, and the spacecraft structure supporting the instruments. A significant improvement in the scalar residuals between the pairs of magnetometers is demonstrated, confirming the high performance of these instruments. The results presented here, including the characterisation of a Sundriven disturbance field, form the basis of the correction of the magnetic vector measurements from Swarm which is applied to the Swarm Level 1b magnetic data.
Introduction
In November 2013 the European Space Agency (ESA) launched the three Swarm satellites, named Alpha, Bravo, and Charlie, with the objective to provide the best ever survey of the geomagnetic field and its temporal evolution (FriisChristensen et al. 2006). Each spacecraft carries an Absolute Scalar Magnetometer (ASM) for measuring Earth’s magnetic field intensity, a Vector Fluxgate Magnetometer (VFM) measuring the direction and strength of the magnetic field, and a threehead Star TRacker (STR) mounted close to the VFM to obtain the attitude needed to transform the vector readings to an Earthfixed coordinate frame. Time and position are provided by an onboard GPS receiver. The payload also includes instruments to measure plasma and electric field parameters as well as nongravitational acceleration.
One of the purposes of the scalar magnetometer (ASM) is to provide the necessary absolute magnetic data to calibrate the vector magnetometer (VFM). For this an approach similar to that adopted for the previous satellite missions Ørsted and CHAMP was foreseen (c.f. Olsen 2003; Yin and Lühr 2011) since those missions carried equivalent instrumentation. However, soon after launch of Swarm it became clear that the magnetic field vector measurements on all three spacecraft were contaminated by unforeseen disturbances which could not be captured by the traditional inflight calibration methods referred to above. Furthermore, the disturbances show systematic variation which could impact or map into scientific investigations based on Swarm magnetic data. The light blue symbols in Fig. 1 show time series of the scalar residuals, which are the difference, \(\Delta F = \vec {B}_\mathrm {VFM}  F_\mathrm {ASM}\), between the modulus of the VFM data, \(\vec {B}_\mathrm {VFM}\), and the magnetic intensity measurements, \(F_\mathrm {ASM}\), taken by the ASM instrument. Based on experience with Ørsted and CHAMP scalar residuals with subnanotesla level were expected (rms value well below 0.5 nT), while for Swarm the scatter of the residuals was observed to reach several nT, resulting in an rms value approaching 1 nT, but crucially showing a very clear local time dependence. A task force was therefore established to investigate and mitigate the effect.
Detailed investigations of the scalar residuals \(\Delta F\) and of the ASM and VFM measurements separately indicated that:

the vector readings of the VFM are affected by a disturbance vector field;

the scalar readings of the ASM are much less, if at all, affected.
Consequently, the task force concluded to pursue models which assume the magnetic disturbance to be affecting the VFM measurements only. Plotting \(\Delta F\) as a function of the Sun incidence angles with respect to the spacecraft, reveals systematic features of the disturbance, as shown in Fig. 2. At the start of “Characterisation and calibration with scalar residuals” section we provide detailed definitions of the two Sun incidence angles \(\alpha\) and \(\beta\). This supports the hypothesis that a magnetic source in the vicinity of the VFM magnetometer, with strength and direction depending on the direction to the Sun (as seen from the spacecraft), is responsible. We refer to such a disturbance field vector that depends on the direction to the Sun, as \(\delta \vec {B}_\mathrm {Sun}\).
The purpose of this article is to document the details of inflight calibration of the Swarm magnetometer package, including an empirical determination and removal of the Sundriven vector disturbance field \(\delta \vec {B}_\mathrm {Sun}\), based on a mitigation approach proposed by Vincent Lesur (Lesur et al. 2015).
“Characterisation and calibration with scalar residuals” section describes the parameterisation of the model of the Sundriven disturbance—in following referred to as the characterisation of the disturbance field—and of the calibration of the VFM instrument, by which means determination of its intrinsic scale factors and their dependence on time and temperature, and determination of the sensoraxis nonorthogonalities. We document the adopted Iteratively Reweighted Least Squared (IRLS) estimation approach that includes a truncated singular value decomposition (SVD) approach to solving the inverse problem. The results obtained for Swarm Alpha, based on data covering the period from launch (22 November 2013) until end of June 2015 (i.e. 19 months), are presented in “Results of model estimation for Swarm Alpha ” section. Application of the scheme to data from the satellites Bravo and Charlie resulted in similar levels of residual improvement and statistics, and the estimates of the Sundriven disturbance \(\delta \vec {B}_\mathrm {Sun}\) show generally similar behaviour and structural features as found for Swarm Alpha, although there are also some differences. Finally, “Conclusions” section summarises the findings and provides perspectives regarding further improvements of the method.
Characterisation and calibration with scalar residuals
The Sun incidence angles \(\alpha\) and \(\beta\) are crucial in our approach to characterise the scalar residual. To clarify, in Fig. 3 we illustrate the definition of these angles with respect to the spacecraft and the Sun position. \(\alpha\) is the azimuth in the spacecraft x–z plane (nominally the orbit plane), and \(\beta\) is the “elevation” out of the x–z plane positive towards left (looking in the nominal flight direction; i.e. positive opposite the spacecraft yaxis). Examples of values for \(\alpha\) and \(\beta\) for particular Sun positions are:

\(\beta = +90^\circ\): Sun directly from \(y\) (i.e. from the left during nominal flight)

\(\beta = 90^\circ\): Sun directly from \(+y\) (i.e. from the right)

\(\beta = 0^\circ\), \(\alpha = 0^\circ\): Sun directly from \(+x\) (i.e. from the front)

\(\beta = 0^\circ\), \(\alpha = +90^\circ\): Sun directly from \(z\) (above)

\(\beta = 0^\circ\), \(\alpha = +180^\circ\): Sun directly from \(x\) (i.e. from the back—slightly above the boom)
Considering how these angles vary over orbits of the Swarm spacecraft during nominal flight, we find that \(\alpha\) varies rapidly from \(360^\circ\) down to \(0^\circ\) within one orbit (i.e. within \({\approx }90\) min), while \(\beta\) varies slowly up and down typically by \({\approx }1.25^\circ\) in one day (for Alpha and Charlie, \(1.20^\circ\) for Bravo).
Although the observed scalar residuals clearly vary with the Sun incidence angles \(\alpha\) and \(\beta\) (see Fig. 2), there is no direct mapping of \(\Delta F\) in terms of these parameters. This is a consequence of the scalar residuals \(\Delta F \approx \delta \vec {B}_\mathrm {Sun} \cdot \vec {b}_0\) being the projection of the magnetic disturbance vector \(\delta \vec {B}_\mathrm {Sun}\), onto the unit vector \(\vec {b}_0\) of the ambient magnetic field direction (Earth’s main field). The former is oriented relative to the spacecraft, while the latter is oriented relative to Earth, which results in the variations with the spacecraft local time (captured by \(\beta\)) as seen in Fig. 2. The spacecraft local time changes by 12 hours (corresponding to a change in \(\beta\) by \(180^\circ\)) within approximately \(4\frac{1}{2}\) months.
To account for the projection on to the ambient field, we consider a vector magnetic disturbance \(\delta \vec {B}_\mathrm {Sun}(\alpha , \beta )\), with each component depending individually on the Sun incidence angles. Mathematically, we describe each component of the disturbance field vector by a spherical harmonic expansion in \(\alpha\) and \(\beta\) i.e. we consider three independent spherical harmonic expansions in all.
This model characterising the Sundriven disturbance is coestimated together with a model of the temporal evolution of the VFM sensitivity and an adjustment of the preflight estimated nonorthogonality angles of the VFM sensor. For this we perform a scalar calibration via a least squares fit, minimising the discrepancy (\(\Delta F\)) between the fully calibrated and corrected measurements from the ASM and the modulus of the vector measurements from the VFM after our model has been applied. Huber weights are used iteratively to eliminate the effect of anomalous measurements (“outliers”) on the estimated models.
Model parameterisation
As outlined above, our model characterising the Sundriven disturbance vector \(\delta \vec {B}_\mathrm {Sun}\) consists of three spherical harmonic expansions up to degree and order 25, one for each of the magnetic field components in the VFM magnetometer frame, with the position of the Sun with respect to the spacecraft parameterised by the Sun incidence angles \(\alpha\) and \(\beta\). It takes the form
where \(\vec {u}_n^m\) and \(\vec {v}_n^m\) are the spherical harmonic expansion coefficients, with one component for each component of the disturbance field, and \(P_n^m\) are the Schmidt seminormalised Legendre functions. Note that \(\delta \vec {B}_\mathrm {Sun}\) includes static terms (\(n = m = 0\)) that describe a static (i.e. independent of the Sun position) disturbance vector. The disturbance field vector \(\delta \vec {B}_\mathrm {Sun}\) is thus described by \(3 \times 26^2 = 2028\) model coefficients.
The model for rescaling the vector measurements and taking into account any small adjustment of the nonorthogonality of the VFM sensors, which is required in order to obtain the fully calibrated and corrected vector field measurements \(\vec {B}_\mathrm {VFM}\), now takes the form
where \(\vec {B}_{{\rm pre}\text{}{\rm flight}}\) are the VFM measurements calibrated using the preflight parameters and corrected for the preflight determined stray fields as described in TøffnerClausen (2015). \(\underline{\underline{S}}\) is a \(3 \times 3\) diagonal scaling matrix with elements
where \(s^{{\rm B}\text{}{\rm spline}}(t)\) is a quadratic Bspline in time with 3month knot separation (common for all three components of the magnetic field), and \(s_{j,\mathrm {Tsensor}}, j=13\) is an adjustment of the preflight estimated dependency of the VFM sensitivity on its sensor temperature, \(T_\mathrm {sensor}\), for each sensor axis j. \(s_{j,\beta }\) is an empirical scaling parameter and \(\beta\) the Sun incidence angle, as defined above. The choice of quadratic Bsplines with 3month knot separation is made to allow sufficient flexibility of the model; the exact choice of Bspline knot times is not crucial as very similar results are obtained with other, similar parameterisations. The estimated Bsplines exhibit very moderate accelerations (in the case of the full model, see Fig. 6), and it may be possible to simplify the parameterisation of the time dependence in future models, e.g. to an exponential saturation in time as this is the expected behaviour of the VFM instrument sensitivity; however, an exponential model is illconditioned on the time span of data used here.
\(\underline{\underline{P}}\) is the nonorthogonality matrix that makes small adjustments to the preflight estimated nonorthogonalities of the VFM sensor (cf. Olsen 2003)
Our inflight calibration model comprises 18 parameters in all; together with the 2028 parameters describing \(\delta \vec {B}_\mathrm {Sun}\) this results in 2046 model parameters to be estimated, as listed in Table 1.
Estimation of model parameters: inversion and regularisation
In order to estimate the 2046 model parameters from the scalar residuals we need to solve a nonlinear inverse problem. The nonlinearity arises from the treatment of nonorthogonalities (Olsen 2003).
The forward relationship between the vector of the scalar residuals, \(\mathbf {d}\) (\(d_i = \Delta F_i\), the scalar residual of the ith data point), and the model parameter vector \(\mathbf {m}\), may therefore be written in the form
where \(\mathbf {g}(\mathbf {m})\) is a nonlinear function of the model parameters and \(\mathbf {e}\) is a small remainder that cannot be explained by the model, which we seek to minimise.
Linearisation of this problem is straightforward. A regularised, iteratively reweighted, least squares solution to the inverse problem is then obtained using the algorithm
where at the kth iteration, \(\underline{\underline{\mathbf {G}}}_k= \left. {\partial \mathbf {g}(\mathbf {m}) \over \partial \mathbf {m}} \right _{\mathbf {m}=\mathbf {m}_k}\), is the appropriate Jacobian matrix, \(\underline{\underline{\mathbf {R}}}\) is a regularisation matrix discussed in detail below, and \(\underline{\underline{\mathbf {W}}}_{k}\) is a (Huber) weighting matrix.
\(\underline{\underline{\mathbf {W}}}_{k}\) is updated at each iteration and consists of diagonal elements
\(^kd_i\) is the scalar residual of the ith data point using model vector \(\mathbf {m}_k\) and
being a (robust) estimate of the standard deviation of the residuals at iteration k. We set \(c = 2\), slightly higher than the value of 1.5 usually chosen, in order to ensure that the less numerous polar data are not overly downweighted in the determination of the calibration parameters.
It turns out that the full set of 2046 parameters is not needed to obtain good results and low data misfit, which is confirmed by inspection of the eigenvalues of the matrix \((\underline{\underline{\mathbf {G}}}_k^T \underline{\underline{\mathbf {W}}}_{k} \underline{\underline{\mathbf {G}}}_k + \lambda \underline{\underline{\mathbf {R}}})\), as presented in Fig. 4 for Swarm Alpha. The magnitudes of the sorted eigenvalues (in order of decreasing magnitude) exhibit a distinct drop around 750–800 degrees of freedom, indicating the smaller eigenvalues contribute little to the solution. The inversion of this matrix was therefore finally performed using a truncated singular value decomposition (TSVD) procedure, retaining only 750 degrees of freedom.
A regularisation matrix \(\underline{\underline{\mathbf {R}}}\) is also included to help stabilise the inversion. This is necessary because the Swarm satellites operate in a tightly controlled attitude orientation which leads to a poor excitation of the VFM instrument along the axis perpendicular to the orbit plane (the east–west direction corresponding to the yaxis of the VFM sensor). Consequently, the parameters related to the yaxis are poorly determined in a scalar calibration. The regularisation matrix \(\underline{\underline{\mathbf {R}}}\) is therefore defined so that it acts on the parameters \(s_{2,\mathrm {Tsensor}}\), \(s_{2,\beta }\), \(u_1\), and \(u_3\) to force \(s_{2,\mathrm {Tsensor}} \simeq \left( s_{1,\mathrm {Tsensor}} + s_{3,\mathrm {Tsensor}}\right) /2\) (to reflect the physical properties of the VFM sensor) and also to minimise the norms \(s_{2,\beta }^2\) and \(u_1^2 + u_3^2\). \(\lambda\) is chosen to be sufficiently large to effectively impose the regularisation on the estimated model. Note that no regularisation is directly imposed on \(\delta \vec {B}_\mathrm {Sun}\) but use of truncated SVD during the inversion automatically acts to suppresses structure in regions that are not well constrained by the input data.
The starting model for the inversions is “unity”, i.e. \(\underline{\underline{P}} = \underline{\underline{S}} =\underline{\underline{I}}\), where \(\underline{\underline{I}}\) is the identity matrix, and \(\vec {u}_n^m = \vec {v}_n^m = \vec {0}\). The inversions typically converge within 25 iterations.
Results of model estimation for Swarm Alpha
The model described above is estimated for Swarm Alpha using data from the beginning of the mission (22 November 2013) until June 2015. Figure 1 shows the final scalar residuals, i.e. the residuals after application of the model (after “calibration and correction”) of the VFM measurements (in green), as a function of time together with the residuals of the uncorrected but rescaled vector field measurements, i.e. \(\vec {B}_\mathrm {VFM} + \delta \vec {B}_\mathrm {Sun}\), in light blue; these data illustrate what can be achieved with the traditional scalar calibration methods. Note the excellent reduction of the scalar residuals achieved by the model; the Huber weighted rms of the residuals drops from 963 to 168 pT. Table 2 provides the corresponding numbers for Bravo and Charlie.
Figure 5 shows normal distribution plots for the scalar residuals. The top plot shows the distributions of all data for uncorrected (red) and fully corrected data (green) and demonstrates a transition from a nonGaussian to Gaussian residual distribution when applying the model. The bottom plots show the distributions of the data split into 3month periods, uncorrected to the left and corrected to the right. These also demonstrate the elimination of systematic and nonGaussian effects.
Table 3 lists the estimated \(s_{\mathrm{Tsensor}}\) and \(s_\beta\) parameters and the nonorthogonality values for all three Swarm satellites together with their estimated preflight values for the VFM instrument itself for reference. I.e. the table shows the adjustments applied in order to reduce the scalar residuals to the level indicated above.
Table 4 shows the increase in the weighted rms of the scalar residuals when omitting individual parts of the model—a full reestimation of the remaining model parameters is carried out for each table entry. Particularly the omission of the nonorthogonalities drastically increases the misfit—the power (the meansquare) is more than doubled. Due to the stable attitude of the Swarm satellites, the small x–z nonorthogonality angle, \(u_2\), is equivalent to first order to a small, relative timeshift between the ASM and VFM measurements—1 arcsecond corresponds roughly to a 3ms timeshift, and it has been discussed whether it would be more reasonable to introduce such timeshifts rather than adjusting the preflight estimated nonorthogonalities. However, the variations in the \(u_2\) angles estimated by this model would imply timeshifts varying from \(3\) ms to \(+13\) ms for the individual satellites which, to the authors, seems quite unlikely.
The temporal evolution of the scaling of the vector field measurements, \(s^{{\rm B}\text{}{\rm spline}}\), is shown in Fig. 6 for the various test models listed in Table 4. The full model, shown in red, shows a smooth behaviour in time, as expected from an instrument design perspective. The blue curve shows the model without \(s_\beta\); this exhibits some small oscillations, whereas the light brown (no \(s_{\mathrm{Tsensor}}\)) and green (no \(\delta \vec {B}_\mathrm {Sun}\)) curves show much higher level of oscillations indicating they are inadequate to capture the behaviour of the measurements. The elimination of the oscillations in the full model is a good indicator of the validity of this model. The magenta curve shows the model without nonorthogonalities; this is rather close to the curve of the full model and indicates the decoupling of the nonorthogonalities from any longterm temporal effect of the measurement disturbances and instruments.
Maps of the three components of the estimated disturbance fields from the full model as function of Sun incidence angles \(\alpha\) (abscissa) and \(\beta\) (ordinate) are given in Figs. 7, 8, and 9 for Swarm Alpha, Bravo, and Charlie, respectively. During nominal flight, the Sun incidence angles traverse these plots horizontally from right to left and move up or down in \(\beta\) as the orbit plane moves through local time. The Suninduced disturbance is observed to have temporal characteristics that are observed in the plots as horizontally stretched features, and these are attributed to thermal capacitance: The Suninduced disturbance exhibits characteristic warmup and cooldown effects, i.e. the disturbance increases when the spacecraft is exposed to the Sun, and decreases when the Sun exposure terminates. The time constants for these effects are up to tens of minutes (corresponding to several tens of degrees in the \(\alpha\) angle). This effect is captured by the spherical harmonic model expansion of \(\delta \vec {B}_\text {Sun}\) and yields the horizontally stretched features in Figs. 7, 8, and 9. Note also the regions of nightside data (eclipse), the circled areas to the left of the figures, which generally show less disturbance; this is not imposed by the model or any regularisation; rather, it is simply a result of the data itself and thus another indicator of the ability of the model to describe the observed disturbances. The plots also show both the similarities and the differences in \(\delta \vec {B}_\mathrm {Sun}\) between the three satellites.
Conclusions
We have established a predominantly empirical model for the calibration and correction of the magnetic vector field measurements of the three Swarm spacecraft. The model is based on detailed studies of the observed scalar residuals between the measurements of the absolute scalar magnetometer, ASM, and the modulus of the measurements of the vector field magnetometer, VFM. The model has proven to be quite robust as more data are incorporated into the estimation of the model parameters, although the ambiguity of determining vector disturbances from a pure scalar calibration affects the estimated correction vectors; these corrections do change slightly (by a few tenths of a nT) as more data are added.
The estimated models reduce the scalar differences between the Swarm magnetometers to generally below 0.5 nT with rms values well below 200 pT for all three satellites and have been in operational use since April 2015 to produce corrected Swarm Level 1b magnetic field vector data (as of version 0401).
Future evolutions of the model presented here are foreseen to include changing the model of the temporal evolution of the VFM sensitivity from Bsplines to an exponentially decaying function. Analysis of \(\delta \vec {B}_\mathrm {Sun}\) also indicates that this vector is generally confined to a few, distinct directions which may be incorporated in future models. Finally, it may be possible to model the effect of the thermal capacitance using appropriate temporal filter functions which would lead to a significant reduction of the number of parameters of the model.
Data availability
The estimated disturbance vectors, \(\delta \vec {B}_\mathrm {Sun}\), are included in the operational Level 1b magnetic Swarm data products as dB_Sun.
Uncorrected data are available at ftp://swarmdiss.eo.esa.int/Advanced/ (login required, access can be requested via https://earth.esa.int/Swarm).
References
FriisChristensen E, Lühr H, Hulot G (2006) Swarm: a constellation to study the Earth’s magnetic field. Earth Planets Space 58:351–358
Lesur V, Rother M, Wardinski I, Schachtschneider R, Hamoudi M, Chambodut A (2015) Parent magnetic field models for the IGRF12 GFZcandidates. Earth Planets Space. doi:10.1186/s4062301502396
Olsen N et al (2003) Calibration of the Ørsted vector magnetometer. Earth Planets Space 55:11–18
TøffnerClausen L (2015), Swarm Level 1b processor algorithms. Esa doc. swrsdscsy0002, National Space Institute, DTU Space, Copenhagen
Yin F, Lühr H (2011) Recalibration of the CHAMP satellite magnetic field measurements. Meas Sci Technol 22(5):055,101
Authors’ contributions
LTC carried out the inflight scalar calibration and characterisation, analysed the results, and led the writing of this manuscript. VL proposed the model for the Suninduced vector disturbance, \(\delta \vec {B}_\mathrm {Sun}\), and made the first estimations using this model. NiO and CF supported the entire project with many discussions, suggestions, and source code. All authors read and approved the final manuscript.
Author information
Affiliations
Corresponding author
Ethics declarations
Acknowledgements
We would like to thank ESA for establishing and providing support to the ASMVFM Task Force with the aim of investigating the source of the scalar residuals observed in the Swarm magnetic measurements and developing a correction scheme. We would also like to thank this Task Force for its work in characterising the behaviour of the magnetic disturbance and for many fruitful discussions and inputs for this work. In particular, we would like to thank Peter Brauer from the VFM instrument team for detailed discussions on the modelling and on the characteristics of the VFM instruments. Two anonymous reviewers are thanked for their comments that helped to improve the clarity of the manuscript. This paper is the IPGP contribution 3761.
Competing interests
The authors declare that they have no competing interests.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
TøffnerClausen, L., Lesur, V., Olsen, N. et al. Inflight scalar calibration and characterisation of the Swarm magnetometry package. Earth Planet Sp 68, 129 (2016). https://doi.org/10.1186/s4062301605016
Received:
Accepted:
Published:
Keywords
 Geomagnetism
 Magnetometer
 Instrument calibration
 Satellite
 Swarm