Geomagnetic Excursions: Gothenburg & Laschamp Field Flux

Executive Summary & Theoretical Thesis

Dynamo Instability and Dipole Collapse Dynamics

The geodynamo operates as a non-linear, self-exciting magnetohydrodynamic (MHD) system governed by the convective motion of a low-viscosity, high-conductivity liquid iron-nickel alloy within the terrestrial outer core. Secular variation of this field is dominated by an axial dipole geometry that accounts for over eighty percent of the surface field energy during stable chron intervals.

A geomagnetic excursion does not constitute an aborted reversal in the kinematic sense. Rather, it represents an acute catastrophic bifurcation where the dominant dipolar mode suffers sudden convective decoupling. During an excursion, the virtual axial dipole moment (VADM) of the planet collapses toward a residual threshold between five and ten percent of its time-averaged baseline. The convective energy is not destroyed; instead, it is violently redistributed into non-dipolar, higher-degree spherical harmonics—chiefly quadrupolar ($l=2$) and octupolar ($l=3$) configurations.

This transformation fundamentally alters the topological structure of the terrestrial magnetosphere, eroding the protective magnetopause from an equilibrium standoff distance of approximately ten to twelve Earth radii ($R_e$) down to an ultra-compressed envelope terminating inside four Earth radii ($R_e$), exposing the upper stratosphere to direct solar wind impingement. This deep-seated reorganization is fundamentally linked to magnetohydrodynamic geodynamo collapse.

[ AXIAL DIPOLE STABLE STATE ]
           |  VADM ~ 8e22 A*m^2
           |  Standoff Distance: 10-12 Re
           v
     (Turbulent Convective Decoupling)
           |
           v
      [ MULTIPOLAR EXCURSION PHASE ]
           |  VADM < 1e22 A*m^2 (5-10% Baseline)
           |  Standoff Distance: < 4 Re
           v
   [ DIELECTRIC BREAKDOWN & ATMOSPHERIC IONIZATION ]

Thresholds of Geomagnetic Cutoff Rigidity Loss

The primary consequence of this transition from dipolar symmetry to chaotic multipolar topology is the catastrophic reduction of global geomagnetic-cutoff-rigidity. Under normal dipolar conditions, the terrestrial magnetic vector acts as a momentum-dependent momentum-filter for charged cosmic radiation. As the axial dipole moment deteriorates, this deflection mechanism degrades across all latitudes.

Primary galactic cosmic rays (GCRs)—predominantly high-energy protons and alpha particles—together with relativistic solar energetic particles (SEPs), cease to be funneled exclusively toward the polar cusps. The collapse in magnetic rigidity permits high-flux fluxes of gigaelectronvolt (GeV) particles to penetrate equatorial and mid-latitude atmospheres.

The consequences for atmospheric chemistry and electrical conductivity are catastrophic: high-altitude ionization rates increase by orders of magnitude, collapsing the stratospheric ozone layer through catalytic odd-nitrogen ($NO_x$) and odd-hydrogen ($HO_x$) cycles. This cascade allows biocidal ultraviolet-B (UV-B) radiation to reach the surface while concurrently driving the dielectric breakdown of the upper atmosphere. This degradation manifests in widespread alterations to the global electric circuit and high-energy secondary particle showers reaching sea-level strata.

Chronostratigraphic Synchronicity: 41.4 ka and 12.3 ka Anchors

The Late Pleistocene is characterized by two significant excursion events that serve as chronostratigraphic markers: the Laschamp excursion, anchored at approximately 41.4 ka BP, and the Gothenburg excursion, situated at approximately 12.3 ka BP. Far from being isolated geomagnetic phenomena, these excursions correspond with severe global paleoenvironmental stress, abrupt climatic reconfigurations, and hominin cultural turnover.

The Laschamp excursion marks an era of intense biospheric stress, coinciding with the demise of Western Eurasian Neanderthal lineages, the rapid emergence of Aurignacian subterranean adaptations, and pronounced alterations in planetary atmospheric dynamics.

Conversely, the Gothenburg excursion, positioned precisely within the termination phase of the Allerød interstadial and the onset of the Younger Dryas cold reversal, implicates a magnetospheric collapse in the triggering of global climatic destabilization. This boundary aligns with the mechanical investigations documented within Younger Dryas boundary mechanisms.

By establishing a chronostratigraphic framework that correlates geodynamo disruptions with shifts in cosmogenic isotope production rates, paleomagnetic field records provide an empirical foundation for investigating how geomagnetic field collapse shapes earthly evolutionary trajectories.

Historical Lineage & Experimental Precedents

The Volcanic Record: Bonhommet’s Discovery at Chaîne des Puys

The empirical foundation of geomagnetic excursions originated not in continuous marine sedimentary cores, but within the volcanic fields of the Massif Central in south-central France. In the late 1960s, French geophysicists Norbert Bonhommet and Jean Babkine analyzed the remanent magnetization preserved within the trachybasaltic and trachyandesitic lava flows of the Chaîne des Puys volcanic chain.

At the Laschamp and Olby volcanic vents, thermal remanent magnetization (TRM) vectors within the basaltic matrix revealed magnetic inclinations and declinations directed oppositely to the modern geomagnetic field vector, exhibiting full polarity reversal vectors within rocks geologically constrained to the Late Pleistocene.

Prior to Bonhommet and Babkine’s 1967 discovery, geochronological orthodoxy adhered strictly to the Brunhes-Matuyama binary model established by Allan Cox, Richard Doell, and Brent Dalrymple. Under this framework, the Brunhes Normal Chron, spanning from 780 ka BP to the present, was assumed to be an uninterrupted epoch of homogeneous normal polarity.

The presence of reversed remanence within lavas dated by potassium-argon ($^{40}\text{K}$-$^{40}\text{Ar}$) and thermoluminescence to ~40 ka BP disrupted this paradigm. The Laschamp lavas provided physical evidence that the Earth’s field could execute ultra-rapid, transient deviations toward full reversal before returning to the ambient polarity state, inaugurating the study of excursion dynamics in geophysics.

Sedimentological Verification and the Gothenburg Flip Controversy

Following the verification of the Laschamp event, Swedish paleogeophysicist Nils-Axel Mörner identified anomalous paleomagnetic signatures within a high-resolution, Late Weichselian marine and lacustrine sediment sequence retrieved from Gothenburg, Sweden.

Mörner documented a rapid flip in both inclination and declination, terminating precisely at the boundary between the Allerød interstadial and the Younger Dryas stadial—a horizon he dated via bio-stratigraphy, varve counting, and radiocarbon methodologies to approximately 12,350 yr BP.

Mörner proclaimed this interval the “Gothenburg Magnetic Excursion” and later hypothesized that it represented an omnipresent, global stratigraphic marker where the field’s virtual geomagnetic pole (VGP) wandered into the Southern Hemisphere.

The Gothenburg excursion triggered vigorous academic debate. Initial attempts to replicate the Gothenburg profile in contiguous lacustrine systems across northern Europe yielded inconsistent results. Skeptics argued that Mörner’s “flip” arose from sedimentological artifacts, such as:

• Post-depositional remanent magnetization (PDRM) lock-in distortions,
• Mechanical shear deformation within soft clay during coring operations, or
• Local hydrodynamics and slumping.

However, subsequent high-resolution testing using continuous u-channel sediment analysis and relative paleointensity (RPI) tracking across un-sheared marine cores in the Skagerrak, the Norwegian Sea, and lacustrine basins in North America confirmed significant paleomagnetic instability at the terminal Pleistocene. The consensus recognizes the Gothenburg event not as a sustained, globally homogeneous directional flip, but as a short-lived collapse in paleointensity characterized by localized multipolar drift.

[ CONVENTIONAL BRUNHES CHRON ]
   ----------------------------------------- (780 ka - Present: Uniform Normal)
                     |
       (Discovered Volcanic Anomalies)
                     |
                     v
   [ LASCHAMP EXCURSION IDENTIFIED ]
   (Bonhommet & Babkine, 1967: ~41.4 ka BP)
   - Reversed TRM in Massif Central lavas
   - Transient field collapse paradigm
                     |
       (High-Resolution Marine Coring)
                     |
                     v
   [ GOTHENBURG CONTROVERSY & VALIDATION ]
   (Mörner, 1971/1977: ~12.3 ka BP)
   - Intense directional deviations in clays
   - Replicated via global Relative Paleointensity (RPI) stacks

Evolution of Superconducting Rock Magnetometry (SQUID)

The quantification of these short-lived, transient geomagnetic excursions evolved in tandem with advancements in rock magnetic instrumentation. Early investigations relied upon astatic and balanced fluxgate spinner magnetometers. While functional for measuring high-coercivity, high-intensity volcanic basalt cores, these systems lacked the sensitivity to resolve weak, sub-millennial magnetic declination and inclination paths within unconsolidated, organic-poor sediments.

The introduction of the three-axis Superconducting Quantum Interference Device (SQUID) rock magnetometer revolutionized the field. Operating at cryogenic temperatures sustained by liquid helium, SQUID magnetometers measure weak magnetic flux variations with noise floors approaching $10^{-12} \text{ Am}^2$.

SQUID systems permitted continuous, automated pass-through measurements of sediment u-channels, resolving directional shifts and alternating-field demagnetization profiles at sub-centimeter scales.

These analytical advances made it possible to deconvolve true geomagnetic field shifts from lithological variations, providing the observational basis required to demonstrate that the Laschamp and Gothenburg excursions were true geodynamic phenomena rather than sedimentary artifacts.

Mathematical Formalism & Physical Mechanics

Magnetohydrodynamic (MHD) Induction and Helicity Decay

The operational physics of the geodynamo is formulated through the coupling of the Navier-Stokes equations for rotating, stratified fluid motion and Maxwell’s equations within the magnetohydrodynamic (MHD) approximation. The temporal evolution of the magnetic field $\mathbf{B}$ within the electrically conducting outer core is governed by the magnetic induction equation:

$$\frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{u} \times \mathbf{B}) + \eta \nabla^2 \mathbf{B}$$

where $\mathbf{u}$ defines the convective velocity field of the iron-nickel fluid, and $\eta = (\mu_0 \sigma)^{-1}$ represents the magnetic diffusivity of the medium, with $\mu_0$ as the magnetic permeability and $\sigma$ as electrical conductivity. The first term on the right-hand side represents the advective generation of magnetic fields through shear and inductive fluid motion, while the second term accounts for Ohmic dissipation.

In steady-state dynamo operation, the inductive stretching of magnetic lines of force balances Ohmic decay. This process relies on the alpha-effect ($\alpha$), where helical convective plumes—organized by Coriolis forces into column-like structures parallel to the planetary rotation axis—twist toroidal field lines into poloidal loops. The mathematical generation of this kinetic helicity is defined by:

$$H_k = \langle \mathbf{u} \cdot (\nabla \times \mathbf{u}) \rangle$$

During the onset of a geomagnetic excursion, this organized convective framework breaks down. Perturbations in thermal or chemical core-mantle boundary (CMB) heat flux suppress the macroscale column-like flow.

When the local kinetic helicity decays below a critical threshold, the production of poloidal magnetic energy via the $\alpha$-effect is halted. Consequently, the diffusion term $\eta \nabla^2 \mathbf{B}$ dominates the induction equation. The primary axial dipole field undergoes turbulent dissipation over magnetic diffusion timescales:

$$\tau_\eta \sim \frac{L^2}{\pi^2 \eta}$$

where $L$ is the characteristic length scale of the outer core convective cell.

+-------------------------------------------------------------------------+
|                  CORE CONVECTIVE HELICITY STABILITY                     |
|                                                                         |
|   d/dt [B] = Curl(u x B) + eta * Del^2(B)                               |
|   Poloidal field sustained via Alpha-Omega dynamo mechanisms            |
+-------------------------------------------------------------------------+
                                    |
                                    | Heat flux perturbation at CMB;
                                    | Helicity H_k falls below critical threshold.
                                    v
+-------------------------------------------------------------------------+
|                  OHMIC DIFFUSION & DIPOLAR COLLAPSE                     |
|                                                                         |
|   Advective generation term drops; Diffusion term eta*Del^2(B) dominates|
|   Axial Dipole (g_1^0) collapses by 90-95%                              |
+-------------------------------------------------------------------------+
                                    |
                                    | Energy cascades into higher-order harmonics;
                                    | Non-dipolar terms dominate surface topology.
                                    v
+-------------------------------------------------------------------------+
|                  MULTIPOLAR TRANSITION & IONOSPHERIC EXPOSURE           |
|                                                                         |
|   Surface field defined by chaotic multi-polar nodes                    |
|   Störmer cutoff rigidity drops: R_c -> 0 at mid-latitudes              |
+-------------------------------------------------------------------------+

Dipolar vs. Multipolar Energy Partitioning

The planetary magnetic potential $V$ at or above the terrestrial surface is mathematically parameterized using a spherical harmonic expansion in spherical coordinates $(r, \theta, \phi)$:

$$V(r, \theta, \phi) = a \sum_{l=1}^{\infty} \sum_{m=0}^{l} \left( \frac{a}{r} \right)^{l+1} \left[ g_l^m \cos(m\phi) + h_l^m \sin(m\phi) \right] P_l^m(\cos\theta)$$

where $a$ is the mean terrestrial radius, $g_l^m$ and $h_l^m$ are the time-varying Gauss coefficients of degree $l$ and order $m$, and $P_l^m(\cos\theta)$ are the Schmidt semi-normalized associated Legendre polynomials. The degree $l=1$ coefficients ($g_1^0, g_1^1, h_1^1$) define the dipolar field, while $l=2$ defines the quadrupolar, and $l=3$ defines the octupolar components. The total magnetic energy density integrated across the core surface is measured via the Mauersberger-Lowes spectrum:

$$R_l = (l + 1) \sum_{m=0}^{l} \left[ (g_l^m)^2 + (h_l^m)^2 \right]$$

During stable chron configurations, the energy spectrum is concentrated in $R_1$, with the axial dipole term $g_1^0$ outstripping higher-degree terms by more than an order of magnitude. In the Laschamp and Gothenburg excursions, the energy distribution across the Mauersberger-Lowes spectrum shifts.

The $g_1^0$ coefficient drops toward zero, while the quadrupolar terms ($g_2^0, g_2^1, h_2^1, g_2^2, h_2^2$) and octupolar terms ($g_3^m, h_3^m$) retain their amplitudes or amplify through non-linear energy cascading.

The planetary field morphology transitions from an axial dipole to a complex multipolar field. Terrestrial surface topology ceases to feature two discrete polar regions; instead, it develops transient, localized magnetic poles wandering erratically across equatorial and temperate latitudes.

Atmospheric Cosmic Ray Cascades and Dielectric Disruption

When the terrestrial dipole field collapses into a multipolar state, the planetary barrier to incoming high-energy galactic cosmic rays (GCRs) and solar energetic particles (SEPs) deteriorates. Protons with kinetic energies $E_k > 1 \text{ GeV}$ collide with nitrogen and oxygen nuclei within the upper mesosphere and stratosphere, initiating extensive air showers governed by high-energy nuclear spallation processes:

$$p + ^{14}\text{N} \longrightarrow p, n, \pi^\pm, \pi^0, K^\pm, ^{10}\text{Be}, ^{14}\text{C}, ^{36}\text{Cl} + \dots$$

Neutral pions ($\pi^0$) decay into high-energy gamma photons ($\pi^0 \to 2\gamma$), initiating electromagnetic cascades. Charged pions ($\pi^\pm$) decay into muons and neutrinos ($\pi^+ \to \mu^+ + \nu_\mu$), producing penetrating secondary radiation that reaches the planetary surface.

Concurrently, intense ionization of the middle atmosphere by relativistic cascades generates high concentrations of nitric oxide ($NO$) and hydroxyl radicals ($OH$) via dissociation reactions:

$$\text{N}_2 + e^- \longrightarrow 2\text{N}(^2\text{D}, ^4\text{S}) + e^-$$ $$\text{N}(^2\text{D}) + \text{O}_2 \longrightarrow \text{NO} + \text{O}$$

These radical species catalyze ozone destruction cycles:

$$\text{NO} + \text{O}_3 \longrightarrow \text{NO}_2 + \text{O}_2$$ $$\text{NO}_2 + \text{O} \longrightarrow \text{NO} + \text{O}_2$$

The net result is the systematic thinning of the stratospheric ozone column by twenty to forty percent globally, with polar and mid-latitude depletions exceeding seventy percent during concurrent coronal mass ejections.

This ion production increases upper-atmospheric electrical conductivity, collapsing the vertical dielectric field of the global electric circuit and driving significant charge exchange down into the troposphere.

Empirical Evidence & Observational Data

Radionuclide Deposition Profiles: Beryllium-10 and Carbon-14 Spikes

Because cosmogenic radionuclides are generated through the spallation of atmospheric nuclei by cosmic ray cascades, their historical deposition rates act as an inverse proxy for paleomagnetic field intensity. The beryllium-10-radionuclide ($^{10}\text{Be}$), characterized by a half-life of $1.387 \times 10^6$ years, is rapidly scavenged by atmospheric aerosols and deposited within polar ice sheets within one to two years of its production, preventing significant oceanic reservoir mixing.

High-resolution analysis of ice cores retrieved from the Greenland Ice Core Project (GRIP), GISP2, and the Antarctic EPICA Dome C and Vostok platforms reveals a prominent anomaly in $^{10}\text{Be}$ deposition flux precisely centered at $41.4 \pm 0.5 \text{ ka BP}$.

During the peak collapse phase of the Laschamp excursion, the $^{10}\text{Be}$ flux rises by more than a factor of two relative to Holocene background averages. Radiocarbon ($^{14}\text{C}$) profiles calibrated across Cariaco Basin marine varves and speleothem records confirm this dynamic, recording a synchronous surge in atmospheric $^{14}\text{C}$ concentration ($\Delta^{14}\text{C}$) that defies explanation by carbon-cycle ocean circulation dynamics alone.

[ HIGH-ALTITUDE SPALLATION CASCADE ]
            p + N/O ---> 10Be + 14C + Spallation Products
                                |
             +------------------+------------------+
             |                                     |
             v                                     v
     [ 10Be METEORIC FLUX ]               [ 14C ATMOSPHERIC SURGE ]
  - Rapid aerosol scavenging (~1-2 yr)   - Integrated into carbon cycle
  - Preserved in polar ice strata        - Sequestered in sub-fossil wood
  - GRIP / EPICA Dome C spikes           - IntCal anomalies & Cariaco varves

Sub-fossil Arboreal Records: The Kauri Tree Atmospheric Benchmark

Sediment and ice core records are often constrained by depositional smoothing and lock-in delays. In 2021, an international team led by Alan Cooper utilized sub-fossil New Zealand swamp kauri (Agathis australis) to construct a calendar-year resolved record of the Laschamp excursion. Preserved in anaerobic peat bogs for over forty millennia, these massive sub-fossil trunks captured atmospheric composition at annual rings throughout the excursion’s duration.

The kauri data revealed an acute event preceding the directional minimum of the Laschamp excursion, designated the “Adams Event” in honor of science communicator Douglas Adams. The analysis documented a continuous, precise tracking of the atmospheric $\Delta^{14}\text{C}$ surge across a 1,700-year window from $42.3\text{ ka to }41.6\text{ ka BP}$.

The data revealed that the most intense atmospheric degradation and ecological stress occurred not during the absolute directional minimum of the field, but throughout the preceding transitional phase. During this period, the dipole field collapsed to under six percent of its normal strength while the magnetic axis drifted across equatorial latitudes, destabilizing weather patterns and driving shifts in the Pacific Intertropical Convergence Zone (ITCZ).

Marine and Lacustrine Sedimentary Relative Paleointensity (RPI)

Beyond volcanic and arboreal records, deep-sea and lacustrine sediment cores provide geographically distributed paleointensity records over deep time. Because natural remanent magnetization (NRM) in sediment varies with lithological concentration and grain-size variations, relative paleointensity (RPI) is calculated by normalizing the NRM against laboratory-induced artificial magnetizations:

• Anhysteretic Remanent Magnetization (ARM), or
• Isothermal Remanent Magnetization (IRM).

The PISO-1500 paleointensity stack (Channell et al., 2009), which integrates high-sedimentation cores from across the Atlantic, Indian, and Pacific oceans, corroborates the global extent of the Laschamp excursion as a deep relative paleointensity minimum.

Similarly, high-resolution continuous measurements from the Black Sea, Lake Baikal, the Skagerrak, and the Western Mediterranean reveal that the Gothenburg excursion at ~12.3 ka BP manifests as a coherent RPI low.

These uniform RPI minimums demonstrate that despite regional differences in paleomagnetic field directions—caused by localized multipolar drift—the overall energy of the geodynamo was globally depressed during both horizons.

Metaphysical Implications & Unified Synthesis

Ionospheric Cavity Perturbations and Schumann Resonances

The planetary electromagnetic environment is framed by the terrestrial waveguide bounded by the conductive lithosphere and the lower ionospheric D-layer (at an altitude of $\sim 60\text{ km}$). Within this cavity, continuous broadband electromagnetic radiation generated by global lightning discharges excites extremely low-frequency (ELF) transverse electromagnetic standing waves: the Schumann resonances. The idealized eigenfrequencies for a lossless, perfectly conducting spherical shell cavity are derived from the planetary radius $R_e$ and the speed of light $c$:

$$f_n = \frac{c}{2\pi R_e} \sqrt{n(n + 1)}$$

producing characteristic fundamental frequencies at $7.83\text{ Hz}, 14.3\text{ Hz}, 20.8\text{ Hz}$, and $27.3\text{ Hz}$.

During a geomagnetic field collapse, this resonant architecture undergoes significant perturbation. The surge in cosmic ray cascades drives high-density ionization down into the mesosphere and upper stratosphere, lowering the effective ionospheric boundary layer from $60\text{ km}$ to below $30\text{ km}$.

This severe boundary compression radically alters cavity volume and boundary-layer surface impedance:

$$Z_s = \sqrt{\frac{i \omega \mu}{\sigma + i \omega \epsilon}}$$

The influx of charges increases conductivity $\sigma$ at lower altitudes, accelerating the attenuation rate (lowering cavity quality factor $Q$) and shifting the fundamental resonant eigenmodes upward.

Because terrestrial biological neuro-electric architectures—including mammalian electroencephalographic (EEG) alpha and theta bands—operate within these exact ELF bandwidths, this ionospheric restructuring during excursions represents a direct mechanism of environmental biophysical stress.

[ STANDARD DIPOLAR CAVITY ]
+---------------------------------------+  Ionospheric D-Layer (~60 km)
|                                       |
|  Cavity Eigenmode: f_1 = 7.83 Hz      |  Transverse EM Standing Waves
|  Low dielectric leakage / High Q      |
+---------------------------------------+  Terrestrial Lithosphere
                   |
                   | Cosmic Ray / SEP Influx Ionizes Mesosphere;
                   | Stratospheric dielectric breakdown occurs.
                   v
       [ EXCURSION COMPRESSED CAVITY ]
+---------------------------------------+
|  Ionospheric Boundary Depressed (<30 km) 
|  f_n Eigenmodes shifted upward        |  Extreme Dissipation / Low Q
|  Enhanced ELF Attenuation Rates       |
+---------------------------------------+  Terrestrial Lithosphere

Paleolithic Cave Substrata: Subterranean Retreats as Radiation Refugia

The chronological intersection between the peak of the Laschamp excursion ($41.4\text{ ka BP}$) and the flowering of Upper Paleolithic cave art across Eurasia is one of the most intriguing convergences in paleoanthropology.

Throughout the Franco-Cantabrian region (exemplified by Chauvet, El Castillo, and Fumane), human populations shifted their spatial habits, venturing hundreds of meters into subterranean karstic labyrinths. Concurrently, archaeological assemblages document a sudden surge in the systematic, heavy processing of red ochre (iron-oxide rich hematite, $\text{Fe}_2\text{O}_3$).

While traditional narratives interpret these subterranean networks through symbolic or ritualistic frameworks, field physics offers an ecological hypothesis: human troglodytism during the Laschamp collapse served as a survival mechanism against extreme cosmic ray flux increases and elevated UV-B exposure.

With the stratospheric ozone layer depleted and high-energy secondary particle showers ionizing the troposphere, surface environments were subjected to elevated mutational and biological pressures.

Deep limestone caverns acted as natural shielding structures, attenuating the ionizing secondary muon flux and providing subterranean radiation refugia.

Concurrently, the extensive topical application of micronized red ochre served a functional purpose as a mineral sunscreen and barrier against cutaneous cellular degradation, mitigating the effects of solar UV-B radiation.

Cataclysmic Archaeoastronomical Realignment Traditions

Beyond the immediate physical and biological impacts, the visual and electromagnetic sky underwent dramatic transformations during these field collapses.

With the dipolar magnetic shield diminished, the planetary auroral oval—typically restricted to geomagnetic latitudes higher than $65^\circ$—broadened and descended toward the equator. During persistent solar energetic particle events, auroral curtains and high-energy synchrotron radiation discharges would have illuminated the night sky globally.

These auroral displays were not the diffuse glows observed today, but dynamic electromagnetic plasma filaments organized along the chaotic field lines of the dominant quadrupolar and octupolar nodes.

These conditions left deep impressions on ancestral memory, preserved in ancient archaeoastronomical traditions and oral mythologies. Traditions collected across geographically disparate indigenous cultures recount historical epochs dominated by “falling skies,” celestial serpents, and intense fire falling from the heavens, often linked to the demise of ancient civilizations or shifts in the stars.

These motifs mirror the real physical phenomena of an excursion: violent coronal mass ejections directly impacting the unshielded atmosphere, accompanied by the radical migration of the virtual geomagnetic pole.

The disorientation of seasonal migratory animal herds—dependent on magnetic inclination vectors for navigation—would have heightened the cultural impression of cosmic collapse, leading to realignments in ancient mythic cosmologies. For complementary analysis of these celestial tracking shifts, see archaeoastronomy and axial precession.

Frequently Asked Questions

Excursion vs. Reversal Mechanics

A geomagnetic excursion and a full geomagnetic polarity reversal differ fundamentally in duration, topological evolution, and final outcome.

A full polarity reversal—such as the Brunhes-Matuyama reversal at 780 ka BP—involves a prolonged collapse of the axial dipole field followed by the dynamic restructuring of the outer core’s global helicity. In a full reversal, the geodynamo reconstructs its axial poloidal geometry with its magnetic vector permanently inverted ($180^\circ$ directional shift) relative to its prior state. This reversed state remains stable for hundreds of thousands to millions of years.

In contrast, a geomagnetic excursion represents a transient instability of the outer core’s convective helicity that fails to achieve a stable inversion. The dipole moment collapses rapidly—dropping by 90% or more—and the virtual geomagnetic pole (VGP) wanders erratically across temperate and equatorial latitudes.

However, this multipolar state is unstable. The geodynamo recovers its original polarity within several centuries to a few millennia without completing a permanent polarity inversion. An excursion can be conceptualized as an uncompleted or aborted reversal where the core’s macroscale convective forces restore the baseline dipolar symmetry before a reversed state can stabilize.

+--------------------------------------------------------------------------+
|                     GEODYNAMO BIFURCATION TAXONOMY                       |
+--------------------------------------------------------------------------+
|  METRIC             | FULL POLARITY REVERSAL    | GEOMAGNETIC EXCURSION  |
+---------------------+---------------------------+------------------------+
|  Dipole Recovery    | Permanent 180° inversion  | Recovers initial state |
|  Time Horizon       | 10^3 to 10^4 years        | 10^2 to 10^3 years     |
|  Chron Impact       | Terminates/starts Chron   | Sub-chron micro-event  |
|  Core State         | Full convective inversion | Localized helicity dip |
+---------------------+---------------------------+------------------------+

Gothenburg Excursion and Younger Dryas Synchronicity

The chronological convergence of the Gothenburg magnetic excursion at $\sim 12.3\text{ ka BP}$ with the Allerød-Younger Dryas boundary is a subject of active research in Quaternary science. While mainstream climatic paradigms attribute the abrupt cooling of the Younger Dryas primarily to the shutdown of the Atlantic Meridional Overturning Circulation (AMOC) via glacial meltwater pulses from Lake Agassiz, paleomagnetic evidence highlights the role of concurrent magnetospheric collapse.

The rapid weakening of the geodynamo to approximately 15–25% of modern baseline levels during the Gothenburg event coincided with an abrupt increase in atmospheric cosmogenic isotope production and a severe compression of the magnetopause.

This weakened magnetospheric shielding allowed solar wind momentum and coronal mass ejection plasma to interact directly with the upper atmosphere. This interaction catalyzed dramatic stratospheric ozone depletion, altered the planetary electric circuit, and perturbed Arctic cloud condensation nuclei formation.

Rather than acting as isolated triggers, the Gothenburg magnetic excursion and high-latitude glaciological destabilizations likely operated in feedback with one another, amplifying the rapid descent into Younger Dryas stadial conditions.

Modern Field Weakening and Multi-Pole Emergence

Instrumental measurements over the past two centuries—beginning with Carl Friedrich Gauss’s absolute geomagnetic intensity observations in the 1830s and continuing through modern satellite arrays like the ESA Swarm mission—confirm that the Earth’s dipolar magnetic field is weakening at an accelerating rate. Globally, the dipole moment has degraded by approximately nine to ten percent over the last 150 years, an attrition rate an order of magnitude faster than typical secular variation.

This decay is concentrated in the South Atlantic Anomaly (SAA), a massive region extending from southern Africa across the South Atlantic to South America. Within this zone, the geomagnetic field intensity is low enough that satellites experience computer system disruptions caused by direct exposure to ionizing Van Allen belt radiation.

MHD core-flow inversions reveal that this anomaly is driven by reverse-flux patches emerging at the core-mantle boundary beneath southern Africa, signaling localized disruptions in convective helicity.

While it remains uncertain whether these modern observations indicate the early stages of a full geomagnetic excursion or a transient secular fluctuation, the observed physics mirrors the initial phases of the Laschamp and Gothenburg collapses: the attenuation of the primary $g_1^0$ axial dipole and the emergence of chaotic multipolar harmonics.

Core Academic Bibliography

• Bonhommet, N., & Babkine, J. (1967). Sur la présence d’aimantations inversées dans la Chaîne des Puys. Comptes Rendus de l’Académie des Sciences, 264, 92–94.
• Channell, J. E. T., Xuan, C., & Hodell, D. A. (2009). Stacking paleointensity and oxygen isotope data for the last 1.5 Myr (PISO-1500). Earth and Planetary Science Letters, 283(1-4), 14–23.
• Cooper, A., Turney, C. S. M., Palmer, J., et al. (2021). A global environmental crisis 42,000 years ago. Science, 371(6531), 811–818.
• Mörner, N.-A. (1977). The Gothenburg Magnetic Excursion. Quaternary Research, 7(3), 413–427.
• Muscheler, R., Beer, J., Kubik, P. W., & Synal, H.-A. (2005). Geomagnetic field intensity during the last 60,000 years based on $^{10}\text{Be}$ and $^{14}\text{C}$ from the GRIP ice core and Cariaco Basin sediments. Quaternary Science Reviews, 24(16-17), 1849–1860.