the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 21 May 2026
The direct hand of the magnetosphere in controlling small-scale auroral plasma turbulence: Introducing the renormalization group
Abstract. Studies of high-latitude plasma turbulence in Earth's upper atmosphere fundamentally focus on the differential response of electrons and ions to strong external electric fields generated during geomagnetic storms. Because ions in the E region are heavy and highly collisional, they remain largely tied to the neutral gas, whereas magnetized electrons undergo rapid E×B-drift. Microscopic polarization electric fields are generated when the relative drift velocity between these streaming electrons and the background ions exceeds the local ion-acoustic speed, triggering two-stream plasma instabilities, producing Farley Buneman waves. We propose a new theory that explicitly considers thousands, or millions, of such waves being excited inside a limited volume of space around aurorae, subject to the renormalization group. The resulting theory constitutes an effective field-theory for Farley-Buneman turbulence in the Martin-Siggia-Rose formalism. At the core of this theory is a statistical description of Farley-Buneman waves, where we allow each individual wave to produce a polarization electric field. We treat the sum total of these "micro-fields" that occur inside a turbulent volume as a stochastic variable, or simply noise. That noise, now a thermodynamic property, becomes the basis for anomalous diffusion, and an effective diffusion tensor, and we recover the expression for Bohm diffusion. In support of this theory, we present a large statistical analysis of how auroral electrojet turbulence responds to magnetospheric driving, revealing a clear tendency for the observed number density of turbulent waves to scale linearly with driving power, matching the predictions made by our field theory's overdamped equations of motion. Crucially, the effective field theory offers closed-form calculations of macroscopic transport relations that are uniquely suitable for sub-grid parameterization in space weather modeling, mimicking the success of stochastic parameterization in hydrodynamic Earth-system models. The derivation of these parameterized equations demonstrates how Bohm diffusion arises from a statistical-mechanical treatment of turbulent volumes, at the expense of an explicit treatment of the turbulent cascade. The equations should be investigated further, and in future, they may model the evolution of momentum and energy in numerical treatments of global magnetohydrodynamic circulation, below the scale-sizes normally considered accessible to fast, predictive models.
Competing interests: At least one of the (co-)authors is a member of the editorial board of Annales Geophysicae.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.- Preprint
(7474 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (extended)
-
RC1: 'Comment on egusphere-2026-2345', Anonymous Referee #1, 18 Jun 2026
reply
-
AC1: 'Reply on RC1', Magnus Ivarsen, 12 Jul 2026
reply
We thank the referee for a thorough report. We have revised the manuscript, addressing the many important points brought up by the referee. This includes changes of title, abstract, Section 2.1, and conclusion; we have added an Appendix C (presenting the complete stochastic closure), Appendix D (detailed RG classification with operator table, flow equation, convergent integral), and Appendix E [kinetic-simulation confrontation, as the referee requested, as well as the results from a reduced (internal) Monte-Carlo test]. Section 2.1.2 is now restructured to present the physical interpretation (retaining the pedagogical derivation of the theory) while deferring the formal operator classification to Appendix D's in-text table. We have added a panel to Figure 3 (d), which shows Doppler distribution statistics, and we have addressed the timescales involved in the turbulence itself (empirically determined to die out in 1-2 seconds, meriting the theory's second central motivation, namely FB waves being ephemeral, in addition to their being amplitude-saturated).
The abstract now reads: "In Earth's E-region ionosphere, the differential magnetization of electrons and ions allows strong auroral electric fields to drive Hall currents that, above the ion-acoustic speed, break the plasma up into Farley-Buneman turbulence. We construct an effective field theory for the macroscopic transport arising from this turbulence, treating the aggregate polarization micro-fields of the saturated wave structures as a stochastic variable in the Martin-Siggia-Rose formalism. A renormalization group analysis classifies the macroscopic dynamics (dynamical exponent z=2, infrared-stable Gaussian fixed point), while a stochastic closure fixes the transport coefficients at the Bohm order of magnitude. The theory is confronted with three-dimensional kinetic simulations from the literature, tested internally by a reduced Monte Carlo simulation, and exposed to a direct empirical test. Here, we aggregate 272 million ICEBEAR radar echo detections concurrent with 530,000 electric-field wave spectra from the Arase spacecraft: above a critical magnetospheric wave power, the radar echo detection rate rises linearly with driving power (exponent 0.98 ± 0.09), while echo Doppler speeds remain saturated near the ion-acoustic speed, which selects the amplitude-saturated, echo count-carried response predicted by the effective field theory. The resulting closed-form transport law brings to space weather modeling the sub-grid stochastic parameterization of unresolved turbulence that Earth-system models rely on for analogous scale-separation problems, at the cost of an explicit turbulent cascade."
What follows is an outline of our revison.
POINT 1: Referee: "[The replacement of the collective effect of many polarization fields by a stochastic variable] is not a small technical step... Which physical terms are kept? Which are neglected? What is the ordering argument?"
Response. The revision no longer presents this step as derived from kinetic theory. It is presented as a stochastic closure with its assumptions stated. Section 2.1.1 names the closure and its ordering (epsilon = tau_c / T_driver << 1, with tau_c ~ 1 s and T_driver of order minutes), lists the kept terms (mean advection, damping, acoustic stiffness, Gaussian micro-field statistics) and the dropped ones (three-wave phase correlations, hence the explicit spectral energy flux). Appendix C constructs the action term by term: bare propagators, multiplicative-noise vertex, vanishing Ito tadpole (C_0 = 0), Bourret-Kraichnan closure, and Kubo integral in both correlation regimes. The Gaussian step is the random phase approximation; see also the internal Monte Carlo test of the theory (Appendix E).
POINT 2: Referee: "If the calculation is a phenomenological stochastic coarse-graining, it can still be useful, but it should be presented as such. If it is meant to be a genuine RG derivation, then the RG flow and its physical content have to be made explicit."
Response. The two possibilities are separated. The transport coefficient comes from the closure (Appendix C opens by stating that this part is not an RG flow). The universality statement comes from the RG, given in Appendix D: tracked fields and couplings, noise covariance, coarse-graining operation, operator classification table (to which Section 2.1.2 now refers for the formal exponents, presenting only the physical consequences in the main text), z = 2 selected by marginality of the damping; see the new Appendix D.
POINT 3: Referee: "There should be a numerical test of whether this behavior actually emerges from nonlinear Farley-Buneman dynamics... If new simulations are not included, then the manuscript still has to confront the existing nonlinear FB simulation literature in a much more direct way."
Response. No new kinetic simulations were performed; this is stated. Performing new, state-of-the-art PIC simulations of electrojet turbulence is an arduous task that needs dedicated supercomputing resources. However, the new Appendix E confronts the theory row by row with the 3-D PIC simulations of Oppenheim & Dimant (2013) and Oppenheim et al. (2026): phase-speed clamping, saturation amplitudes, post-growth spectra, anisotropy, transport and conductance modification, heating, structure lifetimes, ensemble statistics, response scaling (Table E1). A reduced Monte-Carlo (code, seed, documentation archived on Zenodo) test demonstrates Gaussian convergence of the summed micro-field, dynamical emergence of the thresholded linear response without the slow-manifold approximation, discrimination between the clamped response (beta = 1.01) and the unclamped alternative (beta = 2.08), and the two-branch Kubo diffusivity with its tau_c-independent mixing-length plateau; all presented in the new Appendix E, where we walk a balance between the two reigning schools of FB theory (St-Maurice et al. v. Dimant et al.).
POINT 4: Referee: "A thresholded linear relation can appear in many driven dissipative systems. It is not, by itself, a unique signature of RG physics, Bohm diffusion, or phase-locking dynamics."
Response. Agreed. The revised paper no longer asks the observation to validate the chain. Each link with the literature and statistical mechanical standards carries its own support (Appendices C, D, E). The observation is assigned one discriminating role, stated as necessary rather than sufficient: the measured beta = 0.98 +/- 0.09 selects the amplitude-clamped, count-carried response and rejects the unclamped alternative (expected beta ~ 2; Appendix E, Figure E1). Simpler interpretations are addressed under Points 5 and 6.
POINT 5: Referee: "The observed correlation does not yet establish direct control. It may reflect common driving... If this cannot be shown, then the observational result should be framed more modestly as consistency with magnetosphere-ionosphere coupling."
Response. We accept the logic and adopt the framing the referee proposes. The title no longer asserts direct control. The Results section states what the correlation can and cannot establish: covariance with a common driver is not excluded. The directional evidence is the lag structure: the ionospheric response in the second case study of our paper (Figures 2 and A1) trails the magnetospheric signal by 3 m 20 s, with the lower-band chorus channel shifted a further 20 s; a common large-scale driver imposes no such frequency-ordered delay, while wave-particle interaction and propagation physics do (see Ivarsen et al. 2025 for a similar argument regarding the commonality of the driver). On sample independence: the captions now state that one-second intervals are temporally clustered and not independent, that the bootstrap interval understates event-level uncertainty, and that the scaling conclusion rests on the two-decade span of the bin medians. The full multivariate control analysis (partial correlation controlling for SME, substorm phase, MLT, precipitation, geometry, with event-level block bootstrap) is deferred to follow-up work. The present claims are qualified and calibrated, and do not depend on its outcome.
POINT 6: Referee: "A larger number of echoes does not automatically mean a larger physical number density of turbulent waves."
Response. Three arguments support the identification. First, the meteor-trail echoes are a null control: the same radar, processing, propagation, and geometry yield a flat detection rate across the activity range (Figure 3a, below the elbow), so detectability changes cannot manufacture the trend. Second, saturation clamps per-structure amplitude; Figure 3d shows the Doppler distribution unchanged when the detection threshold is doubled (267 M vs 101 M echoes; black and red histograms).
POINT 7: Referee: "Bohm diffusion should not be presented as derived from first principles unless the kinetic, numerical, or observational evidence is much stronger."
Response. In the revised paper, the Bohm diffusion result is not presented in this way. The argument rests on two independent estimates: the Einstein-relation closure with T_eff ~ T_e (natural because C_s^2 = k_B(T_e + T_i)/m_i pins the fluctuation energy to ambient temperatures), and the mixing-length value dv_par * ell_c of Appendix C, which is independent of the uncertain tau_c (verified by tracer dispersion in the Monte-Carlo test, Appendix E) and shares Bohm's structural form (T_e/eB = v_th * rho_e). Neither is a kinetic derivation; this is stated. Scale range: L < L_c ~ 350 m at C_s = 500 m/s. Empirical support: the fully turbulent Doppler speeds in the new Figure~3d track the heated ion acoustic speed (approaching 700 m/s under strong driving; St-Maurice et al., 2023), the saturated speed migrating with electron temperature. The faster, narrow-spectrum population near 1000~m/s consists of weakly turbulent near-threshold eigenmodes that ride the ion drift (same reference); these lie outside the saturated ensemble the theory coarse-grains, and the revised Discussion states this domain boundary explicitly.
POINT 8: Referee: "The key missing issue is the measurable phase variable... What would falsify this interpretation?"
Response. The interpretation is a candidate interpretation throughout, deserving cautious treatment. The measurable phase is the Doppler velocity of the saturated waves relative to the driver. The falsifier: the interpretation requires Doppler organization around C_s confined to driving above Delta_E_c; failure of either half falsifies it. Figure 3d presents the above-threshold half (organization around C_s, robust to doubled detection threshold); the below-threshold conditioning on the Arase-concurrent subset is stated as a remaining test. The sentence claiming the pattern is "cleanly describable" by Josephson dynamics has been rewritten; the abstract now makes no synchronization claim; the conclusion's two-interpretation structure is retained and sharpened.
References:
Oppenheim, M. M., and Y. S. Dimant. ‘Kinetic Simulations of 3‐D Farley‐Buneman Turbulence and Anomalous Electron Heating’. Journal of Geophysical Research: Space Physics 118, no. 3 (2013): 1306–18. https://doi.org/10.1002/jgra.50196.
Oppenheim, Meers, Y. Dimant, R. Koontaweepunya, A. Q. Green, and S. Evans. ‘First Kinetic 3-D Simulations of the High-Latitude Electrojet Spanning an Entire Turbulent Flux Tube’. Geophysical Research Letters 53, no. 13 (2026): e2025GL121042. https://doi.org/10.1029/2025GL121042.
Ivarsen, Magnus F., Jean-Pierre St-Maurice, Glenn C. Hussey, et al. ‘Eastward Transients in the Dayside Ionosphere. I. Electrodynamics on Closed Field Lines’. Physical Review E 112, no. 4 (2025): 045204. https://doi.org/10.1103/r6bv-pzlq.
St-Maurice, Jean-Pierre, Devin Huyghebaert, Magnus F. Ivarsen, and Glenn C. Hussey. ‘Narrow Width Farley-Buneman Spectra Above 100 Km Altitude’. Journal of Geophysical Research: Space Physics 128, no. 10 (2023): e2022JA031191. https://doi.org/10.1029/2022JA031191.-
RC2: 'Reply on AC1', Anonymous Referee #1, 14 Jul 2026
reply
I thank the authors for their detailed response and for the substantial revision of the manuscript. The revised version is clearly improved. In particular, the authors have made useful changes to the title, abstract, theoretical framing, and conclusion. I also appreciate that they now separate the stochastic closure from the RG classification more clearly, and that they have added further discussion of the connection with previous kinetic simulations and reduced numerical tests.However, I do not think that all of the main concerns are fully resolved yet. The manuscript is certainly more careful than before, but several key steps in the argument still remain partly phenomenological. The stochastic closure is now described more explicitly, which is helpful, but it is still an assumed closure rather than a derivation from the full E-region Farley–Buneman plasma dynamics. Similarly, the RG analysis appears to classify the reduced stochastic model after this closure has been introduced; it does not by itself establish the physical transport coefficient.I also remain concerned that the observational scaling is still asked to carry too much theoretical weight. The relation between Arase wave power and ICEBEAR echo occurrence is interesting, but the common-driving issue has not been fully resolved. The authors acknowledge this point and defer a full multivariate control analysis to future work. In my view, this means that the causal language and the interpretation of direct magnetospheric control should remain cautious.Another important point is the use of radar echo detection rate as a proxy for the number density or filling factor of saturated turbulent structures. The revised manuscript gives a clearer physical argument for this approximation, but I am not fully convinced that it is demonstrated. Coherent radar echo occurrence can still depend on aspect angle, SNR threshold, scattering volume, altitude, propagation effects, plasma density, and detectability. This proxy may be reasonable, but it should be presented as a working observational proxy rather than as a fully established physical number density.I also think the Bohm-order transport and Adler/Josephson interpretation should be treated with more caution. The revised manuscript has softened some language, which is welcome. Still, the Bohm result should be described as a Bohm-order coefficient arising within the proposed closure, not as a kinetic derivation of Bohm diffusion for Farley–Buneman turbulence. Likewise, the Adler/Josephson interpretation remains an interesting analogy or candidate interpretation, but I do not think it is yet demonstrated physically. In particular, the measurable phase variable and falsifiable observational signatures need to be clearer if this interpretation is kept as a central part of the paper.Overall, I appreciate the authors’ effort and I think the revised manuscript is moving in a better direction. The observational dataset is valuable, and the theoretical idea may be productive. However, I would still recommend major revision before acceptance. The manuscript should more clearly distinguish between what is derived, what is assumed through closure, what is supported by observations, and what remains a speculative or phenomenological interpretation.Citation: https://doi.org/
10.5194/egusphere-2026-2345-RC2
-
RC2: 'Reply on AC1', Anonymous Referee #1, 14 Jul 2026
reply
-
AC1: 'Reply on RC1', Magnus Ivarsen, 12 Jul 2026
reply
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 129 | 60 | 16 | 205 | 6 | 7 |
- HTML: 129
- PDF: 60
- XML: 16
- Total: 205
- BibTeX: 6
- EndNote: 7
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
Major comments
The manuscript addresses an important problem in auroral plasma physics: how small-scale E-region Farley–Buneman turbulence, which is usually not resolved in large-scale space-weather models, may be represented through an effective transport description. The observational material is valuable, and the attempt to connect radar observations, magnetospheric wave activity, stochastic transport, and effective-field-theory language is certainly interesting.
However, the paper currently puts a very large theoretical interpretation on top of an observational result. The observations appear to provide the final empirical piece of the story, but many of the intermediate physical and mathematical steps are not yet shown with the same level of care. Before the paper can claim an RG/effective-field-theory interpretation of Farley–Buneman turbulence, the authors have to show that the whole chain of reasoning is valid, not only that the final observed trend is compatible with it.
My major concerns are as follows.
1. The derivation from E-region plasma physics to the stochastic field theory is not sufficiently established.
The central theoretical step in the manuscript is the replacement of the collective effect of many Farley–Buneman polarization fields by a stochastic variable. This stochastic variable then becomes the basis for the MSR/effective-field-theory formulation. This is not a small technical step. It is one of the main pillars of the paper.
At present, I do not see that this step has been derived convincingly from the actual E-region plasma system. The E-region plasma is collisional, partially magnetized, altitude dependent, and strongly driven. It includes magnetized electrons, weakly magnetized or demagnetized ions, ion-neutral collisions, polarization electric fields, nonlinear saturation, electron thermal effects, finite-altitude structure, and radar-scale irregularities. These details are not secondary if the paper is claiming a physical theory of Farley–Buneman turbulence.
The authors need to show the route from the collisional two-fluid or kinetic equations to the proposed stochastic action. Which physical terms are kept? Which terms are neglected? What is the ordering argument? What is the small parameter, if there is one? Why is the random-phase or Gaussian-noise closure justified? Why should the sum of polarization micro-fields have exactly the stochastic character assumed in the theory?
Without this derivation, the theory remains a plausible and interesting phenomenological construction, but not yet a demonstrated consequence of Farley–Buneman plasma dynamics.
2. The RG structure is not yet shown with enough precision.
The manuscript uses renormalization-group language, MSR formalism, dynamic scaling, fixed points, effective field theory, anomalous diffusion, and universality. These words carry precise meanings. At present, it is not clear enough what is actually being renormalized.
The authors state that their method is not a conventional Wilsonian momentum-shell calculation and that the scale separation is mainly temporal. This clarification is useful, but it also raises the central question: in what precise sense is this an RG calculation, rather than a stochastic closure followed by scaling analysis?
The paper needs to identify the variables, couplings, noise covariance, coarse-graining operation, running parameters, fixed point, and the quantities that are universal versus model-dependent. It also needs to show how the effective diffusion tensor follows from the coarse-graining, rather than being inferred from the form of the reduced equation.
This is not just a matter of language. If the calculation is a phenomenological stochastic coarse-graining, it can still be useful, but it should be presented as such. If it is meant to be a genuine RG derivation, then the RG flow and its physical content have to be made explicit. The present version does not yet seperate these two possibilities clearly enough.
3. Numerical verification of the proposed mechanism is missing.
The proposed mechanism is strong. The manuscript argues that saturated Farley–Buneman structures generate polarization micro-fields, that the ensemble of these fields can be treated as stochastic noise, and that coarse-graining then produces effective diffusion, anomalous resistivity, Bohm-like transport, and a thresholded macroscopic response. This sequence cannot be established by observational scaling alone.
There should be a numerical test of whether this behavior actually emerges from nonlinear Farley–Buneman dynamics. This could be through nonlinear fluid, hybrid, kinetic, or reduced simulations. The important issue is not simply whether the model can reproduce the observed echo trend, but whether the plasma dynamics themselves generate the proposed coarse-grained transport law.
The numerical evidence should address whether the polarization-field statistics have the assumed stochastic form, whether the coarse-grained diffusion tensor has the proposed anisotropic structure, whether Bohm-like scaling appears without being imposed, and whether a thresholded linear response emerges from the nonlinear saturation process.
If new simulations are not included, then the manuscript still has to confront the existing nonlinear Farley–Buneman simulation literature in a much more direct way. The proposed transport coefficients and scaling arguments should be compared against those results. Otherwise, the proposed theory remains an interpretation of the observations, not a demonstrated mechanism.
4. The observational result does not yet validate the full theoretical chain.
The observational result is interesting, but it is only the final empirical link in a much longer theoretical chain. A statistical relation between Arase wave power and ICEBEAR echo occurrence does not by itself validate the sequence from stochastic polarization fields to RG coarse-graining, effective diffusion, Bohm transport, and Adler/Josephson-type dynamics.
A thresholded linear relation can appear in many driven dissipative systems. It is not, by itself, a unique signature of RG physics, Bohm diffusion, or phase-locking dynamics. The manuscript needs to explain why this observed scaling specifically supports the mechanism proposed here, rather than a simpler interpretation based on stronger geomagnetic activity, precipitation, conductance, electric fields, or ordinary threshold behavior of the Farley–Buneman instability.
At the moment, the paper sometimes reads as if agreement with the observed scaling confirms the whole theoretical construction. I do not think that follows. The authors need to demonstrate that the intermediate steps are physically necessary, not only compatible with the final observational trend.
5. Common geomagnetic driving has not been ruled out.
The title and interpretation imply direct magnetospheric control of small-scale auroral turbulence. This is a strong causal statement. But Arase wave power, ICEBEAR echo occurrence, SME/AE activity, precipitation, conductance, electric fields, MLT, season, and substorm phase can all respond to the same larger geomagnetic driver.
For this reason, the observed correlation does not yet establish direct control. It may reflect common driving. This issue is especially important because the observational relation is being used to support a broad theoretical interpretation.
The authors need to show that magnetospheric wave power explains the radar response beyond what is already explained by SME/AE, substorm phase, MLT, precipitation, conductance, radar geometry, and preferably local or proxy electric-field conditions. If this cannot be shown, then the observational result should be framed more modestly as consistency with magnetosphere-ionosphere coupling, not as evidence for the specific RG/effective-field-theory mechanism proposed in the paper.
The uncertainty analysis also needs more care. The one-second intervals and individual echoes are not independent samples. Auroral activity is temporally clustered, and radar echoes inside an active interval are strongly correlated. Event-level, night-level, or storm-level uncertainty would likely be more realistic than treating short-cadence samples as independent.
6. Radar echo count is not yet justified as turbulent-wave number density.
A major step in the manuscript is the interpretation of ICEBEAR echo occurrence or echo count as a proxy for the number density or proliferation of Farley–Buneman turbulent structures. This step is not sufficiently demonstrated.
Coherent radar echoes depend on Bragg scattering geometry, magnetic aspect angle, propagation effects, SNR thresholding, altitude, plasma density, irregularity amplitude, scattering volume, and the detection/classification method. A larger number of echoes does not automatically mean a larger physical number density of turbulent waves. It may also mean better viewing geometry, stronger scattering amplitude, different propagation conditions, or changes in detectability.
This point directly affects the bridge between the observations and the theory. If echo count is not a reliable proxy for turbulent-structure density, then the comparison with the proposed field theory becomes much weaker. The manuscript needs to clearly seperate radar echo occurrence, backscatter power, irregularity amplitude, turbulent energy, scattering volume, and actual number density of Farley–Buneman structures.
7. The Bohm-diffusion interpretation requires independent justification.
The recovery of Bohm-like diffusion is one of the main theoretical claims of the manuscript. However, the derivation seems to rely on an effective-temperature or generalized Einstein-relation closure in a strongly driven, non-equilibrium plasma. This is a major physical assumption, not a minor detail.
The authors need to explain why such a closure should hold for saturated Farley–Buneman turbulence in the auroral E region. What is the physical meaning of the effective temperature? Why can it be identified with, or approximated by, the electron temperature? What part of the turbulent spectrum is being thermalized? What is the connection between the stochastic polarization fields and the resulting transport coefficient? Over what scale range should Bohm-like transport be expected?
The manuscript acknowledges some phenomenological aspects of this closure, but the broader claims still read stronger than the support provided. Bohm diffusion should not be presented as derived from first principles unless the kinetic, numerical, or observational evidence is much stronger.
8. The Adler/Josephson interpretation is still an analogy unless a measurable phase variable is identified.
The Adler equation and Josephson-junction analogy are interesting, but they are not yet demonstrated physically. A thresholded linear response is not enough to establish Josephson-like synchronization or phase slipping in auroral Farley–Buneman turbulence.
The key missing issue is the measurable phase variable. What plasma phase is locking or slipping? How would this phase be observed in radar data, spacecraft data, or numerical simulation? What evidence would distinguish genuine Adler/Josephson-type dynamics from a more ordinary thresholded response of a driven dissipative instability? What would falsify this interpretation?
Until these questions are answered, the Josephson/Adler language should remain an analogy rather than a central physical conclusion. The manuscript itself admits uncertainty about true macroscopic phase locking, but the abstract, discussion, and conclusion still give this interpretation too much strenght.
Overall assessment
The paper should not be rejected, because the observational data are valuable and the theoretical direction may be productive. However, the manuscript currently asks the reader to accept too many intermediate steps without enough demonstration.
Before making such a strong claim, the authors have to show that the proposed stochastic/RG framework follows from Farley–Buneman plasma physics, that the RG procedure is mathematically well defined, that the mechanism can be verified numerically or at least against existing simulations, that the observational scaling is not simply a common-driving effect, that radar echo count can be interpreted as turbulent-structure density, and that the Bohm and Adler/Josephson interpretations have independent physical support.
Until these points are answered, the manuscript remains an interesting interpretation of auroral turbulence, but not yet a fully established theory of it.