 Full paper
 Open Access
 Published:
Applications for CryoSat2 satellite magnetic data in studies of Earth’s core field variations
Earth, Planets and Space volume 73, Article number: 73 (2021)
Abstract
We use 20 years of continuous magnetic field measurements from the Ørsted, CHAMP and Swarm satellite missions, supplemented by calibrated platform magnetometer data from the CryoSat2 satellite, to study time variations of the Earth’s core field at satellite altitude and at the core–mantle boundary (CMB). From the satellite data we derive composite time series of the core field secular variation (SV) with 4month cadence, at 300 globally distributed Geomagnetic Virtual Observatories (GVO). A previous gap in the GVO series between 2010 and 2014 is successfully filled using CryoSat2, and subdecadal variations are identified during this period. Tests showed that similar subdecadal SV patterns were obtained from the CryoSat2 data regardless of whether IGRF13 or CHAOS6x9 was used in their calibration. Cryosat2 radial field SV series at nonpolar latitudes have a mean standard deviation level compared to smoothing spline fits of 3.5 nT/yr compared to 1.8 nT/yr for CHAMP and 0.9 nT/yr for Swarm. GVO radial SV series display regional fluctuations with 5–10 years duration and amplitudes reaching 20 nT/yr, most notably at low latitudes over Indonesia (2014), over South America and the South Atlantic (2007, 2011 and 2014), and over the central Pacific (2017). Applying the Subtractive Optimally Localized Averages (SOLA) method, we also map the radial SV at the CMB as a collection of locally averaged SV estimates. We demonstrate that using 2year windows of CryoSat2 data, it is possible to reliably estimate the SV and its time derivative, the secular acceleration (SA), at the CMB, with a spatial resolution, corresponding to spherical harmonic degree 10. Along the CMB geographic equator, we find strong SA features with amplitude \(\pm 2.5\mu \mathrm{T}/\mathrm{yr}^2\) under Indonesia from 2011–2014, under central America from 2015 to 2019, and sequences of SA with alternating sign under the Atlantic during 2004–2019. We find that platform magnetometer data from CryoSat2 make a valuable contribution to the emerging picture of subdecadal core field variations. Using 1year windows of data from the Swarm satellites, we show that it is possible to study SA changes at low latitudes on timescales down to 1 year, with spatial resolution corresponding to spherical harmonic degree 10. We find strong positive and negative SA features appearing sidebyside in the Pacific in 2017, and thereafter drift westward.
Introduction
The main part of the Earth’s magnetic field is generated by motions in the electrical conducting liquid outer core, in a process known as the geodynamo. This magnetic field, termed the core field, exhibits both spatial and temporal changes over a broad range of scales. Magnetic measurements from satellites have increased the recovery of smallscale features of this field and revealed rapid changes in its temporal behavior (e.g., Alken et al. 2020a; Finlay et al. 2020; Baerenzung et al. 2020; Ropp et al. 2020; Sabaka et al. 2020). From ground and space magnetic observations, variations in the first and second time derivatives of the field, termed the secular variation (SV) and acceleration (SA), respectively, may now be resolved down to periods of about 1 to 2 years (Lesur et al. 2010, 2017; Ropp et al. 2020).
Studies have revealed oscillating SA pulselike field features at the core–mantle boundary (CMB) focused in the region around the geographical equator (Chulliat et al. 2010; Chulliat and Maus 2014; Chulliat et al. 2015; Sabaka et al. 2018; Alken et al. 2020a). The interpretation and geophysical mechanisms responsible for driving such distinctive behavior in the SA signal remains under debate (e.g., Gillet 2019; Buffett and Matsui 2019; Aubert and Finlay 2019; Gerick et al. 2020), as is the connection to abrupt changes in the SV observed at ground observatories (Mandea et al. 2010). The secular acceleration must be characterized with care, paying attention to those spatial and temporal scales that are well resolved, as its observed spatial spectra at the core surface is blue, showing increasing power with spherical harmonic degree, and its observed temporal spectra seems to be rather flat, meaning that there could be important unresolved fast variations (Christensen et al. 2012; Bouligand et al. 2016; Lesur et al. 2017; Gillet 2019). In this respect, assessing the limitations of the information obtained from measurements by analyzing their resolving power is crucial when aiming to investigate the SA signal at small length scales and short timescales.
Magnetic field measurements from lowEarth orbiting (LEO) satellites provide global field observations, which have proved important for mapping the spatial structures of the core field signal (e.g. Olsen and Stolle 2012). Beginning with the launch of the Danish Ørsted (1999–2014) satellite, the German CHAMP (2000–2010) satellite and the European Swarm (2013) satellite trio, satellites have provided high quality magnetic field measurements, enabling global investigations into the spatiotemporal variations of field. Unfortunately, the CHAMP mission ended in September 2010 and since reliable vector measurements from the Ørsted satellite extend only up to 2006, there is a gap from 2010 and 2014 in the satellite magnetic records (Finlay et al. 2016). However, other satellite missions, not dedicated to measuring the magnetic field, offer a possibility to fill in this gap adding information about the field. In particular, the CryoSat2 (2010) mission, intended for measuring polar ice thickness, carries three platform magnetometers for navigational purposes. Calibrated CryoSat2 measurements, where vector fluxgate magnetometer readings have been transformed into reliable magnetic field vector outputs, from August 2010 to December 2018 have recently become available (Olsen et al. 2020), such that there are now in total 20 years of continuous satellite measurements. Regarding the data from CryoSat2, it is clearly important to assess the quality and limitations of the calibrated platform magnetometer measurements and to test what contribution they can make to the study of core field variations.
The standard approach to using satellite magnetometer data for core field studies is to construct spherical harmonic (SH) field models by leastsquares inversion methods (see e.g., Langel 1987). In such global models Bsplines are often used to parameterize the model timedependence. Typically this necessitates temporal regularization which modifies the timedependence of the harmonics in a nonuniform manner. An undesirable consequence is that, for the higher harmonics, the first time derivative effectively becomes a time average over an increasingly long interval, rather than an estimate of the instantaneous secular variation (Olsen et al. 2009). CryoSat2 data has already been used in the construction of such timedependent spherical harmonic field models by Alken et al. (2020a), Finlay et al. (2020) and Kloss et al. (2021). From the viewpoint of accessing more detailed information on the spatial and temporal structure of the field, it is also of interest to look at alternative techniques for studying secular variation that can complement the traditional SH approach. In this paper, we focus on two local methods for studying core field variations as recorded in satellite measurements, with a focus on assessing the quality and resolving ability of CryoSat2 magnetic data.
In a first assessment of the quality of the CryoSat2 data and its ability to map the SV field at satellite altitude together with CHAMP and Swarm we use the Geomagnetic Virtual Observatory (GVO) technique proposed by Mandea and Olsen (2006) and Olsen and Mandea (2007). This technique involves computing time series of field estimates at specified target locations at satellite altitude, from satellite measurements taken nearby. We apply the processing algorithm recently developed to derive 4monthly Swarm GVO data series (Hammer et al. 2021) and derive time series on a global grid of 300 GVOs. This network of GVOs allows field changes at satellite altitude to be investigated globally at fixed locations. The GVOs provide a useful compression of satellite magnetic measurements and are a convenient dataset for workers wishing to use constraints from satellite measurements for studies of core dynamics. The GVO dataset involves a series of independent local constraints that can be separately assessed rather than the inherently global constraints provided by SH models. In addition the GVOs have well understood error covariances that can be assigned by methods similar to those used with ground observatories. GVOs have already been used by a number of groups for studies of core dynamics (e.g., Whaler and Beggan 2015; Barrois et al. 2018; Domingos et al. 2019).
In the second part of this study, we directly map the radial field SV at the core–mantle boundary, using the technique of Subtractive Optimally Localized Averages (SOLA) that was adapted to geomagnetism by Hammer and Finlay (2019). The SOLA technique can be used to compute estimates of the radial field SV directly at the CMB, based on local spatial averages the SV field centered on target locations of interest and timeaverages over chosen time windows. By collecting many individual SOLA estimates on a grid at the CMB, the SV field can be mapped on regional or global scales. An important foundation of the SOLA technique is that, for noisefree data and a linear forward problem, any linear combination of the data provides a specific average of the true model. With noisy data, a variance is ascribed to this spatial average value such that a tradeoff between resolution and variance arises (Oldenburg 1984; Parker 1994). The SOLA technique readily provides information on the resolution offered by a given set of magnetic field observations, in the form of averaging kernels, as well as estimates of the variance of the locally averaged field. We compare SOLAbased maps of SV and SA estimates derived from CryoSat2 and Swarm data, in order to asses the suitability of the CryoSat2 data for mapping these fields at the CMB. Demonstrating the usefulness of the CryoSat2 data, we then take advantage of this data to map the time evolution of the SA field along the geographic equator at the CMB from 2001 to 2019.
We wish to emphasize that both the GVO and the SOLA methods can result in patterns of secular variation different from those seen in the CHAOS field model, despite the fact they use similar data selection schemes and processing steps that involve the same magnetospheric field model. In the GVO and SOLA methods data close to a location of interest are effectively used to determine localized field or SV estimates. In contrast, estimation of the coefficients of truncated spherical harmonic expansions, in models such as CHAOS, is an inherently global procedure that aims at finding the best possible global model. In the GVO method only data from within a 700km radius of a target location is used to determine the local potential. In the SOLA method measurements far from the target location have essentially no influence on the estimated CMB SV because the data kernels, based on Green’s functions for Laplace’s equation under Neumann boundary conditions, have decreasing sensitivity far away from the target location.
Moreover, the CHAOS model involves temporal regularization whereby one minimizes global norms based on time derivatives of the CMB radial field, integrated over the entire timespan of the model. In the GVO method there is no temporal smoothing beyond the choice of 4month data windows and use of annual differences to produce SV series. In the SOLA method we use 1 or 2year windows to estimate the SV and then annual differences to estimate SA. The SOLA method involves a tradeoff parameter specifying the balance between the spatial resolution and the variance of each local estimate, each estimate being timeaveraged over a 1 or 2year time window; there is no global regularization over longer time spans. Both methods presented here therefore constitute a localized compression of information contained within the satellite data on the potential field near the location of interest within the specified time window. The GVOs and SOLA can thus give a different picture of the core field evolution compared to the sphericalharmonic based field model that involve global temporal regularization, especially regarding rapid changes at short wavelengths when the data quality is high (see Sect. 4.2).
Section 2 describes the satellite measurements used in this study, including the CryoSat2 platform magnetometer data, including how these have been selected and processed. Section 3 presents the GVO technique and results concerning global GVO time series. Section 4 presents the SOLA technique and results of applying this to estimate SV and SA at the CMB. Conclusions and perspectives are given in Sect. 5.
Data selection
We use satellite vector magnetic field measurements from the Ørsted satellite between July 2000 and December 2005, from CHAMP taking L3 magnetic data between July 2000 and September 2010, and from the Swarm trio taking Level 1b MAGL data, version 0505/0506, between January 2014 and April 2020. Most importantly for this study we make use of platform magnetometer data from the CryoSat2 mission, taking calibrated vector measurements with a sampling rate of 4s from the FGM1 magnetometer dataset, version 3, between August 2010 to December 2018. This dataset has had extensive corrections applied for disturbances fields, and was calibrated using a reference field model—for full details see Olsen et al. (2020).
From the Ørsted, CHAMP and Swarm measurements we produced two data sets: dataset #1 used in the GVO application taking a 15 s subsampling of the vector field measurements from Ørsted, CHAMP and Swarm, while taking every 4th measurement from the CryoSat2 dataset (i.e., a 16s subsampling); dataset #2 used in the SOLA application taking a 5s subsampling of the vector field measurements from Ørsted, CHAMP and Swarm, while taking every element of the CryoSat2 dataset with its 4 s sampling rate. Field measurements having gross data outliers for which the vector field components deviated more than 500 nT from the CHAOS7.2 internal field model predictions (Finlay et al. 2020) were rejected. For both data sets we apply a dark geomagnetically quiettime selection criteria scheme. See Table 1 for full details of the selection requirements for the datasets used in the GVO and SOLA applications. In both cases we required the sun to be at least \(10^{\circ }\) below the horizon, adding restrictions on the geomagnetic activity index \(K_p\) and the change in the ring current index (see Olsen et al. 2014), as well as constraints on the merging electric field at magnetopause, and on the magnitudes of the \(B_Y\) and \(B_Z\) components of the interplanetary magnetic field (e.g., Finlay et al. 2020; Ritter et al. 2004). We used minute values of the IMF components and solar wind speed from the OMNI database, http://omniweb.gsfc.nasa.gov, computing twohourly means prior to the time of the considered datum (Finlay et al. 2016).
Since our focus is the core field, we have applied corrections to the data for the lithospheric and external fields. For the lithospheric field model, we used the LCS1 model (Olsen et al. 2017); the precise choice of lithospheric field is not crucial for studies of the SV and SA. For the solarquiet ionospheric field and associated induced fields we used the CIY4 model (Sabaka et al. 2018). For the magnetospheric field and related induced fields we used the CHAOS7 model (Finlay et al. 2020). These models were chosen as they are well established and compatible with the data selection criteria described above.
As noted in Table 1, for the GVO application we use the sums and differences of the magnetic field measurements. It has been shown by Olsen (2015) and Sabaka et al. (2018), that taking differences of the satellite measurements alongtrack and East–West (between Swarm satellites A and C) helps the recovery of the smallscale core field, as this reduces the impact of correlated errors caused by unmodeled largescale external fields. Here we follow such an approach by taking differences of the measurements, but we also include alongtrack and East–West sums of the measurements in order to ensure sufficient constraint on the larger wavelengths of the field (Sabaka et al. 2013; Hammer 2018). We denote the magnetic vector measurements by \(B_k(\mathbf {r})\), where k is any of the three given vector component of the field, \(\Delta d_k\) and \(\Sigma d_k\) denote measurement differences and sums of this particular component, respectively. Here the alongtrack (AT) and East–West (EW) data differences are denoted by \(\Delta d_k=(\Delta d_k^{\mathrm{AT}},\Delta d_k^{\mathrm{EW}})\), and the data sums by \(\Sigma d_k=(\Sigma d_k^{\mathrm{AT}},\Sigma d_k^{\mathrm{EW}})\). The alongtrack data differences are calculated using the 15 s differences \(\Delta d_k^{\mathrm{AT}} = [B_k(\mathbf {r},t)  B_k(\mathbf {r}+\delta \mathbf {r},t+15s)]\). The alongtrack sums were calculated as \(\Sigma d_k^{AT} = [B_k(\mathbf {r},t) + B_k(\mathbf {r}+\delta \mathbf {r},t+15s)]/2\). For Swarm, East–West differences were calculated as \(\Delta d_k^{\mathrm{EW}} = [B_k^{\mathrm{SWA}}(\mathbf {r}_1,t_1)  B_k^{\mathrm{SWC}}(\mathbf {r}_2,t_2)]\) having an East–West orbit separation between the Swarm Alpha (SWA) and Charlie (SWC) satellites of \(\approx 1.4^{\circ }\) corresponding to 155 km at the equator (Olsen 2015). The East–West sums were calculated as \(\Sigma d_k^{\mathrm{EW}} = [B_k^{\mathrm{SWA}}(\mathbf {r}_1,t_1) + B_k^{\mathrm{SWC}}(\mathbf {r}_2,t_2)]/2\). For a particular orbit of Swarm Alpha the corresponding Swarm Charlie data were selected to be those closest in colatitude with the condition that \(\vert \Delta t\vert =\vert t_1t_2\vert <50s\) (Olsen 2015).
Application I: geomagnetic virtual observatories
4monthly core field GVOs and secular variation estimates
In the first application, we compute Geomagnetic Virtual Observatory time series derived from dataset #1. Of particular interest is the quality of the GVO series obtained from CryoSat2 data compared with similar series obtained from CHAMP and Swarm data. The time series consist of estimates of the geocentric spherical polar vector components of the magnetic field at specified target points, referred to as GVOs (Mandea and Olsen 2006; Olsen and Mandea 2007). Here we use the same algorithm described in detail by Hammer et al. (2021) (see also http://www.spacecenter.dk/files/magneticmodels/GVO/GVO_Product_Algorithm.pdf) to produce the Swarm GVO product, and derive global grids of 300 uniformly distributed GVO time series each having 4month cadence. The GVOs are located in an approximately equal area grid based on the sphere computed using the algorithm of Leopardi (2006).
For each GVO in the grid we take data from within a cylinder of horizontal radius \(r_{cyl}=700\,\hbox {km}\). The GVOs have the spherical polar coordinates \(\mathbf {r}_{GVO}=(r,\theta ,\phi )\), and are placed at fixed altitudes \(r=r_a+h_{GVO}\) where \(h_{GVO}\) is the height above the Earth’s mean spherical radius \(r_a=6371.2\,\hbox {km}\). For the CHAMP, CryoSat2 and Swarm missions the GVO altitudes were chosen as \(h_{GVO}=370\,\hbox {km}\), \(h_{GVO}=727\,\hbox {km}\) and \(h_{GVO}=490\,\hbox {km}\), respectively, such that the GVOs are located at approximately the mean orbital altitude for each mission during the time interval considered.
The input measurements of dataset #1 are provided in an EarthCenteredEarthFixed (ECEF) coordinate frame by the spherical polar components \(\mathbf {B}^{obs}=(B_r,B_{\theta },B_{\phi })\). From the extracted vector field measurements surrounding each GVO, within a radius of 700 km and within a 4month time window, a residual magnetic field, \(\delta \mathbf {B}\), is first computed by subtracting off estimates of the main field and noncore fields:
where the field estimates removed are: a) \(\mathbf {B}^{MF}\) the internal field for SH degrees \(n \in [1,13]\) as given by IGRF13 (Alken et al. 2020b), b) \(\mathbf {B}^{lith}\) the static internal field for SH degrees \(n \in [14,185]\) as given by the LCS1 model (Olsen et al. 2017), c) \(\mathbf {B}^{mag}\) the magnetospheric and associated induced field as given by the CHAOS7.2, model (Finlay et al. 2020), d) \(\mathbf {B}^{iono}\) the ionospheric and associated induced field as given by the CIY4 model (Sabaka et al. 2018). Note that estimates of the main field from IGRF13 (with linear time dependence over 5year intervals) are subtracted here. At a later stage, main field estimates, again from IGRF13, are added back for the specified GVO times and positions. Removal of a main field at this stage allows for a more effective prewhitening of the data, such that Huber weights, used in the robust GVO estimation scheme, can be well determined. We emphasize that this approach still allows us to capture departures from the subtracted main field when required by the satellite data.
Of the remaining residual field, we are interested in the core field part of that signal. Although we have removed estimates of the external fields and their associated Earthinduced counterpart, contributions from noncore sources remain in the residual field. A particular concern is contamination from fields caused by rapidly varying ionospheric currents, for example polar electrojet currents at high latitudes in the Elayer, and possibly also signatures of Flayer currents at mid and low latitudes. Further work is needed on these aspects. We attempt to mitigate leakage of fieldaligned currents by removing toroidal field estimates (Sabaka et al. 2010) obtained by performing an epochbyepoch spherical harmonic analysis performed on the global network of GVOs (Hammer et al. 2021), see below for further details.
Next, the residual magnetic field data and their positions are transformed from the spherical system to a righthanded local topocentric Cartesian system (x, y, z) with origin at the GVO target location. At the GVO location and only at this location, x points towards geographic south, y points towards east and z points radially upwards (Hammer et al. 2021). Assuming that the satellite measurements are made in a source free region, the residual magnetic field, \(\delta \mathbf {B}\), is a Laplacian potential field. In the local Cartesian coordinate system the magnetic scalar potential, V, can be expanded as a sum of polynomials of the form \(C_{abc} x^a y^b z^c\) (Backus et al. 1996). Here we use this expansion out to cubic terms:
The forward problem linking the vector of GVO model coefficients, \(\mathbf {m}={[C_{100},C_{010},..., C_{111}]}^T\), with the data vector \(\mathbf {d}^{vec}\) containing the residual field components, \(\delta \mathbf {B}\), can be written
where \(\underline{\underline{\mathbf {G}}}^{vec}\) is a design matrix derived from the spatial derivatives of Eq. (2). As noted above, instead of using residual vector field components to compute the potential, we use sums and differences of the residual field vector components such that the data vector is \(\mathbf {d}=[\Delta d_{x}^{vec},\Delta d_{y}^{vec},\Delta d_{z}^{vec},\Sigma d_{x}^{vec},\Sigma d_{y}^{vec},\Sigma d_{z}^{vec}]^T\), where \(\Delta\) and \(\Sigma\) denotes the differences and sums of the computed residual field as described in Sect. 2. The relevant design matrix is then constructed as \(\underline{\underline{\mathbf {G}}}=[\Delta G_x^{vec};\Delta G_y^{vec};\Delta G_z^{vec};\Sigma G_x^{vec};\Sigma G_y^{vec};\Sigma G_z^{vec})]\) where \(\Delta G_{k}^{vec} =[G_{k}^{vec}(\mathbf {r}_1)  G_{k}^{vec}(\mathbf {r}_2)]\) and \(\Sigma G_{k}^{vec} = [G_{k}^{vec}(\mathbf {r}_1) +G_{k}^{vec}(\mathbf {r}_2)]/2\), where \(k=(x,y,z)\).
To determine the GVO model coefficients, we use a robust iterativelyreweighted leastsquares inversion scheme, based on a diagonal weight matrix consisting of robust (Huber) weights for each entry in the data vector (e.g., Constable 1988). We also include a downweighting factor of 1/2 for the Swarm satellites Alpha and Charlie accounting for the fact that these two satellites fly sidebyside and therefore do not provide completely independent measurements. A minimum number of 30 data points were required to computing the inversion. Using the resulting coefficients of the potential for a given GVO target location and time, derived from the associated sums and differences satellite data, a prediction for the mean residual field at the GVO target point and epoch can be computed as \(\delta \mathbf {B}_{GVO}(x,y,z)=\nabla V(0,0,0)=(C_{100},C_{010},C_{001})\).
Moving back to the vector components in spherical polar coordinates, \(\delta B_{GVO,r}=\delta B_{GVO,z}\), \(\delta B_{GVO,\theta }=\delta B_{GVO,x}\), \(\delta B_{GVO,\phi }=\delta B_{GVO,y}\). We then add back the IGRF13 main field predictions for the given target point and epoch, \(\mathbf {B}_{GVO}^{MF}(\mathbf {r}_{GVO},t)\) to obtain
The above procedure is then repeated for each epoch and each component to obtain time series of GVO estimates of the vector magnetic field at the GVO target locations.
The GVO method assumes that the residual field in Eq.(1), is a potential field. However, because the satellite measurements are made in the ionospheric Fregion, in situ currents can cause nonpotential fields to leak into the estimated potential (Olsen and Mandea 2007). Therefore, in a final postprocessing step, we carry out an epochby epoch spherical harmonic analysis of our 4monthly GVOs, estimating external and toroidal field contributions to SH degree 13 or to degree 6 if fewer than 300 GVOs are available. For epochs having an insufficient number of GVOs available to ensure a stable solution, the external and toroidal coefficients were obtained via a linear interpolation between nearby epochs. The external and toroidal field estimates, reaching a level of \(\pm 15\)nT at high latitudes, are then removed to obtain the estimated Core Field GVO time series (Hammer et al. 2021).
Secular variation at a given GVO location is computed using annual differences between values at time \(t+6\) months and at time \(t6\) months. We have chosen to take annual differences, as this helps to avoid annual noncore signals that may persists in the GVO series despite our best efforts in reducing such contamination.
Figure 1 presents the number of 4monthly Core Field GVO estimates during the past 20 years. The maximum possible number of GVOs per epoch is 300. A strong dip in the number of GVOs is seen during 20022004 due to increased solar activity that meant there were fewer data meeting our selection criteria. As noted above, if fewer than 30 measurements are available within a GVO target cylinder during a given 4month window, we were unable to reliably determine GVO estimates. The remaining epochs from 20042020 are well covered with only few epochs having less than 250 GVOs available. We find that using CryoSat2 data, we are able to provide between 200 and 300 GVO estimates at all times during the gap between the end of the CHAMP mission and the start of the Swarm mission.
GVO results: global time series of secular variation from 2002 to 2020
Table 2 presents the rootmeansquare (rms) and mean of the residuals between the contributing satellite data (sums and differences) and GVO model predictions, summing over all GVOs for a given vector component and a given region (polar or nonpolar). Here we defined polar to be polward of \(\pm 54^{\circ }\) geographic latitude. The polar rms values for both sums and differences are higher than the nonpolar, and the CHAMP values are slightly higher than the Swarm values. The CryoSat2 values are seen to be higher for all components but not unreasonable, given they are derived from platform magnetometer data. The nonpolar rms values for all components are below 2 nT during both CHAMP and Swarm times. The CryoSat2 GVO’s rms values are as expected larger, and especially for the alongtrack differences, indicating that alongtrack correlated noise is less dominant, due to the presence of other noise sources in platform magnetometer data.
The CryoSat2 magnetometer data have been cleaned from known platform signals and calibrated as described by Olsen et al. (2020). This calibration relies on computing residuals with respect to a reference magnetic field which was taken from the CHAOS6x9 field model (Finlay et al. 2016). Here, we carried out an experiment to verify that the GVO secular variation signals obtain from CryoSat2 data are independent of the main field model used for data calibration. To do this we computed CryoSat2 GVO estimates using two different datasets, the first being the official CryoSat2 dataset calibrated using the CHAOS6x9 field model and the second a test version of the CryoSat2 data calibrated instead using IGRF13 (Alken et al. 2020b). Having estimated global grids of GVO series for each dataset, to each series we fit cubic smoothing splines, with a knot spacing at every 4 months, and with the smoothing parameter determined using a GCV (generalized crossvalidation) approach (Green and Silverman 1993). Figure 2 presents SV series for the three field components at three example GVOs, with colatitude/longitude \((36^{\circ },77^{\circ })\) (top), \((72^{\circ },110^{\circ })\) (Center) and \((84^{\circ },137^{\circ })\) (bottom). The SV time series derived from the CHAOS6x9 calibrated dataset are shown with red dots (spline fit in the red line) and similar series from the IGRF13 calibrated dataset are shown with blue dots (spline fit in the blue line). SV model predictions from the CHAOS6x9 model up to SH degree 16 (in green), and IGRF13 (in black), are also shown for reference. There are certainly some differences in the GVO SV estimates derived using the CHAOS6x9 and IGRF13 calibrated data. Since the CryoSat2 calibration relies on computing residuals with respect to a reference field model (here either CHAOS or IGRF), we do expect such differences, especially since the IGRF model assumes a crude piecewise constant SV field. Nevertheless, the IGRF calibrated GVO SV estimates are clearly closer to the CHAOS6x9 SV predictions than to the IGRF13 SV predictions, by an rms of 1.2nT/yr for the horizontal components, and 3.9nT/yr for the radial component. This indicates that the CryoSat2 data does possess an SV signal regardless of the model used to calibrate them.
We find that the CHAOS6x9 and IGRF13 calibrated CryoSat2 GVO SV series show similar subdecadal changes, neither of which exactly match those seen in CHAOS6x9. The overall similarity of the IGRF and CHAOS calibrated SV series give us confidence that the subdecadal secular variation seen in CryoSat2 data is independent of the field model used in its calibration. In addition, we note a strong acceleration in the radial SV component from 2014 to 2018 in both sets of GVOs (top plot) located above Northwest Siberia, a change which is only partially captured by the piecewise constant IGRF13 SV model.
In Table 3, we report the mean rms differences between the GVO SV estimates and GCV spline fits, separated by component (\(B_r,B_\theta ,B_\phi\) and also the intensity F) and by polar and nonpolar regions. These numbers give an indication of the scatter in the SV datasets and allow the quality of the GVO SV series obtained from CHAMP, Swarm and CryoSat2 (both CHAOS6x9 and IGRF13 calibrated versions) to be compared. Similar results are seen for the CHAOS6x9 and IGRF13 model calibrated GVO SV series, with rms differences between GCV fits and intensity SV data, taken the over all GVOs, of \(3.5\,\hbox {nT}/\hbox {yr}\) and \(3.4\,\hbox {nT}/\hbox {yr}\), respectively. As expected, similar numbers for the CHAMP and especially the Swarm derived GVOs, are much smaller being \(1.8\,\hbox {nT}/\hbox {yr}\) and \(0.9\,\hbox {nT}/\hbox {yr}\), respectively. This indicates the lower scatter in the CHAMP and Swarm GVO series.
Figure 3 goes beyond detailed comparisons at a few example locations and presents a global map of CryoSat2 GVO time series for the radial field SV, showing both the official CHAOS6x9 calibrated dataset (red dots) and the IGRF13 calibrated test dataset (blue dots). The scale is shown in the bottom left corner, with the yaxis being 20 nT/yr and the xaxis going from 2011 to 2018. Visual inspection of Fig. 3 supports the results of Table 3, and shows similar scatter levels with rms differences of less than \(1.0\hbox {nT}/\hbox {yr}\) in all three components of the SV series derived from the CHAOS6x9 and IGRF13 calibrations.
With the availability of CryoSat2 magnetic data, it is now possible to use GVOs to study secular variation globally over the past 20 years. In order to illustrate the quality of information they can provide, in Fig. 4 we first present comparisons of GVO SV series, from CHAMP, CryoSat2 and Swarm, to annual differences of Revised Monthly Means (rmm) (Olsen et al. 2014) from example high quality ground observatories, Kourou in South America (top plots), from Novosibirsk observatory in Siberia (middle plots) and from Honolulu ground observatory in the central Pacific (bottom plots). Each plot shows the spherical polar components of the annual differences of revised monthly means (black dots) computed from the ground observatory data (Olsen et al. 2014), and GVO SV time series derived from CHAMP (purple dots), CHAOS6x9 calibrated CryoSat2 (blue dots) and Swarm (red dots) data mapped to ground level. Here the GVO time series for each satellite mission have been mapped to ground by subtracting the field difference as given by the CHAOS7.2 model for SH degrees 120, between ground and GVO altitudes which are close to the mean orbital altitude for each mission. At all three locations the GVOs derived from CryoSat2 data have more scatter – this is particularly noticeable in the \(\theta\) and \(\phi\)components for the examples shown, whereas the scatter in the radial component is closer to the level seen in the CHAMP and Swarm derived GVOs. Table 3 indicates that at nonpolar latitudes the scatter in the r and \(\theta\)components is generally similar, 3.4 nT/year compared to 3.8 nT/year. Notice that the scatter in the ground observatory rmm’s is also enhanced in the horizontal components. The good agreement between the independently estimated CryoSat2 and Swarmderived GVOs at overlapping epochs from 2014 to 2018, is particularly evident in the r and \(\phi\)components. SV variations are coherent in both phase and amplitude between the CryoSat2 GVO time series and the ground observatory records, thus confirming that the CryoSat2 GVOs are able to track the same field changes as observed by ground observatory records on timescales of 1 year and longer.
Next, in Fig. 5 we present a global map of the SV of the radial field component over the past 18 years. Here, to ease visualization, GVO SV series from CHAMP (covering 20022010), CHAOS6x9 calibrated CryoSat2 (covering 20102014) and Swarm (covering 20142020) have been mapped to a common altitude of 700 km, again using the CHAOS7.2 field model, and combined into one composite time series. This allows for easy visual investigation of global patterns of subdecadal SV. We find that regions at low latitudes display strong subdecadal variations, however, not simultaneously at all longitudes. For instance, we observe a change of slope in the radial SV field occurring over the south Atlantic region around 2007 and again in 2014, over Indonesia around 2014, and in the Pacific region centered in 2017. Some of these variations are characterized by distinct “\(\Lambda\)” and “V”shaped behavior occurring over time spans of 5–10 years and locally confined to specific regions, but otherwise reminiscence of the often discussed phenomenon of geomagnetic jerks (Mandea et al. 2010). The availability of CryoSat2 magnetic field data plays a key role in permitting a continuous coverage from satellitebased time series without a gap between 2010 and 2013, thus allowing the study of global patterns in the timevarying core field secular variation over the past 20 years.
Application II: SOLA
SOLA method for local estimation of CMB radial field SV
We now move on to investigate the behavior of the core field not at satellite altitude but down at the core–mantle boundary (CMB), on the edge of the region where it originates. To do this we use the Green’s functions of the Neumann boundary value problem that links the magnetic field at satellite altitude to the radial field at the CMB. Following Hammer and Finlay (2019), a localized estimate, \(\widehat{B}_r\), of the radial magnetic field at a target location and time, \((\mathbf {r}_0,t_0)\), at the CMB can be computed as a localized spatial average around the target location timeaveraged over a specified interval. Because the CMB radial magnetic field is linearly related to the spherical polar components of the vector field at satellite altitude, we can write \(\widehat{B}_r\) as a weighted linear combination of the satellite magnetic measurements with weights \(q_n\) (Backus and Gilbert 1968, 1970; Hammer and Finlay 2019)
where \(d_n\) are satellite magnetic measurements \((n=1,...,N)\) within a specified time window and \(q_n\) are weighting coefficients to be determined. Here the data \(d_{n}\), at positions \(\mathbf {r}_n\) and times \(t_n\), are taken from dataset #2 and for simplicity we consider using only observations of the radial component of the field. Corrections for the lithospheric field for SH degrees \(n \in [14,185]\) as given by the LCS1 model (Olsen et al. 2017), for the magnetospheric and associated induced fields as given by the CHAOS7.2 model (Finlay et al. 2020), and for the ionospheric and associated induced fields as given by the CIY4 model, (Sabaka et al. 2018) are removed from the observations in a preprocessing step. The radial magnetic field measurements, \(d_{n}(\mathbf {r}_n,t_n)\), are then related to the radial magnetic field \(B_r(\mathbf {r}',t_n)\), integrated over the CMB, by (Gubbins and Roberts 1983):
where the surface element is \(dS'=\mathrm{sin}\theta 'd\theta 'd\phi '\). The data kernel \(G_{r}(\mathbf {r}_n, \mathbf {r}')\) is the radial derivative with respect to \(\mathbf {r}\), of the Green’s functions for the exterior Neumann boundary value problem (e.g., Gubbins and Roberts 1983; Barton 1989):
where \(h_n=r'/r_n\) and \(r'\) is the CMB radius, \(f_n=R_n/r_n\) where \(R_n=\sqrt{r_n^2+r'^22r_nr'\zeta _n}\) and \(\zeta _n = \mathrm{cos}\, \gamma _n=\mathrm{cos}\,\theta _n \, \mathrm{cos}\,\theta '+\mathrm{sin}\,\theta _n \mathrm{sin}\,\theta '\, \mathrm{cos(\phi _n\phi ')}\), where \(\gamma _n\) is the angular distance between a measurement at position \((\theta _n,\phi _n)\) and a position on the CMB \((\theta ',\phi ')\). The data kernel describes how a particular measurement samples the CMB radial field; radial magnetic measurements sample the CMB radial field most strongly directly below the measurement position. Regarding the timedependence, we use a first order Taylor expansion around a reference time \(t_0\), such that
The time difference, \(\Delta t_n = t_nt_0\), is computed with respect to the target time, \(t_0\). Inserting Eq.(8) into Eq.(5), we obtain
where \(\mathcal {K}(\mathbf {r}_0, {t_0}, \mathbf {r}')\) and \(\dot{{\mathcal {K}}}(\mathbf {r}_0, t_0, \mathbf {r}')\) are spatial averaging kernels for the CMB field and secular variation, respectively, constructed from the weighting coefficients and the data kernels
By varying the weight coefficients, \(q_n\), the shape of the averaging kernels change. Notice that time differences \(\Delta t_n\), between the measurement times and the target time, are effectively additional temporal weights applied to the kernel \({\mathcal {K}}\) in order to obtain \(\dot{{\mathcal {K}}}\).
In order to obtain estimates of the secular variation of the radial field on the CMB, at the target location and time \(\widehat{\dot{B}}_r(\mathbf {r}_0,t_0)\), we minimize the following objective function:
where \(\lambda\) is a tradeoff parameter (units of \([\mathrm{nT}^{1}]\)), \(\mathbf {q}\) is vector of the weighting coefficients, \(\underline{\underline{\mathbf {E}}}\) is the data error covariance matrix which we define below and \(\dot{{\mathcal {T}}}\) is an SV target kernel that we choose to be a Fisher distribution on the sphere (Fisher 1953):
where \(\gamma _0\) is the angular distance on the CMB between the target position \((\theta _0,\phi _0)\) and another position \((\theta ',\phi ')\). On the basis of tests carried out by Hammer and Finlay (2019), we initially set \(\kappa =600\) corresponding to a target kernel width of \(15^{\circ }\); this is narrower than can be achieved for \(\dot{{\mathcal {K}}}\) with the available data, but it avoids excessive ringing associated with taking a Dirac delta function as the target kernel. When computing SOLA estimates for a given time window, we select a subset (\(n=1,...,N\)) of the measurements. Using this data subset the data error covariance matrix \(\mathbf {E}\) is defined as follows. Using all available measurements (\(m=1,...,M\)), for each satellite mission within 2 degree bins of quasidipole (QD) latitude (Richmond 1995), we first derived robust data error variances as a function of QD latitude:
where \(\epsilon _m\) are residuals with respect to predictions of the CHAOS7.2 internal field model for SH degrees \(n \in [1,13]\), \(\mu\) are robust mean residuals within the considered bin and \(w_m\) are Huber weights (e.g., Constable 1988) for the data within each bin. Here we use QD coordinates, as this is appropriate for characterizing processes related to unmodeled ionospheric currents which we consider to be a likely source of contamination, especially at high latitudes. Figure 6 presents the resulting QDlatitudedependent error estimates \(\sigma (\theta _{QD})\) for the radial field component used in this study, comparing the values for the Ørsted, CHAMP, CryoSat2 and Swarm datasets. When computing SOLA estimates for a specified time window of, e.g., 2 years, we select a data subset of dataset #2. Using this data subset of N measurements, a specific data error covariance matrix E is computed. Diagonal elements of this data error covariance matrix \(\mathbf {E}\) are finally defined as:
where \(w_n\) are additional robust (Huber) weights determined a priori for each datum (\(n=1,...,N\)), based on their residual to CHAOS7.2, in order to account for the expected longtailed error distribution. Offdiagonal elements of \(\mathbf {E}\) are set to zero.
In addition to minimizing Eq.(11), we simultaneously impose the following constraint:
where the first term is in practice very small when estimating the SV, since it is minimized in Eq. (11). This constraint ensures that a valid averaging kernel is obtained.
Discretization of integrals over the CMB was carried out using Lebedev quadrature (Lebedev and Laikov 1999) and the system of equations was solved for the coefficients, \(q_n\), using a Lagrange multiplier method, see Hammer and Finlay (2019) for further details.
Once SOLA estimates of the CMB radial field SV at the chosen target location and epoch are obtained, by minimizing Eq. (11) subject Eq. (15), we are able to easily appraise them based on (i) their averaging kernel width, which we define as the angular distance between the points at which the averaging kernel first reaches zero amplitude moving away from its maximum value, and (ii) the variance of the SOLA estimate which we computed as:
By changing the parameter \(\lambda\) of Eq.(11), a range of solutions can be computed which describes a tradeoff between having an averaging kernel width as small as possible and the variance of the estimate being as small as possible (Parker 1977). Below we discuss the effect of changing \(\lambda\) on our results.
Results: SOLA estimates of CMB SV and SA from Ørsted, CHAMP, CryoSat2 and Swarm data
We begin by first comparing SOLA estimates for the CMB radial field SV obtained using separate data subsets from the Swarm and CryoSat2 missions, respectively. First, Swarm and CryoSat2 data subsets are extracted from the main dataset #2 described in Sect. 2, so that each cover the same 2year time window from 2015.0 to 2017.0 Next, in order to obtain data subsets with suitable spatial and temporal coverage, we considered bins surrounding each point in an approximately equaldistance grid at satellite altitude of \(\approx 2.5^{\circ }\) spacing, based on the partitioning algorithm of Leopardi (2006), and randomly sampled one datapoint from each bin, resetting the bins every 2 months. Data subsets spanning the full 2year window from 2015.0 to 2017.0 were produced by accumulating these 2monthly globallydistributed subsets. The resulting Swarm and CryoSat2 data subsets spanning 2015.0 to 2017.0 consisted of 62469 and 54685 radial field observations, respectively.
In Fig. 7, we compare maps collecting SOLA CMB radial field SV estimates centered on epoch 2016.0, derived using the CryoSat2 and Swarm data subsets spanning 2015.0 to 2017.0. To ensure that we obtained SOLA estimates of comparable resolution, we first computed SOLA SV estimates using the Swarm data subset and taking \(\lambda =3 \times 10^{3}\ \mathrm{nT}^{1}\). This resulted in wellbehaved averaging kernels with widths \(\approx 38^{\circ }\). Next, we used these averaging kernels as the target kernels in order to derive similar estimates using the Swarm and CryoSat2 data subsets, thus effectively seeking Swarm and CryoSat2 SV estimates with the same spatial resolution. The maps in the top row of Fig. 7 show the resulting global collections of SOLA SV estimates, obtained on a \(1^{\circ }\) grid of target locations at the CMB, based on the CryoSat2 (left plot) and Swarm (right plot) data subsets, respectively.
Maps collecting the formal standard deviations for each SOLA estimate (derived using Eq.16) and their averaging kernel widths are presented in middle and bottom rows of Fig. 7. The standard deviations of the Swarmbased SOLA SV estimates are fairly homogeneous with values of \(0.30.4\mu \hbox {T}/\hbox {yr}\); those for CryoSat2 SOLA are somewhat larger, being in the range \(\approx 0.61.8\mu \hbox {T}/\hbox {yr}\). We note that the CryoSat2 error estimates are slightly lower at higher latitudes. The same is true for the Swarm map, but the variations in errors estimates are in that case much smaller. Kernel widths in both cases are also fairly homogeneous except at auroral latitudes where distinct behavior of the kernels are found related to the data error estimates having increased amplitude, as seen in Fig. 6.
As seen from the kernel widths in the bottom row of Fig. 7, very similar resolution has been obtained (by construction) in the Swarm and CryoSat2 SV maps. The same field features are clearly identified in both maps. For instance, we notice high latitude SV patches in the northern hemisphere which have been associated with a high latitude jet of core flow (Livermore et al. 2017), and there is increased amplitude of SV over the hemisphere centered on the Atlantic in comparison with the Pacific hemisphere. This first test gives us confidence that the CryoSat2 measurements can indeed be used to reliably map SV features at the CMB on a timescale of 2 years, and with a spatial resolution down to \(\approx 38^{\circ }\) degrees.
Next, we go further and investigate the second time derivative or secular acceleration (SA) of the radial field at the CMB, which is of great interest for investigating the dynamics of the core (e.g., Chulliat et al. 2010; Finlay et al. 2016; ChiDurán et al. 2020). Again we first compare maps based on CryoSat2 and Swarm data. To obtain SA estimates we initially use the accumulated change between SV estimates 2 years apart. In particular, in Fig. 8 we show the SA in 2017 based on the difference between SOLA SV estimates in 2016.0 and 2018.0. Note, that this is not an instantaneous secular acceleration but a centered difference in SV estimates 2 years apart, each based on 2 years of data. To study the SA, we computed SOLA SV estimates from Swarm data taking \(\lambda =1 \times 10^{2}\ \mathrm{nT}^{1}\), and used the resulting associated averaging kernels as the target kernels for the SOLA estimates from both Swarm and CryoSat2 data.
Figure 8 presents global grids of SOLA CMB radial field SA estimates centered on 2017.0, again with a \(1^{\circ }\) spacing, derived from CryoSat2 (left plot) and Swarm (right plot) data subsets. Here the map is centered on the Pacific region where there has been interesting SA activity during the past 6 years (Finlay et al. 2016, 2020). As for the SV maps, we find error estimates, computed assuming the contributing SV estimates have independent errors, to be fairly homogeneous, with values ranging between \(0.110.27\mu \mathrm{T}/\mathrm{yr}^2\) for the estimates derived using CryoSat2 data and \(0.060.08\mu \mathrm{T}/\mathrm{yr}^2\) for the estimates derived using Swarm data. In both cases kernel widths are close to \(\approx 42^{\circ }\), except in the auroral region. Comparing the CryoSat2 and Swarmbased SA maps, similar features can be observed. In particular, this is the case for the features seen under Asia and Indonesia. A distinctive feature reproduced in both maps is the sequence of intense patches of SA at low latitudes in a localized region below central America, having amplitudes of approximately \(1.9\pm 0.3\mu \mathrm{T}/\mathrm{yr}^2\) and \(1.7\pm 0.1\mu \mathrm{T}/\mathrm{yr}^2\) for the CryoSat2 and Swarm maps, respectively. The location and amplitude of these features are similar in the two maps, confirming that CryoSat2 data can be used to track such SA structures. Because the SOLA estimates are local averages, distant high latitude measurements, where ionospheric electrical currents may be prominent even during dark quiet times, will have little influence on such low latitude SV and SA estimates. Finally, we note a strong patch under the Bering Sea of amplitude approximately \(1.4\pm 0.3\mu \mathrm{T}/\mathrm{yr}^2\) and \(1.0\pm 0.1\mu \mathrm{T}/\mathrm{yr}^2\) for the CryoSat2 and Swarm maps, respectively.
These initial investigations of the CMB radial field SV and SA using the SOLA technique indicate, as seen earlier in the GVO time series, that low latitude regions experience significant subdecadal core field variations. We are therefore motivated to study the field timedependence in the equatorial region in more detail. We do this in Fig. 9 by constructing timelongitude (TL) plots of SOLA CMB radial field SA estimates along the geographic equator from 2002 to 2019, centered on the Pacific. We again compute our SA estimates based on differences of SOLA SV estimates 2 years apart, each derived from 2year data windows, and sliding the windows in 2month steps. Here we use radial field data from the Ørsted, CHAMP, CryoSat2 and Swarm satellites.
The top left plot in Fig. 9 presents estimates based on averaging kernels obtained from Swarm data using \(\lambda =3 \times 10^{2}\ \mathrm{nT}^{1}\). These have associated error estimates between \(0.020.08\mu \mathrm{T}/\mathrm{yr}^2\) and kernel widths \(\approx 50^{\circ }\). Their resolution is lower than that shown in Fig. 8 and corresponds to approximately SH degree 8. For comparison the SA predicted by the CHAOS7.2 model for SH degrees \(n \in [1,8]\) are shown on the top right plot. Note here that the CHAOS7.2 model makes use of uncalibrated vector magnetic data from CryoSat2 between 2010 and 2014, and coestimates magnetometer calibration parameters (Finlay et al. 2020). The SOLA and CHAOS TL plots in the top row of Fig. 9 show largely the same SA features, illustrating the convergence of the two techniques at long wavelengths of the SA and when constructing SOLA SA estimates from SV differences between consecutive 2year time windows. For instance, the evolution of the SA features observed in Fig. 8 under the central Americas, can be identified ranging from longitudes \(240^{\circ }\) to \(320^{\circ }\) centered on 2017. In addition, we find strong SA patches in the CryoSat2 data around 2013 at longitudes \(70^{\circ }\) to \(160^{\circ }\) and \(280^{\circ }\) to \(320^{\circ }\). Notice, that there seems to be a sign changing sequence occurring at longitudes \(240^{\circ }\) to \(320^{\circ }\) going from 2005 to 2019.
Next, we increase the spatial resolution by instead deriving SV estimates using \(\lambda =1 \times 10^{2}\ \mathrm{nT}^{1}\), which leads to slightly larger error estimates in the range \(0.10.5\mu \mathrm{T}/\mathrm{yr}^2\) and kernel widths \(\approx 42^{\circ }\), i.e., similar to Fig. 8. This is shown in the bottom left plot while the CHAOS7.2 model predictions for SH degrees \(n \in [1,10]\), matching approximately the kernel width, are shown on the bottom right plot. Although the SOLA TLplot looks somewhat noisier than the CHAOS plot, similar coherent evolving structures having higher amplitudes can clearly be identified. The noisier appearance in the interval 20102014 likely indicates the limitations of the CryoSat2 data, but they clearly provide useful information during this period.
With data from the Swarm mission, it is possible to go further and also increase the temporal resolution of the SA by taking 1year differences of SOLA SV estimates derived from 1year data windows and sliding in 1month steps. The result of applying this procedure to obtain SA estimates on the geographic equator between 2015.0 and 2019.5 is shown in the left plot of Fig. 10. These SOLA estimates have associated errors of \(0.30.6\mu \mathrm{T}/\mathrm{yr}^2\) and kernel widths \(\approx 42^{\circ }\). The right plot shows similar CHAOS7.2 model predictions for SH degrees \(n \in [1,10]\). Both TLplots shows the similar largescale features, for instance, the features under central America from longitudes \(240^{\circ }\) to \(320^{\circ }\), which are elongated compared with Fig. 9 due to the change in scale of the yaxis (time). However the SOLA results show significantly more timedependence, revealing features that were smoothed out by the temporal regularization of CHAOS7.2. Changes of sign in the SA within about 1 year can be observed.
Particularly interesting is the appearance in the Pacific region around 2017, at longitudes \(150^{\circ }\) to \(220^{\circ }\), of sidebyside positive and negative intense SA features, that have subsequently drifted westwards. This SA change coincides with the peak in the radial SV field observed in the Pacific region during Swarm time seen in the GVO map (Fig. 5). We note the presence of features in Fig. 10 that appear to drift rapidly both eastwards and westwards, for example from \(160^{\circ }\) East in 2015 to \(220^{\circ }\) East in 2017. Such rapidly drifting behavior of low latitude SA patches is difficult to explain in terms of simple core flow advection processes. They may instead be a signature of wave propagation close to the core surface. A range of possible candidates for fast waves in the core have recently been described, some requiring only a strong magnetic field and rotation (Aubert and Finlay 2019; Gerick et al. 2020) while others rely on the presence of a possible stratified layer at the top of the core (Buffett and Matsui 2019). Though tempting, it may be dangerous to interpret such features that are at the limit of the present spatial resolution and temporal resolution (Gillet 2019). It will be important to assess whether such features remain coherent in the future, as the resolution of the SA increases.
Conclusions
In this article, we have studied global patterns and subdecadal changes in geomagnetic secular variation during the past 20 years. We have shown that continuous coverage of magnetic field measurements from lowEarth orbiting satellite missions is now available during this period, provided one takes advantage of calibrated platform magnetometer data from the CryoSat2 satellite. Using vector magnetic field measurements from Ørsted, CHAMP, Cryosat2 and Swarm we have constructed Geomagnetic Virtual Observatory (GVO) time series that track subdecadal changes in the core field secular variation at satellite altitude, and we have used the Subtractive Optimally Localized Averages (SOLA) technique to study the secular variation of the radial field and its time changes down at the core–mantle boundary. These are local methods whereby field measurements in the vicinity of the location of interest are combined so as to estimate the field at that point; measurements far away from the site of interest, have little or no influence on the field estimates.
Using the GVO method, we derived composite time series of geomagnetic secular variation, spanning nearly 20 years, on a global grid of 300 GVOs. GVO time series derived from IGRF13 and CHAOS6x9 calibrated CryoSat2 data, show similar subdecadal SV features, and comparable levels of scatter. We found a scatter level for radial field SV of \(3.5\,\hbox {nT}/\hbox {yr}\) for the CryoSat2 GVOs compared with \(2.5\,\hbox {nT}/\hbox {yr}\) for CHAMP and \(1.0\,\hbox {nT}/\hbox {yr}\) for Swarm. Comparing GVOs with overlapping epochs from 2014 to 2018, derived from calibrated CryoSat2 data and Swarm data, we find similar subdecadal SV changes, thus confirming the possibility of using CryoSat2 for core field studies. In our 20yearlong composite GVO records, we observe fluctuations in the radial SV field of up to 20 nT/yr in amplitude occurring at low latitudes over time periods of 5–10 years. For instance, we see a rapid change of slope in the radial field SV over Indonesia around 2014, over South America and the South Atlantic region around 2007, 2011 and 2014, and in the central Pacific around 2017. Some of these events have previously been discussed in ground observatory records (Brown et al. 2013; Torta et al. 2015). They have the distinct “\(\Lambda\)” and “V”shapes that are often associated with geomagnetic jerks, but as indicated in earlier investigations (Olsen and Mandea 2007, 2008; Chulliat and Maus 2014) we find these events are localized rather than global features.
Using the SOLA technique, we mapped the radial field SV directly at the core–mantle boundary using satellite data, computing spatially localized averages of the SV timeaveraged over specified windows. Taking differences between consecutive SV estimates we obtained estimates for the secular acceleration (SA) at the CMB. Using only CryoSat2 measurements we are able to successfully map the SA at the CMB, down to spatial averaging widths of \(\approx 42^{\circ }\), corresponding approximately to SH degree 10 or length scales of 2500 km at the CMB. Comparing SV and also SA field maps at the CMB, derived from CryoSat2 and Swarm measurements, the same features can be identified having similar amplitudes and latitude/longitude extent.
In time–longitude plots of SOLAbased radial field SA estimates along the geographic equator at the CMB we find strong SA features, with amplitudes of \(\pm 2.5\mu \mathrm{T}/\mathrm{yr}^2\), under Indonesia at longitudes \(70^{\circ }\) to \(160^{\circ }\) from 2011 to 2014 during CryoSat2 time, under central America at longitudes \(240^{\circ }\) to \(320^{\circ }\) from 2015 to 2019, and in sequences with alternating signs of radial SA under South America and the South Atlantic region at longitudes \(240^{\circ }\) to \(360^{\circ }\) during 20042019. The imaged SA features around 2013 at longitudes \(70^{\circ }\) to \(160^{\circ }\) and \(280^{\circ }\) to \(360^{\circ }\) demonstrates the usefulness of the CryoSat2 measurements. Our results lend further support to a sign changing sequence of SA, for length scales down to 2500 km, at longitudes \(240^{\circ }\) to \(360^{\circ }\) from 2005 to 2019; this has been noticed in previous studies (Chulliat et al. 2015; Alken et al. 2020a). We have shown it is possible to increase the temporal resolution of SA estimates during Swarm era, compared to that seen for example in the CHAOS7 model, by computing SA estimates from the differences of consecutive SOLA SV estimates derived using 1yearly time windows. We find similar coherent structures as seen in TL plots constructed using 2yearly time window, but also see that changes of sign in the SA occur within 1 year. In the central Pacific region at longitudes \(150^{\circ }\) to \(220^{\circ }\) we find strong positive and negative SA features appearing sidebyside in late 2017 that drift westwards until 2020. The results presented in Figs. 9 and 10, demonstrate that estimates of core field SA different from those found in the CHAOS model can be obtained despite employing similar data selection and external field modeling schemes. This is due to the important role played by the model parameterization and regularization in the SA recovered in the CHAOS model, especially for the small lengthscales and fast time changes which are towards the limits of what can be reliably resolved from the data.
The rapid fluctuations of the core magnetic field described in this study are likely caused by time variations in the motions of the liquid metal outer core. In particular, changes in secular acceleration patterns at low latitudes provide constraints on the equatorial dynamics of the outer core (e.g. Aubert and Finlay 2019; Kloss and Finlay 2019). This is a topic of active research with various possible phenomenon recently proposed including equatorially trapped MAC waves in a stratified layer close to the core surface (Buffett and Matsui 2019; ChiDurán et al. 2020) or the equatorial focusing of hydrodynamic waves driven by turbulent convection deep within the core (Gerick et al. 2020; Aubert and Finlay 2019).
It is undeniable that much core dynamics occurs on timescales either much longer, or much shorter, than can be can be resolved using the available satellite and ground observations (Gillet 2019). However, with 20 years of lowEarth orbit satellite measurements of the vector magnetic field, and with tools similar to those presented here, it is now possible to probe and characterize core field changes with increasing spatial and temporal resolution. We have shown here that platform magnetometer data can contribute to this activity, provided they are appropriately calibrated.
Availability of data and materials
GVO datasets derived from the Swarm, CHAMP and CryoSat2 missions including composite GVO data series, along with software and additional documentation are available from the Swarm DISC GVO project webpage, https://www.space.dtu.dk/english/research/projects/projectdescriptions/geomagneticvirtualobservatories. The CryoSat2 datasets are available at swarmdiss.eo.esa.int in folder #CryoSat2.
Abbreviations
 AT:

Alongtrack
 CHAMP:

CHAllenging Minisatellite Payload
 CHAOS:

CHAMP, Ørsted, and Swarm field model
 CI:

Comprehensive Inversion Model
 CMB:

Core–mantle boundary
 ECEF:

EarthCenteredEarthFixed reference frame
 \(\scriptsize E_m\) :

Merging electric field at the magnetopause
 EW:

East–West
 FGM:

Fluxgate magnetometer
 GCV:

Generalized crossvalidation
 GVO:

Geomagnetic Virtual Observatory
 IGRF:

International Geomagnetic Reference Field
 IMF:

Interplanetary Magnetic Field
 Kp:

K planetary index
 LCS:

Lithospheric model from CHAMP and Swarm
 LEO:

LowEarth orbiting
 nT:

NanoTesla
 QD:

Quasidipole
 RC:

RingCurrent index
 rms:

Rootmean square
 RMM:

Revised monthly mean
 SA:

Secular acceleration
 SH:

Spherical harmonic
 SOLA:

Subtractive Optimally Localized Averages
 SV:

Secular variation
 TL:

Time–longitude
 yr:

year
References
Alken P, Olsen N, Finlay CC (2020a) Coestimation of geomagnetic field and inorbit fluxgate magnetometer calibration parameters. Earth Planets Space 72:1–32. https://doi.org/10.1186/s40623020011639
Alken P et al (2020b) International Geomagnetic Reference Field: the thirteenth generation. Earth Planets Space, in press,. https://doi.org/10.1186/s4062302001288x
Aubert J, Finlay CC (2019) Geomagnetic jerks and rapid hydromagnetic waves focusing at Earth’s core surface. Nat Geosci 12(5):393–398
Backus G, Gilbert F (1968) The resolving power of gross earth data. Geophys J Int 16(2):169–205
Backus G, Gilbert F (1970) Uniqueness in the inversion of inaccurate gross earth data. Philos Trans Ro Soc Lond A Math Phys Eng Sci 266(1173):123–192
Backus G, Parker R, Constable C (1996) Foundations of Geomagnetism. Cambridge Univ. Press, New York
Baerenzung J, Holschneider M, Wicht J, Lesur V, Sanchez S (2020) The kalmag model as a candidate for IGRF13. Earth Planets Space 72(1):1–13. https://doi.org/10.1186/s4062302001295y
Barrois O, Hammer MD, Finlay CC, Martin Y, Gillet N (2018) Assimilation of ground and satellite magnetic measurements: inference of core surface magnetic and velocity field changes. Geophys J Int 215:695–712
Barton G (1989) Elements of Green’s functions and propagation: potentials, diffusion, and waves. Oxford University Press, Oxford
Bouligand C, Gillet N, Jault D, Schaeffer N, Fournier A, Aubert J (2016) Frequency spectrum of the geomagnetic field harmonic coefficients from dynamo simulations. Geophys Suppl Monthly Notices R Astron Soc 207(2):1142–1157
Brown W, Mound J, Livermore P (2013) Jerks abound: an analysis of geomagnetic observatory data from 1957 to 2008. Phys Earth Planet Int 223:62–76
Buffett B, Matsui H (2019) Equatorially trapped waves in Earth’s core. Geophys J Int 218(2):1210–1225
ChiDurán R, Avery MS, Knezek N, Buffett BA (2020) Decomposition of geomagnetic secular acceleration into traveling waves using complex empirical orthogonal functions. Geophys Res Lett. e2020GL087940
Christensen U, Wardinski I, Lesur V (2012) Timescales of geomagnetic secular acceleration in satellite field models and geodynamo models. Geophys J Int 190(1):243–254
Chulliat A, Maus S (2014) Geomagnetic secular acceleration, jerks, and a localized standing wave at the core surface from 2000 to 2010. J Geophys Res Solid Earth 119(3):1531–1543
Chulliat A, Thébault E, Hulot G (2010) Core field acceleration pulse as a common cause of the 2003 and 2007 geomagnetic jerks. Geophys Res Lett. https://doi.org/10.1029/2009GL042019
Chulliat A, Alken P, Maus S (2015) Fast equatorial waves propagating at the top of the Earth’s core. Geophys Res Lett 42(9):3321–3329
Constable C (1988) Parameter estimation in nonGaussian noise. Geophys J Int 94(1):131–142
Domingos J, Pais MA, Jault D, Mandea M (2019) Temporal resolution of internal magnetic field modes from satellite data. Earth Planets Space 71(1):1–17. https://doi.org/10.1186/s4062301809835
Finlay CC, Olsen N, Kotsiaros S, Gillet N, TøffnerClausen L (2016) Recent geomagnetic secular variation from Swarm and ground observatories as estimated in the CHAOS6 geomagnetic field model. Earth Planets Space 68(1):1–18. https://doi.org/10.1186/s4062301604861
Finlay CC, Kloss C, Olsen N, Hammer MD, TøffnerClausen L, Grayver A, Kuvshinov A (2020) The CHAOS7 geomagnetic field model and observed changes in the South Atlantic Anomaly. Earth Planets Space 72(1):1–31. https://doi.org/10.1186/s40623020012529
Fisher R (1953) Dispersion on a sphere. Proc R Soc Lond A Math Phys Eng Sci 217:295–305
Gerick F, Jault D, Noir J (2020) Fast quasigeostrophic magnetocoriolis modes in the earth’s core. Res Lett Geophys. https://doi.org/10.1029/2020GL090803
Gillet N (2019) Spatial and temporal changes of the geomagnetic field: insights from forward and inverse core field models, pp 115–132, Special Publications of the International Union of Geodesy and Geophysics, Cambridge University Press, https://doi.org/10.1017/9781108290135.010
Green PJ, Silverman BW (1993) Nonparametric regression and generalized linear models: a roughness penalty approach. Chapman and Hall, London
Gubbins D, Roberts N (1983) Use of the frozen flux approximation in the interpretation of archaeomagnetic and palaeomagnetic data. Geophys J Int 73(3):675–687
Hammer MD (2018) Local estimation of the Earth’s core magnetic field, Ph.D. thesis, Technical University of Denmark
Hammer MD, Finlay CC (2019) Local averages of the coremantle boundary magnetic field from satellite observations. Geophys J Int 216(3):1901–1918
Hammer MD, Cox G, Brown W, Beggan CD, Finlay CC (2021) Geomagnetic Virtual Observatories: monitoring geomagnetic secular variation with the Swarm satellites. Earth Planets Space. https://doi.org/10.1186/s40623021013579
Kloss C, Finlay CC (2019) Timedependent lowlatitude core flow and geomagnetic field acceleration pulses. Geophys J Int 217(1):140–168
Kloss C, Finlay CC, Olsen N (2021) Coestimating geomagnetic field and calibration parameters: modelling Earth’s magnetic field with platform magnetometer data. Earth Planets Space. https://doi.org/10.1186/s40623020013517
Langel RA (1987) The main field. In: Jacobs JA (ed) Geomagnetism, vol 1. Academic Press, Cambridge, pp 249–5512
Lebedev VI, Laikov D (1999) A quadrature formula for the sphere of the 131st algebraic order of accuracy, In: Doklady. Mathematics, vol. 59, pp. 477–481, MAIK Nauka/Interperiodica
Leopardi P (2006) A partition of the unit sphere into regions of equal area and small diameter. Electron Trans Numer Anal 25(12):309–327
Lesur V, Wardinski I, Hamoudi M, Rother M (2010) The second generation of the GFZ reference internal magnetic model: GRIMM2. Earth Planets Space 62(10):6. https://doi.org/10.5047/eps.2010.07.007
Lesur V, Wardinski I, Baerenzung J, Holschneider M (2017) On the frequency spectra of the core magnetic field Gauss coefficients. Phys Earth Planet Int
Livermore PW, Hollerbach R, Finlay CC (2017) An accelerating highlatitude jet in Earth’s core. Nat Geosci 10(1):62
Mandea M, Olsen N (2006) A new approach to directly determine the secular variation from magnetic satellite observations. Geophys Res Lett. https://doi.org/10.1029/2006GL026616
Mandea M, Holme R, Pais A, Pinheiro K, Jackson A, Verbanac G (2010) Geomagnetic jerks: rapid core field variations and core dynamics. Space Sci Rev 155(1–4):147–175
Oldenburg DW (1984) An introduction to linear inverse theory. IEEE Trans Geosci Remote Sens GE–22(6):665–674
Olsen N, Mandea M (2007) Investigation of a secular variation impulse using satellite data: the 2003 geomagnetic jerk. Earth Planet Sci Lett 255(1):94–105
Olsen N, Mandea M (2008) Rapidly changing flows in the Earth’s core. Nat Geosci 1(6):390
Olsen N, Stolle C (2012) Satellite geomagnetism. Annu Rev Earth Planet Sci 40:441–465
Olsen N, Mandea M, Sabaka TJ, TøffnerClausen L (2009) CHAOS2—a geomagnetic field model derived from one decade of continuous satellite data. Geophys J Int 179(3):1477–1487
Olsen N, Lühr H, Finlay CC, Sabaka TJ, Michaelis I, Rauberg J, TøffnerClausen L (2014) The CHAOS4 geomagnetic field model. Geophys J Int 197(2):815–827
Olsen N et al (2015) The Swarm Initial Field Model for the 2014 geomagnetic field. Geophys Res Lett 42(4):1092–1098
Olsen N, Ravat D, Finlay CC, Kother LK (2017) LCS1: a highresolution global model of the lithospheric magnetic field derived from CHAMP and Swarm satellite observations. Geophys J Int 211(3):1461–1477
Olsen N, Albini G, Bouffard J, Parrinello T, TøffnerClausen L (2020) Magnetic observations from CryoSat2: calibration and processing of satellite platform magnetometer data. Earth Planets Space 72. https://doi.org/10.1186/s40623020011719
Parker RL (1977) Understanding inverse theory. Annu Rev Earth Planet Sci 5(1):35–64
Parker RL (1994) Geophysical Inverse Theory. Princeton University Press, Princeton
Richmond AD (1995) Ionospheric electrodynamics using magnetic apex coordinates. J Geomagn Geoelectr 47(2):191–212
Ritter P, Lühr H, Maus S, Viljanen A (2004) Highlatitude ionospheric currents during very quiet times: their characteristics and predictability. Annales Geophysicae 22:2001–2014
Ropp G, Lesur V, Baerenzung J, Holschneider M (2020) Sequential modelling of the Earth’s core magnetic field. Earth Planets Space 72(1):1–15. https://doi.org/10.1186/s40623020012301
Sabaka TJ, Hulot G, Olsen N (2010) Mathematical properties relevant to geomagnetic field modeling, In: Handbook of Geomathematics, pp 503–538, Springer
Sabaka TJ, TøffnerClausen L, Olsen N (2013) Use of the comprehensive inversion method for Swarm satellite data analysis. Earth Planets Space 65(11):1201–1222. https://doi.org/10.5047/eps.2013.09.007
Sabaka TJ, TøffnerClausen L, Olsen N, Finlay CC (2018) A comprehensive model of the Earth’s magnetic field determined from 4 years of Swarm Satellite Observations. Earth Planets Space 70(1):130. https://doi.org/10.1186/s4062301808963
Sabaka TJ, TøffnerClausen L, Olsen N, Finlay CC (2020) CM6: a comprehensive geomagnetic field model derived from both CHAMP and Swarm satellite observations. Earth Planets Space 72(1):1–24. https://doi.org/10.1186/s40623020012105
Torta JM, PavonCarrasco FJ, Marsal S, Finlay CC (2015) Evidence for a new geomagnetic jerk in 2014. Geophys Res Lett 42(19):7933–7940
Whaler KA, Beggan CD (2015) Derivation and use of core surface flows for forecasting secular variation. J Geophys Res Solid Earth 120(3):1400–1414
Acknowledgements
We thank the GFZ German Research Centre for Geoscience for providing access to the CHAMP MAGL3 data and the European Space Agency (ESA) for providing access to the CryoSat2 and the Swarm L1b data. Highresolution 1min OMNI data were provided by the Space Physics Data Facility (SPDF), NASA Goddard Space Flight Center. We like to thank two anonymous reviewers for comments that helped improve the manuscript.
Funding
MDH and CCF were partly funded by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 772561). The study has also been partly supported by the Swarm+ 4D Deep Earth: Core project, ESA contract no. 4000127193/19/NL/IA and as part of Swarm DISC activities, funded by ESA contract no. 4000109587.
Author information
Affiliations
Contributions
MDH developed the GVO data processing and modeling scheme and drafted the manuscript. MDH and CCF developed the SOLA modeling scheme. NiO calibrated and prepared the CryoSat2 data and developed the original version of the GVO modeling scheme. All authors contributed to the design of the study. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hammer, M.D., Finlay, C.C. & Olsen, N. Applications for CryoSat2 satellite magnetic data in studies of Earth’s core field variations. Earth Planets Space 73, 73 (2021). https://doi.org/10.1186/s40623021013659
Received:
Accepted:
Published:
Keywords
 Geomagnetism
 Secular variation
 Secular acceleration
 Earth’s core
 Platform magnetometer
 CryoSat2 satellite Swarm satellite constellation