Aims. Cosmogenic isotopes provide the only quantitative proxy for analyzing the long-term solar variability over a centennial timescale. While essential progress has been achieved in both ...measurements and modeling of the cosmogenic proxy, uncertainties still remain in the determination of the geomagnetic dipole moment evolution. Here we aim at improving the reconstruction of solar activity over the past nine millennia using a multi-proxy approach. Methods. We used records of the 14C and 10Be cosmogenic isotopes, current numerical models of the isotope production and transport in Earth’s atmosphere, and available geomagnetic field reconstructions, including a new reconstruction relying on an updated archeo- and paleointensity database. The obtained series were analyzed using the singular spectrum analysis (SSA) method to study the millennial-scale trends. Results. A new reconstruction of the geomagnetic dipole field moment, referred to as GMAG.9k, is built for the last nine millennia. New reconstructions of solar activity covering the last nine millennia, quantified in terms of sunspot numbers, are presented and analyzed. A conservative list of grand minima and maxima is also provided. Conclusions. The primary components of the reconstructed solar activity, as determined using the SSA method, are different for the series that are based on 14C and 10Be. This shows that these primary components can only be ascribed to long-term changes in the terrestrial system and not to the Sun. These components have therefore been removed from the reconstructed series. In contrast, the secondary SSA components of the reconstructed solar activity are found to be dominated by a common ≈2400-year quasi-periodicity, the so-called Hallstatt cycle, in both the 14C and 10Be based series. This Hallstatt cycle thus appears to be related to solar activity. Finally, we show that the grand minima and maxima occurred intermittently over the studied period, with clustering near lows and highs of the Hallstatt cycle, respectively.
Abstract
The
Swarm
mission was selected as the 5th mission in ESA’s Earth Explorer Programme in 2004. The mission will provide the best ever survey of the geomagnetic field and its temporal evolution ...that will lead to new insights into the Earth system by improving our understanding of the Earth’s interior and its effect on Geospace, the vast region around the Earth where electrodynamic processes are influenced by the Earth’s magnetic field. Scheduled for launch in 2010, the mission will comprise a constellation of three satellites, with two spacecraft flying sideby- side at lower altitude (450 km initial altitude), thereby measuring the East-West gradient of the magnetic field, and the third one flying at higher altitude (530 km). High-precision and high-resolution measurements of the strength, direction and variation of the magnetic field, complemented by precise navigation, accelerometer and electric field measurements, will provide the necessary observations that are required to separate and model the various sources of the geomagnetic field. This results in a unique “view” inside the Earth from space to study the composition and processes of its interior. It also allows analysing the Sun’s influence within the Earth system. In addition practical applications in many different areas, such as space weather, radiation hazards, navigation and resource management, will benefit from the
Swarm
concept.
This paper presents a compilation of intensity data covering the past 10 millennia (ArcheoInt). This compilation, which upgrades the one of Korte et al. (2005), contains 3648 data and incorporates ...additional intensity and directional data sets. A large majority of these data (∼87%) were acquired on archeological artifacts, and the remaining ∼13% correspond to data obtained from volcanic products. The present compilation also includes important metadata for evaluating the intensity data quality and providing a foundation to guide improved selection criteria. We show that ∼50% of the data set fulfill reasonable reliability standards which take into account the anisotropic nature of most studied objects (potsherds), the stability of the magnetization, and the data dispersion. The temporal and geographical distributions of this sub–data set are similar to those of the main data set, with ∼72% of the data dated from the past three millennia and ∼76% obtained from western Eurasia. Approximately half of the selected intensity data are associated with at least an inclination value. To constrain the axial and full dipole evolution over the past three millennia requires that we avoid any overrepresentation of the western Eurasian data. We introduce a first‐order regional weighting scheme based on the definition of eight widely distributed regions of 30° width within which the selected data are numerous enough. The regional curves of virtual axial dipole moments (VADM) and of mixed VADM‐virtual dipole moments (VDM) averaged over sliding windows of 200 years and 500 years testify for strong contributions from either equatorial dipole or nondipole components. The computation of global VADM and mixed VADM/VDM variation curves, assuming an equal weight for each region, yields a dipole evolution marked by a distinct minimum around 0 B.C./A.D. followed by a maximum around the third‐fourth century A.D. A second minimum is present around the eighth century A.D. This variation pattern is compatible with the one deduced from earlier, more sophisticated analysis based on the inversion of both intensity and directional data. In particular, there is a good agreement among all VADMs and dipole moment estimates over the historical period, which further strengthens the validity of our weighting scheme.
We present a new empirical model of quiet‐time F‐region ionospheric currents and associated magnetic fields. This model is designed to accurately represent these currents and fields at low and mid ...latitudes. For each individual Swarm satellite, the preprocessed data is represented as a non‐potential toroidal magnetic field using the Mie representation in a thin‐shell and spherical harmonic expansions. This approach allows to fully separate spatial and climatological variations as well as to assess the robustness of the model with respect to both measurement errors and data sampling with local time. The obtained model describes the toroidal magnetic fields and the associated radial poloidal electric currents at two distinct altitudes in the ionosphere F region. Clear signatures of low‐ and mid‐latitude interhemispheric field‐aligned currents (IHFACs) are identified. The model reproduces well‐known characteristics of the climatology of IHFACs and provides new insights, for example, on the average daily variations of IHFACs during winter in the Northern Hemisphere. It also well recovers the variations of IHFACs with longitude. The potential driving mechanisms of these variations, such as longitudinal variations of the main field and modulation by upward propagating atmospheric tides, are discussed. The new model can be used to analyze the relationship between atmospheric tides and IHFACs. It can also be used to investigate the connection between the magnetic fields and electric currents from the ionospheric E and F regions in order to improve the separation of these fields as well as our understanding of the overall ionospheric electric current system.
Key Points
A new Swarm‐based model of low‐ and mid‐latitude ionospheric F‐region magnetic fields and electric currents is presented
The model provides a detailed picture of the climatology of low‐ and mid‐latitude interhemispheric field‐aligned currents
The model provides constraints on the relationship between interhemispheric field‐aligned currents and upward propagating atmospheric tides
Data-based modeling of the magnetic field originating in the Earth’s ionosphere is challenging due to the multiple timescales involved and the small spatial scales of some of the current systems, ...especially the equatorial electrojet (EEJ) that flows along the magnetic dip equator. The Dedicated Ionospheric Field Inversion (DIFI) algorithm inverts a combination of Swarm satellite and ground observatory data at mid- to low latitudes and provides models of the solar-quiet (Sq) and EEJ magnetic fields on the ground and at satellite altitude. The basis functions of these models are spherical harmonics in quasi-dipole coordinates and Fourier series describing the 24-, 12-, 8- and 6-h periodicities, as well as the annual and semiannual variations. A 1-D conductivity model of the Earth and a 2-D conductivity model of the oceans and continents are used to separate the primary ionospheric field from its induced counterpart. First results from the DIFI algorithm confirm several well-known features of the seasonal variability and westward drift speed of the Sq current systems. They also reveal a peculiar seasonal variability of the Sq field in the Southern hemisphere and a longitudinal variability reminiscent of the EEJ wave-4 structure in the same hemisphere. These observations suggest that the Sq and EEJ currents might be electrically coupled, but only for some seasons and longitudes and more so in the Southern hemisphere than in the Northern hemisphere.
We present a new application of a multispacecraft method which provides estimates of the full electric current density vector in the ionospheric low‐ and mid‐latitude F region. The method uses the ...three satellites of the Swarm constellation in configurations when they were close to each other. The current density is calculated inside triangular prisms defined by the satellite positions. The technique is inspired by similar approaches such as the curlometer. Here, we propose an alternative mathematical treatment, which involves the use of the well‐known curl‐B technique and a least squares minimization. The proposed formalism provides a rigorous way to propagate the errors through the problem and special care is taken to validate the algorithm and assess the relevance of the derived currents. Such techniques are powerful tools to study the F‐region ionosphere currents as they naturally exclude all the other contributions outside the considered volumes. The F region is known to host currents with a complex spatiotemporal variability. Besides the intrinsic interest of this system itself, understanding the dynamics of these currents is also critical for global geomagnetic field modeling as they are an important source of error in the models. In this paper, we present the full algorithm and test the method on both synthetic and real Swarm satellite data on February 15, 2014. Several sources of error are investigated. Additionally, the results reveal expected characteristics of F‐region interhemispheric field‐aligned currents. Some new characteristics of F‐region currents, potentially associated with dynamo and pressure currents, are also observed.
Key Points
A new application of a multispacecraft method to derive the three components of the current density in the F‐region ionosphere is presented
An assessment of the effect of biases and errors on the current density estimates is provided
Application of the method on February 15, 2014 reveals both expected and new characteristics of F‐region currents
Using observatory data, we report the detection of a geomagnetic jerk in 2007, which we relate to a jump in the second derivative of the geomagnetic field previously noted in satellite data. Although ...not of worldwide extent, this jerk is very intense in the South Atlantic region. Using the CHAOS‐2 model, we show that both this jerk and the previous 2003 jerk are caused by a single core field acceleration pulse reaching its maximum power near 2006.0. This pulse seems to be simultaneously occurring in several regions of the core surface where it corresponds to dominant n = 5 and 6 spherical harmonic modes. Geometrical attenuation explains why the 2003 and 2007 jerks are local and not fully synchronized at the Earth's surface. Our results suggest that this core field acceleration pulse is the relevant phenomenon to be investigated from the point of view of core dynamics, rather than the jerks themselves.
SUMMARY
We study predictions of reversals of Earth’s axial magnetic dipole field that are based solely on the dipole’s intensity. The prediction strategy is, roughly, that once the dipole intensity ...drops below a threshold, then the field will continue to decrease and a reversal (or a major excursion) will occur. We first present a rigorous definition of an intensity threshold-based prediction strategy and then describe a mathematical and numerical framework to investigate its validity and robustness in view of the data being limited. We apply threshold-based predictions to a hierarchy of numerical models, ranging from simple scalar models to 3-D geodynamos. We find that the skill of threshold-based predictions varies across the model hierarchy. The differences in skill can be explained by differences in how reversals occur: if the field decreases towards a reversal slowly (in a sense made precise in this paper), the skill is high, and if the field decreases quickly, the skill is low. Such a property could be used as an additional criterion to identify which models qualify as Earth-like. Applying threshold-based predictions to Virtual Axial Dipole Moment palaeomagnetic reconstructions (PADM2M and Sint-2000) covering the last two million years, reveals a moderate skill of threshold-based predictions for Earth’s dynamo. Besides all of their limitations, threshold-based predictions suggests that no reversal is to be expected within the next 10 kyr. Most importantly, however, we show that considering an intensity threshold for identifying upcoming reversals is intrinsically limited by the dynamic behaviour of Earth’s magnetic field.
Aims. The Sun shows strong variability in its magnetic activity, from Grand minima to Grand maxima, but the nature of the variability is not fully understood, mostly because of the insufficient ...length of the directly observed solar activity records and of uncertainties related to long-term reconstructions. Here we present a new adjustment-free reconstruction of solar activity over three millennia and study its different modes. Methods. We present a new adjustment-free, physical reconstruction of solar activity over the past three millennia, using the latest verified carbon cycle, 14C production, and archeomagnetic field models. This great improvement allowed us to study different modes of solar activity at an unprecedented level of details. Results. The distribution of solar activity is clearly bi-modal, implying the existence of distinct modes of activity. The main regular activity mode corresponds to moderate activity that varies in a relatively narrow band between sunspot numbers 20 and 67. The existence of a separate Grand minimum mode with reduced solar activity, which cannot be explained by random fluctuations of the regular mode, is confirmed at a high confidence level. The possible existence of a separate Grand maximum mode is also suggested, but the statistics is too low to reach a confident conclusion. Conclusions. The Sun is shown to operate in distinct modes – a main general mode, a Grand minimum mode corresponding to an inactive Sun, and a possible Grand maximum mode corresponding to an unusually active Sun. These results provide important constraints for both dynamo models of Sun-like stars and investigations of possible solar influence on Earth’s climate.
In December 2019, the 13th revision of the International Geomagnetic Reference Field (IGRF) was released by the International Association of Geomagnetism and Aeronomy (IAGA) Division V Working Group ...V-MOD. This revision comprises two new spherical harmonic main field models for epochs 2015.0 (DGRF-2015) and 2020.0 (IGRF-2020) and a model of the predicted secular variation for the interval 2020.0 to 2025.0 (SV-2020-2025). The models were produced from candidates submitted by fifteen international teams. These teams were led by the British Geological Survey (UK), China Earthquake Administration (China), Universidad Complutense de Madrid (Spain), University of Colorado Boulder (USA), Technical University of Denmark (Denmark), GFZ German Research Centre for Geosciences (Germany), Institut de physique du globe de Paris (France), Institut des Sciences de la Terre (France), Pushkov Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation (Russia), Kyoto University (Japan), University of Leeds (UK), Max Planck Institute for Solar System Research (Germany), NASA Goddard Space Flight Center (USA), University of Potsdam (Germany), and Université de Strasbourg (France). The candidate models were evaluated individually and compared to all other candidates as well to the mean, median and a robust Huber-weighted model of all candidates. These analyses were used to identify, for example, the variation between the Gauss coefficients or the geographical regions where the candidate models strongly differed. The majority of candidates were sufficiently close that the differences can be explained primarily by individual modeling methodologies and data selection strategies. None of the candidates were so different as to warrant their exclusion from the final IGRF-13. The IAGA V-MOD task force thus voted for two approaches: the median of the Gauss coefficients of the candidates for the DGRF-2015 and IGRF-2020 models and the robust Huber-weighted model for the predictive SV-2020-2025. In this paper, we document the evaluation of the candidate models and provide details of the approach used to derive the final IGRF-13 products. We also perform a retrospective analysis of the IGRF-12 SV candidates over their performance period (2015–2020). Our findings suggest that forecasting secular variation can benefit from combining physics-based core modeling with satellite observations.