Composite Two-Component Turbulence Model

• The composite turbulence model is defined as a mix of parallel (slab) and perpendicular (2D) fluctuations, where the slab component increases with solar activity.
• The methodology relies on spectral decomposition of magnetic-field data using Taylor mapping and directional spectral analysis to estimate anisotropy and energy fractions.
• Solar-cycle variability modulates the slab fraction, affecting particle scattering and cosmic-ray transport through anisotropic MHD turbulence dynamics.

Searching arXiv for recent and foundational papers on the composite two-component turbulence model in the solar wind and related transport applications. arXiv search: "(Fa et al., 19 May 2025) Solar-cycle Variability of Composite Geometry in the Solar Wind Turbulence" The composite two-component turbulence model, in heliospheric plasma physics, represents solar-wind turbulence as a superposition of a slab component and a two-dimensional component defined relative to the local mean magnetic field $B_0$. In the slab component, wavevectors are predominantly parallel to $B_0$; in the 2D component, wavevectors are predominantly perpendicular to $B_0$. A long-baseline Wind analysis covering 1995 January to 2023 December found that the inertial-range turbulence at $\sim 1$ AU is dominated by the 2D component, with about $80\%$ of the turbulence energy in that mode, while the slab fraction increases with solar activity as measured by sunspot number (Fa et al., 19 May 2025). The same decomposition also underpins contemporary models of particle scattering, perpendicular diffusion, and self-consistent solar-wind turbulence transport (Shalchi, 2016, Wiengarten et al., 2016).

1. Geometric definition and spectral decomposition

In its standard solar-wind form, the model is field-aligned and axisymmetric. The slab component is defined by wavevectors parallel to the mean field, with $k_x=k_y=0$, so that $k_\parallel=k_z\neq 0$ and $k_\perp=0$. The 2D component is defined by wavevectors perpendicular to the mean field, with $k_\parallel=k_z=0$ and $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$ (Fa et al., 19 May 2025).

The corresponding spectral-tensor representation separates the two geometries cleanly. For slab turbulence,

$B_0$0

For 2D turbulence,

$B_0$1

A common inertial-range assumption is

$B_0$2

with a Kolmogorov-like discussion often taking $B_0$3. The slab fraction is then

$B_0$4

This decomposition is not merely geometric. It is a dynamical partition between parallel-propagating, field-aligned fluctuations and a perpendicular cascade that carries most of the inertial-range energy at $B_0$5 AU. In the broader heliospheric literature, the same distinction also appears as a wave-like versus quasi-2D split in fully three-dimensional transport models (Wiengarten et al., 2016).

2. Observational inference at 1 AU

The principal long-term observational implementation uses Wind magnetic-field and plasma data at $B_0$6 AU, measured from 1995 January to 2023 December with 3 s resolution by MFI and 3DP (Fa et al., 19 May 2025). The analysis is carried out in 30-minute intervals, partially overlapping, with each interval required to have $B_0$7 missing data and no consecutive gaps; small gaps are linearly interpolated. Pre-whitening and post-darkening are applied to stabilize spectral estimation, and only intervals with clear inertial-range power laws are retained.

The analysis focuses on 21 frequencies between $B_0$8 and $B_0$9 Hz, corresponding to timescales of about $B_0$0 to $B_0$1 s, with a reference frequency $B_0$2 Hz and $B_0$3. The field-to-flow angle $B_0$4 is binned into nine $B_0$5 bins from $B_0$6 to $B_0$7. Taylor’s hypothesis is used to map temporal to spatial spectra,

$B_0$8

Two complementary tests are used to infer the slab fraction. The spectrum ratio test employs directional spectra $B_0$9 and $\sim 1$0, and fits the anisotropy ratio

$\sim 1$1

as a function of $\sim 1$2. The spectrum anisotropy test instead fits the $\sim 1$3-dependence of the normalized total spectrum at $\sim 1$4, allowing the total amplitude $\sim 1$5 to vary with angle rather than assuming it constant. Empirically, $\sim 1$6 increases by a factor of about $\sim 1$7 from $\sim 1$8 to $\sim 1$9 (Fa et al., 19 May 2025).

At fixed $80\%$0, the composite model gives

$80\%$1

$80\%$2

Hence

$80\%$3

The limiting cases are diagnostically useful: $80\%$4 for pure slab and $80\%$5 for pure 2D.

3. Baseline partition and solar-cycle variability

Across the full 1995–2023 interval, the Wind analysis reports a geometric-mean anisotropy ratio $80\%$6, with annual values in the range $80\%$7 to $80\%$8, confirming significant spectral anisotropy (Fa et al., 19 May 2025). The two inference methods yield different absolute slab fractions, but both imply persistent 2D dominance.

Estimate	Overall slab fraction	Annual range
Ratio test $80\%$9	$k_x=k_y=0$0	$k_x=k_y=0$1–$k_x=k_y=0$2
Anisotropy test $k_x=k_y=0$3	$k_x=k_y=0$4	$k_x=k_y=0$5–$k_x=k_y=0$6
Average $k_x=k_y=0$7	$k_x=k_y=0$8	2D fraction $k_x=k_y=0$9

The central empirical result is the solar-cycle dependence of $k_\parallel=k_z\neq 0$0. The slab fraction rises with annual sunspot number $k_\parallel=k_z\neq 0$1, with Pearson correlations

$k_\parallel=k_z\neq 0$2

The paper reports the best-fit parameterizations

$k_\parallel=k_z\neq 0$3

The dependence is sublinear, and the 2D component remains dominant even at solar maximum.

The year-to-year behavior follows the same pattern. At solar minimum in 2018, $k_\parallel=k_z\neq 0$4, $k_\parallel=k_z\neq 0$5, and $k_\parallel=k_z\neq 0$6. During solar maxima and rise phases, slab fractions are larger: $k_\parallel=k_z\neq 0$7 reaches $k_\parallel=k_z\neq 0$8 in 1999 and $k_\parallel=k_z\neq 0$9 in 2023, while $k_\perp=0$0 reaches $k_\perp=0$1 in 2002. The mean slab fraction $k_\perp=0$2 lies around $k_\perp=0$3–$k_\perp=0$4 in those active years (Fa et al., 19 May 2025).

The same study associates slab enhancement with stronger interplanetary magnetic field magnitude and larger Alfvén speed,

$k_\perp=0$5

This suggests that stronger mean-field and Alfvénic conditions during rise phases favor a larger parallel-propagating contribution without overturning the overall predominance of the perpendicular cascade.

4. Spectral anisotropy and theoretical interpretation

In the Wind data, the directional anisotropy ratio $k_\perp=0$6 increases with $k_\perp=0$7, peaks near $k_\perp=0$8–$k_\perp=0$9, and then decreases slightly; the total spectrum $k_\parallel=k_z=0$0 increases by a factor of about $k_\parallel=k_z=0$1 from $k_\parallel=k_z=0$2 to $k_\parallel=k_z=0$3 (Fa et al., 19 May 2025). Both behaviors are consistent with a mixed slab+2D population rather than either pure limit.

The theoretical context is standard anisotropic MHD turbulence. Models emphasizing critical balance or quasi-perpendicular cascades predict $k_\parallel=k_z=0$4, which is consistent with the observed 2D dominance at $k_\parallel=k_z=0$5 AU. The measured solar-cycle modulation of $k_\parallel=k_z=0$6 adds a second layer: the anisotropic cascade persists, but the relative amount of field-aligned, wave-like content is not constant across the cycle.

A generalized solar-wind turbulence transport model extends this picture by evolving a low-frequency quasi-2D component and a high-frequency, parallel-propagating wave-like component self-consistently with the background MHD solar wind (Wiengarten et al., 2016). In that formulation, both components undergo quasi-perpendicular nonlinear cascades, but they are driven differently: stream shear contributes at low frequencies, while pickup ions can inject high-frequency wave-like power. This broader framework places the slab+2D decomposition within a dynamical, heliocentric transport model rather than treating it as a purely local spectral fit.

5. Particle transport and heliospheric applications

The composite partition matters because slab and 2D fluctuations enter transport theory differently. In quasi-linear treatments, pitch-angle diffusion and parallel mean free paths depend strongly on the slab spectrum $k_\parallel=k_z=0$7, while cross-field diffusion is influenced by the 2D spectrum $k_\parallel=k_z=0$8 and field-line random walk (Fa et al., 19 May 2025). A larger slab fraction therefore implies more efficient parallel scattering and can reduce $k_\parallel=k_z=0$9.

This dependence is sharpened by the nonlinear transport analysis of Shalchi. In two-component turbulence, slab modes do not explicitly contribute to the late-time perpendicular diffusion coefficient $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$0, but they do contribute implicitly by reducing $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$1 through a decorrelation factor $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$2 or, in an approximate form, through an effective $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$3 damping term in the denominator of the transport integral (Shalchi, 2016). In that sense, the slab fraction affects both parallel and perpendicular transport, even when only the 2D spectrum appears explicitly in the late-time $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$4 expression.

Intermittency introduces an additional layer. Test-particle simulations in static, quasi-three-dimensional composite turbulence at $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$5 AU show that large-scale intermittency can produce an extended phase of subdiffusive parallel transport during which cross-field transport diffusion dominates, and can drive the ratio $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$6 close to $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$7 (Alouani-Bibi et al., 2014). This does not replace the slab+2D partition; rather, it indicates that geometry and intermittency act together.

For practical modeling at $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$8 AU, the solar-cycle study proposes a direct parameterization: $k_\perp=(k_x^2+k_y^2)^{1/2}\neq 0$9 with either

$B_0$00

or

$B_0$01

and $B_0$02 when $B_0$03 is unavailable (Fa et al., 19 May 2025). This yields a cycle-aware geometry for cosmic-ray and SEP transport calculations.

6. Robustness, limitations, and scope

The observational case for the model at $B_0$04 AU is strengthened by the scale of the dataset: 29 years of Wind measurements, millions of 30-minute intervals, 3 s cadence, inertial-range quality control, geometric averaging within $B_0$05-bins, and two independent geometry tests (Fa et al., 19 May 2025). Within that framework, the conclusion that the 2D component dominates, while the slab fraction increases with solar activity, is robust.

The main limitations are also explicit. The inference is based on single-spacecraft sampling, Taylor’s hypothesis, axisymmetry, and field-aligned coordinates. Taylor mapping can be challenged in intervals with strong compressibility or slow flow. Uncertainties in the inertial-range slope $B_0$06 propagate into $B_0$07, and the study reports correlation coefficients but no formal confidence intervals for the slab-fraction fits. Fast- versus slow-wind dependence is not segmented in that analysis, even though prior work has found stronger slab fractions in fast wind (Fa et al., 19 May 2025).

The phrase “composite two-component turbulence model” is used in multiple research areas, but in heliospheric turbulence it has a specific meaning: a slab+2D, or wave-like+quasi-2D, decomposition tied to anisotropic MHD turbulence, solar-wind observations, and particle transport. In that domain, the model has evolved from a local spectral diagnostic into a cycle-dependent, transport-relevant parameterization that links solar activity, spectral anisotropy, and the slab/2D energy partition (Fa et al., 19 May 2025, Wiengarten et al., 2016).