Abstract.Numerical simulations of the geodynamo have successfully represented many observable characteristics of the geomagnetic field, yielding insight into the fundamental processes that generate magnetic fields in the Earth's core. Because of limited spatial resolution, however, the diffusivities in numerical dynamo models are much larger than those in the Earth's core, and consequently, questions remain about how realistic these models are. The typical strategy used to address this issue has been to continue to increase the resolution of these quasi-laminar models with increasing computational resources, thus pushing them toward more realistic parameter regimes. We assess which methods are most promising for the next generation of supercomputers, which will offer access to O(10 6 ) processor cores for large problems. Here we report performance and accuracy benchmarks from 15 dynamo codes that employ a range of numerical and parallelization methods. Computational performance is assessed on the basis of weak and strong scaling behavior up to 16,384 processor cores. Extrapolations of our weak scaling results indicate that dynamo codes that employ two-or three-dimensional domain decompositions can perform efficiently on up to ∼ 10 6 processor cores, paving the way for more realistic simulations in the next model generation.

Satellite observations reveal short pulses in the second time derivative of the geomagnetic field. We seek to interpret these signals using complex empirical orthogonal functions (CEOFs). This methodology decomposes the signal into traveling waves, permitting estimates for the period, angular wave number, and phase velocity. We recover CEOFs from the CHAOS‐6 model, focusing on three geographic regions with strong secular acceleration. Two regions are confined to the equator, while the third is located under Alaska. We find evidence for both eastward and westward traveling waves with periods between 7 and 20 years. There is also evidence for weaker standing waves with complex spatial patterns. Two of the three regions have waves that are compatible with predictions for waves in a stratified fluid. Our results yield estimates for the structure of fluid stratification at the top of the core.

No abstract