pith. machine review for the scientific record. sign in

arxiv: 2604.18553 · v1 · submitted 2026-04-20 · ⚛️ physics.space-ph · astro-ph.EP· physics.geo-ph· physics.plasm-ph

Recognition: unknown

On the curlometer measurement of field-aligned and perpendicular currents in low Earth orbit: Swarm observations and whole geospace simulations

R Gajewski , RT Desai , B Hnat , D Lin , MW Dunlop , M Fillion , G Hulot , Shreedevi P R
show 8 more authors
M-T Walach E Panov J-M Leger T Jager D Fischer W Magnes JA Blake T Etchells
Authors on Pith no claims yet

Pith reviewed 2026-05-10 02:39 UTC · model grok-4.3

classification ⚛️ physics.space-ph astro-ph.EPphysics.geo-phphysics.plasm-ph
keywords curlometerfield-aligned currentsSwarmtetrahedral configurationnumerical instabilitymagnetosphere-ionospherecurrent densitymeso-scale
0
0 comments X

The pith

Time-shifted magnetic measurements from spacecraft tetrahedra diverge from true field-aligned currents even at scales of hundreds of kilometers.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tests the curlometer method on both Swarm satellite data and full geospace simulations to measure field-aligned and perpendicular currents. It shows that field-aligned currents change on short time scales, so estimates made from positions separated in time do not match the actual current density even when the tetrahedron spans hundreds of kilometers. Poorly shaped tetrahedra also produce false perpendicular current signals because small errors get amplified in the calculation. Quality checks on the tetrahedron geometry can reduce these errors and allow usable field-aligned current values when one face stays aligned with the local magnetic field. The results emphasize that true simultaneous four-point measurements are needed to study how currents connect the magnetosphere and ionosphere.

Core claim

Even at meso-scales of hundreds of kilometres, time-shifted FAC estimates can diverge significantly from ground truth, and poor tetrahedral configurations produce spurious perpendicular currents due to numerical instability in the inversion process. This can be mitigated using appropriate quality metrics and high-quality FAC reconstructions still achieved with a tetrahedral face well-aligned to the local magnetic field.

What carries the argument

The curlometer technique that derives current density from the curl of the magnetic field sampled at four points forming a tetrahedron.

Load-bearing premise

The magnetic field measurements from the four points can be treated as simultaneous and the current density is approximately uniform across the tetrahedron volume at the scales examined.

What would settle it

A direct comparison in which time-shifted curlometer estimates from tetrahedra of hundreds of kilometers match the known current density in the simulations would falsify the reported divergence.

Figures

Figures reproduced from arXiv: 2604.18553 by B Hnat, D Fischer, D Lin, E Panov, G Hulot, JA Blake, J-M Leger, M Fillion, M-T Walach, MW Dunlop, R Gajewski, RT Desai, Shreedevi P R, T Etchells, T Jager, W Magnes.

Figure 1
Figure 1. Figure 1: Illustration of the curlometer tetrahedron: Nodes denoted A, B, C and D represent spacecraft [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Left: Cartesian components of the magnetic field residuals (SM) for Swarm spacecraft A B and C collected for an auroral passing over the northern hemisphere on 17 June 2014. Right: The corresponding current density estimates obtained using the curlometer decomposed into the field-aligned current (FAC) component and the remaining perpendicular magnitude, for the spacecraft tetrahedron configurations ABCCp (… view at source ↗
Figure 3
Figure 3. Figure 3: Left: Part of the Swarm trajectory (17 June 2014) used to evaluate correlations between estimates, shown in the ySM-zSM plane (top) and xSM-ySM plane (bottom), expressed in terms of the Earth’s radius RE. Nodes Cp and Cpp denote the positions of spacecraft C delayed by 22 and 26 s respectively. The distance between the barycentre positions of the two tetrahedra in the enlarged view is around 7.5 km. The Ri… view at source ↗
Figure 4
Figure 4. Figure 4: Left: Part of the Swarm trajectory (17 June 2014) used to evaluate correlations between estimates, shown in the ySM-zSM plane (top) and xSM-ySM plane (bottom), expressed in terms of the Earth’s radius RE. Nodes Cp and Ap denote the positions of spacecraft C and A delayed by 22 s respectively. The distance between the barycentre positions of the two tetrahedra in the enlarged view is around 35 km. The Right… view at source ↗
Figure 5
Figure 5. Figure 5: Radial, R, meridional, θ, and azimuthal, ϕ components of the simulated magnetic field perturbations (on a geographic grid) corresponding to the conditions on 23 March 2024 at 400 km altitude in the northern hemisphere. The plot shows the evolution of the perturbations and the trajectories of Swarm spacecraft from 29 May 2014, 13:31:40-13:44:10 UT, used to sample them. 9 [PITH_FULL_IMAGE:figures/full_fig_p… view at source ↗
Figure 6
Figure 6. Figure 6: Field-aligned (top) and the magnitude of the perpendicular components (second from top) of [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Field-aligned and perpendicular components of the current density estimates (top) for the 29 [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
read the original abstract

Measuring field-aligned currents (FACs) using magnetic field observations provides a powerful means to probe the multi-scale interactions between the magnetosphere, ionosphere and thermosphere. In this study, we apply the curlometer technique to Swarm spacecraft observations and to simulations of the coupled magnetosphere-ionosphere system. We begin by correlating current density curlometer estimates derived from Swarm tetrahedra with varying spatial scales and barycentre locations. This confirms an apparent departure from stationarity for FACs at spatio-temporal scales below 100 km where measurements appear highly uncorrelated. We then analyse simulated magnetic perturbations, where true four-point measurements are available. This shows how, even at meso-scales of hundreds of kilometres, time-shifted FAC estimates can diverge significantly from this ground truth. In both observational and simulated data we find poor tetrahedral configurations can produce spurious perpendicular currents due to numerical instability in the inversion process. This can be mitigated using appropriate quality metrics and high-quality FAC reconstructions still achieved with a tetrahedral face well-aligned to the local magnetic field. These results highlight the dynamic nature of FACs at large as well as small scales, and underscore the substantial advantages of true four-point observations for their accurate analysis.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The paper applies the curlometer technique to Swarm constellation magnetic field data and to coupled magnetosphere-ionosphere simulations. It reports that FAC estimates from tetrahedra become highly uncorrelated at scales below 100 km, that time-shifted four-point estimates diverge from simulation ground-truth currents even at meso-scales of hundreds of km, and that poorly conditioned tetrahedra generate spurious perpendicular currents via numerical instability in the inversion; these artifacts can be reduced by quality metrics and by aligning a tetrahedron face with the local magnetic field. The work concludes that true simultaneous four-point sampling offers substantial advantages over time-shifted reconstructions.

Significance. If the central comparisons hold, the manuscript supplies direct, simulation-validated evidence that the simultaneity and uniformity assumptions underlying the curlometer break at observationally relevant scales. The use of independent ground-truth current density from the geospace model is a clear strength, as is the explicit link between inversion-matrix condition number and spurious perpendicular currents. The practical guidance on quality metrics and tetrahedron orientation will be useful to the community analyzing multi-spacecraft current measurements in low Earth orbit.

major comments (2)
  1. [Simulation analysis section] The central claim that time-shifted FAC estimates diverge significantly from ground truth even at meso-scales rests on the simulation comparison; however, the quantitative metric (e.g., RMS difference, correlation coefficient, or fractional error) used to establish 'significant' divergence is not stated, making it difficult to judge the scale dependence.
  2. [Swarm observations and correlation analysis] The Swarm tetrahedron correlation results are load-bearing for the non-stationarity conclusion; the manuscript must specify the exact tetrahedron selection criteria, barycentre filtering, and any post-hoc time-shift handling to allow assessment of possible selection bias.
minor comments (3)
  1. [Abstract] The abstract refers to 'whole geospace simulations' without naming the model or its key resolution and boundary conditions; this information should appear in the first paragraph of the methods.
  2. [Throughout] Notation for the current-density vector components (J∥ vs J⊥) is occasionally inconsistent between text and figure labels; a single, explicit definition early in the paper would improve clarity.
  3. [Figure captions] Figure captions describing tetrahedron geometries should include the condition number of the inversion matrix or the face-normal angle to B for each example shown.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and positive review, including the recommendation for minor revision. We address each major comment below and will revise the manuscript to incorporate the requested clarifications and quantitative details.

read point-by-point responses

  1. Referee: [Simulation analysis section] The central claim that time-shifted FAC estimates diverge significantly from ground truth even at meso-scales rests on the simulation comparison; however, the quantitative metric (e.g., RMS difference, correlation coefficient, or fractional error) used to establish 'significant' divergence is not stated, making it difficult to judge the scale dependence.

    Authors: We agree that explicit quantitative metrics are needed to support the claim of divergence. The simulation analysis section currently relies on side-by-side visual comparisons of time-shifted curlometer FACs against the model's ground-truth current densities, along with qualitative descriptions of increasing mismatch at smaller scales. In the revised manuscript we will add the Pearson correlation coefficient and normalized RMS difference (RMS divided by the mean absolute current density) between the time-shifted estimates and the simulation truth, computed and plotted as functions of tetrahedron scale. These metrics will be included in a new panel or supplementary table to allow readers to assess the scale dependence quantitatively. revision: yes

  2. Referee: [Swarm observations and correlation analysis] The Swarm tetrahedron correlation results are load-bearing for the non-stationarity conclusion; the manuscript must specify the exact tetrahedron selection criteria, barycentre filtering, and any post-hoc time-shift handling to allow assessment of possible selection bias.

    Authors: We acknowledge that the current description of the Swarm analysis is insufficiently detailed for full reproducibility and bias assessment. The manuscript states that we correlate curlometer estimates from tetrahedra with varying spatial scales and barycentre locations, but does not list the precise thresholds. In the revised version we will expand the methods section to specify: (i) the tetrahedron selection criteria (maximum spacecraft separation, minimum eigenvalue ratio of the geometry matrix, and quality factor threshold); (ii) the barycentre filtering (latitude/longitude bounds and altitude range); and (iii) confirmation that all estimates use the four simultaneous Swarm measurements with no additional post-hoc time shifting. These additions will enable readers to evaluate potential selection effects on the reported decorrelation below 100 km. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper establishes its claims through direct empirical comparisons between curlometer outputs and independent references: Swarm tetrahedron correlations at varying scales, plus simulated magnetic perturbations where the geospace model supplies separate ground-truth current density for simultaneous four-point sampling. Divergence of time-shifted estimates and spurious perpendicular currents from poor tetrahedra are diagnosed via the inversion matrix condition number and quality metrics, which are standard numerical diagnostics rather than fitted parameters or self-defined quantities. No load-bearing step reduces by construction to an input, self-citation chain, ansatz, or renamed empirical pattern; the derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The curlometer rests on the linear approximation of the magnetic field gradient and the assumption that displacement currents are negligible at these frequencies. No new free parameters or invented entities are introduced.

axioms (2)
  • standard math Magnetic field can be linearly interpolated across the tetrahedron volume
    Standard assumption of the curlometer technique invoked when computing the curl from four-point measurements.
  • domain assumption Displacement current is negligible compared with conduction current at the frequencies of interest
    Implicit in applying Ampere's law without the displacement term to low-frequency magnetospheric currents.

pith-pipeline@v0.9.0 · 5581 in / 1393 out tokens · 36110 ms · 2026-05-10T02:39:32.665203+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 30 canonical work pages

  1. [1]

    Aschelhoug & Co., Christiania,

    Birkeland, K.,The Norwegian Aurora Polaris Expedition, 1902–1903, 998 pp., H. Aschelhoug & Co., Christiania,

    1902

  2. [2]

    Hesse, The substorm current wedge in mhd simulations,Journal of Geophysical Research: Space Physics,118, 3364–3376, doi:10.1002/jgra.50187,

    Birn, J., and M. Hesse, The substorm current wedge in mhd simulations,Journal of Geophysical Research: Space Physics,118, 3364–3376, doi:10.1002/jgra.50187,

    work page doi:10.1002/jgra.50187

  3. [3]

    Desai, R.T., M.P.Freeman, J.P.Eastwood, J.W.B.Eggington, M.O.Archer, Y.Y.Shprits, etal., Inter- planetary shock-induced magnetopause motion: Comparison between theory and global magnetohydro- dynamic simulations,Geophysical Research Letters,48, e2021GL092,554, doi:10.1029/2021GL092554,

    work page doi:10.1029/2021gl092554

  4. [4]

    Southwood, K.-H

    Dunlop, M., D. Southwood, K.-H. Glassmeier, and F. Neubauer, Analysis of multipoint magnetometer data,Advances in Space Research,8(9), 273–277, doi:https://doi.org/10.1016/0273-1177(88)90141-X,

    work page doi:10.1016/0273-1177(88)90141-x

  5. [5]

    W.,Ionospheric Multi-Spacecraft Analysis Tools: Approaches for Deriving Ionospheric Pa- rameters, 1 ed., Springer International Publishing, doi:10.1007/978-3-030-26732-2,

    Dunlop, M. W.,Ionospheric Multi-Spacecraft Analysis Tools: Approaches for Deriving Ionospheric Pa- rameters, 1 ed., Springer International Publishing, doi:10.1007/978-3-030-26732-2,

    work page doi:10.1007/978-3-030-26732-2

  6. [6]

    Dunlop, M. W., J. Y. Yang, Y. Y. Yang, H. Lühr, and J. B. Cao, Multi-spacecraft current estimates at swarm, inIonospheric Multi-Spacecraft Analysis Tools,ISSI Scientific Report Series, vol. 17, edited by M. W. Dunlop and H. Lühr, Springer, Cham, doi:10.1007/978-3-030-26732-2_5,

    work page doi:10.1007/978-3-030-26732-2_5

  7. [7]

    W., H.-S

    Dunlop, M. W., H.-S. Fu, C. Shen, X. Tan, X.-C. Dong, Y.-Y. Yang, P. Robert, and C. P. Escoubet, Curlometer and gradient techniques: past and future applications,Frontiers in Astronomy and Space Sciences,Volume 11 - 2024, doi:10.3389/fspas.2024.1422341,

    work page doi:10.3389/fspas.2024.1422341 2024

  8. [8]

    W., et al., Multispacecraft current estimates at swarm,Journal of Geophysical Research: Space Physics,120(10), 8307–8316, doi:https://doi.org/10.1002/2015JA021707,

    Dunlop, M. W., et al., Multispacecraft current estimates at swarm,Journal of Geophysical Research: Space Physics,120(10), 8307–8316, doi:https://doi.org/10.1002/2015JA021707,

    work page doi:10.1002/2015ja021707

  9. [9]

    Ebihara, Y., and T. Tanaka, Generation mechanism of region 1 field-aligned current and energy transfer from solar wind to polar ionosphere,Reviews of Modern Plasma Physics,8, 20, doi: 10.1007/s41614-024-00154-7,

    work page doi:10.1007/s41614-024-00154-7

  10. [10]

    Hulot, P

    Fillion, M., G. Hulot, P. Alken, A. Chulliat, and P. Vigneron, Multispacecraft current density estimates in the low- and mid-latitude f-region ionosphere using the swarm constellation,Journal of Geophysical Research: Space Physics,126, e2020JA028,872, doi:10.1029/2020JA028872,

    work page doi:10.1029/2020ja028872

  11. [11]

    Lester, A

    Forsyth, C., M. Lester, A. Fazakerley, C. Owen, and A. Walsh, On the effect of line current width and relative position on the multi-spacecraft curlometer technique,Planetary and Space Science,59(7), 598–605, doi:https://doi.org/10.1016/j.pss.2009.12.007, cross-Scale Coupling in Plasmas,

    work page doi:10.1016/j.pss.2009.12.007 2009

  12. [12]

    Forsyth, C., I. J. Rae, I. R. Mann, and I. P. Pakhotin, Identifying intervals of temporally invariant field-aligned currents from swarm: Assessing the validity of single-spacecraft methods,Journal of Geophysical Research: Space Physics,122(3), 3411–3419, doi:https://doi.org/10.1002/2016JA023708,

    work page doi:10.1002/2016ja023708

  13. [13]

    Gary, J. B., L. J. Zanetti, B. J. Anderson, T. A. Potemra, J. H. Clemmons, J. D. Wilmingham, and J. R. Sharber, Identification of auroral oval boundaries from in situ magnetic field measurements,Journal of Geophysical Research: Space Physics,103(A3), 4187–4197, doi:10.1029/97JA02814,

    work page doi:10.1029/97ja02814

  14. [14]

    J., and H

    Hallinan, T. J., and H. C. Stenbaek-Nielsen, The connection between auroral acceleration and auroral morphology,Physics and Chemistry of the Earth, Part C: Solar, Terrestrial & Planetary Science, 26(1–3), 169–177, doi:10.1016/S1464-1917(00)00104-5,

    work page doi:10.1016/s1464-1917(00)00104-5 1917

  15. [15]

    Chapman, and N

    Hnat, B., S. Chapman, and N. Watkins, Topology of turbulence within collisionless plasma reconnection, Scientific Reports,13, 18,665, doi:10.1038/s41598-023-45650-x,

    work page doi:10.1038/s41598-023-45650-x

  16. [16]

    W., and J

    Hockney, R. W., and J. W. Eastwood,Computer simulation using particles, Bristol: Hilger, 1988,

    1988

  17. [17]

    Liu, J., V.Angelopoulos, X.Chu, X.Zhou, andC.Yue, Substormcurrentwedgecompositionbywedgelets, Geophysical Research Letters,42, 1669–1676, doi:10.1002/2015GL063289,

    work page doi:10.1002/2015gl063289

  18. [18]

    Lühr, H., and Y.-L. Zhou, Small-scale and mesoscale field-aligned auroral current structures: their spatial andtemporalcharacteristicsdeducedbytheswarmconstellation,Annales Geophysicae,43(2), 447–468, doi:10.5194/angeo-43-447-2025,

    work page doi:10.5194/angeo-43-447-2025 2025

  19. [19]

    14 Lühr, H., J. Park, J. W. Gjerloev, J. Rauberg, I. Michaelis, J. M. G. Merayo, and P. Brauer, Field-aligned currents’ scale analysis performed with the swarm constellation,Geophysical Research Letters,42, 1–8, doi:10.1002/2014GL062453,

    work page doi:10.1002/2014gl062453

  20. [20]

    L., Kinetic alfvén waves and auroral particle acceleration: a review,Reviews of Modern Plasma Physics,7(1), 6, doi:10.1007/s41614-022-00111-2,

    Lysak, R. L., Kinetic alfvén waves and auroral particle acceleration: a review,Reviews of Modern Plasma Physics,7(1), 6, doi:10.1007/s41614-022-00111-2,

    work page doi:10.1007/s41614-022-00111-2

  21. [21]

    G., and J

    Merkin, V. G., and J. G. Lyon, Effects of the low-latitude ionospheric boundary condition on the global magnetosphere,Journal of Geophysical Research,115(A10), A10,202, doi:10.1029/2010JA015461,

    work page doi:10.1029/2010ja015461

  22. [22]

    Merkin, V. G., H. Arnold, S. Bao, J. Garretson, D. Lin, J. G. Lyon, et al., Jhuapl/kaiju: Mage-1.0.5, doi:10.5281/zenodo.16818682,

    work page doi:10.5281/zenodo.16818682

  23. [23]

    Milan, S.E., L.B.N.Clausen, J.C.Coxon, andetal., Overviewofsolarwind–magnetosphere–ionosphere– atmosphere coupling and the generation of magnetospheric currents,Space Science Reviews,206, 547– 573, doi:10.1007/s11214-017-0333-0,

    work page doi:10.1007/s11214-017-0333-0

  24. [24]

    Panov, E. V., W. Baumjohann, R. Nakamura, J. M. Weygand, B. L. Giles, C. T. Russell, et al., Continent- wide r1/r2 current system and ohmic losses by broad dipolarization-injection fronts,Journal of Geo- physical Research: Space Physics,124, 4064–4082, doi:10.1029/2019JA026521,

    work page doi:10.1029/2019ja026521

  25. [25]

    Lühr, and J

    Ritter, P., H. Lühr, and J. Rauberg, Determining field-aligned currents with the swarm constellation mission,Earth, Planets and Space,65, 1285–1294, doi:10.5047/eps.2013.09.006,

    work page doi:10.5047/eps.2013.09.006 2013

  26. [26]

    Sorathia, K. A., et al., Multiscale magnetosphere-ionosphere coupling during stormtime: A case study of thedawnsidecurrentwedge,Journal of Geophysical Research: Space Physics,128(11), e2023JA031,594, doi:10.1029/2023JA031594,

    work page doi:10.1029/2023ja031594

  27. [27]

    Tian, S., and et al., Evidence for alfvén waves powering auroral arc via a static electric potential drop, Nature Communications,17(1), 297, doi:10.1038/s41467-025-65819-4,

    work page doi:10.1038/s41467-025-65819-4

  28. [28]

    Toffoletto, F. R., S. Sazykin, R. W. Spiro, and R. A. Wolf, Inner magnetospheric modeling with the rice convection model,Space Science Reviews,107(1), 175–196, doi:10.1023/A:1025532008047,

    work page doi:10.1023/a:1025532008047

  29. [29]

    R.,Computational Methods for Inverse Problems, Society for Industrial and Applied Mathe- matics, doi:10.1137/1.9780898717570,

    Vogel, C. R.,Computational Methods for Inverse Problems, Society for Industrial and Applied Mathe- matics, doi:10.1137/1.9780898717570,

    work page doi:10.1137/1.9780898717570

  30. [30]

    Blagau, and L

    Vogt, J., A. Blagau, and L. Pick, Robust adaptive spacecraft array derivative analysis,Earth and Space Science,7, e2019EA000,953, doi:10.1029/2019EA000953,

    work page doi:10.1029/2019ea000953

  31. [31]

    15 Wang, W., D. Lin, K. Pham, and V. Merkin, Mage simulations of thermosphere and ionosphere re- sponses to subauroral polarization streams (saps), in2025 International Conference on Electromagnet- ics in Advanced Applications (ICEAA), pp. 1309–1312, Palermo, Italy, doi:10.1109/ICEAA65662.2025. 11305966,

    work page doi:10.1109/iceaa65662.2025 2025

  32. [32]

    Zhang, B., K. A. Sorathia, J. G. Lyon, V. G. Merkin, J. S. Garretson, and M. Wiltberger, Gamera: A three-dimensional finite-volume mhd solver for non-orthogonal curvilinear geometries,The Astrophys- ical Journal Supplement Series,244(1), 20, doi:10.3847/1538-4365/ab3a4c,

    work page doi:10.3847/1538-4365/ab3a4c

  33. [33]

    Zhou, Y.-L., and H. Lühr, Properties of large-amplitude kilometer-scale field-aligned currents at au- roral latitudes, as derived from swarm satellites,Annales Geophysicae,43, 667–686, doi:10.5194/ angeo-43-667-2025,

    2025