Recent geomagnetic secular variation from Swarm and ground observatories as estimated in the CHAOS-6 geomagnetic field model

We use more than 2 years of magnetic data from the Swarm mission, and monthly means from 160 ground observatories as available in March 2016, to update the CHAOS time-dependent geomagnetic field model. The new model, CHAOS-6, provides information on time variations of the core-generated part of the Earth’s magnetic field between 1999.0 and 2016.5. We present details of the secular variation (SV) and secular acceleration (SA) from CHAOS-6 at Earth’s surface and downward continued to the core surface. At Earth’s surface, we find evidence for positive acceleration of the field intensity in 2015 over a broad area around longitude 90°E that is also seen at ground observatories such as Novosibirsk. At the core surface, we are able to map the SV up to at least degree 16. The radial field SA at the core surface in 2015 is found to be largest at low latitudes under the India–South-East Asia region, under the region of northern South America, and at high northern latitudes under Alaska and Siberia. Surprisingly, there is also evidence for significant SA in the central Pacific region, for example near Hawaii where radial field SA is observed on either side of a jerk in 2014. On the other hand, little SV or SA has occurred over the past 17 years in the southern polar region. Inverting for a quasi-geostrophic core flow that accounts for this SV, we obtain a prominent planetary-scale, anti-cyclonic, gyre centred on the Atlantic hemisphere. We also find oscillations of non-axisymmetric, azimuthal, jets at low latitudes, for example close to 40°W, that may be responsible for localized SA oscillations. In addition to scalar data from Ørsted, CHAMP, SAC-C and Swarm, and vector data from Ørsted, CHAMP and Swarm, CHAOS-6 benefits from the inclusion of along-track differences of scalar and vector field data from both CHAMP and the three Swarm satellites, as well as east–west differences between the lower pair of Swarm satellites, Alpha and Charlie. Moreover, ground observatory SV estimates are fit to a Huber-weighted rms level of 3.1 nT/year for the eastward components and 3.8 and 3.7 nT/year for the vertical and southward components. We also present an update of the CHAOS high-degree lithospheric field, making use of along-track differences of CHAMP scalar and vector field data to produce a new static field model that agrees well with the MF7 field model out to degree 110.


Introduction
The Earth's intrinsic magnetic field is gradually changing as a result of motional induction and Ohmic dissipation processes taking place within its metallic core. This phenomenon, called geomagnetic secular variation (SV), has been well documented, but poorly understood, for centuries (e.g. Gellibrand 1635; Hansteen 1819). In 1980, MAGSAT provided the first truly global set of vector field observations. Combined with novel regularized inversion techniques, this enabled the structure of field at the core surface to be estimated for the first time with some confidence (Langel et al. 1980;Shure et al. 1982). Unfortunately, the MAGSAT mission lasted less than a year, so inferences concerning SV were limited. It has only been in the past decade, thanks to the Ørsted and CHAMP missions, that it has become possible to map the largescale patterns of the SV directly at the core surface (Lesur et al. 2008;Olsen et al. 2006). A picture has emerged of gradual (decadal) variations in SV punctuated by localized pulses of secular acceleration (SA) on shorter interannual timescales Olsen et al. 2014). SA pulses provide an unexpected new window into the dynamics of the core, and we are still in the early stages of their study. We presently lack detailed knowledge of their morphology and their time dependence, and our understanding is severely limited by the relatively short time window for which there has been global monitoring from space.
A new opportunity for studying SV is today provided by the Swarm mission. Launched on 22 November 2013, it consists of three dedicated low-Earth-orbit satellites, each simultaneously measuring the near-Earth magnetic field. After more than 2 years of operation, Swarm data are starting to provide valuable new constraints on the time-varying SV. In this article, we present investigations of SV as observed by the Swarm satellites, and at ground observatories, in 2015 as part of a new time-dependent geomagnetic field model, called CHAOS-6, that also includes data from previous satellite missions (Ørsted, CHAMP, and SAC-C).
CHAOS-6 is the latest generation of the CHAOS series of global geomagnetic field models developed by Olsen et al. (2006Olsen et al. ( , 2009Olsen et al. ( , 2010Olsen et al. ( , 2014. Ten months of Swarm data (up to September 2014) were included in the previous version, CHAOS-5 , a model that was primarily designed for producing candidate field models for IGRF-12. With more than 2 years of Swarm data now available, given there have been advances in the use of spatial differences (gradients) in field modelling (Kotsiaros et al. 2014;Olsen et al. 2015), and because a geomagnetic jerk happened in 2014 (Torta et al. 2015), there is now a clear need to update the CHAOS model series and particularly its timedependent part, as CHAOS-6.
The CHAOS model series aims to estimate the internal geomagnetic field at the Earth's surface with high resolution in space and time. It is derived primarily from magnetic satellite data, although ground-based activity indices and observatory monthly means are also used. It includes a parameterization of the quiet-time, near-Earth magnetospheric field, but there is no explicit representation of the ionospheric field or fields due to magnetosphere-ionosphere coupling currents. A limitation of the CHAOS models series is that its validity is restricted to after 1999, when the Ørsted satellite was launched.
Other models with continuous time dependence are available for studying core field variations on longer timescales, although these typically have lower resolution in both space and time. The gufm1 model (Jackson et al. 2000) is the definitive source for the historical field from 1590 to 1990. A more recent alternative spanning 1840 to 2010 is the COV-OBS model (Gillet et al. 2013) [see also Gillet et al. (2015a) for a version updated to 2015]. These models contain only a small amount of satellite data and are predominantly constrained by observatory data during the twentieth and twenty-first centuries. A rather different approach to modelling the recent field is provided by the comprehensive model series (Sabaka et al. 2004(Sabaka et al. , 2015. The latest versions CM4 and CM5 cover 1960-2002and 2000, and involve simultaneous estimation of fields from a large number of sources including quiet-time ionospheric currents and magnetosphere-ionosphere coupling currents. This requires a much larger number of free parameters than is the case in the CHAOS models. A series of models of similar complexity to the CHAOS model, but derived only from CHAMP and ground observatory data, is the GRIMM series of models (Lesur et al. 2008(Lesur et al. , 2010. The latter study is particularly interesting because it proposed field models with time dependence controlled by co-estimated core flows.
The main purpose of this article is to present CHAOS-6, providing a reference for users regarding its construction. In the "Data" section, we detail the input data, while the model parameterization and estimation scheme are described in the "Field modelling" section. Model results and related discussion are presented in the "Results and discussion" section, including the fit to Swarm and CHAMP satellite data as well as ground observatory SV in the sections "Fit to satellite data" and "Fit to secular variation estimates from ground observatories" , respectively. The field and SV at Earth's surface are described in the sections "Power spectra of field, SV and SA at Earth's surface" and "Time changes in magnetic intensity at Earth's surface" . The lithospheric field part of CHAOS-6 is described in the section "CHAOS-6h and the high-degree lithospheric field". The field, SV and SA at the core surface are described in the section "Secular variation and acceleration at Earth's core surface". In the section "An interpretation based on quasi-geostrophic core flows", we present for epoch 2015.0 a quasi-geostrophic flow derived from the CHAOS-6 time-dependent field and SV. A summary and perspectives are offered in the "Conclusions" section.

Data
The database of magnetic observations used to construct CHAOS-6 is essentially an extension of that used by  to construct the CHAOS-5 model in September 2014. The ground observatory vector field data have been updated as available in March 2016 (see the section "Ground observatory data"), while vector and scalar field data from the Swarm constellation up to 30 March 2016 have been included. The data selection criteria for satellite data at high latitudes have also been slightly altered compared with previous versions of the CHAOS model; further details are given in the "Satellite data" section.
A major improvement compared with CHAOS-5 is the inclusion of field spatial differences (i.e. approximate gradients) as data, along-track for both CHAMP and Swarm Alpha, Bravo and Charlie, and also east-west between Swarm's lower satellite pair Alpha and Charlie. Along-track field differences approximate north-south gradients in non-polar regions (Kotsiaros et al. 2014;Sabaka et al. 2015), while the east-west differences provide information on the longitudinal gradient of the field . We constructed along-track field differences from data points from the same satellite track, separated by 15 s. East-west difference data were found by searching Swarm L1b 1 Hz magnetic data for a datum from Swarm Charlie, with latitude closest to a selected Swarm Alpha datum, within a maximum time difference of 50 s. This procedure typically resulted in time shifts of about 10 s between the contributing Alpha and Charlie data.

Ground observatory data
Annual differences of revised observatory monthly means  between January 1997 and December 2015 provide crucial constraints on the SV at fixed points on the Earth's surface. Revised monthly means were derived from the hourly mean values of 160 observatories (for locations and IAGA codes see Fig. 1) that have been quality controlled, checking for trends, spikes and other errors (Macmillan and Olsen 2013). Quasi-definitive data (Peltier and Chulliat 2010;Clarke et al. 2013) were used when possible, for times when definitive data were not yet available; these quasi-definitive data are vital for determining up-to-date secular variation and for comparisons with the latest data from the Swarm mission. Starting from hourly mean values, estimates of the ionospheric (plus induced) field from the CM4 model (Sabaka et al. 2004) and the large-scale magnetospheric (plus induced) field from a preliminary CHAOS-type model, CHAOS-6pre, were subtracted. Then a Huber-weighted monthly mean, including all local times and all disturbance levels, is computed. Taking annual differences, this procedure resulted in 23,466 vector field triples of SV estimates.

Satellite data
The basic Ørsted, CHAMP and SAC-C dataset, sampled at a rate of 1 datum per 60 s, is the same as that employed in the earlier CHAOS-4 ) and  ) models. One difference compared with previous models is that CHAMP vector data were used only when attitude information from both star cameras was available.
Regarding the Swarm data, we used the Swarm Level-1b data product Mag-L, taking the latest available baseline 0408/09 data files in March 2016, for more than 28 months from 26 November 2013 to 30 March 2016. During this time, Swarm Alpha and Charlie descended from 514 to 445 km altitude, and Swarm Bravo, after being pushed up to 531 km altitude, has descended to 503 km. The local time of the ascending node of the three Swarm satellites has passed through more than two-anda half 24-hour cycles, and Swarm Bravo is now separated from Swarm Alpha and Charlie by about 3 h of local time. The nominal 1 Hz data were sub-sampled at 60-s intervals unless no vector field magnetometer (VFM) or star tracker (STR) data were available. In addition, we rejected known disturbed days (for example when satellite manoeuvres occurred) and excluded gross outliers for which the vector field components were more than 500 nT (and the scalar field more than 100 nT) from the predictions of a preliminary field model, CHAOS-6pre. In contrast to the case for CHAOS-5, no rescaling of Swarm vector data was necessary to ensure compatibility with the scalar data , since the L1b baseline 0408/09 data calibration includes a co-estimated sundriven disturbance that reduces rms scalar differences between the ASM and VFM measurements to under 200 pT . After November 2014, calibration of the vector magnetometer on Swarm Charlie has been carried out using scalar field values mapped over from Swarm Alpha; this is necessary due to total failure of the absolute scalar magnetometers on Swarm Charlie.
Following experience with the Swarm Initial Field Model (SIFM, see Olsen et al. 2015) and in preliminary experiments for CHAOS-6, we use different selection criteria for different classes of data: For vector field data, we adopt the same quiet-time, dark, selection criteria that were used for earlier versions of the CHAOS model series, namely (1) sun at least 10 • below the horizon, (2) strength of the field due to the magnetospheric ring current, estimated using the RC index , was required to change by at most 2 nT/h, (3) it was required that the geomagnetic activity index Kp ≤ 2 0 for quasi-dipole (Richmond 1995) latitudes equatorward of ±55 • .
As for earlier versions of the CHAOS model series, only scalar intensity data were used poleward of ±55 • quasi-dipole latitude, and these were selected only when the merging electric field at the magnetopause (averaged over the preceding hour) E m ≤ 0.8 mV/m. In CHAOS-6, E m was calculated using 1-min values of the interplanetary magnetic field (IMF) and solar wind speed from OMNIWeb (http://omniweb.gsfc.nasa.gov), in contrast to earlier versions of the CHAOS model where 5-min mean values were used. In addition, an additional selection criterion that IMF B Z > 0 was introduced in CHAOS-6, motivated by a desire to avoid as far as possible disturbances related to the sub-storm auroral electrojet that are especially prominent when IMF B Z < 0. Scalar data were also used at lower latitudes when attitude data were not available.
In CHAOS-6, we also make use of along-track and east-west differences of scalar data. As for the Swarm Initial Field model, SIFM , scalar field differences were used at all latitudes and for all local times (including sunlit conditions, but excluding dayside equatorial region <±10 • quasi-dipole latitudes), with slightly relaxed quiet-time criteria (RC index required to change by at most 3 nT/h and Kp ≤ 3 0 ). The same selection criteria as for scalar data regarding E m and IMF B Z were applied to scalar field differences at polar latitudes. Scalar data have the advantage of not being directly perturbed by the field-aligned currents that are a major contribution to the unmodelled external fields, particularly at polar latitudes. Olsen et al. (2015) found that including spatial differences of scalar data helped to improve the quality of both lithospheric field and secular variation models.
Along-track differences of vector data from the single satellite mission CHAMP and both along-track and eastwest vector differences from the Swarm mission were also employed. For the vector field differences, we used the same selection criteria as for the vector data itself i.e. only data from dark (sun at least 10 • below the horizon), non-polar (equatorward of ±55 • quasi-dipole latitude) regions when the RC index changed by at most 2 nT/h and Kp ≤ 2 0 were selected.
For the low-degree part of CHAOS-6, called CHAOS-6l, 3 × 920,871 vector data, 942,303 scalar data, 1,793,294 along-track scalar differences, 424,003 east-west scalar differences, 3 × 403,382 along-track vector differences and 3 × 92,842 east-west vector differences were used. The reason for the much larger number of scalar differences, compared with the number of scalar data, is that scalar differences were included for all local times (not just dark regions) and because their quiet-time selection criteria were less strict. As in previous versions of the CHAOS model, all satellite data were also weighted proportional to sin θ, where θ is geographic co-latitude, in order to simulate an equal-area distribution.
Although iteratively reweighting of the data is performed during the model estimation (to implement a robust measure of misfit based on a Huber distribution of errors, see the section "Field modelling"), we also employed an a priori error budget to account for the differences between the satellites. Regarding the scalar field, we assumed an a prior error estimate of 2.5 nT for Ørsted, CHAMP and SAC-C and 2.2 nT for Swarm. An isotropic pointing error estimate of 5 arc seconds was assumed for Swarm, 10 arc seconds for CHAMP (when both star cameras are available) and an anisotropic pointing error of 10 arc seconds and, respectively, 40 (60) arc seconds for after (before) 22 January 2000 for Ørsted. Note that these error estimates include the expected impact of unmodelled fields, which often dominate over instrumental errors.

Model parameterization
The model parameterization for CHAOS-6 follows closely that of CHAOS-5 and CHAOS-4. Since the focus of this article is the time-dependent internal field, we explicitly describe only this part of the model. See Olsen et al. (2014) for a more detailed account of the CHAOS field modelling scheme, including the external model.
represented as the gradient of the scalar potential where a = 6371.2 km is a reference radius, (r, θ , φ) are geographic coordinates and P m n (cos θ) are the Schmidt semi-normalized associated Legendre functions of degree n and order m. Note that we follow the usual geomagnetic convention and refer to B as the magnetic field, though it is strictly the magnetic flux density. In the vacuum, it is related to the magnetic field H by B = µ 0 H where µ 0 is the permeability of free space. g m n (t), h m n (t) are time-dependent Gauss coefficients that are further expanded in a basis of sixth-order B-splines (De Boor 2001) such that and similarly for h m n (t), where k g m n are the spline coefficients estimated for each Gauss coefficient, K = 6 (sixthorder B-splines), B k are the spline basis functions, and we use a 6-month knot spacing with fivefold repeated knots at the endpoints t = 1997.1 and t = 2016.6.
In addition to a time-dependent internal field, we also estimate a static internal field above degree 20. The lowdegree part of CHAOS-6 model that is the focus of this section, CHAOS-6l, was estimated using a maximum degree of 80 (in contrast the high-degree part of the CHAOS-6 model, CHAOS-6h was estimated using a maximum degree of 120-see the section "CHAOS-6h: estimation of the high-degree lithospheric field").
Regarding the external field, as in earlier CHAOS models, we use a representation of fields due to near-Earth magnetospheric sources, e.g. the magnetospheric ring current, in the solar magnetic (SM) coordinate system (up to n = 2, with time dependence for n = 1 parameterized by the external and induced parts of the RC index) and of fields due to remote magnetospheric sources, e.g. magnetotail and magnetopause currents, in geocentric solar magnetospheric (GSM) coordinates (also up to n = 2, but restricted to order m = 0). Additional offset parameters in bins of width 5 days, respectively, 30 days are included for the degree 1 SM terms, for orders m = 0 , respectively m = 1.
We also co-estimate the Euler angles needed to describe the rotation between the vector magnetometer frame and the star imager frame. For Ørsted, this yields two sets of Euler angles (one for the period before 22 January 2000 when the onboard software of the star imager was updated and one for the period after that date), while for CHAMP and each Swarm satellite we solve for Euler angles in bins of 10 days.

Model estimation
Model parameters were estimated using an iteratively reweighted least-squares algorithm making use of Huber weights. Regularization of temporal variations was also included. Specifically, we minimized the cost function where m is the model vector, C is the data error covariance matrix which includes anisotropic errors due to attitude uncertainty (Holme and Bloxham 1996) and (3) The vector of residuals e comprises differences between data and model predictions It involves vector and scalar data, denoted by d obs , and the associated model predictions d mod = Gm, where G is the design matrix for the forward model. For scalar data, G is the forward operator linearized around the present model. The data in CHAOS-6 also include along-track and east-west vector and scalar field differences, denoted by �d obs = d obs (r 2 , t 2 ) − d obs (r 1 , t 1 ) . The associated model predictions are Further details of the implementation of along-track and acrosstrack differences in field modelling are described by Kotsiaros et al. (2014) and Olsen et al. (2015Olsen et al. ( , 2016.
In deriving the CHAOS-6l time-dependent field model, we estimated 28,766 model parameters from 7,481,013 observations. The final model was obtained after 9 iterations.

CHAOS-6h: estimation of the high-degree lithospheric field
The final version of CHAOS-6 was obtained by taking the model coefficients from CHAOS-6l (as described above) and replacing the static field Gauss coefficients above degree n = 24 with the static field coefficients from the CHAOS-6h model, truncated at degree n = 110. CHAOS-6h is a new dedicated model of the high-degree lithospheric field. As for the CHAOS-4h model , that provided the high-degree static field in both CHAOS-4 and CHAOS-5, it was derived only using low-altitude, solar minimum CHAMP data, from August 2008 to September 2010.
In addition to scalar and vector field data, CHAOS-6h makes use of along-track scalar and vector field differences from CHAMP. The data selection criteria for the vector and scalar field data are the same as for CHAOS-6l. However, for CHAOS-6h, identical selection criteria are used for both scalar and vector field differences. Data are selected only if Kp ≤ 3 0 , and |dD st /dt| ≤ 3 nT/h. Both night and dayside data are selected, excluding the dayside equatorial region (<±10 • quasi-dipole latitudes). Regarding the CHAOS-6h model parameterization, a static field up to n = 120 was estimated, with a time-dependent internal field for n ≤ 16 described by a third-order Taylor expansion (quadratic SV). The same bin lengths as in CHAOS-6l were used for the RC baseline correction terms and for the instrument alignment calibration parameters (Euler angles). As for CHAOS-4h ), we applied regularization above degree n = 90 by minimizing the L2 norm of B r at Earth's surface.
In all, 15,636 model parameters were estimated from 3,306,074 CHAMP observations when deriving CHAOS-6h.

Fit to satellite data
The fit of the CHAOS-6l field model to scalar and vector satellite data is generally similar (within 0.15 nT) to the fits achieved by CHAOS-5. For example, the Huberweighted rms misfit between CHAOS-6l and non-polar scalar field data is 2.12 nT for CHAMP in comparison with 2.20nT, 2.18nT and 2.19nT, respectively, for Swarm Alpha, Beta and Charlie. The misfit to the B φ vector component is 2.54 nT for CHAMP and 2.50, 2.47 and 2.50 nT, respectively, for Swarm Alpha, Beta and Charlie.
Field difference data were not included in earlier CHAOS models; CHAOS-6 is the first field model to be derived using along-track spatial differences of vector data from both CHAMP and Swarm. Figure 2 presents histograms of residuals for the vector field differences. Comparing Swarm along-track and east-west differences, the along-track differences (involving measurements made on the same orbit by the same instrument) have Huber-weighted rms misfits of 0.27, 0.27 and 0.34 nT for the radial, north-south and east-west components, compared with 0.47, 0.50 and 0.57 nT for the eastwest vector field differences between Swarm Alpha and Charlie. Despite involving measurements from different satellites, we conclude the east-west vector field differences from Swarm are reliable and an internal field model is able to fit them to a weighted rms level of approximately 0.5 nT. Of course no east-west differences were possible with CHAMP, but we can compare the alongtrack differences. We find Huber-weighted rms misfits of 0.36, 0.36 and 0.40 nT for along-track differences of the radial, north-south and east-west vector field components from CHAMP. The along-track differences of Swarm vector data thus have generally smaller residuals than similar differences constructed with CHAMP data. This augurs well for the future of the Swarm mission as the satellites descend.

Fit to secular variation estimates from ground observatories
The highest quality records of geomagnetic secular variation and its time variability come from ground magnetic observatories, where absolute calibrations are routinely carried out. If we are to use CHAOS-6 (which is primarily determined by fitting satellite data) to study secular variation, it is essential that it also fits the available ground observatory data well. We find Huber-weighted rms misfits of CHAOS-6 to annual differences of ground observatory revised monthly means of 3.80, 3.65 and 3.07 nT/year, respectively, for the radial, north-south and east-west components. Examples of ground observatory secular variation time series, along with CHAOS-6 model predictions, are presented in Fig. 3. The top row shows examples of some complete series spanning 1999 to 2016 from well-established observatories at Honolulu (HON) in the central Pacific, at Dorbes (DOU) in Europe and at Alice Springs (ASP) in Australia. We find that CHAOS-6 provides a good description of the time-dependent SV in all these locations. There are noticeable sub-decadal changes in the SV trends even in the central Pacific where SV is often considered to be less intense. As pointed out by Torta et al. (2015), a geomagnetic jerk (characterized by a 'V' shape in the SV as the SA changes sign) clearly occurred in 2014. This jerk event is generally well captured by the CHAOS-6 model.
In addition to presenting examples from well-established observatories, Fig. 3 also shows shorter time series at three recently established observatories in remote locations, at Gan in the southern Maldives, at King Edward Point (KEP) in South Georgia and at Tristan da Cunha (TDC) in the mid-Atlantic. Note that these plots are zoomed in compared to the previous plots and they cover only the shorter time interval of 2010-2015. CHAOS-6 again satisfactorily fits the data from these newer observatories, even when sharp changes in SV are observed, for example in dB φ /dt at Tristan da Cunha in 2014. Although the fit to KEP in Fig. 3 is visually less impressive than that at TDC, note that it is for the northsouth field component, while the (typically quieter) eastwest component is presented for TDC. The rms weighted residuals for dB r /dt, dB θ /dt, dB φ /dt are, respectively, 1.2, 1.6, 0.9 nT/year for TDC and 1.65, 2.40, 2.09 nT/year for KEP.

Power spectra of field, SV and SA at Earth's surface
In Fig. 4, we present the Lowes-Mauersberger spherical harmonic power spectra for the vector field, its first time derivative (SV) and its second time derivative (SA) at the Earth's surface in 2015. The spectra for the field itself decreases steadily until approximately degree 14, after which it begins to level out. The change from a negative (decreasing) slope to a positive (increasing) slope, which indicates that lithospheric sources are certainly dominating, does not take place until degree 18. At the Earth's surface, the spectrum of the SV also decreases with degree; the slope begins to level out about degree 19, indicative of the noise floor being reached. In contrast to the field and the SV, the SA spectra converges at the surface for CHAOS-6 in 2015, with essentially zero power remaining by its truncation degree 20. This is a consequence of the model regularization, that forces the SA towards zero at the model endpoints and minimizes time changes in the SA throughout, which is stronger at higher degree. The low values of the SA spectrum at high degrees should thus not be taken as indicative of a detection limit for the SA which would be related to the noise spectrum; the detection limit can only be properly assessed in unregularized inversions. The SA power spectrum shows weak peaks at degrees 3, 5, 7 and 9 in 2015. Given the surface spectra are well behaved and not diverging, the entire time-dependent part of the CHAOS-6 model (up to spherical harmonic degree 20) can legitimately be used to map and investigate timedependent secular variation at the Earth's surface.   Fig. 2 Histograms of satellite data residuals from CHAOS-6l, for vector field spatial differences (gradients). Blue colours are radial component differences, red/brown colours are north-south component differences, and green colours are east-west component differences. Top A comparison of Swarm along-track and east-west differences residuals. Darker colours are the along-track differences and brighter colours the east-west differences. Bottom A comparison of alongtrack Swarm residuals (darker colours) and along-track CHAMP residuals (brighter colours, with dots). The histograms have been normalized so that each has the same integrated area

Time changes in magnetic intensity at Earth's surface
It is well known that the magnetic field intensity F at Earth's surface is changing, with the South Atlantic Anomaly growing in size and moving westwards. With CHAOS-6, it is possible to map both trends and accelerations in the field intensity at Earth's surface. In Fig. 5, we present maps of F, dF /dt and d 2 F /dt 2 at Earth's surface in 2015. We calculate dF /dt in 2015 from F in 2015.5 minus F in 2014.5 and similarly d 2 F /dt 2 in 2015 from dF /dt in 2015.5 minus dF /dt in 2014.5.
We find that the field is presently strengthening in general in the eastern hemisphere and weakening in the western hemisphere. This is partly a consequence of the low-intensity South Atlantic Anomaly moving to the west, bringing lower field strengths with it, while stronger fields replace it in the east as it moves away. But in addition between 1999 and 2016 the maxima of field intensity over North America have clearly decreased in amplitude, while the field intensity maxima over Northern Asia has grown. A movie of the evolution of F is available at www. spacecenter.dk/files/magnetic-models/CHAOS-6.
Examining the acceleration of the intensity (d 2 F /dt 2 ), we find a strong positive acceleration is now taking place in the east, in a broad longitudinal sector from 30 • to 120 • E. There is also a notable patch of negative acceleration in field intensity around South-West Africa and a negative acceleration taking place close to Alaska and in the northern Pacific region. Considering a time series of such field intensity acceleration maps from 2000 to 2015 (also available as a movie at www.spacecenter.dk/

Fig. 3
Example fits of CHAOS-6 to secular variation recorded at ground observatories. Top Examples of relatively complete records, dB r /dt at Honolulu, USA (left), dB θ /dt at Dourbes, Belgium (middle) and dB φ /dt at Alice Springs, Australia (right); bottom zooms in to shorter records from recently established stations, dB r /dt at Gan, Maldives (left), dB θ /dt at King Edward Point, South Georgia (middle), and dB φ /dt at Tristan da Cunha (right). Black dots are annual differences of revised monthly means, and red lines are the SV predictions from CHAOS-6 files/magnetic-models/CHAOS-6), we find the intensity acceleration changes dramatically on sub-decadal timescales. For example, a series of prominent oscillations are observed west of southern Africa. The field, the SV and the SA downward continued to the core surface are presented later in the section "Secular variation and acceleration at Earth's core surface" and in Fig. 9.
In order to test the above inferences concerning field intensity changes made using the CHAOS-6 field model, in Fig. 6 we present time series of field intensity changes, based on annual differences of ground observatory revised monthly means (F is in this case calculated from the revised month values of B r , B θ , B φ ). A positive intensity acceleration in 2015 is clearly seen at Novosibirsk, Russia, at eastern longitudes in the northern hemisphere and is also evident although weaker at Niemegk, Germany, and at Learmouth, Australia. A relatively long-term negative acceleration is evident in the rate of field intensity decrease observed in Alaska. Overall, we are satisfied that CHAOS-6 adequately explains the observed trends and accelerations of the recent geomagnetic field intensity.

CHAOS-6h and the high-degree lithospheric field
Turning to the higher degree static field in CHAOS-6 (from CHAOS-6h, see the section "CHAOS-6h: estimation of the high-degree lithospheric field"), Fig. 7 presents a map of the lithospheric part of the radial field (degrees 15-110) along with the power spectrum and degree correlation at the Earth's surface in comparison with MF7 (Maus 2010) and CHAOS-4. CHAOS-6 agrees with MF7 much better than CHAOS-4 whose power spectra begin to show deviations above degree 83, when the degree correlation also drops below 0.85. In contrast, the spectrum for CHAOS-6 remains close to that of MF7 up to degree 110, and only after this, does its degree correlation fall below 0.85. We therefore consider the static field in CHAOS-6 to be reliable up degree 110 and recommend its use to this degree.
The map in Fig. 7 shows the radial field plotted at the Earth's surface considering degrees 16 to 110. The map Units are µT = 10 3 nT, nT/ year and nT/year 2 , respectively. Map projection is Hammer-Aitoff displays well-localized anomalies, especially over the continents. Over the oceans, short-wavelength northsouth linear features are visible, despite their relatively low amplitude. Differences to MF7 are mostly with respect to these features. Some differences, especially around auroral electrojet latitudes south of Australia, are possibly due to disturbed tracks, but it may also be that MF7 lacks some along-track power due to the filtering applied during its construction. It will be interesting to see how this part of the signal develops in future models constructed from Swarm data, as the lower pair of satellites descend and they are better able to resolve the shortwavelength east-west field gradients.
Statistics regarding the Huber-weighted mean and rms misfits of CHAOS-6h to the CHAMP field and field difference data used to construct it are presented in Table 1.

Secular variation and acceleration at Earth's core surface
In order to study the origin of secular variation, it is necessary to downward continue the field to the outer edge of its source region in the core. We carry out the downward continuation, assuming that there are no current sources in the mantle on the timescale of observable secular variation, so the field continues to be described by a potential. The resulting spectra for the field, SV and SA at the core surface in 2015 are presented in Fig. 8.
Above degree 13, we see an upward trend in the field spectrum that we attribute to lithospheric sources. We therefore choose to present maps of the field at the core surface only to degree 13. The SV spectra increases rapidly with degree at first, but levels out above degree 9. It starts to increase more rapidly again above degree 18; plotting maps of the SV at the core surface, we see this is associated with an increase in disorganized noise in maps. We therefore believe the SV in CHAOS-6 is satisfactory out at least to degree 16, and possibly even as far as degree 18. Turning to the SA spectrum, in CHAOS-6 this converges at high degree at the core surface due to the applied regularization. In 2015 (relatively close to the model endpoint), regularization starts to dominate the solution already above degree 9. We nonetheless choose to present the SA at the core surface also to degree 16, since some information on rapid field changes is possible up to this degree, particularly for epochs more distant from the model endpoints. Maps of the radial field to degree 13 as well as radial SV and radial SA to degree 16 at the core surface in 2015 are presented in Fig. 9. Movies showing the time changes of such maps are available at www.spacecenter.dk/files/magnetic-models/CHAOS-6.
We find that regions of intense radial SV at the core surface occur close to edges of patches of strong radial field that can be seen to drift when examining a sequence of maps (or a movie) of the radial field between 1999 and 2016. Intense SV in 2015 is observed to lie in a broad band equatorward of 30 • latitude between longitudes 100 • E and 90 • W. There is also a well-localized negative-positive-negative series of three patches of radial SV visible under Alaska and Siberia; this appears to be a consequence of a very rapid westward movement of the intense high-latitude radial field patches. The SV is also generally large in the longitudinal sector from 60 • to 120 • E, particularly in the northern hemisphere.
Regarding the radial field SA at the core surface in 2015, the most prominent features are a positive-negative pair under India-South-East Asia, a series of strong radial SA patches of alternating sign in the region under northern South America, and a positive-negative pair at high northern latitudes under Alaska-Siberia, that is linked to the evolution of the high-latitude SV patches described above.
In both the radial SV and SA, there is a striking absence of structure in the southern polar region (see also the discussion in Holme et al. 2011;Olsen et al. 2014). Although the Pacific region shows lower amplitude radial SV (again see Holme et al. 2011;Olsen et al. 2014), we note that in 2015 there is strong radial SA in the central Pacific, consistent with the aftermath of the jerk observed in 2014 at Hawaii (see Fig. 3, top left). Although the flows driving SV may be weaker in this region, they nonetheless seem to undergo similar time variations.
Earlier versions of the CHAOS model  as well as independent models based on CHAMP and DMSP data (Chulliat and Maus 2014;Chulliat et al. 2015) have demonstrated that the SA undergoes dramatic changes on sub-decadal timescales, notably exhibiting a series of pulses in amplitude. In Fig. 10, power spectra of the SA at the core surface for a number of epochs, and the L2 norm of the SA at the core surface (e.g. , calculated for different spherical harmonic truncation levels, are presented for CHAOS-6. The applied regularization forces the SA spectra to decay at high degree and it begins to have an influence already between degree 10 and 12, especially close to the model endpoints. We find peaks in the SA norm, indicating pulses of SA, for all the investigated truncation levels, at around 2006, 2009.5 and 2013. The exact time of the pulses depends on the chosen truncation level of the SA, which was usually set to degree 6, 8 or 9 in earlier studies. The relative sizes of the pulses also change with the chosen truncation level. As is also evident from the associated power spectra, the 2006 pulse displayed more power at high degrees (10-15), while the 2013 pulse has relatively more power at lower odd degrees 5, 7, 9. Although each pulse has a different spectral signature, there is always enhanced power in the band of degrees from 5 to 7. Maps and movies of the radial SA at the core surface also show recurring oscillations at particular  Present limitations in our ability to infer the highdegree SA in 2015 are also illustrated in Fig. 10. The power spectra of the SA in 2015 drops rapidly above spherical harmonic degree 9. Looking at the SA norm versus time, we see that this is a consequence of the imposed model end constraints which force the SA towards zero in 2016. The end constraints have less influence on the lower degrees (for example, see the SA norm truncated at n = 6), but a longer time span of data is certainly required in order to better determine the highdegree (n > 9) SA in 2015.

An interpretation based on quasi-geostrophic core flows
One possible interpretation of the observed secular variation is in terms of rotation-dominated (or quasigeostrophic) flows of liquid metal in the outer core. An estimate of the responsible flow may be obtained by inverting the magnetic induction equation evaluated at the surface of the core, where u is the core surface flow, ∇ H · is the horizontal divergence operator and where we have neglected magnetic diffusion on the decadal and shorter timescales that are of interest here (see Finlay et al. (2016) for a discussion of the effects of diffusion on longer timescales).
Here, we present a quasi-geostrophic solution for u obtained using the inversion method of Gillet et al. (2015b), taking as input the CHAOS-6 internal field to degree 13 and its SV to degree 16, evaluated at 1-year intervals between 1999.0 and 2016.0. We impose a columnar flow constraint at the core surface that follows from quasi-geostrophy and incompressibility in the outer core volume (Amit and Olson 2004) and also force the flow to be equatorially symmetric, consistent with core motions that are to leading order axially invariant (Pais and Jault 2008) so that, The core surface flow is expanded into toroidal and poloidal parts where r is the position vector and T and S are toroidal and poloidal scalars that are further expanded using a Schmidt semi-normalized spherical harmonic basis, up to degree and order 28. We consider in (5) temporally correlated SV model errors arising from the interaction of the flow with temporally correlated, but unresolved, small-scale field from degrees 14 to 30. An iterative scheme is employed, updating at each step the flow model covariance matrix using information from an ensemble of solutions. CHAOS-6 does not provide covariance information for the input SV, so we adopt simple diagonal covariances for the SV observation errors. These are deduced from the errors provided by the COV-OBS.x1 field model (see Fig. 4 in Gillet et al. 2015a), with a fit to the SV uncertainties in 2010 extrapolated to degree 16. Further details of our flow inversion scheme may be u φ (θ, φ) = u φ (90 − θ, φ) and u θ (θ, φ) = −u θ (90 − θ, φ). Red colours indicate times of SA pulses; blue colours indicate times when SA is minimum. Black and grey colours are times approaching the model endpoints, where SA at high degree is strongly influenced by the imposed model end constraints. Bottom Quadratic norm of the SA power versus time for different truncation levels, n = 6 (black), n = 9 (blue), n = 12 (green), and n = 15 (red) found in Gillet et al. (2015b). The flow models presented here go beyond those presented by Gillet et al. (2015b) in using SV to higher degree, and in focusing on explaining rapid field variations during the past 17 years when highquality satellite data have been available. A map of the resulting quasi-geostrophic flow in 2015, truncated at degree 16, is presented in Fig. 11. Here the green lines follow imaginary tracers in the flow, with the thickness of the line indicating the strength of the flow. At degree 16, the kinetic energy of the ensemble average flow is greater than 50 % of the kinetic energy of any of the ensemble realizations-see Gillet et al. (2015b), section 4.1 for a discussion of how the ensemble can be used to characterize the reliability of the inferred flow. As shown in Fig. 11, the flow is dominated by an anticyclonic, planetary-scale, eccentric gyre consisting of equatorward flow around 100 • E, that then meanders westward flow in a belt around 20 • -30 • N and S of the equator, and then flows poleward again around 90 • W , before closing with intense westward flow at high latitudes around 65 • -75 • N and S, close to the tangent cylinder that circumscribes the inner core. Broadly similar planetary gyres are found many recent flow inversions (e.g. Amit and Pais 2013;Aubert 2015;Gillet et al. 2015b;Baerenzung et al. 2016). The planetary gyre obtained here is, by construction, equatorially symmetric. Using and 2013.5, corresponding to times of large azimuthal flow acceleration, in between maxima and minima of u φ the high-resolution SV from CHAOS-6, we are able to obtain more detail regarding the small-scale structure of the gyre; the flow in Fig. 11 is presented to degree 16, while, for example, Gillet et al. (2015b) presented flows only to degree 14. We find the flow within the centre of the gyre is surprisingly quiescent, for example in the vicinity of the South Atlantic reverse flux patches. This is despite these lying in the Atlantic hemisphere, which is typically considered to be an active region characterized by high-amplitude SV.
In agreement with the findings of Gillet et al. (2015b), we find a series of prominent non-axisymmetric azimuthal (i.e. east-west, or u φ ) jets close to the equator. We find these jets undergo time-dependent oscillations at some locations, for example at 40 • W at the equator-see Fig. 12. At this location, strong oscillations of the radial field SA are seen at the core surface. We find that pulses of SA i at a particular location correspond to times of large acceleration in u φ , occurring between maxima and minima of u φ , for example at 40 • W at the equator, where SA extrema occurred in 2005.8, 2009 and 2013.5. Fig . 11 shows that azimuthal flows at low latitude are dominated by their non-axisymmetric part; their amplitude is significantly larger than that of the axisymmetric motions that are often interpreted as torsional Alfvén waves (Gillet et al. 2010(Gillet et al. , 2015b in the same sub-decadal period range. Time-longitude plots of u φ at the equator do not show coherent propagation in longitude, but rather standing oscillatory features, with enhanced amplitude at particular locations. Interpretation of quasi-geostrophic flows at low latitudes requires pause for thought. Quasi-geostrophic models in a thin-shell (β -plane) geometry, as is relevant for the atmosphere and oceans, are known to break down at the equator. However, the outer core is a thick shell and recent tests of the quasi-geostrophic approximation in this geometry (comparing inertial modes in quasi-geostrophic models against full 3D solutions) show encouraging agreement, even for equatorially confined modes (Canet et al. 2014;Labbé et al. 2015). Further work is needed to better understand the dynamics of the low-latitude non-axisymmetric jets. For example: what drives such motions, and does the non-axisymmetric Lorentz force play an important role in producing the observed oscillations?
At this stage, it is important to recognize that other hypotheses are possible regarding the nature of the core flows. For example, there is presently a debate concerning whether a stratified layer may exist close to the core surface (e.g. Buffett 2014; Buffett et al. 2016;Chulliat et al. 2015;, inhomogeneous boundary conditions may force departures from equatorial symmetry (e.g. Amit and Pais 2013) or large scales may for some reason dominate the flow (e.g. Bloxham 1988;Whaler and Beggan 2015). Nonetheless, the primary flow structures identified here, in particular equatorward flow in both the northern and southern hemispheres around 100 • E and time-dependent non-axisymmetric westward flow at low latitudes, are sufficient to reproduce the observed rapid field changes, within the uncertainties due to the unresolved small-scale field.

Conclusions
In this article, we have presented the CHAOS-6 field model and used it to analyse recent patterns of geomagnetic secular variation. CHAOS-6 includes more than 2 years of Swarm data and the latest ground observatory magnetic measurements as available in March 2016, along with data from previous satellite missions, and it provides information on geomagnetic secular variation between 1999.0 and 2016.5. It is the first member of the CHAOS field model series to use spatial field differences as data, utilizing along-track differences from both the Swarm and CHAMP satellites and east-west differences between Swarm Alpha and Charlie.
At Earth's surface, we find large-scale patterns of secular acceleration that change on short, sub-decadal, timescales. A geomagnetic jerk that occurred in 2014 is visible in Australia, and in the central Pacific, as well as in Europe. Transient accelerations are also seen in the strengthening and weakening of the field intensity; there has recently been a notable positive acceleration of the field intensity in the Asian longitude sector. CHAOS-6 captures the secular variation at the core surface up to at least spherical harmonic degree 16. Looking at the time derivative of this secular variation, the secular acceleration, we find that it has been dominated by a series of pulses, seen most clearly at low latitudes in the Atlantic sector and also at longitudes close to 100 • E. Inverting the secular variation for a quasi-geostrophic core flow, we find the dominant time-averaged feature is a planetary-scale gyre structure that flows equatorward around 100 • E , then westward at mid-to low latitudes and then poleward around 90 • W, closing with intense westward flow at high latitudes close to the tangent cylinder. Rapid fluctuations are evident in the eastern, equatorward, limb of the gyre. In addition, the quasi-geostrophic flows show prominent oscillations of non-axisymmetric azimuthal jets at low latitudes that provide a possible explanation for localized, oscillatory SA pulses observed in this region, for example near to 40 • W under northern South America.
Longer time series of Swarm data are needed to test and extend the preliminary results reported here for the secular variation and secular acceleration in 2015. The relatively long timescales involved, even for rapid secular acceleration pulses, mean that long-term monitoring