  The flow of liquid metal inside the Earth's core produces the geomagnetic
field and its time variations. Understanding the variability of those deep
currents is crucial to improve the forecast of geomagnetic field variations,
which affect human spacial and aeronautic activities. Moreover, it may provide
relevant information on the core dynamics. The main goal of this study is to
extract and characterize the leading variability modes of core flows over
centennial periods, and to assess their statistical robustness. To this end, we
use flows that we invert from two geomagnetic field models (gufm1 and COV-OBS),
and apply Principal Component Analysis and Singular Value Decomposition of
coupled fields. The quasi geostrophic (QG) flows inverted from both geomagnetic
field models show similar features. However, COV-OBS has a less energetic mean
and larger time variability. The statistical significance of flow components is
tested from analyses performed on subareas of the whole domain. Bootstrapping
methods are also used to extract significant flow features required by both
gufm1 and COV-OBS. Three main empirical circulation modes emerge,
simultaneously constrained by both geomagnetic field models and expected to be
robust against the particular a priori used to build them (large scale QG
dynamics). Mode 1 exhibits three large vortices at medium/high latitudes, with
opposite circulation under the Atlantic and the Pacific hemispheres. Mode 2
interestingly accounts for most of the variations of the Earth's core angular
momentum. In this mode, the regions close to the tangent cylinder and to the
equator are correlated, and oscillate with a period between 80 and 90 years.
Each of these two modes is energetic enough to alter the mean flow, sometimes
reinforcing the eccentric gyre, and other times breaking it up into smaller
circulations. The three main circulation modes added together to the mean flow
account for about 70% of the flows variability, 90% of the root mean square
total velocities, and 95% of the secular variation induced by the total flows.
Direct physical interpretation of the computed modes is not straightforward.
Nonethe-less, similarities found between the two first modes and time/spatial
features identified in different studies of core dynamics, suggest that our
approach can help to pinpoint the relevant physical processes inside the core
on centennial timescales.
